bioRxivBioinformaticsPreprint
Evaluating Nanopore read trimming and quality filtering in bacterial genomics
Achmad Dimas Cahyaning FurqoniD, Ryan R. WickiD, Leah Wendy RobertsiD, Michael B. HalliD
- Citations
- 0
Abstract
Oxford Nanopore Technologies (ONT) sequencing is widely used in bacterial genomics, particularly for reference-based variant calling and de novo assembly. Both analyses routinely start with quality control (QC) preprocessing, but as ONT basecalling accuracy has improved, its necessity and downstream trade-offs remain poorly defined. We benchmarked 17 preprocessing combinations, built from three adapter trimmers (Dorado, Porechop_ABI, Barbell) and three quality-based processing tools (Filtlong, Chopper, SeqKit) with untrimmed and unprocessed controls, plus Fastplong, which performs both. Each combination was evaluated across 13 bacterial isolates at 100x and 20x depth under the HAC and SUP basecalling models. Variant calling was assessed by F1 score. Assemblies were scored on sequence accuracy, contiguity, replicon recovery and residual adapter contamination, then aggregated into a weighted arithmetic mean using weights from a survey of the microbial genomics community. With SUP basecalling, raw reads with no trimming or quality-based processing produced accurate SNP calls and outperformed several pipelines that applied both operations; adapter trimming alone left SNP accuracy unchanged. Subsetting reads by combined length and quality with Filtlong gave the highest SNP F1 scores at both depths, except with Barbell at 20x, where it ranked second lowest. As is often the case in benchmarking studies, assembly results involved trade-offs, and the best tool depends on which aspects of performance matter most for a given use. Filtlong produced the lowest error rates, but every Filtlong pipeline lost three small plasmids at both depths. Barbell trimming raised error rates and reduced contiguity at 20x, while skipping adapter trimming left residual adapter and barcode sequences in assemblies. Under the community-derived weights, Porechop_ABI followed by Chopper had the highest overall score at both depths, without winning any single metric. QC tool choice should therefore be guided by the downstream analysis goal rather than applied indiscriminately across all workflows.
- organism: bacteria (other)
- assay: whole-genome sequencing
- research focus: variant calling
- approach: benchmark study
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Read trimming has minimal effect on bacterial SNP calling accuracy
Stephen J. Bush·Microbial GenomicsPeer-reviewed·Dec 1, 2020
- Prowler: A novel trimming algorithm for Oxford Nanopore sequence data
Simon Lee, Loan T Nguyen, … Elizabeth M. Ross·BioinformaticsPeer-reviewed·Sep 2, 2021
- How low can you go? Short-read polishing of Oxford Nanopore bacterial genome assemblies
George Spyro Bouras, Louise M. Judd, … Ryan R. Wick·Microbial GenomicsPeer-reviewed·Jun 4, 2024
- Performance of neural network basecalling tools for Oxford Nanopore sequencing
Ryan R. Wick, Louise M. Judd, Kathryn E. Holt·Genome biologyPeer-reviewed·Jun 24, 2019
- Comparison of R9.4.1/Kit10 and R10/Kit12 Oxford Nanopore flowcells and chemistries in bacterial genome reconstruction
Nicholas D. Sanderson, Natalia Kapel, … Nicole E Stoesser·Microbial GenomicsPeer-reviewed·Jan 10, 2023
- Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing
N. D. Sanderson, K. Hopkins, … N. Stoesser·Microbial GenomicsPeer-reviewed·May 7, 2024