Abstract

BackgroundSingle cell RNA sequencing (scRNA-seq) brings unprecedented opportunities for mapping the heterogeneity of complex cellular environments such as bone marrow, and provides insight into many cellular processes. Single cell RNA-seq, however, has a far larger fraction of missing data reported as zeros (dropouts) than traditional bulk RNA-seq. This makes difficult not only the clustering of cells, but also the assignment of the resulting clusters into predefined cell types based on known molecular signatures, such as the expression of characteristic cell surface markers. ResultsWe present a computational tool for processing single cell RNA-seq data that uses a voting algorithm to identify cells based on approval votes received by known molecular markers. Using a stochastic procedure that accounts for biases due to dropout errors and imbalances in the number of known molecular signatures for different cell types, the method computes the statistical significance of the final approval score and automatically assigns a cell type to clusters without an expert curator. We demonstrate the utility of the tool in the analysis of eight samples of bone marrow from the Human Cell Atlas. The tool provides a systematic identification of cell types in bone marrow based on a recently-published manually-curated cell marker database [1], and incorporates a suite of visualization tools that can be overlaid on a t-SNE representation. The software is freely available as a python package at https://github.com/sdomanskyi/DigitalCellSorter ConclusionsThis methodology assures that extensive marker to cell type matching information is taken into account in a systematic way when assigning cell clusters to cell types. Moreover, the method allows for a high throughput processing of multiple scRNA-seq datasets, since it does not involve an expert curator, and it can be applied recursively to obtain cell sub-types. The software is designed to allow the user to substitute the marker to cell type matching information and apply the methodology to different cellular environments.