Abstract:Single cell RNA sequencing (scRNA-seq) is a gene expression profiling technique that is presently revolutionizing the study of complex cellular systems in the biological sciences. Existing scRNA-seq methods suffer from sub-optimal target recovery leading to inaccurate measurements including many false negatives. The resulting 'zero-inflated' data may confound data interpretation and visualization. Since cells have coherent phenotypes defined by conserved molecular circuitries (i.e. multiple gene products working together) and since similar cells utilize similar circuits, information about each expression value or 'node' in a multi-cell, multi-gene scRNA-seq data set is expected to also be predictable from other nodes in the data set. Based on this logic, several approaches have been proposed to impute missing values in a data set by extracting information from its non-zero measurements. In this study, we apply non-negative matrix factorization to a selection of published scRNA-seq data sets followed by multiplication of the factor matrices to generate idealized 'completed' model versions of the data. From the model matrices, we recommend new values where original measurements are likely to be inaccurate and where 'zero' measurements are predicted to be false negatives. The resulting imputed data model predicts novel type markers and expression patterns that match orthogonal measurements and field literature better than those obtained from pre-imputation data or alternative imputation strategies.

Factorization-based Imputation of Expression in Single-cell Transcriptomic Analysis (FIESTA) recovers Gene-Cell-State Relationships