Abstract:A flexible statistical framework is developed for the analysis of read counts from RNA-Seq gene expression studies. It provides the ability to analyse complex experiments involving multiple treatment conditions and blocking variables while still taking full account of biological variation. Biological variation between RNA samples is estimated separately from the technical variation associated with sequencing technologies. Novel empirical Bayes methods allow each gene to have its own specific variability, even when there are relatively few biological replicates from which to estimate such variability. The pipeline is implemented in the edgeR package of the Bioconductor project. A case study analysis of carcinoma data demonstrates the ability of generalized linear model methods (GLMs) to detect differential expression in a paired design, and even to detect tumour-specific expression changes. The case study demonstrates the need to allow for gene-specific variability, rather than assuming a common dispersion across genes or a fixed relationship between abundance and variability. Genewise dispersions de-prioritize genes with inconsistent results and allow the main analysis to focus on changes that are consistent between biological replicates. Parallel computational approaches are developed to make non-linear model fitting faster and more reliable, making the application of GLMs to genomic data more convenient and practical. Simulations demonstrate the ability of adjusted profile likelihood estimators to return accurate estimators of biological variability in complex situations. When variation is gene-specific, empirical Bayes estimators provide an advantageous compromise between the extremes of assuming common dispersion or separate genewise dispersion. The methods developed here can also be applied to count data arising from DNA-Seq applications, including ChIP-Seq for epigenetic marks and DNA methylation analyses.

EBSeq-HMM: a Bayesian Approach for Identifying Gene-Expression Changes in Ordered RNA-seq Experiments

A Bayesian model selection approach for identifying differentially expressed transcripts from RNA-Seq data

Estimation of Isoform Expression in Rna-Seq Data Using A Hierarchical Bayesian Model

A Robust Empirical Bayesian Method For Detecting Differentially Expressed Genes

Differential expression analysis for paired RNA-seq data

Identifying differentially expressed transcripts from RNA-seq data with biological variation

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

A two-step strategy for detecting differential gene expression in cDNA microarray data

Impulse Model-Based Differential Expression Analysis of Time Course Sequencing Data

Rseqdiff: Detecting Differential Isoform Expression from RNA-Seq Data Using Hierarchical Likelihood Ratio Test.

Identifying Differentially Spliced Genes from Two Groups of Rna-Seq Samples

An Empirical Bayes Change-Point Model for Identifying 3′ and 5′ Alternative Splicing by Next-Generation RNA Sequencing

A Closed-Form Solution to the 2-Sample Problem for Quantifying Changes in Gene Expression using Bayes Factors

Bayesian Analysis of RNA-Seq Data Using a Family of Negative Binomial Models.

BADER: Bayesian analysis of differential expression in RNA sequencing data

HBA-DEALS: accurate and simultaneous identification of differential expression and splicing using hierarchical Bayesian analysis

Non-parametric Bayesian modelling of digital gene expression data

DEIsoM: a Hierarchical Bayesian Model for Identifying Differentially Expressed Isoforms Using Biological Replicates

Group Heteroscedasticity - A Silent Saboteur of Power and False Discovery in RNA-Seq Differential Expression