Complex Genotype Mixtures Analyzed by Deep Sequencing in Two Different Regions of Hepatitis B Virus

This study assesses the presence and outcome of genotype mixtures in the polymerase/surface and X/preCore regions of the HBV genome in patients with chronic hepatitis B virus (HBV) infection. Thirty samples from ten chronic hepatitis B patients were included. The polymerase/surface and X/preCore regions were analyzed by deep sequencing (UDPS) in the first available sample at diagnosis, a pre-treatment sample, and a sample while under treatment. HBV genotype was determined by phylogenesis. Quasispecies complexity was evaluated by mutation frequency and nucleotide diversity. The polymerase/surface and X/preCore regions were validated for genotyping from 113 GenBank reference sequences. UDPS yielded a median of 10,960 sequences per sample (IQR 16,645) in the polymerase/surface region and 11_,595 sequences per sample (IQR 14,682) in X/preCore. Genotype mixtures were more common in X/preCore (90%) than in polymerase/surface (30%) (p<0.001). On X/preCore genotyping, all samples were genotype A, whereas polymerase/surface yielded genotypes A (80%), D (16.7%), and F (3.3%) (p = 0.036). Genotype changes in polymerase/surface were observed in four patients during natural quasispecies dynamics and in two patients during treatment. There were no genotype changes in X/preCore. Quasispecies complexity was higher in X/preCore than in polymerase/surface (p = 0.004).The results provide evidence of genotype mixtures and differential genotype proportions in the polymerase/surface and X/preCore regions. The genotype dynamics in HBV infection and the different patterns of quasispecies complexity in the HBV genome suggest a new paradigm for HBV genotype classification.     

Dynamic changes of hepatitis B virus polymerase gene including YMDD motif in lamivudine-treated patients with chronic hepatitis B

The mutation of YMDD motif of hepatitis B virus (HBV) polymerase gene is the most frequent cause in HBV resistant to lamivudine. The aim of the study was to investigate variation features of HBV polymerase gene in chronic hepatitis B (CHB) patients before and after lamivudine treatment. From the serum samples of five CHB patients before and after 12 months of lamivudine treatment, HBV polymerase gene was amplificated and positive DNA fragments were cloned into JM105 competent cell. Twenty positive clones of every sample were checked with mismatched polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and YMDD variants were sequenced. Among five patients after 12 months of lamivudine treatment, M552I mutations in two patients with HBV DNA rebounding and D553G mutation in one non-responder were detected except two patients with negative HBV DNA consecutively. In summary, D553G mutation is probably one of the reasons that caused non-responders during lamivudine treatment. The mutations of YMDD motif occurred after lamivudine treatment are caused by the induction of lamivudine.     

Detection of mixed populations of wild-type and YMDD hepatitis B variants by pyrosequencing in acutely and chronically infected patients

ABSTRACT: BACKGROUND: Lamivudine (LAM) is associated with the highest known rate of resistance mutations amongnucleotide analogs used to treat chronic hepatitis B virus (HBV) infection. Despite this, LAMcontinues in widespread use, especially in combination therapies. The primary LAMresistance mutation (rtM204V/I) occurs in the YMDD motif of HBV polymerase. The aim ofthis study was to characterize Brazilian HBV isolates from acute and chronic cases by directsequencing, and to identify HBV quasispecies in the YMDD motif using a pyrosequencingmethod capable of detecting single-nucleotide polymorphisms. HBV DNA from serumsamples of 20 individuals with acute HBV infection and 44 with chronic infectionundergoing antiviral therapies containing LAM were analyzed by direct sequencing andpyrosequencing methods. RESULTS: Phylogenic analyses of direct-sequenced isolates showed the expected genotypes (A, D andF) for the Brazilian population in both acute and chronic infections. However, withingenotype A isolates, subgenotype A2 was more frequently detected in acute cases than inchronic cases (P = 0.012). As expected, none of the individuals with acute hepatitis B hadLAM-resistant isolates as a dominant virus population, whether detected by direct sequencingor pyrosequencing. However, pyrosequencing analyses showed that 45% of isolates (9/20)had minor subpopulations (4-17%) of LAM-resistant isolates. Among chronic patientsundergoing LAM treatment, YMDD mutants were frequently found as a dominant viruspopulation. In cases where wild-type virus was the dominant population, subpopulations ofYMDD variants were usually found, demonstrating the complexity of HBV quasispecies. CONCLUSIONS: YMDD variants were frequently detected as a minor population in acute HBV infection. Theoccurrence of pre-existing variants may lead to a high frequency of resistant mutants duringantiviral therapy in the chronic phase. In chronic infection, detection of YMDD variantsbefore virological or biochemical breakthrough might contribute to making better therapychoices and thus improving treatment outcome.     

Ultra-Deep Pyrosequencing (UDPS) Data Treatment to Study Amplicon HCV Minor Variants

We have investigated the reliability and reproducibility of HCV viral quasispecies quantification by ultra-deep pyrosequencing (UDPS) methods. Our study has been divided in two parts. First of all, by UDPS sequencing of clone mixes samples we have established the global noise level of UDPS and fine tuned a data treatment workflow previously optimized for HBV sequence analysis. Secondly, we have studied the reproducibility of the methodology by comparing 5 amplicons from two patient samples on three massive sequencing platforms (FLX+, FLX and Junior) after applying the error filters developed from the clonal/control study. After noise filtering the UDPS results, the three replicates showed the same 12 polymorphic sites above 0.7%, with a mean CV of 4.86%. Two polymorphic sites below 0.6% were identified by two replicates and one replicate respectively. A total of 25, 23 and 26 haplotypes were detected by GS-Junior, GS-FLX and GS-FLX+. The observed CVs for the normalized Shannon entropy (Sn), the mutation frequency (Mf), and the nucleotidic diversity (Pi) were 1.46%, 3.96% and 3.78%. The mean absolute difference in the two patients (5 amplicons each), in the GS-FLX and GS-FLX+, were 1.46%, 3.96% and 3.78% for Sn, Mf and Pi. No false polymorphic site was observed above 0.5%.Our results indicate that UDPS is an optimal alternative to molecular cloning for quantitative study of HCV viral quasispecies populations, both in complexity and composition. We propose an UDPS data treatment workflow for amplicons from the RNA viral quasispecies which, at a sequencing depth of at least 10,000 reads per strand, enables to obtain sequences and frequencies of consensus haplotypes above 0.5% abundance with no erroneous mutations, with high confidence, resistant mutants as minor variants at the level of 1%, with high confidence that variants are not missed, and highly confident measures of quasispecies complexity.     

YVDD mutation of hepatitis B Virus, a dominant Lamivudine-resistant type in Guangzhou, South China

The epidemiological effects of native and mutated YMDD motif in the HBV genome under the selective pressure of lamivudine were investigated. YMDD wild and mutation motif in HBV genome were detected by flow through reverse dot blots (FT-RDB) with KaiPuTM DNA HybriMax Rapid Hybridization Machine based on the principle of "Flow-through hybridization" and by the traditional Reverse Dot Blot assay. Sera from 1021 suspected lamivudine-resistant chronic HBV carders after more than 8 months of lamivudine therapy and the corresponding archived sera were collected and assayed. We found 35.94% were single type infections with 8.03% YMDD, 7.93% YIDD and 19.98% YVDD. It was also found that 64.06% were mixed infections including 1.96% YMDD and YIDD, 51.62% YMDD and YVDD, 1.96% YIDD and YVDD, 8.52% YMDD, YIDD and YVDD. The levels of infections containing YVDD motif reached 82.08%. The pretreatment infectious status were: YMDD single infection was 36.93%; YIDD single infection was 6.07%; YVDD single infection was 17.04%; YMDD and YIDD mixed infection was 0.97%; YMDD and YVDD mixed infection was33.99%; YIDD and YVDD mixed infection was 0.98%; YMDD, YIDD and YVDD mixed infection was 4.02%. Infections containing YVDD motif were only 56.03%. The 34.32% mutation rate of YMDD motif to YVDD was significantly higher than the 10.97% of YMDD to YIDD (U=10.98, P<0.05), as estimated by Mann-Whitney U-test for non-parametric data. HBV containing YVDD motif might have an evolutionary ascendancy and become the dominant type under the selective pressure of lamivudine.     

The naturally occurring YMDD mutation among patients chronically infected HBV and untreated with lamivudine: a systematic review and meta-analysis

Background

Several recent reports have demonstrated that tyrosine (Y)-methionine (M)-aspartic acid (D)-aspartic acid (D) (YMDD) motif mutations can naturally occur in chronic HBV patients without antiviral treatment such as lamivudine therapy. This paper aims to assess the overall spontaneous incidence and related risk factors of YMDD-motif mutations among lamivudine-naïve chronic HBV carriers, so as to provide some clue for clinical treatment of hepatitis B.

Methodology/Principal Findings

Chinese and English literatures were searched for studies reporting natural YMDD mutations among untreated chronic HBV patients from 2001 to 2010. The incidence estimates were summarized and analyzed by meta-analyses. Forty-seven eligible articles from eight countries were selected in this review (13 in English and 34 in Chinese). The pooled incidence of YMDD-motif mutation among untreated chronic HBV patients from eight countries was 12.21% (95% CI: 9.69%–14.95%). China had an incidence of 13.38% (95% CI: 10.90%–16.07%) and seven other countries had an incidence of 9.90% (95% CI: 3.28%–19.55%), respectively. Lamivudine therapy would increase the risk of mutations 5.23 times higher than the untreated patients. A higher HBV DNA copy number was associated with increased incidence of natural YMDD mutation. No significant difference was found in YMDD mutation incidence between groups of different gender, age, HBeAg status, patients'' ALT (alanine aminotransferase) level, and between the groups of HBV genotype B and C.

Conclusions

The YMDD-motif mutations can occur spontaneously with a relatively high incidence in CHB patients untreated with lamivudine. These mutations might be the consequence of accumulated base mismatch due to the nature of viral polymerase. More fundamental and clinical studies are needed to clarify the influence of YMDD mutations in hepatitis B progression and antiviral treatment.     

未经抗病毒治疗慢乙肝患者YMDD变异检测及HBV基因分型

目的 :探讨未接受拉米呋啶及干扰素抗病毒治疗的慢性乙型肝炎患者中YMDD是否存在变异及HBV的基因亚型。方法 :采用PCR荧光分子信标探针检测技术及PCR微板核酸杂交 ELISA技术进行HBVYMDD变异株 (YIDD YVDD)及HBV基因分型检测。结果 :2 9例慢性乙肝患者HBVYMDD变异阳性 6例 ,阳性率为 2 0 6%。 5例变异阳性中 3例为YIDD阳性 ,2例为YVDD阳性 ,6例均为YMDD野生株和变异株同时存在 ;HBV基因亚型以C、B、D亚型为主 ,分别占41 3% ( 1 2 2 9)、2 4 1 % ( 7 2 9)、2 0 6% ( 6 2 9)。结论 :在未经抗病毒治疗的慢乙肝患者中存在YMDD变异株 ,临床上在应用拉米呋啶对慢性乙型肝炎患者进行抗病毒治疗前 ,应对患者是否存在YMDD变异加以重视 。     

Clinical Application of Estimating Hepatitis B Virus Quasispecies Complexity by Massive Sequencing: Correlation between Natural Evolution and On-Treatment Evolution

Aim

To evaluate HBV quasispecies (QA) complexity in the preCore/Core regions in relation to HBeAg status, and explore QA changes under natural evolution and nucleoside analogue (NUC) treatment.

Methods

Ultra-deep pyrosequencing of HBV preCore/Core regions in 30 sequential samples (baseline [diagnosis], treatment-free, and treatment-nonresponse) from 10 retrospectively selected patients grouped according to HBeAg status over time: HBeAg+ (N = 4), HBeAg- (N = 2), and fluctuating HBeAg (transient seroreversion/seroconversion pattern) (N = 4). QA complexity was defined by Shannon entropy, mutation frequency, nucleotide diversity, and mutation frequency of amino acids (MfAA) in preCore and Core.

Results

The QA was less complex in HBeAg+ than in HBeAg- or fluctuating HBeAg. High complexity in preCore was associated with decreased viral replication (preCore MfAA negatively correlated with HBV-DNA, p = 0.005). QA complexity in the treatment-free period negatively correlated with values seen during treatment. Specific variants were mainly selected in the Core region in HBeAg- and fluctuating HBeAg patients, suggesting higher immune pressure than in HBeAg+.

Conclusions

The negative correlation between QA natural evolution and on-treatment evolution indicates the importance of pre-treatment QA study to predict QA changes in NUC nonresponders. Study of QA complexity could be useful for managing HBV infection.     

Mismatch repair-independent increase in spontaneous mutagenesis in yeast lacking non-essential subunits of DNA polymerase ε

Yeast DNA polymerase ε (Pol ε) is a highly accurate and processive enzyme that participates in nuclear DNA replication of the leading strand template. In addition to a large subunit (Pol2) harboring the polymerase and proofreading exonuclease active sites, Pol ε also has one essential subunit (Dpb2) and two smaller, non-essential subunits (Dpb3 and Dpb4) whose functions are not fully understood. To probe the functions of Dpb3 and Dpb4, here we investigate the consequences of their absence on the biochemical properties of Pol ε in vitro and on genome stability in vivo. The fidelity of DNA synthesis in vitro by purified Pol2/Dpb2, i.e. lacking Dpb3 and Dpb4, is comparable to the four-subunit Pol ε holoenzyme. Nonetheless, deletion of DPB3 and DPB4 elevates spontaneous frameshift and base substitution rates in vivo, to the same extent as the loss of Pol ε proofreading activity in a pol2-4 strain. In contrast to pol2-4, however, the dpb3Δdpb4Δ does not lead to a synergistic increase of mutation rates with defects in DNA mismatch repair. The increased mutation rate in dpb3Δdpb4Δ strains is partly dependent on REV3, as well as the proofreading capacity of Pol δ. Finally, biochemical studies demonstrate that the absence of Dpb3 and Dpb4 destabilizes the interaction between Pol ε and the template DNA during processive DNA synthesis and during processive 3' to 5'exonucleolytic degradation of DNA. Collectively, these data suggest a model wherein Dpb3 and Dpb4 do not directly influence replication fidelity per se, but rather contribute to normal replication fork progression. In their absence, a defective replisome may more frequently leave gaps on the leading strand that are eventually filled by Pol ζ or Pol δ, in a post-replication process that generates errors not corrected by the DNA mismatch repair system.     

Multilayered defense in HLA-B51-associated HIV viral control

Polymorphism in the HLA region of a chromosome is the major source of host genetic variability in HIV-1 outcome, but there is limited understanding of the mechanisms underlying the beneficial effect of protective class I alleles such as HLA-B57, -B27, and -B51. Taking advantage of a unique cohort infected with clade B' HIV-1 through contaminated blood, in which many variables such as the length of infection, the infecting viral strain, and host genetic background are controlled, we performed a comprehensive study to understand HLA-B51-associated HIV-1 control. We focused on the T cell responses against three dominant HLA-B51-restricted epitopes: Gag327-345(NI9) NANPDCKTI, Pol743-751(LI9) LPPVVAKEI, and Pol283-289(TI8) TAFTIPSI. Mutations in all three dominant epitopes were significantly associated with HLA-B51 in the cohort. A clear hierarchy in selection of epitope mutations was observed through epitope sequencing. L743I in position 1 of epitope LI9 was seen in most B51(+) individuals, followed by V289X in position 8 of the TI8, and then, A328S, in position 2 of the NI9 epitope, was also seen in some B51(+) individuals. Good control of viral load and higher CD4(+) counts were significantly associated with at least one detectable T cell response to unmutated epitopes, whereas lower CD4(+) counts and higher viral loads were observed in patients who had developed escape mutations in all three epitopes or who lacked T cell responses specific to these epitope(s). We propose that patients with HLA-B51 benefit from having multiple layers of effective defense against the development of immune escape mutations.     

The Glycine-Rich Motif of Pyrococcus abyssi DNA Polymerase D Is Critical for Protein Stability

A glycine-rich motif described as being involved in human polymerase δ proliferating cell nuclear antigen (PCNA) binding has also been identified in all euryarchaeal DNA polymerase D (Pol D) family members. We redefined the motif as the (G)-PYF box. In the present study, Pol D (G)-PYF box motif mutants from Pyrococcus abyssi were generated to investigate its role in functional interactions with the cognate PCNA. We demonstrated that this motif is not essential for interactions between PabPol D (P. abyssi Pol D) and PCNA, using surface plasmon resonance and primer extension studies. Interestingly, the (G)-PYF box is located in a hydrophobic region close to the active site. The (G)-PYF box mutants exhibited altered DNA binding properties. In addition, the thermal stability of all mutants was reduced compared to that of wild type, and this effect could be attributed to increased exposure of the hydrophobic region. These studies suggest that the (G)-PYF box motif mediates intersubunit interactions and that it may be crucial for the thermostability of PabPol D.     

Association of hepatitis B virus polymerase with promyelocytic leukemia nuclear bodies mediated by the S100 family protein p11

Hepatitis B virus (HBV) polymerase (Pol) interacts with cellular chaperone proteins and thereby performs multiple functions necessary for viral replication. Yeast two-hybrid analysis was applied to identify additional cellular targets required for HBV Pol function. HBV Pol interacted with S100A10 (p11), a Ca(2+)-modulated protein previously shown to bind to annexin II. The interaction between HBV Pol and p11 was confirmed by co-immunoprecipitation of the two proteins synthesized either in vitro or in transfected cells and by inhibition of the DNA polymerase activity of HBV Pol by p11. Immunofluorescence analysis of transfected human cell lines revealed that, although most HBV Pol and p11 was restricted to the cytoplasm, a small proportion of each protein colocalized as nuclear speckles; HBV Pol was not detected in the nucleus in the absence of p11. The HBV Pol-p11 nuclear speckles coincided with nuclear bodies containing the promyelocytic leukemia protein PML. Furthermore, the association of HBV Pol-p11 with PML was increased by exposure of cells to EGTA and inhibited by valinomycin. These results suggest a role for p11 in modulation of HBV Pol function and implicate PML nuclear bodies and intracellular Ca(2+) in viral replication.