cellchat-与其他单细胞分析工具的接口

与其他单细胞分析工具包的接口

从Seurat或SingleCellExperiment对象创建CellChat对象

从CellChat版本0.5.0开始,USERS可以从Seurat或SingleCellExperiment对象创建新的CellChat对象。。如果输入是Seurat或SingleCellExperiment对象,则默认情况下将使用该对象中的元数据,并且USER必须提供group.by以定义单元格组。例如,group.by =“ ident”用于Seurat对象中的默认单元格标识。

请通过查阅文档中的示例以createCellChat获取详细信息help(createCellChat)

注意:如果用户加载先前计算的CellChat对象(版本<0.5.0),请通过 updateCellChat

在单单元格对象之间转换(Seurat,SingleCellExperiment和andata对象)

有几种出色的工具可在Seurat对象,SingleCellExperiment对象和andata对象之间进行转换。这是两个示例:一个是Satija Lab的Seurat R软件包,另一个是Theis Lab的zellkonverter R软件包

根据从Seurat和Scanpy工具中提取的数据矩阵创建CellChat对象

下面我们展示了如何从其他现有的单细胞分析工具包(包括SeuratScanpy)中将CellChat输入文件提取为数据矩阵。

CellChat中需要的数据输入

CellChat需要两个用户输入:一个是细胞的基因表达数据,另一个是用户分配的细胞标签(即基于标签的模式)或单细胞数据的低维表示形式(即无标签)模式)。对于后者,CellChat通过基于低维空间或伪时间轨迹空间中的单元格距离构建共享的邻居图来自动对单元格进行分组。

资料格式

对于基因表达数据矩阵,基因应在具有行名的行中以及单元格在具有名的列中。需要进行归一化的数据作为CellChat分析的输入,例如,库大小归一化,然后以伪计数1进行对数转换。如果用户提供了计数数据,我们将提供一个normalizeData函数来计算库大小。

对于单元组信息,需要具有行名的数据框作为CellChat的输入。

从Seurat V3对象中提取CellChat输入文件

可以从Seurat对象获得归一化的计数数据和细胞组信息,方法是:

data.input <- GetAssayData(seurat_object, assay = "RNA", slot = "data") # normalized data matrix
labels <- Idents(seurat_object)
meta <- data.frame(group = labels, row.names = names(labels)) # create a dataframe of the cell labels

从Seurat V2对象中提取CellChat输入文件

可以从Seurat对象获得归一化的计数数据和细胞组信息,方法是:

data.input <- seurat_object@data # normalized data matrix
labels <- seurat_object@idents
meta <- data.frame(group = labels, row.names = names(labels)) # create a dataframe of the cell labels

从Scanpy对象中提取CellChat输入文件

anndata提供了一个python类,可用于存储单单元格数据。此数据格式还用于在scanpy软件包中进行存储。我们首先使用reticulate包将数据读入R中,以导入andata模块。

library(reticulate)
ad <- import("anndata", convert = FALSE)
ad_object <- ad$read_h5ad("scanpy_object.h5ad")
# access normalized data matrix
data.input <- t(py_to_r(ad_object$X))
rownames(data.input) <- rownames(py_to_r(ad_object$var))
colnames(data.input) <- rownames(py_to_r(ad_object$obs))
# access meta data
meta.data <- py_to_r(ad_object$obs)
meta <- meta.data

提取所需的CellChat输入文件后,然后创建一个CellChat对象并开始分析。

使用数据矩阵作为输入创建一个CellChat对象

cellchat <- createCellChat(object = data.input, meta = meta, group.by = "labels")

将单元格信息添加到对象的插槽中

如果在创建CellChat对象时未添加单元格mata信息,则USERS也可以稍后使用添加该信息addMeta,并使用设置默认单元格标识setIdent

cellchat <- addMeta(cellchat, meta = meta, meta.name = "labels")
cellchat <- setIdent(cellchat, ident.use = "labels") # set "labels" as default cell identity
levels(cellchat@idents) # show factor levels of the cell labels
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,826评论 6 506
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,968评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,234评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,562评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,611评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,482评论 1 302
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,271评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,166评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,608评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,814评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,926评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,644评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,249评论 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,866评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,991评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,063评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,871评论 2 354

推荐阅读更多精彩内容