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.

Differential expression analysis for paired RNA-seq data

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

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

Differential Expression Analysis for RNAseq Using Poisson Mixed Models.

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

A Two-Part Mixed Model for Differential Expression Analysis in Single-Cell High-Throughput Gene Expression Data.

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

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

A comparison of methods for differential expression analysis of RNA-seq data

Unit-Free and Robust Detection of Differential Expression from RNA-Seq Data

Differential expression analysis on RNA-seq count data based on penalized matrix decomposition

A Unified Model for Joint Normalization and Differential Gene Expression Detection in RNA-Seq Data.

Detecting Differential Expression from RNA-seq Data with Expression Measurement Uncertainty

A Unified Model for Differential Expression Analysis of RNA-seq Data Via L1-Penalized Linear Regression

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

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

A Method Based on Differential Entropy-Like Function for Detecting Differentially Expressed Genes Across Multiple Conditions in RNA-Seq Studies

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

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

Network Module-Based Model in the Differential Expression Analysis for RNA-seq.

Three Differential Expression Analysis Methods for RNA Sequencing: limma, EdgeR, DESeq2