AimsThe dynamics of resistance-associated mutations under combination therapy were explored.MethodsA total of 46 patients were classified into adefovir (n=14) and entecavir (n=32) groups. In the adefovir (ADV) group, six patients receiving combined therapy were DNA-positive after more than 3 years of therapy. Ultra-deep pyrosequencing was used to analyze the dynamics of multi-drugs resistance mutations.ResultsAt baseline, all 46 treatment-naïve patients harbored rtA181V/T substitutions (1.2%-4.6%) and rtN236T substitutions (1.6%-6.1%). In the ADV group, eight patients with long-term treatment were consecutively HBV DNA-positive for more than 3 years. During treatment, the rtA181T resistance-associated site appeared with increasing frequency in six of eight patients (NOs. 1-6), and two patients (NOs.4 and 8) carrying the rtA181T resistance mutations increasingly showed high levels of rtN236T. One patient (NO. 8) experienced virological breakthrough. Other known pre-existing mutations showed no dynamic fluctuations, including in rtA194T, rtP177G, rtF249A, and rtD263E. In addition to the common substitutions, some previously unknown amino acid substitutions, such as rtD134N, rtL145M/S, rtF151Y/L, rtR153Q, and rtS223A, should be further studied.ConclusionsHBV-resistance substitutions conferring to nucleoside analogs are present at baseline. The dynamics of the HBV RT-region quasispecies variation are heterogeneous and complex.