Sequencing of SARS-CoV-2 in local transmission cases through Oxford Nanopore MinION platform from Karachi Pakistan.

INTRODUCTION: The first case of severe acute respiratory syndrome 2 (SARS-CoV-2) was imported to Pakistan in February 2020, since then 8,260 deaths have been witnessed. The virus has been constantly mutating and local transmission cases from different countries vary due to host dependent viral adaptation. Many distinct clusters of variant SARS-CoV-2 have been defined globally. In this study, the epidemiology of SARS-CoV-2 was studied and locally transmitted SARS-CoV-2 isolates from Karachi were sequenced to compared and identify any possible variants. METHODOLOGY: The real time PCR was performed on nasopharyngeal specimen to confirm SARS-CoV-2 with Orf 1ab and E gene as targets. The virus isolates were sequenced through oxford nanopore technology MinION platform. Isolates from the first and second wave of COVID-19 outbreak in Karachi were compared. RESULTS: The overall positivity rate for PCR was 26.24% with the highest number of positive cases in June. Approximately, 37.45% PCR positive subjects aged between 19-40 years. All the isolates belonged to GH clade and shared missense mutation D614G in spike protein linked to increased transmission rate worldwide. Another spike protein mutation A222V coexisted with D614G in the virus from the second wave of COVID-19. CONCLUSIONS: Based on the present findings it is suggested that the locally transmitted virus from Karachi varies from those reported from other parts of Pakistan. Slight variability was also observed between viruses from the first and second wave. Variability in any potential vaccine target may result in failed trials, therefore information on any local viral variants is always useful for effective vaccine design and/or selection.

Analyzing Modern Biomolecules: The Revolution of Nucleic-Acid Sequencing - Review.

Recent developments have revolutionized the study of biomolecules. Among them are molecular markers, amplification and sequencing of nucleic acids. The latter is classified into three generations. The first allows to sequence small DNA fragments. The second one increases throughput, reducing turnaround and pricing, and is therefore more convenient to sequence full genomes and transcriptomes. The third generation is currently pushing technology to its limits, being able to sequence single molecules, without previous amplification, which was previously impossible. Besides, this represents a new revolution, allowing researchers to directly sequence RNA without previous retrotranscription. These technologies are having a significant impact on different areas, such as medicine, agronomy, ecology and biotechnology. Additionally, the study of biomolecules is revealing interesting evolutionary information. That includes deciphering what makes us human, including phenomena like non-coding RNA expansion. All this is redefining the concept of gene and transcript. Basic analyses and applications are now facilitated with new genome editing tools, such as CRISPR. All these developments, in general, and nucleic-acid sequencing, in particular, are opening a new exciting era of biomolecule analyses and applications, including personalized medicine, and diagnosis and prevention of diseases for humans and other animals.

Whole-Genome Sequencing of Pathogens in Saliva : A Target-Enrichment Approach for SARS-CoV-2.

Outbreak analysis and transmission surveillance of viruses can be performed via whole-genome sequencing after viral isolation. Such techniques have recently been applied to characterize and monitor SARS-CoV-2 , the etiological agent of the COVID-19 pandemic. However, the isolation and culture of SARS-CoV-2 is time consuming and requires biosafety level 3 containment, which is not ideal for many resource-constrained settings. An alternate method, bait capture allows target enrichment and sequencing of the entire SARS-CoV-2 genome eliminating the need for viral culture. This method uses a set of hybridization probes known as "baits" that span the genome and provide sensitive, accurate, and minimal off-target hybridization. Baits can be designed to detect any known virus or bacteria in a wide variety of specimen types, including oral secretions. The bait capture method presented herein allows the whole genome of SARS-CoV-2 in saliva to be sequenced without the need to culture and provides an outline of bait design and bioinformatic analysis to guide a bioinformatician.

Implementation of whole genome sequencing for tuberculosis diagnostics in a low-middle income, high MDR-TB burden country.

Whole genome sequencing (WGS) is revolutionary for diagnostics of TB and its mutations associated with drug-resistances, but its uptake in low- and middle-income countries is hindered by concerns of implementation feasibility. Here, we provide a proof of concept for its successful implementation in such a setting. WGS was implemented in the Kyrgyz Republic. We estimated needs of up to 55 TB-WGS per week and chose the MiSeq platform (Illumina, USA) because of its capacity of up to 60 TB-WGS per week. The project's timeline was completed in 93-weeks. Costs of large equipment and accompanying costs were 222,065 USD and 8462 USD, respectively. The first 174 WGS costed 277 USD per sequence, but this was skewed by training inefficiencies. Based on real prices and presuming optimal utilization of WGS capacities, WGS costs could drop to 167 and 141 USD per WGS using MiSeq Reagent Kits v2 (500-cycles) and v3 (600-cycles), respectively. Five trainings were required to prepare the staff for autonomous WGS which cost 48,250 USD. External assessment confirmed excellent performance of WGS by the Kyrgyz laboratory in an interlaboratory comparison of 30 M. tuberculosis genomes showing complete agreeance of results.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples.

BACKGROUND: Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS. METHODS: Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance. RESULTS: At 100 parasites/µL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published. CONCLUSIONS: Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

Genome-Wide Identification of Host-Segregating Single-Nucleotide Polymorphisms for Source Attribution of Clinical Campylobacter coli Isolates.

Campylobacter is among the most common causes of gastroenteritis worldwide. Campylobacter jejuni and Campylobacter coli are the most common species causing human disease. DNA sequence-based methods for strain characterization have focused largely on C. jejuni, responsible for 80 to 90% of infections, meaning that C. coli epidemiology has lagged behind. Here, we have analyzed the genome of 450 C. coli isolates to determine genetic markers that can discriminate isolates sampled from 3 major reservoir hosts (chickens, cattle, and pigs). These markers then were applied to identify the source of infection of 147 C. coli strains from French clinical cases. Using STRUCTURE software, 259 potential host-segregating markers were revealed by probabilistic characterization of single-nucleotide polymorphism (SNP) frequency variation in strain collections from three different hosts. These SNPs were found in 41 genes or intergenic regions, mostly coding for proteins involved in motility and membrane functions. Source attribution of clinical isolates based on the differential presence of these markers confirmed chickens as the most common source of C. coli infection in France.IMPORTANCE Genome-wide and source attribution studies based on Campylobacter species have shown their importance for the understanding of foodborne infections. Although the use of multilocus sequence typing based on 7 genes from C. jejuni is a powerful method to structure populations, when applied to C. coli, results have not clearly demonstrated its robustness. Therefore, we aim to provide more accurate data based on the identification of single-nucleotide polymorphisms. Results from this study reveal an important number of host-segregating SNPs, found in proteins involved in motility, membrane functions, or DNA repair systems. These findings offer new, interesting opportunities for further study of C. coli adaptation to its environment. Additionally, the results demonstrate that poultry is potentially the main reservoir of C. coli in France.

Genopo: a nanopore sequencing analysis toolkit for portable Android devices.

The advent of portable nanopore sequencing devices has enabled DNA and RNA sequencing to be performed in the field or the clinic. However, advances in in situ genomics require parallel development of portable, offline solutions for the computational analysis of sequencing data. Here we introduce Genopo, a mobile toolkit for nanopore sequencing analysis. Genopo compacts popular bioinformatics tools to an Android application, enabling fully portable computation. To demonstrate its utility for in situ genome analysis, we use Genopo to determine the complete genome sequence of the human coronavirus SARS-CoV-2 in nine patient isolates sequenced on a nanopore device, with Genopo executing this workflow in less than 30 min per sample on a range of popular smartphones. We further show how Genopo can be used to profile DNA methylation in a human genome sample, illustrating a flexible, efficient architecture that is suitable to run many popular bioinformatics tools and accommodate small or large genomes. As the first ever smartphone application for nanopore sequencing analysis, Genopo enables the genomics community to harness this cheap, ubiquitous computational resource.

Comparison of commercially available whole-genome sequencing kits for variant detection in circulating cell-free DNA.

Circulating cell-free DNA (ccfDNA) has great potential for non-invasive diagnosis, prognosis and monitoring treatment of disease. However, a sensitive and specific whole-genome sequencing (WGS) method is required to identify novel genetic variations (i.e., SNVs, CNVs and INDELS) on ccfDNA that can be used as clinical biomarkers. In this article, five WGS methods were compared: ThruPLEX Plasma-seq, QIAseq cfDNA All-in-One, NEXTFLEX Cell Free DNA-seq, Accel-NGS 2 S PCR FREE DNA and Accel-NGS 2 S PLUS DNA. The Accel PCR-free kit did not produce enough material for sequencing. The other kits had significant common number of SNVs, INDELs and CNVs and showed similar results for SNVs and CNVs. The detection of variants and genomic signatures depends more upon the type of plasma sample rather than the WGS method used. Accel detected several variants not observed by the other kits. ThruPLEX seemed to identify more low-abundant SNVs and SNV signatures were similar to signatures observed with the QIAseq kit. Accel and NEXTFLEX had similar CNV and SNV signatures. These results demonstrate the importance of establishing a standardized workflow for identifying non-invasive candidate biomarkers. Moreover, the combination of variants discovered in ccfDNA using WGS has the potential to identify enrichment pathways, while the analysis of signatures could identify new subgroups of patients.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples.

BACKGROUND: Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined. METHODS: DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA. RESULTS: MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples. CONCLUSION: This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

The complete costs of genome sequencing: a microcosting study in cancer and rare diseases from a single center in the United Kingdom.

PURPOSE: The translation of genome sequencing into routine health care has been slow, partly because of concerns about affordability. The aspirational cost of sequencing a genome is $1000, but there is little evidence to support this estimate. We estimate the cost of using genome sequencing in routine clinical care in patients with cancer or rare diseases. METHODS: We performed a microcosting study of Illumina-based genome sequencing in a UK National Health Service laboratory processing 399 samples/year. Cost data were collected for all steps in the sequencing pathway, including bioinformatics analysis and reporting of results. Sensitivity analysis identified key cost drivers. RESULTS: Genome sequencing costs £6841 per cancer case (comprising matched tumor and germline samples) and £7050 per rare disease case (three samples). The consumables used during sequencing are the most expensive component of testing (68-72% of the total cost). Equipment costs are higher for rare disease cases, whereas consumable and staff costs are slightly higher for cancer cases. CONCLUSION: The cost of genome sequencing is underestimated if only sequencing costs are considered, and likely surpasses $1000/genome in a single laboratory. This aspirational sequencing cost will likely only be achieved if consumable costs are considerably reduced and sequencing is performed at scale.

Prediction of Minimum Inhibitory Concentrations of Antimicrobials for Neisseria gonorrhoeae Using Whole-Genome Sequencing.

Bacterial whole-genome sequencing is now increasingly available to researchers, reference laboratories and individual healthcare institutions. It can be possible to predict antimicrobial minimum inhibitory concentrations (MICs) for Neisseria gonorrhoeae using sequencing data, for many antimicrobials within one or two MIC doubling dilutions of the phenotypic value. With emerging rapid sequencing technologies, it may be possible in future to predict antimicrobial resistance faster than existing culture-based methods. Sequencing also provides insights into the genetic mechanisms underlying antimicrobial resistance, their spread in time and space, as well as the molecular epidemiology of the gonococcal strains.

Operational models and criteria for incorporating microbial whole genome sequencing in hospital microbiology - A systematic literature review.

BACKGROUND: Microbial whole genome sequencing (WGS) has many advantages over standard microbiological methods. However, it is not yet widely implemented in routine hospital diagnostics due to notable challenges. OBJECTIVES: The aim was to extract managerial, financial and clinical criteria supporting the decision to implement WGS in routine diagnostic microbiology, across different operational models of implementation in the hospital setting. METHODS: This was a systematic review of literature identified through PubMed and Web of Science. English literature studies discussing the applications of microbial WGS without limitation on publication date were eligible. A narrative approach for categorization and synthesis of the sources identified was adopted. RESULTS: A total of 98 sources were included. Four main alternative operational models for incorporating WGS in clinical microbiology laboratories were identified: full in-house sequencing and analysis, full outsourcing of sequencing and analysis and two hybrid models combining in-house/outsourcing of the sequencing and analysis components. Six main criteria (and multiple related sub-criteria) for WGS implementation emerged from our review and included cost (e.g. the availability of resources for capital and operational investment); manpower (e.g. the ability to provide training programmes or recruit trained personnel), laboratory infrastructure (e.g. the availability of supplies and consumables or sequencing platforms), bioinformatics requirements (e.g. the availability of valid analysis tools); computational infrastructure (e.g. the availability of storage space or data safety arrangements); and quality control (e.g. the existence of standardized procedures). CONCLUSIONS: The decision to incorporate WGS in routine diagnostics involves multiple, sometimes competing, criteria and sub-criteria. Mapping these criteria systematically is an essential stage in developing policies for adoption of this technology, e.g. using a multicriteria decision tool. Future research that will prioritize criteria and sub-criteria that were identified in our review in the context of operational models will inform decision-making at clinical and managerial levels with respect to effective implementation of WGS for routine use. Beyond WGS, similar decision-making challenges are expected with respect to future integration of clinical metagenomics.

Whole genome sequence of Vibrio cholerae directly from dried spotted filter paper.

BACKGROUND: Global estimates for cholera annually approximate 4 million cases worldwide with 95,000 deaths. Recent outbreaks, including Haiti and Yemen, are reminders that cholera is still a global health concern. Cholera outbreaks can rapidly induce high death tolls by overwhelming the capacity of health facilities, especially in remote areas or areas of civil unrest. Recent studies demonstrated that stool specimens preserved on filter paper facilitate molecular analysis of Vibrio cholerae in resource limited settings. Specimens preserved in a rapid, low-cost, safe and sustainable manner for sequencing provides previously unavailable data about circulating cholera strains. This may ultimately contribute new information to shape public policy response on cholera control and elimination. METHODOLOGY/PRINCIPAL FINDINGS: Whole genome sequencing (WGS) recovered close to a complete sequence of the V. cholerae O1 genome with satisfactory genome coverage from stool specimens enriched in alkaline peptone water (APW) and V. cholerae culture isolates, both spotted on filter paper. The minimum concentration of V. cholerae DNA sufficient to produce quality genomic information was 0.02 ng/µL. The genomic data confirmed the presence or absence of genes of epidemiological interest, including cholera toxin and pilus loci. WGS identified a variety of diarrheal pathogens from APW-enriched specimen spotted filter paper, highlighting the potential for this technique to explore the gut microbiome, potentially identifying co-infections, which may impact the severity of disease. WGS demonstrated that these specimens fit within the current global cholera phylogenetic tree, identifying the strains as the 7th pandemic El Tor. CONCLUSIONS: WGS results allowed for mapping of short reads from APW-enriched specimen and culture isolate spotted filter papers. This provided valuable molecular epidemiological sequence information on V. cholerae strains from remote, low-resource settings. These results identified the presence of co-infecting pathogens while providing rare insight into the specific V. cholerae strains causing outbreaks in cholera-endemic areas.

Tuberculosis at the animal-human interface in the Ugandan cattle corridor using a third-generation sequencing platform: a cross-sectional analysis study.

INTRODUCTION: Mycobacterial diseases are a significant source of disease burden, with Mycobacterium tuberculosis being the most common infectious cause of death worldwide. Given this, the emergence of antibiotic resistance in these species is of particular interest. By examining the epidemiology of mycobacteria in humans and cattle in an area of intense human-animal contact (the Ugandan cattle corridor [UCC]) and using a novel whole-genome sequencing technique to analyse organism diversity, this study will explore the role bidirectional transmission of mycobacteria plays in the local ecology, as well the significance of zoonotic Mycobacterium bovis in the human population. METHODS AND ANALYSIS: This ongoing study includes both a cross-sectional analysis of the UCC mycobacteria-positive population and novel laboratory-based research focused on differentiating the species causing M. tuberculosis complex-linked disease. We will use a third-generation sequencing platform (the MinION sequencer from Oxford Nanopore Technology), comparing data from a sample subset to the Illumina platform as a means of measuring viability of the MinION platform in this specific setting. Our full sample set will be sequenced on Illumina and the data will be used to perform epidemiological and phylogenetic analyses. ETHICS AND DISSEMINATION: Ethical approval was obtained from both the University of Minnesota IRB committee and the Ugandan National Council for Science and Technology-Research Ethics Committee. Samples were obtained after patients signed an informed consent indicating samples could be retained and used for research purposes. All samples are deidentified, with only basic demographic and geographic information being retained per national tuberculosis (TB) recording guidelines. Significant drug resistance results will be referred back to the local TB control officer to inform patient care. Final results of the study will be submitted to infectious disease-specific journals and will be submitted to the annual Infectious Diseases Society of America (IDSA) meeting.

Sequencing of human genomes with nanopore technology.

Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. Here we evaluate the potential of the Oxford Nanopore Technologies (ONT) MinION long-read sequencer for routine WGS by sequencing the reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome and severe immune dysregulation. We develop and apply a novel reference panel-free analytical method to infer and then exploit phase information which improves single-nucleotide variant (SNV) calling performance from otherwise modest levels. In the clinical sample, we identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they lie on the same paternal haplotype. Whilst consensus SNV-calling error rates from ONT data remain substantially higher than those from short-read methods, we demonstrate the substantial benefits of analytical innovation. Ongoing improvements to base-calling and SNV-calling methodology must continue for nanopore sequencing to establish itself as a primary method for clinical WGS.

Sequencing-based methods and resources to study antimicrobial resistance.

Antimicrobial resistance extracts high morbidity, mortality and economic costs yearly by rendering bacteria immune to antibiotics. Identifying and understanding antimicrobial resistance are imperative for clinical practice to treat resistant infections and for public health efforts to limit the spread of resistance. Technologies such as next-generation sequencing are expanding our abilities to detect and study antimicrobial resistance. This Review provides a detailed overview of antimicrobial resistance identification and characterization methods, from traditional antimicrobial susceptibility testing to recent deep-learning methods. We focus on sequencing-based resistance discovery and discuss tools and databases used in antimicrobial resistance studies.

Preimplantation genetic diagnosis and screening (PGD/S) using a semiconductor sequencing platform.

BACKGROUND: Recent advances in semiconductor sequencing platform (SSP) have provided new methods for preimplantation genetic diagnosis/screening (PGD/S). The present study aimed to evaluate the applicability and efficiency of SSP in PGD/S. METHODS: The artificial positive single-cell-like DNAs and normal single-cell samples were chosen to test our semiconductor sequencing platform for preimplantation genetic diagnosis/screening (SSP-PGD/S) method with two widely used whole-genome amplification (WGA) kits. A total of 557 single blastomeres were collected from in vitro fertilization (IVF) couples, and their WGA products were processed and analyzed by our SSP-PGD/S method in comparison with array comparative genomic hybridization (array-CGH). RESULTS: Our SSP-PGD/S method indicated high compatibilities with two commercial WGA kits. For 557 single blastomeres, our method with four million reads in average could detect 24-chromosome aneuploidies as well as microdeletion/microduplication of the size over 4 Mb, providing 100% consistent conclusion with array-CGH method in the classification of whether it was transplantable. CONCLUSIONS: Our studies suggested that SSP-PGD/S represents a valuable alternative to array-CGH and brought PGD/S into a new era of more rapid, accurate, and economic.

Applying Rapid Whole-Genome Sequencing To Predict Phenotypic Antimicrobial Susceptibility Testing Results among Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates.

Standard antimicrobial susceptibility testing (AST) approaches lead to delays in the selection of optimal antimicrobial therapy. Here, we sought to determine the accuracy of antimicrobial resistance (AMR) determinants identified by Nanopore whole-genome sequencing in predicting AST results. Using a cohort of 40 clinical isolates (21 carbapenemase-producing carbapenem-resistant Klebsiella pneumoniae, 10 non-carbapenemase-producing carbapenem-resistant K. pneumoniae, and 9 carbapenem-susceptible K. pneumoniae isolates), three separate sequencing and analysis pipelines were performed, as follows: (i) a real-time Nanopore analysis approach identifying acquired AMR genes, (ii) an assembly-based Nanopore approach identifying acquired AMR genes and chromosomal mutations, and (iii) an approach using short-read correction of Nanopore assemblies. The short-read correction of Nanopore assemblies served as the reference standard to determine the accuracy of Nanopore sequencing results. With the real-time analysis approach, full annotation of acquired AMR genes occurred within 8 h from subcultured isolates. Assemblies sufficient for full resistance gene and single-nucleotide polymorphism annotation were available within 14 h from subcultured isolates. The overall agreement of genotypic results and anticipated AST results for the 40 K. pneumoniae isolates was 77% (range, 30% to 100%) and 92% (range, 80% to 100%) for the real-time approach and the assembly approach, respectively. Evaluating the patients contributing the 40 isolates, the real-time approach and assembly approach could shorten the median time to effective antibiotic therapy by 20 h and 26 h, respectively, compared to standard AST. Nanopore sequencing offers a rapid approach to both accurately identify resistance mechanisms and to predict AST results for K. pneumoniae isolates. Bioinformatics improvements enabling real-time alignment, coupled with rapid extraction and library preparation, will further enhance the accuracy and workflow of the Nanopore real-time approach.

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.