Abstract:BackgroundRNA-sequencing (RNA-Seq) has become a powerful technology to characterize gene expression profiles because it is more accurate and comprehensive than microarrays. Although statistical methods that have been developed for microarray data can be applied to RNA-Seq data, they are not ideal due to the discrete nature of RNA-Seq data. The Poisson distribution and negative binomial distribution are commonly used to model count data. Recently, Witten (Annals Appl Stat 5:2493–2518, 2011) proposed a Poisson linear discriminant analysis for RNA-Seq data. The Poisson assumption may not be as appropriate as the negative binomial distribution when biological replicates are available and in the presence of overdispersion (i.e., when the variance is larger than or equal to the mean). However, it is more complicated to model negative binomial variables because they involve a dispersion parameter that needs to be estimated.ResultsIn this paper, we propose a negative binomial linear discriminant analysis for RNA-Seq data. By Bayes’ rule, we construct the classifier by fitting a negative binomial model, and propose some plug-in rules to estimate the unknown parameters in the classifier. The relationship between the negative binomial classifier and the Poisson classifier is explored, with a numerical investigation of the impact of dispersion on the discriminant score. Simulation results show the superiority of our proposed method. We also analyze two real RNA-Seq data sets to demonstrate the advantages of our method in real-world applications.ConclusionsWe have developed a new classifier using the negative binomial model for RNA-seq data classification. Our simulation results show that our proposed classifier has a better performance than existing works. The proposed classifier can serve as an effective tool for classifying RNA-seq data. Based on the comparison results, we have provided some guidelines for scientists to decide which method should be used in the discriminant analysis of RNA-Seq data. R code is available at http://www.comp.hkbu.edu.hk/~xwan/NBLDA.Ror https://github.com/yangchadam/NBLDA

PLNseq: a multivariate Poisson lognormal distribution for high-throughput matched RNA-sequencing read count data.

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

NBLDA: negative binomial linear discriminant analysis for RNA-Seq data

Gene regulatory network in single cells based on the Poisson log-normal model

PennSeq: Accurate Isoform-Specific Gene Expression Quantification in RNA-Seq by Modeling Non-Uniform Read Distribution

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

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

Estimating graphical models for count data with applications to single-cell gene network

Robust estimation of isoform expression with RNA-Seq data

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

Sparse Inverse Covariance Estimation for High-throughput microRNA Sequencing Data in the Poisson Log-Normal Graphical Model

Rseqnp: a Non-Parametric Approach for Detecting Differential Expression and Splicing from RNA-Seq Data

A sparse negative binomial mixture model for clustering RNA-seq count data

PDEGEM: Modeling non-uniform read distribution in RNA-Seq data

Nonparametric clustering of RNA-sequencing data

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

Joint Between-Sample Normalization and Differential Expression Detection Through ℓ0-Regularized Regression

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

Finite mixtures of matrix-variate Poisson-log normal distributions for three-way count data

clrDV: a differential variability test for RNA-Seq data based on the skew-normal distribution

Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model