The impact of read length on quantification of differentially expressed genes and splice junction detection

Abstract

BackgroundThe initial next-generation sequencing technologies produced reads of 25 or 36 bp, and only from a single-end of the library sequence. Currently, it is possible to reliably produce 300 bp paired-end sequences for RNA expression analysis. While read lengths have consistently increased, people have assumed that longer reads are more informative and that paired-end reads produce better results than single-end reads. We used paired-end 101 bp reads and trimmed them to simulate different read lengths, and also separated the pairs to produce single-end reads. For each read length and paired status, we evaluated differential expression levels between two standard samples and compared the results to those obtained by qPCR.ResultsWe found that, with the exception of 25 bp reads, there is little difference for the detection of differential expression regardless of the read length. Once single-end reads are at a length of 50 bp, the results do not change substantially for any level up to, and including, 100 bp paired-end. However, splice junction detection significantly improves as the read length increases with 100 bp paired-end showing the best performance. We performed the same analysis on two ENCODE samples and found consistent results confirming that our conclusions have broad application.ConclusionsA researcher could save substantial resources by using 50 bp single-end reads for differential expression analysis instead of using longer reads. However, splicing detection is unquestionably improved by paired-end and longer reads. Therefore, an appropriate read length should be used based on the final goal of the study.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-015-0697-y) contains supplementary material, which is available to authorized users.

Hypothalamic transcriptomic alterations in male and female California mice (Peromyscus californicus) developmentally exposed to bisphenol A or ethinyl estradiol

“…And it has been demonstrated that using single‐end (SE) reads for transcriptome assembly and quantification has relatively minor effects on gene‐level DE estimates, as described previously (Gonzalez and Joly ; Chhangawala et al. ). Thus, we chose to use SE data for this work.…”

Transcriptome Analysis of Piperlongumine-Treated Human Pancreatic Cancer Cells Reveals Involvement of Oxidative Stress and Endoplasmic Reticulum Stress Pathways

“…Six barcoded libraries were created, which were pooled together and sequenced in one full lane of the flow cell on an Illumina HiSeq 2500. Single-end reads of 50 base pair (bp), which are sufficient to reliably identify mRNA expression differences relative to paired-end longer reads, 13 were obtained. TopHat 14 was used to map reads to the reference genome (Build 38 and Ensembl 78) for each sample (max mismatches 15 and maximum insertion and deletion lengths 18).…”

Dual RNA-Sequencing to Elucidate the Plant-Pathogen Duel

“…If a high quality and well annotated reference sequence is available, increasing read length above 50 bp is unnecessary for accurate detection of differential expression (Chhangawala et al, 2015). Similarly, sequencing of paired-end instead of single-end reads does not significantly affect detection of differential expression in these cases (Chhangawala et al, 2015).…”

“…If a high quality and well annotated reference sequence is available, increasing read length above 50 bp is unnecessary for accurate detection of differential expression (Chhangawala et al, 2015). Similarly, sequencing of paired-end instead of single-end reads does not significantly affect detection of differential expression in these cases (Chhangawala et al, 2015). Conversely, when studying organisms with less well defined reference sequences, sequencing of longer paired-end reads increases the accuracy of splice junction detection (Chhangawala et al, 2015).…”

“…Similarly, sequencing of paired-end instead of single-end reads does not significantly affect detection of differential expression in these cases (Chhangawala et al, 2015). Conversely, when studying organisms with less well defined reference sequences, sequencing of longer paired-end reads increases the accuracy of splice junction detection (Chhangawala et al, 2015). When no reference sequence is available, it is often assumed that longer reads equate to increased accuracy for de novo assembly.…”