Statistical models for RNA-seq data derived from a two-condition 48-replicate experiment

Marek Gierliński, Christian Cole, Pietà Schofield, Nicholas J Schurch, Alexander Sherstnev, Vijender Singh, Nicola Wrobel, Karim Gharbi, Gordon Simpson, Tom Owen-Hughes, Mark Blaxter, Geoffrey J Barton

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

MOTIVATION: High-throughput RNA sequencing (RNA-seq) is now the standard method to determine differential gene expression. Identifying differentially expressed genes crucially depends on estimates of read count variability. These estimates are typically based on statistical models such as the negative binomial distribution, which is employed by the tools edgeR, DESeq and cuffdiff. Until now, the validity of these models has usually been tested on either low-replicate RNA-seq data or simulations.

RESULTS: A 48-replicate RNA-seq experiment in yeast was performed and data tested against theoretical models. The observed gene read counts were consistent with both log-normal and negative binomial distributions, while the mean-variance relation followed the line of constant dispersion parameter of ~0.01. The high-replicate data also allowed for strict quality control and screening of "bad" replicates, which can drastically affect the gene read-count distribution.

AVAILABILITY: RNA-seq data have been submitted to ENA archive with project ID PRJEB5348.

CONTACT: [email protected].

Original languageEnglish
Publication statusPublished - 23 Jul 2015


Dive into the research topics of 'Statistical models for RNA-seq data derived from a two-condition 48-replicate experiment'. Together they form a unique fingerprint.

Cite this