Next Article in Journal
Convalescent Plasma Transfusion for the Treatment of COVID-19 in Adults: A Global Perspective
Next Article in Special Issue
Continuing Orthohantavirus Circulation in Deer Mice in Western Montana
Previous Article in Journal
In Vivo Models to Study the Pathogenesis of Extra-Respiratory Complications of Influenza A Virus Infection
Previous Article in Special Issue
Identification of Novel Antiviral Compounds Targeting Entry of Hantaviruses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea

1
BK21 Graduate Program, Department of Biomedical Sciences, Korea University College of Medicine, Seoul 02841, Korea
2
Department of Microbiology, Korea University College of Medicine, Seoul 02841, Korea
3
Division of High-Risk Pathogens, Bureau of Infectious Diseases Diagnosis Control, Korea Disease Control and Prevention Agency, Cheongju 28726, Korea
4
Department of Microbiology, College of Medicine, Hallym University, Chuncheon 24252, Korea
5
Department of Preventive Medicine, Army Headquarters, Gyeryong 32800, Korea
6
Gangwon Institute of Health and Environment, Infectious Disease Research Department, Chuncheon 24203, Korea
7
Force Health Protection and Preventive Medicine, MEDDAC-Korea, 65th Medical Brigade, Unit 15281, APO AP 96271-5281, USA
8
Institute of Medical Research, College of Medicine, Hallym University, Chuncheon 24252, Korea
*
Author to whom correspondence should be addressed.
Submission received: 12 April 2021 / Revised: 1 May 2021 / Accepted: 4 May 2021 / Published: 6 May 2021
(This article belongs to the Special Issue Hantavirus)

Abstract

:
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.

1. Introduction

Orthohantaviruses (family, Hantaviridae; order, Bunyavirales) are enveloped negative-sense 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]. Orthohantaviruses are the etiological agents of hemorrhagic fever with renal syndrome (HFRS), with case fatality rates varying from <1% to 15% in Eurasia [3]. HFRS is mainly caused by Hantaan virus (HTNV), Seoul virus, Dobrava virus, and Puumala virus [4]. The natural reservoir hosts of hantaviruses include rodents (Rodentia), bats (Chiroptera), and insectivores (Soricomorpha) [5,6,7]. Hantaviruses are transmitted to humans when aerosolized infectious particles from the saliva, urine, and feces of infected rodents are inhaled through the respiratory tract [8,9].
Next-generation sequencing (NGS) technologies have been applied to a wide range of fields in virology, including metagenomics of the virome, whole-genome sequencing for population-level surveillance, the tracking of the spread of infectious agents, the development of vaccines, and the identification of putative pathogens [10,11,12,13]. Whole-genome sequencing plays a critical role in understanding the characteristics and transmission dynamics during the outbreaks of emerging or re-emerging viruses that pose a public health threat [14,15,16,17]. However, a significant limitation of whole-genome sequencing for the identification of viral sequences is the ultra-low copy number of target genomes in clinical or environmental sources. To detect low-titer viruses, researchers have developed various NGS approaches, including sequence-independent, single-primer amplification (SISPA); specific oligonucleotide probe-mediated enrichment of target viral nucleic acids; amplicon-based NGS; and small-RNA deep sequencing [18,19,20,21].
MinION (Oxford Nanopore Technologies) is a palm-sized portable sequencer that is smaller and cheaper than conventional sequencing platforms; it allows real-time sequencing in field situations or hospitals [22,23]. The nanopore sequencing has allowed researchers to reveal the dynamic genetic diversity and geographic distribution of the rabies virus in field [24]. Previous studies have used nanopore sequencing in point-of-care testing for the real-time identification of infectious agents, including severe acute respiratory syndrome coronavirus 2, Chikungunya virus (CHIKV), Ebola virus (EBOV), Zika virus (ZIKV), and hepatitis C virus in clinical samples [25,26,27]. However, to our knowledge, nanopore-based NGS has not been used to conduct whole-genome sequencing of orthohantaviruses directly from the lung tissue of a natural host, Apodemus agrarius (the striped field mouse).
In this study, we evaluated whether the MinION device could be used for the whole-genome sequencing of HTNV independently of virus culture and determined whether it could rapidly and accurately generate genomic sequences for the identification of genetic diversity at the variants. The epidemiological surveillance of small mammals showed the genetic diversity and geographical prevalence of HTNV in HFRS-endemic areas in the Republic of Korea (ROK). This study provides important insights into the application of nanopore sequencing for genomic-based diagnosis and analysis of orthohantaviruses from clinical and animal specimens in the field.

2. Materials and Methods

2.1. Ethics Statement

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.

2.2. Sample Collection

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 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.

2.3. Mitochondrial DNA Analysis

Total DNA was extracted from the liver tissues using a High Pure PCR Template Preparation Kit (Roche, Basel, Switzerland). Mitochondrial cytochrome b (Cyt b) gene-specific 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).

2.4. Indirect Immunofluorescence Antibody Test

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 anti-mouse (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.

2.5. Reverse Transcription-Polymerase Chain Reaction

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′, Han-L-F2 (inner): 5′-TGC WGA TGC HAC IAA RTG GTC-3′ and Han-L-R2 (inner): 5′-GCR TCR TCW GAR TGR TGD GCA A-3′ for the L segment; G2F1 (outer): 5′-TGG GTG CAA GTG C-3′, G2-2 (outer): 5′-ACA TGC TGT ACA GCC TGT GCC-3′, G2-1 (inner): 5′-TGG GCT GCA AGT GCA TCA GAG-3′, G2-4 (inner): 5′-ATG GAT TAC AAC CCC AGC TCG-3′ for the M segment; HTN-S6 (outer): 5′-AGC TCI GGA TCC ATI TCA TC-3′ and OSQ84 (outer): 5′-ATC TTA CAT CCT TTG TCG TCC C-3′, HTN-S4 (inner): 5′-GAI IGI TGT CCA CCA ACA TG-3′ and OSQ85 (inner): 5′-AGT TGT CCA CAG CCT CCT TT-3′ for the S segment. 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 product using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) on automated sequencer (ABI 3730XL DNA Analyzer, Applied Biosystems).

2.6. Quantitative 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].

2.7. Multiplex PCR-Based Enrichment

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.

2.8. Library Preparation for Sequencing on MinION Platform

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.

2.9. Library Preparation for Sequencing on MiSeq Platform

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 of HTNV L, M, and S segments were determined by conventional nested RT-PCR. The genomic sequences of HTNV are described in Table S1.

2.10. Phylogenetic 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-likelihood (ML) method in MEGA 7.0 [31]. Topologies were evaluated using bootstrap analysis of 1000 iterations.

3. Results

3.1. Trapping Sites of Small Mammals from the ROK

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).

3.2. Serological and Molecular Prevalence of HTNV Collected from the ROK

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.

3.3. HTNV RNA Loads and Determination of Viral Copy Number from A. agrarius Lung Tissues

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 104 to 105 copies/μL in the evaluated rodents. Aa19-36 and Aa19-236 showed Ct values ranging from 29.4 to 29.7 (102 to 103 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 101 to 102 copies/μL of viral RNA. Aa19-144 and Aa19-152 consisted of 1 to 101 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.

3.4. Multiplex PCR-Based NGS of HTNV Using MinION or MiSeq Platforms

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.

3.5. Phylogeographic Analysis

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 (Munhye-ri) (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.

4. Discussion

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 102 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 MinION 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://0-www-mdpi-com.brum.beds.ac.uk/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..

Author Contributions

Conceptualization: K.P. (Kyungmin Park), S.-H.L., and J.-W.S.; Methodology: K.P. (Kyungmin Park), S.-H.L., J.K., J.L., G.-Y.L., K.P. (Kkothanahreum Park), S.B., and W.-K.K.; Validation: J.K., J.L., G.-Y.L., S.C., S.H.L., J.S.N., and W.-K.K.; Formal analysis: K.P. (Kyungmin Park) and S.-H.L.; Investigation: K.P. (Kyungmin Park)and S.-H.L.; Resources: K.P. (Kyungmin Park), S.-H.L., Y.-J.K., Y.-S.K., H.-C.K., and T.A.K.; Data curation: K.P. (Kyungmin Park) and S.-H.L.; Writing—original draft preparation: K.P. (Kyungmin Park) and S.-H.L.; Writing—review and editing: J.-W.S.; Visualization: K.P. (Kyungmin Park) and S.-H.L.; Supervision: J.-W.S.; Funding acquisition: J.-W.S. All authors have read and approved the published version of the manuscript.

Funding

This study was supported by the Research Program to Solve Social Issues of the National Research Foundation of Korea (NRF) funded by the Ministry of Science and Information and Communication Technology (ICT; NRF-2017M3A9E4061992) and the Agency for Defense Development (UE202026GD). Partial funding was provided by the Research Program to Solve Social Issues of the NRF funded by the Ministry of Science and ICT (NRF-2019R1I1A2A01060902); the Armed Forces Health Surveillance Branch, Global Emerging Infections Surveillance and Response System (AFHSB-GEIS); Silver Spring, MD (ProMIS ID #P0025-2016-2018-ME); and the 65th Medical Brigade, Camp Humphreys, Republic of Korea.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Korea University Institutional Animal Care and Use Committee (KUIACUC #2016-49).

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data generated for this publication have been included in the current manuscript.

Acknowledgments

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vaheri, A.; Henttonen, H.; Voutilainen, L.; Mustonen, J.; Sironen, T.; Vapalahti, O. Hantavirus infections in Europe and their impact on public health. Revi. Med. Virol. 2013, 23, 35–49. [Google Scholar] [CrossRef]
  2. Laenen, L.; Vergote, V.; Calisher, C.H.; Klempa, B.; Klingstrom, J.; Kuhn, J.H.; Maes, P. Hantaviridae: Current Classification and Future Perspectives. Viruses 2019, 11, 788. [Google Scholar] [CrossRef] [Green Version]
  3. Peters, C.J.; Simpson, G.L.; Levy, H. Spectrum of hantavirus infection: Hemorrhagic fever with renal syndrome and hantavirus pulmonary syndrome. Annu. Rev. Med. 1999, 50, 531–545. [Google Scholar] [CrossRef]
  4. Gizzi, M.; Delaere, B.; Weynand, B.; Clement, J.; Maes, P.; Vergote, V.; Laenen, L.; Hjelle, B.; Verroken, A.; Dive, A. Another case of “European hantavirus pulmonary syndrome” with severe lung, prior to kidney, involvement, and diagnosed by viral inclusions in lung macrophages. Eur. J. Clin. Microbiol. Infect. Dis. 2013, 32, 1341–1345. [Google Scholar] [CrossRef] [PubMed]
  5. Yanagihara, R.; Gu, S.H.; Arai, S.; Kang, H.J.; Song, J.W. Hantaviruses: Rediscovery and new beginnings. Virus Res. 2014, 187, 6–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Guo, W.P.; Lin, X.D.; Wang, W.; Tian, J.H.; Cong, M.L.; Zhang, H.L.; Wang, M.R.; Zhou, R.H.; Wang, J.B.; Li, M.H.; et al. Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents. PLoS Pathog. 2013, 9, e1003159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Bennett, S.N.; Gu, S.H.; Kang, H.J.; Arai, S.; Yanagihara, R. Reconstructing the evolutionary origins and phylogeography of hantaviruses. Trends Microbiol. 2014, 22, 473–482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Forbes, K.M.; Sironen, T.; Plyusnin, A. Hantavirus maintenance and transmission in reservoir host populations. Curr. Opin. Virol. 2018, 28, 1–6. [Google Scholar] [CrossRef] [Green Version]
  9. Tsai, T.F. Hemorrhagic fever with renal syndrome: Mode of transmission to humans. Lab. Anim. Sci 1987, 37, 428–430. [Google Scholar]
  10. Radford, A.D.; Chapman, D.; Dixon, L.; Chantrey, J.; Darby, A.C.; Hall, N. Application of next-generation sequencing technologies in virology. J. Gen. Virol. 2012, 93, 1853. [Google Scholar] [CrossRef] [Green Version]
  11. Victoria, J.G.; Kapoor, A.; Dupuis, K.; Schnurr, D.P.; Delwart, E.L. Rapid identification of known and new RNA viruses from animal tissues. PLoS Pathog. 2008, 4, e1000163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Briand, F.-X.; Schmitz, A.; Ogor, K.; Le Prioux, A.; Guillou-Cloarec, C.; Guillemoto, C.; Allée, C.; Le Bras, M.-O.; Hirchaud, E.; Quenault, H. Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: Phylogenetic analyses and markers for zoonotic potential. Eurosurveillance 2017, 22, 30473. [Google Scholar] [CrossRef] [Green Version]
  13. Harvala, H.; Frampton, D.; Grant, P.; Raffle, J.; Ferns, R.B.; Kozlakidis, Z.; Kellam, P.; Pillay, D.; Hayward, A.; Nastouli, E. Emergence of a novel subclade of influenza A (H3N2) virus in London, December 2016 to January 2017. Eurosurveillance 2017, 22, 30466. [Google Scholar] [CrossRef]
  14. Smith, G.J.; Vijaykrishna, D.; Bahl, J.; Lycett, S.J.; Worobey, M.; Pybus, O.G.; Ma, S.K.; Cheung, C.L.; Raghwani, J.; Bhatt, S. Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic. Nature 2009, 459, 1122–1125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gire, S.K.; Goba, A.; Andersen, K.G.; Sealfon, R.S.; Park, D.J.; Kanneh, L.; Jalloh, S.; Momoh, M.; Fullah, M.; Dudas, G. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science 2014, 345, 1369–1372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Neverov, A.; Chumakov, K. Massively parallel sequencing for monitoring genetic consistency and quality control of live viral vaccines. Proc. Natl. Acad. Sci. USA 2010, 107, 20063–20068. [Google Scholar] [CrossRef] [Green Version]
  17. Munnink, B.B.O.; Münger, E.; Nieuwenhuijse, D.F.; Kohl, R.; Van Der Linden, A.; Schapendonk, C.; Van Der Jeugd, H.; Kik, M.; Rijks, J.; Reusken, C. Genomic monitoring to understand the emergence and spread of Usutu virus in the Netherlands, 2016–2018. Sci. Rep. 2020, 10, 1–10. [Google Scholar]
  18. Kim, W.-K.; Kim, J.-A.; Song, D.H.; Lee, D.; Kim, Y.C.; Lee, S.-Y.; Lee, S.-H.; No, J.S.; Kim, J.H.; Kho, J.H. Phylogeographic analysis of hemorrhagic fever with renal syndrome patients using multiplex PCR-based next generation sequencing. Sci. Rep. 2016, 6, 1–8. [Google Scholar] [CrossRef]
  19. Song, D.H.; Kim, W.K.; Gu, S.H.; Lee, D.; Kim, J.A.; No, J.S.; Lee, S.H.; Wiley, M.R.; Palacios, G.; Song, J.W.; et al. Sequence-Independent, Single-Primer Amplification Next-Generation Sequencing of Hantaan Virus Cell Culture-Based Isolates. Am. J. Trop. Med. Hyg. 2017, 96, 389–394. [Google Scholar] [CrossRef] [PubMed]
  20. No, J.S.; Kim, W.-K.; Cho, S.; Lee, S.-H.; Kim, J.-A.; Lee, D.; Song, D.H.; Gu, S.H.; Jeong, S.T.; Wiley, M.R. Comparison of targeted next-generation sequencing for whole-genome sequencing of Hantaan orthohantavirus in Apodemus agrarius lung tissues. Sci. Rep. 2019, 9, 1–9. [Google Scholar] [CrossRef]
  21. Wang, F.; Sun, Y.; Ruan, J.; Chen, R.; Chen, X.; Chen, C.; Kreuze, J.F.; Fei, Z.; Zhu, X.; Gao, S. Using small RNA deep sequencing data to detect human viruses. BioMed Res. Int. 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
  22. McNaughton, A.L.; Roberts, H.E.; Bonsall, D.; de Cesare, M.; Mokaya, J.; Lumley, S.F.; Golubchik, T.; Piazza, P.; Martin, J.B.; de Lara, C. Illumina and Nanopore methods for whole genome sequencing of hepatitis B virus (HBV). Sci. Rep. 2019, 9, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Faria, N.R.; da Silva Azevedo, R.d.S.; Kraemer, M.U.; Souza, R.; Cunha, M.S.; Hill, S.C.; Thézé, J.; Bonsall, M.B.; Bowden, T.A.; Rissanen, I. Zika virus in the Americas: Early epidemiological and genetic findings. Science 2016, 352, 345–349. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gigante, C.M.; Yale, G.; Condori, R.E.; Costa, N.C.; Long, N.V.; Minh, P.Q.; Chuong, V.D.; Tho, N.D.; Thanh, N.T.; Thin, N.X. Portable Rabies Virus Sequencing in Canine Rabies Endemic Countries Using the Oxford Nanopore MinION. Viruses 2020, 12, 1255. [Google Scholar] [CrossRef] [PubMed]
  25. Paden, C.R.; Tao, Y.; Queen, K.; Zhang, J.; Li, Y.; Uehara, A.; Tong, S. Rapid, sensitive, full-genome sequencing of severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 2020, 26, 2401. [Google Scholar] [CrossRef]
  26. Quick, J.; Loman, N.J.; Duraffour, S.; Simpson, J.T.; Severi, E.; Cowley, L.; Bore, J.A.; Koundouno, R.; Dudas, G.; Mikhail, A.; et al. Real-time, portable genome sequencing for Ebola surveillance. Nature 2016, 530, 228–232. [Google Scholar] [CrossRef] [Green Version]
  27. Greninger, A.L.; Naccache, S.N.; Federman, S.; Yu, G.; Mbala, P.; Bres, V.; Stryke, D.; Bouquet, J.; Somasekar, S.; Linnen, J.M. Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis. Genome Med. 2015, 7, 1–13. [Google Scholar] [CrossRef] [Green Version]
  28. Irwin, D.M.; Kocher, T.D.; Wilson, A.C. Evolution of the cytochrome b gene of mammals. J. Mol. Evol. 1991, 32, 128–144. [Google Scholar] [CrossRef]
  29. Lee, G.-Y.; Kim, W.-K.; Park, K.; Lee, S.-H.; Hwang, J.; No, J.S.; Cho, S.; Lee, D.; Song, D.-H.; Gu, S.H.; et al. Phylogeographic diversity and hybrid zone of Hantaan orthohantavirus collected in Gangwon Province, Republic of Korea. PLOS Negl. Trop. Dis. 2020, 14, e0008714. [Google Scholar] [CrossRef]
  30. No, J.S.; Kim, W.K.; Kim, J.A.; Lee, S.H.; Lee, S.Y.; Kim, J.H.; Kho, J.H.; Lee, D.; Song, D.H.; Gu, S.H.; et al. Detection of Hantaan virus RNA from anti-Hantaan virus IgG seronegative rodents in an area of high endemicity in Republic of Korea. Microbiol. Immunol. 2016, 60, 268–271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Kabwe, E.; Davidyuk, Y.; Shamsutdinov, A.; Garanina, E.; Martynova, E.; Kitaeva, K.; Malisheni, M.; Isaeva, G.; Savitskaya, T.; Urbanowicz, R.A.; et al. Orthohantaviruses, Emerging Zoonotic Pathogens. Pathogens 2020, 9, 775. [Google Scholar] [CrossRef]
  33. Korea Disease Control and Prevention Agency (KDCA). Infectious Disease Portal. Available online: http://www.cdc.go.kr/npt/biz/npp/ist/bass/bassDissStatsMain.do (accessed on 10 March 2021).
  34. Kim, J.A.; Kim, W.K.; No, J.S.; Lee, S.H.; Lee, S.Y.; Kim, J.H.; Kho, J.H.; Lee, D.; Song, D.H.; Gu, S.H.; et al. Genetic Diversity and Reassortment of Hantaan Virus Tripartite RNA Genomes in Nature, the Republic of Korea. PLoS Negl. Trop. Dis. 2016, 10, e0004650. [Google Scholar] [CrossRef] [PubMed]
  35. Ryou, J.; Lee, H.I.; Yoo, Y.J.; Noh, Y.T.; Yun, S.M.; Kim, S.Y.; Shin, E.H.; Han, M.G.; Ju, Y.R. Prevalence of hantavirus infection in wild rodents from five provinces in Korea, 2007. J. Wild Dis. 2011, 47, 427–432. [Google Scholar] [CrossRef] [PubMed]
  36. Joshi, Y.P.; Kim, E.H.; Cheong, H.K. 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. BMC Infect. Dis. 2017, 17, 406. [Google Scholar] [CrossRef] [Green Version]
  37. Song, J.W.; Moon, S.S.; Gu, S.H.; Song, K.J.; Baek, L.J.; Kim, H.C.; Kijek, T.; O’Guinn, M.L.; Lee, J.S.; Turell, M.J.; et al. Hemorrhagic fever with renal syndrome in 4 US soldiers, South Korea, 2005. Emerg Infect. Dis. 2009, 15, 1833–1836. [Google Scholar] [CrossRef]
  38. Kim, W.K.; No, J.S.; Lee, D.; Jung, J.; Park, H.; Yi, Y.; Kim, J.A.; Lee, S.H.; Kim, Y.; Park, S.; et al. Active Targeted Surveillance to Identify Sites of Emergence of Hantavirus. Clin. Infect. Dis. 2020, 70, 464–473. [Google Scholar] [CrossRef]
  39. Kim, W.-K.; Cho, S.; Lee, S.-H.; No, J.S.; Lee, G.-Y.; Park, K.; Lee, D.; Jeong, S.T.; Song, J.-W. Genomic Epidemiology and Active Surveillance to Investigate Outbreaks of Hantaviruses. Front. Cell. Infect. Microbiol. 2021, 10, 803. [Google Scholar] [CrossRef]
  40. Domingo, E.; Martin, V.; Perales, C.; Grande-Perez, A.; Garcia-Arriaza, J.; Arias, A. Viruses as quasispecies: Biological implications. Quasispecies 2006, 299, 51–82. [Google Scholar]
  41. Vignuzzi, M.; Stone, J.K.; Arnold, J.J.; Cameron, C.E.; Andino, R. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 2006, 439, 344–348. [Google Scholar] [CrossRef]
  42. Gardy, J.L.; Loman, N.J. Towards a genomics-informed, real-time, global pathogen surveillance system. Nat. Rev. Genet 2018, 19, 9–20. [Google Scholar] [CrossRef]
  43. Oh, M.-d.; Park, W.B.; Park, S.-W.; Choe, P.G.; Bang, J.H.; Song, K.-H.; Kim, E.S.; Kim, H.B.; Kim, N.J. Middle East respiratory syndrome: What we learned from the 2015 outbreak in the Republic of Korea. Korean J. Intern. Med. 2018, 33, 233. [Google Scholar] [CrossRef]
  44. Ladner, J.T.; Beitzel, B.; Chain, P.S.; Davenport, M.G.; Donaldson, E.; Frieman, M.; Kugelman, J.; Kuhn, J.H.; O’Rear, J.; Sabeti, P.C. Standards for sequencing viral genomes in the era of high-throughput sequencing. Am. Soc. Microbiol. 2014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Macalalad, A.R.; Zody, M.C.; Charlebois, P.; Lennon, N.J.; Newman, R.M.; Malboeuf, C.M.; Ryan, E.M.; Boutwell, C.L.; Power, K.A.; Brackney, D.E. Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data. PLoS Comput. Biol. 2012, 8, e1002417. [Google Scholar] [CrossRef]
  46. Wang, C.; Mitsuya, Y.; Gharizadeh, B.; Ronaghi, M.; Shafer, R.W. Characterization of mutation spectra with ultra-deep pyrosequencing: Application to HIV-1 drug resistance. Genome Res. 2007, 17, 1195–1201. [Google Scholar] [CrossRef] [Green Version]
  47. Munnink, B.B.O.; Nieuwenhuijse, D.F.; Stein, M.; O’Toole, Á.; Haverkate, M.; Mollers, M.; Kamga, S.K.; Schapendonk, C.; Pronk, M.; Lexmond, P. Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands. Nat. Med. 2020, 26, 1405–1410. [Google Scholar] [CrossRef]
  48. Illumina Technical Note: Sequencing. Quality scores for Next-Generation Sequencing. Assessing sequencing accuracy using Phred quality scoring. Available online: https://www.illumina.com/documents/products/technotes/technote_Q-Scores.pdf (accessed on 8 March 2021).
  49. Jain, M.; Tyson, J.R.; Loose, M.; Ip, C.L.C.; Eccles, D.A.; O’Grady, J.; Malla, S.; Leggett, R.M.; Wallerman, O.; Jansen, H.J.; et al. MinION Analysis and Reference Consortium: Phase 2 data release and analysis of R9.0 chemistry. F1000Research 2017, 6, 760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Tyson, J.R.; O’Neil, N.J.; Jain, M.; Olsen, H.E.; Hieter, P.; Snutch, T.P. MinION-based long-read sequencing and assembly extends the Caenorhabditis elegans reference genome. Genome Res. 2018, 28, 266–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Stancu, M.C.; Van Roosmalen, M.J.; Renkens, I.; Nieboer, M.M.; Middelkamp, S.; De Ligt, J.; Pregno, G.; Giachino, D.; Mandrile, G.; Valle-Inclan, J.E. Mapping and phasing of structural variation in patient genomes using nanopore sequencing. Nat. Commun. 2017, 8, 1–13. [Google Scholar]
  52. Wick, R.R.; Judd, L.M.; Holt, K.E. Performance of neural network basecalling tools for Oxford Nanopore sequencing. Genome Biol. 2019, 20, 129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Jain, M.; Koren, S.; Miga, K.H.; Quick, J.; Rand, A.C.; Sasani, T.A.; Tyson, J.R.; Beggs, A.D.; Dilthey, A.T.; Fiddes, I.T.; et al. Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat. Biotechnol. 2018, 36, 338–345. [Google Scholar] [CrossRef] [Green Version]
  54. Loman, N.J.; Quick, J.; Simpson, J.T. A complete bacterial genome assembled de novo using only nanopore sequencing data. Nat. Methods 2015, 12, 733–735. [Google Scholar] [CrossRef] [PubMed]
  55. Koren, S.; Walenz, B.P.; Berlin, K.; Miller, J.R.; Bergman, N.H.; Phillippy, A.M. Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017, 27, 722–736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Muyangwa, M.; Martynova, E.V.; Khaiboullina, S.F.; Morzunov, S.P.; Rizvanov, A.A. Hantaviral Proteins: Structure, Functions, and Role in Hantavirus Infection. Front Microbiol. 2015, 6, 1326. [Google Scholar] [CrossRef] [PubMed]
  57. Spiropoulou, C.F.; Goldsmith, C.S.; Shoemaker, T.R.; Peters, C.J.; Compans, R.W. Sin Nombre virus glycoprotein trafficking. Virology 2003, 308, 48–63. [Google Scholar] [CrossRef] [Green Version]
  58. Grubaugh, N.D.; Ladner, J.T.; Lemey, P.; Pybus, O.G.; Rambaut, A.; Holmes, E.C.; Andersen, K.G. Tracking virus outbreaks in the twenty-first century. Nat Microbiol. 2019, 4, 10–19. [Google Scholar] [CrossRef] [PubMed]
  59. Xu, Y.; Lewandowski, K.; Lumley, S.; Pullan, S.; Vipond, R.; Carroll, M.; Foster, D.; Matthews, P.C.; Peto, T.; Crook, D. Detection of viral pathogens with multiplex nanopore MinION sequencing: Be careful with cross-Talk. Front. Microbiol. 2018, 9, 2225. [Google Scholar] [CrossRef]
  60. Quick, J.; Grubaugh, N.D.; Pullan, S.T.; Claro, I.M.; Smith, A.D.; Gangavarapu, K.; Oliveira, G.; Robles-Sikisaka, R.; Rogers, T.F.; Beutler, N.A.; et al. Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples. Nat. Protoc. 2017, 12, 1261–1276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
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.
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.
Viruses 13 00847 g001
Figure 2. Geographical location of the trapping sites of Hantaan virus (HTNV) collected from Gyeonggi and Gangwon Provinces, the Republic of Korea. The geographic map shows different trapping areas where small mammals were captured at Gyeonggi and Gangwon Provinces in 2019. 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.
Figure 2. Geographical location of the trapping sites of Hantaan virus (HTNV) collected from Gyeonggi and Gangwon Provinces, the Republic of Korea. The geographic map shows different trapping areas where small mammals were captured at Gyeonggi and Gangwon Provinces in 2019. 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.
Viruses 13 00847 g002
Figure 3. Average depth of coverage for Hantaan virus (HTNV) obtained from MinION sequencing at different cumulated running times. The illustration shows the average depth of coverage generated from amplicon-based MinION sequencing based on (A) HTNV L segment, (B) M segment, and (C) S segment genomes by using different cumulated running times (2, 4, 8, and 16 h).
Figure 3. Average depth of coverage for Hantaan virus (HTNV) obtained from MinION sequencing at different cumulated running times. The illustration shows the average depth of coverage generated from amplicon-based MinION sequencing based on (A) HTNV L segment, (B) M segment, and (C) S segment genomes by using different cumulated running times (2, 4, 8, and 16 h).
Viruses 13 00847 g003
Figure 4. Phylogenetic analysis of Hantaan virus (HTNV) whole-genome sequences by using amplicon next-generation sequencing (NGS) and MiSeq and MinION platforms. Phylogenetic trees were generated using maximum-likelihood method and bootstrap of 1000 iterations based on (A) HTNV L segment (1–6533 nt), (B) HTNV M segment (1–3616 nt), and (C) HTNV S segment (1–1696 nt). The numbers at each node are bootstrap probabilities, as determined for 1000 iterations. The newly obtained HTNV are shown in bold lettering and designated by color indicating the specific sites: blue, Jangam-ri of Pocheon; green, Dosin-ri of Yeoncheon; brown, Sang-ri of Yeoncheon; orange, Gwanu-ri of Cheorwon; red, Dochang-ri of Cheorwon; and violet, Pungsan-ri of Hwacheon. The newly acquired HTNV genomic sequences in this study are described in Table S1. The following other HTNV genome sequences were used in the analysis: HTNV Aa03-387 (L segment, KT934958; M segment, KT934992; S segment, KT935026), Aa05-331 (L segment, KT934962; M segment, KT934996; S segment, KT935030), Aa09-410 (L segment, KU207177; M segment, KU207185; S segment, KU207193), Aa09-948 (L segment, KT934966; M segment, KT935000; S segment, KT935034), Aa10-434 (L segment, KT934970; M segment, KT935004; S segment, KT935038), Aa10-518 (L segment, KT934971; M segment, KT935005; S segment, KT935039), Aa14-266 (L segment, KT934979; M segment, KT935013; S segment, KT935047), Aa14-272 (L segment, KT934980; M segment, KT935014; S segment, KT935048), Aa14-362 (L segment, KT934981; M segment, KT935015; S segment, KT935049), Aa14-368 (L segment, KT934982; M segment, KT935016; S segment, KT935050), Aa15-56 (L segment, KU207179; M segment, KU207187; S segment, KU207195), Aa15-58 (L segment, KU207180; M segment, KU207188; S segment, KU207196), Aa15-82 (L segment, MT012572; M segment, MT012560; S segment, MT012548), Aa15-84 (L segment, MT012573; M segment, MT012561; S segment, MT012549), Aa17-353 (L segment, MT012575; M segment, MT012563; S segment, MT012551), Aa17-367 (L segment, MT012576; M segment, MT012564; S segment, MT012552), Aa17-421 (L segment, MT012577; M segment, MT012565; S segment, MT012553), Aa17-422 (L segment, MT012578; M segment, MT012566; S segment, MT012554), Aa18-164 (L segment, MT012579; M segment, MT012567; S segment, MT012555), Aa18-179 (L segment, MT012580; M segment, MT012568; S segment, MT012556), HTNV 76-118 (L segment, NC_005222; M segment, M14627; S segment, M14626), HV004 (L segment, JQ083393; M segment, JQ083394; S segment, JQ093395).
Figure 4. Phylogenetic analysis of Hantaan virus (HTNV) whole-genome sequences by using amplicon next-generation sequencing (NGS) and MiSeq and MinION platforms. Phylogenetic trees were generated using maximum-likelihood method and bootstrap of 1000 iterations based on (A) HTNV L segment (1–6533 nt), (B) HTNV M segment (1–3616 nt), and (C) HTNV S segment (1–1696 nt). The numbers at each node are bootstrap probabilities, as determined for 1000 iterations. The newly obtained HTNV are shown in bold lettering and designated by color indicating the specific sites: blue, Jangam-ri of Pocheon; green, Dosin-ri of Yeoncheon; brown, Sang-ri of Yeoncheon; orange, Gwanu-ri of Cheorwon; red, Dochang-ri of Cheorwon; and violet, Pungsan-ri of Hwacheon. The newly acquired HTNV genomic sequences in this study are described in Table S1. The following other HTNV genome sequences were used in the analysis: HTNV Aa03-387 (L segment, KT934958; M segment, KT934992; S segment, KT935026), Aa05-331 (L segment, KT934962; M segment, KT934996; S segment, KT935030), Aa09-410 (L segment, KU207177; M segment, KU207185; S segment, KU207193), Aa09-948 (L segment, KT934966; M segment, KT935000; S segment, KT935034), Aa10-434 (L segment, KT934970; M segment, KT935004; S segment, KT935038), Aa10-518 (L segment, KT934971; M segment, KT935005; S segment, KT935039), Aa14-266 (L segment, KT934979; M segment, KT935013; S segment, KT935047), Aa14-272 (L segment, KT934980; M segment, KT935014; S segment, KT935048), Aa14-362 (L segment, KT934981; M segment, KT935015; S segment, KT935049), Aa14-368 (L segment, KT934982; M segment, KT935016; S segment, KT935050), Aa15-56 (L segment, KU207179; M segment, KU207187; S segment, KU207195), Aa15-58 (L segment, KU207180; M segment, KU207188; S segment, KU207196), Aa15-82 (L segment, MT012572; M segment, MT012560; S segment, MT012548), Aa15-84 (L segment, MT012573; M segment, MT012561; S segment, MT012549), Aa17-353 (L segment, MT012575; M segment, MT012563; S segment, MT012551), Aa17-367 (L segment, MT012576; M segment, MT012564; S segment, MT012552), Aa17-421 (L segment, MT012577; M segment, MT012565; S segment, MT012553), Aa17-422 (L segment, MT012578; M segment, MT012566; S segment, MT012554), Aa18-164 (L segment, MT012579; M segment, MT012567; S segment, MT012555), Aa18-179 (L segment, MT012580; M segment, MT012568; S segment, MT012556), HTNV 76-118 (L segment, NC_005222; M segment, M14627; S segment, M14626), HV004 (L segment, JQ083393; M segment, JQ083394; S segment, JQ093395).
Viruses 13 00847 g004
Table 1. Trapping results of small mammals collected from Gyeonggi and Gangwon Provinces in 2019.
Table 1. Trapping results of small mammals collected from Gyeonggi and Gangwon Provinces in 2019.
SpeciesSmall Mammal Trapping LocationTotal (%)
PajuPocheonPyeongtaekYeoncheonCheorwonChuncheonHwacheon
Apodemus agrarius71841481053846303 (86.6)
Crocidura lasiura-a--174113 (3.7)
Crocidura shantungensis----1113 (0.9)
Mus musculus--7-9-117 (4.9)
Myodes regulus2--11--4 (1.1)
Micromys minutus--2-12-5 (1.4)
Microtus fortis----1--1 (0.3)
Rattus norvegicus--3----3 (0.9)
Tscherskia triton------11 (0.3)
Total91853501254550350 (100)
a: no collection.
Table 2. Serological and molecular prevalence of Hantaan virus (HTNV) in Apodemus agrarius captured from Gyeonggi and Gangwon Provinces in 2019.
Table 2. Serological and molecular prevalence of Hantaan virus (HTNV) in Apodemus agrarius captured from Gyeonggi and Gangwon Provinces in 2019.
Trapping SiteNumber of Captured
A. agrarius
Seropositivity for
Anti-HTNV IgG (%)
HTNV RNA
Positivity (%)
MaleFemaleTotalMaleFemaleTotal
Paju70/20/50/7-a--
Pocheon181/9 (11.1)0/91/18 (5.6)1/1 (100)-1/1 (100)
Pyeongtaek410/230/180/41---
Yeoncheon482/22 (9.1)3/26 (11.5)5/48 (10.4)1/2 (50)1/3 (33.3)2/5 (40)
Cheorwon1055/60 (8.3)2/45 (4.4)7/105 (6.7)4/5 (80)2/2 (100)6/7 (85.7)
Chuncheon380/181/20 (5)1/38 (2.6)-0/10/1
Hwacheon464/20 (20)2/26 (7.7)6/46 (13.0)3/4 (75)0/23/6 (50)
Total30312/154 (7.8)8/149 (5.4)20/303 (6.6)9/12 (75)3/8 (37.5)12/20 (60)
a: seronegative sample was not analyzed using RT-PCR.
Table 3. Results of MinION sequencing based on viral RNA copy number of Hantaan virus (HTNV) in the lung tissues of rodents collected from Gyeonggi and Gangwon Provinces, ROK, in 2019.
Table 3. Results of MinION sequencing based on viral RNA copy number of Hantaan virus (HTNV) in the lung tissues of rodents collected from Gyeonggi and Gangwon Provinces, ROK, in 2019.
Viral RNA
Copy Number
(Copies/µL)
SampleSiteAnti-HTNV
IgG Titer
Nested
RT-PCR
Ct
Value
HTNV Genomes, % Coverage
(Minimum Read Depth, >10×)
L SegmentM SegmentS Segment
104 to 105Aa19-233Cheorwon1:64 aPos20.596.898.5100
Aa19-278Cheorwon1:32 aPos22.292.998.4100
102 to 103Aa19-236Cheorwon1:32 aPos29.495.398.4100
Aa19-36Cheorwon1:512 bPos29.791.798.4100
10 to 102Aa19-89Yeoncheon1:64 aPos31.895.998.4100
Aa19-57Pocheon1:8096 bPos32.795.898.4100
Aa19-167Yeoncheon1:256 aPos34.095.996.4100
Aa19-38Cheorwon1:16,384 bPos34.492.398.4100
Aa19-153Hwacheon1:4096 bPos34.592.798.4100
1 to 10Aa19-152Hwacheon1:512 aPos36.482.697.092.1
Aa19-144Hwacheon1:256 bPos37.388.896.495.1
0 to 1Aa19-39Cheorwon1:512 bPos40.065.490.692.2
HTNV: Hantaan virus; IgG, immunoglobulin G; nested RT-PCR, nested reverse transcription-polymerase chain reaction; Ct, cycle threshold; Aa, Apodemus agrarius; Pos, Positive; Neg, Negative; a, IFA test was performed from heart fluids; b, IFA test was performed from sera.
Table 4. Sequencing error and consensus accuracy of amplicon-based MinION sequencing for Hantaan virus.
Table 4. Sequencing error and consensus accuracy of amplicon-based MinION sequencing for Hantaan virus.
Viral RNA
Copy Number
(Copies/µL)
SampleSiteCt
Value
L SegmentM SegmentS Segment
InsertionDeletionMismatchAccuracy aInsertionDeletionMismatchAccuracyInsertionDeletionMismatchAccuracy
104 to 105Aa19-233Cheorwon20.5031299.805299.901499.8
Aa19-278Cheorwon22.204999.906299.9010100
102 to 103Aa19-236Cheorwon29.405999.906399.901199.9
Aa19-36Cheorwon29.7031299.805299.901199.9
10 to 102Aa19-89Yeoncheon31.803699.905699.801199.9
Aa19-57Pocheon32.704999.905499.901299.9
Aa19-167Yeoncheon34.004899.904499.9010100
Aa19-38Cheorwon34.4031299.805799.801299.9
Aa19-153Hwacheon34.5051499.804399.901399.8
1 to 10Aa19-152Hwacheon36.4051899.7021199.701799.6
Aa19-144Hwacheon37.3041499.8041199.701799.6
0 to 1Aa19-39Cheorwon40.0025999.1023699.0012898.3
TotalAverage (%)03.815.299.804.47.699.8014.799.7
Ct, cycle threshold; Aa, Apodemus agrarius; a, Percentage of nucleotide similarity between Oxford nanopore and Illumina MiSeq platforms.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Park, K.; Lee, S.-H.; Kim, J.; Lee, J.; Lee, G.-Y.; Cho, S.; Lee, S.H.; Park, K.; No, J.S.; Budhathoki, S.; et al. Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea. Viruses 2021, 13, 847. https://0-doi-org.brum.beds.ac.uk/10.3390/v13050847

AMA Style

Park K, Lee S-H, Kim J, Lee J, Lee G-Y, Cho S, Lee SH, Park K, No JS, Budhathoki S, et al. Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea. Viruses. 2021; 13(5):847. https://0-doi-org.brum.beds.ac.uk/10.3390/v13050847

Chicago/Turabian Style

Park, Kyungmin, Seung-Ho Lee, Jongwoo Kim, Jingyeong Lee, Geum-Young Lee, Seungchan Cho, Seung Ho Lee, Kkothanahreum Park, Jin Sun No, Shailesh Budhathoki, and et al. 2021. "Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea" Viruses 13, no. 5: 847. https://0-doi-org.brum.beds.ac.uk/10.3390/v13050847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop