{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,6]],"date-time":"2026-07-06T22:33:40Z","timestamp":1783377220776,"version":"3.54.6"},"reference-count":48,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2023,4,11]],"date-time":"2023-04-11T00:00:00Z","timestamp":1681171200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100008658","name":"German Cancer Research Center","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100008658","id-type":"DOI","asserted-by":"publisher"}]},{"name":"German Federal Ministry of Education and Research"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,5,4]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Factor analysis is a widely used tool for unsupervised dimensionality reduction of high-throughput datasets in molecular biology, with recently proposed extensions designed specifically for spatial transcriptomics data. However, these methods expect (count) matrices as data input and are therefore not directly applicable to single molecule resolution data, which are in the form of coordinate lists annotated with genes and provide insight into subcellular spatial expression patterns. To address this, we here propose FISHFactor, a probabilistic factor model that combines the benefits of spatial, non-negative factor analysis with a Poisson point process likelihood to explicitly model and account for the nature of single molecule resolution data. In addition, FISHFactor shares information across a potentially large number of cells in a common weight matrix, allowing consistent interpretation of factors across cells and yielding improved latent variable estimates.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We compare FISHFactor to existing methods that rely on aggregating information through spatial binning and cannot combine information from multiple cells and show that our method leads to more accurate results on simulated data. We show that our method is scalable and can be readily applied to large datasets. Finally, we demonstrate on a real dataset that FISHFactor is able to identify major subcellular expression patterns and spatial gene clusters in a data-driven manner.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>The model implementation, data simulation and experiment scripts are available under https:\/\/www.github.com\/bioFAM\/FISHFactor.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad183","type":"journal-article","created":{"date-parts":[[2023,4,11]],"date-time":"2023-04-11T11:06:58Z","timestamp":1681211218000},"source":"Crossref","is-referenced-by-count":18,"title":["FISHFactor: a probabilistic factor model for spatial transcriptomics data with subcellular resolution"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-7865-1749","authenticated-orcid":false,"given":"Florin C","family":"Walter","sequence":"first","affiliation":[{"name":"Division of Computational Genomics and System Genetics, German Cancer Research Center (DKFZ) , Heidelberg 69120, Germany"},{"name":"Genome Biology Unit, European Molecular Biology Laboratory (EMBL) , Heidelberg 69117, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8818-7193","authenticated-orcid":false,"given":"Oliver","family":"Stegle","sequence":"additional","affiliation":[{"name":"Division of Computational Genomics and System Genetics, German Cancer Research Center (DKFZ) , Heidelberg 69120, Germany"},{"name":"Genome Biology Unit, European Molecular Biology Laboratory (EMBL) , Heidelberg 69117, Germany"},{"name":"Cellular Genetics Programme, Wellcome Sanger Institute, Wellcome Genome Campus , Hinxton, Cambridge CB10 1SA, United Kingdom"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8397-3515","authenticated-orcid":false,"given":"Britta","family":"Velten","sequence":"additional","affiliation":[{"name":"Division of Computational Genomics and System Genetics, German Cancer Research Center (DKFZ) , Heidelberg 69120, Germany"},{"name":"Cellular Genetics Programme, Wellcome Sanger Institute, Wellcome Genome Campus , Hinxton, Cambridge CB10 1SA, United Kingdom"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2023,4,11]]},"reference":[{"key":"2023051312281836000_btad183-B1","doi-asserted-by":"crossref","first-page":"e8124","DOI":"10.15252\/msb.20178124","article-title":"Multi-Omics factor analysis\u2014a framework for unsupervised integration of multi-omics data sets","volume":"14","author":"Argelaguet","year":"2018","journal-title":"Mol Syst Biol"},{"key":"2023051312281836000_btad183-B2","doi-asserted-by":"crossref","first-page":"111","DOI":"10.1186\/s13059-020-02015-1","article-title":"MOFA+: a statistical framework for comprehensive integration of multi-modal single-cell data","volume":"21","author":"Argelaguet","year":"2020","journal-title":"Genome Biol"},{"key":"2023051312281836000_btad183-B3","doi-asserted-by":"crossref","first-page":"2419","DOI":"10.1038\/s41467-018-04724-5","article-title":"Spatial maps of prostate cancer transcriptomes reveal an unexplored landscape of heterogeneity","volume":"9","author":"Berglund","year":"2018","journal-title":"Nat Commun"},{"key":"2023051312281836000_btad183-B4","first-page":"1","article-title":"Pyro: deep universal probabilistic programming","volume":"20","author":"Bingham","year":"2019","journal-title":"J Mach Learn Res"},{"key":"2023051312281836000_btad183-B5","doi-asserted-by":"crossref","first-page":"4164","DOI":"10.1073\/pnas.0308531101","article-title":"Metagenes and molecular pattern discovery using matrix factorization","volume":"101","author":"Brunet","year":"2004","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023051312281836000_btad183-B6","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1038\/nrm3918","article-title":"In the right place at the right time: visualizing and understanding mRNA localization","volume":"16","author":"Buxbaum","year":"2015","journal-title":"Nat Rev Mol Cell Biol"},{"key":"2023051312281836000_btad183-B7","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1207\/s15327906mbr0102_10","article-title":"The scree test for the number of factors","volume":"1","author":"Cattell","year":"1966","journal-title":"Multivariate Behav Res"},{"key":"2023051312281836000_btad183-B8","doi-asserted-by":"crossref","first-page":"aaa6090","DOI":"10.1126\/science.aaa6090","article-title":"Spatially resolved, highly multiplexed RNA profiling in single cells","volume":"348","author":"Chen","year":"2015","journal-title":"Science"},{"key":"2023051312281836000_btad183-B9","doi-asserted-by":"crossref","first-page":"e22","DOI":"10.1093\/nar\/gkx1206","article-title":"Efficient in situ barcode sequencing using padlock probe-based BaristaSeq","volume":"46","author":"Chen","year":"2018","journal-title":"Nucleic Acids Res"},{"key":"2023051312281836000_btad183-B10","doi-asserted-by":"crossref","first-page":"932","DOI":"10.1038\/s41592-018-0175-z","article-title":"Spatial organization of the somatosensory cortex revealed by osmFISH","volume":"15","author":"Codeluppi","year":"2018","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B11","volume-title":"Advances in Neural Information Processing Systems","author":"Dugas","year":"2000"},{"key":"2023051312281836000_btad183-B12","first-page":"10445","article-title":"Temporal alignment and latent Gaussian process factor inference in population spike trains","volume":"31","author":"Duncker","year":"2018","journal-title":"Advances in Neural Information Processing Systems"},{"key":"2023051312281836000_btad183-B13","doi-asserted-by":"crossref","first-page":"1153","DOI":"10.1038\/nmeth.4500","article-title":"Profiling the transcriptome with RNA SPOTs","volume":"14","author":"Eng","year":"2017","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B14","doi-asserted-by":"crossref","first-page":"235","DOI":"10.1038\/s41586-019-1049-y","article-title":"Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH+","volume":"568","author":"Eng","year":"2019","journal-title":"Nature"},{"key":"2023051312281836000_btad183-B15","doi-asserted-by":"crossref","first-page":"585","DOI":"10.1126\/science.280.5363.585","article-title":"Visualization of single RNA transcripts in situ","volume":"280","author":"Femino","year":"1998","journal-title":"Science"},{"key":"2023051312281836000_btad183-B16","first-page":"7576","article-title":"GPyTorch: blackbox matrix-matrix Gaussian process inference with GPU acceleration","volume":"31","author":"Gardner","year":"2018","journal-title":"Advances in Neural Information Processing Systems"},{"key":"2023051312281836000_btad183-B17","first-page":"351","author":"Hensman","year":"2015"},{"key":"2023051312281836000_btad183-B18","first-page":"1303","article-title":"Stochastic variational inference","volume":"14","author":"Hoffman","year":"2013","journal-title":"J Mach Learn Res"},{"key":"2023051312281836000_btad183-B19","doi-asserted-by":"crossref","first-page":"857","DOI":"10.1038\/nmeth.2563","article-title":"In situ sequencing for RNA analysis in preserved tissue and cells","volume":"10","author":"Ke","year":"2013","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B20","author":"Kingma","year":"2014"},{"key":"2023051312281836000_btad183-B21","doi-asserted-by":"crossref","first-page":"788","DOI":"10.1038\/44565","article-title":"Learning the parts of objects by non-negative matrix factorization","volume":"401","author":"Lee","year":"1999","journal-title":"Nature"},{"key":"2023051312281836000_btad183-B22","doi-asserted-by":"crossref","first-page":"442","DOI":"10.1038\/nprot.2014.191","article-title":"Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues","volume":"10","author":"Lee","year":"2015","journal-title":"Nat Protoc"},{"key":"2023051312281836000_btad183-B23","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1002\/nav.3800260304","article-title":"Simulation of nonhomogeneous poisson processes by thinning","volume":"26","author":"Lewis","year":"1979","journal-title":"Naval Res Logistics"},{"key":"2023051312281836000_btad183-B24","doi-asserted-by":"crossref","first-page":"e10108","DOI":"10.15252\/msb.202010108","article-title":"Joint cell segmentation and cell type annotation for spatial transcriptomics","volume":"17","author":"Littman","year":"2021","journal-title":"Mol Syst Biol"},{"key":"2023051312281836000_btad183-B25","doi-asserted-by":"crossref","first-page":"743","DOI":"10.1038\/nmeth.2069","article-title":"Single-cell systems biology by super-resolution imaging and combinatorial labeling","volume":"9","author":"Lubeck","year":"2012","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B26","doi-asserted-by":"crossref","first-page":"360","DOI":"10.1038\/nmeth.2892","article-title":"Single-cell in situ RNA profiling by sequential hybridization","volume":"11","author":"Lubeck","year":"2014","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B27","doi-asserted-by":"crossref","first-page":"1743","DOI":"10.1038\/nprot.2013.109","article-title":"Single-molecule mRNA detection and counting in mammalian tissue","volume":"8","author":"Lyubimova","year":"2013","journal-title":"Nat Protoc"},{"key":"2023051312281836000_btad183-B28","author":"Mah","year":"2022"},{"key":"2023051312281836000_btad183-B29","doi-asserted-by":"crossref","first-page":"308","DOI":"10.1038\/s41587-021-01182-1","article-title":"Spatial components of molecular tissue biology","volume":"40","author":"Palla","year":"2022","journal-title":"Nat Biotechnol"},{"key":"2023051312281836000_btad183-B30","doi-asserted-by":"crossref","first-page":"3545","DOI":"10.1038\/s41467-021-23807-4","article-title":"Cell segmentation-free inference of cell types from in situ transcriptomics data","volume":"12","author":"Park","year":"2021","journal-title":"Nat Commun"},{"key":"2023051312281836000_btad183-B31","doi-asserted-by":"crossref","first-page":"1859","DOI":"10.1111\/febs.15572","article-title":"Spage2vec: unsupervised representation of localized spatial gene expression signatures","volume":"288","author":"Partel","year":"2021","journal-title":"FEBS J"},{"key":"2023051312281836000_btad183-B32","first-page":"2825","article-title":"Scikit-learn: machine learning in python","volume":"12","author":"Pedregosa","year":"2011","journal-title":"J Mach Learn Res"},{"key":"2023051312281836000_btad183-B33","doi-asserted-by":"crossref","first-page":"345","DOI":"10.1038\/s41587-021-01044-w","article-title":"Cell segmentation in imaging-based spatial transcriptomics","volume":"40","author":"Petukhov","year":"2022","journal-title":"Nat Biotechnol"},{"key":"2023051312281836000_btad183-B34","doi-asserted-by":"crossref","first-page":"101","DOI":"10.1038\/s41592-019-0631-4","article-title":"Probabilistic cell typing enables fine mapping of closely related cell types in situ","volume":"17","author":"Qian","year":"2020","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B35","doi-asserted-by":"crossref","first-page":"877","DOI":"10.1038\/nmeth.1253","article-title":"Imaging individual mRNA molecules using multiple singly labeled probes","volume":"5","author":"Raj","year":"2008","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B36","doi-asserted-by":"crossref","first-page":"211","DOI":"10.1038\/s41586-021-03634-9","article-title":"Exploring tissue architecture using spatial transcriptomics","volume":"596","author":"Rao","year":"2021","journal-title":"Nature"},{"key":"2023051312281836000_btad183-B37","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1038\/s41467-017-02554-5","article-title":"A general and flexible method for signal extraction from single-cell RNA-seq data","volume":"9","author":"Risso","year":"2018","journal-title":"Nat Commun"},{"key":"2023051312281836000_btad183-B38","doi-asserted-by":"crossref","first-page":"1463","DOI":"10.1126\/science.aaw1219","article-title":"Slide-seq: a scalable technology for measuring genome-wide expression at high spatial resolution","volume":"363","author":"Rodriques","year":"2019","journal-title":"Science"},{"key":"2023051312281836000_btad183-B39","doi-asserted-by":"crossref","DOI":"10.1002\/9781118575574","volume-title":"Multivariate Density Estimation: Theory, Practice, and Visualization","author":"Scott","year":"2015"},{"key":"2023051312281836000_btad183-B40","doi-asserted-by":"crossref","first-page":"78","DOI":"10.1126\/science.aaf2403","article-title":"Visualization and analysis of gene expression in tissue sections by spatial transcriptomics","volume":"353","author":"St\u00e5hl","year":"2016","journal-title":"Science"},{"key":"2023051312281836000_btad183-B41","doi-asserted-by":"crossref","first-page":"790","DOI":"10.1016\/j.tig.2018.07.003","article-title":"Enter the matrix: factorization uncovers knowledge from omics","volume":"34","author":"Stein-O\u2019Brien","year":"2018","journal-title":"Trends Genet"},{"key":"2023051312281836000_btad183-B42","doi-asserted-by":"crossref","first-page":"229","DOI":"10.1038\/s41592-022-01687-w","article-title":"Nonnegative spatial factorization applied to spatial genomics","volume":"20","author":"Townes","year":"2023","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B43","doi-asserted-by":"crossref","first-page":"295","DOI":"10.1186\/s13059-019-1861-6","article-title":"Feature selection and dimension reduction for single-cell RNA-Seq based on a multinomial model","volume":"20","author":"Townes","year":"2019","journal-title":"Genome Biol"},{"key":"2023051312281836000_btad183-B44","doi-asserted-by":"crossref","first-page":"179","DOI":"10.1038\/s41592-021-01343-9","article-title":"Identifying temporal and spatial patterns of variation from multimodal data using MEFISTO","volume":"19","author":"Velten","year":"2022","journal-title":"Nat Methods"},{"key":"2023051312281836000_btad183-B45","doi-asserted-by":"crossref","first-page":"eaat5691","DOI":"10.1126\/science.aat5691","article-title":"Three-dimensional intact-tissue sequencing of single-cell transcriptional states","volume":"361","author":"Wang","year":"2018","journal-title":"Science"},{"key":"2023051312281836000_btad183-B46","doi-asserted-by":"crossref","first-page":"515","DOI":"10.1093\/biostatistics\/kxp008","article-title":"A penalized matrix decomposition, with applications to sparse principal components and canonical correlation analysis","volume":"10","author":"Witten","year":"2009","journal-title":"Biostatistics"},{"key":"2023051312281836000_btad183-B47","doi-asserted-by":"crossref","first-page":"19490","DOI":"10.1073\/pnas.1912459116","article-title":"Spatial transcriptome profiling by MERFISH reveals subcellular RNA compartmentalization and cell cycle-dependent gene expression","volume":"116","author":"Xia","year":"2019","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023051312281836000_btad183-B48","first-page":"1881","article-title":"Gaussian-process factor analysis for low-dimensional single-trial analysis of neural population activity","volume":"21","author":"Yu","year":"2008","journal-title":"Advances in Neural Information Processing Systems"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad183\/49841450\/btad183.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/5\/btad183\/50310876\/btad183.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/5\/btad183\/50310876\/btad183.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,12,10]],"date-time":"2023-12-10T09:37:40Z","timestamp":1702201060000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad183\/7114027"}},"subtitle":[],"editor":[{"given":"Christina","family":"Kendziorski","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,4,11]]},"references-count":48,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2023,5,4]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad183","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2021.11.04.467354","asserted-by":"object"}]},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,5,1]]},"published":{"date-parts":[[2023,4,11]]},"article-number":"btad183"}}