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

Characterization of BCP/PreC/C region quasispecies in treatment-naive patients with different phases of HBV infection using next-generation sequencing

BACKGROUND

To analysis of quasispecies (QS) changes and high-frequency mutations in the BCP/PreC/C region of patients at different phases of hepatitis B virus (HBV) infection and provides novel biomarkers for the diagnosis of chronic hepatitis B (CHB) patients.

METHODS

With the application of next-generation sequencing technology, we were able to sequence the HBV BCP/PreC/C regions in 40 patients, each at different phases of the HBV infection. The heterogeneity of QS and the frequency of mutations were calculated using MEGA 7 software.

RESULTS

Our results show that the complexity and diversity of the BCP/PreC/C QS in HBeAg-positive CHB patients are significantly higher than those in HBeAg-positive chronic infection patients, while HBeAg-negative chronic infection patients had significantly higher QS complexity and diversity than HBeAg-negative CHB patients. In addition, HBeAg-negative patients showed reduced complexity but increased diversity compared with HBeAg-positive patients. Receiver operating characteristic curves showed that G1764A, C2102T, dN and complexity of QS could be used as potential biomarkers for diagnosing HBeAg-positive CHB, while the A2189C, dS and complexity of QS could be used as potential biomarkers for diagnosing HBeAg-negative chronic hepatitis. Finally, our study also found that G1896A and A2159G may be hotspot mutations affecting HBeAg seroconversion.

CONCLUSION

Our research elucidates the evolution of HBV by analyzing QS heterogeneity and mutation patterns, offering novel serum biomarkers for enhancing clinical diagnosis and disease prognosis. This comprehensive approach sheds light on the intricate dynamics of HBV progression and paves the way for more precise medical interventions.

The central nervous system is a potential reservoir and possible origin of drug resistance in hepatitis B infection

Background

The significance of hepatitis B virus (HBV) in cerebrospinal fluid (CSF) is unclear.

Methods

Synchronous serum and CSF samples were collected from 13 patients. HBV DNA, full-length genome, quasispecies, phylogenetic tree, compartmentalization and mutation of the reverse transcriptase (RT) region were performed based on PCR and sequencing methods.

Results

HBV DNA was detected in the CSF of 3 antiviral-naïve individuals and 1 individual after successful antiviral therapy. Complete full-length HBV genomes were isolated from the CSF of 5 individuals, including 2 with undetectable serum HBV DNA. Ten individuals exhibited distinct CSF-serum quasispecies, 8 harbored independent CSF-serum genetic compartmentalization and phylogenetic trees, and 5 lamivudine/entecavir-associated resistance mutations only in the CSF. The frequencies of rtL180M and rtM204I/V mutations in both serum and CSF were higher in HIV-HBV-coinfected individuals than in the HBV-monoinfected ones (serum: rtL180M: 3.9% vs. 0, P = 0.004; rtM204I/V: 21.3% vs. 0, P < 0.001; CSF: rtL180M: 7.6% vs. 0, P = 0.026; rtM204I/V 7.6% vs. 1.6%, P = 0.097).

Conclusion

CSF is a potential HBV reservoir, and HBV in CSF harbors distinct evolution and mutation characteristics from those in serum. HIV infection increases the possibility of HBV rtL180M and rtM204I/V mutations in both serum and CSF.

Despite low viral titer in saliva samples, Sanger-based SARS-CoV-2 spike gene sequencing is highly applicable for the variant identification

BACKGROUND

This study aimed to compare the performance of Sanger-based SARS-CoV-2 spike gene sequencing and Next Generation Sequencing (NGS)-based full-genome sequencing for variant identification in saliva samples with low viral titer.

METHODS

Using 241 stocked saliva samples collected from confirmed COVID-19 patients between November 2020 and March 2022 in Hiroshima, SARS-CoV-2 spike gene sequencing (nt22735-nt23532) was performed by nested RT-PCR and Sanger platform using in-house primers. The same samples underwent full-genome sequencing by NGS using Illumina NextSeq2000.

RESULTS

Among 241 samples, 147 were amplified by both the Sanger and the Illumina NextSeq2000 NGS, 86 by Sanger only, and 8 were not amplified at all. The overall amplification rates of Illumina NextSeq2000 NGS and Sanger were 61% and 96.7%, respectively. At low viral titer (< 103 copies/mL), Illumina NextSeq2000 NGS provided 19.2% amplification, while Sanger was 89.7% (p < 0.0001). Both platforms identified 38 wild type, 54 Alpha variants, 84 Delta variants, and 57 Omicron variants.

CONCLUSIONS

Our study provided evidence to expand the capacity of Sanger-based SARS-CoV-2 spike gene sequencing for variants identification over full-genome by Illumina NextSeq2000 NGS for mass screening. Therefore, the feasible and simple Sanger-based SARS-CoV-2 spike gene sequencing is practical for the initial variants screening, which might reduce the gap between the rapid evolution of SARS-CoV-2 and its molecular surveillance.

Using machine learning models to predict HBeAg seroconversion in CHB patients receiving pegylated interferon-α monotherapy

BACKGROUND AND OBJECTIVE

Though there are many advantages of pegylated interferon-α (PegIFN-α) treatment to chronic hepatitis B (CHB) patients, the response rate of PegIFN-α is only 30 ~ 40%. Therefore, it is important to explore predictors at baseline and establish models to improve the response rate of PegIFN-α.

METHODS

We randomly divided 260 HBeAg-positive CHB patients who were not previously treated and received PegIFN-α monotherapy (180 μg/week) into a training dataset (70%) and testing dataset (30%). The intersect features were extracted from 50 routine laboratory variables using the recursive feature elimination method algorithm, Boruta algorithm, and Least Absolute Shrinkage and Selection Operator Regression algorithm in the training dataset. After that, based on the intersect features, eight machine learning models including Logistic Regression, k-Nearest Neighbors, Support Vector Machine, Decision Tree, Random Forest, Gradient Boosting, Extreme Gradient Boosting (XGBoost), and Naïve Bayes were applied to evaluate HBeAg seroconversion in HBeAg-positive CHB patients receiving PegIFN-α monotherapy in the training dataset and testing dataset.

RESULTS

XGBoost model showed the best performance, which had largest AUROC (0.900, 95% CI: 0.85-0.95 and 0.910, 95% CI: 0.84-0.98, in training dataset and testing dataset, respectively), and the best calibration curve performance to predict HBeAg seroconversion. The importance of XGBoost model indicated that treatment time contributed greatest to HBeAg seroconversion, followed by HBV DNA(log), HBeAg, HBeAb, HBcAb, ALT, triglyceride, and ALP.

CONCLUSIONS

XGBoost model based on common laboratory variables had good performance in predicting HBeAg seroconversion in HBeAg-positive CHB patients receiving PegIFN-α monotherapy.

Precision therapy for three Chinese families with maturity-onset diabetes of the young (MODY12)

Maturity-onset diabetes of the young (MODY) is rare monogenic diabetes. However, MODY is often undiagnosed or misdiagnosed. In this study, we aimed to investigate the pathogenic gene for diabetes and provide precise treatment for diabetes patients in three families. Three families with suspected MODY were enrolled and screened for germline mutations using Whole exome sequencing (WES). Candidate pathogenic variants were validated in other family members and non-related healthy controls. Three heterozygous missense mutations in the ABCC8 gene (NM_001287174), c.1555 C>T (p.R519C), c.3706 A>G (p.I1236V), and c.2885 C>T (p.S962L) were found in families A, B, and C, respectively. All mutation sites cosegregated with diabetes, were predicted to be harmful by bioinformatics and were not found in non-related healthy controls. Two probands (onset ages, 8 and 12 years) were sensitive to glimepiride. However, an insufficient dose (2 mg/day) led to ketoacidosis. When the dosage of glimepiride was increased to 4 mg/day, blood sugar remained under control. A dose of 4 mg glimepiride daily also effectively controlled blood sugar in an adult patient 25-year-old. In addition, all patients were sensitive to liraglutide, which could control blood sugar better. These data suggest that ABCC8 was the pathogenic gene in three families with diabetes. Glimepiride (2 mg/day) was not effective in controlling blood sugar in children with ABCC8 mutations, however, 4 mg/daily glimepiride was effective in both adults and children. Moreover, liraglutide was effective in controlling blood sugar in both adults and children with ABCC8 mutations.

Characterization of mutations in the reverse transcriptase region of hepatitis B virus in treated and untreated chronic hepatitis B patients

BACKGROUND

The reverse transcriptase (RT) region of the hepatitis B virus (HBV) is the target of antiviral treatment. However, the discrepancy in RT mutations between nucleos(t)ide analogue (NA)-treated and -untreated chronic hepatitis B (CHB) patients is un clear.

METHODS

Serum samples were collected from 119 NA-treated and 135 NA-untreated patients. The sampling time was decided by the clinician. Full-length HBV RT regions were amplified using nest polymerase chain reaction. The mutations within the RT region were analysed by direct sequencing.

RESULTS

The incidence of RT mutations in treated patients was higher than that in untreated patients (p<0.05). The classic drug-resistant mutations were detected in 44.5% (53/119) of treated patients, which was significantly higher than in untreated patients (6.7% [9/135]) (p<0.05). The non-classical mutations showed their complexity and diversity in both patient groups. Multiple mutations (three or more) were more frequent in treated patients than in untreated patients (p<0.05). Several novel mutations might be related to NA resistance.

CONCLUSIONS

The selection pressures of NAs accelerated the development of RT mutations, especially within the functional domain. Mutations in the RT region occurred not only at classical sites, but also at other non-classical sites, which might be related to drug resistance and/or viral replication. The biological function and fitness of HBV isolates harbouring these novel mutations need further in vitro and in vivo verification experiments.

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.

Detection of S-HBsAg Mutations in Patients with Hematologic Malignancies

Multiple studies of hepatitis B virus (HBV) genetic variability and its relationship with the disease pathogenesis are currently ongoing, stemming from growing evidence of the clinical significance of HBV mutations. It is becoming increasingly evident that patients with hematologic malignancies may be particularly prone to a higher frequency of such mutations. The present report is the first extensive study of the prevalence of escape mutations in S-HBsAg, performed using isolates from 59 patients from hospital hematology departments with diagnoses of leukemia (n = 32), lymphoma (n = 20), multiple myeloma (n = 3), and non-tumor blood diseases (n = 4). The isolates were serologically examined for the presence of HBV markers and sequenced using either next-generation sequencing (NGS) or Sanger sequencing. Occult hepatitis B was found in 5.1% of cases. Genetic analysis of the region corresponding to S-HBsAg demonstrated an exceptionally high mutation frequency in patients with leukemias (93.4%) and lymphomas (85.0%), along with the prominent mutation heterogeneity. Additionally, more than 15 mutations in one sample were found in patients with leukemias (6.3% of cases) and lymphomas (5.0% of cases). Most of the mutations were clinically significant. The study analyzes the mutation profile of HBV in different oncohematological diseases and the frequency of individual mutations. The data strongly suggest that the NGS method, capable of detecting minor populations of HBV mutations, provides a diagnostic advantage, lays the foundation for the development of screening methods, and allows for the study of the virological and pathogenetic aspects of hepatitis B.

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.

Global Occurrence of Clinically Relevant Hepatitis B Virus Variants as Found by Analysis of Publicly Available Sequencing Data

Several viral factors impact the natural course of hepatitis B virus (HBV) infection, the sensitivity of diagnostic tests, or treatment response to interferon-α and nucleos(t)ide analogues. These factors include the viral genotype and serotype but also mutations affecting the HBV surface antigen, basal core promoter/pre-core region, or reverse transcriptase. However, a comprehensive overview of the distribution of HBV variants between HBV genotypes or different geographical locations is lacking. To address this, we performed an in silico analysis of publicly available HBV full-length genome sequences. We found that not only the serotype frequency but also the majority of clinically relevant mutations are primarily associated with specific genotypes. Distinct mutations enriched in certain world regions are not explained by the local genotype distribution. Two HBV variants previously identified to confer resistance to the nucleotide analogue tenofovir in vitro were not identified, questioning their translational relevance. In summary, our work elucidates the differences in the clinical manifestation of HBV infection observed between genotypes and geographical locations and furthermore helps identify suitable diagnostic tests and therapies.

Prevalence of Potential Resistance Related Variants Among Chinese Chronic Hepatitis B Patients Not Receiving Nucleos(T)ide Analogues

Background and Aims

Potential drug resistance (DR) related variants in the hepatitis B virus (HBV) reverse transcriptase (RT) region may be associated with the effectiveness of antiviral drugs and disease progression. The aim of this study was to investigate the prevalence and clinical characteristics of potential DR-related variants in Chinese CHB patients not receiving nucleos(t)ide analogues (NAs).

Patients and Methods

Two hundred and six untreated CHB patients from Huzhou Central Hospital in eastern China were recruited for this study. The serum DNA was extracted and the HBV RT region was amplified using nest polymerase chain reaction (nest-PCR). The 42 potential DR-related variants were analyzed by direct sequencing.

Results

Among these CHB patients, HBV genotype B and genotype C were identified in 121 (58.7%) and 85 (41.3%) patients, respectively. Potential DR-related variants were detected in 42.7% (88/206) of patients. Primary and secondary DR variants were found in 7.3% (15/206) of patients, including rtL80I/V, rtI169T, rtV173L rtL180M, rtA181T/V, rtM204I/V, and rtN236T. The variants at rt53, rt82, rt221, rt233, rt237, and rt256 were specific for genotype B, and those at rt38, rt84, rt126, rt139, rt153, rt191, rt214, rt238, and rt242 were specific for genotype C. Moreover, the variation frequency in the A-B interdomain (3.96%) was significantly higher than that in the functional domains (1.17%) and non-A-B interdomains (1.11%). Multivariate logistic regression analysis showed that lower HBV-DNA load (<10⁶ IU/mL) was an independent factor associated with potential DR-related variants in untreated CHB patients (P <0.05).

Conclusion

Potential DR-related variants were frequent and complex in untreated Chinese CHB patients. Furthermore, the variants may contribute to decreased serum HBV-DNA loads. However, the effects of potential DR-related variants on the antiviral therapy and liver disease progression require further study.

The modulation of HBsAg level by sI126T is affected by additional amino acid substitutions in the S region of HBV

The hepatitis B surface antigen (HBsAg) is a vital serum marker for hepatitis B virus (HBV) infection. Amino acid (AA) substitutions in small hepatitis B surface protein (SHBs) are known to affect HBsAg level. However, how the genetic backbones of SHBs sequences would affect the roles of a specific AA substitution on HBsAg level remains unclear. In this study, we found that sI126 had a very high substitution detection rate of 17.54% (40/228) in untreated chronic hepatitis B cohort with subgenotype C2 HBV infection. Among different substitution types at sI126, the sI126T (N = 28) was found to be associated with significantly lower serum HBsAg level. Clone sequencing revealed that sI126T-harboring SHBs sequences had varied genetic backbones with zero to nine additional AA substitutions. Thus, we constructed 24 HBsAg expression plasmids harboring sI126T without (plasmid 1, P1) or with (P2-P24) additional AA substitution(s) and studied them in the HepG2 cells. The HBsAg levels were determined by both ELISA and Western blot. In vitro experiments showed that P1 significantly reduced HBsAg level and its secretion (p < .05), however, P2-P24 showed various extracellular and intracellular HBsAg levels. No significant differences were detected among the HBsAg mRNA levels of nine representative mutant plasmids. Our findings suggest that the modulation of HBsAg level by sI126T is affected by additional AA substitution(s) in the S region of HBV. The effects of AA combination substitutions in SHBs sequences on HBsAg levels are worthwhile for more attentions in terms of HBV biology and its clinical application.