{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,19]],"date-time":"2026-06-19T02:57:11Z","timestamp":1781837831361,"version":"3.54.5"},"reference-count":43,"publisher":"Oxford University Press (OUP)","issue":"15","license":[{"start":{"date-parts":[[2018,12,7]],"date-time":"2018-12-07T00:00:00Z","timestamp":1544140800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["11571349"],"award-info":[{"award-number":["11571349"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["91630314"],"award-info":[{"award-number":["91630314"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["81673833"],"award-info":[{"award-number":["81673833"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Strategic Priority Research Program of Chinese Academy of Sciences","award":["XDB13050000"],"award-info":[{"award-number":["XDB13050000"]}]},{"name":"National Center for Mathematics and Interdisciplinary Sciences of Chinese Academy of Sciences"},{"name":"LSC of Chinese Academy of Sciences"},{"DOI":"10.13039\/501100004739","name":"Youth Innovation Promotion Association of Chinese Academy of Sciences","doi-asserted-by":"crossref","id":[{"id":"10.13039\/501100004739","id-type":"DOI","asserted-by":"crossref"}]},{"name":"Mathematical Biosciences Institute"},{"name":"MBI"},{"DOI":"10.13039\/100006928","name":"Ohio State University","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100006928","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["DMS 1440386"],"award-info":[{"award-number":["DMS 1440386"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,8,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Visualizing and reconstructing cell developmental trajectories intrinsically embedded in high-dimensional expression profiles of single-cell RNA sequencing (scRNA-seq) snapshot data are computationally intriguing, but challenging.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We propose DensityPath, an algorithm allowing (i) visualization of the intrinsic structure of scRNA-seq data on an embedded 2-d space and (ii) reconstruction of an optimal cell state-transition path on the density landscape. DensityPath powerfully handles high dimensionality and heterogeneity of scRNA-seq data by (i) revealing the intrinsic structures of data, while adopting a non-linear dimension reduction algorithm, termed elastic embedding, which can preserve both local and global structures of the data; and (ii) extracting the topological features of high-density, level-set clusters from a single-cell multimodal density landscape of transcriptional heterogeneity, as the representative cell states. DensityPath reconstructs the optimal cell state-transition path by finding the geodesic minimum spanning tree of representative cell states on the density landscape, establishing a least action path with the minimum-transition-energy of cell fate decisions. We demonstrate that DensityPath can ably reconstruct complex trajectories of cell development, e.g. those with multiple bifurcating and trifurcating branches, while maintaining computational efficiency. Moreover, DensityPath has high accuracy for pseudotime calculation and branch assignment on real scRNA-seq, as well as simulated datasets. DensityPath is robust to parameter choices, as well as permutations of data.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>DensityPath software is available at https:\/\/github.com\/ucasdp\/DensityPath.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Supplementary information<\/jats:title>\n                    <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty1009","type":"journal-article","created":{"date-parts":[[2018,12,6]],"date-time":"2018-12-06T17:37:40Z","timestamp":1544117860000},"page":"2593-2601","source":"Crossref","is-referenced-by-count":31,"title":["DensityPath: an algorithm to visualize and reconstruct cell state-transition path on density landscape for single-cell RNA sequencing data"],"prefix":"10.1093","volume":"35","author":[{"given":"Ziwei","family":"Chen","sequence":"first","affiliation":[{"name":"NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing"},{"name":"University of Chinese Academy of Sciences, Beijing"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Shaokun","family":"An","sequence":"additional","affiliation":[{"name":"NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing"},{"name":"University of Chinese Academy of Sciences, Beijing"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiangqi","family":"Bai","sequence":"additional","affiliation":[{"name":"NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing"},{"name":"University of Chinese Academy of Sciences, Beijing"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Fuzhou","family":"Gong","sequence":"additional","affiliation":[{"name":"NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing"},{"name":"University of Chinese Academy of Sciences, Beijing"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Liang","family":"Ma","sequence":"additional","affiliation":[{"name":"University of Chinese Academy of Sciences, Beijing"},{"name":"Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Lin","family":"Wan","sequence":"additional","affiliation":[{"name":"NCMIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing"},{"name":"University of Chinese Academy of Sciences, Beijing"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2018,12,7]]},"reference":[{"key":"2023062803480124900_bty1009-B1","doi-asserted-by":"crossref","first-page":"714","DOI":"10.1016\/j.cell.2014.04.005","article-title":"Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development","volume":"157","author":"Bendall","year":"2014","journal-title":"Cell"},{"key":"2023062803480124900_bty1009-B2","doi-asserted-by":"crossref","first-page":"999","DOI":"10.1016\/j.jmva.2005.05.004","article-title":"Kernel estimation of density level sets","volume":"97","author":"Cadre","year":"2006","journal-title":"J. Multivar. Anal"},{"key":"2023062803480124900_bty1009-B3","first-page":"167","author":"Carreira-Perpi\u00f1\u00e1n","year":"2010"},{"key":"2023062803480124900_bty1009-B4","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1016\/j.mam.2017.10.004","article-title":"Constructing cell lineages from single-cell transcriptomes","volume":"59","author":"Chen","year":"2018","journal-title":"Mol. Aspects Med"},{"key":"2023062803480124900_bty1009-B5","doi-asserted-by":"crossref","first-page":"161","DOI":"10.1080\/24709360.2017.1396742","article-title":"A tutorial on kernel density estimation and recent advances","volume":"1","author":"Chen","year":"2017","journal-title":"Biostat. Epidemiol"},{"key":"2023062803480124900_bty1009-B6","doi-asserted-by":"crossref","first-page":"790","DOI":"10.1109\/34.400568","article-title":"Mean Shift, Mode Seeking, and Clustering","volume":"17","author":"Cheng","year":"1995","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell"},{"key":"2023062803480124900_bty1009-B7","doi-asserted-by":"crossref","first-page":"185","DOI":"10.1016\/j.cels.2018.05.017","article-title":"Generalizable and scalable visualization of single-cell data using neural networks","volume":"7","author":"Cho","year":"2018","journal-title":"Cell Syst"},{"key":"2023062803480124900_bty1009-B8","doi-asserted-by":"crossref","first-page":"603","DOI":"10.1109\/34.1000236","article-title":"Mean shift: a robust approach toward feature space analysis","volume":"24","author":"Comaniciu","year":"2002","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell"},{"key":"2023062803480124900_bty1009-B9","doi-asserted-by":"crossref","first-page":"I102","DOI":"10.1093\/bioinformatics\/btx232","article-title":"HopLand: single-cell pseudotime recovery using continuous Hopfield network-based modeling of Waddington\u2019s epigenetic landscape","volume":"33","author":"Guo","year":"2017","journal-title":"Bioinformatics"},{"key":"2023062803480124900_bty1009-B10","doi-asserted-by":"crossref","first-page":"845","DOI":"10.1038\/nmeth.3971","article-title":"Diffusion pseudotime robustly reconstructs lineage branching","volume":"13","author":"Haghverdi","year":"2016","journal-title":"Nat. Methods"},{"key":"2023062803480124900_bty1009-B11","volume-title":"Clustering Algorithms.","author":"Hartigan","year":"1975"},{"key":"2023062803480124900_bty1009-B12","doi-asserted-by":"crossref","DOI":"10.1007\/978-0-387-84858-7","volume-title":"The Elements of Statistical Learning: Data Mining, inference and Prediction","author":"Hastie","year":"2009"},{"key":"2023062803480124900_bty1009-B13","doi-asserted-by":"crossref","first-page":"e98679","DOI":"10.1371\/journal.pone.0098679","article-title":"ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software","volume":"9","author":"Jacomy","year":"2014","journal-title":"PLoS One"},{"key":"2023062803480124900_bty1009-B14","doi-asserted-by":"crossref","first-page":"e117","DOI":"10.1093\/nar\/gkw430","article-title":"TSCAN: pseudo-time reconstruction and evaluation in single-cell RNA-seq analysis","volume":"44","author":"Ji","year":"2016","journal-title":"Nucleic Acids Res"},{"key":"2023062803480124900_bty1009-B15","doi-asserted-by":"crossref","first-page":"2077","DOI":"10.1093\/bioinformatics\/bty058","article-title":"scEpath: energy landscape-based inference of transition probabilities and cellular trajectories from single-cell transcriptomic data","volume":"34","author":"Jin","year":"2018","journal-title":"Bioinformatics"},{"key":"2023062803480124900_bty1009-B16","first-page":"1839","author":"Kim","year":"2016"},{"key":"2023062803480124900_bty1009-B17","doi-asserted-by":"crossref","first-page":"17","DOI":"10.1242\/dev.133058","article-title":"Understanding development and stem cells using single cell-based analyses of gene expression","volume":"144","author":"Kumar","year":"2017","journal-title":"Development"},{"key":"2023062803480124900_bty1009-B18","doi-asserted-by":"crossref","first-page":"32","DOI":"10.1016\/j.coisb.2018.02.009","article-title":"Exploring intermediate cell states through the lens of single cells","volume":"9","author":"MacLean","year":"2018","journal-title":"Curr. Opin. Syst. Biol"},{"key":"2023062803480124900_bty1009-B19","doi-asserted-by":"crossref","first-page":"2227","DOI":"10.1109\/TPAMI.2016.2635657","article-title":"Principal graph and structure learning based on reversed graph embedding","volume":"39","author":"Mao","year":"2017","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell"},{"key":"2023062803480124900_bty1009-B20","doi-asserted-by":"crossref","first-page":"E5643","DOI":"10.1073\/pnas.1408993111","article-title":"Bifurcation analysis of single-cell gene expression data reveals epigenetic landscape","volume":"111","author":"Marco","year":"2014","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062803480124900_bty1009-B21","article-title":"Visualizing transitions and structure for high dimensional data exploration","author":"Moon","year":"2017"},{"key":"2023062803480124900_bty1009-B22","volume-title":"Level Set Methods and Dynamic Implicit Surfaces. Applied Mathematical Sciences","author":"Osher","year":"2002"},{"key":"2023062803480124900_bty1009-B23","doi-asserted-by":"crossref","first-page":"1663","DOI":"10.1016\/j.cell.2015.11.013","article-title":"Transcriptional heterogeneity and lineage commitment in myeloid progenitors","volume":"163","author":"Paul","year":"2015","journal-title":"Cell"},{"key":"2023062803480124900_bty1009-B24","doi-asserted-by":"crossref","first-page":"1012","DOI":"10.1016\/j.cell.2016.03.023","article-title":"Single-cell RNA-Seq reveals lineage and X chromosome dynamics in human preimplantation embryos","volume":"165","author":"Petropoulos","year":"2016","journal-title":"Cell"},{"key":"2023062803480124900_bty1009-B25","doi-asserted-by":"crossref","first-page":"979","DOI":"10.1038\/nmeth.4402","article-title":"Reversed graph embedding resolves complex single-cell trajectories","volume":"14","author":"Qiu","year":"2017","journal-title":"Nat. Methods"},{"key":"2023062803480124900_bty1009-B26","doi-asserted-by":"crossref","first-page":"846","DOI":"10.1080\/01621459.1971.10482356","article-title":"Objective criteria for the evaluation of clustering methods","volume":"66","author":"Rand","year":"1971","journal-title":"J. Am. Stat. Assoc"},{"key":"2023062803480124900_bty1009-B27","doi-asserted-by":"crossref","first-page":"551","DOI":"10.1038\/nbt.3854","article-title":"Single-cell topological RNA-seq analysis reveals insights into cellular differentiation and development","volume":"35","author":"Rizvi","year":"2017","journal-title":"Nat. Biotechnol"},{"key":"2023062803480124900_bty1009-B28","doi-asserted-by":"crossref","DOI":"10.1101\/191056","article-title":"Reconstruction of developmental landscapes by optimal-transport analysis of single-cell gene expression sheds light on cellular reprogramming","author":"Schiebinger","year":"2017"},{"key":"2023062803480124900_bty1009-B29","doi-asserted-by":"crossref","first-page":"637","DOI":"10.1038\/nbt.3569","article-title":"Wishbone identifies bifurcating developmental trajectories from single-cell data","volume":"34","author":"Setty","year":"2016","journal-title":"Nat. Biotechnol"},{"key":"2023062803480124900_bty1009-B30","doi-asserted-by":"crossref","first-page":"683","DOI":"10.1111\/j.2517-6161.1991.tb01857.x","article-title":"A reliable data-based bandwidth selection method for kernel density-estimation","volume":"53","author":"Sheather","year":"1991","journal-title":"J. R. Stat. Soc. Series B Stat. Methodol"},{"key":"2023062803480124900_bty1009-B31","first-page":"91","author":"Singh","year":"2007"},{"key":"2023062803480124900_bty1009-B32","doi-asserted-by":"crossref","first-page":"331","DOI":"10.1038\/nature21350","article-title":"Scaling single-cell genomics from phenomenology to mechanism","volume":"541","author":"Tanay","year":"2017","journal-title":"Nature"},{"key":"2023062803480124900_bty1009-B33","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":"2023062803480124900_bty1009-B34","first-page":"477","author":"Vladymyrov","year":"2013"},{"key":"2023062803480124900_bty1009-B35","first-page":"345","article-title":"Partial-hessian strategies for fast learning of nonlinear embeddings","volume-title":"Proceedings of the 29th International Conference on Machine Learning (ICML-12)","author":"Vladymyrov","year":"2012"},{"key":"2023062803480124900_bty1009-B36","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1080\/00018732.2015.1037068","article-title":"Landscape and flux theory of non-equilibrium dynamical systems with application to biology","volume":"64","author":"Wang","year":"2015","journal-title":"Adv. Phys"},{"key":"2023062803480124900_bty1009-B37","doi-asserted-by":"crossref","first-page":"8257","DOI":"10.1073\/pnas.1017017108","article-title":"Quantifying the Waddington landscape and biological paths for development and differentiation","volume":"108","author":"Wang","year":"2011","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062803480124900_bty1009-B38","doi-asserted-by":"crossref","first-page":"501","DOI":"10.1146\/annurev-statistics-031017-100045","article-title":"Topological data analysis","volume":"5","author":"Wasserman","year":"2018","journal-title":"Annu. Rev. Stat. Appl"},{"key":"2023062803480124900_bty1009-B39","volume-title":"All of Nonparametric Statistics","author":"Wassermann","year":"2006"},{"key":"2023062803480124900_bty1009-B40","doi-asserted-by":"crossref","first-page":"E2467","DOI":"10.1073\/pnas.1714723115","article-title":"Fundamental limits on dynamic inference from single-cell snapshots","volume":"115","author":"Weinreb","year":"2018","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062803480124900_bty1009-B41","doi-asserted-by":"crossref","first-page":"1246","DOI":"10.1093\/bioinformatics\/btx792","article-title":"SPRING: a kinetic interface for visualizing high dimensional single-cell expression data","volume":"34","author":"Weinreb","year":"2018","journal-title":"Bioinformatics"},{"key":"2023062803480124900_bty1009-B42","doi-asserted-by":"crossref","first-page":"1665","DOI":"10.1214\/aoms\/1177696810","article-title":"On choosing a delta-sequence","volume":"41","author":"Woodroofe","year":"1970","journal-title":"Ann. Math. Stat"},{"key":"2023062803480124900_bty1009-B43","article-title":"Topslam: Waddington landscape recovery for single cell experiments","author":"Zwiessele","year":"2017"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/15\/2593\/50722416\/bioinformatics_35_15_2593.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/15\/2593\/50722416\/bioinformatics_35_15_2593.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,12]],"date-time":"2024-07-12T23:03:53Z","timestamp":1720825433000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/15\/2593\/5233001"}},"subtitle":[],"editor":[{"given":"Bonnie","family":"Berger","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2018,12,7]]},"references-count":43,"journal-issue":{"issue":"15","published-print":{"date-parts":[[2019,8,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty1009","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/276311","asserted-by":"object"}]},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2019,8,1]]},"published":{"date-parts":[[2018,12,7]]}}}