{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,18]],"date-time":"2026-07-18T07:27:37Z","timestamp":1784359657159,"version":"3.55.0"},"reference-count":51,"publisher":"Oxford University Press (OUP)","issue":"8","license":[{"start":{"date-parts":[[2023,7,27]],"date-time":"2023-07-27T00:00:00Z","timestamp":1690416000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,8,1]]},"abstract":"Abstract<\/jats:title>\n \n Motivation<\/jats:title>\n Single-cell RNA-sequencing (scRNA-seq) has enabled the molecular profiling of thousands to millions of cells simultaneously in biologically heterogenous samples. Currently, the common practice in scRNA-seq is to determine cell type labels through unsupervised clustering and the examination of cluster-specific genes. However, even small differences in analysis and parameter choosing can greatly alter clustering results and thus impose great influence on which cell types are identified. Existing methods largely focus on determining the optimal number of robust clusters, which can be problematic for identifying cells of extremely low abundance due to their subtle contributions toward overall patterns of gene expression.<\/jats:p>\n <\/jats:sec>\n \n Results<\/jats:title>\n Here, we present a carefully designed framework, SCISSORS, which accurately profiles subclusters within broad cluster(s) for the identification of rare cell types in scRNA-seq data. SCISSORS employs silhouette scoring for the estimation of heterogeneity of clusters and reveals rare cells in heterogenous clusters by a multi-step semi-supervised reclustering process. Additionally, SCISSORS provides a method for the identification of marker genes of high specificity to the cell type. SCISSORS is wrapped around the popular Seurat R package and can be easily integrated into existing Seurat pipelines.<\/jats:p>\n <\/jats:sec>\n \n Availability and implementation<\/jats:title>\n SCISSORS, including source code and vignettes, are freely available at https:\/\/github.com\/jr-leary7\/SCISSORS.<\/jats:p>\n <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad449","type":"journal-article","created":{"date-parts":[[2023,7,27]],"date-time":"2023-07-27T15:12:52Z","timestamp":1690470772000},"source":"Crossref","is-referenced-by-count":24,"title":["S<\/u>ub-C<\/u>luster I<\/u>dentification through S<\/u>emi-S<\/u>upervised O<\/u>ptimization of R<\/u>are-Cell S<\/u>ilhouettes (SCISSORS) in single-cell RNA-sequencing"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0009-0004-8821-3269","authenticated-orcid":false,"given":"Jack R","family":"Leary","sequence":"first","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Biostatistics, University of Florida , Gainesville, FL 32603, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yi","family":"Xu","sequence":"additional","affiliation":[{"name":"Department of Pharmacology, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ashley B","family":"Morrison","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Chong","family":"Jin","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Emily C","family":"Shen","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Peyton C","family":"Kuhlers","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Biostatistics, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ye","family":"Su","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Naim U","family":"Rashid","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Biostatistics, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jen Jen","family":"Yeh","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Pharmacology, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Surgery, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8684-5762","authenticated-orcid":false,"given":"Xianlu Laura","family":"Peng","sequence":"additional","affiliation":[{"name":"Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"},{"name":"Department of Pharmacology, University of North Carolina at Chapel Hill , Chapel Hill, NC 27599, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2023,7,27]]},"reference":[{"key":"2023081001055488500_btad449-B1","doi-asserted-by":"crossref","first-page":"420","DOI":"10.1007\/3-540-44503-X_27","volume-title":"Database theory\u2014ICDT 2001","author":"Aggarwal","year":"2001"},{"key":"2023081001055488500_btad449-B2","doi-asserted-by":"crossref","first-page":"114","DOI":"10.1016\/j.mam.2017.07.002","article-title":"Identifying cell populations with scRNASeq","volume":"59","author":"Andrews","year":"2018","journal-title":"Mol Aspects Med"},{"key":"2023081001055488500_btad449-B3","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":"2023081001055488500_btad449-B4","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1186\/s13059-016-0927-y","article-title":"Design and computational analysis of single-cell RNA-sequencing experiments","volume":"17","author":"Bacher","year":"2016","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B5","doi-asserted-by":"crossref","first-page":"47","DOI":"10.1038\/nature16965","article-title":"Genomic analyses identify molecular subtypes of pancreatic cancer","volume":"531","author":"Bailey","year":"2016","journal-title":"Nature"},{"key":"2023081001055488500_btad449-B6","first-page":"1345","article-title":"Clustering on the unit hypersphere using von Mises-Fisher distributions","volume":"6","author":"Banerjee","year":"2005","journal-title":"J Mach Learn Res"},{"key":"2023081001055488500_btad449-B7","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1016\/j.cels.2016.08.011","article-title":"A single-cell transcriptomic map of the human and mouse pancreas reveals inter- and intra-cell population structure","volume":"3","author":"Baron","year":"2016","journal-title":"Cell Syst"},{"key":"2023081001055488500_btad449-B8","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":"J Stat Mech"},{"key":"2023081001055488500_btad449-B9","doi-asserted-by":"crossref","first-page":"5","DOI":"10.1016\/0167-9473(96)88919-5","article-title":"Probabilistic models in cluster analysis","volume":"23","author":"Bock","year":"1996","journal-title":"Comput Stat Data Anal"},{"key":"2023081001055488500_btad449-B10","doi-asserted-by":"crossref","first-page":"231","DOI":"10.1038\/s41588-019-0566-9","article-title":"Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution","volume":"52","author":"Chan-Seng-Yue","year":"2020","journal-title":"Nat Genet"},{"key":"2023081001055488500_btad449-B11","doi-asserted-by":"crossref","first-page":"giz121","DOI":"10.1093\/gigascience\/giz121","article-title":"IKAP-Identifying K mAjor cell Population groups in single-cell RNA-sequencing analysis","volume":"8","author":"Chen","year":"2019","journal-title":"Gigascience"},{"key":"2023081001055488500_btad449-B12","doi-asserted-by":"crossref","first-page":"500","DOI":"10.1038\/nm.2344","article-title":"Subtypes of pancreatic ductal adenocarcinoma and their differing responses to therapy","volume":"17","author":"Collisson","year":"2011","journal-title":"Nat Med"},{"key":"2023081001055488500_btad449-B13","doi-asserted-by":"crossref","first-page":"158","DOI":"10.1186\/s12859-020-3482-1","article-title":"GiniClust3: a fast and memory-efficient tool for rare cell type identification","volume":"21","author":"Dong","year":"2020","journal-title":"BMC Bioinformatics"},{"key":"2023081001055488500_btad449-B14","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":"Duo","year":"2018","journal-title":"F1000Res"},{"key":"2023081001055488500_btad449-B15","doi-asserted-by":"crossref","first-page":"1102","DOI":"10.1158\/2159-8290.CD-19-0094","article-title":"Cross-species single-cell analysis of pancreatic ductal adenocarcinoma reveals antigen-presenting cancer-associated fibroblasts","volume":"9","author":"Elyada","year":"2019","journal-title":"Cancer Discov"},{"key":"2023081001055488500_btad449-B16","first-page":"226","author":"Ester","year":"1996"},{"key":"2023081001055488500_btad449-B17","doi-asserted-by":"crossref","first-page":"e94","DOI":"10.1093\/nar\/gkaa582","article-title":"Variance-adjusted Mahalanobis (VAM): a fast and accurate method for cell-specific gene set scoring","volume":"48","author":"Frost","year":"2020","journal-title":"Nucleic Acids Res"},{"key":"2023081001055488500_btad449-B18","doi-asserted-by":"crossref","first-page":"251","DOI":"10.1038\/nature14966","article-title":"Single-cell messenger RNA sequencing reveals rare intestinal cell types","volume":"525","author":"Gr\u00fcn","year":"2015","journal-title":"Nature"},{"key":"2023081001055488500_btad449-B19","doi-asserted-by":"crossref","first-page":"296","DOI":"10.1186\/s13059-019-1874-1","article-title":"Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression","volume":"20","author":"Hafemeister","year":"2019","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B20","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":"2023081001055488500_btad449-B21","first-page":"100","article-title":"Algorithm as 136: a K-means clustering algorithm","volume":"28","author":"Hartigan","year":"1979","journal-title":"J R Stat Soc"},{"key":"2023081001055488500_btad449-B22","doi-asserted-by":"crossref","first-page":"193","DOI":"10.1007\/BF01908075","article-title":"Comparing partitions","volume":"2","author":"Hubert","year":"1985","journal-title":"J Classif"},{"key":"2023081001055488500_btad449-B23","doi-asserted-by":"crossref","first-page":"1159","DOI":"10.1093\/bioinformatics\/btz704","article-title":"Spectrum: fast density-aware spectral clustering for single and multi-omic data","volume":"36","author":"John","year":"2020","journal-title":"Bioinformatics"},{"key":"2023081001055488500_btad449-B24","doi-asserted-by":"crossref","first-page":"12214","DOI":"10.1096\/fj.202000363RR","article-title":"Basal-like and classical cells coexist in pancreatic cancer revealed by single-cell analysis on biopsy-derived pancreatic cancer organoids from the classical subtype","volume":"34","author":"Juiz","year":"2020","journal-title":"FASEB J"},{"key":"2023081001055488500_btad449-B25","doi-asserted-by":"crossref","first-page":"2035","DOI":"10.3389\/fimmu.2019.02035","article-title":"Human monocyte subsets and phenotypes in major chronic inflammatory diseases","volume":"10","author":"Kapellos","year":"2019","journal-title":"Front Immunol"},{"key":"2023081001055488500_btad449-B26","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":"2023081001055488500_btad449-B27","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1186\/s13059-020-1926-6","article-title":"Eleven grand challenges in single-cell data science","volume":"21","author":"L\u00e4hnemann","year":"2020","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B28","doi-asserted-by":"crossref","first-page":"3155","DOI":"10.1038\/s41467-020-16904-3","article-title":"An entropy-based metric for assessing the purity of single cell populations","volume":"11","author":"Liu","year":"2020","journal-title":"Nat Commun"},{"key":"2023081001055488500_btad449-B29","doi-asserted-by":"crossref","first-page":"232","DOI":"10.1186\/s13059-021-02445-5","article-title":"MultiK: an automated tool to determine optimal cluster numbers in single-cell RNA sequencing data","volume":"22","author":"Liu","year":"2021","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B30","doi-asserted-by":"crossref","first-page":"1168","DOI":"10.1038\/ng.3398","article-title":"Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma","volume":"47","author":"Moffitt","year":"2015","journal-title":"Nat Genet"},{"key":"2023081001055488500_btad449-B31","doi-asserted-by":"crossref","first-page":"333","DOI":"10.1038\/s41587-019-0392-8","article-title":"Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas","volume":"38","author":"Moncada","year":"2020","journal-title":"Nat Biotechnol"},{"key":"2023081001055488500_btad449-B32","doi-asserted-by":"crossref","first-page":"3217","DOI":"10.1093\/bioinformatics\/bty316","article-title":"CONICS integrates scRNA-seq with DNA sequencing to map gene expression to tumor Sub-clones","volume":"34","author":"M\u00fcller","year":"2018","journal-title":"Bioinformatics"},{"key":"2023081001055488500_btad449-B33","doi-asserted-by":"crossref","first-page":"579","DOI":"10.1084\/jem.20162024","article-title":"Distinct populations of inflammatory fibroblasts and myofibroblasts in pancreatic cancer","volume":"214","author":"\u00d6hlund","year":"2017","journal-title":"J Exp Med"},{"key":"2023081001055488500_btad449-B34","doi-asserted-by":"crossref","first-page":"4901","DOI":"10.1158\/1078-0432.CCR-19-3724","article-title":"GATA6 expression distinguishes classical and basal-like subtypes in advanced pancreatic cancer","volume":"26","author":"O'Kane","year":"2020","journal-title":"Clin Cancer Res"},{"key":"2023081001055488500_btad449-B35","doi-asserted-by":"crossref","first-page":"1761","DOI":"10.3389\/fimmu.2019.01761","article-title":"A novel, five-marker alternative to CD16-CD14 gating to identify the three human monocyte subsets","volume":"10","author":"Ong","year":"2019","journal-title":"Front Immunol"},{"key":"2023081001055488500_btad449-B36","first-page":"2825","article-title":"Scikit-learn: machine learning in Python","volume":"12","author":"Pedregosa","year":"2011","journal-title":"JMLR"},{"key":"2023081001055488500_btad449-B37","doi-asserted-by":"crossref","first-page":"4729","DOI":"10.1038\/s41467-019-12517-7","article-title":"De novo compartment deconvolution and weight estimation of tumor samples using DECODER","volume":"10","author":"Peng","year":"2019","journal-title":"Nat Commun"},{"key":"2023081001055488500_btad449-B38","doi-asserted-by":"crossref","first-page":"135","DOI":"10.1214\/aos\/1176345339","article-title":"Strong consistency of K-means clustering","volume":"9","author":"Pollard","year":"1981","journal-title":"Ann Statist"},{"key":"2023081001055488500_btad449-B39","doi-asserted-by":"crossref","first-page":"185","DOI":"10.1016\/j.ccell.2017.07.007","article-title":"Integrated genomic characterization of pancreatic ductal adenocarcinoma","volume":"32","author":"Raphael","year":"2017","journal-title":"Cancer Cell"},{"key":"2023081001055488500_btad449-B40","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1158\/1078-0432.CCR-19-1467","article-title":"Purity Independent Subtyping of Tumors (PurIST), a clinically robust, single-sample classifier for tumor subtyping in pancreatic cancer","volume":"26","author":"Rashid","year":"2020","journal-title":"Clin Cancer Res"},{"key":"2023081001055488500_btad449-B41","doi-asserted-by":"crossref","first-page":"53","DOI":"10.1016\/0377-0427(87)90125-7","article-title":"Silhouettes: a graphical aid to the interpretation and validation of cluster analysis","volume":"20","author":"Rousseeuw","year":"1987","journal-title":"J Comput Appl Math"},{"key":"2023081001055488500_btad449-B42","author":"Satija","year":"2019"},{"key":"2023081001055488500_btad449-B43","doi-asserted-by":"crossref","first-page":"1888","DOI":"10.1016\/j.cell.2019.05.031","article-title":"Comprehensive integration of single-cell data","volume":"177","author":"Stuart","year":"2019","journal-title":"Cell"},{"key":"2023081001055488500_btad449-B44","doi-asserted-by":"crossref","first-page":"377","DOI":"10.1038\/nmeth.1315","article-title":"mRNA-Seq whole-transcriptome analysis of a single cell","volume":"6","author":"Tang","year":"2009","journal-title":"Nat Methods"},{"key":"2023081001055488500_btad449-B45","doi-asserted-by":"crossref","first-page":"5233","DOI":"10.1038\/s41598-019-41695-z","article-title":"From Louvain to Leiden: guaranteeing well-connected communities","volume":"9","author":"Traag","year":"2019","journal-title":"Sci Rep"},{"key":"2023081001055488500_btad449-B46","doi-asserted-by":"crossref","first-page":"236","DOI":"10.1080\/01621459.1963.10500845","article-title":"Hierarchical grouping to optimize an objective function","volume":"58","author":"Ward","year":"1963","journal-title":"J Am Stat Assoc"},{"key":"2023081001055488500_btad449-B47","doi-asserted-by":"crossref","first-page":"142","DOI":"10.1186\/s13059-019-1739-7","article-title":"CellSIUS provides sensitive and specific detection of rare cell populations from complex single-cell RNA-seq data","volume":"20","author":"Wegmann","year":"2019","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B48","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1186\/s13059-017-1382-0","article-title":"SCANPY: large-scale single-cell gene expression data analysis","volume":"19","author":"Wolf","year":"2018","journal-title":"Genome Biol"},{"key":"2023081001055488500_btad449-B49","doi-asserted-by":"crossref","first-page":"1269","DOI":"10.1093\/bioinformatics\/bty793","article-title":"SAFE-clustering: single-cell aggregated (from ensemble) clustering for single-cell RNA-seq data","volume":"35","author":"Yang","year":"2019","journal-title":"Bioinformatics"},{"key":"2023081001055488500_btad449-B50","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":"2023081001055488500_btad449-B51","doi-asserted-by":"crossref","first-page":"1317","DOI":"10.1016\/j.immuni.2019.03.009","article-title":"Single-cell transcriptomics of human and mouse lung cancers reveals conserved myeloid populations across individuals and species","volume":"50","author":"Zilionis","year":"2019","journal-title":"Immunity"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad449\/50975968\/btad449.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/8\/btad449\/51078655\/btad449.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/8\/btad449\/51078655\/btad449.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,8,10]],"date-time":"2023-08-10T01:06:52Z","timestamp":1691629612000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad449\/7232228"}},"subtitle":[],"editor":[{"given":"Christina","family":"Kendziorski","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,7,27]]},"references-count":51,"journal-issue":{"issue":"8","published-print":{"date-parts":[[2023,8,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad449","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,8,1]]},"published":{"date-parts":[[2023,7,27]]},"article-number":"btad449"}}