Skip to main content

METHODS article

Front. Microbiol., 13 November 2018
Sec. Infectious Agents and Disease

Whole Genome Sequencing of Influenza A and B Viruses With the MinION Sequencer in the Clinical Setting: A Pilot Study

\r\nKazuo Imai,,Kazuo Imai1,2,3Kaku TamuraKaku Tamura3Tomomi TanigakiTomomi Tanigaki4Mari TakizawaMari Takizawa5Eiko NakayamaEiko Nakayama5Takahiko TaniguchiTakahiko Taniguchi5Misako OkamotoMisako Okamoto5Yasumasa NishiyamaYasumasa Nishiyama4Norihito Tarumoto,Norihito Tarumoto1,2Kotaro Mitsutake,Kotaro Mitsutake2,6Takashi Murakami,Takashi Murakami2,7Shigefumi Maesaki,Shigefumi Maesaki1,2Takuya Maeda,*Takuya Maeda2,7*
  • 1Department of Infectious Disease and Infection Control, Saitama Medical University, Saitama, Japan
  • 2Center for Clinical Infectious Diseases and Research, Saitama Medical University, Saitama, Japan
  • 3Department of Infectious Diseases, Self-Defense Forces Central Hospital, Japan Ground Self-Defense Forces, Tokyo, Japan
  • 4NBC Counter Medical Unit, Japan Ground Self-Defense Forces, Tokyo, Japan
  • 5Camp Asaka Medical Office, Japan Ground Self-Defense Forces, Tokyo, Japan
  • 6Department of Infectious Disease and Infection Control, Saitama Medical University International Medical Center, Saitama, Japan
  • 7Department of Microbiology, Saitama Medical University, Saitama, Japan

Introduction: Whole genome sequencing (WGS) of influenza viruses is important for preparing vaccines and coping with newly emerging viruses. However, WGS is difficult to perform using conventional next-generation sequencers in developing countries, where facilities are often inadequate. In this study, we developed a high-throughput WGS method for influenza viruses in clinical specimens with the MinION portable sequencer.

Methods: Whole genomes of influenza A and B viruses were amplified by multiplex RT-PCR from 13 clinical specimens collected in Tokyo, Japan. Barcode tags for multiplex MinION sequencing were added with each multiplex RT-PCR amplicon by nested PCR with custom barcoded primers. All barcoded amplicons were mixed and multiplex sequencing using the MinION sequencer with 1D2 sequencing kit. In addition, multiplex RT-PCR amplicons generated from each clinical specimen were sequenced using the Illumina MiSeq platform to validate the performance of MinION sequencer. The accuracy, recall, and precision rates of MinION sequencing were calculated by comparing the results of variant calling in the Illumina MiSeq platform and MinION sequencer.

Results: Whole genomes of influenza A and B viruses were successfully amplified by multiplex RT-PCR from 13 clinical samples. We identified 6 samples as influenza type A virus H3N2 subtype and 7 as influenza B virus Yamagata lineage using the Illumina MiSeq platform. The overall accuracy, recall, and precision rates of the MinION sequencer were, respectively 99.95%, 89.41%, and 97.88% from 1D reads and 99.97%, 93.28%, and 99.86% from 1D2 reads.

Conclusion: We developed a novel WGS method for influenza A and B viruses. It is necessary to improve read accuracy and analytical tools in order to better utilize the MinION sequencer for real-time monitoring of genetic rearrangements and for evaluation of newly emerging viruses.

Introduction

Influenza viruses, which are enveloped, single-stranded negative-sense RNA viruses, belong to the Orthomyxoviridae family and cause respiratory diseases in humans around the world. While there are four types of seasonal influenza viruses – types A, B, C, and Thogotovirus – influenza A and B viruses in particular can cause seasonal epidemics that have high rates of morbidity and mortality (Taubenberger and Morens, 2008). The genomes of influenza A and B viruses are composed of the following eight RNA segments: segment 1, RNA polymerase PB2 unit; segment 2, RNA polymerase PB1 unit; segment 3, RNA polymerase PA unit; segment 4, hemagglutinin (HA); segment 5, nucleoprotein (NP); segment 6, neuraminidase (NA); segment 7, matrix (M); and segment 8, non-structural protein (NS).

The genomes of influenza viruses frequently develop genetic changes and have a high level of genetic diversity, which contributes to their ability to evade the human immune system (Grenfell et al., 2004). To cope with the emergence of new influenza virus strains and select appropriate vaccine strains, considerable effort is being devoted to monitoring the genetic rearrangements of these viruses (Ghedin et al., 2005; Shu and McCauley, 2017). A simple amplification method using RT-PCR for influenza virus has been developed and is contributing to the efficient surveillance of influenza viruses via next-generation sequencing (NGS) (Zhou et al., 2009; Zhou et al., 2014; Zhao et al., 2016). NGS technology is used in a wide range of applications, including population surveillance for emerging strains and the subsequent creation of new vaccines, outbreak research in communities, and tracking of nosocomial infections (Briand et al., 2017; Harvala et al., 2017; Houlihan et al., 2018; MacFadden et al., 2018). However, amplicon sequencing using NGS platforms such as 454 (Roche Life Science, Penzberg, Germany), the Illumina platform (Illumina, San Diego, CA, United States), and Ion Torrent technology (Thermo Fisher Scientific, Waltham, MA, United States) has high initial and running costs and is time-consuming. Therefore, it is challenging to establish facilities with an NGS platform for the surveillance of influenza viruses in developing countries (Vemula et al., 2016).

The MinION nanopore sequencer, recently developed by Oxford Nanopore Technologies (ONT; Oxford, United Kingdom), is a palm-size portable sequencer that can read long-read single DNA molecules in real time and is regarded as third-generation DNA sequencing technology. The most notable feature of the MinION sequencer is its portability, which makes it different from other NGS platforms. For this reason, the MinION sequencer has been used as an alternative means of whole genome sequencing (WGS) and can be used for genomic surveillance and infection control for certain pathogens, taking advantage of its high portability, and its quick and easy sample preparation. However, its high sequence error rate remains a challenge (Quick et al., 2015, 2016; Hoenen et al., 2016). In 2017, ONT released the new 1D2 method, which has higher single-read accuracy than the standard 1D method1. The WGS method for influenza A virus using the MinION sequencer has been previously reported (Wang et al., 2015). However, the method still has some practical problems, such as a high per-sample cost because of the single plex sequencing, applicability to only influenza A virus, and unknown performance in clinical settings.

In this study, we developed, a quick and accurate multiplex WGS method for influenza A and B viruses using clinical specimens, the MinION nanopore sequencer, and the 1D2 method.

Materials and Methods

Clinical Samples and Ethics

Clinical samples were collected from suspected influenza patients at the Japan Self-Defense Forces Central Hospital and Medical Office of Camp Asaka, Tokyo, Japan, from December 2017 to April 2018. All patients were diagnosed based on commercially available rapid diagnostic tests for influenza using nasopharyngeal swab specimens, according to the manufacturer’s instructions (QuickNavi-Flu 2; Denka Seiken, Niigata, Japan). In particular, the samples obtained were composed of two specimens, which were separated for subsequent tests. One swab was used for the rapid diagnosis of influenza and the other for the subsequent RNA extraction of influenza virus. Each swab specimen was stored at -80°C until just before use. A diagnosis of influenza was ultimately confirmed via RT-PCR according to previously reported protocols for influenza viruses A (M30F2/08 and M264R3/08 for all subtypes) and B (Bvf224 and Bvr507 for Victoria lineage, and Bvf226 and BYr613 for Yamagata lineage). Both primer sets are recommended for the detection of influenza viruses by the World Health Organization [World Health Organization (WHO), 2017]. To evaluate our procedures in clinical practice, 13 nasopharyngeal swab samples were confirmed as influenza-positive (6 = influenza A virus and 7 = influenza B virus Yamagata lineage) and used in our study.

This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Ethics Committee of the Self-Defense Forces Central Hospital (No. 29-012). All samples were collected after written informed consent was obtained.

RNA Extraction and DNase Treatment

Total RNA was extracted using a High Pure Viral RNA Kit (Roche Life Science) according to the manufacturer’s instructions. Additionally, each eluted RNA was treated with 4U of TURBO DNase (Thermo Fisher Scientific) for 30 min at 37°C to remove the human genomic DNA and purified again by Agencourt RNAClean XP (Beckman Coulter, Brea, CA, United States) according to the manufacturer’s instructions.

Multiplex RT- PCR and Purification

Supernatants of each purified RNA sample (5 μL) were then used for multiplex RT-PCR reactions. In this scheme, viral RNA segments of influenza A (PB2 = 2.3 kb, PB1 = 2.3 kb, PA = 2.2 kb, HA = 1.8 kb, NP = 1.6 kb, NA = 1.4 kb, M = 1.0 kb, and NS = 0.9 kb) were simultaneously amplified using two primers (Uni12, Uni13; Supplementary Table S1) targeting the highly conserved sequence of viral RNA termini, which are commonly present at the ends of the PB2, PB1, PA, HA, NP, NA, M, and NS genome segments, following protocols previously reported for influenza A virus (Zhou et al., 2009). For influenza B (PB2 = 2.3 kb, PB1 = 2.3 kb, PA = 2.3 kb, HA = 1.8 kb, NP = 1.8 kb, NA = 1.5 kb, M = 1.1 kb, and NS = 1.0 kb), multiplex amplification was simultaneously performed with multiple primer sets targeting the corresponding conserved sequence of each viral RNA termini in accordance with a previous report (Supplementary Table S1; Zhou et al., 2014). Multiplex RT-PCR amplicons were generated by the SuperScriptTM III One-Step RT-PCR System with PlatinumTM Taq High Fidelity DNA Polymerase (Thermo Fisher Scientific). PCR products were analyzed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States) with an Agilent DNA 7500 Kit (Agilent Technologies). Each PCR product was purified again using Agencourt AMPure XP according to the manufacturer’s instructions.

Sequencing of Multiplex RT-PCR Products Using the Illumina MiSeq Sequencer

The complete influenza virus genome was analyzed with multiplex RT-PCR amplicon sequencing using the Illumina MiSeq platform. Paired-end libraries for MiSeq platform were prepared using Nextera XT DNA Library Prep Kit (Illumina) and Nextera XT Index Kit (Illumina), and sequencing was performed using a 300-cycle (2 × 150-bp paired-end) MiSeq v2 reagent kit (Illumina) via MiSeq platform according to the manufacturer’s protocols.

Analysis of Illumina MiSeq Platform Data

Low-quality reads in demultiplexed data were removed and 3′ terminal nucleotides were trimmed with SeqKit and mapped with BWA-MEM (v 0.7.15) using the reference sequences for influenza A (GenBank; CY259943-50) and influenza B viruses (GenBank; MH233733-40). Single nucleotide variant calling from mapped data was obtained by Samtools, Picard, and GATK via “HaplotypeCaller” commands (Cornish and Guda, 2015).

Sequencing of PCR Products Using the MinION Sequencer

To enable use of multiplex MinION sequencing at reduced cost with shorter analysis times, each specific barcode tag for the MinION was added at the start of each reverse primer for eight viral RNA segments (influenza A, Uni13/Inf1-Bc; influenza B, UniR-Bc) (Supplementary Table S1). An additional nested PCR assay with custom barcoded primers was performed consecutively to match the sequence reads to the correct clinical specimens. Briefly, 1 μL of 10 ng from each purified PCR amplicon of eight viral RNA segments was used as template DNA. The nested PCR assays were performed using KAPA HiFi HotStart ReadyMix PCR Kit (Kapa Biosystems, Wilmington, MA, United States) and thermal cycling was carried out under the following conditions: 94°C for 5 min, followed by 18 cycles at 98°C for 10 s, 60°C for 30 s, 72°C for 120 s, and a final extension at 72°C for 5 min. The nested PCR products with each ONT-specific barcode tag were analyzed using the Agilent 2100 Bioanalyzer with an Agilent DNA 7500 Kit. Each PCR product was consequently re-purified with Agencourt AMPure XP according to the manufacturer’s instructions. The nested PCR amplicons (300 ng each) were measured by Qubit 3.0 Fluorometer (Thermo Fisher Scientific) and simultaneously processed for library preparation using the SQK-LSK308 1D2 Sequencing Kit (ONT) according to the manufacturer’s instructions. After the MinION Platform QC run, the mixed DNA library derived from 13 clinical samples was loaded into a single MinION Flow Cell (FLO-MIN107 R9.5 Version) and the “NC_48Hr_sequencing_FLO-MIN107_SQK-LSK308_plus_basecaller” protocol was initiated using MinKNOW software (v 1.4.1, ONT).

Analysis of MinION Data

Local basecalling was performed using MinKNOW automatically in real time. 1D2 reads were generated from FAST5 reads by Albacore (v 2.2.1). Low-quality reads and < 500- and > 3000-length reads were removed by SeqKit (v 0.8.1). Demultiplexing and adapter trimming was performed using Porechop (v 0.2.2), and all ONT-barcode data were eventually collected and mapped by Minimap2 (v 0.7.15) using the reference sequences for influenza A H3N2 (GenBank; CY259943-50) and influenza B Yamagata lineage (GenBank; MH233733-40). Single nucleotide variant calling from mapped data was obtained by Samtools and BCFtools (v 1.5.0) via the mpileup command, and variant calling under 10× read depth was filtered. Mapped data were visualized by IGV software (v 2.3.8) and analyzed by Samtools and Qualimap (v 2.2.1).

Results

Multiplex RT-PCR for Clinical Samples

All eight segments of influenza gene from each nasopharyngeal swab samples were successfully amplified with multiplex RT-PCR. Based on the results of RT-PCR amplicon sequencing using the Illumina MiSeq, 6 samples were identified as influenza A H3N2 subtype and 7 samples were identified as influenza B Yamagata lineage. The complete influenza virus genomes were deposited in the Global Initiative on Sharing All Influenza Data database (GISAID: EPI_ISL_313902 - 313914).

Results of MinION Sequencing, Mapped for Reference Genomes

Amplicons with ONT-specific barcodes were successfully amplified by nested RT-PCR. Finally, a total of 2,120,657 1D reads (total 3.2 Gbases) were generated within 48 h after the start of MinION sequencing, and a total of 71,066 1D2 reads (total 101 Mbases) were generated by Albacore. Among the total collected reads, 1.87% ± 1.42% and 5.97% ± 2.83% were classified into ONT-barcodes after barcode demultiplexing as 1D and 1D2, respectively. Tables 1, 2, and Supplementary Table S2 show the results of mapped data for reference sequences by Minimap2, and Figures 1, 2 show the depth of coverage for the reference sequences per sample. The mean coverage per sample for influenza A virus was 1,578 and 89 and that for influenza B virus was 6,545 and 86 for 1D and 1D2, respectively. The raw sequence data of the Illumina MiSeq and MinION sequencers were deposited in the Sequence Read Archive (SRA: PRJNA491222).

TABLE 1
www.frontiersin.org

TABLE 1. Results of whole genome sequencing of influenza A virus using the MinION sequencer.

TABLE 2
www.frontiersin.org

TABLE 2. Results of whole genome sequencing of influenza B virus by MinION sequencer.

FIGURE 1
www.frontiersin.org

FIGURE 1. Coverage of MinION sequencing reads for a reference genome of influenza A virus. The upper graph shows the results of 1D reads and the lower graph shows results of 1D2 reads.

FIGURE 2
www.frontiersin.org

FIGURE 2. Coverage of MinION sequencing reads for a reference genome of influenza B virus. The upper graph shows the results of 1D reads and the lower graph shows results of 1D2 reads.

Comparison With the Illumina MiSeq Platform

Tables 1, 2, and Supplementary Table S3 show the nucleotide sequence identity between influenza gene segments sequenced using the MinION and Illumina MiSeq sequencers. The total accuracy rate of nucleotide sequence identity for influenza A virus when using the MinION sequencer was 99.94 and 99.98% and that for influenza B virus was 99.95 and 99.98% for 1D and 1D2, respectively. Variant calling from 1D2 rather than from 1D reads decreased the numbers of false positives and false negatives when using the MinION sequencer. The total recall and precision rates were 89.41% and 97.88% from 1D reads and 93.28 and 99.86% from 1D2 reads, respectively (Supplementary Table S4). The accuracy of the read level of MinION sequencing, which was calculated based on the mapped data for consensus influenza genome sequences generated from Illumina MiSeq reads, was 88.74% for 1D reads and 92.81% for 1D2 reads (Supplementary Table S5).

Discussion

In this pilot study, we proposed a high-throughput MinION sequencing protocol that can be performed both in the field and in the clinical setting. We demonstrated that multiplex MinION sequencing could achieve variant calling easily and rapidly from 1D2 reads generated by combining multiplex RT-PCR methods for both influenza A and B viruses. In particular, sample processing and library preparation processes can be completed within 12 h without expensive or large experimental devices. Furthermore, the sequencing data can be analyzed in real time using a laptop computer. Moreover, samples can be prepared and sequenced without storage or transfer to a fully equipped laboratory. Therefore, our method would have the potential to allow real-time monitoring of both genetic rearrangements and the emergence of new influenza virus strains.

For accurate surveillance and vaccine subpopulation selection, a highly precise sequencing and variant calling method is always required. In particular, the recall and precision rate of the current MinION sequence technology were inferior to those of the Illumina MiSeq platform, and these limitations were similar to findings from other studies on genotyping using the MinION sequencer (Wang et al., 2015; Cornelis et al., 2017; Lang et al., 2018; Runtuwene et al., 2018). Thus, our method is currently unsuitable to replace the Illumina platform. The high error rate and number of coverages of the MinION sequencer could affect the accuracy of variant calling because the read level accuracy rate was estimated to be about 88.74% for 1D and 92.81% for 1D2. In fact, false-positive results can be significantly reduced via high-quality 1D2 reads (precision rate: 97.88% from 1D reads vs. 99.86% from 1D2 reads). However, 1D2 reads would unfortunately still include false negatives due to single nucleotide polymorphisms (recall rate: 89.41% from 1D reads vs. 93.28% from 1D2 reads) because they would include deletion errors more frequently around the homopolymer sequence and have lower coverage than 1D reads. Further improvements in sequence accuracy of the MinION sequencer would make it possible to perform surveillance and vaccine subpopulation selection.

In this study, we used the clinically applicable amplicon sequencing protocol, to analyze influenza virus genome sequencing. Previously, Greninger et al. (2015) demonstrated that metagenomic sequencing analysis with the MinION sequencer could detect unknown emerging viruses without prior information and would be desirable for point-of-care genomic sequencing even in the developing world. However, its sequencing scheme still has lower sensitivity and specificity levels compared with the amplicon sequencing protocol of clinical samples, because of contamination from the large amount of host genome fragments in clinical samples (Rosseel et al., 2012; Quick et al., 2017). Moreover, there is currently no practical preparation kit for adding barcode tags in multiple-sample processing based on metagenomic sequencing with the 1D2 sequencing kit, even though multiplex sequencing is an essential technology that reduces the cost and time of analysis. Amplicon sequencing can solve one of these problems precisely by reducing contamination from host genome fragments via the added multiplex RT-PCR step. Furthermore, it is possible to simply add a barcode tag using nested PCR, thus enabling multiple sample processing. Of course, the DNA fragmentation step, including mechanical or chemical pretreatment, is required for the Illumina MiSeq platform because of its read length of up to 600 bp. However, a fragmentation step and additional device for amplicon sequencing are not necessary when using the MinION sequencer. This is because it can obtain long read data over 100 kb, which is another one of its advantageous characteristics (Jain et al., 2018). Unfortunately, even with high-fidelity DNA polymerase reactions, amplification by multiple RT-PCR and nested PCR could be a drawback as it may cause PCR-generated errors (Brandariz-Fontes et al., 2015; Potapov and Ong, 2017).

Multiplex sequencing with introduced barcode tag nucleotides would also contribute to improving the accuracy of pairing 1D sequence reads. In this study, we used nested PCR protocols to add 30-base ONT-specific barcodes against multiplex RT-PCR products of eight viral RNA segments and analyzed multiple samples simultaneously. In fact, in 1D2 sequencing using PCR amplicons, false-positive pairings between different samples could occur because 1D2 reads are generated with aberrant in silico chimeric reads from 1D reads. However, the addition of 30-base ONT-specific barcodes can effectively decrease the risk of false 1D2 pairing of strands.

Avian influenza types, such as influenza A virus subtypes H5N1 and H7N9, are likely to emerge with the potential to cause pandemics in developing countries (Oshitani et al., 2008; Fedson, 2009). Furthermore, recent meta-analyses showed that morbidity and mortality are higher in patients with influenza viruses in developing countries than in those in developed countries (Nair et al., 2011; Ortiz et al., 2013). The WHO indicated an urgent need for strengthening local responses to emerging influenza viruses [World Health Organization (WHO), 2005]. Whole genome analyses of influenza viruses using NGS is important for the surveillance of novel strains of influenza virus with genetic rearrangements and antiviral resistance. Compared with other NGS technologies, the MinION sequencer has excellent portability in addition to reduced initial installation costs, running costs, and analysis times. In the near future, the low-priced Flongle Flow Cell sequencing device2 – an adapter for MinION flow cells that has 126 channels compared with 512 channels in the current MinION – will become available and provide further versatility. Therefore, this novel procedure has the potential to be applicable as a surveillance technology for emerging influenza strains and for use even in situations where research facilities are inadequate.

Conclusion

We developed novel high-throughput MinION sequencing methods for real-time sequencing of influenza A and B viruses using a portable sequencing setup with quick and easy sample preparation. Further improvement of read level accuracy and analytical tools is necessary to achieve better utilization of the MinION sequencer for both real-time monitoring of genetic rearrangements and evaluation of newly emerging viruses.

Author Contributions

KT, YN, MO, and TMa designed the research. KI, ToT, MT, EN, TaT, and NT performed the research. KM, SM, and TMu provided scientific guidance. KI and TMa prepared the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Mr. Makio, Mr. Takao, Mr. Ogasawara, Mr. Takano, and Mr. Miyazawa, medical technologists of the Japan-Ground Self-Defense Forces, for their technical assistance with the NGS experiments, and Mrs. Kobayashi and Mrs. Saito, nurses of the Japan-Ground Self-Defense Forces, for their assistance in collecting samples.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02748/full#supplementary-material

Footnotes

  1. ^https://nanoporetech.com/
  2. ^https://nanoporetech.com/products/flongle

References

Brandariz-Fontes, C., Camacho-Sanchez, M., Vila, C., Vega-Pla, J. L., Rico, C., and Leonard, J. A. (2015). Effect of the enzyme and PCR conditions on the quality of high-throughput DNA sequencing results. Sci. Rep. 5:8056. doi: 10.1038/srep08056

PubMed Abstract | CrossRef Full Text | Google Scholar

Briand, F. X., Schmitz, A., Ogor, K., Le Prioux, A., Guillou-Cloarec, C., Guillemoto, C., et al. (2017). Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: phylogenetic analyses and markers for zoonotic potential. Euro. Surveill. 22:30473. doi: 10.2807/1560-7917.ES.2017.22.9.30473

PubMed Abstract | CrossRef Full Text | Google Scholar

Cornelis, S., Gansemans, Y., Deleye, L., Deforce, D., and Van Nieuwerburgh, F. (2017). Forensic SNP genotyping using nanopore MinION sequencing. Sci. Rep. 7:41759. doi: 10.1038/srep41759

PubMed Abstract | CrossRef Full Text | Google Scholar

Cornish, A., and Guda, C. (2015). A comparison of variant calling pipelines using genome in a bottle as a reference. Biomed. Res. Int. 2015:456479. doi: 10.1155/2015/456479

PubMed Abstract | CrossRef Full Text | Google Scholar

Fedson, D. S. (2009). Meeting the challenge of influenza pandemic preparedness in developing countries. Emerg. Infect. Dis. 15, 365–371. doi: 10.3201/eid1503.080857

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghedin, E., Sengamalay, N. A., Shumway, M., Zaborsky, J., Feldblyum, T., Subbu, V., et al. (2005). Large-scale sequencing of human influenza reveals the dynamic nature of viral genome evolution. Nature 437, 1162–1166. doi: 10.1038/nature04239

PubMed Abstract | CrossRef Full Text | Google Scholar

Grenfell, B. T., Pybus, O. G., Gog, J. R., Wood, J. L., Daly, J. M., Mumford, J. A., et al. (2004). Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303, 327–332. doi: 10.1126/science.1090727

PubMed Abstract | CrossRef Full Text | Google Scholar

Greninger, A. L., Naccache, S. N., Federman, S., Yu, G., Mbala, P., Bres, V., et al. (2015). Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis. Genome Med. 7:99. doi: 10.1186/s13073-015-0220-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Harvala, H., Frampton, D., Grant, P., Raffle, J., Ferns, R. B., Kozlakidis, Z., et al. (2017). Emergence of a novel subclade of influenza A(H3N2) virus in London, December 2016 to January 2017. Euro. Surveill. 22:30466. doi: 10.2807/1560-7917.ES.2017.22.8.30466

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoenen, T., Groseth, A., Rosenke, K., Fischer, R. J., Hoenen, A., Judson, S. D., et al. (2016). Nanopore sequencing as a rapidly deployable ebola outbreak tool. Emerg. Infect. Dis. 22, 331–334. doi: 10.3201/eid2202.151796

PubMed Abstract | CrossRef Full Text | Google Scholar

Houlihan, C., Frampton, D., Ferns, R. B., Raffle, J., Grant, P., Reidy, M., et al. (2018). The use of whole genome sequencing in the investigation of a nosocomial influenza virus outbreak. J. Infect. Dis. 218, 1485–1489. doi: 10.1093/infdis/jiy335

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, M., Koren, S., Miga, K. H., Quick, J., Rand, A. C., Sasani, T. A., et al. (2018). Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat. Biotechnol. 36, 338–345. doi: 10.1038/nbt.4060

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, K., Surendranath, V., Quenzel, P., Schofl, G., Schmidt, A. H., and Lange, V. (2018). Full-length HLA class I genotyping with the MinION nanopore sequencer. Methods Mol. Biol. 1802, 155–162. doi: 10.1007/978-1-4939-8546-3_10

PubMed Abstract | CrossRef Full Text | Google Scholar

MacFadden, D. R., Mcgeer, A., Athey, T., Perusini, S., Olsha, R., Li, A., et al. (2018). Use of genome sequencing to define institutional influenza outbreaks, Toronto, Ontario, Canada, 2014-15. Emerg. Infect. Dis. 24, 492–497. doi: 10.3201/eid2403.171499

PubMed Abstract | CrossRef Full Text | Google Scholar

Nair, H., Brooks, W. A., Katz, M., Roca, A., Berkley, J. A., Madhi, S. A., et al. (2011). Global burden of respiratory infections due to seasonal influenza in young children: a systematic review and meta-analysis. Lancet 378, 1917–1930. doi: 10.1016/S0140-6736(11)61051-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ortiz, J. R., Jacob, S. T., and West, T. E. (2013). Clinical care for severe influenza and other severe illness in resource-limited settings: the need for evidence and guidelines. Influenza Other Respir. Viruses 7(Suppl. 2), 87–92. doi: 10.1111/irv.12086

PubMed Abstract | CrossRef Full Text | Google Scholar

Oshitani, H., Kamigaki, T., and Suzuki, A. (2008). Major issues and challenges of influenza pandemic preparedness in developing countries. Emerg. Infect. Dis. 14, 875–880. doi: 10.3201/eid1406.070839

PubMed Abstract | CrossRef Full Text | Google Scholar

Potapov, V., and Ong, J. L. (2017). Examining sources of error in PCR by single-molecule sequencing. PLoS One 12:e0169774. doi: 10.1371/journal.pone.0169774

PubMed Abstract | CrossRef Full Text | Google Scholar

Quick, J., Ashton, P., Calus, S., Chatt, C., Gossain, S., Hawker, J., et al. (2015). Rapid draft sequencing and real-time nanopore sequencing in a hospital outbreak of Salmonella. Genome Biol. 16:114. doi: 10.1186/s13059-015-0677-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Quick, J., Grubaugh, N. D., Pullan, S. T., Claro, I. M., Smith, A. D., Gangavarapu, K., et al. (2017). Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples. Nat. Protoc. 12, 1261–1276. doi: 10.1038/nprot.2017.066

PubMed Abstract | CrossRef Full Text | Google Scholar

Quick, J., Loman, N. J., Duraffour, S., Simpson, J. T., Severi, E., Cowley, L., et al. (2016). Real-time, portable genome sequencing for Ebola surveillance. Nature 530, 228–232. doi: 10.1038/nature16996

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosseel, T., Scheuch, M., Hoper, D., De Regge, N., Caij, A. B., Vandenbussche, F., et al. (2012). DNase SISPA-next generation sequencing confirms Schmallenberg virus in Belgian field samples and identifies genetic variation in Europe. PLoS One 7:e41967. doi: 10.1371/journal.pone.0041967

PubMed Abstract | CrossRef Full Text | Google Scholar

Runtuwene, L. R., Tuda, J. S. B., Mongan, A. E., Makalowski, W., Frith, M. C., Imwong, M., et al. (2018). Nanopore sequencing of drug-resistance-associated genes in malaria parasites, Plasmodium falciparum. Sci. Rep. 8:8286. doi: 10.1038/s41598-018-26334-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, Y., and McCauley, J. (2017). GISAID: global initiative on sharing all influenza data - from vision to reality. Euro. Surveill. 22:30494. doi: 10.2807/1560-7917.ES.2017.22.13.30494

PubMed Abstract | CrossRef Full Text | Google Scholar

Taubenberger, J. K., and Morens, D. M. (2008). The pathology of influenza virus infections. Annu. Rev. Pathol. 3, 499–522. doi: 10.1146/annurev.pathmechdis.3.121806.154316

CrossRef Full Text | Google Scholar

Vemula, S. V., Zhao, J., Liu, J., Wang, X., Biswas, S., and Hewlett, I. (2016). Current approaches for diagnosis of influenza virus infections in humans. Viruses 8:96. doi: 10.3390/v8040096

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Moore, N. E., Deng, Y. M., Eccles, D. A., and Hall, R. J. (2015). MinION nanopore sequencing of an influenza genome. Front. Microbiol. 6:766. doi: 10.3389/fmicb.2015.00766

CrossRef Full Text | Google Scholar

World Health Organization (WHO) (2005). utbreak Communication : Best Practices for Communicating with the Public During an Outbreak: Report of the WHO Expert Consultation on Outbreak Communications Held in Singapore, 21-23 September 2004. Geneva: WHO/CDS.

World Health Organization (WHO) (2017). WHO Information for the Molecular Detection of Influenza Viruses - Update. Geneva: WHO.

Zhao, J., Liu, J., Vemula, S. V., Lin, C., Tan, J., Ragupathy, V., et al. (2016). Sensitive detection and simultaneous discrimination of influenza A and B viruses in nasopharyngeal swabs in a single assay using next-generation sequencing-based diagnostics. PLoS One 11:e0163175. doi: 10.1371/journal.pone.0163175

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, B., Donnelly, M. E., Scholes, D. T., St George, K., Hatta, M., Kawaoka, Y., et al. (2009). Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses. J. Virol. 83, 10309–10313. doi: 10.1128/JVI.01109-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, B., Lin, X., Wang, W., Halpin, R. A., Bera, J., Stockwell, T. B., et al. (2014). Universal influenza B virus genomic amplification facilitates sequencing, diagnostics, and reverse genetics. J. Clin. Microbiol. 52, 1330–1337. doi: 10.1128/JCM.03265-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: influenza virus, MinION, MiSeq, whole genome sequencing, Japan

Citation: Imai K, Tamura K, Tanigaki T, Takizawa M, Nakayama E, Taniguchi T, Okamoto M, Nishiyama Y, Tarumoto N, Mitsutake K, Murakami T, Maesaki S and Maeda T (2018) Whole Genome Sequencing of Influenza A and B Viruses With the MinION Sequencer in the Clinical Setting: A Pilot Study. Front. Microbiol. 9:2748. doi: 10.3389/fmicb.2018.02748

Received: 01 August 2018; Accepted: 26 October 2018;
Published: 13 November 2018.

Edited by:

Leonard Peruski, Centers for Disease Control and Prevention (CDC), United States

Reviewed by:

Alfonso Benítez-Páez, Instituto de Agroquímica y Tecnología de Alimentos (IATA), Spain
Zheng Wang, Yale University, United States

Copyright © 2018 Imai, Tamura, Tanigaki, Takizawa, Nakayama, Taniguchi, Okamoto, Nishiyama, Tarumoto, Mitsutake, Murakami, Maesaki and Maeda. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Takuya Maeda, t_maeda@saitama-med.ac.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.