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.

Application of the fractional stable distributions for approximation of gene expression profiles

An Efficient and Assumption-Free Method to Approximate Protein Level Distribution in the Two-States Gene Expression Model.

Analytical distributions for stochastic gene expression

Estimation distribution of algorithm for fuzzy clustering gene expression data

Cancer Classification Using Entropy Analysis in Fractional Fourier Domain of Gene Expression Profile

Estimation of Distribution Algorithm Using Variety of Information

Probabilistic estimation of microarray data reliability and underlying gene expression

An Extended Fractional Kalman Filter for Inferring Gene Regulatory Networks Using Time-Series Data

Exact distributions for stochastic models of gene expression with arbitrary regulation

A Fokker-Planck approach to the study of robustness in gene expression

Fuzzy Entropy Clustering by Searching Local Border Points for the Analysis of Gene Expression Data

A Generalized Expression Of Multi-Target Probability Density In The Framework Of Fisst

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

Differential analysis in Transcriptomic: The strength of randomly picking 'reference' genes

Exact protein distributions for stochastic models of gene expression using partitioning of Poisson processes

Approximation of State Variables for Discrete-Time Stochastic Genetic Regulatory Networks with Leakage, Distributed, and Probabilistic Measurement Delays: a Robust Stability Problem.

On discrete approximations of stable distributions

Moment-Based Approximations of Probability Mass Functions with Applications Involving Order Statistics

An Uncertainty Quantification Approach to the Study of Gene Expression Robustness.

Sufficient Direction Factor Model and Its Application to Gene Expression Quantitative Trait Loci Discovery

A Bayesian Determination of Threshold for Identifying Differentially Expressed Genes in Microarray Experiments