Software for pre-processing Illumina next-generation sequencing short read sequences

Abstract

BackgroundWhen compared to Sanger sequencing technology, next-generation sequencing (NGS) technologies are hindered by shorter sequence read length, higher base-call error rate, non-uniform coverage, and platform-specific sequencing artifacts. These characteristics lower the quality of their downstream analyses, e.g. de novo and reference-based assembly, by introducing sequencing artifacts and errors that may contribute to incorrect interpretation of data. Although many tools have been developed for quality control and pre-processing of NGS data, none of them provide flexible and comprehensive trimming options in conjunction with parallel processing to expedite pre-processing of large NGS datasets.MethodsWe developed ngsShoRT (next-generation sequencing Short Reads Trimmer), a flexible and comprehensive open-source software package written in Perl that provides a set of algorithms commonly used for pre-processing NGS short read sequences. We compared the features and performance of ngsShoRT with existing tools: CutAdapt, NGS QC Toolkit and Trimmomatic. We also compared the effects of using pre-processed short read sequences generated by different algorithms on de novo and reference-based assembly for three different genomes: Caenorhabditis elegans, Saccharomyces cerevisiae S288c, and Escherichia coli O157 H7.ResultsSeveral combinations of ngsShoRT algorithms were tested on publicly available Illumina GA II, HiSeq 2000, and MiSeq eukaryotic and bacteria genomic short read sequences with the focus on removing sequencing artifacts and low-quality reads and/or bases. Our results show that across three organisms and three sequencing platforms, trimming improved the mean quality scores of trimmed sequences. Using trimmed sequences for de novo and reference-based assembly improved assembly quality as well as assembler performance. In general, ngsShoRT outperformed comparable trimming tools in terms of trimming speed and improvement of de novo and reference-based assembly as measured by assembly contiguity and correctness.ConclusionsTrimming of short read sequences can improve the quality of de novo and reference-based assembly and assembler performance. The parallel processing capability of ngsShoRT reduces trimming time and improves the memory efficiency when dealing with large datasets. We recommend combining sequencing artifacts removal, and quality score based read filtering and base trimming as the most consistent method for improving sequence quality and downstream assemblies.ngsShoRT source code, user guide and tutorial are available at http://research.bioinformatics.udel.edu/genomics/ngsShoRT/. ngsShoRT can be incorporated as a pre-processing step in genome and transcriptome assembly projects.

Spike-in Genomic DNA for Validating Performance of Metagenomics Workflows

“…Sequence bases with quality score below 30 were discarded from both 3 and 5 end using cutadapt version 1.13 [14]. Remaining reads were further filtered by BWA-MEM version 0.7.13-r1126, paired-end mode [15] and samtools version 1.3.1 with htslib 1.3.1 [16] to remove reads mapped to human genome (GRCh38) and Phix genome.…”

The 3'-5' exoribonuclease Dis3 regulates the expression of specific microRNAs inDrosophilawing imaginal discs

“…Adaptors were removed using Cutadapt. 46 The remaining reads were mapped to miRBase release19 using Bowtie 2. 47 Mapped reads were normalized per million reads (RPKM).…”

Integrated tRNA, transcript, and protein profiles in response to steroid hormone signaling

“…Small RNA library preparation and sequencing was performed by the Genomic Services Laboratory at HudsonAlpha. Small RNAsequencing reads and adaptor sequences were trimmed using the ngsShort toolkit (Chen et al 2014). Sequencing reads were subsequently mapped to the Dm6 release of the Drosophila genome (Hoskins et al 2015) using bowtie1.0.0 (Langmead 2010), allowing one mismatch and randomly assigning multimapping reads a single alignment to the highest scoring region.…”