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.

Robustified MANOVA with Applications in Detecting Differentially Expressed Genes from Oligonucleotide Arrays

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

An Empirical Bayes Method for Robust Variance Estimation in Detecting DEGs Using Microarray Data.

Shrinkage-based regularization tests for high-dimensional data with application to gene set analysis

Power enhancement via multivariate outlier testing with gene expression arrays

Model-based Analysis of Oligonucleotide Arrays: Expression Index Computation and Outlier Detection.

Identifying differentially expressed genes in human acute leukemia and mouse brain microarray datasets utilizing QTModel

Novel Statistical Framework to Identify Differentially Expressed Genes Allowing Transcriptomic Background Differences

Detecting False Expression Signals in High-Density Oligonucleotide Arrays by an in Silico Approach.

Identification of differentially expressed genes with multivariate outlier analysis.

A Tractable Probabilistic Model for Affymetrix Probe-Level Analysis Across Multiple Chips.

Hotelling's T2 multivariate profiling for detecting differential expression in microarrays

Robust detection and genotyping of single feature polymorphisms from gene expression data.

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

A new framework for identifying differentially expressed genes

MAGE: Monte Carlo method for Aberrant Gene Expression

An imputation approach for oligonucleotide microarrays

A Statistical Framework for the Design of Microarray Experiments and Effective Detection of Differential Gene Expression

Robustly detecting differential expression in RNA sequencing data using observation weights

Identifying Differentially Expressed Genes Based on Probe Level Data for GeneChip Arrays.

A robust RUV-testing procedure via gamma-divergence