Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads

Abstract

BackgroundAdapter trimming is a prerequisite step for analyzing next-generation sequencing (NGS) data when the reads are longer than the target DNA/RNA fragments. Although typically used in small RNA sequencing, adapter trimming is also used widely in other applications, such as genome DNA sequencing and transcriptome RNA/cDNA sequencing, where fragments shorter than a read are sometimes obtained because of the limitations of NGS protocols. For the newly emerged Nextera long mate-pair (LMP) protocol, junction adapters are located in the middle of all properly constructed fragments; hence, adapter trimming is essential to gain the correct paired reads. However, our investigations have shown that few adapter trimming tools meet both efficiency and accuracy requirements simultaneously. The performances of these tools can be even worse for paired-end and/or mate-pair sequencing.ResultsTo improve the efficiency of adapter trimming, we devised a novel algorithm, the bit-masked k-difference matching algorithm, which has O(kn) expected time with O(m) space, where k is the maximum number of differences allowed, n is the read length, and m is the adapter length. This algorithm makes it possible to fully enumerate all candidates that meet a specified threshold, e.g. error ratio, within a short period of time. To improve the accuracy of this algorithm, we designed a simple and easy-to-explain statistical scoring scheme to evaluate candidates in the pattern matching step. We also devised scoring schemes to fully exploit the paired-end/mate-pair information when it is applicable. All these features have been implemented in an industry-standard tool named Skewer (https://sourceforge.net/projects/skewer). Experiments on simulated data, real data of small RNA sequencing, paired-end RNA sequencing, and Nextera LMP sequencing showed that Skewer outperforms all other similar tools that have the same utility. Further, Skewer is considerably faster than other tools that have comparative accuracies; namely, one times faster for single-end sequencing, more than 12 times faster for paired-end sequencing, and 49% faster for LMP sequencing.ConclusionsSkewer achieved as yet unmatched accuracies for adapter trimming with low time bound.

An Evolutionarily Conserved Transcriptional Activator-Repressor Module Controls Expression of Genes for D-Galacturonic Acid Utilization inAspergillus niger

“…We refer to A. niger gene IDs based on the most up-to-date and accurate annotation of the A. niger NRRL3 genome (http://genome.fungalgenomics.ca/). The RNA-seq reads were cleaned by correcting sequencing errors with Rcorrector (Song and Florea 2015), trimming sequencing adapters and low quality sequences with Skewer (Jiang et al 2014), and removing ribosomal RNA with SortMeRNA (Kopylova et al 2012). The cleaned reads were mapped to NRRL3 transcripts and counted with Salmon (Patro et al 2016), and the read counts were analyzed for differences in transcript expression between genotypes with DESeq2 (Love et al 2014).…”

Plants Release Precursors of Histone Deacetylase Inhibitors to Suppress Growth of Competitors

“…We trimmed reads for adapter sequences using the software skewer (v0.1.120) (Jiang et al, 2014). Reads were then quality-filtered using SHORE (v0.9.0) (Ossowski et al, 2008), and reads trimmed to <30 bp due to low quality were discarded.…”

Venus flytrap carnivorous lifestyle builds on herbivore defense strategies

“…Passed data sets were quality trimmed using skewer (version 0.1.67; -Q 30; -q 30; -l 75; -m pe) (Jiang et al 2014) The transcriptome was assembled using Trinity (release 2013-02-16; -jaccard_clip; -min_kmer_cov 2; -path_reinforcement_distance 75) (Grabherr et al 2011). The assembly was screened for artificial fusion events caused by low-complexity regions or highly similar UTRs.…”