Evaluation of a Next Generation Sequencing Assay for Hepatitis B Antiviral Drug Resistance on the Oxford Nanopore System ======================================================================================================================== * Michael Payne * Gordon Ritchie * Tanya Lawson * Matthew Young * Willson Jang * Aleksandra Stefanovic * Marc G. Romney * Nancy Matic * Christopher F. Lowe ## Abstract **Background** Next-generation sequencing (NGS) for Hepatitis B virus (HBV) antiviral resistance (AVR) testing is a highly sensitive diagnostic method, able to detect low-level mutant subpopulations. Our clinical virology laboratory previously transitioned from DNA hybridization (INNO-LiPA) to NGS, initially with the GS Junior System and subsequently the MiSeq. The Oxford Nanopore Technology (ONT) sequencing system was evaluated for HBV resistance testing, with regards to sequencing accuracy and turn-around time. **Methods** We performed amplicon sequencing of the HBV polymerase gene from patient plasma and external quality assessment (EQA) samples on the MiSeq Reagent Nano Kit v2 and GridION ONT with R10.4.1 flowcells. Mutational analysis and genotyping were performed by DeepChekAssay-HBV (version 2.0). **Results** A total of 49 patient samples and 15 EQA samples were tested on both the MiSeq and ONT. There was high agreement for both patient and EQA samples between the MiSeq and ONT systems, with regards to total drug resistance mutations detected and total patient sample agreement, 68/70 (97%) and 47/49 (96%), respectively. **Conclusion** The ONT NGS platform provided accurate HBV AVR results, with improved turn- around times. Sequencing error rates at AVR codons were below 1%. Key Words * Hepatitis B virus * Viral resistance testing * Next generation sequencing * Molecular methods ## 1. Background Hepatitis B virus (HBV) is an important cause of chronic hepatitis and hepatocellular carcinoma (HCC) worldwide. Globally, there are an estimated 296 million people living with HBV, along 1.1. million deaths annually (1). Treatment of chronic HBV can decrease the risk of progression to cirrhosis and HCC. Identification of the HBV genotype, along with antiviral resistance (AVR) mutations, can help inform the management of chronic HBV infection (2–4). Newer antiviral agents (Entecavir and Tenofovir) have a higher barrier to AVR mutations and viral breakthrough; however, AVR testing is useful to monitor for emergence of mutations associated with decreased virologic response, particularly for patients with previous exposure to HBV antiviral therapies (3, 5). HBV AVR line probe (LiPA) testing is limited by its inability to detect novel mutations, hybridization failures, and decreased sensitivity with low viral loads (<1000 IU/mL) (6, 7). In 2014, testing in our clinical laboratory was transitioned to next-generation sequencing (NGS) with the GS Junior System (454 Life Sciences, Branford, CT) and in 2017, the MiSeq (Illumina) platform. MiSeq offered increased sequencing depth; however, a shorter read length did not include all AVR mutation sites (8). The Oxford Nanopore Technology (ONT) sequencing system (GridION ONT) with longer read capability and updated R10.4.1 flowcells has reported improved error rates, allowing for more accurate detection of variant subpopulations (9). In addition, this technology has less hands-on library preparation and sequencing run time, allowing for improved workflow in clinical laboratories. Our study aim was to optimize HBV AVR testing using NGS on the ONT platform, along with validation of this test for clinical use. ## 2. Materials and Methods We used guidance from molecular testing references and draft recommendations from the US FDA to design this NGS validation study (10–12). Parameters evaluated included: Limit of detection, inclusivity/accuracy, reproducibility/precision and cross-contamination. DNA was extracted from samples using the MagNA Pure LC 2.0 (Roche Diagnostics, Mannheim, Germany) for MiSeq and the MagNA Pure 24 (Roche) for ONT. For both ONT and MiSeq, PCR primers amplified an 820 base pair (bp) region of the HBV polymerase gene on the Roche LightCycler® 480. Primer sequence (5’ → 3’) forward (ILF) CGT GGT GGA CTT CTC TCA ATT TTC and reverse (ILR) AGA AAG GCC TTG TAA GTT GGC GA. For MiSeq, the second nested PCR reaction mixture consisted of the first round PCR product, and fusion indexed and barcoded primers. The nested PCR final amplicon was 364 bp. ONT AVR sequencing was performed using the GridION with R10.4.1 flow cells and the SQK-NBD114.24 library kit (Oxford Nanopore Technologies). Guppy (v. 6.4.6) was utilized for sequence base-calling. For MiSeq and ONT, mutational analysis and genotyping of FASTQ files was performed with DeepChek® HBV v.2.0 (ABL SA Group). Interpretation of clinically significant AVR mutations was based on the 2017 EASL Clinical Practice Guidelines (2). The 2nd World Health Organization (WHO) International Standard for Hepatitis B Virus DNA was used for evaluation of the limit of detection for the initial HBV PCR. Serial dilutions of the WHO Standard were tested with 5 replicates at: 1000, 500, 250, 125, 62.5, and 30 IU/mL. Probit regression was performed to determine limit of detection. A total of 49 patient samples and 15 external quality assurance (EQA) samples (QCMD, Glasgow, United Kingdom) were selected for AVR testing by NGS, to include a broad range of AVR mutations, at varying viral loads and percentage mutation frequency. There were 51 patient and 35 EQA samples available for genotype testing. The ATCC 45020D plasmid DNA was utilized as a control in all runs to determine error rates for base calling. To evaluate reproducibility/precision, a well characterized patient sample (genotype C; viral load of 9.9 × 108 IU/mL) was utilized, with mutations at approximately 70%, 51%, 26%, and 7%, at codons M204I, V173L, L180M, and T184S, respectively. The patient sample was diluted in Basematrix plasma (Seracare, Milford, MA) to 10,000, 4000, 1000, 250 and 125 IU/mL. These were tested in 3-5 total replicates over 2 separate NGS runs to determine the reproducibility and sensitivity of NGS at detecting different percent mutations at various viral loads. In order to assess cross-contamination risk, each run contained a negative control. ## 3. Results There was a high agreement rate between MiSeq and ONT with regards to total mutations detected and patient sample agreement, 68/70 (97%) and 47/49 (96%), respectively. All EASL mutations were included in the study, with the exceptions of M250V, N236T and T184G (Table 1) (2). View this table: [Table 1.](http://medrxiv.org/content/early/2024/10/02/2024.10.01.24314730/T1) Table 1. Hepatitis B Antiviral Resistance Mutations Identified by Next Generation Sequencing There were two patient samples with a mutation detected only by MiSeq (Table 2). The discrepant mutation results were found with low viral load samples (<2000 IU/mL). For one sample, the mutant subpopulation was detected at a low frequency (4.3%, M204I) only on MiSeq and would not have changed the genotypic resistance interpretation as a M204V mutation was also detected at a high frequency on both MiSeq and ONT platforms. For the second sample, the A181T mutation was found on MiSeq at a higher frequency (38%). QCMD samples with known HBV AVR mutations were tested by MiSeq and ONT. There was 80% (12/15) agreement for both platforms, with 3 samples having a L80I/V mutation identified only by the ONT platform. This is a known limitation of the MiSeq assay, as the sequenced amplicon is shorter and does not cover that particular mutation. View this table: [Table 2.](http://medrxiv.org/content/early/2024/10/02/2024.10.01.24314730/T2) Table 2. Discordant Result Analysis There were 35 EQA samples tested with the following genotypes: A, B, C, D, F. All were identified accurately by MiSeq and ONT. A total of 51 patient samples were tested for genotype by both MiSeq and ONT, and all results were concordant. These were composed of HBV genotypes: A (9), B (12), C (22), D (5), E (2), and A/G (1). ONT showed reproducible sequence variant detection for replicates with high viral loads and/or high-level sequence variants (Table 3). Sequence variant detection showed increased variability with low viral loads samples (<1000 IU/mL), particularly with low-level frequency subpopulations (<10%). The limit of detection for the initial PCR was 119 IU/mL. From the plasmid control, average base calling error rates at AVR codons were: V173 (0.74%), L180 (0.84%), M204 (0.36%) and M250 (0.42%). There was no evidence of sample cross- contamination risk, with each runs negative control testing negative. View this table: [Table 3.](http://medrxiv.org/content/early/2024/10/02/2024.10.01.24314730/T3) Table 3. Reproducibility in minority variant analysis with decreasing viral loads The ONT assay PCR and library preparation time was approximately 30% less than the MiSeq assay (5 versus 7 hours for one sample with controls); however, the time required for library preparation is also dependent on the number of samples included with each sequencing run. Sequencing time is markedly reduced with the ONT platform: approximately 15-30 minutes for ONT as compared to 23 hours for MiSeq assay. The time required for bioinformatic analysis for each sample was comparable between the assays (∼15 minutes). ## 4. Discussion This study describes the clinical validation of a NGS HBV AVR assay using the ONT platform. Results for ONT were highly concordant with MiSeq, with 100% (51/51) accuracy for genotype and 96% (47/49) accuracy for AVR testing on patient samples. Two samples had an additional AVR mutation detected by MiSeq which were not detected by ONT. These additional mutations were detected only on MiSeq from samples with low viral loads (<2000 IU/mL). The increased detection of low-level subpopulations by the MiSeq assay could be a result of the increased read depth per run with the MiSeq platform (∼12,500) as compared to ONT (∼8000). Additionally, the MiSeq assay used a more sensitive nested PCR for amplicon sequencing, which may explain improved detection of AVR mutations in low viral load samples. The nested PCR step was discontinued for the ONT assay, as it was more labour intensive, and had increased risks for sample cross-contamination. In addition, for samples with HBV viral load <500 IU/mL, our laboratory routinely reports a qualifying comment regarding decreased sensitivity for detecting mutant subpopulations of <10%, which is consistent with the reproducibility results from Table 3. ONT was able to improve PCR/library preparation and sequencing TAT, which has been previously demonstrated (13). Using the 45020D plasmid DNA control, which has a known sequence, ONT sequencing error rates were found to be less than 1% (0.36-0.84%). This is consistent with error rates reported for the R10.4.1 flow cells (9). Sequencing error rates are minimized by further quality filtering by Q score and forward/reverse read bias performed during bioinformatic analysis. This allows for confidence in reporting mutations found at a frequency of 3.5% or greater. Although resistance to antiviral agents is uncommon with entecavir and tenofovir, particularly in treatment-naïve patients, accessibility to HBV AVR is required for clinical follow up of patients with virological breakthrough to guide treatment (2–4). While *de novo* resistance to entecavir is rare, resistance can develop over time with chronic entecavir treatment: 0.9% at ≥5 years for treatment naïve and 20.1% for nucleos/tide analogue experienced (14). Limited data exist regarding tenofovir resistance, though suspected mutations have been reported which were not described in the EASL guidance (2, 5). In addition to tenofovir, novel antiviral agents aiming for functional cure are in the pipeline, and AVR sequencing will be needed to provide ongoing surveillance for any potential novel mutations associated with drug resistance (15, 16). There are practical barriers to the introduction of NGS in clinical laboratories, including cost, technical expertise and data analysis/interpretation (17, 18). As a clinical laboratory without a dedicated bioinformatician, we utilized a commercial platform to analyze FASTQ sequencing data and to remove the barrier relating to bioinformatics. Laboratories can also consider other commercial bioinformatic tools, which may increasingly become available as more laboratories launch NGS testing programs (19). With developing infrastructure and demand for clinical sequencing, laboratories may also consider the feasibility of in-house developed bioinformatic pipelines to enable customization of sequence interpretation, particularly for novel drugs/mutations (20). Limitations of our study include, despite attempts to include a broad array of AVR mutations and genotypes (patient and EQA samples), less frequently encountered AVR mutations and genotypes could not be included. In addition, the clinical relevance of identifying low-level subpopulations of AVR mutations through NGS is still unclear, and further research is required to determine relevant thresholds for identification of low-level subpopulations. In conclusion, a single real-time PCR followed by NGS with the ONT platform was found to be highly accurate for the detection of HBV AVR mutations and HBV genotype testing in a clinical laboratory. This was consistent across a wide variety of HBV genotypes and resistance mutations. Benefits of the ONT HBV assay included increased read length, decreased sample processing time, and improved TAT. ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## Conflict of Interest No relevant conflict of interests to declare. ## Funding This research did not receive any specific grant from funding agencies in the public, commercial or not-for-profit sectors. ## Authorship Conceptualization, MP, GR, NM, CFL; Data curation, MP, GR, TL, MY, WJ; Formal analysis, MP, GR; Sequencing and analysis, GR; Funding acquisition, N/A; Investigation, MP, GR, AS, MGR, NM, CFL; Methodology, GR, TL, MY; Project administration, WJ, MGR, NM, CFL; Supervision, NM, CFL; Visualization, MP; Writing – original draft, MP; Writing – review & editing, MP, GR, TL, MY, WJ, AS, MGR, NM, CFL. * Received October 1, 2024. * Revision received October 1, 2024. * Accepted October 2, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.World Health Organization. 2024. Guidelines for the prevention, diagnosis, care and treatment for people with chronic hepatitis B infection. 2. 2.European Association for the Study of the Liver. 2017. EASL 2017 Clinical Practice Guidelines on the management of hepatitis B virus infection. J Hepatol 67:370–398. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2017.03.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28427875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F10%2F02%2F2024.10.01.24314730.atom) 3. 3.Coffin CS, Fung SK, Alvarez F, Cooper CL, Doucette KE, Fournier C, Kelly E, Ko HH, Ma MM, Martin SR, Osiowy C, Ramji A, Tam E, Villeneuve JP. 2018. Management of Hepatitis B Virus Infection: 2018 Guidelines from the Canadian Association for the Study of Liver Disease and Association of Medical Microbiology and Infectious Disease Canada. Can liver J 1:156–217. 4. 4.Terrault NA, Lok ASF, McMahon BJ, Chang K-M, Hwang JP, Jonas MM, Brown RSJ, Bzowej NH, Wong JB. 2018. Update on Prevention, Diagnosis, and Treatment of Chronic Hepatitis B: AASLD 2018 Hepatitis B Guidance. Clin liver Dis 12:33–34. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F10%2F02%2F2024.10.01.24314730.atom) 5. 5.Mokaya J, McNaughton AL, Bester PA, Goedhals D, Barnes E, Marsden BD, Matthews PC. 2020. Hepatitis B virus resistance to tenofovir: fact or fiction? A systematic literature review and structural analysis of drug resistance mechanisms. Wellcome open Res. England. 6. 6.Solmone M, Vincenti D, Prosperi MCF, Bruselles A, Ippolito G, Capobianchi MR. 2009. Use of Massively Parallel Ultradeep Pyrosequencing To Characterize the Genetic Diversity of Hepatitis B Virus in Drug-Resistant and Drug-Naive Patients and To Detect Minor Variants in Reverse Transcriptase and Hepatitis B S Antigen. J Virol 83:1718– 1726. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjk6IjgzLzQvMTcxOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzEwLzAyLzIwMjQuMTAuMDEuMjQzMTQ3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 7. 7.Qutub MO, Germer JJ, Rebers SPH, Mandrekar JN, Beld MGHM, Yao JDC. 2006. Simplified PCR protocols for INNO-LiPA HBV Genotyping and INNO-LiPA HBV PreCore assays. J Clin Virol Off Publ Pan Am Soc Clin Virol 37:218–221. 8. 8.Lowe CL, Merrick L, Harrigan PR, Mazzulli T, Sherlock CH, Ritchie G. 2016. Next-generation sequencing for hepatitis B genotype and resistance testing in a clinical microbiology laboratory. J Clin Microbiol 54:127–133. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjg6IjU0LzEvMTI3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMTAvMDIvMjAyNC4xMC4wMS4yNDMxNDczMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. 9.Kim BY, Gellert HR, Church SH, Suvorov A, Anderson SS, Barmina O, Beskid SG, Comeault AA, Crown KN, Diamond SE, Dorus S, Fujichika T, Hemker JA, Hrcek J, Kankare M, Katoh T, Magnacca KN, Martin RA, Matsunaga T, Medeiros MJ, Miller DE, Pitnick S, Simoni S, Steenwinkel TE, Schiffer M, Syed ZA, Takahashi A, Wei KH-C, Yokoyama T, Eisen MB, Kopp A, Matute D, Obbard DJ, O’Grady PM, Price DK, Toda MJ, Werner T, Petrov DA. 2023. Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life. bioRxiv Prepr Serv Biol. United States. 10. 10.FDA US. 2016. Infectious disease next generation sequencing based diagnostic devices: microbial identification and detection of antimicrobial resistance and virulence markers: draft guidance for industry and Food and Drug Administration staff. U.S. Food and Drug Administration, Silver Spring, MD. 11. 11.Burd EM. 2010. Validation of Laboratory-Developed Molecular Assays for Infectious Diseases. Clin Microbiol Rev 23:550–576. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjg6IjIzLzMvNTUwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMTAvMDIvMjAyNC4xMC4wMS4yNDMxNDczMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. 12.Gargis AS, Kalman L, Lubin M. 2016. Assuring the Quality of Next-Generation Sequencing in Clinical Microbiology and Public Health Laboratories. J Clin Microbiol 54:2857–2865. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjEwOiI1NC8xMi8yODU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMTAvMDIvMjAyNC4xMC4wMS4yNDMxNDczMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. 13.Colson P, Boschi C, Bengone-Abogourin JG, Brechard L, Motte A, Allemand I. 2021. Concurrent Nanopore Next-Generation Sequencing of Hepatitis B and Delta Virus Genomes Directly From Patient Plasma. Ann Lab Med. Korea (South). 14. 14.Lumley SF, Delphin M, Mokaya JF, Tan CCS, Martyn E, Anderson M, Li KC, Waddilove E, Sukali G, Downs LO, Said K, Okanda D, Campbell C, Harriss E, Shimakawa Y, Matthews PC. 2024. A systematic review and meta-analysis of the risk of hepatitis B virus (HBV) resistance in people treated with entecavir or tenofovir. J Clin Virol Off Publ Pan Am Soc Clin Virol 174:105711. 15. 15.Wong GLH, Gane E, Lok ASF. 2022. How to achieve functional cure of HBV: Stopping NUCs, adding interferon or new drug development? J Hepatol 76:1249–1262. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2021.11.024&link_type=DOI) 16. 16.Yuen M-F, Lim S-G, Plesniak R, Tsuji K, Janssen HLA, Pojoga C, Gadano A, Popescu CP, Stepanova T, Asselah T, Diaconescu G, Yim HJ, Heo J, Janczewska E, Wong A, Idriz N, Imamura M, Rizzardini G, Takaguchi K, Andreone P, Arbune M, Hou J, Park SJ, Vata A, Cremer J, Elston R, Lukić T, Quinn G, Maynard L, Kendrick S, Plein H, Campbell F, Paff M, Theodore D. 2022. Efficacy and Safety of Bepirovirsen in Chronic Hepatitis B Infection. N Engl J Med 387:1957–1968. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2210027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36346079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F10%2F02%2F2024.10.01.24314730.atom) 17. 17.Boers SA, Jansen R, Hays JP. 2019. Understanding and overcoming the pitfalls and biases of next-generation sequencing (NGS) methods for use in the routine clinical microbiological diagnostic laboratory. Eur J Clin Microbiol Infect Dis Off Publ Eur Soc Clin Microbiol 38:1059–1070. 18. 18.Greninger AL. 2018. The challenge of diagnostic metagenomics. Expert Rev Mol Diagn 18:605–615. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/14737159.2018.1487292&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F10%2F02%2F2024.10.01.24314730.atom) 19. 19.Ringlander J, Andersson ME, Prakash K, Larsson SB, Lindh M. 2022. Deep sequencing of hepatitis B virus using Ion Torrent fusion primer method. J Virol Methods 299:114315. 20. 20.Wang H. 2024. Practical updates in clinical antiviral resistance testing. J Clin Microbiol 62:e0072823.