key: cord-1041906-fe136t5f
authors: Park, Kyungmin; Lee, Seung-Ho; Kim, Jongwoo; Lee, Jingyeong; Lee, Geum-Young; Cho, Seungchan; Lee, Seung Ho; Park, Kkothanahreum; No, Jin Sun; Budhathoki, Shailesh; Kim, Yu-Jin; Kim, Young-Su; Kim, Heung-Chul; Klein, Terry A.; Kim, Won-Keun; Song, Jin-Won
title: Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea
date: 2021-05-06
journal: Viruses
DOI: 10.3390/v13050847
sha: 01f65a3fd099d657e66094c0f85126818097af71
doc_id: 1041906
cord_uid: fe136t5f

Whole-genome sequencing of infectious agents enables the identification and characterization of emerging viruses. The MinION device is a portable sequencer that allows real-time sequencing in fields or hospitals. Hantaan orthohantavirus (Hantaan virus, HTNV), harbored by Apodemus agrarius, causes hemorrhagic fever with renal syndrome (HFRS) and poses a critical public health threat worldwide. In this study, we aimed to evaluate the feasibility of using nanopore sequencing for whole-genome sequencing of HTNV from samples having different viral copy numbers. Amplicon-based next-generation sequencing was performed in A. agrarius lung tissues collected from the Republic of Korea. Genomic sequences of HTNV were analyzed based on the viral RNA copy numbers. Amplicon-based nanopore sequencing provided nearly full-length genomic sequences of HTNV and showed sufficient read depth for phylogenetic analysis after 8 h of sequencing. The average identity of the HTNV genome sequences for the nanopore sequencer compared to those of generated from Illumina MiSeq revealed 99.8% (L and M segments) and 99.7% (S segment) identities, respectively. This study highlights the potential of the portable nanopore sequencer for rapid generation of accurate genomic sequences of HTNV for quicker decision making in point-of-care testing of HFRS patients during a hantavirus outbreak.

Orthohantaviruses (family, Hantaviridae; order, Bunyavirales) are enveloped negativesense single-stranded RNA viruses [1] . The viral genome consists of tripartite RNA segments (large, medium, and small) encoding an RNA-dependent RNA polymerase (RdRp), two surface glycoproteins (Gn and Gc), and a nucleocapsid (N) protein, respectively [2] .

Rodents and shrews were handled in accordance with the ethical guidelines of the Korea University Institutional Animal Care and Use Committee (KUIACUC #2016-49). Small mammals were euthanized via cardiac puncture under isoflurane anesthesia. Necropsy and experiments were performed in a biosafety level 3 laboratory at Korea University.

Field trappings were performed using Sherman live traps (8 cm × 9 cm × 23 cm; H. B. Sherman, Tallahassee, FL, USA) in Gyeonggi and Gangwon Provinces, ROK, in 2019. For each day, a total of 100 traps was placed in unmanaged grasses, herbaceous vegetations, farmlands, mountains, fields, and forests at intervals of 1-5 m for 1-3 days. The captured rodents and shrews were sequentially numbered, placed in secure containers, and transported to the biosafety level 3 facility at Korea University in accordance with the safety guidelines. Sera and multiple tissue samples were aseptically obtained from the Viruses 2021, 13, 847 3 of 16 small mammals. The sera and heart fluids were separated by centrifugation for 5 min at 4 • C. The tissue samples were stored at −80 • C until use.

Total DNA was extracted from the liver tissues using a High Pure PCR Template Preparation Kit (Roche, Basel, Switzerland). Mitochondrial cytochrome b (Cyt b) genespecific PCR was performed using universal primers forward: 5 -CGA AGC TTG ATA TGA AAA ACC ATC GTT G-3 and reverse: 5 -CTG GTT TAC AAG ACC AGA GTA AT-3 [28] . Polymerase chain reaction (PCR) was performed in 25 µL reaction mixtures, containing 2.5 mM dNTP (ELPIS Biotech, Daejeon, South Korea) and 1 U of Super-Therm Taq polymerase (JMR Holdings, London, UK). First cycling conditions consisted of an initial denaturation at 95 • C for 5 min followed by 8 cycles with denaturation at 94 • C for 1 min, annealing at 50 • C for 1 min, and elongation at 72 • C for 2 min. Second cycling was conducted by 30 cycles of denaturation at 94 • C for 1 min, annealing at 55 • C for 1 min, and elongation at 72 • C for 2 min, and final cycle at 72 • C for 10 min in a thermal cycler (ProFlex PCR System, Life Technology, CA, USA). The obtained sequences were deposited in GenBank (accession numbers: MW796110-MW796121).

The sera or heart fluids were diluted 1:32 and 1:2, respectively, and placed in duplicate wells of acetone-fixed Vero E6 cells infected with HTNV. The plate was incubated at 37 • C for 30 min, then the cells were washed with Phosphate-Buffered Saline (PBS) and distilled water (D.W.). The slides were treated with fluorescein isothiocyanate-conjugated antimouse (for rodents) or mouse/rat (for shrews) immunoglobulin G (IgG) antibodies (MP Bio, CA, USA). Subsequently, the slides were incubated at 37 • C for 30 min, washed with PBS and D.W., and examined under a fluorescence microscope (Axio Scope; Zeiss, Berlin, Germany). Virus-specific fluorescence is indicative of HTNV infection.

Total RNA was extracted from the lung tissues of the rodents using TRI Reagent Solution (Ambion, Austin, Texas, USA) according to the manufacturer's instructions. cDNA was synthesized from 1 µg of total RNA by using the High Capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA, USA) with a random hexamer and OSM55 (5 -TAG TAG ACT CC-3 ). PCR conditions and primer sequences for amplification of HTNV genomes have been previously described [29] . Oligonucleotide primer sequences followed up Han-L-F1 (outer): 5 -ATG TAY GTB AGT GCW GAT GC-3 and Han-L-R1 (outer):  5 -AAC CAD TCW GTY CCR TCA TC-3 First-and second-round PCR were performed in 25 µL reaction mixture containing 2.5 mM dNTP (ELPIS Biotech, Daejeon, South Korea), 1U of Super-Therm Taq DNA polymerase (JMR Holdings, London, UK), and 1.5 µg of cDNA with 10 pM of each primer. Cycling condition consisted of initial denaturation at 94 • C for 5 min, followed by 6 cycles of denaturation at 94 • C for 30 sec, annealing at 37 • C for 40 sec, and elongation at 72 • C for 1 min. Second cycling was performed by 32 cycles of denaturation at 94 • C for 30 sec, annealing at 42 • C for 40 sec, and elongation at 72 • C for 1 min, and final cycle at 72 • C for 5 min (ProFlex PCR System, Life Technology, CA, USA). PCR products were purified using MinElute PCR purification kit (Qiagen, Hilden, Germany). Amplicons were sequenced in forward and reversed directions of each PCR 

Quantitative PCR (qPCR) was performed using SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) on a Quantstudio 5 Flex Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) using the following reaction steps: denaturation at 95 • C for 10 min, followed by 40 cycles at 95 • C for 15 s and 60 • C for 1 min. The forward (5 -TTA  TTG TGC TCT TCA TGG TTG C-3 ) and reverse primer (5 -CAT CCC CTA AGT GGA AGT TGT C-3 ) sequences of the HTNV S segment were used [30] . The viral copy number of the HTNV genomes was calculated using a linear regression curve described in a previous study [20] .

cDNA was amplified using HTNV-specific primer mixtures using Solg 2X Uh-Taq PCR Smart Mix (Solgent, Daejeon, Republic of Korea) according to the manufacturer's instructions. The PCR conditions and primer information were described in a previous study [18] . The enrichment was performed in 25 µL reaction mixture containing 12.5 µL 2X Uh premix, 2.0 µL of each primer mixture, 10.5 µL of D.W., and 1.0 µL of DNA template. Multiplex PCR was performed as follows: a cycle at 95 • C for 15 min, followed by 40 cycles or 25 cycles at 95 • C for 20 s, 50 • C for 40 s, and 72 • C for 1 min, and a cycle at 72 • C for 3 min.

The library was prepared using Ligation Sequencing Kit (SQK-LSK109) and Native Barcoding Kit (EXP-NBD104 and NBD 114) according to the standard protocols (Oxford Nanopore Technologies, London, UK). Amplified libraries were purified using AMPure XP beads (Beckman Coulter, CA, USA). Barcoded libraries were pooled, ligated to sequencing adapters, and sequenced using the MinION by using FLO-MIN106 (R9.4) flow cell for 16 h. Real-time basecalling was performed using Guppy (v3.0.3; MinKNOW software) in the MK1C device. Raw data were demultiplexed, and the adaptor sequences were trimmed. Only reads with a default Q-score of 7 were included in the subsequent analysis to improve the reliability of the data. The filtered sequences were combined into a single FASTA file for each barcode by using Porechop (v.9.0). The reads were mapped to the reference genome sequence of HTNV 76-118, and the consensus sequences were extracted using the CLC Genomics Workbench (v7.5.2; Qiagen, Hilden, Germany) [20] . The consensus sequences were manually polished using the error-correction method. An indel, insertion, or deletion in the open reading frame (ORF) was considered as a nanopore sequencing error. If the indels occurred only in the consensus sequence, the correction was performed based on the reference sequences. The uncovered sequences of HTNV tripartite genomes were determined using conventional nested RT-PCR. The workflow and time-span overview of the nanopore sequencing used in this study are shown in Figure 1 . The genomic sequences of HTNV are described in Table S1.

DNA libraries were prepared using the TruSeq Nano DNA LT sample preparation kit (Illumina, San Diego, USA) according to the manufacturer's instructions and a previously described method [18] . NGS was performed on a MiSeq benchtop sequencer platform (Illumina, San Diego, USA) with 2 × 150 bp by using a MiSeq reagent kit v2 (Illumina, San Diego, USA). Total reads were trimmed to the adaptor and index sequences, and low-quality reads were filtered using CLC Genomics Workbench version 7.5.2 (CLC Bio, Cambridge, MA, USA). The filtered reads were mapped to the tripartite genome sequence of HTNV 76-118, and consensus sequences were extracted. The uncovered genomic sequences Table S1 .

Viruses 2021, 13, x FOR PEER REVIEW 5 of 19 Figure 1 . Workflow overview of the amplicon-based nanopore sequencing for Hantaan virus (HTNV) used in this study. Total RNA was extracted from rodent lung tissues, and cDNA was synthesized from 1 μg of total RNA with a random hexamer and OSM55. cDNA was amplified using HTNV-specific primer mixtures for library preparation. Barcoded libraries were pooled, ligated to sequencing adapter, and sequenced using the MinION device with FLO-MIN106 (R9.4) flow cell for 8 h. The basecalled sequences were combined into a single FASTA file, and the reads were mapped to the reference genome sequences. The consensus sequences were extracted for downstream analysis.

DNA libraries were prepared using the TruSeq Nano DNA LT sample preparation kit (Illumina, San Diego, USA) according to the manufacturer's instructions and a previously described method [18] . NGS was performed on a MiSeq benchtop sequencer platform (Illumina, San Diego, USA) with 2 × 150 bp by using a MiSeq reagent kit v2 (Illumina, San Diego, USA). Total reads were trimmed to the adaptor and index sequences, and lowquality reads were filtered using CLC Genomics Workbench version 7.5.2 (CLC Bio, Figure 1 . Workflow overview of the amplicon-based nanopore sequencing for Hantaan virus (HTNV) used in this study. Total RNA was extracted from rodent lung tissues, and cDNA was synthesized from 1 µg of total RNA with a random hexamer and OSM55. cDNA was amplified using HTNVspecific primer mixtures for library preparation. Barcoded libraries were pooled, ligated to sequencing adapter, and sequenced using the MinION device with FLO-MIN106 (R9.4) flow cell for 8 h. The basecalled sequences were combined into a single FASTA file, and the reads were mapped to the reference genome sequences. The consensus sequences were extracted for downstream analysis.

The genomic sequences of the HTNV L, M, and S segments were aligned using the Clustal W algorithm in Lasergene version 5 (DNASTAR Inc., Madison, WI, USA). Phylogenetic analysis was performed using the best-fit TN93+G+I (for L segment) and T92+G+I (for M and S segments) substitution models of evolution by using the maximum-Viruses 2021, 13, 847 6 of 16 likelihood (ML) method in MEGA 7.0 [31] . Topologies were evaluated using bootstrap analysis of 1000 iterations.

Small mammals were collected from multiple sites, including Paju, Pocheon, Pyeongtaek, Yeoncheon, Cheorwon, Chuncheon, and Hwacheon, in Gyeonggi and Gangwon Provinces, ROK, in 2019 ( Figure 2 ). Geographic information for the sites is shown in Table S2 . The captured animals included 303 A. agrarius, 13 Crocidura lasiura, 3 C. shantungensis, 17 Mus musculus, four Myodes regulus, 5 Micromys minutus, 1 Microtus fortis, 3 Rattus norvegicus, and 1 Tscherskia triton (Table 1) . A. agrarius (86.6% of all individuals) was the most frequently collected species at all trapping sites. Phylogenetic analysis of the mitochondrial Cyt b gene showed that the captured rodents formed a genetic lineage with A. agrarius ( Figure S1 ).

quences of HTNV L, M, and S segments were determined by conventional nested RT-PCR. The genomic sequences of HTNV are described in Table S1 .

The genomic sequences of the HTNV L, M, and S segments were aligned using the Clustal W algorithm in Lasergene version 5 (DNASTAR Inc. WI, USA). Phylogenetic analysis was performed using the best-fit TN93+G+I (for L segment) and T92+G+I (for M and S segments) substitution models of evolution by using the maximum-likelihood (ML) method in MEGA 7.0 [31] . Topologies were evaluated using bootstrap analysis of 1000 iterations.

Small mammals were collected from multiple sites, including Paju, Pocheon, Pyeongtaek, Yeoncheon, Cheorwon, Chuncheon, and Hwacheon, in Gyeonggi and Gangwon Provinces, ROK, in 2019 ( Figure 2 ). Geographic information for the sites is shown in Figure S1 ). The colored circles indicate the HTNV RNA positive sites: Pocheon, blue (Jangam-ri); Yeoncheon, green (Dosin-ri) and brown (Sang-ri); Cheorwon, orange (Gwanu-ri) and red (Dochang-ri); Hwacheon-gun, violet (Pungsan-ri). The black circles represent the locations where no HTNV RNA The colored circles indicate the HTNV RNA positive sites: Pocheon, blue (Jangam-ri); Yeoncheon, green (Dosin-ri) and brown (Sang-ri); Cheorwon, orange (Gwanu-ri) and red (Dochang-ri); Hwacheon-gun, violet (Pungsan-ri). The black circles represent the locations where no HTNV RNA was detected: Paju (Jeogam-ri); Pocheon (Jail-ri); Yeoncheon (Hyunga-ri, Dongi-ri, Ohkye-ri, Yangwon-ri, and Bugok-ri); Cheorwon (Eupnae-ri, Munhye-ri, and Chungyang-ri); and Chuncheon (Gamjeong-ri, Sinchon-ri, Jinae-ri, Balsan-ri, and Yulmun-ri). A Quantum Geographical Information System (QGIS) 3.10 for Mac was used to create the map, which was modified in Adobe Illustrator CC 2019. 

A total of 20/303 (6.6%) rodents were seropositive, of which 1/18 (5.6%), 5/48 (10.4%), 7/105 (6.7%), 1/38 (2.6%), and 6/46 (13.0%) were detected in Pocheon, Yeoncheon, Cheorwon, Chuncheon, and Hwacheon, respectively; no seropositive rodents were detected in Paju and Pyeongtaek ( Table 2 ). The seropositive samples were examined for hantaviral RNA by using HTNV-specific RT-PCR. In total, 12/20 (60%) seropositive A. agrarius harbored HTNV RNA, including 1/1 (100%) in Pocheon, 2/5 (40%) in Yeoncheon, 6/7 (85.7%) in Cheorwon, and 3/6 (50%) in Hwacheon. Hantaviral antibodies were not detected in the captured rodents and shrews except for A. agrarius. 

Viral RNA load was measured in the lung tissues of the captured rodents (Aa19-36, Aa19-38, Aa19-39, Aa19-57, Aa19-89, Aa19-144, Aa19-152, Aa19-153, Aa19-167, Aa19-233, Aa19-236, and Aa19-278; Table 3 ). Aa19-233 and Aa19-278 harbored the highest amount of HTNV RNA (Ct value, 20.5-22.2), corresponding to 10 4 to 10 5 copies/µL in the evaluated rodents. Aa19-36 and Aa19-236 showed Ct values ranging from 29.4 to 29.7 (10 2 to 10 3 copies/µL of viral RNA). The Ct values of the rodents Aa19-38, Aa19-57, Aa19-89, Aa19-153, and Aa19-167 ranged from 31.8 to 34.5, indicating that the samples contained 10 1 to 10 2 copies/µL of viral RNA. Aa19-144 and Aa19-152 consisted of 1 to 10 1 copies/µL of HTNV RNA with Ct values ranging from 36.4 to 37.3. In Aa19-39, HTNV RNA was not detectable by qPCR. 

In the MinION sequencing, 12 nearly whole-genome sequences of HTNV were obtained from the lung tissues of A. agrarius captured from Gyeonggi and Gangwon Provinces, ROK. The average depth of coverage for HTNV genomes was calculated using the mapped reads corresponding to sequencing running times (2, 4, 8, and 16 h; Figure 3 ). After 8 h, MinION sequencing produced nearly complete coverage of the HTNV tripartite genomes with over 1000× average read depth for each segment. The sequencing of Aa19-39 showed a low depth of coverage for the HTNV L, M, and S segments in all the tested specimens owing to the lowest viral load. The genome coverage of HTNV was 65.4-96.8% for the L segment, 90.6-98.5% for the M segment, and 92.1-100% for the S segment ( Table 3) . The uncovered genome sequences of the HTNV L, M, and S segments were determined using conventional RT-PCR. The 3 and 5 terminal sequences were empirically replaced owing to the conserved region of the family Hantaviridae. Raw reads obtained after 8 h of sequencing on the MinION platform were extracted and evaluated for indel (insertion and deletion) and mismatch sequences ( Table 4 ). The mean number of single base indels in the homopolymer regions was 3.8 for L segment, 4.4 for M segment, and 1 for S segment. Insertion error was not detected in any of the samples. The total accuracy rates of the HTNV genome sequences obtained using the MinION platform compared to those generated by MiSeq were 99.1-99.9% for the L segment, 99.0-99.9% for the M segment, and 98.3-100% for the S segment. Notably, nucleotide mismatch was the most frequent in the three observed error patterns.

All amplified samples were sequenced in the MiSeq platform to confirm the consensus accuracy of the nanopore sequencing. The coverage of the genome sequences of 12 HTNV strains was 98.2-100% for the L segment, 93.8-100% for the M segment, and 93.8-100% for the S segment (Table S3 ). The average number of mapped reads and the depth of coverage are shown in Table S4 . The uncovered sequences of the HTNV tripartite genomes were determined using conventional RT-PCR. The genomic sequences of the 3 and 5 ends were empirically filled with the conserved sequence of the family Hantaviridae.

All amplified samples were sequenced in the MiSeq platform to confirm the consensus accuracy of the nanopore sequencing. The coverage of the genome sequences of 12 HTNV strains was 98.2-100% for the L segment, 93.8-100% for the M segment, and 93.8-100% for the S segment (Table S3 ). The average number of mapped reads and the depth of coverage are shown in Table S4 . The uncovered sequences of the HTNV tripartite genomes were determined using conventional RT-PCR. The genomic sequences of the 3 and 5 ends were empirically filled with the conserved sequence of the family Hantaviridae. 

A total of 12 full-length genomic sequences of HTNV were obtained from multiple regions of Gyeonggi and Gangwon Provinces, ROK. Seven whole-genome sequences of HTNV were recovered from new trapping sites of HFRS endemic areas, including Yeoncheon (Sang-ri), Cheorwon (Dochang-ri), and Hwacheon (Pungsan-ri). The phylogenies of the HTNV L, M, and S segments from Aa19-36, Aa19-38, and Aa19-39 in Cheorwon (Dochang-ri) formed a genetic cluster with the strains collected from Cheorwon (Munhyeri) (Figure 4) . The phylogenetic pattern of the HTNV tripartite genomes from Aa19-57 shared a common ancestor with the HTNV collected in Pocheon. Phylogenetically, the L and M segments of HTNV from Aa19-89 clustered with the HTNV from Cheorwon (Gwanu-ri), whereas the S segment was grouped with the HTNV from Pocheon. The tripartite segments of Aa19-144, Aa19-152, and Aa19-153 in Hwacheon (Pungsan-ri) showed a genetic lineage with HTNV in Hwacheon (Sanyang-ri). The L, M, and S segments of HTNV from Aa19-167 in Yeoncheon (Sang-ri) were phylogenetically grouped with HTNV in Yeoncheon (Tonghyeon-ri). The tripartite HTNV genomes of Aa19-233, Aa19-236, and Aa19-278 from Cheorwon (Gwanu-ri) formed a homologous genetic group with HTNV in Gwanu-ri, Cheorwon. 

Rodent-borne orthohantaviruses pose a critical public health threat worldwide [32] . In the ROK, approximately 400 HFRS cases are reported annually, with a mean case fatality rate of 1-4% [33] . HFRS is highly endemic in Gyeonggi and Gangwon Provinces, ROK, and affects both military personnel and civilians [29, [34] [35] [36] . Phylogenetic analyses of hantaviruses delineated a clear molecular epidemiological relationship between HFRS patients and their reservoir hosts, indicating the most putative site of infection [37] . In 2005, HTNV genomic sequences obtained from four US soldier patients diagnosed with HFRS in the ROK showed a phylogenetic link with the viral sequences acquired by targeted rodent trapping in military training locations where the patients had exercised. In our previous studies, multiplex PCR-based NGS was used to obtain nearly whole-genome sequences of HTNV from four ROK and US soldiers during 2013-2015 [18] . Recently, active surveillance with targeted small mammal trapping in the assumed regions of hantavirus outbreaks associated with HFRS enabled the phylogeographical analyses of HTNV within a short distance (approximately 5 km) to track the precise infectious agents and potential exposed sites [38] . Establishing a database of viral genomic sequences to clarify the epidemiological association of patients with infectious sources at a higher resolution can contribute to the development of control strategies of hantavirus outbreaks to prevent endemic transmission [39] . In this study, 12 whole-genome sequences of HTNV newly recovered from A. agrarius collected from Gyeonggi and Gangwon Provinces enhanced the resolution of the phylogeographic map of orthohantaviruses. The phylogeny of the HTNV tripartite genomes showed well-supported phylogeographic lineages and genetic diversity in the restricted regions. Furthermore, seven full-length genomic sequences of HTNV were acquired in new trapping sites, including Yeoncheon (Sang-ri), Cheorwon (Dochang-ri), and Hwacheon (Pungsan-ri). Our findings highlight that disease risk assessment and vigilance should be conducted to prevent potential hantavirus outbreaks to humans in HFRS endemic areas.

High-throughput sequencing technology offers a robust platform to investigate the genetic diversity across viral populations, which is essential for understanding viral evolution and pathogenesis [40] [41] [42] [43] . The accuracy of consensus sequences depends on the background error profiling of sequencing technology and the required level of sensitivity [44] . A high level of sequencing coverage is required for population level characterization [45] . Wang et al. revealed that the identification of minor variants at the virus population level required a minimum depth of 400× and 1000× coverage to recognize variants found at a 1% and 0.5% SNP frequency with 99.999% confidence owing to the polymerase-based error and need for amplification, respectively [46] . The targeted amplification of the viral genomes is necessary to achieve these requirements. Three distinct approaches, including SISPA, RNA access, and amplicon-based methods, have been reported for whole-genome sequencing of HTNV from patient and rodent specimens [18] [19] [20] . Of these methods, amplicon-based NGS demonstrated the most sensitive approach to attain nearly complete genome coverage of HTNV and provided the most read depth by using only 10 2 viral RNA copies/µL from the lung tissues of wild reservoir rodents [20] . In this study, amplicon-based nanopore sequencing showed nearly full genome coverage of HTNV except for Aa19-39, with an average depth of 1321× for the L segment; 3029× for the M segment; and 3766× for the S segment, after 8 h of sequencing. At the consensus level, the mean accuracy of the HTNV genome sequences using the MinION platform demonstrated 99.8 (L and M segments) and 99.7% (S segment) identities compared to those using Illumina MiSeq. Our data suggest that the amplicon-based method is the most suitable for obtaining high-quality whole-genome sequences of HTNV in the nanopore platform.

The accuracy of genome sequences generated from raw reads plays an important role in metadata analysis for epidemiological surveillance [47] . Nanopore technology provides low-quality reads compared to the Illumina MiSeq platform, which offers 90% of bases with a basecalling accuracy higher than a Q-score of 30 (99.9%) [48] . Currently, the MinION device has error rates of approximately 15% in single reads by using the flow cell (R9 version) [49] [50] [51] [52] [53] . Nanopore-based sequencing led to the generation of indel errors that cannot be resolved naturally in homopolymer regions, even with high read depth and coverage [52, 54, 55] . Single nucleotide indels in the coding regions might induce frameshift mutations. The indel errors in homopolymer sites are a significant limitation to the reliability of sequence accuracy in nanopore technologies. In this study, manual correction of errors at homopolymer sites was performed in the ORF regions to obtain the highest quality consensus sequences of HTNV tripartite genomes. The numbers of nucleotide deletions in HTNV genomes were 2-5 for L segment, 2-6 for M segment, and 1 for S segment. The proteins generated from each segment of HTNV genome are essential for viral function [56, 57] . The likelihood of a frameshift mutation naturally occurring in the coding regions of proteins in a viable virus is extremely low. Given these observations, we suggest that manual correction should be appropriately performed in homopolymer sites to resolve indel errors within the coding regions of the HTNV tripartite genomes.

The rapid and accurate generation of full-length viral sequence data is important for improving the approaches to clinical diagnosis (including point-of-care diagnostics), patient-stratified management, epidemiological surveillance, and long-term development of cure strategies [58] . Although the Illumina MiSeq platforms can produce massive and high-quality data, their application for clinical diagnosis and public health has been limited by sequencer size, slowness, and library preparation complexity [59] . In contrast, the Min-ION device, a compact portable sequencer, offers the application of real-time sequencing in point-of-care tests for infectious diseases, including identification of CHIKV, EBOV, and ZIKV from clinical blood samples [26, 59, 60] . The major advantages of the MinION system are the reasonable price, the ability to produce real-time data, and portability, which would make it an appropriate device to be used in a hospital or field for rapid diagnosis [24] . In previous studies, a genomic surveillance system based on phylogenetic analyses provided an effective prevention strategy against emerging orthohantavirus outbreaks for mitigating HFRS incidences, such as the restriction of field activities and environmental cleaning to eliminate rodents in endemic areas [37, 39] . These surveillance systems depend on the rapid and accurate sequencing of infectious agents from clinical samples. In this study, we showed that amplicon-based nanopore sequencing provides nearly whole-genome sequences for the identification of the phylogeographic diversity of HTNV within 8 h of sequencing. Our study indicates that amplicon-based sequencing with the MinION sequencer can offer the potential for quicker diagnosis and facilitation of rapid decision making in point-of-care testing of HFRS patients in clinical or field situations.

In conclusion, we report 12 complete HTNV sequences newly obtained from A. agrarius collected from HFRS-endemic areas of Gyeonggi and Gangwon Provinces, ROK. These results improve the phylogeographic maps for monitoring and responding to hantaviral outbreaks in endemic regions. We found that amplicon-based nanopore sequencing provides rapid and accurate genomic sequences of HTNV for identifying infectious sources and their genetic diversity at the strain level. Our findings suggest the potential for the application of MinION sequencing for the rapid analyses of HTNV from clinical and animal specimens in fields during an endemic outbreak. This study provides significant insights into nanopore sequencing for genomic-based diagnosis and control strategies for hantavirus outbreaks in highly HFRS-endemic areas.

Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/v13050847/s1. Figure S1 , Phylogenetic analysis based on the mitochondrial cytochrome b gene of Apodemus spp; Table S1, GenBank accession numbers of Hantaan virus sequences obtained in this study; Table S2 , Topography and GPS coordinates of the trapping sites; Table S3 , Results of MiSeq sequencing from rodent lung tissues by viral RNA copy number of Hantaan virus (HTNV) collected in Gyeonggi and Gangwon Provinces, ROK, 2019; Table S4 , Summary of mapped reads and depth of coverages for Hantaan virus by amplicon-based NGS in MiSeq platform..

Hantavirus infections in Europe and their impact on public health

Current Classification and Future Perspectives. Viruses

Spectrum of hantavirus infection: Hemorrhagic fever with renal syndrome and hantavirus pulmonary syndrome

Another case of "European hantavirus pulmonary syndrome" with severe lung, prior to kidney, involvement, and diagnosed by viral inclusions in lung macrophages

Hantaviruses: Rediscovery and new beginnings

Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents

Reconstructing the evolutionary origins and phylogeography of hantaviruses

Hantavirus maintenance and transmission in reservoir host populations

Hemorrhagic fever with renal syndrome: Mode of transmission to humans

Application of next-generation sequencing technologies in virology

Rapid identification of known and new RNA viruses from animal tissues

Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: Phylogenetic analyses and markers for zoonotic potential

Emergence of a novel subclade of influenza A (H3N2) virus in London

Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic

Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak

Massively parallel sequencing for monitoring genetic consistency and quality control of live viral vaccines

Genomic monitoring to understand the emergence and spread of Usutu virus in the Netherlands

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

Sequence-Independent, Single-Primer Amplification Next-Generation Sequencing of Hantaan Virus Cell Culture-Based Isolates

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

Using small RNA deep sequencing data to detect human viruses

Illumina and Nanopore methods for whole genome sequencing of hepatitis B virus (HBV)

Rissanen, I. Zika virus in the Americas: Early epidemiological and genetic findings

Portable Rabies Virus Sequencing in Canine Rabies Endemic Countries Using the Oxford Nanopore MinION

Rapid, sensitive, full-genome sequencing of severe acute respiratory syndrome coronavirus 2. Emerg

Real-time, portable genome sequencing for Ebola surveillance

Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis

Evolution of the cytochrome b gene of mammals

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

Detection of Hantaan virus RNA from anti-Hantaan virus IgG seronegative rodents in an area of high endemicity in Republic of Korea

MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets

Korea Disease Control and Prevention Agency (KDCA)

Genetic Diversity and Reassortment of Hantaan Virus Tripartite RNA Genomes in Nature, the Republic of Korea

Prevalence of hantavirus infection in wild rodents from five provinces in Korea

The influence of climatic factors on the development of hemorrhagic fever with renal syndrome and leptospirosis during the peak season in Korea: An ecologic study

Hemorrhagic fever with renal syndrome in 4 US soldiers

Active Targeted Surveillance to Identify Sites of Emergence of Hantavirus

Genomic Epidemiology and Active Surveillance to Investigate Outbreaks of Hantaviruses

Viruses as quasispecies: Biological implications

Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population

Towards a genomics-informed, real-time, global pathogen surveillance system

Middle East respiratory syndrome: What we learned from the 2015 outbreak in the Republic of Korea

Standards for sequencing viral genomes in the era of high-throughput sequencing

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data

Characterization of mutation spectra with ultra-deep pyrosequencing: Application to HIV-1 drug resistance

Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands

Illumina Technical Note: Sequencing. Quality scores for Next-Generation Sequencing. Assessing sequencing accuracy using Phred quality scoring

MinION Analysis and Reference Consortium: Phase 2 data release and analysis of R9.0 chemistry

MinION-based long-read sequencing and assembly extends the Caenorhabditis elegans reference genome

Mapping and phasing of structural variation in patient genomes using nanopore sequencing

Performance of neural network basecalling tools for Oxford Nanopore sequencing

Nanopore sequencing and assembly of a human genome with ultra-long reads

A complete bacterial genome assembled de novo using only nanopore sequencing data

Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation

Hantaviral Proteins: Structure, Functions, and Role in Hantavirus Infection. Front Microbiol

Sin Nombre virus glycoprotein trafficking

Tracking virus outbreaks in the twenty-first century

Detection of viral pathogens with multiplex nanopore MinION sequencing: Be careful with cross-Talk

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

We thank Man-Seong Park in Korea University for supporting the experiments.

The authors declare no conflict of interest.