Enrichment by hybridisation of long DNA fragments for Nanopore sequencing

Advances in genotypic antimicrobialresistance testing: a comprehensive review

Antimicrobial resistance (AMR) represents a substantial threat to global public health, complicating the treatment of common infections and leading to prolonged illness and escalated healthcare expenses. To effectively combat AMR, timely and accurate detection is crucial for AMR surveillance and individual-based therapy. Phenotypic antibiotic resistance testing (AST) has long been considered the gold standard in clinical applications, serving as the foundation for clinical AMR diagnosis and optimized therapy. It has significantly contributed to ensuring patients' health and the development of novel antimicrobials. Despite advancements in automated culture-based AST technologies, inherent limitations impede the widespread use of phenotypic AST in AMR surveillance. Genotypic AST technologies offer a promising alternative option, exhibiting advantages of rapidity, high sensitivity, and specificity. With the continuous advancement and expanding applications of genotypic AST technologies, such as microfluidics, mass spectrometry, and high-resolution melting curve analysis, new vigor has been injected into the development and clinical implementation of genotypic AST technologies. In this narrative review, we discuss the principles, applications, and advancements of emerging genotypic AST methods in clinical settings. The comprehensive review aims to highlight the significant scientific potential of emerging genotypic AST technologies in clinical AMR diagnosis, providing insights to enhance existing methods and explore novel approaches.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Nano-Strainer: A workflow for the identification of single-copy nuclear loci for plant systematic studies, using target capture kits and Oxford Nanopore long reads

In modern plant systematics, target enrichment enables simultaneous analysis of hundreds of genes. However, when dealing with reticulate or polyploidization histories, few markers may suffice, but often are required to be single-copy, a condition that is not necessarily met with commercial capture kits. Also, large genome sizes can render target capture ineffective, so that amplicon sequencing would be preferable; however, knowledge about suitable loci is often missing. Here, we present a comprehensive workflow for the identification of putative single-copy nuclear markers in a genus of interest, by mining a small dataset from target capture using a few representative taxa. The proposed pipeline assesses sequence variability contained in the data from targeted loci and assigns reads to their respective genes, via a combined BLAST/clustering procedure. Cluster consensus sequences are then examined based on four pre-defined criteria presumably indicative for absence of paralogy. This is done by calculating four specialized indices; loci are ranked according to their performance in these indices, and top-scoring loci are considered putatively single- or low copy. The approach can be applied to any probe set. As it relies on long reads, the present contribution also provides template workflows for processing Nanopore-based target capture data. Obtained markers are further tested and then entered into amplicon sequencing. For the detection of possibly remaining paralogy in these data, which might occur in groups with rampant paralogy, we also employ the long-read assembly tool canu. In diploid representatives of the young Compositae genus Leucanthemum, characterized by high levels of polyploidy, our approach resulted in successful amplification of 13 loci. Modifications to remove traces of paralogy were made in seven of these. A species tree from the markers correctly reproduced main relationships in the genus, however, at low resolution. The presented workflow has the potential to valuably support phylogenetic research, for example in polyploid plant groups.

Rapid Real-time Squiggle Classification for Read until using RawMap

ReadUntil enables Oxford Nanopore Technology's (ONT) sequencers to selectively sequence reads of target species in real-time. This enables efficient microbial enrichment for applications such as microbial abundance estimation and is particularly beneficial for metagenomic samples with a very high fraction of non-target reads (> 99% can be human reads). However, read-until requires a fast and accurate software filter that analyzes a short prefix of a read and determines if it belongs to a microbe of interest (target) or not. The baseline Read Until pipeline uses a deep neural network-based basecaller called Guppy and is slow and inaccurate for this task (~60% of bases sequenced are unclassified). We present RawMap, an efficient CPU-only microbial species-agnostic Read Until classifier for filtering non-target human reads in the squiggle space. RawMap uses a Support Vector Machine (SVM), which is trained to distinguish human from microbe using non-linear and non-stationary characteristics of ONT's squiggle output (continuous electrical signals). Compared to the baseline Read Until pipeline, RawMap is a 1327X faster classifier and significantly improves the sequencing time and cost, and compute time savings. We show that RawMap augmented pipelines reduce sequencing time and cost by ~24% and computing cost by 22%. Additionally, since RawMap is agnostic to microbial species, it can also classify microbial species it is not trained on. We also discuss how RawMap may be used as an alternative to the RT-PCR test for viral load quantification of SARS-CoV-2.

Nanopore Sequencing for Mycobacterium tuberculosis: a Critical Review of the Literature, New Developments, and Future Opportunities

The next-generation, short-read sequencing technologies that generate comprehensive, whole-genome data with single nucleotide resolution have already advanced tuberculosis diagnosis, treatment, surveillance, and source investigation. Their high costs, tedious and lengthy processes, and large equipment remain major hurdles for research use in high tuberculosis burden countries and implementation into routine care. The portable next-generation sequencing devices developed by Oxford Nanopore Technologies (ONT) are attractive alternatives due to their long-read sequence capability, compact low-cost hardware, and continued improvements in accuracy and throughput. A systematic review of the published literature demonstrated limited uptake of ONT sequencing in tuberculosis research and clinical care. Of the 12 eligible articles presenting ONT sequencing data on at least one Mycobacterium tuberculosis sample, four addressed software development for long-read ONT sequencing data with potential applications for M. tuberculosis. Only eight studies presented results of ONT sequencing of M. tuberculosis, of which five performed whole-genome and three did targeted sequencing. Based on these findings, we summarize the standard processes, reflect on the current limitations of ONT sequencing technology, and the research needed to overcome the main hurdles. The low capital cost, portable nature and continued improvement in the performance of ONT sequencing make it an attractive option for sequencing for research and clinical care, but limited data are available on its application in the tuberculosis field. Important research investment is needed to unleash the full potential of ONT sequencing for tuberculosis research and care.

GRAde: a long-read sequencing approach to efficiently identifying the CYP11B1/CYP11B2 chimeric form in patients with glucocorticoid-remediable aldosteronism

BACKGROUND

Glucocorticoid-remediable aldosteronism (GRA) is a form of heritable hypertension caused by a chimeric fusion resulting from unequal crossing over between 11β-hydroxylase (CYP11B1) and aldosterone synthase (CYP11B2), which are two genes with similar sequences. Different crossover patterns of the CYP11B1 and CYP11B2 chimeric genes may be associated with a variety of clinical presentations. It is therefore necessary to develop an efficient approach for identifying the differences between the hybrid genes of a patient with GRA.

RESULTS

We developed a long-read analysis pipeline named GRAde (GRA deciphering), which utilizes the nonidentical bases in the CYP11B1 and CYP11B2 genomic sequences to identify and visualize the chimeric form. We sequenced the polymerase chain reaction (PCR) products of the CYP11B1/CYP11B2 chimeric gene from 36 patients with GRA using the Nanopore MinION device and analyzed the sequences using GRAde. Crossover events were identified for 30 out of the 36 samples. The crossover sites appeared in the region exhibiting high sequence similarity between CYP11B1 and CYP11B2, and 53.3% of the cases were identified as having a gene conversion in intron 2. More importantly, there were six cases for whom the PCR products indicated a chimeric gene, but the GRAde results revealed no crossover pattern. The crossover regions were further verified by Sanger sequencing analysis.

CONCLUSIONS

PCR-based target enrichment followed by long-read sequencing is an efficient and precise approach to dissecting complex genomic regions, such as those involved in GRA mutations, which could be directly applied to clinical diagnosis. The scripts of GRAde are available at https://github.com/hsu-binfo/GRAde .

Long range PCR-based deep sequencing for haplotype determination in mixed HCMV infections

BACKGROUND

Short read sequencing has been used extensively to decipher the genome diversity of human cytomegalovirus (HCMV) strains, but falls short to reveal individual genomes in mixed HCMV strain populations. Novel third-generation sequencing platforms offer an extended read length and promise to resolve how distant polymorphic sites along individual genomes are linked. In the present study, we established a long amplicon PacBio sequencing workflow to identify the absolute and relative quantities of unique HCMV haplotypes spanning over multiple hypervariable sites in mixtures. Initial validation of this approach was performed with defined HCMV DNA templates derived from cell-culture enriched viruses and was further tested for its suitability on patient samples carrying mixed HCMV infections.

RESULTS

Total substitution and indel error rate of mapped reads ranged from 0.17 to 0.43% depending on the stringency of quality trimming. Artificial HCMV DNA mixtures were correctly determined down to 1% abundance of the minor DNA source when the total HCMV DNA input was 4 × 104 copies/ml. PCR products of up to 7.7 kb and a GC content < 55% were efficiently generated when DNA was directly isolated from patient samples. In a single sample, up to three distinct haplotypes were identified showing varying relative frequencies. Alignments of distinct haplotype sequences within patient samples showed uneven distribution of sequence diversity, interspersed by long identical stretches. Moreover, diversity estimation at single polymorphic regions as assessed by short amplicon sequencing may markedly underestimate the overall diversity of mixed haplotype populations.

CONCLUSIONS

Quantitative haplotype determination by long amplicon sequencing provides a novel approach for HCMV strain characterisation in mixed infected samples which can be scaled up to cover the majority of the genome by multi-amplicon panels. This will substantially improve our understanding of intra-host HCMV strain diversity and its dynamic behaviour.

Assessment of Viral Targeted Sequence Capture Using Nanopore Sequencing Directly from Clinical Samples

Shotgun metagenomic sequencing (SMg) enables the simultaneous detection and characterization of viruses in human, animal and environmental samples. However, lack of sensitivity still poses a challenge and may lead to poor detection and data acquisition for detailed analysis. To improve sensitivity, we assessed a broad scope targeted sequence capture (TSC) panel (ViroCap) in both human and animal samples. Moreover, we adjusted TSC for the Oxford Nanopore MinION and compared the performance to an SMg approach. TSC on the Illumina NextSeq served as the gold standard. Overall, TSC increased the viral read count significantly in challenging human samples, with the highest genome coverage achieved using the TSC on the MinION. TSC also improved the genome coverage and sequencing depth in clinically relevant viruses in the animal samples, such as influenza A virus. However, SMg was shown to be adequate for characterizing a highly diverse animal virome. TSC on the MinION was comparable to the NextSeq and can provide a valuable alternative, offering longer reads, portability and lower initial cost. Developing new viral enrichment approaches to detect and characterize significant human and animal viruses is essential for the One Health Initiative.

Whole-Genome Sequencing of Human Enteroviruses from Clinical Samples by Nanopore Direct RNA Sequencing

Enteroviruses are small RNA viruses that affect millions of people each year by causing an important burden of disease with a broad spectrum of symptoms. In routine diagnostic laboratories, enteroviruses are identified by PCR-based methods, often combined with partial sequencing for genotyping. In this proof-of-principle study, we assessed direct RNA sequencing (DRS) using nanopore sequencing technology for fast whole-genome sequencing of viruses directly from clinical samples. The approach was complemented by sequencing the corresponding viral cDNA via Illumina MiSeq sequencing. DRS of total RNA extracted from three different enterovirus-positive stool samples produced long RNA fragments, covering between 59% and 99.6% of the most similar reference genome sequences. The identification of the enterovirus sequences in the samples was confirmed by short-read cDNA sequencing. Sequence identity between DRS and Illumina MiSeq enterovirus consensus sequences ranged between 94% and 97%. Here, we show that nanopore DRS can be used to correctly identify enterovirus genotypes from patient stool samples with high viral load and that the approach also provides rich metatranscriptomic information on sample composition for all life domains.

Metagenomic Nanopore Sequencing of Influenza Virus Direct from Clinical Respiratory Samples

Influenza is a major global public health threat as a result of its highly pathogenic variants, large zoonotic reservoir, and pandemic potential. Metagenomic viral sequencing offers the potential for a diagnostic test for influenza virus which also provides insights on transmission, evolution, and drug resistance and simultaneously detects other viruses. We therefore set out to apply the Oxford Nanopore Technologies sequencing method to metagenomic sequencing of respiratory samples. We generated influenza virus reads down to a limit of detection of 102 to 103 genome copies/ml in pooled samples, observing a strong relationship between the viral titer and the proportion of influenza virus reads (P = 4.7 × 10-5). Applying our methods to clinical throat swabs, we generated influenza virus reads for 27/27 samples with mid-to-high viral titers (cycle threshold [CT ] values, <30) and 6/13 samples with low viral titers (CT values, 30 to 40). No false-positive reads were generated from 10 influenza virus-negative samples. Thus, Nanopore sequencing operated with 83% sensitivity (95% confidence interval [CI], 67 to 93%) and 100% specificity (95% CI, 69 to 100%) compared to the current diagnostic standard. Coverage of full-length virus was dependent on sample composition, being negatively influenced by increased host and bacterial reads. However, at high influenza virus titers, we were able to reconstruct >99% complete sequences for all eight gene segments. We also detected a human coronavirus coinfection in one clinical sample. While further optimization is required to improve sensitivity, this approach shows promise for the Nanopore platform to be used in the diagnosis and genetic analysis of influenza virus and other respiratory viruses.

Nanopore Targeted Sequencing for Rapid Gene Mutations Detection in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) clinical settings cannot do without molecular testing to confirm or rule out predictive biomarkers for prognostic stratification, in order to initiate or withhold targeted therapy. Next generation sequencing offers the advantage of the simultaneous investigation of numerous genes, but these methods remain expensive and time consuming. In this context, we present a nanopore-based assay for rapid (24 h) sequencing of six genes (NPM1, FLT3, CEBPA, TP53, IDH1 and IDH2) that are recurrently mutated in AML. The study included 22 AML patients at diagnosis; all data were compared with the results of S5 sequencing, and discordant variants were validated by Sanger sequencing. Nanopore approach showed substantial advantages in terms of speed and low cost. Furthermore, the ability to generate long reads allows a more accurate detection of longer FLT3 internal tandem duplications and phasing double CEBPA mutations. In conclusion, we propose a cheap, rapid workflow that can potentially enable all basic molecular biology laboratories to perform detailed targeted gene sequencing analysis in AML patients, in order to define their prognosis and the appropriate treatment.

Human Herpesvirus Sequencing in the Genomic Era: The Growing Ranks of the Herpetic Legion

The nine human herpesviruses are some of the most ubiquitous pathogens worldwide, causing life-long latent infection in a variety of different tissues. Human herpesviruses range from mild childhood infections to known tumour viruses and 'trolls of transplantation'. Epstein-Barr virus was the first human herpesvirus to have its whole genome sequenced; GenBank now includes thousands of herpesvirus genomes. This review will cover some of the recent advances in our understanding of herpesvirus diversity and disease that have come about as a result of new sequencing technologies, such as target enrichment and long-read sequencing. It will also look at the problem of resolving mixed-genotype infections, whether with short or long-read sequencing methods; and conclude with some thoughts on the future of the field as herpesvirus population genomics becomes a reality.

Long-fragment targeted capture for long-read sequencing of plastomes

PREMISE

Third-generation sequencing methods generate significantly longer reads than those produced using alternative sequencing methods. This provides increased possibilities for the study of biodiversity, phylogeography, and population genetics. We developed a protocol for in-solution enrichment hybridization capture of long DNA fragments applicable to complete plastid genomes.

METHODS AND RESULTS

The protocol uses cost-effective in-house probes developed via long-range PCR and was used in six non-model monocot species (Poaceae: African rice, pearl millet, fonio; and three palm species). DNA was extracted from fresh and silica gel-dried leaves. Our protocol successfully captured long-read plastome fragments (3151 bp median on average), with an enrichment rate ranging from 15% to 98%. DNA extracted from silica gel-dried leaves led to low-quality plastome assemblies when compared to DNA extracted from fresh tissue.

CONCLUSIONS

Our protocol could also be generalized to capture long sequences from specific nuclear fragments.

A novel CRISPR/Cas9 associated technology for sequence-specific nucleic acid enrichment

Massively parallel sequencing technologies have made it possible to generate large quantities of sequence data. However, as research-associated information is transferred into clinical practice, cost and throughput constraints generally require sequence-specific targeted analyses. Therefore, sample enrichment methods have been developed to meet the needs of clinical sequencing applications. However, current amplification and hybrid capture enrichment methods are limited in the contiguous length of sequences for which they are able to enrich. PCR based amplification also loses methylation data and other native DNA features. We have developed a novel technology (Negative Enrichment) where we demonstrate targeting long (>10 kb) genomic regions of interest. We use the specificity of CRISPR-Cas9 single guide RNA (Cas9/sgRNA) complexes to define 5' and 3' termini of sequence-specific loci in genomic DNA, targeting 10 to 36 kb regions. The complexes were found to provide protection from exonucleases, by protecting the targeted sequences from degradation, resulting in enriched, double-strand, non-amplified target sequences suitable for next-generation sequencing library preparation or other downstream analyses.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.

Rapid Sequencing of Multiple RNA Viruses in Their Native Form

Long-read nanopore sequencing by a MinION device offers the unique possibility to directly sequence native RNA. We combined an enzymatic poly-A tailing reaction with the native RNA sequencing to (i) sequence complex population of single-stranded (ss)RNA viruses in parallel, (ii) detect genome, subgenomic mRNA/mRNA simultaneously, (iii) detect a complex transcriptomic architecture without the need for assembly, (iv) enable real-time detection. Using this protocol, positive-ssRNA, negative-ssRNA, with/without a poly(A)-tail, segmented/non-segmented genomes were mixed and sequenced in parallel. Mapping of the generated sequences on the reference genomes showed 100% length recovery with up to 97% identity. This work provides a proof of principle and the validity of this strategy, opening up a wide range of applications to study RNA viruses.

Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations

BACKGROUND

Human cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies.

RESULTS

The MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions.

CONCLUSIONS

Third generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis.

Influenza A/H3N2 virus infection in immunocompromised ferrets and emergence of antiviral resistance

Influenza viruses can cause severe life threatening infections in high-risk patients, including young children, the elderly and patients with compromised immunity due to underlying medical conditions or immunosuppressive treatment. The impaired immunity of these patients causes prolonged virus infection and combined with antiviral treatment facilitates the emergence of viruses with resistance mutations. The diverse nature of their immune status makes them a challenging group to study the impact of influenza virus infection and the efficacy of antiviral therapy. Immunocompromised ferrets may represent a suitable animal model to assess influenza virus infection and antiviral treatment strategies in immunocompromised hosts. Here, ferrets were given a daily oral solution of mycophenolate mofetil, tacrolimus and prednisolone sodium phosphate to suppress their immune system. Groups of immunocompromised and immunocompetent ferrets were inoculated with an A/H3N2 influenza virus and were subsequently treated with Oseltamivir or left untreated. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) was performed on the throat and nose specimens to study virus replication during the course of infection. All immunocompromised ferrets had prolonged presence of viral RNA and a higher total amount of virus shedding compared to the immunocompetent ferrets. Although Oseltamivir reduced the total amount of virus shedding from the nose and throat of treated ferrets, it also resulted in the emergence of the neuraminidase R292K resistance substitution in all these animals, as determined by mutation specific RT-PCR and next-generation sequencing. No additional mutations that could be associated with the emergence of the R292K resistance mutation were detected. The immunocompromised ferret model can be used to study A/H3N2 virus shedding and is a promising model to study new antiviral strategies and the emergence of antiviral resistance in immunocompromised hosts.

Capture Hybridization of Long-Range DNA Fragments for High-Throughput Sequencing

Capture hybridization coupled with high-throughput sequencing (HTS) has become one of the most popular approaches to address some scientific problems not only for fundamental evolution but also for ecology and human disease in recent years. However, the technical problem of limited probe capture ability affects its widespread application. Here, we propose to capture hybridize long-range DNA fragments for HTS (termed LR-LCH). We provide a case of three amphibian samples to examine LR-LCH with 2 kb libraries and comparison of standard capture hybridization with 480 bp libraries. Capture sensitivity increased from an average 13.57% of standard capture hybridization to an average 19.80% of LR-LCH; capture efficiency also increased from an average 72.56% of standard capture hybridization to an average 97.71% of LR-LCH. These indicate that longer fragments in the library generally contain both relatively variable regions and relatively conservative regions. The divergent parts of target DNA are enriched along with conservative parts of DNA sequence that effectively captured during hybridization. We present a protocol that allows users to overcome the low capture sensitivity problem for high divergent regions.

Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches

This study aimed to assess the feasibility of using the Oxford Nanopore Technologies (ONT) MinION long-read sequencer in reconstructing fully closed plasmid sequences from eight Enterobacteriaceae isolates of six different species with plasmid populations of varying complexity. Species represented were Escherichia coli, Klebsiella pneumoniae, Citrobacter freundii, Enterobacter cloacae, Serratia marcescens and Klebsiella oxytoca, with plasmid populations ranging from 1–11 plasmids with sizes of 2–330 kb. Isolates were sequenced using Illumina (short-read) and ONT’s MinION (long-read) platforms, and compared with fully resolved PacBio (long-read) sequence assemblies for the same isolates. We compared the performance of different assembly approaches including SPAdes, plasmidSPAdes, hybridSPAdes, Canu, Canu+Pilon (canuPilon) and npScarf in recovering the plasmid structures of these isolates by comparing with the gold-standard PacBio reference sequences. Overall, canuPilon provided consistently good quality assemblies both in terms of assembly statistics (N50, number of contigs) and assembly accuracy [presence of single nucleotide polymorphisms (SNPs)/indels with respect to the reference sequence]. For plasmid reconstruction, Canu recovered 70 % of the plasmids in complete contigs, and combining three assembly approaches (Canu or canuPilon, hybridSPAdes and plasmidSPAdes) resulted in a total 78 % recovery rate for all the plasmids. The analysis demonstrated the potential of using MinION sequencing technology to resolve important plasmid structures in Enterobacteriaceae species independent of and in conjunction with Illumina sequencing data. A consensus assembly derived from several assembly approaches could present significant benefit in accurately resolving the greatest number of plasmid structures.

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

Genome sequencing has become a powerful tool for studying emerging infectious diseases; however, genome sequencing directly from clinical samples (i.e., without isolation and culture) remains challenging for viruses such as Zika, for which metagenomic sequencing methods may generate insufficient numbers of viral reads. Here we present a protocol for generating coding-sequence-complete genomes, comprising an online primer design tool, a novel multiplex PCR enrichment protocol, optimized library preparation methods for the portable MinION sequencer (Oxford Nanopore Technologies) and the Illumina range of instruments, and a bioinformatics pipeline for generating consensus sequences. The MinION protocol does not require an Internet connection for analysis, making it suitable for field applications with limited connectivity. Our method relies on multiplex PCR for targeted enrichment of viral genomes from samples containing as few as 50 genome copies per reaction. Viral consensus sequences can be achieved in 1-2 d by starting with clinical samples and following a simple laboratory workflow. This method has been successfully used by several groups studying Zika virus evolution and is facilitating an understanding of the spread of the virus in the Americas. The protocol can be used to sequence other viral genomes using the online Primal Scheme primer designer software. It is suitable for sequencing either RNA or DNA viruses in the field during outbreaks or as an inexpensive, convenient method for use in the lab.