{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,12]],"date-time":"2026-05-12T10:13:48Z","timestamp":1778580828914,"version":"3.51.4"},"reference-count":28,"publisher":"Oxford University Press (OUP)","issue":"13","license":[{"start":{"date-parts":[[2016,12,28]],"date-time":"2016-12-28T00:00:00Z","timestamp":1482883200000},"content-version":"vor","delay-in-days":1024,"URL":"http:\/\/creativecommons.org\/licenses\/by\/3.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2014,7,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>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.<\/jats:p>\n               <jats:p>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.<\/jats:p>\n               <jats:p>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 .<\/jats:p>\n               <jats:p>Contact: \u00a0fbuettner.phys@gmail.com<\/jats:p>\n               <jats:p>Supplementary information: Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btu134","type":"journal-article","created":{"date-parts":[[2014,3,12]],"date-time":"2014-03-12T00:16:49Z","timestamp":1394583409000},"page":"1867-1875","source":"Crossref","is-referenced-by-count":21,"title":["Probabilistic PCA of censored data: accounting for uncertainties in the visualization of high-throughput single-cell qPCR data"],"prefix":"10.1093","volume":"30","author":[{"given":"Florian","family":"Buettner","sequence":"first","affiliation":[{"name":"1 Institute of Computational Biology, Helmholtz-Zentrum M\u00fcnchen, 85764 Neuherberg, Germany, 2 Department of Haematology, University of Cambridge, Cambridge Institute for Medical Research and Wellcome Trust & MRC Cambridge Stem Cell Institute, Cambridge CB2 0XY, UK and 3 Department of Mathematics, TU M\u00fcnchen, 85748 Garching, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Victoria","family":"Moignard","sequence":"additional","affiliation":[{"name":"1 Institute of Computational Biology, Helmholtz-Zentrum M\u00fcnchen, 85764 Neuherberg, Germany, 2 Department of Haematology, University of Cambridge, Cambridge Institute for Medical Research and Wellcome Trust & MRC Cambridge Stem Cell Institute, Cambridge CB2 0XY, UK and 3 Department of Mathematics, TU M\u00fcnchen, 85748 Garching, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Berthold","family":"G\u00f6ttgens","sequence":"additional","affiliation":[{"name":"1 Institute of Computational Biology, Helmholtz-Zentrum M\u00fcnchen, 85764 Neuherberg, Germany, 2 Department of Haematology, University of Cambridge, Cambridge Institute for Medical Research and Wellcome Trust & MRC Cambridge Stem Cell Institute, Cambridge CB2 0XY, UK and 3 Department of Mathematics, TU M\u00fcnchen, 85748 Garching, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Fabian J.","family":"Theis","sequence":"additional","affiliation":[{"name":"1 Institute of Computational Biology, Helmholtz-Zentrum M\u00fcnchen, 85764 Neuherberg, Germany, 2 Department of Haematology, University of Cambridge, Cambridge Institute for Medical Research and Wellcome Trust & MRC Cambridge Stem Cell Institute, Cambridge CB2 0XY, UK and 3 Department of Mathematics, TU M\u00fcnchen, 85748 Garching, Germany"},{"name":"1 Institute of Computational Biology, Helmholtz-Zentrum M\u00fcnchen, 85764 Neuherberg, Germany, 2 Department of Haematology, University of Cambridge, Cambridge Institute for Medical Research and Wellcome Trust & MRC Cambridge Stem Cell Institute, Cambridge CB2 0XY, UK and 3 Department of Mathematics, TU M\u00fcnchen, 85748 Garching, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2014,3,10]]},"reference":[{"key":"2023012711111797000_btu134-B1","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1038\/nbt.2594","article-title":"viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia","volume":"31","author":"Amir","year":"2013","journal-title":"Nat. Biotechnol."},{"key":"2023012711111797000_btu134-B2","doi-asserted-by":"crossref","first-page":"R106","DOI":"10.1186\/gb-2010-11-10-r106","article-title":"Differential expression analysis for sequence count data","volume":"11","author":"Anders","year":"2010","journal-title":"Genome Biol."},{"key":"2023012711111797000_btu134-B3","doi-asserted-by":"crossref","first-page":"e46423","DOI":"10.1371\/journal.pone.0046423","article-title":"Novel statistical approaches for non-normal censored immunological data: analysis of cytokine and gene expression data","volume":"7","author":"Ballenberger","year":"2012","journal-title":"PLoS One"},{"key":"2023012711111797000_btu134-B4","volume-title":"Pattern Recognition and Machine Learning (Information Science and Statistics)","author":"Bishop","year":"2006"},{"key":"2023012711111797000_btu134-B5","doi-asserted-by":"crossref","first-page":"1093","DOI":"10.1038\/nmeth.2645","article-title":"Accounting for technical noise in single-cell RNA-seq experiments","volume":"10","author":"Brennecke","year":"2013","journal-title":"Nat. Methods"},{"key":"2023012711111797000_btu134-B6","doi-asserted-by":"crossref","first-page":"i626","DOI":"10.1093\/bioinformatics\/bts385","article-title":"A novel approach for resolving differences in single-cell gene expression patterns from zygote to blastocyst","volume":"28","author":"Buettner","year":"2012","journal-title":"Bioinformatics"},{"key":"2023012711111797000_btu134-B7","doi-asserted-by":"crossref","first-page":"1120","DOI":"10.1038\/nbt.2038","article-title":"Single-cell dissection of transcriptional heterogeneity in human colon tumors","volume":"29","author":"Dalerba","year":"2011","journal-title":"Nat. Biotechnol."},{"key":"2023012711111797000_btu134-B8","doi-asserted-by":"crossref","first-page":"133","DOI":"10.1016\/j.jim.2013.03.002","article-title":"Highly multiplexed quantitation of gene expression on single cells","volume":"391","author":"Dominguez","year":"2013","journal-title":"J. Immunol. Methods"},{"key":"2023012711111797000_btu134-B9","article-title":"Application guidance: single-cell data analysis","author":"Fluidigm Corporation","year":"2012"},{"key":"2023012711111797000_btu134-B10","first-page":"522","article-title":"Style-based inverse kinematics","volume-title":"ACM Transactions on Graphics (TOG)","author":"Grochow","year":"2004"},{"key":"2023012711111797000_btu134-B11","doi-asserted-by":"crossref","first-page":"675","DOI":"10.1016\/j.devcel.2010.02.012","article-title":"Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst","volume":"18","author":"Guo","year":"2010","journal-title":"Dev. Cell"},{"key":"2023012711111797000_btu134-B12","first-page":"1783","article-title":"Probabilistic non-linear principal component analysis with gaussian process latent variable models","volume":"6","author":"Lawrence","year":"2005","journal-title":"J. Mach. Learn. Res."},{"key":"2023012711111797000_btu134-B13","first-page":"56","article-title":"Extensions of the informative vector machine","volume-title":"Deterministic and Statistical Methods in Machine Learning, Volume 3635 of Lecture Notes in Computer Science","author":"Lawrence","year":"2005"},{"key":"2023012711111797000_btu134-B14","article-title":"Gaussian process latent variable models for visualisation of high dimensional data","volume-title":"NIPS","author":"Lawrence","year":"2004"},{"key":"2023012711111797000_btu134-B15","doi-asserted-by":"crossref","first-page":"513","DOI":"10.1145\/1143844.1143909","article-title":"Local distance preservation in the gp-lvm through back constraints","volume-title":"ICML","author":"Lawrence","year":"2006"},{"key":"2023012711111797000_btu134-B16","doi-asserted-by":"crossref","first-page":"1691","DOI":"10.1289\/ehp.7199","article-title":"Epidemiologic evaluation of measurement data in the presence of detection limits","volume":"112","author":"Lubin","year":"2004","journal-title":"Environ. Health Perspect."},{"key":"2023012711111797000_btu134-B17","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1093\/bioinformatics\/bts714","article-title":"Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments","volume":"29","author":"McDavid","year":"2013","journal-title":"Bioinformatics"},{"key":"2023012711111797000_btu134-B18","first-page":"362","article-title":"Expectation propagation for approximate bayesian inference","volume-title":"Proceedings of the Seventeenth Conference Annual Conference on Uncertainty in Artificial Intelligence (UAI-01)","author":"Minka","year":"2001"},{"key":"2023012711111797000_btu134-B19","doi-asserted-by":"crossref","first-page":"363","DOI":"10.1038\/ncb2709","article-title":"Characterization of transcriptional networks in blood stem and progenitor cells using high-throughput single-cell gene expression analysis","volume":"15","author":"Moignard","year":"2013","journal-title":"Nat. Cell Biol."},{"key":"2023012711111797000_btu134-B20","article-title":"Netlab: algorithms for pattern recognition","volume-title":"Advances in Pattern Recognition","author":"Nabney","year":"2001"},{"key":"2023012711111797000_btu134-B21","doi-asserted-by":"crossref","first-page":"631","DOI":"10.1016\/j.cell.2008.01.025","article-title":"Hematopoiesis: an evolving paradigm for stem cell biology","volume":"132","author":"Orkin","year":"2008","journal-title":"Cell"},{"key":"2023012711111797000_btu134-B22","doi-asserted-by":"crossref","first-page":"287","DOI":"10.1038\/ncb2442","article-title":"Inferring rules of lineage commitment in haematopoiesis","volume":"14","author":"Pina","year":"2012","journal-title":"Nat. Cell Biol."},{"key":"2023012711111797000_btu134-B23","article-title":"Gaussian processes for machine learning","volume-title":"Adaptive Computation and Machine Learning","author":"Rasmussen","year":"2006"},{"key":"2023012711111797000_btu134-B24","doi-asserted-by":"crossref","first-page":"503","DOI":"10.1038\/nmeth.1338","article-title":"Quantitative analysis of gene expression in a single cell by qPCR","volume":"6","author":"Taniguchi","year":"2009","journal-title":"Nat. Methods."},{"key":"2023012711111797000_btu134-B25","doi-asserted-by":"crossref","first-page":"2501","DOI":"10.1093\/molbev\/msr077","article-title":"Complex principal component and correlation structure of 16 yeast genomic variables","volume":"28","author":"Theis","year":"2011","journal-title":"Mol. Biol. Evol."},{"key":"2023012711111797000_btu134-B26","doi-asserted-by":"crossref","first-page":"59","DOI":"10.1186\/1471-2172-9-59","article-title":"Evaluation of regression methods when immunological measurements are constrained by detection limits","volume":"9","author":"Uh","year":"2008","journal-title":"BMC Immunol."},{"key":"2023012711111797000_btu134-B27","first-page":"2579","article-title":"Visualizing data using t-SNE","volume":"9","author":"van der Maaten","year":"2008","journal-title":"J. Mach. Learn. Res."},{"key":"2023012711111797000_btu134-B28","doi-asserted-by":"crossref","first-page":"1131","DOI":"10.1038\/nsmb.2660","article-title":"Single-cell rna-seq profiling of human preimplantation embryos and embryonic stem cells","volume":"20","author":"Yan","year":"2013","journal-title":"Nat. Struct. Mol. Biol."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/30\/13\/1867\/48923654\/bioinformatics_30_13_1867.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/30\/13\/1867\/48923654\/bioinformatics_30_13_1867.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,1,27]],"date-time":"2023-01-27T11:51:21Z","timestamp":1674820281000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/30\/13\/1867\/2422238"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2014,3,10]]},"references-count":28,"journal-issue":{"issue":"13","published-print":{"date-parts":[[2014,7,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btu134","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2014,7,1]]},"published":{"date-parts":[[2014,3,10]]}}}