Abstract:Abstract Background High throughput RNA sequencing is a powerful approach to study gene expression. Due to the complex multiple-steps protocols in data acquisition, extreme deviation of a sample from samples of the same treatment group may occur due to technical variation or true biological differences. The high-dimensionality of the data with few biological replicates make it challenging to accurately detect those samples, and this issue is not well studied in the literature currently. Robust statistics is a family of theories and techniques aim to detect the outliers by first fitting the majority of the data and then flagging data points that deviate from it. Robust statistics have been widely used in multivariate data analysis for outlier detection in chemometrics and engineering. Here we apply robust statistics on RNA-seq data analysis. Results We report the use of two robust principal component analysis (rPCA) methods, PcaHubert and PcaGrid , to detect outlier samples in multiple simulated and real biological RNA-seq data sets with positive control outlier samples. PcaGrid achieved 100% sensitivity and 100% specificity in all the tests using positive control outliers with varying degrees of divergence. We applied rPCA methods and classical principal component analysis (cPCA) on an RNA-Seq data set profiling gene expression of the external granule layer in the cerebellum of control and conditional SnoN knockout mice. Both rPCA methods detected the same two outlier samples but cPCA failed to detect any. We performed differentially expressed gene detection before and after outlier removal as well as with and without batch effect modeling. We validated gene expression changes using quantitative reverse transcription PCR and used the result as reference to compare the performance of eight different data analysis strategies. Removing outliers without batch effect modeling performed the best in term of detecting biologically relevant differentially expressed genes. Conclusions rPCA implemented in the PcaGrid function is an accurate and objective method to detect outlier samples. It is well suited for high-dimensional data with small sample sizes like RNA-seq data. Outlier removal can significantly improve the performance of differential gene detection and downstream functional analysis.

Identifying Relevant Covariates in RNA-seq Analysis by Pseudo-Variable Augmentation

Variance component score test for time-course gene set analysis of longitudinal RNA-seq data

Robust principal component analysis for accurate outlier sample detection in RNA-Seq data

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

AutoGeneS: Automatic Gene Selection Using Multi-Objective Optimization for RNA-seq Deconvolution

Machine learning-assisted identification of factors contributing to the technical variability between bulk and single-cell RNA-seq experiments

RNAseqCovarImpute: a multiple imputation procedure that outperforms complete case and single imputation differential expression analysis

Group Heteroscedasticity - A Silent Saboteur of Power and False Discovery in RNA-Seq Differential Expression

Rank-based Bayesian variable selection for genome-wide transcriptomic analyses

V-SVA: an R Shiny application for detecting and annotating hidden sources of variation in single-cell RNA-seq data

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

Bayesian Nonparametric Variable Selection as an Exploratory Tool for Finding Genes that Matter

Comparison of Parametric and Machine Methods for Variable Selection in Simulated Genetic Analysis Workshop 19 Data.

Online FDR Control for RNAseq Data

Accounting for technical noise in differential expression analysis of single-cell RNA sequencing data

A censored-Poisson model based approach to the analysis of RNA-seq data

Variable selection method for the identification of epistatic models.

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

The impact of package selection and versioning on single-cell RNA-seq analysis

Controlled Variable Selection from Summary Statistics Only? A Solution via GhostKnockoffs and Penalized Regression