key: cord-0786748-yayaq81u authors: Ciuffreda, Laura; Rodríguez-Pérez, Héctor; Flores, Carlos title: Nanopore sequencing and its application to the study of microbial communities date: 2021-03-07 journal: Comput Struct Biotechnol J DOI: 10.1016/j.csbj.2021.02.020 sha: c682cdf534bc2a2a2c45b68d5005929e57b3b308 doc_id: 786748 cord_uid: yayaq81u Since its introduction, nanopore sequencing has enhanced our ability to study complex microbial samples through the possibility to sequence long reads in real time using inexpensive and portable technologies. The use of long reads has allowed to address several previously unsolved issues in the field, such as the resolution of complex genomic structures, and facilitated the access to metagenome assembled genomes (MAGs). Furthermore, the low cost and portability of platforms together with the development of rapid protocols and analysis pipelines have featured nanopore technology as an attractive and ever-growing tool for real-time in-field sequencing for environmental microbial analysis. This review provides an up-to-date summary of the experimental protocols and bioinformatic tools for the study of microbial communities using nanopore sequencing, highlighting the most important and recent research in the field with a major focus on infectious diseases. An overview of the main approaches including targeted and shotgun approaches, metatranscriptomics, epigenomics, and epitranscriptomics is provided, together with an outlook to the major challenges and perspectives over the use of this technology for microbial studies. The study of microbial communities, including bacteria, viruses, archaea and fungi, is crucial to understand important aspects of the environment and/or human health. Over the last decades, there has been an important shift in the way microbial communities have been explored due to the introduction of sequencing technologies. In the late '80s, scientists realized that the world of uncultured microorganisms outsized the cultured world, and the analysis of DNA sequences replaced culturing for the study of complex microbial communities [1, 2] . Since then, the advent of nextgeneration sequencing (NGS) technologies has unquestionably led to a real revolution in the area of microbiology, with major breakthroughs the complete characterization of the human gut microbiome [3] or the identification of novel phyla with undiscovered biology [4] , to name a few. The introduction of third-generation sequencing represented another major turning point in the field, because it opened the possibility for real-time sequencing of long reads. The marketed technologies that currently dominate this field are single-molecule real-time sequencing (SMRT, commercialized by Pacific Biosciences) and nanopore sequencing (NS, commercialized by Oxford Nanopore Technologies (ONT)). ONT nanopore sequencing allows single molecule sequencing based on bioengineered nanopores, which are embedded into an electrically resistant membrane where a voltage is applied. When a single stranded DNA/RNA fragment passes through the nanopore, it causes a change in electrical current through the membrane that is translated into a specific sequence of nucleotides using recurrent neural network (RNN)-based algorithms. The first marketed device to sequence by nanopores, the MinION sequencer, was introduced in 2014 and is a small and inexpensive device (starting pack available for 1000$) [5] . A single flow cell of MinION contains up to 2,048 nanopores which are controlled in groups of 512 via an application specific integrated circuit (ASIC). It can be directly plugged into a portable computer allowing real-time acquisition and analysis of data, making it the first sequencer to enable in-field sample genomic characterization [6] . Later on, other three devices were introduced by ONT, the GridION and the PromethION, which allow parallel running of up to five and 48 flow cells, respectively, and Flongle, a smaller adaptor device for running smaller experiments both on MinION and GridION. Since nanopore sequencing does not need to occur in amplified DNA, very long reads (more than 2 Mb [7] ) can be generated, with no theoretical limit in read length [8] . The possibility of obtaining these very long reads has allowed significant improvements in diverse applications such as de novo genome assembly [9] , and in the deeper characterization of repetitive DNA elements [10, 11] , which could not be resolved relying only on existing short-read NGS technologies. Furthermore, because sequencing is mediated by the translation of an electrical signal into a sequence of nucleotides, nanopore sequencing allows the identification of native base modifications [12] and the direct sequencing of RNA molecules [13] . The main drawback of this technology has been the high error rate [14] . However, this has been continuously improved over the years, with a current modal raw read accuracy of 97% with the latest released flow cell [15] . Furthermore, the continuous development of novel basecallers and bioinformatic tools for read error correction (polishing) and consensus generation has dramatically helped to further improve this aspect. As a consequence, the impact of nanopore read errors on taxonomic classification and other microbial analyses is limited. Sequencing throughput is also very variable. At the moment, around 30 Gb of data can theoretically be generated using a MinION flow cell, corresponding to the sequencing of 25 E. coli genomes with a coverage of 100Â, and up to 200 Gb of data per flow cell using a PromethION system, allowing to sequence a minimum of one human genome with a coverage of around 30Â. Thus, through parallel run of 48 flow cells, the acquisition of around 9600 Gb of DNA/RNA data is theoretically possible at the moment, corresponding to the sequencing of a minimum of 48 human genomes in a single run. Box 1 Basic concepts. Taxonomic profiling: type of analysis aiming at the identification of taxa present in a sample together with their relative abundances. It is typically performed using marker genes which enable the discrimination between different taxa. It answers the question: ''who is there?" Functional analysis: the study of the metabolic and other biological pathways related to the taxa present in a sample. It is usually performed by comparing gene sequences to functional databases. It answers the question: ''what are they doing?" Operational Taxonomic Unit (OTU): a group of closely related individuals which are arranged together based on the similarity of specific sequences (usually the 16S rRNA gene). Basecalling: computational process to assign nucleotides to sequence from the raw electric signal data (squiggle) generated by a nanopore sequencing device. Quality control: set of read filtering steps prior to analysis which usually consist of read-length and read quality filtering. Contig: a series of overlapping sequences used to reconstruct the original DNA sequence of a genomic region. Polishing: it refers to the analytical process aiming to improve the base accuracy of a contig. k-mer: subsequences of length k that are contained within a sequence. k-mer frequencies encompass features that are characteristic of particular organisms and are suitable for taxonomic binning or microbial composition inference. Genome or metagenomic assembly: computational process for reconstructing individual contigs from a genomic or metagenomic dataset that spans complete or nearlycomplete microbial genes and genomes. Assembly graph: graph representation of the final assembly of a genome or metagenome based on read overlap or kmer information and including all possible paths for contig reconstruction. Metataxonomics is defined as the process of characterization of the microbiota through the amplification and sequencing of conserved marker genes followed by the assignment of the generated reads to specific taxonomic levels, and/or the construction of a taxonomic tree [16] . Metataxonomics is commonly differentiated from metagenomics, which is defined as the study of the totality of genomes of the microbiota through shotgun sequencing techniques. An explanation of the basic concepts related to microbial study and bioinformatic analysis is given in Box 1. A general workflow of experimental and bioinformatics steps for targeted-and metagenomics studies is presented in Fig. 1 . The most commonly used marker genes in metataxonomics are the 16S rRNA for bacteria and archaea, and the 18S rRNA for fungi, while there are no marker genes for viruses, even though they are an integral part of the microbiota. The 16S rRNA gene has been historically used in the taxonomic classification of known and new microbial taxa and in phylogenetic analysis [17, 18] . The 16S rRNA gene spans~1500 bp including nine hypervariable segments (V1-V9) flanked by highly conserved regions which can be targeted by universal primers for amplification. The variable regions allow to distinguish between different taxa. A similarity threshold of 98.65% in the 16S rRNA gene sequences has been identified to discriminate between two species [19] . These characteristics make the 16S rRNA gene the gold standard for microbial classification [20] . Compared to the NGS technologies, where only short segments of the gene encompassing one or several hypervariable regions can be targeted, long-read sequencing allows amplification and analysis of the full-length 16S rRNA gene, which provides a more realistic representation of the taxa in a sample [21] . In fact, despite the higher error rate characterizing nanopore sequencing, the increased read length achieved through the full-length 16S rRNA gene sequencing allows species-level classification, improving taxa resolution over previous technologies [22] [23] [24] [25] . Furthermore, the use of clustering-based algorithms (e.g. NanoCLUST [25] ) and other polishing tools allows overcoming the read error rate, with the generation of highly accurate consensus sequences of the 16S rRNA genes that are now able to discriminate between species. ONT allows sequencing with portable devices, opening the possibility for the use of 16S rRNA gene sequencing for rapid infield pathogen detection. In particular, ONT has developed two 16S barcoding kits for rapid sequencing of full-length 16S rRNA genes, which allow simultaneous sequencing of up to 12 or 24 samples ( Table 1 ). The library preparation protocol for both kits has the advantage to be fast (< 2 h) and easy to perform, and consists of an amplification step where barcodes are added to the amplicons, followed by the attachment of adapters necessary to mediate molecule entrance into the nanopores on the flow cell. Alternatively, when multiplexing of more than 24 samples is needed, it is possible to use the standard PCR barcoding amplicon kit, which currently enables simultaneous sequencing of up to 96 samples in a single experiment. A few considerations need to be made before embarking on a 16S rRNA gene sequencing experiment. First of all, only bacterial and archaea communities are identified using this approach, while viruses and fungi are missed. Alternatively, the mycobiome (the fungal microbiota) can be studied using the 18S rRNA gene or the internal transcribed spacer (ITS) as markers. However, ONT kits specifically targeting these regions have not yet been developed, and this may have limited the study of the fungal communities using this technology. Secondly, although the PCR step during library preparation increases the chances to detect low abundant taxa, it is well-known that the PCR introduces biases in taxonomic classification and estimation of relative abundances [26] . In fact, Kai et al. reported that Bifidobacterium is not detected by the ONT 16S rRNA library kit, due to the lack of annealing of the universal primers to the flanking regions of the 16S rRNA gene of this taxon [27] . They further addressed this bias by changing the reverse primer sequence to target all taxa present in their sample [28] . An additional source of bias to consider is the differential number of rrn operons in the genomes of different taxa, which often leads to inaccuracies in the estimation of the abundance profiles. Although algorithms for 16S rRNA gene copy number normalization have often been used to overcome this bias, it has recently been proved that they fail to provide a more reliable picture of the community composition in metataxonomics studies [29] . There are a broad range of bioinformatics pipelines for metagenomics which are also valid for metataxonomic analysis ( Table 2) , including the widely used multi-purpose pipelines based on Operational Taxonomic Unit (OTU) picking and/or Amplicon Sequence Variants (ASV) analysis [30, 31] . Most of these were initially developed to work with short-read data (particularly for those from Illumina) and are not suitable for nanopore read lengths and error profiles, commonly leading to issues such as an overestimation of taxa diversity. Hence, the potential benefits of performing taxonomic classification with full-length 16S rRNA reads have not been extensively explored. Therefore, the rapid changes in the sequencing technology have outpaced the availability of specific tools and benchmark studies of nanopore 16S rRNA reads. Computational methods for nanopore 16S rRNA analysis have been reviewed recently [32] . Depending on the chosen approach to classify the sequences, read classification techniques can be categorized into alignment-based and alignment-free methods. EPI2ME (ONT) is the most extensively used analysis pipeline for nanopore 16S rRNA. It covers end-to-end analysis of nanopore 16S rRNA data in a cloud-based environment and includes demultiplexing, quality filtering and taxonomic assignment using the BLAST tool against the NCBI database. The main drawback of using EPI2ME is its limited possibility to customize workflow parameters, such as reference databases and alignment options. Furthermore, this tool can only be accessed by ONT customers through a web application and the output data format is incompatible with other software for downstream analysis, highlighting the need for the development of alternatives based on available opensource tools. Simpler workflows proceed with the alignment of the input sequences using tools designed to work efficiently with long noisy reads, like minimap2 [33] , against specifically designed 16S rRNA databases which contain only curated 16S rRNA sequences and their taxonomy (Table 3) . Additionally, alignmentfree methods like Centrifuge [34] or Kraken [35, 64] (described in more detail in following sections) have emerged as feasible options for taxonomic classification of nanopore 16S rRNA reads. A recently proposed approach, NanoCLUST, relies on the Uniform Manifold Approximation and Projection (UMAP) algorithm to cluster full-length 16S rRNA reads and then classifies a representative polished sequence from each cluster to deliver abundance profiles at different taxonomic levels [25] . ONT nanopore sequencing has opened the possibility for the first time to sequence and analyse data in real-time at competitive costs. A major application in hospital settings is the rapid diagnostics of infectious diseases to allow prompt patient management and appropriate treatment decisions. A summary of the successful clinical applications of ONT nanopore sequencing in infectious diseases is given in Table 4 . Proof-of-principle studies have been conducted leveraging ONT 16S rRNA targeted sequencing for the diagnosis of bacterial infections. In an early study, Mitsuhashi et al. developed a protocol for rapid characterization of bacterial composition using a mock community, which was then evaluated on a pleural effusion sample from a patient with empyema [36] . The protocol was based on 16S rRNA gene sequencing followed by analysis with BLAST-based searching or Centrifuge classification. They efficiently characterized the mock community at the species level using BLAST against the GenomeSync database, while Centrifuge missed the identification of one of the species. By comparing results at different times during the sequencing run, they concluded that 5 min of sequencing time was enough to obtain a sufficient amount of data to identify all bacteria taxa in the sample, and to reach a sensitivity >90% using BLAST. However, when the rapid protocol was tested on the clinical sample, a longer sequencing time was necessary to identify low abundant taxa. Later studies adopted a similar analysis protocol (replacing BLAST by minimap2) to conduct point-of-care diagnostics, both in developed [37] and in resource-poor countries [38] . In particular, this last study was conducted on cerebrospinal fluid samples from eleven patients with bacterial meningitis in Zambia, where a portable sequencing-based system was set up. Although they could successfully confirm the results from culture-based methods on four samples, two positive samples showed different bacterial compositions using the MinION sequencer, and these results were Table 2 Main long-read bioinformatics tools for targeted and shotgun approaches. Aligners/Alignment-based classifiers BLAST, MEGABLAST [58, 59] Targeted; Shotgun Gold-standard alignment tools for classification of nucleotide and protein sequences. Feature webbased version and multiple implementations for specific purposes. minimap2 [33] Targeted; Shotgun Versatile tool for fast read alignments against large reference databases. Kraken, Kraken2 [35, 64] Targeted; Shotgun Taxonomic classification tool implementing an accurate and fast k-mer matching. KrakenUniq [65] Shotgun Classifier that combines Kraken classification tool with the assessment of the coverage of unique kmers for better recall and precision. Bracken [66] Targeted; Shotgun Relative abundance estimation tool for single-level abundance using Kraken read classification output. Metamaps [69] Shotgun Read assignment and sample composition estimation for nanopore metagenomic datasets. Centrifuge [34] Targeted; Shotgun Read classification based on the Burrows-Wheeler transform (BWT) and the Ferragina-Manzini (FM) index that performs fast classification relying on small pre-computed index databases. Mash [72] Targeted; Shotgun Fast genome and metagenome distance estimation tool that computes distances between sequences using the MinHash algorithm. Long-read assemblers Canu [90] Shotgun Assembly pipeline for long-reads that compute and process read overlaps for the generation of contigs and draft assemblies. miniasm [73] Shotgun Fast OLC-based assembler for long reads that builds assembly graphs from all-vs-all read mappings. wtdbg2 [91] Shotgun De-novo sequence assembler for uncorrected long-reads based on Fuzzy Bruijn graphs to compute contigs. OPERA-MS [95] Shotgun Hybrid metagenomic assembler that first performs a short-read assembly and then maps short and long reads to resolve contiguity of contigs. MetaFlye [96] Shotgun Metagenomic assembler from the Flye package featuring repeat graphs to compute high-quality metagenome assemblies. MetaSPAdes [74] Shotgun Metagenomic assembly module from SPAdes assembler that features a hybrid assembly option. Shotgun Snakemake-based pipeline for assembly and polishing of long genomes from long nanopore reads. Metagenomics workflow for hybrid assembly, differential coverage binning, transcriptomics and pathway analysis. NanoSPC [71] Shotgun Metagenomic analysis pipeline that includes viral and bacterial pathogen identification, genome assembly and variant calling. BusyBee https://ccb-microbe.cs.unisaarland.de/busybee/ Shotgun Web-based metagenomic analysis pipeline for long-reads and contigs that features taxonomic and functional annotation of AMR elements along with a comprehensive visualization of results. not confirmed using any additional method. Furthermore, culturenegative samples were positive using MinION, suggesting higher sensitivity of the sequencing method compared to the traditional culturing. However, no further validation with complementary methods was carried out. For the adoption of the MinION sequencer for rapid diagnostics of infectious diseases in clinical settings, it is necessary for future studies to be conducted with standardized protocols on large sam-ple size, followed by appropriate cross-validation of the results. This issue was partially addressed by Neuenschwander et al. who developed LORCAN, a standardized laboratory protocol and automated pipeline for taxonomic identification of bacterial mixtures based on the ONT 16S rRNA sequencing [39] . The library consists of PCR amplicons from regions of the 16S rRNA gene (length between 500 and 1000 bp) and the pipeline generates consensus sequences that improve the accuracy of taxonomic classification. The workflow was tested on culture isolates and artificial mock communities from read or amplicon mixing. Analysis of samples using LORCAN generated consensus sequences with a sequence identity of 99.6% to their corresponding Sanger-obtained sequences and a turnaround time from raw amplicons to reports of about 8 h. While this workflow has the potential to be used in the clinics, it has not been tested yet on clinical samples, leaving open the question of whether it could truly be adopted in real settings. Metagenomics allows deep investigations of microbial communities, usually encompassing taxonomic analysis, functional profiling and whole-genome assembly. In contrast to targeted approaches, the entire genomic content of a microbial sample is sequenced, providing greater genomic information [40] . In fact, strain-level information can be accurately recovered with the metagenomic approach, and helps to find the associations between phylogeny and function [41] . Furthermore, metagenomics enables the identification of antimicrobial resistance (AMR) genes or other virulence elements that can be assigned to specific pathogens through genome assembly. It also helps to identify viral communities that are missed using targeted sequencing approaches since their genomes lack consensus sequences that are necessary for universal primer attachments during the library preparation stage [42] . Nanopore sequencing has led to important improvements in the metagenomic analysis because reads can span extended areas of the genome. In fact, long reads can resolve complex genomic structures such as repetitive elements [43] , and allow to identify the position and organization of bacterial pathogenicity islands (PAIs) encoding virulence factors [44] . This long-range genomic information has been particularly relevant for de novo and metagenomics assemblies where it enabled to resolve areas of the genomes where short reads failed. For example, ONT sequencing had allowed to localize AMR gene positions in Klebsiella pneumoniae [45] , to reconstruct plasmid genomes in Enterobacteriaceae isolates [46, 47] , and to circularize the genome of Bordetella pertussis isolates [48] , as well as to assemble the human gut microbiome [43, 49] . The accurate and complete reconstruction of genomes is important for the identification of gene-genome associations in functional and evolutionary studies that, for example, aim to shed light on the genomic organization of metabolic pathways or on the horizontal gene transfer. [82] Bioinformatic resources and database for the annotation of antimicrobial resistance genes (AMR) and mutations from genomic sequences. DeepARG-DB [87] Antibiotic resistance genes database generated by a deep-learning prediction algorithm trained with ARG from other sequence collections. MEGARes [83] Hand-curated database containing AMR genes optimized for use with high-throughput sequencing data. [36] [37] [38] RNA sequencing [118] Rapid identification of pathogens and AMR genes Shotgun metagenomics [97, 98] Surveillance of pathogens and AMR in hospital settings Shotgun metagenomics [99] Genomic surveillance for viral outbreaks RNA sequencing [119, 120, 124, 125] Currently, there are four ONT library preparation kits available for metagenomic studies that differ in the input DNA quantity, preparation time, and throughput ( Table 1 ). The Rapid Sequencing Kit allows fast library preparation (~10 min) and requires an input of > 400 ng of high molecular weight (HMW) genomic DNA (> 30 kb). Multiplexing up to 12 samples can be achieved using the Rapid Barcoding Kit. When a lower starting amount of DNA is available (< 10 ng), the alternative Rapid PCR Barcoding Kit can be used, which includes a PCR step for the amplification of the target DNA and for the attachment of barcodes for multiplexing (up to 12 samples). By using this kit, library preparation requires~15 min added to the time for the PCR step, and read length distribution is centred at around 2 kb. These library preparation kits use a transposase for fragmentation of the HMW genomic DNA and the attachment of transposase adapters which can then be used as anchors for sequencing adapters or as binding sites for PCR primers in the Rapid PCR Barcoding Kit. When higher throughput is needed, ONT suggests the use of other two library preparation kits. The Ligation Sequencing Kit, which requires~1 lg of starting HMW DNA and consists of a protocol lasting~60 min, allowing the maximum throughput achievable while retaining the possibility to call base modifications. DNA molecules are nick-repaired and dA tailed, and then sequencing adapters are ligated onto the prepared ends. This kit can be combined with upstream processes such as target enrichment by capture, size selection, or whole genome amplification (when < 1 ng of original DNA is available), and multiplexing of up to 96 samples can be achieved by using the Native Barcoding Expansion kits. When lower starting material is accessible or sample purity is compromised, the PCR Sequencing Kit can be adopted. In this protocol, the original DNA is fragmented, sheared ends are repaired and dA tailed, adapters containing primer binding sites are ligated and an amplification step follows using primers containing tags for ligase-free attachment of rapid sequencing adapters. The mean read length distribution of the protocol is larger than that obtained by the Rapid PCR Sequencing Kit and is limited only by the processivity of the DNA polymerase. This kit requires around 60 min of hands-on time, added to that for the PCR step, which is variable depending on the number of cycles, template length, and polymerase speed. Multiplexing of up to 12 samples is possible by integrating the PCR Barcoding Kit in the protocol, which uses barcoded primers during the PCR step. The sequencing data yield will depend on the experimental aim and on the sample analysed. Since the recommended coverage varies depending on the aim of the study -the recommended minimal depth of coverage is 10Â for detection of taxa, 20Â for taxonomic assignation and AMR gene analysis, and 30Â for genome assembly [50] -the sequencing data necessary to achieve that specific goal is also variable. As an example, to assemble a genome of 3 Mb, 90 Mb of data from that genome is needed. If 1 Gb of data is sequenced, it will be possible to assemble genomes representing up to 9% of the total data. In addition, the presence of host DNA in the sample is also likely to reduce the amount of data related to the metagenome itself, requiring a larger amount of sequencing data. Host DNA depletion protocols can overcome this issue, which is especially relevant in clinical samples (e.g. respiratory specimens and swabs) where up to 95% of reads could be host-derived [51] . Examples of host depletion protocols include saponin, molYsis kits (Molzym, Germany), or kits for rRNA depletion. In particular, Charalampous et al. have recently optimized a metagenomic protocol for the detection of lower respiratory infections which includes a saponin-based depletion step removing up to the 99.99% of host nucleic acids and enables profiling of pathogen and AMR genes within 6 h [52] . They achieved a sensitivity of 96.6% and a limit of detection similar to that of culture-based methods. Determining the taxonomic entities present in a sample from a metagenomic sequencing dataset is a key step in metagenomics studies. The assignment of taxonomic labels to the sequencing reads and the subsequent inference of the composition of a microbial community are increasingly popular research areas due to the growing use of high-throughput technologies demanding more accurate and efficient tools for metagenomic analyses (Table 2) . Generally, long reads enable better taxonomic and functional analysis than short reads due to the higher information content enclosed in the sequence. Yet, most of the widely adopted metagenomic classification tools or pipelines often rely on algorithms built on short reads which, by default, do not scale well with long-read datasets -ranging from 13 kb up to 2 Mb -and/or do not account for the higher error profile of nanopore reads. The inclusion of long-read datasets in benchmarking studies and bioinformatic software updates have provided some guidance for metagenomic tool suitability and performance with nanopore reads [53] . Furthermore, long-read specific tools are being continuously developed [54] , including error correction methods [55] and hybrid approaches to overcome read error-related issues [56, 57] . Traditional read classification methods are based on the detection of similarities between sequencing reads and genomes from known organisms through an initial alignment against databases containing taxonomic information. Of these, the most popular are classification tools based on the classic BLAST algorithm [58] , which remains the gold standard for the taxonomic assignment task. Apart from the classic nucleotide BLAST, more recently developed methods like MEGABLAST [59] provide faster alignments. Alternative methods built upon BLAST or other alignment tools have also been proposed to improve classification results combining sequence alignment of input reads with machine learning techniques for taxonomic resolution at different levels. For example, MEGAN-LR [60] expands the functionality of the interactive metagenomics pipeline featuring long-read approaches for taxonomic and functional analysis. It adopts alignment-based comparisons using LAST aligner [61] to compute frameshiftaware DNA-to-protein alignments and applies a custom lowest common ancestor (LCA) algorithm to resolve taxonomy and deliver classification results. Alignment-based techniques also output useful information for results interpretation such as genomic locations and qualities of alignments. This feature usually comes at the cost of an increase in required computing resources and time when analyzing long-read datasets. Because of this, in the last few years, alignment-free classification methods have become popular for the analysis of short-and long-read datasets. These methods mainly rely on a k-mer based classification against precomputed indexes and guarantee the efficient search and storage of sequence databases. As a result, most of these tools are capable of classifying millions of reads per minute with a relatively small memory footprint and enable the analysis of extensive long-read datasets. A major limitation of k-mer based methods is that they are sensitive to low error rates and may lead to misclassifications when used with error-prone long reads, especially when classifying similar organisms at the species level or organisms that have high sequence identity. However, the continuous chemistry updates and the release of novel basecalling algorithms, such as Bonito (https://github.com/nanoporetech/bonito) [62, 63] , have improved the raw read accuracy, leading to an overall increase in the performance of downstream analysis tools. Despite that, further inspection and post hoc analysis of k-mer based classification outputs have been suggested in order to limit possible misclassifications [58] [59] [60] . An example of an alignment-free classification tool is Kraken [64] , which uses exact k-mer matches for each read against a k-mer-to-LCA records database. These records are generated from sequence databases and indexed in time-efficient data structures that enable faster look-up searches. However, this process is memory-intensive, an issue that was addressed by the development of Kraken2 [35] , characterized by an enhanced database efficiency and improved k-mer-based read analysis. Using the same Kraken k-mer based classification technique, KrakenUniq [65] improves precision and recall by assessing the coverage of unique k-mers of each taxon that is present in the dataset. An additional development leveraging the Kraken classification output is Bracken [66] , a statistical method to compute the abundance estimation of a sample at any given taxonomic level from Kraken/ Kraken2 classifications for each read. Another software, Centrifuge , builds a data structure based on the Ferragina-Manzini index, a technique based on widely used read aligners, i.e. the Burrows-Wheeler Aligner (BWA) algorithm [67] and Bowtie [68] . This data structure provides efficient storage of database sequences and the classification is also performed by k-mer matching against the pre-built index. Metamaps [69] is a long-read specific approach featuring taxonomic assignment and sample composition estimations at strain-level along with an output that includes per read positional and quality information. In addition, a number of pipelines for analysis of metagenomic data have been developed, which integrate previously mentioned tools into easy-to-deploy workflows and generate comprehensive outputs for result interpretation ( Taxonomic classification methods use pre-computed or indexed reference databases (Table 3) . While some tools are designed to work with specific databases, most tools allow a variety of sequence collections to be indexed. Thus, the database choice is important for metagenomic workflows. Popular reference databases include the NCBI RefSeq collection of complete genomes [75] , encompassing both prokaryotic and eukaryotic genomes, and the nt BLAST database built from more than 50 million highquality nucleotide sequences. The GenBank database [76] includes a wider collection of complete genomes albeit with lower quality standards than RefSeq collections. Other databases that are better suited for metataxonomics include GreenGenes [77] , SILVA [78] , RDP [79] for 16S rRNA gene sequencing, and the NCBI RefSeqTargeted Loci Project database containing 16S/23S and 18S/28S rRNA genes from the GenBank database (https://www.ncbi.nlm.nih.gov/ refseq/targetedloci/). These databases contain partial and fulllength 16S rRNA gene sequences providing a more lightweight and comprehensive sequence collection for metataxonomic analysis. Other purpose databases are Prokka [80] , for the annotation of assembled genomes, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [81] , which includes annotated genes and genomes for functional annotation profiling. Regarding AMR-specific databases, CARD [82] and MEGARes [83] contain well-curated gene sequences and annotations for integration in the AMR analysis and detection workflows. The size of microbial genome databases grows exponentially every year and can reach 100s of GBs for some of the large collections, e.g., RefSeq. Despite that, their size increase does not guarantee the successful classification of every generated read of the experiment [84] . False negatives are common due to undiscovered and yet to be sequenced microorganisms. Modifications or updates of databases include changes in the structure of the taxonomy tree and the inclusion of new sequences or resequenced strains. The increase in the number of total sequences added to a database (through updates) is not representative of a gain on species richness. Some of these additions tend to create redundancy of certain genera and species in the database. Hence, the use of databases created at different times and, therefore, containing different sequence content and taxonomy can affect the results and profoundly confound the software benchmarkings [85] . These issues emphasize the importance of performing continuous comparisons and benchmarkings of widely adopted tools using varied testing datasets and databases, necessary for result interpretation and evaluation of computational requirements and performance. In addition to the improvements in the taxonomic analysis when relying on long reads, it has been shown that nanopore sequencing data can produce contiguous and highly accurate assemblies, including full-length plasmids and viruses, without any preprocessing step, such as initial read binning [88] . These assemblies are in the range of 95% of completeness when applied to whole-genome sequencing combined with posterior assembly polishing using error correction tools [55] . Long-read metagenomic assembly has the potential not only to improve contiguity over short-read assemblies but also to enable strain resolution, to sequence novel plasmids and viruses, and to enhance the power of identifying horizontal gene transfer. Despite the advantages over short-read sequencing data, complete contiguous assemblies are still constrained by the relatively high error rate of nanopore reads, and the quality of the metagenome assembly is related to the coverage of the different species present in the sample, which in turn depends on the experiment throughput. In practice, the readlength advantage of nanopore enables nearly complete assemblies even for low abundance strains, provided that they are covered at a minimal level [89] . Metagenomic assembly requires the development of specificpurpose algorithms to overcome the limitations of classic assemblers which assume that the depth of coverage is approximately uniform across the genome. Some long-read assemblers have been used for metagenomic assembly even if not specifically designed for the task ( Table 2) . For example, Canu [90] is one of the first and most popular assemblers for long-read data and generates contigs using an adaptive k-mer weighting strategy to produce an assembly with high coverage long-read data. wtdbg2 [91] adopts a fast all-versus-all read alignment and a layout algorithm based on fuzzy-Bruijn graphs for sequence assembly. A recent development is Raven [92] , a specific long-read assembler that features a faster read overlap step and assembly graph building. While these approaches have been successfully used and benchmarked for metagenomics assembly [20, 93, 94] , novel assemblers that were built specifically for metagenomics datasets have recently been developed. An example is OPERA-MS [95] , a hybrid assembler that leverages the strengths of both short-read and long-read sequencing approaches using first a short-read assembler to create contigs and then the long-read information to create an assembly graph for all genomes that are distinguished in a coverage-based clustering. Another software is MetaFlye [96] , which selects high-frequency kmers from the dataset to detect read overlaps and creates errorprone contigs used for developing the final assembly graph. Nanopore sequencing has been used in a number of clinically relevant metagenomics applications due to the possibility of realtime sequencing (Table 4 ). One example is the rapid species and AMR profiling to guide proper antibiotic treatment, which is of particular importance for the current worldwide AMR threat. A recent study used nanopore sequencing for the characterization of pathogenic bacteria and AMR genes in gut-associated microbial commu-nities in preterm infants [97] . To do this, they used the NanoOK RT software, where reads are aligned to bacteria and AMR databases on the fly while they are generated in the run, resulting in a sequencing turnaround time of 1 h. Furthermore, they successfully linked AMR genes to the harboring pathogens, suggesting that nanopore sequencing coupled with bioinformatic analysis can tailor antibiotic therapies. Similarly, Břinda et al. have recently developed an innovative method, called genomic neighbor typing, which accelerates pathogen detection and AMR typing [98] . The method is based on the assumption that resistance elements are genetically linked to the rest of the genome and that it is possible to define the antimicrobial susceptibility by only inferring the bacterial strains that are present in the sample. It relies on a two-step algorithm where the sequence is first compared to a reference database, and then the most probable phenotype (drug resistance or susceptibility) of the sample is determined from the phenotype of its nearest genomic neighbor. They were able to identify the correct lineage, strain, and antibiotic susceptibility of pneumococcal and gonococcal isolates in < 10 min. They also compared it to the AMR gene-based approach, which required 25 min for single copies to be detected. They further validated the method by successfully confirming resistance in pneumococcus in six sputum samples from patients suffering lower respiratory infections. The main limitation of this method is that strain and AMR detection are based on the information provided by the database. Therefore, it is crucial that the constructed database includes genomic sequences and resistance metadata of the strains encountered in the clinical samples analysed. As a consequence, further applications of this method may include pathogen diagnostics and surveillance, provided that the target microbe and AMR are known. Given that long reads facilitate accessing repetitive elements and structural variants, nanopore sequencing has improved our capacity for de novo assembly of genomes, metagenomes, and plasmids, which in turn allow gathering information on the localization of resistance and virulence factors, such as AMR genes and PAIs. An example of this is a recent study that collected the most comprehensive characterization of opportunistic pathogens and their resistomes colonizing tertiary hospital environments [99] . The workflow consisted of culturing, antibiotic selection, metagenomic sequencing, and OPERA-MS-based genome assembly of environmental samples collected from 179 sites associated with 45 beds. They obtained genomes for 69 species, 16% of them from novel species. Furthermore, they reconstructed plasmids and phages, of which more than 90% were uncharacterized. By doing this, they were able to identify novel associations between AMR genes, characterize chromosomal cassettes and AMR gene combinations in plasmids, and detect the persistence of multidrugresistant organisms in hospital environments over years. This study highlighted the importance of monitoring environmental pathogens in clinical settings, which can be responsible for occasional outbreaks in hospitals. More recently, another research group has introduced a novel long-read assembly-based approach for metagenome-assembled genomes (MAGs), called Lathe [43] . By coupling Lathe to a newly developed experimental protocol for HMW DNA extraction, they assembled seven genomes (out of 12) from a mock sample into single contigs, obtaining circular genomes, while three more genomes were assembled into four or fewer contigs. They validated their protocol on 13 human stool samples and demonstrated that Lathe generated better assembly contiguity than those from short reads and the read cloud assembler (hybrid). Notably, they were able to resolve the circular genome of Prevotella copri, known to be characterized by a high degree of repetitive sequences. By resolving complete circularized genomes, it is possible to optimally study microbial phenomena, such as horizontal gene transfers, and to investigate how inter-strain structural variants may be linked to a specific phenotype of a microbial community [100] . Because of their reduced length (typically within 3-300 kb), viral genomes can be sequenced as a single molecule using nanopore sequencing. Beaulaurier et al. discovered more than 1,800 phage genomes in seawater samples [101] . The analytical method consisted at first of filtering those sequences containing the direct terminal regions, which are characteristic of the virus genome termini in dsDNA-tailed phages. Then, the analysis consisted of a step of dimensionality reduction and clustering, followed by polishing to create high-quality draft phage genomes. By adopting this method, they were able to identify viral microheterogeneity, otherwise very difficult to detect using short-read sequencing. Furthermore, this approach allowed inferring the phage packaging strategy and identifying concatemers of sequences similar to the phage-inducible chromosomal islands, revealing the utility of this approach to identify repeat sequences derived from phage-induced mobile elements. The study of the microbiome can be approached through metatranscriptomics, i.e. the study of the totality of transcripts in a sample. Nanopore sequencing enables obtaining full-length transcripts in a single read, facilitating transcriptome analysis by avoiding the challenging steps necessary for short-read transcriptomics. Furthermore, ONT technologies can directly sequence the RNA molecules eliminating the biases introduced by the reverse transcription or the amplification step, given that all transcripts do not amplify with the same efficiency [13, 102] . In addition, the processes of retrotranscription and amplification erase the epitranscriptomic information, which is known to have a role in modulating transcript activity and stability. Viral RNA genomes can also be sequenced as native RNA molecules or as cDNA after retrotranscription using nanopore sequencing. This is of particular importance for many emerging human viral diseases, such as Ebola, severe acute respiratory syndrome (SARS), and the coronavirus disease 2019 (COVID-19), all of them caused by RNA viruses. ONT offers three main sequencing library preparation kits for the analysis of transcriptome and viral RNA genomes ( Table 1) . Two of them (the Direct cDNA Sequencing Kit and the cDNA PCR Sequencing Kit) are based on a retrotranscription step, followed by either digestion of the RNA strand and ligation of sequencing primers, or by a PCR step with rapid attachment primers when the initial target RNA does not reach the minimum required amount of 100 ng. A third library preparation kit, the Direct RNA Sequencing Kit, is based on a DNA primer annealing and ligation to the RNA strand, followed by an optional retrotranscription step (necessary only to stabilize the RNA molecule) and by the attachment of sequencing adapter at the RNA 3 0 end. This library preparation protocol is faster (< 2 h) because it does not include cDNA synthesis. However, it requires a higher amount of initial RNA (~500 ng) and has lower throughput compared to the cDNA kits. However, ONT is continuously modifying the chemistry of these kits in order to improve current accuracy and throughput. Although ONT successfully enabled gaining insights into eukaryotic messenger RNA [103] and viral RNA with a poly-A tail [104, 105] , the study of the prokaryotic transcriptome has been hindered by the lack of a poly-A tail required for the attachment of primers during library construction. A way to overcome this issue would be by adding a step of polyadenylation of prokaryotic transcriptomes in the experimental protocol to make them recognizable and modifiable by the ONT Direct RNA Sequencing Kit [106] . Another possibility is to design custom adapters to ligate the 3 0 end of the transcripts, the tRNAs, or the rRNAs of interest. For example, Smith et al. designed adapters containing a 20nucleotide overhanging sequence targeting the conserved anti-Shine Dalgarno region, present in prokaryotic 16S rRNA, to study how canonical and non-canonical base modifications affect antimicrobial susceptibility in Escherichia coli strains [107] . A similar approach was employed by Keller et al. to sequence, for the first time, the complete RNA genome of the influenza A virus, by designing adapters targeting the highly conserved genome termini of the virus [108] . Whereas some well-known pipelines for short-read metagenomics, such as MEGAN [30] and MG-RAST [109] , can be used with RNA reads for taxonomic assignment, the availability of specific tools to analyze long-read RNA profiles is limited. Metataxonomic workflows, as described elsewhere [110] , can be performed alternatively by extracting rRNA sequences, such as the small subunits (16S/18S) and large subunits (23S/28S), using specialized software, e.g., METAXA2 [111] . Functional analysis is performed with BLASTx or Magic-BLAST [112] to align the RNA sequences to a protein database in order to assign either Gene Ontology terms (GO) with Blas-t2GO [113] or metabolic pathway annotations according to KEGG [114] . Recent examples of specific tools for long-read transcriptomics are Poreplex (https://github.com/hyeshik/poreplex), a signal-level processor for ONT direct RNA sequencing data that features real-time basecalling, quality filtering, 3 0 adapter trimming, and alignment to reference transcriptomes. Bambu [115] is an R software package for multi-sample transcript discovery and quantification from long-read RNA data. Regarding alignment-free methods, a recent development is isONclust [116] , a tool for de novo transcript reconstruction with a cluster-based approach that accounts for large dataset scaling. The study of the transcriptome for metagenomics can be adopted to address several issues which cannot be tackled by DNA sequencing. In fact, RNA sequencing provides additional information, such as the functionality of AMR genes, allowing to identify a situation where the resistance gene is present but not transcribed, and therefore does not generate a resistant phenotype [106] . Metatranscriptomics is also very useful for the identification of viable pathogens since DNA-based approaches are unable to differentiate between viable and unviable bacterial cells [117] . This approach is particularly important in the detection of food pathogens, where food processing and storage often kill bacteria cells without removing their genomic DNA. Direct RNA sequencing has recently been compared to multiplex real-time PCR amplicon sequencing for this purpose. Results suggest it to be especially applicable to complex microbiomes because it does not require assay customization for specific biohazards when a complete database is used during the bioinformatic analysis step [117] . Other applications of metatranscriptomics using nanopore sequencing also include pathogen detection from clinical samples, which is particularly useful for diseases caused by RNA viruses (Table 4) . For example, nanopore RNA sequencing has been used for differential diagnosis of dengue and chikungunya viruses, two singlestranded positive RNA viruses circulating in the same geographical areas and causing diseases with similar symptomatology [118] . Despite that, transcriptomics is still a very immature application of nanopore sequencing for microbial studies. In recent years, emerging RNA viruses have become a threat to global health, and viral genome sequencing has turned into an essential tool for outbreak identification and monitoring of transmission patterns. Nanopore technology has been demonstrated to be an exceptionally valuable tool for this purpose because it can produce data in real-time, directly in-field and under extreme conditions thanks to inexpensive portable devices. For this reason, nanopore technology has been adopted for genomic surveillance during the Ebola outbreak in West Africa [119] , the Zika outbreak in Brazil and the Americas [120] , and the ongoing COVID-19 pandemic. Since the first identification of a novel coronavirus in December 2019 [121] , thousands of SARS-CoV-2 genomes have been sequenced using the ARTIC [122] or alternative [123] protocols developed for fast sequencing of the virus using nanopore technology, allowing to gather information on virus biology, transmission, and viral dynamics. For example, an important study performed by Fauver et al. coupled genomic data with domestic and international travel patterns in the USA, tracking down the SARS-CoV-2 transmission dynamics in early March 2020 [124] . In this study, nine viral genomes from early cases in Connecticut were sequenced within 24 h and used to build a phylogenetic tree. When compared to other 168 publicly available genomes at that time, seven out of nine genomes clustered into one clade containing sequences from other USA samples, suggesting domestic transmission of SARS-CoV-2 in the USA early in the first wave of the pandemic. This information was further confirmed by estimating the SARS-CoV-2 travel importation risk into Connecticut using airline travel data and epidemiological dynamics in regions where travel routes came from. Nanopore sequencing of SARS-CoV-2 has also been adopted in a prospective genomic surveillance study aiming to identify healthcareassociated infections in a hospital in the UK [125] . Around 1,000 genomes were sequenced within five months and results were compared on the basis of the ward location data of patients or healthcare workers in order to unravel transmission patterns within the hospital. This information was transmitted almost in real-time to the hospital management team and allowed to identify risk factors for transmission in clinical settings. Importantly, this study supports the adoption of combined epidemiological and genomic data for the implementation of infection control measures and highlights the importance of genomic epidemiology to guide decision-making on a local, national and international level. Apart from the above-described protocols, which require a retrotranscription step, direct RNA sequencing has also been used for SARS-CoV-2 sequencing [126] [127] [128] . This method allowed to sequence regions spanning almost the entire viral genome (~30 kb), although the coverage was found to be extremely variable along the genome, ranging from 34Â to >160,000Â, and biased towards the poly-A 3 0 end [128] . The reason for this is the abundance of subgenomic mRNAs carrying these regions and the directional nanopore sequencing from the poly-A 3 0 end. Nonetheless, this study allowed gaining insights into the transcriptome and epitranscriptome of the virus, with the identification of eight major transcripts and 42 positions with 5 0 methyl-cytosine (5mC) modifications. Furthermore, by direct RNA sequencing, Taiaroa et al. were able to estimate the evolutionary rate of the virus, which is important for epidemiological studies [128] . Epigenetic modifications of DNA in bacteria and DNA/RNA in viruses are responsible for several biological functions, such as the regulation of DNA/RNA replication and repair, control of gene expression, and protection from external pathogens [129] . So far, methylation is the only nucleotide modification known in bacterial DNA, with three forms of methylation identified: 5mC, N4methylcytosine (4mC), and N6-methyladenine (6 mA), the latter being the most prevalent form [130] . Each of these types of epimodifications occurs in a highly motif-driven manner, where every occurrence of the motif is methylated. Nanopore sequencing allows direct detection of the native modified bases on the nucleic acid during its passage through the nanopore. In fact, the characteristic ionic current observed when a certain sequence passes through the pore is altered by the presence of a methylated base, generating a distinctive current pattern that can be distinguished from the non-methylated DNA/RNA. While ONT amplificationfree libraries can be easily generated via standard ONT kits, the bottleneck for the study of nucleotide modifications in nanopore sequencing is still the basecalling process, where the presence of multiple new current signals generated by one or multiple methylated bases in the k-mer passing the pore causes a considerable computational challenge [12, 131] . Multiple research groups have tried to face this by developing tools to detect methylated bases. For example, Stoiber et al. presented a method based on the statistical comparison between ionic current signals from native and methylated sequences [132] , which has then evolved into the current ONT Tombo platform (https://github.com/nanoporetech/tombo). Other methods use pre-trained classification models to capture epigenomic modifications such as Nanopolish [131] and SignalAlign [12] , which use Hidden-Markov models, or the recently developed mCaller [133] , Deepsignal [134] , and Deepmod [135] which adopt neural network classifiers. However, these tools are characterised by detection accuracies that vary based on the methylation type and the target motif, and the capability to detect de novo methylated motifs is limited by the training data [136] . A recent study [137] tried to address this issue by generating a large training dataset for de novo methylation typing and mapping of all three forms of DNA methylation and applied it to individual bacteria and mouse gut microbiome samples. In this work, Tourancheau et al. also developed a novel approach for methylation binning of metagenome contigs and demonstrated how methylation patterns may assist in the process of metagenome assembly. Although this method enabled de novo methylation typing and fine mapping, accuracy is still highly dependent on the type and position of the methylated base which remains an issue to be addressed in the future. Nanopore technology has improved many aspects of microbial analysis and has the potential to be adopted routinely in clinical settings in the near future. In fact, many proof-of-concept studies have demonstrated that nanopore sequencing can be adopted for infectious disease diagnostics and for monitoring the human microbiome, which can be a useful tool in clinical medicine. For example, dysbiosis of the lung microbiome has been shown to have a prognostic value for mortality in patients with non-pulmonary sepsis in intensive care units [138] , and nanopore sequencing was proposed to be used as a prognostic tool for real-time monitoring of this dysbiosis. Similarly, nanopore sequencing could be used to monitor changes in the gut microbiome over time, before or after antibiotic treatments [97] , or to assess species engraftment after faecal microbiota transplantation [139] . Although there are many possible applications for nanopore sequencing in the field of metagenomics, there are still numerous challenges that need to be addressed. For example, standardized protocols for microbial characterization are needed in order to use this technology in clinical settings. Furthermore, novel and efficient HMW DNA extraction protocols from microbial samples are required to pro-duce high-quality long reads, and library preparation protocols need to be simplified, especially for in-field and educational settings. For this purpose, VolTRAX has been released by ONT as a system for automated library preparation to provide high reproducibility and portability to the library preparation step. However, the cartridge used by the system can hold up a maximum of ten barcoded samples, which is far too low for experiments requiring multiplexing of 96 samples. Another aspect that needs to be considered is the read error rate. This issue is continuously addressed by ONT through the constant improvement of the flow cell and the release of more rapid and accurate basecalling methods. However, both read and consensus accuracy are limited by the organisms chosen for model training, and they are drastically reduced when the basecaller is used for less frequently sequenced microbial species [140] . This issue can be addressed by developing more custom-trained basecallers, built-on taxonspecific training data so that the users can choose which basecaller most closely matches the organisms present in their samples [140, 141] . Sequencing data throughput by single flow cell is also continuously improving, enabling the detection and sequencing of DNA/RNA from microbes even when present in very low abundance. Furthermore, real-time selective sequencing [142] , or read until, is also becoming popular among nanopore users and consists of extruding specific molecules from the pores, such as host DNA or other non-interesting molecules. It has the potential to increase the efficiency of the run by reducing the time taken to complete an experiment, to enrich the data with the less represented genomes in the sample and, at the same time, to simplify library preparation protocols by eliminating host DNA depletion or target enrichment steps [143, 144] . Regarding the analysis of metagenomic data, improvements are expected in taxonomic analysis and sequence comparison software in order to achieve better resolution of closely related strains and higher classification accuracy, along with the development of efficient indexing techniques for metagenomics databases. Metagenomic assembly software and hybrid techniques using both short and long reads have the potential to enhance the analysis of complex samples improving the detection of unknown organisms and enabling the assembly of mobile elements and resistance genes, which are crucial for the characterization of complex microbial environments such as the human microbiome [43, 145] . Advances in nanopore technologies, such as direct RNA sequencing [13] and the detection of epimodifications [12] , have highlighted the need for novel bioinformatic tools enabling accurate characterization or discovery of transcriptomes and the determination of the type and position of modified bases on a sequence and their functional impact. Furthermore, benchmarking studies of metagenomics tools for the analysis of long reads are lacking, together with long-read metagenomic datasets representing complex microbial communities to be used during software development and tool assessments. Future advancements in metagenomic analysis tools and workflows for long reads will need to follow the quick pace in modifications and updates of the sequencing technology and also account for software efficiency and scalability in order to enable the analysis of sequence data from high-throughput devices such as the GridION and PromethION. This work was supported by the Ministerio de Ciencia e Innovación (RTC-2017-6471-1; AEI/FEDER, UE) and the Instituto de Salud Carlos III (PI14/00844, PI17/00610, FI18/00230), which were co-financed by the European Regional Development Funds 'A way of making Europe' from the European Union; Fundación Canaria Instituto de Investigación Sanitaria de Canarias (PIFUN48/18); Cabildo Insular de Tenerife (CGIEU0000219140 and ''Apuestas científicas del ITER para colaborar en la lucha contra la COVID-19"); and by agreement OA17/008 with the Instituto Tecnológico y de Energías Renovables (ITER) to strengthen scientific and technological education, training, research, development, and innovation in genomics, personalized medicine, and biotechnology. The sponsors had no involvement in the review conceptualization, the manuscript writing and the decision to submit the article for publication. CRediT authorship contribution statement The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Rapid determination of 16S ribosomal RNA sequences for phylogenetic analyses The Analysis of Natural Microbial Populations by Ribosomal RNA Sequences A human gut microbial gene catalogue established by metagenomic sequencing Unusual biology across a group comprising more than 15% of domain Bacteria The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community Oxford Nanopore MinION Sequencing and Genome Assembly BulkVis: a graphical viewer for Oxford nanopore bulk FAST5 files A world of opportunities with nanopore sequencing A complete bacterial genome assembled de novo using only nanopore sequencing data Accurate characterization of expanded tandem repeat length and sequence through whole genome long-read sequencing on PromethION Analysis of short tandem repeat expansions and their methylation state with nanopore sequencing Mapping DNA methylation with high-throughput nanopore sequencing Highly parallel direct RN A sequencing on an array of nanopores From squiggle to basepair: Computational approaches for improving nanopore sequencing read accuracy Oxford Nanopore Technologies. Nanopore sequencing accuracy The vocabulary of microbiome research: a proposal Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases Then and now: use of 16S rDNA gene sequencing for bacterial identification and discovery of novel bacteria in clinical microbiology laboratories Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes Current challenges and best-practice protocols for microbiome analysis Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis Species-level resolution of 16S rRNA gene amplicons sequenced through the MinION TM portable nanopore sequencer Analysis of the mouse gut microbiome using full-length 16S rRNA amplicon sequencing A preliminary study on the potential of Nanopore MinION and Illumina MiSeq 16S rRNA gene sequencing to characterize building-dust microbiomes NanoCLUST: a species-level analysis of 16S rRNA nanopore sequencing data btaa900 PCR-based quantification of taxa-specific abundances in microbial communities: Quantifying and avoiding common pitfalls Rapid bacterial identification by direct PCR amplification of 16S rRNA genes using the MinION TM nanopore sequencer Full-length 16S rRNA gene amplicon analysis of human gut microbiota using MinION TM nanopore sequencing confers species-level resolution 16S rRNA gene copy number normalization does not provide more reliable conclusions in metataxonomic surveys MEGAN analysis of metagenomic data Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2 Computational methods for 16S etabarcoding studies using Nanopore sequencing data Sequence analysis Minimap2 : pairwise alignment for nucleotide sequences Centrifuge: Rapid and sensitive classification of metagenomic sequences Improved metagenomic analysis with Kraken 2 A portable system for rapid bacterial composition analysis using a nanoporebased sequencer and laptop computer Real-time diagnostic analysis of MinION TM -based metagenomic sequencing in clinical microbiology evaluation: a case report Rapid sequencing-based diagnosis of infectious bacterial species from meningitis patients in Zambia A sample-to-report solution for taxonomic identification of cultured bacteria in the clinical setting based on nanopore sequencing Shotgun metagenomics, from sampling to analysis Metagenomics: genomic analysis of microbial communities New dimensions of the virus world discovered through metagenomics Complete, closed bacterial genomes from microbiomes using nanopore sequencing MinION nanopore sequencing identifies the position and structure of a bacterial antibiotic resistance island Completing bacterial genome assemblies with multiplex MinION sequencing Comparison of bacterial genome assembly software for MinION data and their applicability to medical microbiology Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches Complete genome sequences of bordetella pertussis clinical isolate FR5810 and reference strain tohama from combined oxford nanopore and illumina sequencing Improved high-molecular-weight DNA extraction, nanopore sequencing and metagenomic assembly from the human gut microbiome Metagenomic sequencing with Oxford Nanopore 2020 New opportunities for managing acute and chronic lung infections Nanopore metagenomics enables rapid clinical diagnosis of bacterial lower respiratory infection Benchmarking the MinION: evaluating long reads for microbial profiling Bioinformatics of nanopore sequencing A comprehensive evaluation of long read error correction methods A comparative evaluation of hybrid error correction methods for error-prone long reads Benchmarking long-read assemblers for genomic analyses of bacterial pathogens using oxford nanopore sequencing Basic local alignment search tool High speed BLASTN: an accelerated MegaBLAST search tool MEGAN-LR: new algorithms allow accurate binning and easy interactive exploration of metagenomic long reads and contigs Adaptive seeds tame genomic sequence comparison Halcyon: an accurate basecaller exploiting an encoder-decoder model with monotonic attention Pair consensus decoding improves accuracy of neural network basecallers for nanopore sequencing Ultrafast metagenomic sequence classification using exact alignments KrakenUniq: confident and fast metagenomics classification using unique k-mer counts Bracken: estimating species abundance in metagenomics data Fast and accurate short read alignment with Burrows-Wheeler transform Strain-level metagenomic assignment and compositional estimation for long reads with MetaMaps Metagenomics workflow for hybrid assembly, differential coverage binning, transcriptomics and pathway analysis (MUFFIN) NanoSPC: a scalable, portable, cloud compatible viral nanopore metagenomic data processing pipeline Mash: fast genome and metagenome distance estimation using MinHash Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences MetaSPAdes: A new versatile metagenomic assembler RefSeq : an update on prokaryotic genome annotation and curation Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB The SILVA ribosomal RNA gene database project: Improved data processing and webbased tools Ribosomal Database Project: Data and tools for high throughput rRNA analysis Prokka: rapid prokaryotic genome annotation KEGG: integrating viruses and cellular organisms Antibiotic resistome surveillance with the comprehensive antibiotic resistance database MEGARes: An antimicrobial resistance database for high throughput sequencing RefSeq database growth influences the accuracy of k-mer-based lowest common ancestor species identification Duplicates, redundancies and inconsistencies in the primary nucleotide databases: A descriptive study An integrated catalog of reference genes in the human gut microbiome DeepARG: a deep learning approach for predicting antibiotic resistance genes from metagenomic data Ultra-deep, long-read nanopore sequencing of mock microbial community standards Implications of error-prone long-read whole-genome shotgun sequencing on characterizing reference microbiomes Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation Fast and accurate long-read assembly with wtdbg2 Raven: a de novo genome assembler for long reads Compendium of 4,941 rumen metagenome-assembled genomes for rumen microbiome biology and enzyme discovery Assembly methods for nanopore-based metagenomic sequencing: a comparative study Hybrid metagenomic assembly enables high-resolution analysis of resistance L. Ciuffreda, Héctor Rodríguez-Pérez and C metaFlye: scalable long-read metagenome assembly using repeat graphs Rapid MinION profiling of preterm microbiota and antimicrobial-resistant pathogens Rapid inference of antibiotic resistance and susceptibility by genomic neighbour typing Cartography of opportunistic pathogens and antibiotic resistance genes in a tertiary hospital environment Structural variation in the gut microbiome associates with host health Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities Accurate detection of m6A RNA modifications in native RNA sequences Nanopore longread RNAseq reveals widespread transcriptional variation among the surface receptors of individual B cells Multiple long-read sequencing survey of herpes simplex virus dynamic transcriptome Novel splicing and open reading frames revealed by long-read direct RNA sequencing of adenovirus transcripts Evaluating the genome and resistome of extensively drug-resistant Klebsiella pneumoniae using native DNA and RNA Nanopore sequencing Reading canonical and modified nucleobases in 16S ribosomal RNA using nanopore native RNA sequencing Direct RNA sequencing of the coding complete influenza a virus genome MG-RAST, a metagenomics service for analysis of microbial community structure and function Evaluating the potential of direct RNA nanopore sequencing: Metatranscriptomics highlights possible seasonal differences in a marine pelagic crustacean zooplankton community METAXA2: Improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data Magic-BLAST, an accurate RNA-seq aligner for long and short reads Blast2GO: A comprehensive suite for functional analysis in plant genomics New perspectives on genomes, pathways, diseases and drugs De novo clustering of long-read transcriptome data using a greedy, quality value-based algorithm Direct metatranscriptome RNA-seq and multiplex RT-PCR amplicon sequencing on nanopore MinIONpromising strategies for multiplex identification of viable pathogens in food Assessment of metagenomic Nanopore and Illumina sequencing for recovering whole genome sequences of chikungunya and dengue viruses directly from clinical samples Realtime, portable genome sequencing for Ebola surveillance Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples A novel coronavirus from patients with pneumonia in China nCoV-2019 sequencing protocol v3 (LoCost) Rapid and inexpensive wholegenome sequencing of SARS-CoV-2 using 1200 bp tiled amplicons and oxford nanopore rapid barcoding bpaa014 Coast-to-Coast Spread of SARS-CoV-2 during the Early Epidemic in the United States 990-996.e5 Rapid implementation of SARS-CoV-2 sequencing to investigate cases of health-care associated COVID-19: a prospective genomic surveillance study Direct RNA nanopore sequencing of full-length coronavirus genomes provides novel insights into structural variants and enables modification analysis The architecture of SARS-CoV-2 transcriptome Direct RNA sequencing and early evolution of SARS-CoV-2 Epigenetic Gene Regulation in the Bacterial World Deciphering bacterial epigenomes using modern sequencing technologies Detecting DNA cytosine methylation using nanopore sequencing De novo identification of DNA modifications enabled by genome-guided Nanopore signal processing Single-molecule sequencing detection of N6-methyladenine in microbial reference materials DeepSignal: detecting DNA methylation state from Nanopore sequencing reads using deep-learning Detection of DNA base modifications by deep recurrent neural network on Oxford Nanopore sequencing data Opportunities and challenges in long-read sequencing data analysis Discovering and exploiting multiple types of DNA methylation from individual bacteria and microbiome using nanopore sequencing Could lung bacterial dysbiosis predict ICU mortality in patients with extra-pulmonary sepsis? A proof-of-concept study Strand-wise and baitassisted assembly of nearly-full rrn operons applied to assess species engraftment after faecal microbiota transplantation Performance of neural network basecalling tools for Oxford Nanopore sequencing High quality genome assemblies of Mycoplasma bovis using a taxon-specific Bonito basecaller for MinION and Flongle long-read nanopore sequencing Real-time selective sequencing using nanopore technology Targeted nanopore sequencing by real-time mapping of raw electrical signal with UNCALLED Real-Time, Direct Classification of Nanopore Signals with SquiggleNet Rapid resistome mapping using nanopore sequencing