Comparison of Detection Rate and Mutational Pattern of Drug-Resistant Mutations Between a Large Cohort of Genotype B and Genotype C Hepatitis B Virus-Infected Patients in North China

Effect of mutations across reverse transcriptase region on HBV replication and progression of liver diseases in Chinese patients

It was known that mutations in the RT region were mainly related to nucleot(s)ide analogs resistance. Increasing studies indicated that RT mutations were related to advanced liver diseases (ALD) and had effects on HBV replication, but the distribution characteristics of mutations across RT region in the development of liver diseases and the effect of RT mutations on HBV replication were not fully clarified. HBV RT region was direct-sequenced in 1473 chronic HBV-infected patients. Mutation frequencies were analyzed to identify the specific mutations differing between groups classified by genotypes, loads of HBV DNA, or progression of liver diseases. In the range of rt145-rt290, rt145, rt221, rt222, rt267, and rt271 were the genotype-polymorphic sites, while rt238 was the genotype-specific sites. Mutations at rt163, rt173, rt180, rt181, rt184, rt191, rt199, and rt214 were more frequent among patients with C-genotype HBV, while those at rt220, rt225, rt226, rt269, and rt274 were more frequent among patients with B-genotype HBV. RtM204V/I could reduce the HBV DNA loads while rtQ/L267H/R could increase the HBV DNA loads. RtV214A/E/I (OR 3.94, 95% CI 1.09 to 14.26) was an independent risk factor for advanced liver diseases. In summary, the hotspots of mutations were different between B and C genotypes. Besides the effect on the S region, RT mutations had effects on HBV replication by other unknown ways. RtV214A/E/I was found to be an independent risk factor for ALD, suggesting that mutations at rt214 site could be used as a potential virological marker for the liver disease progression.

Molecular Epidemiological Characteristics and Risk Factors for Acquiring HBV Among Li Ethnic in Baisha County, Hainan Island-Subgenotype D3 Was First Discovered in China

The residents of Baisha, a county of Hainan Island, mainly composed of Li ethnic population and relatively closed living environment with its unique geographical location. Our previous study showed that Li ethnic population of Baisha is an endemic center for hepatitis C virus, with significantly higher rates than in other parts of China. However, the epidemiology of HBV in this region remains unclear. Therefore, we conducted a comprehensive epidemiological survey of HBV in Baisha County, including 1,682 Li ethnic residents. The total seropositive rate for HBsAg was 10.2% and was higher than other parts of China. HBV-positive status was associated with the 20–40-year-old group (OR = 1.27, 95%CI 1.04–1.39, P < 0.01) and alcohol consumption (OR = 2.17, 95%CI 1.58–2.99, P < 0.01). Phylogenetic analysis showed that HBV subgenotype D3 was predominant in Baisha County which was first discovered in China, followed by C5, C1, B2, and undetermined subgenotypes which were significantly different from other geographical distribution of main genotypes in China. The most recent common ancestor (tMRCA) of the HBV genotype C in the Li ethnic of Baisha County was 1846 (95%CI: 1739–1932), and Baisha-C5 was earlier than Baisha-C1 and Baisha-C2. Most Baisha-D3 sequences were concentrated in one bundle and unrelated to those D3 genome sequences elsewhere in the world. According to the phylogenetic tree, D3 was introduced into Baisha County in 1884 (95%CI: 1816–1993) and became a local endemic virus. In conclusion, HBV infection in the Li ethnic group is characterized by a high prevalence rate in 20–40-year-old individuals and a unique genotype distribution which were significantly different from other geographical distribution of main genotypes in China, and subgenotype D3 was first discovered in China.

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.

Global prevalence and phylogeny of hepatitis B virus (HBV) drug and vaccine resistance mutations

Vaccination and anti-viral therapy with nucleos(t)ide analogues (NAs) are key approaches to reducing the morbidity, mortality and transmission of hepatitis B virus (HBV) infection. However, the efficacy of these interventions may be reduced by the emergence of drug resistance-associated mutations (RAMs) and/or vaccine escape mutations (VEMs). We have assimilated data on the global prevalence and distribution of HBV RAMs/VEMs from publicly available data and explored the evolution of these mutations. We analysed sequences downloaded from the HBV Database and calculated prevalence of 41 RAMs and 38 VEMs catalogued from published studies. We generated maximum likelihood phylogenetic trees and used treeBreaker to investigate the distribution and estimated the age of selected mutations across tree branches. RAM M204I/V had the highest prevalence, occurring in 3.8% (109/2838) of all HBV sequences in our data set, and a significantly higher rate in genotype C at 5.4% (60/1102, p = 0.0007). VEMs had an overall prevalence of 1.3% (37/2837) and had the highest prevalence in genotype C and in Asia at 2.2% (24/1102; p = 0.002) and 1.6% (34/2109; p = 0.009), respectively. Phylogenetic analysis suggested that RAM/VEMs can arise independently of treatment/vaccine exposure. In conclusion, HBV RAMs/VEMs have been found globally and across genotypes, with the highest prevalence observed in genotype C. Screening for genotype and for resistance-associated mutations may help to improve stratified patient treatment. As NAs and HBV vaccines are increasingly being deployed for HBV prevention and treatment, monitoring for resistance and advocating for better treatment regimens for HBV remains essential.

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.

Molecular cloning and phenotypic analysis of drug-resistance mutants with relevant S-region variants of HBV for a patient during 189-month anti-HBV treatment

BACKGROUND

A unique chronic hepatitis B patient was followed over 189 months of nucleoside/nucleotide analogue (NA) therapies with the analysis of multiple drug-resistance HBV mutants.

METHODS

Clonal sequencing (≥20 clones/sample) was performed on sera sampled at 41 time points, and the phenotypic features of eight representative mutants were analysed.

RESULTS

Lamivudine (LAM)-, adefovir dipivoxil (ADV)-, entecavir (ETV)- and repeat ADV-resistance mutants emerged upon individual sequential NA monotherapy. The efficacy of NA combination rescue therapies ranked as LAM+ADV < ETV+ADV < ETV+ tenofovir disoproxil fumarate (TDF). Specifically, LAM+ADV and ETV+ADV suppressed viral loads to <100 IU/ml for a long period of time, either with or without late stage HBV DNA fluctuations. Furthermore, ETV+TDF suppressed the viral load to <10 IU/ml. During the LAM+ADV and ETV+ADV combination therapies, ETV-resistance mutants dominated at most time points, and multidrug-resistance (MDR) mutants that harboured LAM-, ETV- and ADV-resistance mutations were intermittently detected. Interestingly, the rtA181T-causative sW172stop to sW172non-stop mutation transition was observed at HBV DNA fluctuations. In a phenotypic analysis, two MDR strains had cross-resistance to LAM, ETV and ADV, and a lower susceptibility to TDF (<10-fold decrease compared with the wild-type strain). In contrast, the natural replication capacity was inversely associated with the number of primary resistant mutations which would limit MDR mutant development.

CONCLUSIONS

Taken together, viral drug susceptibility, replication capacity, and perhaps immunological adaptation may play coordinated roles in the fitness of drug-resistance mutants. ETV+TDF therapy is the preferred option for treating chronic hepatitis B patients with multiple drug failure.

Characteristics and evolutionary history of hepatitis B virus quasi-subgenotype B3 in Southeast Asia

To analyze the hepatitis B virus (HBV) quasi-subgenotype B3 characters and molecular evolution in Southeast Asia, 411 serum samples with HBsAg positive were collected from Xishuangbanna, China. After DNA extraction, PCR amplification and sequencing, a total of 183 HBV full-length genomes were obtained. Phylogenetic analysis showed 139 stains (76.0%) were genotype B, 41 strains were genotype C (22.4%) and 3 strains were genotype I (1.6%). Among genotype B, 34 sequences were identified as quasi-subgenotype B3. Quasi-subgenotype B3 sequences from this study and quasi-subgenotype B3 sequences from the GenBank (total of 141 complete genome sequences) were grouped into quasi-subgenotype B3 (B3, formerly B5, Chinese B6 and B7-9). Sixteen peculiar nucleotides distributed in quasi-subgenotype B3 were identified, which were differ from B1, B2, B4 and B5(formerly B6) (nt93 T, nt100C, nt355 G, nt843 T, nt861C, nt912C, nt929 T, nt930 G, nt1023 T, nt1041 T, nt2651C, nt2693 T, nt2970C, nt3054A, nt3087A and nt3171 G). Then Evolutionary dynamics analysis of HBV quasi-subgenotype B3 was performed. The mean rate of nucleotide substitution for HBV quasi-subgenotype B3 was estimated to be around 5.556-5.660 × 10-4 substitutions/site/year. Estimated time to most recent ancestor of quasi-subgenotype B3 was around the 1847-1945(95%HPD), and Yunnan strains might be the parental strains. The Bayesian sky plot showed a steady spreading of HBV quasi-genotype B3 from early of 1940s to 90 s. In summary, HBV quasi-subgenotype B3 infection is prevalent in Southeast Asia based on the current reports and still with a high prevalence rate based on the evolutionary dynamics analysis. Current vaccine and nucleotide analogues might have effective prevention and treatment for HBV quasi-subgenotype B3 based on the rare clinically relevant mutation sites included in quasi-subgenotype B3.

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.

Antiviral drug resistance and hepatitis B: a continuing public health problem

The use of selective antiviral therapy has been very successful in controlling HBV replication in individuals, leading to a reduction in disease progression and mortality. However, the use of first-generation therapies, often the only available option in low-resource settings, can result in a high prevalence of drug resistance mutations. Variants selected by antiviral therapies targeting the viral polymerase can also result in variants in the viral envelope. These variants can allow the virus to escape the host immune response. The effect of antiviral selection pressure on viral variants that may contribute to immune escape requires further investigation.

Trends in hepatitis B virus resistance to nucleoside/nucleotide analogues in North China from 2009-2016: A retrospective study

Nucleos(t)ide analogues (NAs) are widely used in anti-hepatitis B virus (anti-HBV) therapy for effective inhibition of HBV replication. However, HBV resistance to NAs has emerged, resulting in virus reactivation and disease recurrence. Data on the current dynamics of HBV resistance are still rare in China. This study analysed 4491 plasma samples with HBV primary genotypic resistance mutations representative of the general HBV resistance situation in northern China from 2009-2016. We found that entecavir (ETV), representing 57.6% (12 713/22 060) of NA users in North China in 2016, has become the major NA for treating Chinese patients infected with HBV. Despite >50% of M204I/V±L180M among all HBV resistance cases annually and extensive exposure of patients to lamivudine (LAM), telbivudine (LdT) and adefovir dipivoxil (ADV), ETV resistance also showed a dramatically increased incidence, which rose to 17.1% in 2016. Moreover, A181T/V, ETV resistance mutations and multidrug resistance mutations were found more frequently in HBV genotype C compared with genotype B (21.2% vs. 8.5%, 12.4% vs. 7.9% and 5.9% vs. 3.0%, respectively), whereas M204I and N236T were more predominant in genotype B than genotype C (40.3% vs. 20.8% and 11.3% vs. 1.8%, respectively). In conclusion, we report the dynamic changes of HBV NA resistance mutation patterns and the current NA usage profile for anti-HBV treatment in North China over the past 8 years. These data provide valuable information on HBV NA resistance that is an important reference for clinicians to devise more effective treatment regimens for individual patients.

Evolution of drug-resistant mutations in HBV genomes in patients with treatment failure during the past seven years (2010-2016)

The objective of this study was to analyze the prevalence of drug-resistant HBV mutants in patients with treatment failure during the past seven years (2010-2016). 4055 HBV-infected patients who underwent HBV polymerase gene mutation test from 2010 to 2016 were enrolled. The nucleos(t)ide analogues (NAs) resistance mutation positions, including rtL180, rtA181, rtT184, rtS202, rtM204, rtI233, rtN236, rtI169, rtV173, and rtM250 were analyzed. Genotypic resistance mutations were detected in 30.8% (1248/4055) of the patients with treatment failure. Rates of drug-resistant mutations associated with LAM, ADV, ETV, and multidrug were 27.23% (1104/4055), 9.67% (392/4055), 3.69% (150/4055), and 0.79% (32/4055). Among the primary NA-resistant mutations, rtM204I (13.44%, 545/4055) occurred more frequently, followed by rtM204V, rtN236T, rtA181T, and rtA181V. For single-base mutations, rtL180M and rtA181V increased gradually during the past seven years, while rtM204I/V and rtN236T decreased after 2015. The development of drug-resistant mutations positively correlated with the consumption of ETV (r = 0.964, P = 0.002), and weakly correlated with that of LAM (r = 0.679, P = 0.109) and ADV (r = 0.429, P = 0.354). Moreover, single-base mutation rtA181V and multi-base mutations (rtL180M + M204I and rtL180M + M204V + M204I) were more common in HBV genotype C than those in genotype B (1.94% vs. 0.66%, 1.84% vs. 0.16%, 1.02% vs. 0.16%, respectively). NA-related mutations in HBV RT region increased in the past seven years, especially for LAM. Frequencies of rtL180M and rtA181T/V increased gradually in the past seven years, to which we should pay more attention.

Molecular evolution of hepatitis B vaccine escape variants in China, during 2000-2016

Hepatitis B vaccine escape variants are the main threat to hepatitis B virus (HBV) infection in vaccination era worldwide. With 215 genotype B HBV and 313 genotype C HBV vaccine escape variants isolated from China during 2000-2016, we reported that genotype B HBV vaccine escape strains diverged in ∼1997 (95% HPD; 1987-2005), while genotype C HBV vaccine escape strains diverged in ∼1976 (95% HPD; 1955-2003). Additionally, the p-distance of genotype C HBV vaccine escape strains was 0.0291±0.0169, which was significantly higher than that in the genotype B HBV (t=131.02, p<0.05). However, genotype B HBV vaccine escape strains evolved more rapidly than genotype C HBV (2.103×10^-3 vs 1.083×10^-3 substitutions/site/year). Bayesian skyline plot analysis showed that the populations of genotype C HBV vaccine escape strains fluctuated more than those in genotype B HBV. Four sites (A5T/S, L21S, T/A126S and T/N131I/A) and 13 sites (N3S, T5A, G10Q/R/E, L21S, T47K/A/V, L98V/P, I/S126N/V/T, Q129H/R/L, T131P/I/N/A, G145A/R, L175S/F, L213I/S, V224A/G) were found to be under positive selection in genotype B and C HBV vaccine escape strains, respectively. More importantly, N3S, L21S, T47K, L98V, I/S126T and L213I mutations were detected in 1 (2.5%), 1 (2.5%), 1 (2.5%), 3 (7.5%), 1 (2.5%), 1 (2.5%) genotype C HBV infected Chinese younger with neonatal HBV vaccination, respectively. Therefore, our results should be valuable in further understanding the molecular evolution of HBV and providing new ideas for the elimination of HBV infection.

Relationship between pre-existing resistance mutations and HBV genotype in patients with HBV related liver cirrhosis

AIM

To investigate the relationship between hepatitis B virus (HBV)-related drug resistance mutations and HBV genotype in patients with HBV-related liver cirrhosis (LC).

METHODS

Serum samples were collected from 112 patients with HBV-related LC who did not receive prior anti-viral treatment. HBV reverse transcriptase regions were amplified by PCR, and PCR products were directly sequenced.

RESULTS

Among 112 HBV-related LC patients, 13 (11.61%, 13/112) had lamivudine resistance associated mutations, 15 (13.39%, 15/112) had adefovir resistance associated mutations, and 1 had mutations associated with both lamivudine and adefovir resistance. All of these resistance mutations occurred in patients infected with HBV genotype C, and there was no mutation in patients infected with HBV genotype B.

CONCLUSION

Adefovir/lamivudine resistance associated mutations exist in patients with HBV-related LC, and the mutations are related to HBV genotype. Compared with patients infected with HBV genotype B, those infected with HBV genotype C are prone to resistance mutations.