Exponential-Family Embedding With Application to Cell Developmental Trajectories for Single-Cell RNA-Seq Data
Year: 2021
Type: article
Abstract: Scientists often embed cells into a lower-dimensional space when studying single-cell RNA-seq data for improved downstream analyses such as developmental trajectory analyses, but the statistical properties of such nonlinear embedding methods are often not well understood. In this article, we develop the exponential-family SVD (eSVD), a nonlinear embedding method for both cells and genes jointly with respect to a random dot product model using exponential-family distributions. Our estimator uses alternating minimization, which enables us to have a computationally efficient method, prove the identifiability conditions and consistency of our method, and provide statistically principled procedures to tune our method. All these qualities help advance the single-cell embedding literature, and we provide extensive simulations to demonstrate that the eSVD is competitive compared to other embedding methods. We apply the eSVD via Gaussian distributions where the standard deviations are proportional to the means to analyze a single-cell dataset of oligodendrocytes in mouse brains. Using the eSVD estimated embedding, we then investigate the cell developmental trajectories of the oligodendrocytes. While previous results are not able to distinguish the trajectories among the mature oligodendrocyte cell types, our diagnostics and results demonstrate there are two major developmental trajectories that diverge at mature oligodendrocytes. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplementary materials.
Authors: Kevin Z. Lin, Jing Lei, Kathryn Roeder
Institutions: University of Pennsylvania, Carnegie Mellon University
Language: English
Subfield: Molecular Biology
Domain: Life Sciences
Open Access status: green
Awards: G1363864812, Estimation and Inference of Gene Regulatory Networks, G202538826, Genetic Association and Personalized Medicine, G3297477028, CAREER: Modernizing Classical Nonparametric and Multivariate Theory for Large-scale, High-dimensional Data Analysis, Theory and Methods for Large-Scale Multi-Modal Matrix Data, Genetics Association in Schizophrenia and Other Disorders
OpenAlex ID: W3132119066