Abstract:Advances in small RNAs (sRNAs)-related studies have posed a challenge for NGS-related bioinformatics, especially regarding the correct mapping of sRNAs. Depending on the algorithms and scoring matrices on which they are based, aligners are influenced by the characteristics of the dataset and the reference genome. These influences have been studied mainly in eukaryotes and to some extent in prokaryotes. However, in bacteria, the selection of aligners depending on sRNA-seq data associated with outer membrane vesicles (OMVs) and the features of the corresponding bacterial reference genome has not yet been investigated. We selected five aligners: BBmap, Bowtie2, BWA, Minimap2 and Segemehl, known for their generally good performance, to test them in mapping OMV-associated sRNAs from Aliivibrio fischeri to the bacterial reference genome. Significant differences in the performance of the five aligners were observed, resulting in differential recognition of OMV-associated sRNA biotypes in A. fischeri . Our results suggest that aligner(s) should not be arbitrarily selected for this task, which is often done, as this can be detrimental to the biological interpretation of NGS analysis results. Since each aligner has specific advantages and disadvantages, these need to be considered depending on the characteristics of the input OMV sRNAs dataset and the corresponding bacterial reference genome to improve the detection of existing, biologically important OMV sRNAs. Until we learn more about these dependencies, we recommend using at least two, preferably three, aligners that have good metrics for the given dataset/bacterial reference genome. The overlapping results should be considered trustworthy, yet their differences should not be dismissed lightly, but treated carefully in order not to overlook any biologically important OMV sRNA. This can be achieved by applying the intersect-then-combine approach. For the mapping of OMV-associated sRNAs of A. fischeri to the reference genome organized into two circular chromosomes and one circular plasmid, containing copies of sequences with rRNA- and tRNA-related features and no copies of sequences with protein-encoding features, if the aligners are used with their default parameters, we advise avoiding Segemehl, and recommend using the intersect-then-combine approach with BBmap, BWA and Minimap2 to improve the potential for discovery of biologically important OMV-associated sRNAs.

benchNGS : An approach to benchmark short reads alignment tools

NCBench: providing an open, reproducible, transparent, adaptable, and continuous benchmark approach for DNA-sequencing-based variant calling

Systematic benchmark of state-of-the-art variant calling pipelines identifies major factors affecting accuracy of coding sequence variant discovery

A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing

Benchmarking RNA-Seq Aligners at Base-Level and Junction Base-Level Resolution Using the Arabidopsis thaliana Genome

SequenceLab: A Comprehensive Benchmark of Computational Methods for Comparing Genomic Sequences

Who Watches the Watchmen? An Appraisal of Benchmarks for Multiple Sequence Alignment

CGAP-align: a high performance DNA short read alignment tool.

An efficient Burrows-Wheeler transform-based aligner for short read mapping

Benchmarking genome assembly methods on metagenomic sequencing data

Fast and accurate short read alignment with hybrid hash-tree data structure

New Developments Of Alignment-Free Sequence Comparison: Measures, Statistics And Next-Generation Sequencing

Benchmarking of bioinformatics tools for the hybrid de novo assembly of human whole-genome sequencing data

A Re-Sequencing Tool For High Mismatch-Tolerant Short Read Alignment Based On Burrows-Wheeler Transform

Benchmarking long-read aligners and SV callers for structural variation detection in Oxford nanopore sequencing data

NGSCheckMate: software for validating sample identity in next-generation sequencing studies within and across data types

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Benchmarking Recent Computational Tools for DNA-binding Protein Identification

Benchmarking of five NGS mapping tools for the reference alignment of bacterial outer membrane vesicles-associated small RNAs

Benchmarking short and long read polishing tools for nanopore assemblies: achieving near-perfect genomes for outbreak isolates

Benchmarking DNA Methylation Analysis of 14 Alignment Algorithms for Whole Genome Bisulfite Sequencing in Mammals