Probabilistic PCA of censored data: accounting for uncertainties in the visualization of high-throughput single-cell qPCR data

MotivationHigh-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.ResultsWe 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 implementationThe 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.

1. 1
   Diffusion maps for high-dimensional single-cell analysis of differentiation data

   Bioinformatics2015709 citationsDataRank 0.985
2. 2
   Revealing the vectors of cellular identity with single-cell genomics

   Nature Biotechnology2016701 citationsDataRank 0.983
3. 3
   Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments

   Bioinformatics2012441 citationsDataRank 0.914
4. 4
   Characterization of transcriptional networks in blood stem and progenitor cells using high-throughput single-cell gene expression analysis

   Nature Cell Biology2013279 citationsDataRank 0.845
5. 5
   A novel approach for resolving differences in single-cell gene expression patterns from zygote to blastocyst

   Bioinformatics201267 citationsDataRank 2.6