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.

Distortion-free PCA on sample space for highly variable gene detection from single-cell RNA-seq data

Probabilistic PCA of Censored Data: Accounting for Uncertainties in the Visualization of High-Throughput Single-Cell Qpcr Data.

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

CHLPCA: Correntropy-Based Hypergraph Regularized Sparse PCA for Single-Cell Type Identification.

VPAC: Variational Projection for Accurate Clustering of Single-Cell Transcriptomic Data

Preprocessing of Single Cell RNA Sequencing Data Using Correlated Clustering and Projection

Accurate Identification of Single-Cell Types Via Correntropy-Based Sparse PCA Combining Hypergraph and Fusion Similarity

Joint L2,p-norm and Random Walk Graph Constrained PCA for Single-Cell RNA-seq Data.

Joint Modeling of Cellular Heterogeneity and Condition Effects with scPCA in Single-Cell RNA-Seq

K-nearest-neighbors induced topological PCA for single cell RNA-sequence data analysis

A Partitioning PCA Approach for Extracting Cell-Cycle Regulated Genes

Dual Graph regularized PCA based on Different Norm Constraints for Bi-clustering Analysis on Single-cell RNA-seq Data

Robust Principal Component Analysis Based on Hypergraph Regularization for Sample Clustering and Co-Characteristic Gene Selection

Probabilistic Count Matrix Factorization for Single Cell Expression Data Analysis

Identification of Differentially Expressed Genes Using Disjoint Principal Component Analysis Coupled with Genetic Algorithm

Deterministic column subset selection for single-cell RNA-Seq

PCA via joint graph Laplacian and sparse constraint: Identification of differentially expressed genes and sample clustering on gene expression data

scGIR: deciphering cellular heterogeneity via gene ranking in single-cell weighted gene correlation networks

K-Nearest-Neighbors Induced Topological PCA for scRNA Sequence Data Analysis

Feature selection and dimension reduction for single-cell RNA-Seq based on a multinomial model

Partially-independent component analysis for tissue heterogeneity correction in microarray gene expression analysis