Hepatitis B virus drug resistance tools: one sequence, two predictions.

Drug resistance testing, genotype analysis, and the determination of immune and vaccine escape variants support personalized antiviral treatment for hepatitis B patients. As phenotypic drug resistance testing for hepatitis B virus (HBV) is especially labor-intensive, due to the lack of simple cell culture systems, genotypic methods play a very pronounced role. The genetic structure of HBV allows the simultaneous analysis of two different genes by examination of a single region in the genome of HBV. Nevertheless, the overlapping open reading frames of the hepatitis B surface antigen (HBsAg) and the reverse transcriptase (RT) have to be interpreted separately. In diagnostic procedures, standard Sanger type sequencing (mostly performed as a dye-dideoxynucleotide terminator system) is still the most commonly used method. This allows using established techniques for interpreting those types of genetic information. Besides reviewing the foundation of drug resistance interpretation for HBV, different interpretation systems, either commercially available (TRUGENE, Abbott, ViroScore) or free to use (geno2pheno[HBV], HIV-GRADE HBV tool), are compared in this article.

High frequency of complex mutational patterns in lamivudine resistant hepatitis B virus isolates.

The evolution of hepatitis-B virus (HBV) drug resistance is characterized by the emergence of resistance conferring mutations and compensatory mutations. Therefore HBV drug-resistant isolates often harbor multiple mutations in the reverse transcriptase (RT) and corresponding mutations in the hepatitis B surface antigen (HBsAg). In this study mutational patterns of 60 HBV isolates harboring drug resistance mutations rtM204V or rtM204I were retrospectively analyzed. Both mutations, especially mutation rtM204V, were most often accompanied by compensatory mutations rtV173L and rtL180M but also by mutations conferring entecavir (n = 5) or adefovir resistance (n = 4). In addition, 22 (36.7%) drug resistant HBV isolates carried mutations related to immune escape in the HBsAg. In seven cases premature stop codons in HBsAg were detected resulting in the expression of truncated HBsAg. Clonal analysis of these seven quasispecies even disclosed the presence of HBV isolates carrying further stop codons and in one case the occurrence of resistance mutation rtA181T without the stop codon mutation sW172*. Interestingly, only one HBV clone carried the resistance mutations rtM204V and rtA181T. HBV drug resistant isolates frequently harbored HBsAg mutations associated with immune escape or disease progression pointing to a complex interaction of both proteins. HBV genotypic resistance tests based on population sequencing methods seemed to be inappropriate for determining the clinical relevance of stop codons in the HBsAg.

Genotyping hepatitis B virus dual infections using population-based sequence data.

The hepatitis B virus (HBV) is classified into distinct genotypes A-H that are characterized by different progression of hepatitis B and sensitivity to interferon treatment. Previous computational genotyping methods are not robust enough regarding HBV dual infections with different genotypes. The correct classification of HBV sequences into the present genotypes is impaired due to multiple ambiguous sequence positions. We present a computational model that is able to identify and genotype inter- and intragenotype dual infections using population-based sequencing data. Model verification on synthetic data showed 100 % accuracy for intergenotype dual infections and 36.4 % sensitivity in intragenotype dual infections. Screening patient sera (n = 241) revealed eight putative cases of intergenotype dual infection (one A-D, six A-G and one D-G) and four putative cases of intragenotype dual infection (one A-A, two D-D and one E-E). Clonal experiments from the original patient material confirmed three out of three of our predictions. The method has been integrated into geno2pheno([hbv]), an established web-service in clinical use for analysing HBV sequence data. It offers exact and detailed identification of HBV genotypes in patients with dual infections that helps to optimize antiviral therapy regimens. geno2pheno([hbv]) is available under http://www.genafor.org/g2p_hbv/index.php.

Results of the first international HIV-1 coreceptor proficiency panel test.

BACKGROUND: Quality Assurance (QA) programs are essential to evaluate performance in diagnostics laboratories. OBJECTIVES: We present the results from the first QA for HIV-1 genotypic tropism testing, organized and coordinated by the Institute of Virology at the University of Cologne. STUDY DESIGN: 12 cell culture-derived viral strains of different HIV-1 clades from the NIH AIDS Reagent Program were sent to the participants to be processed with their standard diagnostic methods Fasta files containing the V3 region sequence were centrally analyzed at the Institute of Virology, Cologne. All samples were sent in parallel for phenotypic testing. RESULTS: 36 laboratories from 16 countries reported genotypic results. The sequence-generation efficacy was 95.1%, while the phenotypic assays ESTA® and PhenXR only achieved results for 58.3% of the samples. All four X4 samples were identified by 31/36 laboratories, two laboratories amplified 3/4×4 samples, and three detected 2/4×4 samples. There was high concordance among the genotypic and phenotypic results, although differences in FPR values were detected. Most deficiencies in sequence editing did not result in wrong classification of X4 viruses as R5, with the exception of sample NRZ05 by laboratory 38, but in an overestimation of CXCR4 use. CONCLUSIONS: This demonstrates that genotypic tropism prediction is a safe procedure for clinical purposes. As we used homogeneous cell culture samples and all sequence fasta files were centrally analyzed, variations in FPR values can only be attributed to sample preparation, RT-PCR or sequence editing protocols.

Geno2pheno[HCV] - A Web-based Interpretation System to Support Hepatitis C Treatment Decisions in the Era of Direct-Acting Antiviral Agents.

The face of hepatitis C virus (HCV) therapy is changing dramatically. Direct-acting antiviral agents (DAAs) specifically targeting HCV proteins have been developed and entered clinical practice in 2011. However, despite high sustained viral response (SVR) rates of more than 90%, a fraction of patients do not eliminate the virus and in these cases treatment failure has been associated with the selection of drug resistance mutations (RAMs). RAMs may be prevalent prior to the start of treatment, or can be selected under therapy, and furthermore they can persist after cessation of treatment. Additionally, certain DAAs have been approved only for distinct HCV genotypes and may even have subtype specificity. Thus, sequence analysis before start of therapy is instrumental for managing DAA-based treatment strategies. We have created the interpretation system geno2pheno[HCV] (g2p[HCV]) to analyse HCV sequence data with respect to viral subtype and to predict drug resistance. Extensive reviewing and weighting of literature related to HCV drug resistance was performed to create a comprehensive list of drug resistance rules for inhibitors of the HCV protease in non-structural protein 3 (NS3-protease: Boceprevir, Paritaprevir, Simeprevir, Asunaprevir, Grazoprevir and Telaprevir), the NS5A replicase factor (Daclatasvir, Ledipasvir, Elbasvir and Ombitasvir), and the NS5B RNA-dependent RNA polymerase (Dasabuvir and Sofosbuvir). Upon submission of up to eight sequences, g2p[HCV] aligns the input sequences, identifies the genomic region(s), predicts the HCV geno- and subtypes, and generates for each DAA a drug resistance prediction report. g2p[HCV] offers easy-to-use and fast subtype and resistance analysis of HCV sequences, is continuously updated and freely accessible under http://hcv.geno2pheno.org/index.php. The system was partially validated with respect to the NS3-protease inhibitors Boceprevir, Telaprevir and Simeprevir by using data generated with recombinant, phenotypic cell culture assays obtained from patients' virus variants.

Entecavir allows an unexpectedly high residual replication of HBV mutants resistant to lamivudine.

BACKGROUND: Entecavir is an efficient inhibitor of HBV reverse transcriptase (RT) and widely used for therapy of chronic hepatitis B. Entecavir treatment of HBV patients with lamivudine-resistant viral strains, however, often fails, but the mechanism of cross-resistance development is not fully understood. METHODS: Using non-linear regression models, dose-response curves of cloned HBV strains from patients pre-treated with RT inhibitors were established in human hepatoma cell lines after transfection with HBV genomes containing HBV polymerase genes from patient isolates. 50% and 90% inhibitory concentrations (IC50 and IC90) and corresponding antiviral resistance factors (RF50 and RF90) were calculated. RESULTS: The entecavir dose-response curve of lamivudine-resistant HBV RT mutants rtM204 for the replication of HBV decreased less than expected with increasing drug dose. Remarkably, due to the flat dose-response curves, RF90 values against entecavir of samples with rtM204 substitutions were up to 30× higher than their RF50 values. CONCLUSIONS: The unexpectedly high IC90 indicates a strong residual replication capacity of lamivudine-resistant HBV rtM204 variants under entecavir therapy, although IC50 values are initially within the therapeutic range of entecavir. This characteristic favours rapid selection of additional mutants with overt resistance against entecavir. Thus, the current phenotypic resistance assays should include determination of IC90.

Clinical performance of the novel DiaSorin LIAISON(®) XL murex: HBsAg Quant, HCV-Ab, HIV-Ab/Ag assays.

BACKGROUND: The fully automated and closed LIAISON(®)XL platform was developed for reliable detection of infection markers like hepatitis B virus (HBV) surface antigen (HBsAg), hepatitis C virus (HCV) antibodies (Ab) or human immunodeficiency virus (HIV)-Ag/Ab. To date, less is known about the diagnostic performance of this system in direct comparison to the common Abbott ARCHITECT(®) platform. OBJECTIVES: We compared the diagnostic performance and usability of the DiaSorin LIAISON(®)XL with the commonly used Abbott ARCHITECT(®) system. STUDY DESIGN: The qualitative performance of the above mentioned assays was compared in about 500 sera. Quantitative tests were performed for HBsAg-positive samples from patients under therapy (n=289) and in vitro expressed mutants (n=37). For HCV-Ab, a total number of 155 selected samples from patients chronically infected with different HCV genotypes were tested. RESULTS: The concordance between both systems was 99.4% for HBsAg, 98.81% for HCV-Ab, and 99.6% for HIV-Ab/Ag. The quantitative LIAISON(®)XL murex HBsAg assay detected all mutants in comparable amounts to the HBsAg wild type and yielded highly reliable HBsAg kinetics in patients treated with antiviral drugs. Dilution experiments using the 2nd International Standard for HBsAg (WHO) showed a high accuracy of this test. HCV-Ab from patients infected with genotypes 1-3 were equally detected in both systems. Interestingly, S/CO levels of HCV-Ab from patients infected with genotype 3 seem to be relatively low using both systems. CONCLUSIONS: The LIAISON(®)XL platform proved to be an excellent system for diagnostics of HBV, HCV, and HIV with equal performance compared to the ARCHITECT(®) system.

Inferring short-range linkage information from sequencing chromatograms.

Direct Sanger sequencing of viral genome populations yields multiple ambiguous sequence positions. It is not straightforward to derive linkage information from sequencing chromatograms, which in turn hampers the correct interpretation of the sequence data. We present a method for determining the variants existing in a viral quasispecies in the case of two nearby ambiguous sequence positions by exploiting the effect of sequence context-dependent incorporation of dideoxynucleotides. The computational model was trained on data from sequencing chromatograms of clonal variants and was evaluated on two test sets of in vitro mixtures. The approach achieved high accuracies in identifying the mixture components of 97.4% on a test set in which the positions to be analyzed are only one base apart from each other, and of 84.5% on a test set in which the ambiguous positions are separated by three bases. In silico experiments suggest two major limitations of our approach in terms of accuracy. First, due to a basic limitation of Sanger sequencing, it is not possible to reliably detect minor variants with a relative frequency of no more than 10%. Second, the model cannot distinguish between mixtures of two or four clonal variants, if one of two sets of linear constraints is fulfilled. Furthermore, the approach requires repetitive sequencing of all variants that might be present in the mixture to be analyzed. Nevertheless, the effectiveness of our method on the two in vitro test sets shows that short-range linkage information of two ambiguous sequence positions can be inferred from Sanger sequencing chromatograms without any further assumptions on the mixture composition. Additionally, our model provides new insights into the established and widely used Sanger sequencing technology. The source code of our method is made available at http://bioinf.mpi-inf.mpg.de/publications/beggel/linkageinformation.zip.

Anti-HBV treatment induces novel reverse transcriptase mutations with reflective effect on HBV S antigen.

INTRODUCTION: The identification of novel reverse-transcriptase (RT) drug-resistance mutations is critical in predicting the probability of success to anti-HBV treatment. Furthermore, due to HBV-RT/HBsAg gene-overlap, they can have an impact on HBsAg-detection and quantification. METHODS: 356 full-length HBV-RT sequences from 197 drug-naive patients and 159 patients experiencing virological-breakthrough to nucleoside/nucleotide-analogs (NUCs) were analyzed. Mutants and wild-type HBs-antigens were expressed in HuH7-hepatocytes and quantified in cell-supernatants and cell-lysates by Architect HBsAg-assay. RESULTS: Ten novel RT-mutations (rtN53T-rtS78T-rtS85F-rtS135T-rtA181I-rtA200V-rtK212Q-rtL229V/F-rtM309K) correlated with specific NUC-treatments and classical drug-resistance mutations on divergent evolutionary pathways. Some of them reduced RT-binding affinity for anti-HBV drugs and altered S-antigen structure. Indeed, rtS78T (prevalence: 1.1% in drug-naïve and 12.2% in adefovir-failing patients) decreased the RT-affinity for adefovir more than the classical adefovir-resistance mutations rtA181 T/V (WT:-9.63 kcal/mol, rtA181T:-9.30 kcal/mol, rtA181V:-7.96 kcal/mol, rtS78T:-7.37 kcal/mol). Moreover, rtS78T introduced a stop-codon at HBsAg-position 69, and completely abrogated HBsAg-quantification in both supernatants and cell-lysates, indicating an impaired HBsAg-secretion/production. Furthermore, the HBsAg-mutation sP217L, silent in RT, significantly correlated with M204V/I-related virological-breakthrough and increased HBsAg-quantification in cell-lysate. CONCLUSIONS: Mutations beyond those classically known can affect drug-binding affinity of mutated HBV-RT, and may have potential effects on HBsAg. Their cumulative effect on resistance and HBV-pathogenicity indicates the importance of preventing therapeutic failures.

The detection of HBsAg mutants expressed in vitro using two different quantitative HBsAg assays.

BACKGROUND: The quantification of hepatitis-B surface antigen (HBsAg) is useful to identify inactive carriers in chronically HBV infected patients and to predict interferon treatment outcome. OBJECTIVE: To compare two quantitative HBsAg assays for the detection of HBsAg mutants expressed in vitro. STUDY DESIGN: HBsAg mutants (n=35) were expressed in HuH7 cells and the supernatants were tested in two different quantitative HBsAg assays (Architect, Abbott and Elecsys, Roche). RESULTS: All HBsAg mutants were detected by both assays, but in general the results of the Architect system were higher than those of the Elecsys system. The detection of HBsAg mutants in comparison to wild type was similar using both assays. However, HBsAg mutation T123A was under quantified by the Architect, whereas HBsAg mutations P142L, P142S and G145K yielded lower results in the Elecsys system. CONCLUSIONS: Even though HBsAg assays are optimised for the detection of HBsAg mutants, discrepant results were obtained for some HBsAg mutants in two quantitative HBsAg assays. These findings have to be considered when testing samples from one patient with two different quantitative HBsAg assays.

Novel HBsAg markers tightly correlate with occult HBV infection and strongly affect HBsAg detection.

Occult HBV infection (OBI) is a threat for the safety of blood-supply, and has been associated with the onset of HBV-related hepatocellular carcinoma and lymphomagenesis. Nevertheless, genetic markers in HBsAg (particularly in D-genotype, the most common in Europe) significantly associated with OBI in vivo are missing. Thus, the goal of this study is to define: (i) prevalence and clinical profile of OBI among blood-donors; (ii) HBsAg-mutations associated with OBI; (iii) their impact on HBsAg-detection. OBI was searched among 422,278 blood-donors screened by Nucleic-Acid-Testing. Following Taormina-OBI-definition, 26 (0.006%) OBI-patients were identified. Despite viremia <50IU/ml, HBsAg-sequences were obtained for 25/26 patients (24/25 genotype-D). OBI-associated mutations were identified by comparing OBI-HBsAg with that of 82 chronically-infected (genotype-D) patients as control. Twenty HBsAg-mutations significantly correlated for the first time with OBI. By structural analysis, they localized in the major HBV B-cell-epitope, and in HBsAg-capsid interaction region. 14/24 OBI-patients (58.8%) carried in median 3 such mutations (IQR:2.0-6.0) against 0 in chronically-infected patients. By co-variation analysis, correlations were observed for R122P+S167L (phi=0.68, P=0.01), T116N+S143L (phi=0.53, P=0.03), and Y100S+S143L (phi=0.67, p<0.001). Mutants (obtained by site-directed mutagenesis) carrying T116N, T116N+S143L, R122P, R122P+Q101R, or R122P+S167L strongly decreased HBsAg-reactivity (54.9±22.6S/CO, 31.2±12.0S/CO, 6.1±2.4S/CO, 3.0±1.0S/CO and 3.9±1.3S/CO, respectively) compared to wild-type (306.8±64.1S/CO). Even more, Y100S and Y100S+S143L supernatants show no detectable-HBsAg (experiments in quadruplicate). In conclusions, unique HBsAg-mutations in genotype-D, different than those described in genotypes B/C (rarely found in western countries), tightly correlate with OBI, and strongly affect HBsAg-detection. By altering HBV-antigenicity and/or viral-particle maturation, they may affect full-reliability of universal diagnostic-assays for HBsAg-detection.