Low-level persistence of drug resistance mutations in hepatitis B virus-infected subjects with a past history of Lamivudine treatment.

We sought to determine the prevalence of hepatitis B virus (HBV) lamivudine (LAM)-resistant minority variants in subjects who once received LAM but had discontinued it prior to virus sampling. We performed direct PCR Sanger sequencing and ultradeep pyrosequencing (UDPS) of HBV reverse transcriptase (RT) of plasma viruses from 45 LAM-naive subjects and 46 LAM-experienced subjects who had discontinued LAM a median of 24 months earlier. UDPS was performed to a depth of ∼3,000 reads per nucleotide. Minority variants were defined as differences from the Sanger sequence present in ≥0.5% of UDPS reads in a sample. Sanger sequencing identified ≥1 LAM resistance mutations (rtL80I/V, rtM204I, and rtA181T) in samples from 5 (11%) of 46 LAM-experienced and none of 45 LAM-naive subjects (0%; P = 0.06). UDPS detected ≥1 LAM resistance mutations (rtL80I/V, rtV173L, rtL180M, rtA181T, and rtM204I/V) in 10 (22%) of the 46 LAM-experienced subjects, including 5 in whom LAM resistance mutations were not identified by Sanger sequencing. Overall, LAM resistance mutations were more likely to be present in LAM-experienced (10/46, 22%) than LAM-naive subjects (0/45, 0%; P = 0.001). The median time since LAM discontinuation was 12.8 months in the 10 subjects with a LAM resistance mutation compared to 30.5 months in the 36 LAM-experienced subjects without a LAM resistance mutation (P < 0.001). The likelihood of detecting a LAM resistance mutation was significantly increased using UDPS compared to Sanger sequencing and was inversely associated with the time since LAM discontinuation.

A classification model for G-to-A hypermutation in hepatitis B virus ultra-deep pyrosequencing reads.

MOTIVATION: G → A hypermutation is an innate antiviral defense mechanism, mediated by host enzymes, which leads to the mutational impairment of viruses. Sensitive and specific identification of host-mediated G → A hypermutation is a novel sequence analysis challenge, particularly for viral deep sequencing studies. For example, two of the most common hepatitis B virus (HBV) reverse transcriptase (RT) drug-resistance mutations, A181T and M204I, arise from G → A changes and are routinely detected as low-abundance variants in nearly all HBV deep sequencing samples. RESULTS: We developed a classification model using measures of G → A excess and predicted indicators of lethal mutation and applied this model to 325 920 unique deep sequencing reads from plasma virus samples from 45 drug treatment-naïve HBV-infected individuals. The 2.9% of sequence reads that were classified as hypermutated by our model included most of the reads with A181T and/or M204I, indicating the usefulness of this model for distinguishing viral adaptive changes from host-mediated viral editing. AVAILABILITY: Source code and sequence data are available at http://hivdb.stanford.edu/pages/resources.html. CONTACT: ereuman@stanfordalumni.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Similar prevalence of low-abundance drug-resistant variants in treatment-naive patients with genotype 1a and 1b hepatitis C virus infections as determined by ultradeep pyrosequencing.

BACKGROUND AND OBJECTIVES: Hepatitis C virus (HCV) variants that confer resistance to direct-acting-antiviral agents (DAA) have been detected by standard sequencing technology in genotype (G) 1 viruses from DAA-naive patients. It has recently been shown that virological response rates are higher and breakthrough rates are lower in G1b infected patients than in G1a infected patients treated with certain classes of HCV DAAs. It is not known whether this corresponds to a difference in the composition of G1a and G1b HCV quasispecies in regards to the proportion of naturally occurring DAA-resistant variants before treatment. METHODS: We used ultradeep pyrosequencing to determine the prevalence of low-abundance (<25% of the sequence reads) DAA-resistant variants in 191 NS3 and 116 NS5B isolates from 208 DAA-naive G1-infected patients. RESULTS: A total of 3.5 million high-quality reads of ≥ 200 nucleotides were generated. The median coverage depth was 4150x and 4470x per NS3 and NS5B amplicon, respectively. Both G1a and G1b populations showed Shannon entropy distributions, with no difference between G1a and G1b in NS3 or NS5B region at the nucleotide level. A higher number of substitutions that confer resistance to protease inhibitors were observed in G1a isolates (mainly at amino acid 80 of the NS3 region). The prevalence of amino acid substitutions that confer resistance to NS5B non-nucleoside inhibitors was similar in G1a and G1b isolates. The NS5B S282T variant, which confers resistance to the polymerase inhibitors mericitabine and sofosbuvir, was not detected in any sample. CONCLUSION: The quasispecies genetic diversity and prevalence of DAA-resistant variants was similar in G1a and G1b isolates and in both NS3 and NS5B regions, suggesting that this is not a determinant for the higher level of DAA resistance observed across G1a HCV infected patients upon treatment.

Comparison of the Mechanisms of Drug Resistance among HIV, Hepatitis B, and Hepatitis C.

Human immunodeficiency virus (HIV), hepatitis B virus (HBV), and hepatitis C virus (HCV) are the most prevalent deadly chronic viral diseases. HIV is treated by small molecule inhibitors. HBV is treated by immunomodulation and small molecule inhibitors. HCV is currently treated primarily by immunomodulation but many small molecules are in clinical development. Although HIV is a retrovirus, HBV is a double-stranded DNA virus, and HCV is a single-stranded RNA virus, antiviral drug resistance complicates the development of drugs and the successful treatment of each of these viruses. Although their replication cycles, therapeutic targets, and evolutionary mechanisms are different, the fundamental approaches to identifying and characterizing HIV, HBV, and HCV drug resistance are similar. This review describes the evolution of HIV, HBV, and HCV within individuals and populations and the genetic mechanisms associated with drug resistance to each of the antiviral drug classes used for their treatment.

Hepatitis B virus reverse transcriptase sequence variant database for sequence analysis and mutation discovery.

Drug resistance resulting from reverse transcriptase (RT) mutations is one of the main obstacles to successful hepatitis B virus (HBV) therapy. Indeed, HBV treatment guidelines recommend HBV genotypic resistance testing for patients receiving nucleos(t)ide RT inhibitors (N(t)RTIs) who develop virological failure. N(t)RTI-resistance mutations at 10 RT positions have been well characterized in phenotypic studies, however, data are lacking on the relative frequency of these mutations in N(t)RTI-treated and untreated individuals. There are also few published data on the extent of amino acid variation at most of the 344 positions of HBV RT and the extent to which this variation is influenced by N(t)RTI treatment. We retrieved 23,871 HBV RT sequences from GenBank and reviewed the published reports of these sequences to ascertain the number of individuals from whom the sequences were obtained, the N(t)RTI treatments of these individuals, and the year and region of virus sampling. We then used these data to populate a relational database we named HBVrtDB. As of July 2010, HBVrtDB contained 6811 sequences from 3869 individuals reported in 281 references. Among these 3869 individuals, 73% were N(t)RTI-naïve and 27% received one or more N(t)RTIs. Among the 10 well-characterized N(t)RTI-resistance mutations, L80I/V, V173L, L180M, A181T, T184S, S202G and M204I/V were significantly associated with treatment with lamivudine, an l-nucleoside analog, and A181S/T/V and N236T were significantly associated with treatment with adefovir, an acyclic nucleoside phosphonate. A similar analysis of ten additional less well-characterized resistance mutations demonstrated a significant association with N(t)RTI treatment for four of the mutations: L82M, S85A, A200V, and Q215S. We also created an interactive program, HBVseq, to enable users to identify mutations in submitted sequences and retrieve the prevalence of these mutations in HBVrtDB according to genotype and N(t)RTI treatment. HBVrtDB and HBVseq are available at http://hivdb.stanford.edu/HBV/releaseNotes/.

Ultra-deep pyrosequencing of hepatitis B virus quasispecies from nucleoside and nucleotide reverse-transcriptase inhibitor (NRTI)-treated patients and NRTI-naive patients.

The dynamics of emerging nucleoside and nucleotide reverse-transcriptase inhibitor (NRTI) resistance in hepatitis B virus (HBV) are not well understood because standard dideoxynucleotide direct polymerase chain reaction (PCR) sequencing assays detect drug-resistance mutations only after they have become dominant. To obtain insight into NRTI resistance, we used a new sequencing technology to characterize the spectrum of low-prevalence NRTI-resistance mutations in HBV obtained from 20 plasma samples from 11 NRTI-treated patients and 17 plasma samples from 17 NRTI-naive patients, by using standard direct PCR sequencing and ultra-deep pyrosequencing (UDPS). UDPS detected drug-resistance mutations that were not detected by PCR in 10 samples from 5 NRTI-treated patients, including the lamivudine-resistance mutation V173L (in 5 samples), the entecavir-resistance mutations T184S (in 2 samples) and S202G (in 1 sample), the adefovir-resistance mutation N236T (in 1 sample), and the lamivudine and adefovir-resistance mutations V173L, L180M, A181T, and M204V (in 1 sample). G-to-A hypermutation mediated by the apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like family of cytidine deaminases was estimated to be present in 0.6% of reverse-transcriptase genes. Genotype A coinfection was detected by UDPS in each of 3 patients in whom genotype G virus was detected by direct PCR sequencing. UDPS detected low-prevalence HBV variants with NRTI-resistance mutations, G-to-A hypermutation, and low-level dual genotype infection with a sensitivity not previously possible.