key: cord-0833824-i3vp1ved
authors: Lee, Geum-Young; Kim, Won-Keun; Cho, Seungchan; Park, Kyungmin; Kim, Jongwoo; Lee, Seung-Ho; Lee, Jingyeong; Lee, Young-Sun; Kim, Ji Hoon; Byun, Kwan Soo; Song, Jin-Won
title: Genotyping and Molecular Diagnosis of Hepatitis A Virus in Human Clinical Samples Using Multiplex PCR-Based Next-Generation Sequencing
date: 2022-01-04
journal: Microorganisms
DOI: 10.3390/microorganisms10010100
sha: fd528ea0ed0c8c9eb4ec8f7af2d2c5581a11e049
doc_id: 833824
cord_uid: i3vp1ved

Hepatitis A virus (HAV) is a serious threat to public health worldwide. We used multiplex polymerase chain reaction (PCR)-based next-generation sequencing (NGS) to derive information on viral genetic diversity and conduct precise phylogenetic analysis. Four HAV genome sequences were obtained using multiplex PCR-based NGS. HAV whole-genome sequence of one sample was obtained by conventional Sanger sequencing. The HAV strains demonstrated a geographic cluster with sub-genotype IA strains in the Republic of Korea. The phylogenetic pattern of HAV viral protein (VP) 3 region showed no phylogenetic conflict between the whole-genome and partial-genome sequences. The VP3 region in serum and stool samples showed sensitive detection of HAV with differences of quantification that did not exceed <10 copies/μL than the consensus VP4 region using quantitative PCR (qPCR). In conclusion, multiplex PCR-based NGS was implemented to define HAV genotypes using nearly whole-genome sequences obtained directly from hepatitis A patients. The VP3 region might be a potential candidate for tracking the genotypic origin of emerging HAV outbreaks. VP3-specific qPCR was developed for the molecular diagnosis of HAV infection. This study may be useful to predict for the disease management and subsequent development of hepatitis A infection at high risk of severe illness.

Hepatitis A virus (HAV) infections are the critical etiology of viral hepatitis and impose a massive socioeconomic burden worldwide [1] . HAV is mainly transmitted via the fecaloral route through contaminated food and water. The prevalence rate of HAV in different countries varies with hygiene levels, and approximately, 1.5 million people are annually infected with HAV worldwide. In the late 1970s and early 1990s, 85-95% of the population were serologically positive for anti-HAV immunoglobulin in Republic of Korea (ROK), China, Japan, Philippines, and Thailand [2] . The seroprevalence of anti-HAV antibodies in the Korean population rapidly declined from more than 80% in the 1970s to less than 20% HAV (genus, Hepatovirus; family, Picornaviridae) is a non-enveloped virus with a singlestranded 7.5 kb genome that contains 5' and 3' noncoding regions, one open reading frame (ORF), and a poly(A) tract [5] . The 5' untranslated region (UTR) of the HAV RNA harbors an internal ribosome entry site, which directs the ribosomes to the initiation codon of the polyprotein. The single ORF is divided into three functional regions: the P1 region that encodes the capsid polypeptides viral protein 1 (VP1) to VP4 and the P2 (2A, 2B, and 2C) and P3 (3A, 3B, 3C, and 3D) regions that encode nonstructural polypeptides associated with viral replication. Globally spread HAV strains display significant variability in RNA genomic sequences, despite the limited heterogeneity of amino acids [6] . In 1987, a genetic analysis of HAV strains was conducted for seven genotypes including human (I-III and VII) and simian (IV-VI) groups. The HAV strains have six sub-genotypes (IA, IB, IIA, IIB, IIIA, and IIIB) in humans. Phylogenetic and genotypic analyses of HAV were performed using partial-genome sequences, including the VP1/P2A junction [5, 7] , VP1/P2B junction [8, 9] , C termini of VP3 [7] , N terminus of VP1 [8, 10] , entire VP1 [11, 12] , and 5' UTR [13] regions. In clinical microbiology, the partial-genome sequences of HAV exhibited inaccurate assessment of phylogeographic relationships or potential transmission routes in an outbreak [14] . Therefore, the precise genotypic relationship of global HAV strains with partial-genome sequences remains to be investigated.

Genomic epidemiology of HAV plays a critical role in identifying the epidemiologic surveillance, microbial source tracking, and pathogen identification of food-borne outbreak [4] . The genomic epidemiology of emerging HAV in eastern Spain demonstrated a food-borne origin of the virus that caused the outbreak and patients with hepatitis A [15, 16] . Epidemiological analysis of tourists returning from Egypt identified that orange juice was the source of infection in the hepatitis A outbreak [17] . Genomic epidemiology of HAV IB genotype reported the causes of an outbreak associated with homelessness and drug abuse in USA [14] . Epidemiologic investigations of hepatitis A outbreaks revealed that hepatitis A cases had consumed semi-dried tomatoes in Australia, Netherlands, and England [18] [19] [20] . However, identification and sequencing of HAV remains challenging due to the low titer of pathogens in human or food samples.

Next-generation sequencing (NGS) has a broad range of applications, including simultaneous multiple pathogens detection, phylogenetic diversity estimation, and wholegenome sequencing (WGS) [21] [22] [23] [24] [25] [26] [27] [28] . Metagenomic NGS combined with spiked primer enrichment and/or tiling multiplex PCR recovered viral genomes from coronavirus disease 19 patients, revealing the genetic diversity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the transmission of virus lineages [29] . Multiplex polymerase chain reaction (PCR)-based NGS acquired the nearly WGS of Hantaan virus (HTNV) directly from hemorrhagic fever with renal syndrome (HFRS) patient samples and rodent tissues [27] . These studies demonstrated the putative infection site of HTNV by phylogeographic analysis in HFRS-endemic areas. NGS allowed nearly WGS of HAV in frozen berries linked to HAV outbreaks in Italy [30] .

In this study, multiplex PCR-based NGS facilitated nearly WGS of HAV directly from samples collected from patients with hepatitis A. The results demonstrated that HAV VP3 gene might be a potential candidate for tracing the genotypic origin of an outbreak, showing a well-established phylogenic pattern compared with the whole-genome sequences of HAV. Overall, our findings provide significant insights into the genomic epidemiology, genetic identification, and molecular diagnosis of hepatitis A infection. 

Written informed consent was obtained from all participants, and the Institutional Review Board of the Korea University Guro Hospital (2008GR0015) approved this study.

Clinical samples were obtained from the Korea University Guro Hospital, Seoul, Korea. Hepatitis A infection was diagnosed based on laboratory and clinical tests in five patients. Serum and stool samples from patients with anti-HAV immunoglobulin M (IgM) antibodies were used in this study. All clinical samples were collected from patients of 22-38 years in age.

Total RNA was extracted using TRI Reagent LS (Ambion, Austin, TX, USA). Reverse transcription was performed using 1 µg of RNA with a High-Capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA, USA) and random hexamers. Primer sequences used for RT-PCR included HAV VP3-F1 (forward), 5 -GCT TTG GAT CAG GAA GAT TGG A-3 ; HAV VP3-F2 (forward), 5 -ACT CAT TTT ACT TUT TTG GAC ATC CA-3 ; and HAV VP3-R1 (reverse), 5 -CAT GAT AAA GAG GAG CAA AAC ATT CC-3 . The first and second RT-PCRs were performed at 94 • C for 5 min, followed by six cycles at 94 • C for 40 s, 40-45 • C for 1 min, and 72 • C for 1-2 min; 32 cycles at 94 • C for 40 s, 42 • C for 40 s, and 72 • C for 1-2 min; and a final cycle at 72 • C for 5 min (ProFlex PCR System, Life Technology, CA, USA).

The recombinant plasmid DNA of HAV VP3 gene was isolated using the pTOP Blunt V2 vector (Enzynomics Co., Ltd., Seoul, Korea). The concentration of recombinant plasmid DNA was measured by UV absorbance at 260 nm and 280 nm using Nano drop (Thermo Fisher Scientific, CA, USA). Standard curve was obtained from serial 10-fold dilutions of recombinant plasmid DNA ranging from 1 × 10 10 to 1 × 10 copies/µL. The copy number of plasmids per microgram of DNA was calculated using the total number of nucleotides in the plasmid using a previously described formula [31] . The HAV-specific primer set was designed based on the VP3 region: HAV VP3-qF (forward), 5 -ATG AAG ATG CAA GGG CAA A-3 and HAV VP3-qR (reverse), 5 -GGA ATG GAT GTC CAA GTA GTA AA-3 . Viral RNA quantification was compared with the qPCR assay of the VP4 genomic region targeting the HAV genotype IA [32] . qPCR was performed at 95 • C for 10 min, followed by 40 cycles at 95 • C for 15 s and 60 • C for 1 min using a QuantStudio 5 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) with a Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA).

For the whole-genome strategy, we used several sets of primers to amplify overlapping fragments of 200-500 bp in length spanning the entire viral genome (Table S1 ). The cDNA library was enriched using Solg 2× Uh-Taq PCR Smart mix (SolGent, Daejeon, Korea), according to the manufacturer's instructions. The reaction mixture of 25 µL contained 12.5 µL of 2× Uh pre-mix, 2 µL of each primer, 9.5 µL of distilled water, and 1 µL of DNA template. Multiplex PCR was performed at 95 • C for 15 min, followed by 40 and/or 25 cycles at 95 • C for 20 s, 50 • C for 40 s, 72 • C for 1 min, and a final cycle at 72 • C for 3 min.

DNA libraries were prepared using a TruSeq Nano DNA LT Sample Preparation Kit (Illumina, San Diego, CA, USA), according to the manufacturer's instructions. To obtain size-selected amplicons, cDNA templates were mechanically sheared using an M220 focused ultrasonicator (Covaris, Woburn, MA, USA). The cDNA amplicons were prepared by size selection, A-tailing, and ligation with indexes and adaptors. The enrichment reaction contained 5 µL of PCR primers and 20 µL of the enhanced PCR mixture (Illumina, San Diego, CA, USA). Library quality was evaluated using a bioanalyzer with an Agilent DNA 1000 Kit (Agilent Technologies, Santa Clara, CA, USA). NGS was performed using a 2 × 150 bp MiSeq benchtop sequencer (Illumina, San Diego, CA, USA) with an MiSeq Reagent Kit V2 (Illumina, San Diego, CA, USA).

Whole-genome sequences of HAV were aligned using ClustalW in Lasergene 5 (DNAS-TAR, Madison, WI, USA). The phylogenetic tree was generated using the best-fit general time reversible+gamma+invariable model of evolution. Support for topologies was assessed using bootstrapping for 1000 iterations. After model optimization for each data set, the pairwise genetic distance between HAV strains was calculated using MEGA 7.0 [33] .

To compare the phylogenetic patterns in the whole-genome and partial-genome sequences of HAV, two phylogenetic trees were generated using the tanglegram algorithm for matching HAV strains. The auxiliary lines in the center connected the phylogenetic trees. Tanglegrams were generated for all phylogenetic links between the whole-genome and partial-genome sequences using the 'dendextend' package of R [34] .

Five patients exhibited symptoms of acute hepatitis A ( Table 1 ). All patients who tested positive for anti-HAV IgM showed significant elevation in aspartate aminotransferase, alanine aminotransferase, and alkaline phosphatase. All patients manifested typical symptoms compatible with acute viral hepatitis A (e.g., icterus, dark urine, anorexia, and malaise). All five patients manifested hyperbilirubinemia (serum bilirubin levels > 1.2 mg/dl). None of the participants had other chronic liver diseases such as chronic hepatitis B, chronic hepatitis C, autoimmune hepatitis, or primary biliary cholangitis. 

The viral load of the four serum and stool samples was quantified by HAV VP3-specific qPCR ( Figure 1 ). The copy number of HAV RNAs was determined using a linear regression curve with a coefficient of correlation (r) value of 0.993. HAV KUMC 02-1, HAV KUMC 

The viral load of the four serum and stool samples was quantified by HAV VP3-specific qPCR (Figure 1 ). The copy number of HAV RNAs was determined using a linear regression curve with a coefficient of correlation (r) value of 0.993. HAV KUMC 02-1, HAV KUMC 14-1, and HAV KUMC 14-2 showed Ct values of 27.6-31, corresponding to 10 1 -10 2 viral RNA copies/μL in HAV-positive clinical samples, whereas HAV KUMC 04-1 showed a Ct value of 33, corresponding to <10 copies/μL of viral RNA. The threshold values of the HAV VP4 region were assessed at 29.2-35.3 with <10 2 copies/μL of HAV RNA. 

The coverage of HAV genomic sequences was associated with the viral RNA copy number (Table 2) . Nearly whole-genome sequences of HAV (92.3-94.6%) were recovered from HAV KUMC 14-1, HAV KUMC 14-2, and HAV KUMC 02-1 that contained 10-100 viral RNA copies/μL. The HAV coverage rate was relatively low (86.1%) in HAV KUMC 04-1, as it contained only <10 viral RNA copies/μL. HAV KUMC 98-34 was excluded due to insufficient volume, and the whole-genome sequence was obtained by conventional Sanger sequencing.

Newly obtained HAV genome sequences were compared with representative HAV sequences of different sub-genotypes available in the GenBank (Table 3) . HAV strains in ROK showed nucleotide and amino acid similarities with representative HAV strains of the sub-genotype IA at a rate of 93.8-99.8% and 97.1-99.7%, respectively. 

The coverage of HAV genomic sequences was associated with the viral RNA copy number (Table 2) . Nearly whole-genome sequences of HAV (92.3-94.6%) were recovered from HAV KUMC 14-1, HAV KUMC 14-2, and HAV KUMC 02-1 that contained 10-100 viral RNA copies/µL. The HAV coverage rate was relatively low (86.1%) in HAV KUMC 04-1, as it contained only <10 viral RNA copies/µL. HAV KUMC 98-34 was excluded due to insufficient volume, and the whole-genome sequence was obtained by conventional Sanger sequencing.

Newly obtained HAV genome sequences were compared with representative HAV sequences of different sub-genotypes available in the GenBank (Table 3) . HAV strains in ROK showed nucleotide and amino acid similarities with representative HAV strains of the sub-genotype IA at a rate of 93.8-99.8% and 97.1-99.7%, respectively. 

Partial-genome sequences of the 52 HAV strains were assessed for gene homologies ( Figure S1 ). 

The five HAV strains from ROK demonstrated a geographic cluster with the subgenotype IA strains from Japan, China, and Mongolia ( Figure 2 ). HAV KUMC 98-34 and HAV KUMC 04-1 formed a homologous genetic lineage with HAJFF-Kan12 and HA16-0511 from Japan. HAV KUMC 14-1 was phylogenetically grouped with LU38 from China, forming a distinct monophyletic branch with TD51, HD9, MNA09-B1141, and MNA06-2130 from China and Mongolia. HAV KUMC 02-1 was closely associated with HAV strains from Japan, China, and Mongolia, whereas HAV KUMC 14-2 showed the closest relationship with H2 and MNA10-B1355 from China and Mongolia.

Phylogenetic analysis of the nearly whole-genome hepatitis A virus (HAV) sequences in the Republic of Korea (ROK). Nearly whole-genome sequences of HAV from serum and stool samples were obtained by multiplex polymerase chain reaction-based next-generation sequencing. Branch lengths are proportional to the number of nucleotide substitutions. Vertical distances are included for clarity. Phylogenetic trees of HAV were generated by the maximum likelihood method, with 1000 bootstrap iterations. Colored circles indicate specific HAV strains from ROK; red, genotype IA; blue, genotype IB; orange, genotype IIA; green, genotype IIB; pink, genotype IIIA; sky blue, genotype IIIB; and violet, genotype V.

The individual HAV genes were estimated as the probability of intersecting clades compared with the phylogenetic pattern of whole-genome sequences (Figure 3 and Figure  S2 . Nearly whole-genome sequences of HAV from serum and stool samples were obtained by multiplex polymerase chain reaction-based next-generation sequencing. Branch lengths are proportional to the number of nucleotide substitutions. Vertical distances are included for clarity. Phylogenetic trees of HAV were generated by the maximum likelihood method, with 1000 bootstrap iterations. Colored circles indicate specific HAV strains from ROK; red, genotype IA; blue, genotype IB; orange, genotype IIA; green, genotype IIB; pink, genotype IIIA; sky blue, genotype IIIB; and violet, genotype V.

The individual HAV genes were estimated as the probability of intersecting clades compared with the phylogenetic pattern of whole-genome sequences (Figure 3 and Figure  S2 ). The VP0 (VP4/VP2), VP3, and VP1 regions showed phylogenetic clades with 66.1% Whole-genome and partial-genome phylogenies were generated using the maximum likelihood method. Full lines indicate significant phylogenetic conflicts. Grey lines indicate phylogenetic clade patterns falling between two different regions.

The genomic epidemiology of food-borne viruses plays a critical role in identifying and tracking the source of pathogens during outbreaks [35] [36] [37] [38] [39] [40] [41] [42] . Partial-genome sequences are suitable for investigating the presence of a virus and its genotypic diversity but not for identifying the geographical origin or tracing back the source of an outbreak. In 1990, WGS was applied for analyzing the molecular epidemiology of a prototype norovirus [26] and revealed the periodic emergence of globally spreading variants as well as the origin of food-borne outbreaks [36] . The WGS of hepatitis E virus revealed the transmission source in infected animals or humans, associated with food intake, such as of contaminated shellfish, salad, or vegetables [35, 39] . However, the clinical use of WGS using Sanger sequencing and NGS is limited due to the low levels of viral RNA in human samples [21, [43] [44] [45] [46] . Sanger sequencing requires a markedly long time to yield whole-genome sequences from patient samples with ultra-low viral loads. The application of WGS is often limited by the lack of NGS technology and the high diagnostic test cost. However, NGS has revolutionized the acquisition of whole genomes by generating high-yield genomic data from small-concentration samples. For instance, novel NGS methods, such as the shotgun metagenomics-based workflow and target enrichment, were used for identifying SARS-CoV-2 in clinical samples [23, 24] . The nanopore sequencing method was used to obtain nearly whole-genome sequences of HAV from cultured cells [25] . However, NGSbased WGS has not yet been attempted for detecting HAV in patient samples [37, 46, 47] . Here, multiplex PCR-based NGS was developed to ensure high genome coverage of HAV in clinical samples containing a relatively low viral load. Target-enriched NGS is an essential method for successfully obtaining genome sequences of HAV directly from patient samples. Our results suggested that multiplex PCR-based NGS of HAV clinical samples could help investigate the genotypic diversity and tracking of HAV in global outbreaks. 

The genomic epidemiology of food-borne viruses plays a critical role in identifying and tracking the source of pathogens during outbreaks [35] [36] [37] [38] [39] [40] [41] [42] . Partial-genome sequences are suitable for investigating the presence of a virus and its genotypic diversity but not for identifying the geographical origin or tracing back the source of an outbreak. In 1990, WGS was applied for analyzing the molecular epidemiology of a prototype norovirus [26] and revealed the periodic emergence of globally spreading variants as well as the origin of food-borne outbreaks [36] . The WGS of hepatitis E virus revealed the transmission source in infected animals or humans, associated with food intake, such as of contaminated shellfish, salad, or vegetables [35, 39] . However, the clinical use of WGS using Sanger sequencing and NGS is limited due to the low levels of viral RNA in human samples [21, [43] [44] [45] [46] . Sanger sequencing requires a markedly long time to yield whole-genome sequences from patient samples with ultra-low viral loads. The application of WGS is often limited by the lack of NGS technology and the high diagnostic test cost. However, NGS has revolutionized the acquisition of whole genomes by generating high-yield genomic data from small-concentration samples. For instance, novel NGS methods, such as the shotgun metagenomics-based workflow and target enrichment, were used for identifying SARS-CoV-2 in clinical samples [23, 24] . The nanopore sequencing method was used to obtain nearly whole-genome sequences of HAV from cultured cells [25] . However, NGS-based WGS has not yet been attempted for detecting HAV in patient samples [37, 46, 47] . Here, multiplex PCR-based NGS was developed to ensure high genome coverage of HAV in clinical samples containing a relatively low viral load. Target-enriched NGS is an essential method for successfully obtaining genome sequences of HAV directly from patient samples. Our results suggested that multiplex PCR-based NGS of HAV clinical samples could help investigate the genotypic diversity and tracking of HAV in global outbreaks.

Human-to-human transmission of hepatitis A poses a critical public health threat owing to a high proportion of men who have sex with men (MSM) [48] . In 2015-2020, the widespread outbreaks of hepatitis A among MSM have been reported in Taiwan, Berlin, Netherlands, and United States [49] [50] [51] [52] [53] [54] . Genomic epidemiology of MSM-associated HAV outbreaks allows to detect the infection origin, transmission chains linked to outbreaks or sporadic cases, and populations at high risk. In phylogenetic studies, whole-genome sequences are required to interpret fundamental phylogenetic and genotypic identification of HAV in humans, as well as animals and plants [24, 26, 35, 39, 55] . Whole-genome sequencing of the viruses become increasing overall phylogenetic support, while the phylogeny of partial-genome sequences yielded the incomplete lineage sorting in each genomic region [56] [57] [58] . The construction of phylogenetic conflict signals may be associated with evolutionary biological mechanisms such as recombination, host-switching and co-speciation events [6, [59] [60] [61] [62] [63] [64] [65] . A previous phylogenetic analysis of the hepatitis B virus showed that partial-genome sequences were insufficient for defining phylogeographical links owing to phylogenetic branches with low support, despite identifying new genotypes [56] . The partial-genome with high genomic variability of HAV, such as VP1/P2A and VP1/P2B junctions, has been used to identify and classify HAV genotypes worldwide [38, 41, 42, 46, 57, 66] . Here, the accuracy of genotyping to describe the phylogenetic relationship of HAV strains has not been studied using partial-genome sequences. In this study, the partial genomic sequence of VP1/P2A, VP1/P2B, and entire VP1 regions showed 35.6-94.9% phylogenetic conflicts compared to the whole-genome sequences of HAV. These highly variant regions might be useful for categorizing HAV genotypes but not for identifying reliable phylogenetic relationships or infection sources. In contrast, the partial genomic sequence of HAV VP3 region demonstrated a significantly no phylogenetic conflict (100% phylogenetic parallel patterns) compared with that of the whole-genome sequences, indicating that it might be useful for distinguishing viral genotypes and detecting the source of viral infection in emerging HAV outbreaks. Furthermore, the VP3 gene could be used for the rapid molecular diagnosis of hepatitis A by real-time qPCR, as differences in quantification with the consensus VP4 gene. Therefore, these findings would allow the rapid and accurate molecular detection, genotypic identification, and monitoring of emerging HAV. The limitation of our study is the paucity of clinical samples for various HAV strain comparisons, and thus, further research is needed to characterize the VP3 region and evaluate its accuracy in viral surveillance.

Nearly whole-genome sequences of HAV using multiplex PCR-based NGS facilitate definition of the genetic diversity, molecular epidemiology, and origin of outbreaks. The VP3 gene might serve a representative genetic marker for tracking the infectious origin and molecular diagnosis of emerging HAV outbreaks. Thus, this study provides significant insights into the disease control and preparedness of hepatitis A patients at high risk of severe illness.

The following are available online at https://www.mdpi.com/article/ 10.3390/microorganisms10010100/s1, Table S1 : Primers for whole-genome sequencing of hepatitis A virus. Figure S1 : Comparison of nucleotide similarities within various hepatitis A virus (HAV) genome segments at the genotype level. Figure S2 : Tanglegram of hepatitis A virus (HAV) sequences based on whole-(left) and partial-(right) genome phylogenies. Informed Consent Statement: Informed consent was obtained from all participants involved in the study.

The datasets generated for this study can be found in Genbank, under accession number MW405346-MW405350.

An economic analysis of different strategies of immunization against hepatitis A virus in developed countries

Hepatitis A shifting epidemiology in South-East Asia and China

Genotypic shift of the hepatitis A virus and its clinical impact on acute hepatitis A in Korea: A nationwide multicenter study. Scand

Type A viral hepatitis: A summary and update on the molecular virology, epidemiology, pathogenesis and prevention

Genetic Relatedness of Hepatitis-a Virus-Strains Recovered from Different Geographical Regions

Mattion, N. First full-length genomic sequence of a hepatitis A virus isolated in Argentina shows recombination between subgenotypes IA and IB

Molecular epidemiology of human hepatitis A virus defined by an antigen-capture polymerase chain reaction method

A multistate, foodborne outbreak of hepatitis A. National Hepatitis A Investigation Team

Hepatitis a molecular epidemiology in the United States, 1996-1997: Sources of infection and implications of vaccination policy

Direct sequencing of hepatitis A virus strains isolated during an epidemic in France

Molecular evolution of hepatitis A virus: A new classification based on the complete VP1 protein

Sequence analysis of a new hepatitis A virus naturally infecting cynomolgus macaques (Macaca fascicularis)

Genotypes of hepatitis a virus in Turkey: First report and clinical profile of children infected with sub-genotypes IA and IIIA

Molecular Genotyping of Hepatitis A Virus

Human enteric viruses in Coquina clams associated with a large hepatitis A outbreak

Molecular characterization of hepatitis A virus isolates from a transcontinental shellfish-borne outbreak

Major outbreak of hepatitis A associated with orange juice among tourists

A possible outbreak of hepatitis A associated with semi-dried tomatoes

Update: A food-borne outbreak of hepatitis A in the Netherlands related to semi-dried tomatoes in oil

A Multistate Outbreak of Hepatitis A Associated With Semidried Tomatoes in Australia

Implementation of Next-Generation Sequencing for Hepatitis B Virus Resistance Testing and Genotyping in a Clinical Microbiology Laboratory

Next generation sequencing technology: Advances and applications

Clinical Metagenomic Next-Generation Sequencing for Pathogen Detection

Multiple approaches for massively parallel sequencing of SARS-CoV-2 genomes directly from clinical samples

Whole Genome Sequencing of Hepatitis A Virus Using a PCR-Free Single-Molecule Nanopore Sequencing Approach

Sequence and genomic organization of Norwalk virus

Phylogeographic analysis of hemorrhagic fever with renal syndrome patients using multiplex PCR-based next generation sequencing

Multiplex PCR-Based Next-Generation Sequencing and Global Diversity of Seoul Virus in Humans and Rats

Genomic surveillance reveals multiple introductions of SARS-CoV-2 into Northern California

Isolation and Genomic Sequence of Hepatitis A Virus from Mixed Frozen Berries in Italy

Design and Application of Nucleic Acid Standards for Quantitative Detection of Enteric Viruses by Real-Time PCR

Hepatitis A virus subgenotyping based on RT-qPCR assays

MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets

An R package for visualizing, adjusting and comparing trees of hierarchical clustering

Complete genome analysis of avian hepatitis E virus from chicken with hepatic rupture hemorrhage syndrome

A multi-site study of norovirus molecular epidemiology in Australia

An Outbreak of Acute Hepatitis Caused by Genotype IB Hepatitis A Viruses Contaminating the Water Supply in Thailand

High prevalence of antibodies to hepatitis A and E viruses and viremia of hepatitis B, C, and D viruses among apparently healthy populations in Mongolia

Analysis of complete genome sequences of swine hepatitis E virus and possible risk factors for transmission of HEV to humans in Korea

Hepatitis A outbreak associated with a revolving sushi bar in Chiba, Japan: Application of molecular epidemiology

Genetic diversity of hepatitis A virus in China: VP3-VP1-2A genes and evidence of quasispecies distribution in the isolates

Evolution and genetic characterization of hepatitis A virus isolates in China

Comparison of targeted next-generation sequencing for whole-genome sequencing of Hantaan orthohantavirus in Apodemus agrarius lung tissues

Analysis of frozen strawberries involved in a large norovirus gastroenteritis outbreak using next generation sequencing and digital PCR

Phylogeographic diversity and hybrid zone of Hantaan orthohantavirus collected in Gangwon Province

Full-length genomic sequence of subgenotype IIIA hepatitis A virus isolate in Republic of Korea

Full-length sequences of subgenotype IIIA and IIIB hepatitis A virus isolates: Characterization of genotype IIIHAV genornes

Public Health Responses to Person-to-Person Hepatitis A Outbreaks

Hepatitis A Virus Infections-United States

Hepatitis A outbreak predominantly affecting men who have sex with men in Northern Ireland

Hepatitis A outbreak among men who have sex with men (MSM) predominantly linked with the EuroPride, the Netherlands

Hepatitis A Virus Infections Among Men Who Have Sex with Men-Eight US States

Hepatitis A outbreak among MSM in Berlin due to low vaccination coverage: Epidemiology, management, and successful interventions

Hepatitis A Outbreak Among Men Who Have Sex With Men in a Country of Low Endemicity of Hepatitis A Infection

Near-Complete Genome Sequence of a Hepatitis A Subgenotype IB Virus Isolated from Frozen Raspberries

A New Method for Next-Generation Sequencing of the Full Hepatitis B Virus Genome from A Clinical Specimen: Impact for Virus Genotyping

Analysis of the full-length genome of a subgenotype IIIB hepatitis a virus isolate: Primers for broadly reactive PCR and genotypic analysis

Full-length genome characterization and genetic relatedness analysis of hepatitis A virus outbreak strains associated with acute liver failure among children

Evidence of recombination in natural populations of hepatitis A virus

Recombination in hepatitis A virus: Evidence for reproductive isolation of genotypes

Identification of recombination between subgenotypes IA and IB of hepatitis A virus

Bats, Primates, and the Evolutionary Origins and Diversification of Mammalian Gammaherpesviruses

Gene trees in species trees

Deciphering the Origin and Evolution of Hepatitis B Viruses by Means of a Family of Non-enveloped Fish Viruses

We thank Man-Seong Park for providing an Illumina MiSeq instrument at Korea University.

The authors declare no conflict of interest.