Clinical features and viral quasispecies characteristics associated with infection by the hepatitis B virus G145R immune escape mutant. Emerg Microbes Infect. 2017;6:e15. [PMID: 28325923 PMCID: PMC5378923 DOI: 10.1038/emi.2017.2]

Host immunity and HBV S gene mutation in HBsAg-negative HBV-infected patients

Background

Clinically, some patients whose HBsAg becomes negative owing to antiviral therapy or spontaneously still show a low level of HBV DNA persistence in serum. T-lymphocyte subsets, cytokine levels and HBV S gene sequences were analyzed in this study.

Methods

A total of 52 HBsAg-negative and HBV DNA-positive patients(HBsAg-/HBV DNA+ patients), 52 persistently HBsAg-positive patients(HBsAg+/HBV DNA+ patients) and 16 healthy people were evaluated. T-lymphocyte subsets of these patients were detected by flow cytometry, serum cytokines and chemokines were detected by the Luminex technique, and the HBV S region was evaluated by Sanger sequencing. T%, T-lymphocyte, CD8+ and CD4+T lymphocyte were lower in the HBsAg-negative group than in the HC group. Compared with the HBsAg-positive group, the HBsAg-negative group had lower levels in T lymphocyte %, CD8+T lymphocyte %, CD8+T lymphocyte and CD4/CD8. These difference were statistically significant (P<0.05). Serum IFN-γ, IFN-α and FLT-3L levels were significantly higher in the HBsAg-negative group than in the HBsAg-positive group (P<0.05). However, levels of many cytokines related to inflammation (i.e., IL-6, IL-8, IL10, IL-12, IL-17A) were lower in the HBsAg-negative group. Fifty-two HBsAg-negative samples were sequenced, revealing high-frequency amino acid substitution sites in the HBV S protein, including immune escape mutations (i.e., Y100C, S114T, C124Y, P127L, G130R, T131N, M133T, C137S, G145A) and TMD region substitutions (i.e., E2K/R/D, G7D/R, G10D, A17R, F20L/S, L21V, L22V).

Conclusions

According to the results of T-lymphocyte subsets and serum cytokines, it can be deduced that the cellular immune function of HBsAg-negative patients is superior to that of HBsAg-positive patients, with attenuation of liver inflammation. HBsAg-negative patients may show a variety of mutations and amino acid replacement sites at high frequency in the HBV S region, and these mutations may lead to undetectable HBsAg, HBsAg antigenic changes or secretion inhibition.

Molecular and genetic characterization of the hepatitis B virus full-length genome sequences identified in HBsAg-negative blood donors in Ural Federal District

Introduction. The World Health Organization estimates that as of 2019, more than 296 million people were living with chronic hepatitis B virus (HBV) infection. The prevalence of HBsAg-negative, occult form of the disease in blood donors varies depending on the region of the world and the sensitivity of the methods of analysis used. Considering that the genetic diversity of viruses demonstrates space and time variations and taking into account that the genetic profile of isolates in key groups, which may turn into a source of the pathogen spread, is important for forecasting of the epidemiological situation, the attention should be given to identification of HBV genotypes currently circulating among regular blood donors in regions of the Russian Federation. The aim of this work was molecular and genetic characterization of HBV genomes identified in HBsAg-negative blood donors in the Ural Federal District. Materials and methods. The study material was 1400 plasma samples obtained from HBsAg-negative blood donors in Ural Federal District. The study included the testing for HBsAg, anti-HBs IgG and anti-HBcore IgG antibodies, HBV DNA. For all identified HBV DNA containing samples, sequencing and analysis of the nucleotide sequences of the complete HBV genomes were performed. Results. The prevalence of HBV DNA was 4.93%, including 4 (0.28%) cases of false occult hepatitis B. Among anti-HBcore IgG-positive samples, HBV DNA was found in 18.08% of cases, while in persons with detected HBV DNA the anti-HBcore IgG positivity rate was 46.38%. In 8.69% of the isolates, anti-HBs IgG antibodies and viral DNA were detected simultaneously in the absence of anti-HBcore IgG. Based on phylogenetic analysis, HBV subgenotypes distribution in HBsAg-negative blood donors was as follows: D3 53.62%, D2 21.74%, D1 18.84%, C2 5.8%. The high variability in the S, C, P regions of the virus genome in the examined group was shown. In all cases of HBsAg-negative chronic HBV infection identified in blood donors, viral sequences contained at least one amino acid substitution in positions, mutations in which are associated with immune escape. In 3 (4.35%) cases mutations in reverse transcriptase region of P gene that are associated with resistance to the following drugs were identified: lamivudine, telbivudine, entecavir. Mutations in the preCore/Core regions that contribute to the progression of liver disease were also identified. Conclusion. Occult HBsAg-negative chronic HBV infection poses a threat of HBV transmission through transfusion of blood and its components due to the extremely low viral load, which does not allow the virus to be detected using routinely used diagnostic kits. The situation can be exacerbated by the abundance and diversity of virus amino acid substitutions that we have identified, including immune escape mutations, drug resistance mutations, and mutations that contribute to the progression of the disease.

T-Cell Receptor β Chain and B-Cell Receptor Repertoires in Chronic Hepatitis B Patients with Coexisting HBsAg and Anti-HBs

Antibodies in response to antigens are related to the immune repertoire of T- and B-cell receptors. However, some patients with chronic hepatitis B (CHB) have coexisting HBsAg and anti-HBsAg antibodies (anti-HBs) that cannot neutralize HBV. We attempted to investigate the repertoires that produce this response in CHB patients. The T-cell receptor β chain (TRB) and B-cell receptor (BCR) repertoires of peripheral blood genomic DNA were analyzed using MiXCR. T-cell receptor (TCR) cluster analysis was carried out by clusTCR, and motifs prediction was selected by Multiple Em for Motif Elicitation (MEME). A total of 76 subjects were enrolled, including 26 HBsAg and anti-HBs coexisting patients with CHB (DP group), 25 anti-HBs single-positive healthy people (SP group), and 25 CHB patients (CHB group). The clone length of BCR in 39, 90 was significantly different among these groups (p = 0.005, 0.036). The motif “CASSLG” in the DP group was significantly higher than SP and CHB groups and may relate to coexistence, and the motif “GAGPLT” was only shown in the SP group and may relate to anti-HB expression. These provide important insights into vaccine development and CHB treatment.

T-Cell Receptor β Chain and B-Cell Receptor Repertoires in Chronic Hepatitis B Patients with Coexisting HBsAg and Anti-HBs

Antibodies in response to antigens are related to the immune repertoire of T- and B-cell receptors. However, some patients with chronic hepatitis B (CHB) have coexisting HBsAg and anti-HBsAg antibodies (anti-HBs) that cannot neutralize HBV. We attempted to investigate the repertoires that produce this response in CHB patients. The T-cell receptor β chain (TRB) and B-cell receptor (BCR) repertoires of peripheral blood genomic DNA were analyzed using MiXCR. T-cell receptor (TCR) cluster analysis was carried out by clusTCR, and motifs prediction was selected by Multiple Em for Motif Elicitation (MEME). A total of 76 subjects were enrolled, including 26 HBsAg and anti-HBs coexisting patients with CHB (DP group), 25 anti-HBs single-positive healthy people (SP group), and 25 CHB patients (CHB group). The clone length of BCR in 39, 90 was significantly different among these groups (p = 0.005, 0.036). The motif “CASSLG” in the DP group was significantly higher than SP and CHB groups and may relate to coexistence, and the motif “GAGPLT” was only shown in the SP group and may relate to anti-HB expression. These provide important insights into vaccine development and CHB treatment.

Protective Role of Rheumatic Diseases Against Hepatitis B Virus Infection and Human Leukocyte Antigen B27 Highlighted

BACKGROUND

By determining the hepatitis B virus (HBV) surface antigen (HBsAg) positive rate postexposure and HBV-specific antigen/antibody (Ag/Ab) level in patients with rheumatic diseases, we aimed at exploring the rheumatic link to HBV control.

METHODS

Patients who underwent HBV screening in the Ruijin Hospital from 2020 to 2021 were enrolled for the exposure rate estimation. Among antibody to HBV core antigen (HBcAb)-positive patients, we adopted propensity score matching (PSM) to study the impact of rheumatism on HBsAg seroprevalence after exposure. A second PSM evaluated the Ag/Ab differences. We also had HBsAg prevalence in human leukocyte antigen B2 (HLA-B27) tested patients studied.

RESULTS

With 33,989 screened patients, exposure rates remained comparable between rheumatic and non-rheumatic patients: 48.94 vs. 49.86%. PSM first yielded 2,618 balanced pairs. We observed significantly fewer patients with rheumatic diseases in HBsAg positive cases than negative ones (p < 0.001). In the second round, PSM matched 279 pairs, HBsAg (p < 0.001) and HBeAg (p < 0.05) positivity rates were significantly lower in the rheumatic patients, whereas HBsAb positivity rate (p < 0.001) and level (p < 0.01) were significantly higher. Though the value of HBcAb was overall significantly lower (p < 0.001) within the realm of rheumatic diseases, patients with ankylosing spondylitis (AS) demonstrated a significantly higher value than other rheumatic diseases. We saw significantly fewer HBV infections in HLA-B27 positive subjects than in the negative ones (p < 0.001).

CONCLUSION

In this propensity score-matched study, rheumatic patients had an advantage in HBV control. In rheumatic patients, HBcAb levels, together with the beneficial role of HLA-B27, were highlighted.

The Impact of HBV Quasispecies Features on Immune Status in HBsAg+/HBsAb+ Patients With HBV Genotype C Using Next-Generation Sequencing

Background

This study aimed to explore the molecular mechanism of the coexistence of hepatitis B surface antigen (HBsAg) and hepatitis B surface antibody (HBsAb) serological pattern via intensive characterization of HBV s gene in both chronic hepatitis B (CHB) and hepatocellular carcinoma (HCC) patients.

Method

A total of 73 HBsAg+/HBsAb+ patients (CHB = 36, HCC = 37) and 96 HBsAg+/HBsAb− patients (CHB = 47, HCC = 49) were enrolled from 13 medical centers in China. The sequence features were elaborated based on the combination of next-generation sequencing (NGS) and multidimensional bioinformatics analysis.

Results

The 16 high-frequency missense mutations, changes of stop codon mutation, clustering, and random forest models based on quasispecies features demonstrated the significant discrepancy power between HBsAg+/HBsAb+ and HBsAg+/HBsAb− in CHB and HCC, respectively. The immunogenicity for cytotoxic T lymphocyte (CTL) epitope Se and antigenicity for the major hydrophilic region (MHR) were both reduced in HBsAg+/HBsAb+ patients (CTL Se: p < 0.0001; MHR: p = 0.0216). Different mutation patterns were observed between HBsAg+/HBsAb+ patients with CHB and with HCC. Especially, mutations in antigenic epitopes, such as I126S in CHB and I126T in HCC, could impact the conformational structure and alter the antigenicity/immunogenicity of HBsAg.

Conclusion

Based on NGS and bioinformatics analysis, this study indicates for the first time that point mutations and quasispecies diversities of HBV s gene could alter the MHR antigenicity and CTL Se immunogenicity and could contribute to the concurrent HBsAg+/HBsAb+ with different features in HCC and CHB. Our findings might renew the understanding of this special serological profile and benefit the clinical management in HBV-related diseases.

Composition and Dynamics of H1N1 and H7N9 Influenza A Virus Quasispecies in a Co-infected Patient Analyzed by Single Molecule Sequencing Technology

A human co-infected with H1N1 and H7N9 subtypes influenza A virus (IAV) causes a complex infectious disease. The identification of molecular-level variations in composition and dynamics of IAV quasispecies will help to understand the pathogenesis and provide guidance for precision medicine treatment. In this study, using single-molecule real-time sequencing (SMRT) technology, we successfully acquired full-length IAV genomic sequences and quantified their genotypes abundance in serial samples from an 81-year-old male co-infected with H1N1 and H7N9 subtypes IAV. A total of 26 high diversity nucleotide loci was detected, in which the A-G base transversion was the most abundant substitution type (67 and 64%, in H1N1 and H7N9, respectively). Seven significant amino acid variations were detected, such as NA:H275Y and HA: R222K in H1N1 as well as PB2:E627K and NA: K432E in H7N9, which are related to viral drug-resistance or mammalian adaptation. Furtherly, we retrieved 25 H1N1 and 22 H7N9 genomic segment haplotypes from the eight samples based on combining high-diversity nucleotide loci, which provided a more concise overview of viral quasispecies composition and dynamics. Our approach promotes the popularization of viral quasispecies analysis in a complex infectious disease, which will boost the understanding of viral infections, pathogenesis, evolution, and precision medicine.

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.

Investigation of immune escape-associated mutations of hepatitis B virus in patients harboring hepatitis B virus drug-resistance mutations

BACKGROUND

It is unclear whether immune escape-associated mutations in the major hydrophilic region of hepatitis B virus surface antigen (HBsAg) are associated with nucleoside/nucleotide analog resistance.

AIM

To evaluate the association between immune escape-associated mutations and nucleoside/nucleotide analog resistance mutations.

METHODS

In total, 19440 patients with chronic hepatitis B virus infection, who underwent resistance testing at the Fifth Medical Center of Chinese PLA General Hospital between July 2007 and December 2017, were enrolled. As determined by sequence analysis, 6982 patients harbored a virus with resistance mutations and 12458 harbored a virus lacking resistance mutations. Phenotypic analyses were performed to evaluate HBsAg production, replication capacity, and drug-induced viral inhibition of patient-derived drug-resistant mutants with or without the coexistence of sA159V.

RESULTS

The rate of immune escape-associated mutation was significantly higher in 9 of the 39 analyzed mutation sites in patients with resistance mutations than in patients without resistance mutations. In particular, these mutations were sQ101H/K/R, sS114A/L/T, sT118A/K/M/R/S/V, sP120A/L/Q/S/T, sT/I126A/N/P/S, sM133I/L/T, sC137W/Y, sG145A/R, and sA159G/V. Among these, sA159V was detected in 1.95% (136/6982) of patients with resistance mutations and 1.08% (134/12,458) of patients lacking resistance mutations (P < 0.05). The coexistence of sA159V with lamivudine (LAM) and entecavir (ETV)-resistance mutations in the same viral genome was identified during follow-up in some patients with drug resistance. HBsAg production was significantly lower and the replication capacity was significantly higher, without a significant difference in LAM/ETV susceptibility, in sA159V-containing LAM/ETV-resistant mutants than in their sA159V-lacking counterparts.

CONCLUSION

In summary, we observed a close link between the increase in certain immune escape-associated mutations and the development of resistance mutations. sA159V might increase the fitness of LAM/ETV-resistant mutants under environmental pressure in some cases.

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.

Investigation of hepatitis B virus (HBV) rtS78T/sC69* mutation in a large cohort of chronic HBV-infected patients with nucleoside/nucleotide analogue treatment

This study aimed to investigate clinical occurrence and significance of the rtS78T/sC69* mutation of hepatitis B virus (HBV). A total of 22,009 consecutive chronic HBV-infected patients who underwent resistance testing at the Fifth Medical Center of Chinese PLA General Hospital (Original name Beijing 302 Hospital) from 2007 to 2016 were enrolled. Serum samples were collected for sequence analysis of HBV reverse-transcriptase (RT) and S regions. Phenotypic analysis was performed to evaluate the viral replication capacity and drug susceptibility. The rtS78T mutation was detected in 0.83% (182/22,009) of the patients' samples. All mutations simultaneously created a stop codon at sC69 (sC69*). The prevalence of rtS78T/sC69* did not differ significantly between the patients with and without entecavir/tenofovir treatment. Of the 182 mutation-positive samples, 41 (22.5%) were detected with signature drug-resistance mutations to adefovir (n = 26), lamivudine (n = 11), entecavir (n = 3), and lamivudine plus adefovir (n = 1). The HBV DNA and RNA levels of the rtS78T/sC69* mutant were significantly increased compared to the wild-type; while the mutant had undetectable secreted and intracellular HBsAg, and its half maximal effective concentration to lamivudine, adefovir, entecavir, and tenofovir were 3.73-, 1.61-, 4.76-, and 3.71-fold of the wild-type, respectively. Artificial elimination of the rtS78T mutation had a limited effect on the drug susceptibilities. The data obtained in the present study suggested that the emergence of the rtS78T/sC69* mutation was not closely related to entecavir/tenofovir treatment and itself appeared insufficient to confer drug resistance unless it coexisted with signature drug-resistance mutations.

Hepatitis B virus mutation pattern rtL180M+A181C+M204V may contribute to entecavir resistance in clinical practice

Background and Aims: Entecavir (ETV) resistance of hepatitis B virus (HBV) conventionally requires rt184, 202, or 250 mutations plus lamivudine-resistance mutation (rtM204V/I ± L180M). This study aimed to clarify whether rtL180M+A181C+M204V mutations may contribute to HBV ETV resistance.

Methods: Serum samples were collected from 22,009 patients who underwent resistance testing in Beijing 302 Hospital from 2007 to 2016. HBV reverse transcriptase (RT) gene was screened by direct sequencing and verified by clonal sequencing. Phenotypic analysis was performed for evaluating replication capacity and drug susceptibility.

Results: Classical ETV-resistance mutations of HBV were detected in 1252 patients who were receiving ETV therapy. The rtA181C mutation was detected with rtL180M+M204V mutations in 18 lamivudine-experienced ETV-treated patients, and the emergence of the mutations was associated with virological breakthrough or inadequate virological response to ETV. Patient-derived representative rtA181C-containing mutants, rtL180M+A181C+M204V, rtL180M+A181C+M204V+M250V, and rtL180M+A181C+S202G+M204V, exhibited 45.7%, 25.9%, and 25.0% replication capacity and 85.6-, 356.1-, and 307.1-fold decreased susceptibility to ETV respectively compared to the wild-type strain, while the three mutants remained sensitive to tenofovir (TDF). Artificial elimination of rtA181C largely restored the rtL180M+A181C+M204V mutant’s sensitivity to ETV. Molecular modelling of viral RT binding to ETV showed that the rtL180M+A181C+M204V mutant had a less stable conformation compared to rtL180M+M204V mutant. In clinical practice, undetectable serum HBV DNA was achieved in two of five longitudinally followed rtA181C-positive patients who received switching-to TDF therapy, but not in the other three who received add-on adefovir therapy during observation.

Conclusions: Both clinical and experimental data support rtL180M+A181C+M204V as a novel non-classical ETV-resistance mutation pattern.

Application of deep sequencing methods for inferring viral population diversity

The first deep sequencing method was announced in 2005. Due to an increasing number of sequencing data and a reduction in the costs of each sequencing dataset, this innovative technique was soon applied to genetic investigations of viral genome diversity in various viruses, particularly RNA viruses. These deep sequencing findings documented viral epidemiology and evolution and provided high-resolution data on the genetic changes in viral populations. Here, we review deep sequencing platforms that have been applied in viral quasispecies studies. Further, we discuss recent deep sequencing studies on viral inter- and intrahost evolution, drug resistance, and humoral immune selection, especially in emerging and re-emerging viruses. Deep sequencing methods are becoming the standard for providing comprehensive results of viral population diversity, and their applications are discussed.

Identified OAS3 gene variants associated with coexistence of HBsAg and anti-HBs in chronic HBV infection

The underlying mechanism of coexistence of hepatitis B surface antigen (HBsAg) and hepatitis B surface antigen antibody (anti-HBs) is still controversial. To identify the host genetic factors related to this unusual clinical phenomenon, a two-stage study was conducted in the Chinese Han population. In the first stage, we performed a case-control (1:1) age- and gender-matched study of 101 cases with concurrent HBsAg and anti-HBs and 102 controls with negative HBsAg and positive anti-HBs using whole exome sequencing. In the second validation stage, we directly sequence the 16 exons on the OAS3 gene in two dependent cohorts of 48 cases and 200 controls. Although, in the first stage, a genome-wide association study of 58,563 polymorphism variants in 101 cases and 102 controls found no significant loci (P-value ≤ .05/58563), and neither locus achieved a conservative genome-wide significance threshold (P-value ≤ 5e-08), gene-based burden analysis showed that OAS3 gene rare variants were associated with the coexistence of HBsAg and anti-HBs. (P-value = 4.127e-06 ≤ 0.05/6994). A total of 16 rare variants were screened out from 21 cases and 3 controls. In the second validation stage, one case with a stop-gained rare variant was identified. Fisher's exact test of all 149 cases and 302 controls showed that the rare coding sequence mutations were more frequent in cases vs controls (P-value = 7.299e-09, OR = 17.27, 95% CI [5.01-58.72]). Protein-coding rare variations on the OAS3 gene are associated with the coexistence of HBsAg and anti-HBs in patients with chronic HBV infection in Chinese Han population.

Impact of deletions and mutations in Hepatitis B virus envelope proteins on serological profile and clinical evolution

The Hepatitis B virus (HBV) envelope glycoproteins are essential for viral entry into the hepatocyte and are also targets for host immune response. The study of these proteins could allow us to highlight molecular hot points influencing HBV fitness, which would subsequently modify the clinical evolution of the disease, both under anti-viral therapy or without treatment. The present short communication underlines the importance of the high variability in HBV envelope proteins, in regard with the literature and in our hands, for HBV-infected patients either on anti-HBV treatment or not. We report mutations in antigenic areas of S protein, i.e. CD8⁺/CD4⁺ T-cell epitopes and B-cell epitopes in the major hydrophilic region (MHR), such as sI126N and sG145R possibly involved in the rare coexisting Hepatitis B surface Antigen (HBsAg)/anti-HBs serological pattern. We mostly report serial mutations in preS region including preS1 deletion (aa 1-6, 31-71, 38-73, 72-104) and preS2 deletion (aa132-141) in patients with various clinical evolutions. Some of these viral envelope mutations, due to immune selection pressure, may result in a worsening of the hepatic disease.