Antimicrobial resistance prediction and phylogenetic analysis of Neisseria gonorrhoeae isolates using the Oxford Nanopore MinION sequencer

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is common, compromising gonorrhoea treatment internationally. Rapid characterisation of AMR strains could ensure appropriate and personalised treatment, and support identification and investigation of gonorrhoea outbreaks in nearly real-time. Whole-genome sequencing is ideal for investigation of emergence and dissemination of AMR determinants, predicting AMR, in the gonococcal population and spread of AMR strains in the human population. The novel, rapid and revolutionary long-read sequencer MinION is a small hand-held device that generates bacterial genomes within one day. However, accuracy of MinION reads has been suboptimal for many objectives and the MinION has not been evaluated for gonococci. In this first MinION study for gonococci, we show that MinION-derived sequences analysed with existing open-access, web-based sequence analysis tools are not sufficiently accurate to identify key gonococcal AMR determinants. Nevertheless, using an in house-developed CLC Genomics Workbench including de novo assembly and optimised BLAST algorithms, we show that 2D ONT-derived sequences can be used for accurate prediction of decreased susceptibility or resistance to recommended antimicrobials in gonococcal isolates. We also show that the 2D ONT-derived sequences are useful for rapid phylogenomic-based molecular epidemiological investigations, and, in hybrid assemblies with Illumina sequences, for producing contiguous assemblies and finished reference genomes.

amplification steps 23 . Individual DNA molecules can also be sequenced by monitoring their transfer through various types of pores 24,25 , which can potentially result in very long and unbiased sequence reads, because no amplification or chemical reactions are used for the sequencing 26 . Oxford Nanopore Technologies (ONT) (Oxford, UK) have introduced this approach with their single-molecule nanopore genome sequencing device MinION, a TGS platform with unique technology that was commercialized in mid-2015 [26][27][28][29][30] . When both strands are sequenced, a consensus sequence of the molecule can be produced; these consensus reads are termed two-directional reads (1D 2 or 2D ONT reads) and generally have higher accuracy than reads from only a single pass of the molecule (1D ONT reads). The MinION provides several new advantages: small size (10 × 3 × 2 cm), portability, speed, low cost, and direct connection to a laptop through a USB 3.0 interface; the library construction involves simplified methods; no amplification step is required, and data acquisition and analysis occur in real time. The MinION is promising for microbiological applications including describing the microbiome 31 , rapid diagnostics 32 , transmission and surveillance 33,34 , de novo assemblies [35][36][37] , and microbial AMR profiling 38,39 . However, the high error rate for the MinION sequencer 27,28 has limited its ability to compete with existing sequencing technologies and the MinION has not been previously evaluated for gonococci. We evaluated the performance characteristics, ideal sequence analysis (tools and workflow for taxonomy, assembly, assembly improvement ("polishing") and mapping), phylogenomic analysis, and prediction of decreased susceptibility or resistance to recommended therapeutic antimicrobials in N. gonorrhoeae isolates using the Oxford Nanopore MinION sequencer. We sequenced the 2016 WHO gonococcal reference genomes/strains 40 and clinical gonococcal isolates, and evaluated the performance of ONT assemblies and hybrid assemblies including both ONT and Illumina reads. The assemblies were evaluated against the PacBio RS II-sequenced 2016 WHO gonococcal reference genomes (Bioproject PRJEB14020) 40 for sequence variation and to assess whether the assemblies were of sufficient quality for characterisation of AMR determinants and phylogenomic-based molecular epidemiology of gonococci.

N. gonorrhoeae isolates, culture, antimicrobial susceptibility testing and DNA isolation.
Twenty-eight gonococcal isolates, including the previously PacBio-sequenced 2016 WHO reference strains (n = 14) 40,41 and 14 clinical isolates, were examined. The clinical isolates were selected from 140 isolates from the RaDAR-Go project, Switzerland, in 2015-2016 42 . Selection was based on having decreased susceptibility or resistance to ceftriaxone and/or cefixime, resistance to azithromycin or ciprofloxacin, and to represent different AMR phenotypes. All clinical isolates were cultured as part of routine diagnostics, and no human clinical samples or information enabling patient identification was available for this study, so ethical approval was not required. All strains were subsequently cultivated and preserved as described 43 . We determined minimum inhibitory concentrations (MICs) of ceftriaxone, cefixime, azithromycin, spectinomycin, ciprofloxacin, tetracycline, and benzylpenicillin using the Etest (bioMérieux). Resistance breakpoints from the European Committee on Antimicrobial Susceptibility Testing (EUCAST; www.eucast.org/clinical_breakpoints/) were applied. Nitrocefin test (Thermo Fisher Scientific, Wilmington, DE, USA) was used to detect β-lactamase production. Genomic DNA was isolated using the Wizard Genomic DNA Purification Kit (Promega Corporation, Madison, WI, USA). Isolated DNA was quality controlled using fluorometric quantification (Qubit; Thermo Fisher Scientific) and electrophoresis (Tapestation; Agilent, Santa Clara, CA, USA).
Oxford Nanopore Technologies library preparation and MinION sequencing. Three µg of genomic DNA was sheared to an average fragment length of 8 kb with g-TUBES (Covaris, Woburn, WA, USA). Sequencing libraries were prepared according to the 2D library preparation protocol with the SQK-LSK208 2D ligation kit (Oxford Nanopore Technologies), including the DNA repair step with the NEBNext FFPE DNA repair module (New England Biolabs, Ipswich, MA, USA). The sequencing libraries were purified using MyOne C1 beads (Thermo Fisher Scientific), and 6 μl of sequencing library were loaded onto a R9.4 SpotON flow cell and sequenced with the MinION Mk 1B sequencing device (Oxford Nanopore Technologies) for 24 hours, including a top-up with additional 6 µl of DNA library after the first 6 hours of the sequencing run. Base-calling was performed using the Metrichor cloud software and all ONT reads (quality passed 1D and 2D) as well as only 2D ONT reads were then extracted using Poretools 44 for downstream analysis.
Illumina library preparation and sequencing. The sequencing libraries for all clinical isolates (n = 14) were prepared using Nextera XT DNA library preparation kit (Illumina, San Diego, CA, USA) and sequenced on the MiSeq Platform (Illumina), according to manufacturer's instructions, resulting in an average of 967,369 reads with an average read length of 257 bp after quality control and average coverage of 86.5× per base. The raw sequence files for the 2016 WHO gonococcal reference strains 40 were obtained from the European Nucleotide Archive (ENA; Bioproject PRJEB14020). The 2D ONT reads and Illumina reads were species confirmed using the online tool One Codex (www.onecodex.com). All sequenced reads are available from the European Nucleotide Archive with the following accession number: PRJEB25703.
Finally, to evaluate the accuracy of our assemblies, the ONT (corrected and non-corrected) and hybrid assemblies for the WHO reference strains were compared with the finished and closed 2016 WHO gonococcal reference genomes 40 . All 2D ONT reads were aligned to the genome sequence of the respective WHO reference strain, to be able to easily and directly determine the quality and accuracy using BWA-MEM (for mapping and phylogenetics) and Quast (v4.6.0) (for de novo assemblies). The top three assessed ONT de novo assemblies and the best hybrid assemblies were subjects for downstream analysis. Coding sequences (CDS) were annotated using Prokka (v1.12) on the chosen assemblies 54 .
Multiple-sequence alignment and phylogenetics. All Illumina and 2D ONT reads were mapped to the chromosome of FA1090 (GenBank: NC_002946.2) separately using BWA-MEM, the 2D ONT reads were specified using the nanopore option (-x -ont2d). All consensus sequences were merged into a multiple sequence alignment, single-nucleotide polymorphisms (SNPs) were called, and a maximum-likelihood phylogenetic tree based on SNPs was obtained using RAxML (version 8.2.8) 55 . This was also performed on Illumina sequences and 2D ONT sequences separately and the phylogeny was compared in a tanglegram.
Antimicrobial resistance determinants. To detect relevant gonococcal AMR determinants as quick and straight-forwardly as possible for future routine use of MinION sequencing, we evaluated several open-access and user-friendly web-based sequence analysis tools. The top three assessed ONT assemblies and the best hybrid assembly were examined using the Whole Genome Sequence Analysis (WGSA; www.wgsa.net) 21 and the PubMLST N. gonorrhoeae AMR (still in development) scheme (www.pubmlst.org) 56 . The focus was on AMR determinants for ceftriaxone and cefixime (penA, mtrR, penB), azithromycin (23S rDNA, mtrR), ciprofloxacin (gyrA), tetracycline (rpsJ, tet(M) (plasmid-mediated resistance)), and benzylpenicillin (penA, ponA, mtrR, penB, bla TEM (plasmid-mediated resistance)). We only detected the AMR determinants on the plasmids and no additional downstream analysis was performed on the plasmids. Furthermore, the 2D ONT reads were examined in our in house-customised CLC Genomics Workbench workflow to characterise the AMR determinants based on alleles in the N. gonorrhoeae Sequence Typing for Antimicrobial Resistance (NG-STAR) database (www.ngstar.canada.ca) 57 . The NG-STAR database lacked the possibility to run contig assemblies at the time of this study. The workflow employed de novo assembly without mapping and used BLAST with optimised algorithms for the identification, with the highest (%) hit being reported 58 .
The best performing tool of assembly and subsequent detection of AMR determinants was subsequently used to characterise the clinical isolates (n = 14).

Results
Overview of sequenced data. The 14 2016 WHO reference strains 40 and 14 clinical isolates were sequenced using the ONT MinION platform. In these MinION runs, we obtained between 105 Mb and 922 Mb of 2D sequences with longest reads ranging from 26.0 to 58.6 kb and an average read length of 4.6-6.4 kb. The number of passed 2D ONT reads ranged from 17,500 to 147,700 with an average of 33% of all reads. For evaluation of the ONT sequences and further downstream analysis (taxonomy, hybrid assemblies and phylogeny), we used the 2 × 100 bp Illumina HiSeq paired-end reads with >250 depth for each WHO reference genome 40 and 2 × 300 bp Illumina MiSeq paired-end reads with 80-100× depth for each clinical isolate. Taxonomy was defined for all reads using One Codex and 97-99% and 100% of all reads were classified as N. gonorrhoeae using ONT ( Fig. 1) and Illumina MiSeq, respectively.

Comparison of MinION de novo assemblies with the 2016 WHO gonococcal reference genomes.
Assemblies were compared and selected based on the Quast results when comparing to the WHO reference genomes 40 using the following criteria: lowest number of mismatches, misassemblies and number of contigs, and highest fraction of genome coverage (Supplementary Table). We selected and analysed ONT assemblies produced by Canu using only 2D ONT reads and polished with Nanopolish, Miniasm twice error corrected with Racon and polished with Nanopolish, and SMARTdenovo. Furthermore, the MaSuRCA hybrid assemblies (based on the ONT plus Illumina reads) were selected for downstream analysis (Table 1). Generally, ONT assemblies based on only the 2D ONT reads had a higher accuracy compared with assemblies using all reads (1D and 2D ONT reads) with fewer mismatches but very modest to no change in contiguity of the assemblies and fraction of the WHO reference genome covered. All error corrections were made using only 2D ONT reads and statistics such as mismatches and indels improved. Nanopolish had the best performance of the tested error correction tools and reduced the number of indels by up to 13 times. Assemblies produced with hybridSPAdes were more affected by the lower numbers of ONT reads, i.e. when using only 2D ONT reads. The hybrid assembly using SPAdes had less contiguity (up to 5.5 times more contigs) with higher numbers of mismatches (up to 4.4 more mismatches per 100 kb) and introductions of misassemblies (Supplementary Table). In general, the 2D ONT reads were used more efficiently by MaSuRCA for hybrid assembly, with overall improved statistics, especially for contiguity and mismatches.
Although the length of the assemblies based on the ONT sequences did not largely differ from the WHO reference genomes, the number of CDS was vastly different. The 2016 WHO gonococcal reference genomes include 2295-2450 CDS 40 . However, 3210-4712 CDS were identified using Prokka in our selected ONT assemblies. In contrast, Prokka identified 2228-2400 CDS in the MaSuRCA hybrid assemblies, which was more in line with the published WHO reference genomes (Table 1).
Error corrections using Pilon and Circlator, which is not designed for error correction, did not substantially improve the assemblies (Supplementary Table) and the assemblies error corrected with these tools were not used in the downstream analysis.
MinION for molecular epidemiology, including rapid outbreak investigations. The Illumina reads mapped to the PacBio-produced WHO reference genomes (n = 14) 40 with no SNPs detected except for WHO M (one SNP) and a range of 93.3% to 96.4% of the reads mapped to the respective PacBio reference genome. The ONT reads mapped with 2-58 SNPs per genome sequence and 95.2% to 98.3% of the reads mapped. By executing BWA-MEM with the nanopore -x ont2d option, no SNPs between the ONT reads and the respective reference genomes were detected, but we found variation in the number of 2D ONT reads and percentage mapped to the FA1090 reference genome (Fig. 2). Therefore, the BWA-MEM nanopore option was used when mapping all 2D ONT reads to the FA1090 reference genome for the phylogenetic analysis (Figs 3 and 4). By creating the phylogeny with the ONT-produced genomes (n = 28) and Illumina MiSeq genomes (n = 28) separately, we showed that the ONT sequenced libraries produced a phylogenetic tree topology that was comparable with the one using the Illumina dataset (Fig. 4). The main difference in the tree topology was that WHO G and WHO N were closer to the root when Illumina sequences were used for the phylogeny. Moreover, all identical isolates separately sequenced with ONT and Illumina clustered together, except for some (n = 7) of the clinical isolates where the Illumina sequences were more closely related. In four of these cases (11 cervix (C)/12 rectum (R) and 13 (R)/14 urethra (U)), the isolates were from the same patients but from different anatomical sites. For three isolates (5, 6 and 7), the Illumina sequences were also more closely related but all of the separately sequenced datasets were still on the same branch showing that these isolates were very similar (Fig. 3).
WGSA is an easy-to-use web interface where you drag-and-drop the assemblies into the web browser and the characterisation is done in minutes 21 . We used this interface for all chosen assemblies (n = 112) and the analysis of AMR determinants took 15 minutes. We observed good concordance overall with the verified AMR determinants in the reference genomes (Table 2), with the exception of the 23S rDNA macrolide resistance determinant C2611T, which was incorrectly identified in the Canu (polished) and Miniasm (polished) assemblies. No ONT assembly was able to correctly characterise penA in WHO W and none, including the MaSuRCA hybrid assembly, was able to correctly characterise penA in WHO K ( Table 2). The characterisation of plasmid-mediated AMR determinants was inaccurate in all four chosen assemblies for the detection of tet(M). Also, the chosen assemblies failed to detect bla TEM in WHO M, WHO O and WHO V and incorrectly reported presence of bla TEM in WHO G in the Miniasm (polished), SMARTdenovo and MaSuRCA assemblies ( Table 3). Out of the 2D ONT assemblies, the SMARTdenovo assembly performed the best. We also used the PubMLST database, which is another easy-to-use web-based tool, that can mine assembly data for various purposes. The database uses assemblies as input and delivers results in few minutes. The PubMLST database (www.pubmlst.org), using 2D ONT assemblies, was only able to characterise the mtrR AMR determinant and this was only done in the Canu assemblies. However, using the MaSuRCA hybrid assembly PubMLST had 100% concordance with the published reference genomes for detection of penA, mtrR, penB, ponA, 23S rDNA, gyrA, and tet(M). PubMLST does not include rpsJ in its gonococcal AMR characterisation module and also failed to detect the bla TEM gene in WHO M and WHO N ( Table 3).
The NG-STAR database 57 , included in the WHO CC in house-customised CLC Genomics Workbench, gave 100% concordance for the presence of all AMR determinants characterised in NG-STAR 57 with the reference genomes ( Table 2). All the plasmid-mediated AMR determinants ( rDNA, and parC were also correctly characterised. The customised CLC Genomics Workbench was therefore used to characterise the AMR determinants in the clinical isolates using only the 2D ONT reads.  (Table 4).

In house-customised CLC Genomics Workbench for characterisation of AMR determinants in
No isolate displayed spectinomycin resistance or produced β-lactamase. The in house-customised CLC Genomics Workbench including the NG-STAR database 57 correctly characterised all AMR determinants in the isolates  using only the 2D ONT reads and found 12, 1, and 13 of the 14 isolates carrying AMR mutations in penA (cefixime and ceftriaxone resistance determinant), 23S rDNA (azithromycin resistance determinant), and gyrA (ciprofloxacin resistance determinant), respectively. Twelve of the 14 isolates also harboured parC resistance mutations, which contribute to the high MICs in isolates with high-level ciprofloxacin resistance. Specific mutations in mtrR, resulting in an over-expression of the MtrCDE efflux pump, and porB1b (penB AMR determinant), causing a decreased influx of antimicrobials through PorB, were identified in 12 of the 14 isolates. No 16S rDNA resistance mutations (spectinomycin resistance determinant) or plasmid-mediated AMR determinants (high-level resistance to benzylpenicillin and tetracycline) were detected. Resistance mutation in rpsJ (decreased susceptibility/ resistance to tetracycline) was found in all isolates. There was 100% concordance between the AMR determinants detected using the 2D ONT reads only and Illumina reads. Consequently, resistance or decreased susceptibility to extended-spectrum cephalosporins, azithromycin, spectinomycin, ciprofloxacin, tetracycline, and benzylpenicillin could be accurately predicted using solely 2D ONT reads (Table 4).

Discussion
In the present study, we show that gonococcal genomes can be de novo assembled with high accuracy and contiguity running assemblies with MinION 2D ONT reads in combination with Illumina reads using MaSuRCA. These assemblies varied from one to eight contigs and could be further investigated using tools like Circlator 53 to obtain circular and finished gonococcal genomes. Furthermore, it is possible to obtain accurate de novo assemblies for AMR determinant detection by initially performing Illumina-based correction of individual ONT reads and subsequently using the reads for de novo assembly. It was also recently shown that complete hybrid genome assembly for gonococcal isolates can be obtained by combining short-read (Illumina) and long-read (ONT) sequence data 14 . However, this strategy makes the ONT platform dependent on additional Illumina sequencing. Nevertheless, this is another approach for producing accurate finished genomes and, accordingly, an alternative to PacBio sequencing of microbial reference genomes. When relying only on the ONT data, the assemblies contained high numbers of mismatches (27-607 mismatches per 100 kb) and indels (641-1125 indels per 100 kb), using even the best ONT assemblies obtained using SMARTdenovo, Canu, and Miniasm with 2D ONT reads. There are several possible reasons for the high and incorrect number of CDSs in the ONT assemblies using Prokka, including sequencing artefacts, homopolymers and overall high error rate (SNPs, insertions and deletions) in the 2D ONT reads causing the assemblies to contain incorrect internal stop codons and false pseudogenes. Generally, the 2D ONT reads generated more accurate assemblies than the 1D ONT reads and were used in the majority of our downstream analysis. We aimed to obtain highly accurate de novo assemblies to be able to rapidly identify relevant gonococcal AMR determinants to predict AMR, using only user-friendly and rapid online sequence analysis tools such as WGSA (www.wgsa.net) 21 and PubMLST (www.pubmlst.org). Using these online tools, the characterisation of AMR determinants was generally inaccurate using ONT assemblies and, as expected, highly accurate with hybrid assemblies. However, the in house-customised CLC Genomics Workbench workflow, including the NG-STAR database 57 , provided 100% concordance with the AMR determinants of the PacBio sequenced 2016 WHO gonococcal reference genomes 40 using only 2D ONT reads. We can also extract additional AMR determinants and in general genes of interest that are not included in the online tools, because our AMR database is fully customisable. The software workflow is based on extraction of the genes of interest in the customised database from a de novo assembly using BLAST algorithms optimised for each AMR determinant and reporting the highest hit, i.e. not only the 100% hit, in also genes with premature stop codons. This is essential because the ONT de novo assemblies are prone to errors. The main limitations of any CLC Genomics Workbench are that it is commercial (not an open-source online tool) and it has fairly high system requirements (16-32 GB RAM). Nevertheless, it has a simple general user interface for users who are not familiar with command line bioinformatics and is available to all widely used operating systems. Hopefully, open-source tools and public databases that can handle error-prone assemblies will be inspired by this approach and incorporate e.g. NG-STAR data 57 rendering costly software obsolete. NG-STAR can currently not use genome assemblies as input for determining the AMR profiles, which limits its possibilities when performing WGS.  Easy and rapid genome sequencing, using platforms such as the ONT MinION, in combination with algorithms that can appropriately predict AMR profiles using only genetic data could be very valuable in future surveillance of gonococcal AMR and spread of AMR gonococcal strains, nationally and internationally. Briefly, using  Table 3. Detection of plasmid-mediated antimicrobial resistance determinants in the 2016 WHO Neisseria gonorrhoeae reference strains (n = 14) 40,41 using Whole Genome Sequence Analysis (WGSA; www.wgsa.net) 21 and an in house-customised CLC Genomics Workbench. The discrepant results from the published references are highlighted in bold letters.   the customised CLC Genomics Workbench, we showed that a mosaic penA allele or a penA A501 substitution was detected in all clinical isolates expressing decreased susceptibility to extended-spectrum cephalosporins (ESCs; cefixime and ceftriaxone) and predicted decreased susceptibility to ESCs with 100% sensitivity and specificity. However, these AMR determinants can also be found in gonococcal isolates susceptible to ESCs 6,20 , so the specificity of the prediction of decreased susceptibility or resistance to ESCs is expected to be lower in studies that test larger numbers of isolates. Detection of the C2611T mutation in all four alleles of the 23S rRNA gene predicted azithromycin resistance in one clinical isolate (MIC = 4 mg/L), but not in two isolates with low-level azithromycin resistance (MIC = 1 mg/L). Finally, detection of the GyrA S91F mutation predicted the resistance in all (100%) the 13 clinical ciprofloxacin resistant isolates. Molecular tests, including WGS, will probably never completely replace culture-based AMR testing because these methods can only detect known AMR determinants; some gonococcal AMR determinants remain unknown and new AMR determinants are continuously evolving. However, molecular AMR testing can complement culture-based AMR testing. In settings in which molecular diagnostics is replacing culture-based diagnostics, rapid and accurate detection of known AMR determinants would significantly increase the number of samples tested for antimicrobial susceptibility. Ideally, molecular AMR testing should use WGS to identify all known and potentially new AMR determinants that can be used to predict AMR and even the MICs of antimicrobials 60 . The use of 2D ONT reads for phylogenomics of gonococci, e.g. for molecular epidemiological surveillance purposes, would be exceedingly valuable for outbreak investigations and monitoring in nearly real time in local, national and international gonococcal surveillance programmes, which aim to replace traditional, labour intensive and less accurate genotyping techniques such as N. gonorrhoeae multi-antigen sequence typing (NG-MAST) and multi-locus sequence typing (MLST) with WGS techniques 21 . Accordingly, we examined the accuracy of performing phylogenomic analysis using 2D ONT reads mapped to a reference genome using BWA-MEM with the nanopore option and building a phylogeny using the multiple sequence alignment. The ONT data produced a phylogenetic tree topology that was comparable with the one using the Illumina dataset and clustered all related isolates similarly (Figs 3 and 4). However, the number of isolates was limited and the genomic heterogeneity of the strains was high, which might have slightly biased our analysis.
Interestingly, the clustering of isolates was not highly dependent on the number of 2D ONT reads (Fig. 1). This suggests that the read length provides sufficient genome coverage even when a relatively low number of reads are available. For example, the lower number of 2D ONT reads for WHO G (Fig. 1) still provided a 99.96-100% genome fraction (Supplementary Table) and the read length compensated for the low number of reads. For the different WHO reference strains, 89.93-100% fraction of the PacBio-sequenced genomes were covered by the ONT reads (Supplementary Table). Consequently, for some purposes including gonococcal phylogenomic analysis, the MinION sequencer has likely collected enough reads in less than one day, which opens up new possibilities in running WGS to investigate gonorrhoea outbreaks in nearly real time, as well as directly from clinical samples.
The main benefits of using the MinION for genome sequencing included the very rapid turn-around time, high accessibility by connecting the small hand-held device to a laptop, and low cost of the sequencer. The main limitations included the lower accuracy and consistency compared to Illumina and PacBio sequencing reads and the high cost of each sequencing run. Cost remains a major obstacle to the use of MinION outside well-resourced research-focused laboratories. It is essential that the cost of MinION runs is reduced further, for example by using cheaper reagents or multiplexing libraries, and/or having international support and funding for use of MinION, or similar technologies, at national, regional and local levels in less-resourced settings. Furthermore, new bioinformatic tools for analysing the MinION reads are being developed continuously. Future studies should appropriately evaluate, for example, different basecalling tools (https://github.com/rrwick/Basecalling-comparison), assemblers such as Unicycler 61 , and mapping tools, e.g., minimap2 62 and GraphMap 63 .
In conclusion, we show, in the first MinION study for gonococci, that ONT sequences analysed with currently existing open-access, web-based sequence analysis tools are not sufficiently accurate to identify key gonococcal AMR determinants. However, using an appropriate analysis workflow such as an in house-developed CLC Genomics Workbench, we show that 2D ONT sequence data can be used for rapid and accurate identification of AMR determinants in N. gonorrhoeae isolates to predict decreased susceptibility or resistance to recommended therapeutic antimicrobials. We also show that 2D ONT sequence data can be useful for phylogenomics of N. gonorrhoeae, e.g. for molecular epidemiological investigations in nearly real time, and, using 2D ONT-Illumina hybrid assemblies, for producing contiguous assemblies and finished gonococcal reference genomes. The performance of MinION for WGS of N. gonorrhoeae directly from non-cultured nucleic acid amplification test samples now needs to evaluated.