OUP accepted manuscript

Impact of microbiological molecular methodologies on adaptive sampling using nanopore sequencing in metagenomic studies

INTRODUCTION

Metagenomics, the genomic analysis of all species present within a mixed population, is an important tool used for the exploration of microbiomes in clinical and environmental microbiology. Whilst the development of next-generation sequencing, and more recently third generation long-read approaches such as nanopore sequencing, have greatly advanced the study of metagenomics, recovery of unbiased material from microbial populations remains challenging. One promising advancement in genomic sequencing from Oxford Nanopore Technologies (ONT) is adaptive sampling, which enables real-time enrichment or depletion of target sequences. As sequencing technologies continue to develop, and advances such as adaptive sampling become common techniques within the microbiological toolkit, it is essential to evaluate the benefits of such advancements to metagenomic studies, and the impact of methodological choices on research outcomes.

AIM AND METHODS

Given the rapid development of sequencing tools and chemistry, this study aimed to demonstrate the impacts of choice of DNA extraction kit and sequencing chemistry on downstream metagenomic analyses. We first explored the quality and accuracy of 16S rRNA amplicon sequencing for DNA extracted from the ZymoBIOMICS Microbial Community Standard, using a range of commercially available DNA extraction kits to understand the effects of different kit biases on assessment of microbiome composition. We next compared the quality and accuracy of metagenomic analyses for two nanopore-based ligation chemistry kits with differing levels of base-calling error; the older and more error-prone (~ 97% accuracy) LSK109 chemistry, and newer more accurate (~ 99% accuracy) LSK112 Q20 + chemistry. Finally, we assessed the impact of the nanopore sequencing chemistry version on the output of the novel adaptive sampling approach for real-time enrichment of the genome for the yeast Saccharomyces cerevisiae from the microbial community.

RESULTS

Firstly, DNA extraction kit methodology impacted the composition of the yield, with mechanical bead-beating methodologies providing the least biased picture due to efficient lysis of Gram-positive microbes present in the community standard, with differences in bead-beating methodologies also producing variation in composition. Secondly, whilst use of the Q20 + nanopore sequencing kit chemistry improved the base-calling data quality, the resulting metagenomic assemblies were not significantly improved based on common metrics and assembly statistics. Most importantly, we demonstrated the effective application of adaptive sampling for enriching a low-abundance genome within a metagenomic sample. This resulted in a 5-7-fold increase in target enrichment compared to non-adaptive sequencing, despite a reduction in overall sequencing throughput due to strand-rejection processes. Interestingly, no significant differences in adaptive sampling enrichment efficiency were observed between the older and newer ONT sequencing chemistries, suggesting that adaptive sampling performs consistently across different library preparation kits.

CONCLUSION

Our findings underscore the importance of selecting a DNA extraction methodology that minimises bias to ensure an accurate representation of microbial diversity in metagenomic studies. Additionally, despite the improved base-calling accuracy provided by newer Q20 + sequencing chemistry, we demonstrate that even older ONT sequencing chemistries can achieve reliable metagenomic sequencing results, enabling researchers to confidently use these approaches depending on their specific experimental needs. Critically, we highlight the significant potential of ONT's adaptive sampling technology for targeted enrichment of specific genomes within metagenomic samples. This approach offers broad applicability for enriching target organisms or genetic elements (e.g., pathogens or plasmids) or depleting unwanted DNA (e.g., host DNA) in diverse sample types from environmental and clinical studies. However, researchers should carefully weigh the benefits of adaptive sampling against the potential trade-offs in sequencing throughput, particularly for low-abundance targets, where strand rejection can lead to pore blocking. These results provide valuable guidance for optimising adaptive sampling in metagenomic workflows to achieve specific research objectives.

Nanopore adaptive sampling to identify the NLR gene family in melon (Cucumis melo L.)

BACKGROUND

Nanopore adaptive sampling (NAS) offers a promising approach for assessing genetic diversity in targeted genomic regions. Here we designed and validated an experiment to enrich a set of resistance genes in several melon cultivars as a proof of concept.

RESULTS

Using the same reference to guide read acceptance or rejection with NAS, we successfully and accurately reconstructed the 15 regions in two newly assembled ssp. melo genomes and in a third ssp. agrestis cultivar. We obtained fourfold enrichment regardless of the tested samples, but with some variations according to the enriched regions. The accuracy of our assembly was further confirmed by PCR in the agrestis cultivar. We discussed parameters that could influence the enrichment and accuracy of NAS generated assemblies.

CONCLUSIONS

Overall, we demonstrated that NAS is a simple and efficient approach for exploring complex genomic regions, such as clusters of Nucleotide-binding site leucine-rich repeat (NLR) resistance genes. These regions are characterized by containing a high number of copy number variations, presence-absence polymorphisms and repetitive elements. These features make accurate assembly challenging but are crucial to study due to their central role in plant immunity and disease resistance. This approach facilitates resistance gene characterization in a large number of individuals, as required when breeding new cultivars suitable for the agroecological transition.

Nanopore sequencing of protozoa: Decoding biological information on a string of biochemical molecules into human-readable signals

Biological information is encoded in a sequence of biochemical molecules such as nucleic acids and amino acids, and nanopore sequencing is a long-read sequencing technology capable of directly decoding these molecules into human-readable signals. The long reads from nanopore sequencing offer the advantage of obtaining contiguous information, which is particularly beneficial for decoding complex or repetitive regions in a genome. In this study, we investigated the efficacy of nanopore sequencing in decoding biological information from distinctive genomes in metagenomic samples, which pose significant challenges for traditional short-read sequencing technologies. Specifically, we sequenced blood and fecal samples from mice infected with Trypanosoma brucei, a unicellular protozoan known for its hypervariable and dynamic regions that help it evade host immunity. Such characteristics are also prevalent in other host-dependent parasites, such as bacteriophages. The taxonomic classification results showed a high proportion of nanopore reads identified as T. brucei in the infected blood samples, with no significant identification in the control blood samples and fecal samples. Furthermore, metagenomic de novo assembly of these nanopore reads yielded contigs that mapped to the reference genome of T. brucei in the infected blood samples with over 96 % accuracy. This exploratory work demonstrates the potential of nanopore sequencing for the challenging task of classifying and assembling hypervariable and dynamic genomes from metagenomic samples.

Real-time and programmable transcriptome sequencing with PROFIT-seq

The high diversity and complexity of the eukaryotic transcriptome make it difficult to effectively detect specific transcripts of interest. Current targeted RNA sequencing methods often require complex pre-sequencing enrichment steps, which can compromise the comprehensive characterization of the entire transcriptome. Here we describe programmable full-length isoform transcriptome sequencing (PROFIT-seq), a method that enriches target transcripts while maintaining unbiased quantification of the whole transcriptome. PROFIT-seq employs combinatorial reverse transcription to capture polyadenylated, non-polyadenylated and circular RNAs, coupled with a programmable control system that selectively enriches target transcripts during sequencing. This approach achieves over 3-fold increase in effective data yield and reduces the time required for detecting specific pathogens or key mutations by 75%. We applied PROFIT-seq to study colorectal polyp development, revealing the intricate relationship between host immune responses and bacterial infection. PROFIT-seq offers a powerful tool for accurate and efficient sequencing of target transcripts while preserving overall transcriptome quantification, with broad applications in clinical diagnostics and targeted enrichment scenarios.

Advances in Host Depletion and Pathogen Enrichment Methods for Rapid Sequencing-Based Diagnosis of Bloodstream Infection

Bloodstream infection is a major cause of morbidity and death worldwide. Timely and appropriate treatment can reduce mortality among critically ill patients. Current diagnostic methods are too slow to inform precise antibiotic choice, leading to the prescription of empirical antibiotics, which may fail to cover the resistance profile of the pathogen, risking poor patient outcomes. Additionally, overuse of broad-spectrum antibiotics may lead to more resistant organisms, putting further pressure on the dwindling pipeline of antibiotics, and risk transmission of these resistant organisms in the health care environment. Therefore, rapid diagnostics are urgently required to better inform antibiotic choice early in the course of treatment. Sequencing offers great promise in reducing time to microbiological diagnosis; however, the amount of host DNA compared with the pathogen in patient samples presents a significant obstacle. Various host-depletion and bacterial-enrichment strategies have been used in samples, such as saliva, urine, or tissue. However, these methods have yet to be collectively integrated and/or extensively explored for rapid bloodstream infection diagnosis. Although most of these workflows possess individual strengths, their lack of analytical/clinical sensitivity and/or comprehensiveness demands additional improvements or synergistic application. This review provides a distinctive classification system for various methods based on their working principles to guide future research, and discusses their strengths and limitations and explores potential avenues for improvement to assist the reader in workflow selection.

RawHash2: mapping raw nanopore signals using hash-based seeding and adaptive quantization

SUMMARY

Raw nanopore signals can be analyzed while they are being generated, a process known as real-time analysis. Real-time analysis of raw signals is essential to utilize the unique features that nanopore sequencing provides, enabling the early stopping of the sequencing of a read or the entire sequencing run based on the analysis. The state-of-the-art mechanism, RawHash, offers the first hash-based efficient and accurate similarity identification between raw signals and a reference genome by quickly matching their hash values. In this work, we introduce RawHash2, which provides major improvements over RawHash, including more sensitive quantization and chaining algorithms, weighted mapping decisions, frequency filters to reduce ambiguous seed hits, minimizers for hash-based sketching, and support for the R10.4 flow cell version and POD5 and SLOW5 file formats. Compared to RawHash, RawHash2 provides better F1 accuracy (on average by 10.57% and up to 20.25%) and better throughput (on average by 4.0× and up to 9.9×) than RawHash.

AVAILABILITY AND IMPLEMENTATION

RawHash2 is available at https://github.com/CMU-SAFARI/RawHash. We also provide the scripts to fully reproduce our results on our GitHub page.

A case of an Angelman-syndrome caused by an intragenic duplication of UBE3A uncovered by adaptive nanopore sequencing

Adaptive nanopore sequencing as a diagnostic method for imprinting disorders and episignature analysis revealed an intragenic duplication of Exon 6 and 7 in UBE3A (NM_000462.5) in a patient with relatively mild Angelman-like syndrome. In an all-in-one nanopore sequencing analysis DNA hypomethylation of the SNURF:TSS-DMR, known contributing deletions on the maternal allele and point mutations in UBE3A could be ruled out as disease drivers. In contrast, breakpoints and orientation of the tandem duplication could clearly be defined. Segregation analysis in the family showed that the duplication derived de novo in the maternal grandfather. Our study shows the benefits of an all-in-one nanopore sequencing approach for the diagnostics of Angelman syndrome and other imprinting disorders.

ReadCurrent: a VDCNN-based tool for fast and accurate nanopore selective sequencing

Nanopore selective sequencing allows the targeted sequencing of DNA of interest using computational approaches rather than experimental methods such as targeted multiplex polymerase chain reaction or hybridization capture. Compared to sequence-alignment strategies, deep learning (DL) models for classifying target and nontarget DNA provide large speed advantages. However, the relatively low accuracy of these DL-based tools hinders their application in nanopore selective sequencing. Here, we present a DL-based tool named ReadCurrent for nanopore selective sequencing, which takes electric currents as inputs. ReadCurrent employs a modified very deep convolutional neural network (VDCNN) architecture, enabling significantly lower computational costs for training and quicker inference compared to conventional VDCNN. We evaluated the performance of ReadCurrent across 10 nanopore sequencing datasets spanning human, yeasts, bacteria, and viruses. We observed that ReadCurrent achieved a mean accuracy of 98.57% for classification, outperforming four other DL-based selective sequencing methods. In experimental validation that selectively sequenced microbial DNA from human DNA, ReadCurrent achieved an enrichment ratio of 2.85, which was higher than the 2.7 ratio achieved by MinKNOW using the sequence-alignment strategy. In summary, ReadCurrent can rapidly classify target and nontarget DNA with high accuracy, providing an alternative in the toolbox for nanopore selective sequencing. ReadCurrent is available at https://github.com/Ming-Ni-Group/ReadCurrent.

Fast and space-efficient taxonomic classification of long reads with hierarchical interleaved XOR filters

Metagenomic long-read sequencing is gaining popularity for various applications, including pathogen detection and microbiome studies. To analyze the large data created in those studies, software tools need to taxonomically classify the sequenced molecules and estimate the relative abundances of organisms in the sequenced sample. Because of the exponential growth of reference genome databases, the current taxonomic classification methods have large computational requirements. This issue motivated us to develop a new data structure for fast and memory-efficient querying of long reads. Here, we present Taxor as a new tool for long-read metagenomic classification using a hierarchical interleaved XOR filter data structure for indexing and querying large reference genome sets. Taxor implements several k-mer-based approaches, such as syncmers, for pseudoalignment to classify reads and an expectation-maximization algorithm for metagenomic profiling. Our results show that Taxor outperforms state-of-the-art tools regarding precision while having a similar recall for long-read taxonomic classification. Most notably, Taxor reduces the memory requirements and index size by >50% and is among the fastest tools regarding query times. This enables real-time metagenomics analysis with large reference databases on a small laptop in the field.

Biochemical-free enrichment or depletion of RNA classes in real-time during direct RNA sequencing with RISER

The heterogeneous composition of cellular transcriptomes poses a major challenge for detecting weakly expressed RNA classes, as they can be obscured by abundant RNAs. Although biochemical protocols can enrich or deplete specified RNAs, they are time-consuming, expensive and can compromise RNA integrity. Here we introduce RISER, a biochemical-free technology for the real-time enrichment or depletion of RNA classes. RISER performs selective rejection of molecules during direct RNA sequencing by identifying RNA classes directly from nanopore signals with deep learning and communicating with the sequencing hardware in real time. By targeting the dominant messenger and mitochondrial RNA classes for depletion, RISER reduces their respective read counts by more than 85%, resulting in an increase in sequencing depth of 47% on average for long non-coding RNAs. We also apply RISER for the depletion of globin mRNA in whole blood, achieving a decrease in globin reads by more than 90% as well as an increase in non-globin reads by 16% on average. Furthermore, using a GPU or a CPU, RISER is faster than GPU-accelerated basecalling and mapping. RISER's modular and retrainable software and intuitive command-line interface allow easy adaptation to other RNA classes. RISER is available at https://github.com/comprna/riser .

SimReadUntil for benchmarking selective sequencing algorithms on ONT devices

MOTIVATION

The Oxford Nanopore Technologies (ONT) ReadUntil API enables selective sequencing, which aims to selectively favor interesting over uninteresting reads, e.g. to deplete or enrich certain genomic regions. The performance gain depends on the selective sequencing decision-making algorithm (SSDA) which decides whether to reject a read, stop receiving a read, or wait for more data. Since real runs are time-consuming and costly, simulating the ONT sequencer with support for the ReadUntil API is highly beneficial for comparing and optimizing new SSDAs. Existing software like MinKNOW and UNCALLED only return raw signal data, are memory-intensive, require huge and often unavailable multi-fast5 files (≥100GB) and are not clearly documented.

RESULTS

We present the ONT device simulator SimReadUntil that takes a set of full reads as input, distributes them to channels and plays them back in real time including mux scans, channel gaps and blockages, and allows to reject reads as well as stop receiving data from them. Our modified ReadUntil API provides the basecalled reads rather than the raw signal, reducing computational load and focusing on the SSDA rather than on basecalling. Tuning the parameters of tools like ReadFish and ReadBouncer becomes easier because a GPU for basecalling is no longer required. We offer various methods to extract simulation parameters from a sequencing summary file and adapt ReadFish to replicate one of their enrichment experiments. SimReadUntil's gRPC interface allows standardized interaction with a wide range of programming languages.

AVAILABILITY AND IMPLEMENTATION

Code and fully worked examples are available on GitHub (https://github.com/ratschlab/sim_read_until).

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Nanopore adaptive sampling effectively enriches bacterial plasmids

Bacterial plasmids play a major role in the spread of antibiotic resistance genes. However, their characterization via DNA sequencing suffers from the low abundance of plasmid DNA in those samples. Although sample preparation methods can enrich the proportion of plasmid DNA before sequencing, these methods are expensive and laborious, and they might introduce a bias by enriching only for specific plasmid DNA sequences. Nanopore adaptive sampling could overcome these issues by rejecting uninteresting DNA molecules during the sequencing process. In this study, we assess the application of adaptive sampling for the enrichment of low-abundant plasmids in known bacterial isolates using two different adaptive sampling tools. We show that a significant enrichment can be achieved even on expired flow cells. By applying adaptive sampling, we also improve the quality of de novo plasmid assemblies and reduce the sequencing time. However, our experiments also highlight issues with adaptive sampling if target and non-target sequences span similar regions.

IMPORTANCE

Antimicrobial resistance causes millions of deaths every year. Mobile genetic elements like bacterial plasmids are key drivers for the dissemination of antimicrobial resistance genes. This makes the characterization of plasmids via DNA sequencing an important tool for clinical microbiologists. Since plasmids are often underrepresented in bacterial samples, plasmid sequencing can be challenging and laborious. To accelerate the sequencing process, we evaluate nanopore adaptive sampling as an in silico method for the enrichment of low-abundant plasmids. Our results show the potential of this cost-efficient method for future plasmid research but also indicate issues that arise from using reference sequences.

A signal processing and deep learning framework for methylation detection using Oxford Nanopore sequencing

Oxford Nanopore sequencing can detect DNA methylations from ionic current signal of single molecules, offering a unique advantage over conventional methods. Additionally, adaptive sampling, a software-controlled enrichment method for targeted sequencing, allows reduced representation methylation sequencing that can be applied to CpG islands or imprinted regions. Here we present DeepMod2, a comprehensive deep-learning framework for methylation detection using ionic current signal from Nanopore sequencing. DeepMod2 implements both a bidirectional long short-term memory (BiLSTM) model and a Transformer model and can analyze POD5 and FAST5 signal files generated on R9 and R10 flowcells. Additionally, DeepMod2 can run efficiently on central processing unit (CPU) through model pruning and can infer epihaplotypes or haplotype-specific methylation calls from phased reads. We use multiple publicly available and newly generated datasets to evaluate the performance of DeepMod2 under varying scenarios. DeepMod2 has comparable performance to Guppy and Dorado, which are the current state-of-the-art methods from Oxford Nanopore Technologies that remain closed-source. Moreover, we show a high correlation (r = 0.96) between reduced representation and whole-genome Nanopore sequencing. In summary, DeepMod2 is an open-source tool that enables fast and accurate DNA methylation detection from whole-genome or adaptive sequencing data on a diverse range of flowcell types.

RUBICON: a framework for designing efficient deep learning-based genomic basecallers

Nanopore sequencing generates noisy electrical signals that need to be converted into a standard string of DNA nucleotide bases using a computational step called basecalling. The performance of basecalling has critical implications for all later steps in genome analysis. Therefore, there is a need to reduce the computation and memory cost of basecalling while maintaining accuracy. We present RUBICON, a framework to develop efficient hardware-optimized basecallers. We demonstrate the effectiveness of RUBICON by developing RUBICALL, the first hardware-optimized mixed-precision basecaller that performs efficient basecalling, outperforming the state-of-the-art basecallers. We believe RUBICON offers a promising path to develop future hardware-optimized basecallers.

Direct RNA sequencing coupled with adaptive sampling enriches RNAs of interest in the transcriptome

Abundant cellular transcripts occupy most of the sequencing reads in the transcriptome, making it challenging to assay for low-abundant transcripts. Here, we utilize the adaptive sampling function of Oxford Nanopore sequencing to selectively deplete and enrich RNAs of interest without biochemical manipulation before sequencing. Adaptive sampling performed on a pool of in vitro transcribed RNAs resulted in a net increase of 22-30% in the proportion of transcripts of interest in the population. Enriching and depleting different proportions of the Candida albicans transcriptome also resulted in a 11-13.5% increase in the number of reads on target transcripts, with longer and more abundant transcripts being more efficiently depleted. Depleting all currently annotated Candida albicans transcripts did not result in an absolute enrichment of remaining transcripts, although we identified 26 previously unknown transcripts and isoforms, 17 of which are antisense to existing transcripts. Further improvements in the adaptive sampling of RNAs will allow the technology to be widely applied to study RNAs of interest in diverse transcriptomes.

NASCarD (Nanopore Adaptive Sampling with Carrier DNA): A Rapid, PCR-Free Method for SARS-CoV-2 Whole-Genome Sequencing in Clinical Samples

Whole-genome sequencing (WGS) represents the main technology for SARS-CoV-2 lineage characterization in diagnostic laboratories worldwide. The rapid, near-full-length sequencing of the viral genome is commonly enabled by high-throughput sequencing of PCR amplicons derived from cDNA molecules. Here, we present a new approach called NASCarD (Nanopore Adaptive Sampling with Carrier DNA), which allows a low amount of nucleic acids to be sequenced while selectively enriching for sequences of interest, hence limiting the production of non-target sequences. Using COVID-19 positive samples available during the omicron wave, we demonstrate how the method may lead to >99% genome completeness of the SARS-CoV-2 genome sequences within 7 h of sequencing at a competitive cost. The new approach may have applications beyond SARS-CoV-2 sequencing for other DNA or RNA pathogens in clinical samples.

NanoDeep: a deep learning framework for nanopore adaptive sampling on microbial sequencing

Nanopore sequencers can enrich or deplete the targeted DNA molecules in a library by reversing the voltage across individual nanopores. However, it requires substantial computational resources to achieve rapid operations in parallel at read-time sequencing. We present a deep learning framework, NanoDeep, to overcome these limitations by incorporating convolutional neural network and squeeze and excitation. We first showed that the raw squiggle derived from native DNA sequences determines the origin of microbial and human genomes. Then, we demonstrated that NanoDeep successfully classified bacterial reads from the pooled library with human sequence and showed enrichment for bacterial sequence compared with routine nanopore sequencing setting. Further, we showed that NanoDeep improves the sequencing efficiency and preserves the fidelity of bacterial genomes in the mock sample. In addition, NanoDeep performs well in the enrichment of metagenome sequences of gut samples, showing its potential applications in the enrichment of unknown microbiota. Our toolkit is available at https://github.com/lysovosyl/NanoDeep.

Targeted haplotyping in pharmacogenomics using Oxford Nanopore Technologies' adaptive sampling

Pharmacogenomics (PGx) studies the impact of interindividual genomic variation on drug response, allowing the opportunity to tailor the dosing regimen for each patient. Current targeted PGx testing platforms are mainly based on microarray, polymerase chain reaction, or short-read sequencing. Despite demonstrating great value for the identification of single nucleotide variants (SNVs) and insertion/deletions (INDELs), these assays do not permit identification of large structural variants, nor do they allow unambiguous haplotype phasing for star-allele assignment. Here, we used Oxford Nanopore Technologies' adaptive sampling to enrich a panel of 1,036 genes with well-documented PGx relevance extracted from the Pharmacogenomics Knowledge Base (PharmGKB). By evaluating concordance with existing truth sets, we demonstrate accurate variant and star-allele calling for five Genome in a Bottle reference samples. We show that up to three samples can be multiplexed on one PromethION flow cell without a significant drop in variant calling performance, resulting in 99.35% and 99.84% recall and precision for the targeted variants, respectively. This work advances the use of nanopore sequencing in clinical PGx settings.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Beyond assembly: the increasing flexibility of single-molecule sequencing technology

The maturation of high-throughput short-read sequencing technology over the past two decades has shaped the way genomes are studied. Recently, single-molecule, long-read sequencing has emerged as an essential tool in deciphering genome structure and function, including filling gaps in the human reference genome, measuring the epigenome and characterizing splicing variants in the transcriptome. With recent technological developments, these single-molecule technologies have moved beyond genome assembly and are being used in a variety of ways, including to selectively sequence specific loci with long reads, measure chromatin state and protein-DNA binding in order to investigate the dynamics of gene regulation, and rapidly determine copy number variation. These increasingly flexible uses of single-molecule technologies highlight a young and fast-moving part of the field that is leading to a more accessible era of nucleic acid sequencing.

RawHash: enabling fast and accurate real-time analysis of raw nanopore signals for large genomes

Nanopore sequencers generate electrical raw signals in real-time while sequencing long genomic strands. These raw signals can be analyzed as they are generated, providing an opportunity for real-time genome analysis. An important feature of nanopore sequencing, Read Until, can eject strands from sequencers without fully sequencing them, which provides opportunities to computationally reduce the sequencing time and cost. However, existing works utilizing Read Until either (i) require powerful computational resources that may not be available for portable sequencers or (ii) lack scalability for large genomes, rendering them inaccurate or ineffective. We propose RawHash, the first mechanism that can accurately and efficiently perform real-time analysis of nanopore raw signals for large genomes using a hash-based similarity search. To enable this, RawHash ensures the signals corresponding to the same DNA content lead to the same hash value, regardless of the slight variations in these signals. RawHash achieves an accurate hash-based similarity search via an effective quantization of the raw signals such that signals corresponding to the same DNA content have the same quantized value and, subsequently, the same hash value. We evaluate RawHash on three applications: (i) read mapping, (ii) relative abundance estimation, and (iii) contamination analysis. Our evaluations show that RawHash is the only tool that can provide high accuracy and high throughput for analyzing large genomes in real-time. When compared to the state-of-the-art techniques, UNCALLED and Sigmap, RawHash provides (i) 25.8× and 3.4× better average throughput and (ii) significantly better accuracy for large genomes, respectively. Source code is available at https://github.com/CMU-SAFARI/RawHash.

Rapid and precise genotyping of transgene zygosity in mice using an allele-specific method

Precise determination of transgene zygosity is essential for use of transgenic mice in research. Because integration loci of transgenes are usually unknown due to their random insertion, assessment of transgene zygosity remains a challenge. Current zygosity genotyping methods (progeny testing, qPCR, and NGS-computational biology analysis) are time consuming, prone to error or technically challenging. Here, we developed a novel method to determine transgene zygosity requiring no knowledge of transgene insertion loci. This method applies allele-specific restriction enzyme digestion of PCR products (RE/PCR) to rapidly and reliably quantify transgene zygosity. We demonstrate the applicability of this method to three transgenic strains of mice (Atm TgC3001L, Nes-Cre, and Syn1-Cre) harboring a unique restriction enzyme site on either the transgene or its homologous sequence in the mouse genome. This method is as accurate as the gold standard of progeny testing but requires 2 d instead of a month or more. It is also exceedingly more accurate than the most commonly used approach of qPCR quantification. Our novel method represents a significant technical advance in determining transgene zygosities in mice.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

Portable nanopore-sequencing technology: Trends in development and applications

Sequencing technology is the most commonly used technology in molecular biology research and an essential pillar for the development and applications of molecular biology. Since 1977, when the first generation of sequencing technology opened the door to interpreting the genetic code, sequencing technology has been developing for three generations. It has applications in all aspects of life and scientific research, such as disease diagnosis, drug target discovery, pathological research, species protection, and SARS-CoV-2 detection. However, the first- and second-generation sequencing technology relied on fluorescence detection systems and DNA polymerization enzyme systems, which increased the cost of sequencing technology and limited its scope of applications. The third-generation sequencing technology performs PCR-free and single-molecule sequencing, but it still depends on the fluorescence detection device. To break through these limitations, researchers have made arduous efforts to develop a new advanced portable sequencing technology represented by nanopore sequencing. Nanopore technology has the advantages of small size and convenient portability, independent of biochemical reagents, and direct reading using physical methods. This paper reviews the research and development process of nanopore sequencing technology (NST) from the laboratory to commercially viable tools; discusses the main types of nanopore sequencing technologies and their various applications in solving a wide range of real-world problems. In addition, the paper collates the analysis tools necessary for performing different processing tasks in nanopore sequencing. Finally, we highlight the challenges of NST and its future research and application directions.

Efficient real-time selective genome sequencing on resource-constrained devices

BACKGROUND

Third-generation nanopore sequencers offer selective sequencing or "Read Until" that allows genomic reads to be analyzed in real time and abandoned halfway if not belonging to a genomic region of "interest." This selective sequencing opens the door to important applications such as rapid and low-cost genetic tests. The latency in analyzing should be as low as possible for selective sequencing to be effective so that unnecessary reads can be rejected as early as possible. However, existing methods that employ a subsequence dynamic time warping (sDTW) algorithm for this problem are too computationally intensive that a massive workstation with dozens of CPU cores still struggles to keep up with the data rate of a mobile phone-sized MinION sequencer.

RESULTS

In this article, we present Hardware Accelerated Read Until (HARU), a resource-efficient hardware-software codesign-based method that exploits a low-cost and portable heterogeneous multiprocessor system-on-chip platform with on-chip field-programmable gate arrays (FPGA) to accelerate the sDTW-based Read Until algorithm. Experimental results show that HARU on a Xilinx FPGA embedded with a 4-core ARM processor is around 2.5× faster than a highly optimized multithreaded software version (around 85× faster than the existing unoptimized multithreaded software) running on a sophisticated server with a 36-core Intel Xeon processor for a SARS-CoV-2 dataset. The energy consumption of HARU is 2 orders of magnitudes lower than the same application executing on the 36-core server.

CONCLUSIONS

HARU demonstrates that nanopore selective sequencing is possible on resource-constrained devices through rigorous hardware-software optimizations. The source code for the HARU sDTW module is available as open source at https://github.com/beebdev/HARU, and an example application that uses HARU is at https://github.com/beebdev/sigfish-haru.