Statistical Methods for Single-Cell RNA-Seq

单细胞rna测序（scRNA-seq）是最近出现的一种强大的技术，用于研究个体细胞水平上的转录组变异和调控。传统的“散装”RNA-seq从大量细胞中收集RNA，并测量样本中的平均表达。相反，scRNA-seq揭示了细胞间的异质性，为理解发育、分化和疾病病因的生物学过程提供了关键信息。这项新技术在基础研究和临床研究方面的应用不断扩大，但也以其独特的数据特性给分析带来了挑战。这些包括：1)由于少量起始材料和额外的样品制备程序，在存在技术伪影的情况下难以估计分子计数；2)单细胞RNA计数数据缺乏适当的功能聚类方法，这些数据比批量RNA-seq要稀疏得多；3)缺乏对异质性的定量测量和比较。我们建议通过开发一系列用于scRNA-seq数据预处理和分析的新颖统计方法来解决这些挑战。这包括消除RNA捕获和扩增的技术偏差，以获得准确的分子水平计数，识别细胞的功能类型/亚型和可解释的特征群，解释样品和细胞之间的异质性，以及识别差异异质性。在这个项目中开发的所有方法都将作为免费的开源软件实施和发布，以使基因组学研究社区受益。本文建立的概率模型和统计框架将为单细胞ATAC-seq或BS-seq等其他单细胞测序实验的方法学发展奠定基础。

Single-cell RNA-sequencing (scRNA-seq) has emerged very recently as a powerful technology to investigate transcriptomic variation and regulation at the individual cell level. Traditional "bulk" RNA-seq pools RNA from a large number of cells and measures the averaged expressions in a sample. In contrast, scRNA-seq reveals cell to cell heterogeneity, providing critical information to the understanding of biological processes in development, differentiation, and disease etiologies. This new technology leads to an expansion of applications in both basic and clinical research, but also brings challenges in analysis with its unique data characteristics. These include: 1) difficulty in estimating molecule counts with the presence of technical artifacts, due to small amount of starting material and additional sample preparation procedures; 2) lack of appropriate methods for functional clustering for single cell RNA count data, which are much sparser than bulk RNA-seq; 3) lack of a quantitative measure and comparison of heterogeneity. We propose to address these challenges by developing a series of novel statistical methods for scRNA-seq data preprocessing and analyses. This includes removing technical bias in RNA capture and amplification to obtain accurate molecule level counts, identifying functional types/ subtypes of cells and interpretable feature groups, explaining heterogeneities between samples and cells, and identifying differential heterogeneity. All methods developed in this project will be implemented and released as free, open source software to benefit the genomics research community. The probability model and statistical framework established in this proposal will lay a foundation for future methodology development for other single cell sequencing experiments such as single-cell ATAC-seq or BS-seq.