{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"institution":[{"name":"bioRxiv"}],"indexed":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T19:34:18Z","timestamp":1768505658738,"version":"3.49.0"},"posted":{"date-parts":[[2023,1,12]]},"group-title":"Genomics","reference-count":87,"publisher":"openRxiv","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"accepted":{"date-parts":[[2023,1,12]]},"abstract":"<jats:title>ABSTRACT<\/jats:title>\n                <jats:p>Hepatitis B virus (HBV) remains a global public health concern, with over 250 million individuals living with chronic HBV infection (CHB) and no curative therapy currently available. Viral diversity is associated with CHB pathogenesis and immunological control of infection. Improved methods to characterize the viral genome at both the population and intra-host level could aid drug development efforts. Conventionally, HBV sequencing data are aligned to a linear reference genome and only sequences capable of aligning to the reference are captured for analysis. Reference selection has additional consequences, including sample-specific \u2018consensus\u2019 sequence construction. It remains unclear how to select a reference from available sequences and whether a single reference is sufficient for genetic analyses. Using simulated short-read sequencing data generated from full-length publicly available HBV genome sequences and HBV sequencing data from a longitudinally sampled individual with CHB, we investigate alternative graph-based alignment approaches. We demonstrate that using a phylogenetically representative \u2018genome graph\u2019 for alignment, rather than linear reference sequences, avoids issues of reference ambiguity, improves alignment, and facilitates the construction of sample-specific consensus sequences genetically similar to an individual\u2019s infection. Graph-based methods can therefore improve efforts to characterize the genetics of viral pathogens, including HBV, and may have broad implications in host pathogen research.<\/jats:p>","DOI":"10.1101\/2023.01.11.523611","type":"posted-content","created":{"date-parts":[[2023,1,12]],"date-time":"2023-01-12T22:55:18Z","timestamp":1673564118000},"source":"Crossref","is-referenced-by-count":0,"title":["A hepatitis B virus (HBV) sequence variation graph improves sequence alignment and sample-specific consensus sequence construction for genetic analysis of HBV"],"prefix":"10.64898","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-5629-2012","authenticated-orcid":false,"given":"Dylan","family":"Duchen","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Steven","family":"Clipman","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Candelaria","family":"Vergara","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Chloe L.","family":"Thio","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"David L.","family":"Thomas","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Priya","family":"Duggal","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Genevieve L.","family":"Wojcik","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"54368","reference":[{"key":"2023011403100559000_2023.01.11.523611v1.1","unstructured":"WHO. Preventing Perinatal Hepatitis B Virus Transmission \u25a1: A Guide for Introducing and Strengthening Hepatitis B Birth Dose Vaccination. (2015)."},{"key":"2023011403100559000_2023.01.11.523611v1.2","doi-asserted-by":"crossref","first-page":"883","DOI":"10.1016\/S2468-1253(19)30190-6","article-title":"Targets and future direct-acting antiviral approaches to achieve hepatitis B virus cure","volume":"4","year":"2019","journal-title":"Lancet Gastroenterol. Hepatol"},{"key":"2023011403100559000_2023.01.11.523611v1.3","doi-asserted-by":"publisher","DOI":"10.3389\/fgene.2019.00735"},{"key":"2023011403100559000_2023.01.11.523611v1.4","doi-asserted-by":"publisher","DOI":"10.1016\/S0140-6736(14)60220-8"},{"key":"2023011403100559000_2023.01.11.523611v1.5","doi-asserted-by":"crossref","first-page":"10529","DOI":"10.1038\/s41598-019-46609-7","article-title":"Large-scale viral genome analysis identifies novel clinical associations between hepatitis B virus and chronically infected patients","volume":"9","year":"2019","journal-title":"Sci. Rep"},{"key":"2023011403100559000_2023.01.11.523611v1.6","doi-asserted-by":"publisher","DOI":"10.1016\/0016-5085(90)90632-B"},{"key":"2023011403100559000_2023.01.11.523611v1.7","doi-asserted-by":"crossref","unstructured":"G\u00fcnther, S. et al. Type, prevalence, and significance of core promoter\/enhancer II mutations in hepatitis B viruses from immunosuppressed patients with severe liver disease. J. Virol. (1996).","DOI":"10.1128\/jvi.70.12.8318-8331.1996"},{"key":"2023011403100559000_2023.01.11.523611v1.8","doi-asserted-by":"publisher","DOI":"10.1099\/0022-1317-79-2-375"},{"key":"2023011403100559000_2023.01.11.523611v1.9","doi-asserted-by":"publisher","DOI":"10.1086\/592571"},{"key":"2023011403100559000_2023.01.11.523611v1.10","doi-asserted-by":"publisher","DOI":"10.3748\/wjg.15.5761"},{"key":"2023011403100559000_2023.01.11.523611v1.11","first-page":"384","article-title":"Insights From Deep Sequencing of the HBV Genome\u2014 Unique, Tiny, and Misunderstood","volume":"156","year":"2018","journal-title":"Gastroenterology"},{"key":"2023011403100559000_2023.01.11.523611v1.12","doi-asserted-by":"publisher","DOI":"10.1002\/hep.21188"},{"key":"2023011403100559000_2023.01.11.523611v1.13","doi-asserted-by":"publisher","DOI":"10.1002\/jmv.20985"},{"key":"2023011403100559000_2023.01.11.523611v1.14","doi-asserted-by":"crossref","first-page":"1307","DOI":"10.1002\/jmv.23965","article-title":"Molecular epidemiological study of hepatitis B virus genotypes in Southwest, China","volume":"86","year":"2014","journal-title":"J. Med. Virol"},{"key":"2023011403100559000_2023.01.11.523611v1.15","doi-asserted-by":"crossref","first-page":"26","DOI":"10.1016\/j.meegid.2017.10.023","article-title":"Novel HBV recombinants between genotypes B and C in 3\u2032-terminal reverse transcriptase (RT) sequences are associated with enhanced viral DNA load, higher RT point mutation rates and place of birth among Chinese patients","volume":"57","year":"2018","journal-title":"Infect. Genet. Evol"},{"key":"2023011403100559000_2023.01.11.523611v1.16","doi-asserted-by":"publisher","DOI":"10.1128\/JVI.02556-07"},{"key":"2023011403100559000_2023.01.11.523611v1.17","doi-asserted-by":"publisher","DOI":"10.1128\/JVI.00462-09"},{"key":"2023011403100559000_2023.01.11.523611v1.18","doi-asserted-by":"publisher","DOI":"10.1016\/j.ijid.2010.03.020"},{"key":"2023011403100559000_2023.01.11.523611v1.19","doi-asserted-by":"publisher","DOI":"10.1016\/j.meegid.2013.03.021"},{"key":"2023011403100559000_2023.01.11.523611v1.20","doi-asserted-by":"publisher","DOI":"10.1128\/MMBR.05023-11"},{"key":"2023011403100559000_2023.01.11.523611v1.21","doi-asserted-by":"publisher","DOI":"10.1016\/j.coviro.2017.03.013"},{"key":"2023011403100559000_2023.01.11.523611v1.22","doi-asserted-by":"crossref","first-page":"661","DOI":"10.1038\/s41598-017-00694-8","article-title":"Evolution of full-length genomes of HBV quasispecies in sera of patients with a coexistence of HBsAg and anti-HBs antibodies","volume":"7","year":"2017","journal-title":"Sci. Rep"},{"key":"2023011403100559000_2023.01.11.523611v1.23","doi-asserted-by":"publisher","DOI":"10.1128\/JCM.00068-15"},{"key":"2023011403100559000_2023.01.11.523611v1.24","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pgen.1008271"},{"key":"2023011403100559000_2023.01.11.523611v1.25","doi-asserted-by":"publisher","DOI":"10.1128\/JVI.01123-14"},{"key":"2023011403100559000_2023.01.11.523611v1.26","doi-asserted-by":"crossref","first-page":"1064","DOI":"10.1007\/s00535-017-1334-1","article-title":"Deep sequencing analysis of quasispecies in the HBV pre-S region and its association with hepatocellular carcinoma","volume":"52","year":"2017","journal-title":"J. Gastroenterol"},{"key":"2023011403100559000_2023.01.11.523611v1.27","doi-asserted-by":"publisher","DOI":"10.1136\/gutjnl-2012-302408"},{"key":"2023011403100559000_2023.01.11.523611v1.28","doi-asserted-by":"crossref","first-page":"201","DOI":"10.1016\/j.cmi.2015.10.007","article-title":"Increased intrahepatic quasispecies heterogeneity correlates with off-treatment sustained response to nucleos(t)ide analogues in e antigen-positive chronic hepatitis B patients","volume":"22","year":"2016","journal-title":"Clin. Microbiol. Infect"},{"key":"2023011403100559000_2023.01.11.523611v1.29","doi-asserted-by":"publisher","DOI":"10.1136\/gut.2010.226225"},{"key":"2023011403100559000_2023.01.11.523611v1.30","doi-asserted-by":"publisher","DOI":"10.1046\/j.1365-2893.1997.00066.x"},{"key":"2023011403100559000_2023.01.11.523611v1.31","doi-asserted-by":"publisher","DOI":"10.1038\/nbt.4227"},{"key":"2023011403100559000_2023.01.11.523611v1.32","doi-asserted-by":"crossref","first-page":"147","DOI":"10.1007\/s12072-015-9645-x","article-title":"Aligning to the sample-specific reference sequence to optimize the accuracy of next-generation sequencing analysis for hepatitis B virus","volume":"10","year":"2016","journal-title":"Hepatol. Int"},{"key":"2023011403100559000_2023.01.11.523611v1.33","doi-asserted-by":"crossref","first-page":"e1008678","DOI":"10.1371\/journal.pcbi.1008678","article-title":"One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads","volume":"17","year":"2021","journal-title":"PLOS Comput. Biol"},{"key":"2023011403100559000_2023.01.11.523611v1.34","doi-asserted-by":"crossref","first-page":"267","DOI":"10.1186\/s13059-021-02473-1","article-title":"Pandora: nucleotide-resolution bacterial pan-genomics with reference graphs","volume":"22","year":"2021","journal-title":"Genome Biol"},{"key":"2023011403100559000_2023.01.11.523611v1.35","doi-asserted-by":"publisher","DOI":"10.1101\/2022.03.10.483737"},{"key":"2023011403100559000_2023.01.11.523611v1.36","doi-asserted-by":"crossref","first-page":"271","DOI":"10.1099\/jgv.0.001387","article-title":"Analysis of genomic-length HBV sequences to determine genotype and subgenotype reference sequences","volume":"101","year":"2020","journal-title":"J. Gen. Virol"},{"key":"2023011403100559000_2023.01.11.523611v1.37","doi-asserted-by":"crossref","unstructured":"Eizenga, J. M. et al. Pangenome Graphs. Annu. Rev. Genomics Hum. Genet. 21, annurev-genom-120219-080406 (2020).","DOI":"10.1146\/annurev-genom-120219-080406"},{"key":"2023011403100559000_2023.01.11.523611v1.38","doi-asserted-by":"crossref","unstructured":"Eizenga, J. M. et al. Succinct dynamic variation graphs. 1\u20136 (2020).","DOI":"10.1101\/2020.04.23.056317"},{"key":"2023011403100559000_2023.01.11.523611v1.39","doi-asserted-by":"publisher","DOI":"10.1038\/ng.3257"},{"key":"2023011403100559000_2023.01.11.523611v1.40","doi-asserted-by":"crossref","first-page":"5077","DOI":"10.1038\/s41598-022-09035-w","article-title":"The ViReflow pipeline enables user friendly large scale viral consensus genome reconstruction","volume":"12","year":"2022","journal-title":"Sci. Rep"},{"key":"2023011403100559000_2023.01.11.523611v1.41","doi-asserted-by":"publisher","DOI":"10.1038\/nprot.2017.066"},{"key":"2023011403100559000_2023.01.11.523611v1.42","doi-asserted-by":"crossref","first-page":"601","DOI":"10.1186\/s12879-017-2697-x","article-title":"Molecular characterization of hepatitis B virus in Vietnam","volume":"17","year":"2017","journal-title":"BMC Infect. Dis"},{"key":"2023011403100559000_2023.01.11.523611v1.43","doi-asserted-by":"crossref","first-page":"7081","DOI":"10.1038\/s41598-019-43524-9","article-title":"Illumina and Nanopore methods for whole genome sequencing of hepatitis B virus (HBV)","volume":"9","year":"2019","journal-title":"Sci. Rep"},{"key":"2023011403100559000_2023.01.11.523611v1.44","doi-asserted-by":"publisher","DOI":"10.1126\/sciimmunol.aau6905"},{"key":"2023011403100559000_2023.01.11.523611v1.45","doi-asserted-by":"publisher","DOI":"10.14806\/EJ.17.1.200"},{"key":"2023011403100559000_2023.01.11.523611v1.46","unstructured":"Andrews, S. FastQC. Babraham Bioinforma. (2010)."},{"key":"2023011403100559000_2023.01.11.523611v1.47","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/bty630"},{"key":"2023011403100559000_2023.01.11.523611v1.48","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btaa777"},{"key":"2023011403100559000_2023.01.11.523611v1.49","doi-asserted-by":"publisher","DOI":"10.1101\/2022.02.14.4804"},{"key":"2023011403100559000_2023.01.11.523611v1.50","doi-asserted-by":"crossref","unstructured":"Sir\u00e9n, J. et al. Pangenomics enables genotyping of known structural variants in 5202 diverse genomes. Science (80-.). 374, (2021).","DOI":"10.1126\/science.abg8871"},{"key":"2023011403100559000_2023.01.11.523611v1.51","first-page":"400","article-title":"Haplotype-aware graph indexes","volume":"36","year":"2018","journal-title":"Bioinformatics"},{"key":"2023011403100559000_2023.01.11.523611v1.52","doi-asserted-by":"crossref","first-page":"3941","DOI":"10.21105\/joss.03941","article-title":"Rasusa: Randomly subsample sequencing reads to a specified coverage","volume":"7","year":"2022","journal-title":"J. Open Source Softw"},{"key":"2023011403100559000_2023.01.11.523611v1.53","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btp352"},{"key":"2023011403100559000_2023.01.11.523611v1.54","doi-asserted-by":"crossref","first-page":"495","DOI":"10.3390\/genes9100495","article-title":"The Global Hepatitis B Virus Genotype Distribution Approximated from Available Genotyping Data","volume":"9","year":"2018","journal-title":"Genes (Basel)"},{"key":"2023011403100559000_2023.01.11.523611v1.55","doi-asserted-by":"publisher","DOI":"10.1128\/JVI.76.12.5985-5992.2002"},{"key":"2023011403100559000_2023.01.11.523611v1.56","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btp696"},{"key":"2023011403100559000_2023.01.11.523611v1.57","unstructured":"Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv (2013)."},{"key":"2023011403100559000_2023.01.11.523611v1.58","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btv566"},{"key":"2023011403100559000_2023.01.11.523611v1.59","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btv383"},{"key":"2023011403100559000_2023.01.11.523611v1.60","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-020-02135-8"},{"key":"2023011403100559000_2023.01.11.523611v1.61","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-018-1618-7"},{"key":"2023011403100559000_2023.01.11.523611v1.62","doi-asserted-by":"publisher","DOI":"10.1093\/gigascience\/giab016"},{"key":"2023011403100559000_2023.01.11.523611v1.63","unstructured":"Bushnell, B. BBMap: A Fast, Accurate, Splice-Aware Aligner. (2014)."},{"key":"2023011403100559000_2023.01.11.523611v1.64","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btt593"},{"key":"2023011403100559000_2023.01.11.523611v1.65","doi-asserted-by":"publisher","DOI":"10.1093\/nar\/gks918"},{"key":"2023011403100559000_2023.01.11.523611v1.66","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-016-0997-x"},{"key":"2023011403100559000_2023.01.11.523611v1.67","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pone.0191962"},{"key":"2023011403100559000_2023.01.11.523611v1.68","doi-asserted-by":"publisher","DOI":"10.1186\/s40168-017-0283-5"},{"key":"2023011403100559000_2023.01.11.523611v1.69","doi-asserted-by":"publisher","DOI":"10.1101\/gr.215038.116"},{"key":"2023011403100559000_2023.01.11.523611v1.70","doi-asserted-by":"publisher","DOI":"10.1101\/645721"},{"key":"2023011403100559000_2023.01.11.523611v1.71","doi-asserted-by":"crossref","first-page":"217","DOI":"10.3390\/v14020217","article-title":"ViralFlow: A Versatile Automated Workflow for SARS-CoV-2 Genome Assembly, Lineage Assignment, Mutations and Intrahost Variant Detection","volume":"14","year":"2022","journal-title":"Viruses"},{"key":"2023011403100559000_2023.01.11.523611v1.72","doi-asserted-by":"crossref","first-page":"1321","DOI":"10.1038\/s41467-022-28852-1","article-title":"Rescuing low frequency variants within intra-host viral populations directly from Oxford Nanopore sequencing data","volume":"13","year":"2022","journal-title":"Nat. Commun"},{"key":"2023011403100559000_2023.01.11.523611v1.73","doi-asserted-by":"publisher","DOI":"10.1186\/s12864-017-3899-8"},{"key":"2023011403100559000_2023.01.11.523611v1.74","doi-asserted-by":"publisher","DOI":"10.1128\/JVI.00243-16"},{"key":"2023011403100559000_2023.01.11.523611v1.75","doi-asserted-by":"crossref","first-page":"146387","DOI":"10.1016\/j.gene.2022.146387","article-title":"Role of genomics in combating COVID-19 pandemic","volume":"823","year":"2022","journal-title":"Gene"},{"key":"2023011403100559000_2023.01.11.523611v1.76","first-page":"2","article-title":"GISAID: Global initiative on sharing all influenza data \u2013 from vision to reality","volume":"22","year":"2017","journal-title":"Eurosurveillance"},{"key":"2023011403100559000_2023.01.11.523611v1.77"},{"key":"2023011403100559000_2023.01.11.523611v1.78","doi-asserted-by":"publisher","DOI":"10.1056\/NEJMsr1813907"},{"key":"2023011403100559000_2023.01.11.523611v1.79","doi-asserted-by":"publisher","DOI":"10.1101\/2022.03.10.483737"},{"key":"2023011403100559000_2023.01.11.523611v1.80","doi-asserted-by":"publisher","DOI":"10.1101\/2021.02.02.429378"},{"key":"2023011403100559000_2023.01.11.523611v1.81","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-018-1595-x"},{"key":"2023011403100559000_2023.01.11.523611v1.82","doi-asserted-by":"crossref","unstructured":"Monsu, M. & Comin, M. Fast alignment of reads to a variation graph with application to SNP detection. J. Integr. Bioinform. 18, (2021).","DOI":"10.1515\/jib-2021-0032"},{"key":"2023011403100559000_2023.01.11.523611v1.83","doi-asserted-by":"crossref","unstructured":"Baaijens, J. A. et al. Computational graph pangenomics: a tutorial on data structures and their applications. Nat. Comput. 6, (2022).","DOI":"10.1007\/s11047-022-09882-6"},{"key":"2023011403100559000_2023.01.11.523611v1.84","doi-asserted-by":"publisher","DOI":"10.1101\/2022.04.27.489683"},{"key":"2023011403100559000_2023.01.11.523611v1.85","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-020-02229-3"},{"key":"2023011403100559000_2023.01.11.523611v1.86","first-page":"68","article-title":"GEDEVO: An evolutionary graph edit distance algorithm for biological network alignment","volume":"34","year":"2013","journal-title":"OpenAccess Ser. Informatics"},{"key":"2023011403100559000_2023.01.11.523611v1.87","doi-asserted-by":"publisher","DOI":"10.1101\/2022.02.18.481102"}],"container-title":[],"original-title":[],"link":[{"URL":"https:\/\/syndication.highwire.org\/content\/doi\/10.1101\/2023.01.11.523611","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T04:07:25Z","timestamp":1768450045000},"score":1,"resource":{"primary":{"URL":"http:\/\/biorxiv.org\/lookup\/doi\/10.1101\/2023.01.11.523611"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,1,12]]},"references-count":87,"URL":"https:\/\/doi.org\/10.1101\/2023.01.11.523611","relation":{"is-preprint-of":[{"id-type":"doi","id":"10.1371\/journal.pone.0301069","asserted-by":"subject"}]},"subject":[],"published":{"date-parts":[[2023,1,12]]},"subtype":"preprint"}}