Abstract:MOTIVATION:High-throughput single-cell quantitative real-time polymerase chain reaction (qPCR) is a promising technique allowing for new insights in complex cellular processes. However, the PCR reaction can be detected only up to a certain detection limit, whereas failed reactions could be due to low or absent expression, and the true expression level is unknown. Because this censoring can occur for high proportions of the data, it is one of the main challenges when dealing with single-cell qPCR data. Principal component analysis (PCA) is an important tool for visualizing the structure of high-dimensional data as well as for identifying subpopulations of cells. However, to date it is not clear how to perform a PCA of censored data. We present a probabilistic approach that accounts for the censoring and evaluate it for two typical datasets containing single-cell qPCR data.RESULTS:We use the Gaussian process latent variable model framework to account for censoring by introducing an appropriate noise model and allowing a different kernel for each dimension. We evaluate this new approach for two typical qPCR datasets (of mouse embryonic stem cells and blood stem/progenitor cells, respectively) by performing linear and non-linear probabilistic PCA. Taking the censoring into account results in a 2D representation of the data, which better reflects its known structure: in both datasets, our new approach results in a better separation of known cell types and is able to reveal subpopulations in one dataset that could not be resolved using standard PCA.AVAILABILITY AND IMPLEMENTATION:The implementation was based on the existing Gaussian process latent variable model toolbox (https://github.com/SheffieldML/GPmat); extensions for noise models and kernels accounting for censoring are available at http://icb.helmholtz-muenchen.de/censgplvm.

HePPCAT: Probabilistic PCA for Data with Heteroscedastic Noise

HeMPPCAT: Mixtures of Probabilistic Principal Component Analysers for Data with Heteroscedastic Noise

ALPCAH: Sample-wise Heteroscedastic PCA with Tail Singular Value Regularization

Towards a Theoretical Analysis of PCA for Heteroscedastic Data

Streaming Probabilistic PCA for Missing Data with Heteroscedastic Noise

Robust Principal Component Analysis Based on Maximum Correntropy Criterion

Sparse Probabilistic Principal Component Analysis Model for Plant-Wide Process Monitoring

Self-Paced Probabilistic Principal Component Analysis for Data with Outliers

On the Estimation Performance of Generalized Power Method for Heteroscedastic Probabilistic PCA

Inference for Heteroskedastic PCA with Missing Data

Deflated HeteroPCA: Overcoming the curse of ill-conditioning in heteroskedastic PCA

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

High-Dimensional Principal Component Analysis with Heterogeneous Missingness

1 A Probability Extension of PCA to Detect and Diagnose Sensor Faults in Air Handling Units

Generalized probabilistic principal component analysis of correlated data

Ensemble Principal Component Analysis

Improved Algorithms for High-Dimensional Robust Pca

Self-paced Principal Component Analysis

Robust PCA for High Dimensional Data based on Characteristic Transformation

Poisson PCA: Poisson Measurement Error corrected PCA, with Application to Microbiome Data

Functional PCA With Covariate-Dependent Mean and Covariance Structure