{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,8]],"date-time":"2026-05-08T07:47:59Z","timestamp":1778226479158,"version":"3.51.4"},"reference-count":41,"publisher":"Oxford University Press (OUP)","issue":"4","license":[{"start":{"date-parts":[[2023,5,24]],"date-time":"2023-05-24T00:00:00Z","timestamp":1684886400000},"content-version":"vor","delay-in-days":1,"URL":"https:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001809","name":"Natural Science Foundation of China","doi-asserted-by":"publisher","award":["31601072"],"award-info":[{"award-number":["31601072"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"Natural Science Foundation of China","doi-asserted-by":"publisher","award":["62271175"],"award-info":[{"award-number":["62271175"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"Natural Science Foundation of China","doi-asserted-by":"publisher","award":["62072142"],"award-info":[{"award-number":["62072142"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,7,20]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Single-cell RNA sequencing (scRNA-seq) detects whole transcriptome signals for large amounts of individual cells and is powerful for determining cell-to-cell differences and investigating the functional characteristics of various cell types. scRNA-seq datasets are usually sparse and highly noisy. Many steps in the scRNA-seq analysis workflow, including reasonable gene selection, cell clustering and annotation, as well as discovering the underlying biological mechanisms from such datasets, are difficult. In this study, we proposed an scRNA-seq analysis method based on the latent Dirichlet allocation (LDA) model. The LDA model estimates a series of latent variables, i.e. putative functions (PFs), from the input raw cell\u2013gene data. Thus, we incorporated the \u2018cell-function-gene\u2019 three-layer framework into scRNA-seq analysis, as this framework is capable of discovering latent and complex gene expression patterns via a built-in model approach and obtaining biologically meaningful results through a data-driven functional interpretation process. We compared our method with four classic methods on seven benchmark scRNA-seq datasets. The LDA-based method performed best in the cell clustering test in terms of both accuracy and purity. By analysing three complex public datasets, we demonstrated that our method could distinguish cell types with multiple levels of functional specialization, and precisely reconstruct cell development trajectories. Moreover, the LDA-based method accurately identified the representative PFs and the representative genes for the cell types\/cell stages, enabling data-driven cell cluster annotation and functional interpretation. According to the literature, most of the previously reported marker\/functionally relevant genes were recognized.<\/jats:p>","DOI":"10.1093\/bib\/bbad199","type":"journal-article","created":{"date-parts":[[2023,5,25]],"date-time":"2023-05-25T00:09:34Z","timestamp":1684973374000},"source":"Crossref","is-referenced-by-count":19,"title":["An interpretable single-cell RNA sequencing data clustering method based on latent Dirichlet allocation"],"prefix":"10.1093","volume":"24","author":[{"given":"Qi","family":"Yang","sequence":"first","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1799-5529","authenticated-orcid":false,"given":"Zhaochun","family":"Xu","sequence":"additional","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3676-0584","authenticated-orcid":false,"given":"Wenyang","family":"Zhou","sequence":"additional","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8888-8005","authenticated-orcid":false,"given":"Pingping","family":"Wang","sequence":"additional","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Qinghua","family":"Jiang","sequence":"additional","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8826-0269","authenticated-orcid":false,"given":"Liran","family":"Juan","sequence":"additional","affiliation":[{"name":"School of Life Science and Technology, Harbin Institute of Technology , Harbin 150001 , China"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2023,5,23]]},"reference":[{"key":"2023072020034552300_ref1","doi-asserted-by":"crossref","first-page":"1396","DOI":"10.1126\/science.1254257","article-title":"Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma","volume":"344","author":"Patel","year":"2014","journal-title":"Science"},{"key":"2023072020034552300_ref2","doi-asserted-by":"crossref","DOI":"10.1126\/science.aah4573","article-title":"Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors","volume":"356","author":"Villani","year":"2017","journal-title":"Science"},{"key":"2023072020034552300_ref3","doi-asserted-by":"crossref","first-page":"1138","DOI":"10.1126\/science.aaa1934","article-title":"Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq","volume":"347","author":"Zeisel","year":"2015","journal-title":"Science"},{"key":"2023072020034552300_ref4","doi-asserted-by":"crossref","first-page":"381","DOI":"10.1038\/nbt.2859","article-title":"The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells","volume":"32","author":"Trapnell","year":"2014","journal-title":"Nat Biotechnol"},{"key":"2023072020034552300_ref5","doi-asserted-by":"crossref","first-page":"140","DOI":"10.1186\/s12859-016-0984-y","article-title":"pcaReduce: hierarchical clustering of single cell transcriptional profiles","volume":"17","author":"Zurauskiene","year":"2016","journal-title":"BMC Bioinformatics"},{"key":"2023072020034552300_ref6","doi-asserted-by":"crossref","first-page":"483","DOI":"10.1038\/nmeth.4236","article-title":"SC3: consensus clustering of single-cell RNA-seq data","volume":"14","author":"Kiselev","year":"2017","journal-title":"Nat Methods"},{"key":"2023072020034552300_ref7","doi-asserted-by":"crossref","first-page":"379","DOI":"10.1038\/nmeth.4662","article-title":"FateID infers cell fate bias in multipotent progenitors from single-cell RNA-seq data","volume":"15","author":"Herman","year":"2018","journal-title":"Nat Methods"},{"key":"2023072020034552300_ref8","doi-asserted-by":"crossref","first-page":"184","DOI":"10.1016\/j.cell.2015.05.047","article-title":"Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis","volume":"162","author":"Levine","year":"2015","journal-title":"Cell"},{"key":"2023072020034552300_ref9","doi-asserted-by":"crossref","first-page":"411","DOI":"10.1038\/nbt.4096","article-title":"Integrating single-cell transcriptomic data across different conditions, technologies, and species","volume":"36","author":"Butler","year":"2018","journal-title":"Nat Biotechnol"},{"key":"2023072020034552300_ref10","doi-asserted-by":"crossref","first-page":"P10008","DOI":"10.1088\/1742-5468\/2008\/10\/P10008","article-title":"Fast unfolding of communities in large networks","volume":"2008","author":"Blondel","year":"2008","journal-title":"Journal of Statistical Mechanics: Theory and Experiment"},{"key":"2023072020034552300_ref11","doi-asserted-by":"crossref","DOI":"10.1093\/nar\/gkz826","article-title":"Latent cellular analysis robustly reveals subtle diversity in large-scale single-cell RNA-seq data","volume":"47","author":"Cheng","year":"2019","journal-title":"Nucleic Acids Res"},{"key":"2023072020034552300_ref12","doi-asserted-by":"crossref","first-page":"163","DOI":"10.1038\/s41590-018-0276-y","article-title":"Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage","volume":"20","author":"Aran","year":"2019","journal-title":"Nat Immunol"},{"key":"2023072020034552300_ref13","doi-asserted-by":"crossref","first-page":"1095","DOI":"10.1038\/s41587-021-00896-6","article-title":"Gene signature extraction and cell identity recognition at the single-cell level with cell-ID","volume":"39","author":"Cortal","year":"2021","journal-title":"Nat Biotechnol"},{"key":"2023072020034552300_ref14","first-page":"993","article-title":"Latent Dirichlet allocation","volume":"3","author":"Blei","year":"2003","journal-title":"Journal of Machine Learning Research"},{"key":"2023072020034552300_ref15","doi-asserted-by":"crossref","first-page":"150","DOI":"10.1186\/s13059-020-02033-z","article-title":"Approaches for integrating heterogeneous RNA-seq data reveal cross-talk between microbes and genes in asthmatic patients","volume":"21","author":"Spakowicz","year":"2020","journal-title":"Genome Biol"},{"key":"2023072020034552300_ref16","doi-asserted-by":"crossref","first-page":"4543","DOI":"10.1093\/bioinformatics\/btz266","article-title":"Discovering novel mutation signatures by latent Dirichlet allocation with variational Bayes inference","volume":"35","author":"Matsutani","year":"2019","journal-title":"Bioinformatics"},{"key":"2023072020034552300_ref17","doi-asserted-by":"crossref","first-page":"519","DOI":"10.1186\/s12859-018-2530-6","article-title":"A latent allocation model for the analysis of microbial composition and disease","volume":"19","author":"Abe","year":"2018","journal-title":"BMC Bioinformatics"},{"key":"2023072020034552300_ref18","doi-asserted-by":"crossref","first-page":"962","DOI":"10.1186\/s12864-016-3257-2","article-title":"MetaTopics: an integration tool to analyze microbial community profile by topic model","volume":"18","author":"Yan","year":"2017","journal-title":"BMC Genomics"},{"key":"2023072020034552300_ref19","doi-asserted-by":"crossref","first-page":"13738","DOI":"10.1073\/pnas.1608041113","article-title":"Topic modeling for untargeted substructure exploration in metabolomics","volume":"113","author":"Van der Hooft","year":"2016","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2023072020034552300_ref20","doi-asserted-by":"crossref","first-page":"144","DOI":"10.1016\/j.compbiolchem.2014.08.019","article-title":"Exploring the complexity of pathway-drug relationships using latent Dirichlet allocation","volume":"53 Pt A","author":"Pratanwanich","year":"2014","journal-title":"Comput Biol Chem"},{"key":"2023072020034552300_ref21","doi-asserted-by":"crossref","first-page":"11","DOI":"10.1242\/bio.059001","article-title":"Dirichlet process mixture models for single-cell RNA-seq clustering","volume":"11","author":"Adossa","year":"2022","journal-title":"Biology Open"},{"key":"2023072020034552300_ref22","doi-asserted-by":"crossref","first-page":"543","DOI":"10.1007\/s12561-021-09304-8","article-title":"Penalized latent Dirichlet allocation model in single-cell RNA sequencing","volume":"13","author":"Wu","year":"2021","journal-title":"Statistics in Biosciences"},{"key":"2023072020034552300_ref23","doi-asserted-by":"crossref","first-page":"397","DOI":"10.1038\/s41592-019-0367-1","article-title":"cisTopic: cis-regulatory topic modeling on single-cell ATAC-seq data","volume":"16","author":"Bravo Gonz\u00e1lez-Blas","year":"2019","journal-title":"Nat Methods"},{"key":"2023072020034552300_ref24","doi-asserted-by":"crossref","first-page":"538","DOI":"10.1038\/nature25981","article-title":"The cis-regulatory dynamics of embryonic development at single-cell resolution","volume":"555","author":"Cusanovich","year":"2018","journal-title":"Nature"},{"key":"2023072020034552300_ref25","article-title":"Visualizing the structure of RNA-seq expression data using grade of membership models","volume":"13","author":"Dey","year":"2017","journal-title":"PLoS Genet"},{"key":"2023072020034552300_ref26","doi-asserted-by":"crossref","first-page":"363","DOI":"10.1186\/s12859-016-1175-6","article-title":"CellTree: an R\/bioconductor package to infer the hierarchical structure of cell populations from single-cell RNA-seq data","volume":"17","author":"duVerle","year":"2016","journal-title":"BMC Bioinformatics"},{"key":"2023072020034552300_ref27","doi-asserted-by":"crossref","first-page":"479","DOI":"10.1038\/s41592-019-0425-8","article-title":"Benchmarking single cell RNA-sequencing analysis pipelines using mixture control experiments","volume":"16","author":"Tian","year":"2019","journal-title":"Nat Methods"},{"key":"2023072020034552300_ref28","doi-asserted-by":"crossref","first-page":"193","DOI":"10.1007\/BF01908075","article-title":"Comparing partitions","volume":"2","author":"Hubert","year":"1985","journal-title":"Journal of Classification"},{"key":"2023072020034552300_ref29","volume-title":"International Conference on Neural Information Processing Systems","author":"Hoffman","year":"2010"},{"key":"2023072020034552300_ref30","article-title":"Graph embedding and Gaussian mixture variational autoencoder network for end-to-end analysis of single-cell RNA sequencing data","volume":"3","author":"Xu","year":"2023","journal-title":"Cell Rep Methods"},{"key":"2023072020034552300_ref31","doi-asserted-by":"crossref","first-page":"1053","DOI":"10.1038\/s41592-018-0229-2","article-title":"Deep generative modeling for single-cell transcriptomics","volume":"15","author":"Lopez","year":"2018","journal-title":"Nat Methods"},{"key":"2023072020034552300_ref32","doi-asserted-by":"crossref","first-page":"390","DOI":"10.1038\/s41467-018-07931-2","article-title":"Single-cell RNA-seq denoising using a deep count autoencoder","volume":"10","author":"Eraslan","year":"2019","journal-title":"Nat Commun"},{"key":"2023072020034552300_ref33","volume-title":"MALLET: A Machine Learning for Language Toolkit","author":"Mccallum","year":"2002"},{"key":"2023072020034552300_ref34","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1126\/science.aad0501","article-title":"Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq","volume":"352","author":"Tirosh","year":"2016","journal-title":"Science"},{"key":"2023072020034552300_ref35","doi-asserted-by":"crossref","DOI":"10.1126\/science.aay3224","article-title":"A cell atlas of human thymic development defines T cell repertoire formation","volume":"367","author":"Park","year":"2020","journal-title":"Science"},{"key":"2023072020034552300_ref36","doi-asserted-by":"crossref","first-page":"14049","DOI":"10.1038\/ncomms14049","article-title":"Massively parallel digital transcriptional profiling of single cells","volume":"8","author":"Zheng","year":"2017","journal-title":"Nat Commun"},{"key":"2023072020034552300_ref37","doi-asserted-by":"crossref","first-page":"194","DOI":"10.1186\/s13059-019-1795-z","article-title":"A comparison of automatic cell identification methods for single-cell RNA sequencing data","volume":"20","author":"Abdelaal","year":"2019","journal-title":"Genome Biol"},{"key":"2023072020034552300_ref38","doi-asserted-by":"crossref","first-page":"3573","DOI":"10.1016\/j.cell.2021.04.048","article-title":"Integrated analysis of multimodal single-cell data","volume":"184","author":"Hao","year":"2021","journal-title":"Cell"},{"key":"2023072020034552300_ref39","doi-asserted-by":"crossref","first-page":"531","DOI":"10.3390\/genes10070531","article-title":"SCINA: a semi-supervised subtyping algorithm of single cells and bulk samples","volume":"10","author":"Zhang","year":"2019","journal-title":"Genes (Basel)"},{"key":"2023072020034552300_ref40","doi-asserted-by":"crossref","first-page":"220","DOI":"10.1186\/s13059-017-1349-1","article-title":"xCell: digitally portraying the tissue cellular heterogeneity landscape","volume":"18","author":"Aran","year":"2017","journal-title":"Genome Biol"},{"key":"2023072020034552300_ref41","doi-asserted-by":"crossref","first-page":"174","DOI":"10.1186\/s13059-017-1305-0","article-title":"Splatter: simulation of single-cell RNA sequencing data","volume":"18","author":"Zappia","year":"2017","journal-title":"Genome Biol"}],"container-title":["Briefings in Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/24\/4\/bbad199\/50916599\/bbad199.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/24\/4\/bbad199\/50916599\/bbad199.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,7,20]],"date-time":"2023-07-20T20:05:41Z","timestamp":1689883541000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bib\/article\/doi\/10.1093\/bib\/bbad199\/7176310"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,5,23]]},"references-count":41,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2023,7,20]]}},"URL":"https:\/\/doi.org\/10.1093\/bib\/bbad199","relation":{},"ISSN":["1467-5463","1477-4054"],"issn-type":[{"value":"1467-5463","type":"print"},{"value":"1477-4054","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,7]]},"published":{"date-parts":[[2023,5,23]]},"article-number":"bbad199"}}