Abstract:BACKGROUND:The advent of the NGS technologies has permitted profiling of whole-genome transcriptomes (i.e., RNA-Seq) at unprecedented speed and very low cost. RNA-Seq provides a far more precise measurement of transcript levels and their isoforms compared to other methods such as microarrays. A fundamental goal of RNA-Seq is to better identify expression changes between different biological or disease conditions. However, existing methods for detecting differential expression from RNA-Seq count data have not been comprehensively evaluated in large-scale RNA-Seq datasets. Many of them suffer from inflation of type I error and failure in controlling false discovery rate especially in the presence of abnormal high sequence read counts in RNA-Seq experiments.RESULTS:To address these challenges, we propose a powerful and robust tool, termed deGPS, for detecting differential expression in RNA-Seq data. This framework contains new normalization methods based on generalized Poisson distribution modeling sequence count data, followed by permutation-based differential expression tests. We systematically evaluated our new tool in simulated datasets from several large-scale TCGA RNA-Seq projects, unbiased benchmark data from compcodeR package, and real RNA-Seq data from the development transcriptome of Drosophila. deGPS can precisely control type I error and false discovery rate for the detection of differential expression and is robust in the presence of abnormal high sequence read counts in RNA-Seq experiments.CONCLUSIONS:Software implementing our deGPS was released within an R package with parallel computations ( https://github.com/LL-LAB-MCW/deGPS ). deGPS is a powerful and robust tool for data normalization and detecting different expression in RNA-Seq experiments. Beyond RNA-Seq, deGPS has the potential to significantly enhance future data analysis efforts from many other high-throughput platforms such as ChIP-Seq, MBD-Seq and RIP-Seq.

Multi-scale Poisson process approaches for differential expression analysis of high-throughput sequencing data

Differential Expression Analysis for RNAseq Using Poisson Mixed Models.

A flexible model for correlated count data, with application to multi-condition differential expression analyses of single-cell RNA sequencing data

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

An Effective Differential Expression Analysis of Deep-Sequencing Data Based on the Poisson Log-Normal Model.

Differential expression analysis for paired RNA-seq data

Degps is a Powerful Tool for Detecting Differential Expression in RNA-sequencing Studies

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

Using Poisson mixed-effects model to quantify transcript-level gene expression in RNA-Seq.

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

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

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

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

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

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

Robustly detecting differential expression in RNA sequencing data using observation weights

Count-based differential expression analysis of RNA sequencing data using R and Bioconductor

Comparison of differential accessibility analysis strategies for ATAC-seq data

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

Predicting the Number of Bases to Attain Sufficient Coverage in High-Throughput Sequencing Experiments.

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