指数族嵌入及其在单细胞RNA测序数据细胞发育轨迹中的应用

Exponential-Family Embedding With Application to Cell Developmental Trajectories for Single-Cell RNA-Seq Data

Journal of the American Statistical Association · 2021
被引 16
ABS 4

中文导读

提出指数族SVD(eSVD)非线性嵌入方法,同时嵌入细胞和基因,基于随机点积模型和指数族分布,通过交替最小化实现高效估计,证明可识别性和一致性,应用于小鼠脑少突胶质细胞数据,揭示成熟少突胶质细胞的两条主要发育轨迹。

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.

单细胞RNA测序嵌入方法细胞发育轨迹指数族分布统计估计