HBV DNA was quantified by real-time polymerase chain reaction (PCR) using the TaqMan Fast Advanced Master Mix reagent (Thermo Fisher Scientific, MA, USA) on Applied Biosystems StepOne (Thermo Fisher Scientific, MA, USA).
Using PrimeScript One-Step RT-PCR Kit Ver.2, HBV DNA was amplified by nested PCR targeting the overlapping surface-polymerase (SP) or the surface region of the viral genome in case the SP region amplification failed. PCR products were then sequenced by the Sanger sequencing method on SeqStudio Sequence Analyzer (Thermo Fisher Scientific, MA, USA) using the BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific, CA, USA). Then, HBV genotypes were identified by analyzing partial genome sequences with 103 reference strains retrieved from GenBank, using the neighbor-joining method with the Molecular Evolutionary Genetics Analysis software version 10 (Pennsylvania State University, PA, USA). The resulting phylogenetic trees are shown in Supplementary Figs.