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.

Defining pervasive transcription units using chromatin RNA-sequencing data.

Whole-genome Epigenetic Function Annotation Through Sgrna Libraries Synthesized by Controlled Template-Dependent Elongation

Protocol for Quantitative Analysis of RNA 3′-End Processing Induced by Disassociated Subunits Using Chromatin-Associated RNA-seq Data

The USTC co-opts an ancient machinery to drive piRNA transcription in C. elegans

Transient transcriptome sequencing: computational pipeline to quantify genome-wide RNA kinetic parameters and transcriptional enhancer activity

Genome wide full-length transcript analysis using 5' and 3' paired-end-tag next generation sequencing (RNA-PET).

BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution

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

Enhancing Transcriptome Mapping with Rapid PRO-seq Profiling of Nascent RNA

Using Population-Scale Transcriptomic and Genomic Data to Map 3′ UTR Alternative Polyadenylation Quantitative Trait Loci

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

SCAPTURE: a deep learning-embedded pipeline that captures polyadenylation information from 3′ tag-based RNA-seq of single cells

The apparent loss of PRC2 chromatin occupancy as an artefact of RNA depletion

Functional classification of noncoding RNAs associated with distinct histone modifications by PIRCh-seq

Quantitative sequencing using BID-seq uncovers abundant pseudouridines in mammalian mRNA at base resolution

Measuring Nascent Transcripts by Nascent-seq

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

TransRef Enables Accurate Transcriptome Assembly by Redefining Accurate Neo-Splicing Graphs.

Quantitative Profiling of Pseudouridylation Landscape in the Human Transcriptome

A Streamlined and Generalized Analysis of Chromatin ImmunoPrecipitation Paired-End diTag Data