使用monocle包进行pseudotime分析

欢迎关注”生信修炼手册”!

monocle是一个专门用于分析单细胞转录组数据的R包，提供了聚类，pseudotime, 差异分析等多种功能，该项目的网址如下

https://cole-trapnell-lab.github.io/monocle-release/

本文主要介绍使用该R包进行pseudotime分析的步骤，可以分为以下5步

1. 读取数据

monocle包有很多种读取数据的方式，这里只展示了读取Seurat中的对象的方法，代码如下

# 加载需要的R包
library(Seurat)
library(monocle)
# 设置cell ranger输出结果目录
input_dir <- "/scRNA/outs/filtered_gene_bc_matrices/GRCh38/"
# 读取数据
pbmc.data <- Read10X(data.dir = input_dir)
# 创建Seurat中的对象
pbmc <- CreateSeuratObject(raw.data = pbmc.data, project = "10X")
# 将Seurat中的对象转换为monocle识别的对象
cds <- importCDS(pbmc)

在monocle中，用CellDataSet这种class来存储数据，示意如下