SERE: Single-parameter quality control and sample comparison for RNA-Seq

Abstract

BackgroundAssessing the reliability of experimental replicates (or global alterations corresponding to different experimental conditions) is a critical step in analyzing RNA-Seq data. Pearson’s correlation coefficient r has been widely used in the RNA-Seq field even though its statistical characteristics may be poorly suited to the task.ResultsHere we present a single-parameter test procedure for count data, the Simple Error Ratio Estimate (SERE), that can determine whether two RNA-Seq libraries are faithful replicates or globally different. Benchmarking shows that the interpretation of SERE is unambiguous regardless of the total read count or the range of expression differences among bins (exons or genes), a score of 1 indicating faithful replication (i.e., samples are affected only by Poisson variation of individual counts), a score of 0 indicating data duplication, and scores >1 corresponding to true global differences between RNA-Seq libraries. On the contrary the interpretation of Pearson’s r is generally ambiguous and highly dependent on sequencing depth and the range of expression levels inherent to the sample (difference between lowest and highest bin count). Cohen’s simple Kappa results are also ambiguous and are highly dependent on the choice of bins. For quantifying global sample differences SERE performs similarly to a measure based on the negative binomial distribution yet is simpler to compute.ConclusionsSERE can therefore serve as a straightforward and reliable statistical procedure for the global assessment of pairs or large groups of RNA-Seq datasets by a single statistical parameter.

Functional and RNA-Sequencing Analysis Revealed Expression of a Novel Stay-Green Gene from Zoysia japonica (ZjSGR) Caused Chlorophyll Degradation and Accelerated Senescence in Arabidopsis

“…Pearson's correlation coefficient was employed to test the reliability of RNA-seq data (Schulze et al, 2012). Correlation coefficients (approximately close to 1) suggested a strong correlation among the biological replicates in WT and ZjSGR -overexpressing plants (Figure S3).…”

OsPRR37confers an expanded regulation of the diurnal rhythms of the transcriptome and photoperiodic flowering pathways in rice

“…The normalized counts per million (CPM) was calculated according to Cuffdiff in the Cufflinks software (Jiang & Wong, ; Trapnell et al, ). The Pearson's correlation coefficient was used to evaluate the correlation between any biological repeats at the same time point (Schulze, Kanwar, Golzenleuchter, Therneau, & Beutler, ). To eliminate the nonexpressed genes from the rice transcriptome, the genes that showed CPM > 0 in at least three samples in either GL or OsPRR37 ‐OE5 were considered to comprise the entire transcriptome in the present study.…”

AsHSP17, a creeping bentgrass small heat shock protein modulates plant photosynthesis and ABA‐dependent and independent signalling to attenuate plant response to abiotic stress

“…Two biological replicates of the wild type and AsHSP17 transgenic plants were used for cDNA library construction and sequencing. The reliability of RNA‐seq data was confirmed by Pearson's correlation coefficient (Schulze et al ) between each pair of the biological replicates and further visualized in the correlation plot (Fig. S4).…”