{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,30]],"date-time":"2026-03-30T03:42:04Z","timestamp":1774842124111,"version":"3.50.1"},"reference-count":32,"publisher":"Oxford University Press (OUP)","issue":"Supplement_1","license":[{"start":{"date-parts":[[2022,6,27]],"date-time":"2022-06-27T00:00:00Z","timestamp":1656288000000},"content-version":"vor","delay-in-days":3,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2022,6,24]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Improvements in single-cell RNA-seq technologies mean that studies measuring multiple experimental conditions, such as time series, have become more common. At present, few computational methods exist to infer time series-specific transcriptome changes, and such studies have therefore typically used unsupervised pseudotime methods. While these methods identify cell subpopulations and the transitions between them, they are not appropriate for identifying the genes that vary coherently along the time series. In addition, the orderings they estimate are based only on the major sources of variation in the data, which may not correspond to the processes related to the time labels.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>We introduce psupertime, a supervised pseudotime approach based on a regression model, which explicitly uses time-series labels as input. It identifies genes that vary coherently along a time series, in addition to pseudotime values for individual cells, and a classifier that can be used to estimate labels for new data with unknown or differing labels. We show that psupertime outperforms benchmark classifiers in terms of identifying time-varying genes and provides better individual cell orderings than popular unsupervised pseudotime techniques. psupertime is applicable to any single-cell RNA-seq dataset with sequential labels (e.g. principally time series but also drug dosage and disease progression), derived from either experimental design and provides a fast, interpretable tool for targeted identification of genes varying along with specific biological processes.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>R package available at github.com\/wmacnair\/psupertime and code for results reproduction at github.com\/wmacnair\/psupplementary.<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btac227","type":"journal-article","created":{"date-parts":[[2022,4,14]],"date-time":"2022-04-14T11:10:15Z","timestamp":1649934615000},"page":"i290-i298","source":"Crossref","is-referenced-by-count":39,"title":["psupertime: supervised pseudotime analysis for time-series single-cell RNA-seq data"],"prefix":"10.1093","volume":"38","author":[{"given":"Will","family":"Macnair","sequence":"first","affiliation":[{"name":"Institute of Molecular Systems Biology, Department of Biology, ETH Zurich, Zurich 8093, Switzerland"}]},{"given":"Revant","family":"Gupta","sequence":"additional","affiliation":[{"name":"Inner Medicine I, Faculty of Medicine, University of T\u00fcbingen, University Hospital T\u00fcbingen, 72074, Germany"}]},{"given":"Manfred","family":"Claassen","sequence":"additional","affiliation":[{"name":"Inner Medicine I, Faculty of Medicine, University of T\u00fcbingen, University Hospital T\u00fcbingen, 72074, Germany"},{"name":"Department of Computer Science, University of T\u00fcbingen , T\u00fcbingen 72074, Germany"}]}],"member":"286","published-online":{"date-parts":[[2022,6,27]]},"reference":[{"key":"2023041407533932900_","article-title":"Gene set enrichment analysis with topGO","volume":"27, 1\u201326","author":"Alexa","year":"2009","journal-title":"Bioconductor Improv."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"1464","DOI":"10.1002\/sim.4484","article-title":"L1 penalized continuation ratio models for ordinal response prediction using high-dimensional datasets","volume":"31","author":"Archer","year":"2012","journal-title":"Stat. Med."},{"key":"2023041407533932900_","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":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"1408","DOI":"10.1038\/s41587-020-0591-3","article-title":"Generalizing RNA velocity to transient cell states through dynamical modeling","volume":"38","author":"Bergen","year":"2020","journal-title":"Nat. Biotechnol"},{"key":"2023041407533932900_","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":"2023041407533932900_","first-page":"96","author":"Caruana","year":"2008"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"1141","DOI":"10.12688\/f1000research.15666.2","article-title":"A systematic performance evaluation of clustering methods for single-cell RNA-seq data","volume":"7","author":"Du\u00f2","year":"2018","journal-title":"F1000Res"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1016\/j.cell.2017.09.004","article-title":"Single-cell analysis of human pancreas reveals transcriptional signatures of aging and somatic mutation patterns","volume":"171","author":"Enge","year":"2017","journal-title":"Cell"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v033.i01","article-title":"Regularization paths for generalized linear models via coordinate descent","volume":"33","author":"Friedman","year":"2010","journal-title":"J. Stat. Softw."},{"key":"2023041407533932900_","doi-asserted-by":"publisher","first-page":"37","DOI":"10.1016\/j.cels.2017.10.012","article-title":"Unsupervised trajectory analysis of single-cell RNA-seq and imaging data reveals alternative tuft cell origins in the gut","volume":"6","author":"Herring","year":"2018","journal-title":"Cell Syst."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"160109","DOI":"10.1038\/sdata.2016.109","article-title":"An atlas of transcriptional, chromatin accessibility, and surface marker changes in human mesoderm development","volume":"3","author":"Koh","year":"2016","journal-title":"Sci. Data"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"5619","DOI":"10.1038\/ncomms6619","article-title":"An epigenomic roadmap to induced pluripotency reveals DNA methylation as a reprogramming modulator","volume":"5","author":"Lee","year":"2014","journal-title":"Nat. Commun."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"605","DOI":"10.1002\/dvdy.22556","article-title":"Essential role of clusterin in pancreas regeneration","volume":"240","author":"Lee","year":"2011","journal-title":"Dev. Dyn."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"858","DOI":"10.1016\/j.stem.2017.03.007","article-title":"Single-cell RNA-Seq analysis maps development of human germline cells and gonadal niche interactions","volume":"20","author":"Li","year":"2017","journal-title":"Cell Stem Cell"},{"key":"2023041407533932900_","first-page":"2122","article-title":"A step-by-step workflow for low-level analysis of single-cell RNA-seq data with bioconductor","volume":"5","author":"Lun","year":"2016","journal-title":"F1000Res"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"109","DOI":"10.1111\/j.2517-6161.1980.tb01109.x","article-title":"Regression models for ordinal data","volume":"42","author":"McCullagh","year":"1980","journal-title":"J. R. Stat. Soc. Series B Stat. Methodol"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"48","DOI":"10.1002\/stem.755","article-title":"Concise review: pluripotency and the transcriptional inactivation of the female mammalian X chromosome","volume":"30","author":"Minkovsky","year":"2012","journal-title":"Stem Cells"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"224","DOI":"10.1016\/0167-4838(93)90087-8","article-title":"Identification of the characteristic amino-acid sequence for human \u03b1-amylase encoded by the AMY2B gene","volume":"1203","author":"Omichi","year":"1993","journal-title":"Biochim. Biophys. Acta"},{"key":"2023041407533932900_","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":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"1194","DOI":"10.1016\/j.cmet.2017.04.003","article-title":"Deciphering pancreatic islet \u03b2 cell and \u03b1 cell maturation pathways and characteristic features at the single-cell level","volume":"25","author":"Qiu","year":"2017","journal-title":"Cell Metab."},{"key":"2023041407533932900_","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":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"Edger: a bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"547","DOI":"10.1038\/s41587-019-0071-9","article-title":"A comparison of single-cell trajectory inference methods","volume":"37","author":"Saelens","year":"2019","journal-title":"Nat. Biotechnol."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"928","DOI":"10.1016\/j.cell.2019.01.006","article-title":"Optimal-transport analysis of single-cell gene expression identifies developmental trajectories in reprogramming","volume":"176","author":"Schiebinger","year":"2019","journal-title":"Cell"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"7329","DOI":"10.1038\/ncomms8329","article-title":"CD24 tracks divergent pluripotent states in mouse and human cells","volume":"6","author":"Shakiba","year":"2015","journal-title":"Nat. Commun."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"4132","DOI":"10.1016\/j.csbj.2021.07.016","article-title":"Identify differential genes and cell subclusters from time-series scRNA-seq data using scTITANS","volume":"19","author":"Shao","year":"2021","journal-title":"Comput. Struct. Biotechnol. J."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"477","DOI":"10.1186\/s12864-018-4772-0","article-title":"Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics","volume":"19","author":"Street","year":"2018","journal-title":"BMC Genomics"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"e1008205","DOI":"10.1371\/journal.pcbi.1008205","article-title":"Tempora: cell trajectory inference using time-series single-cell RNA sequencing data","volume":"16","author":"Tran","year":"2020","journal-title":"PLoS Comput. Biol."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"391","DOI":"10.1038\/nature18323","article-title":"Dissecting direct reprogramming from fibroblast to neuron using single-cell RNA-seq","volume":"534","author":"Treutlein","year":"2016","journal-title":"Nature"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"692","DOI":"10.1016\/j.ajpath.2013.05.020","article-title":"Differences in the degree of cerulein-induced chronic pancreatitis in C57BL\/6 mouse substrains lead to new insights in identification of potential risk factors in the development of chronic pancreatitis","volume":"183","author":"Ulmasov","year":"2013","journal-title":"Am. J. Pathol."},{"key":"2023041407533932900_","doi-asserted-by":"crossref","DOI":"10.1007\/978-0-387-21706-2","volume-title":"Modern Applied Statistics with S","author":"Venables","year":"2002"},{"key":"2023041407533932900_","doi-asserted-by":"crossref","first-page":"234","DOI":"10.1097\/00006676-200210000-00004","article-title":"Expression of clusterin in human pancreatic cancer","volume":"25","author":"Xie","year":"2002","journal-title":"Pancreas"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/Supplement_1\/i290\/49886767\/btac227.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/Supplement_1\/i290\/49886767\/btac227.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,9,22]],"date-time":"2024-09-22T06:55:42Z","timestamp":1726988142000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/38\/Supplement_1\/i290\/6617492"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,6,24]]},"references-count":32,"journal-issue":{"issue":"Supplement_1","published-print":{"date-parts":[[2022,6,24]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btac227","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2022,7,1]]},"published":{"date-parts":[[2022,6,24]]}}}