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.