首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >CNS图表复现05—免疫细胞亚群再分类

CNS图表复现05—免疫细胞亚群再分类

作者头像
生信技能树jimmy
发布于 2020-09-28 02:48:42
发布于 2020-09-28 02:48:42
2K00
代码可运行
举报
文章被收录于专栏:单细胞天地单细胞天地
运行总次数:0
代码可运行

我们的CNS图表复现之旅已经开始,前面4讲是;

如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。

既然我们根据Cells were defined as non-immune if belonging to a cluster low for PTPRC (gene for CD45)的规则,把细胞区分为了免疫细胞和非免疫细胞,详见:我在单细胞天地的教程:是否是免疫细胞很容易区分那是否是肿瘤细胞呢?

现在我们就可以来复现文章里面的免疫细胞亚群再分类了,如下:

免疫细胞亚群再分类

这幅图超级容易理解,就是13431个免疫细胞,可以继续细分为不同的亚群,如上所述。

首先我们挑选免疫单细胞

因为这个是系列连载教程,所以务必看前面的前面3讲是;

知道我是把全部的2万多个细胞的表达矩阵读入R,变成了Seurat对象,然后判断了细胞是否属于免疫细胞哈。前面的代码核心步骤是:

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
sce@meta.data$immune_annotation <-ifelse(sce@meta.data$RNA_snn_res.0.5  %in% imm ,'immune','non-immune')
# MAke a table 
table(sce@meta.data$immune_annotation)
phe=sce@meta.data
save(phe,file = 'phe-of-immune-or-not.Rdata')

就是把细胞是否属于免疫细胞这个信息(phe=sce@meta.data)保存下来了。后面就可以直接使用啦:

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
rm(list=ls())
options(stringsAsFactors = F)
library(Seurat)
library(ggplot2)
load(file = 'first_sce.Rdata')
sce <- FindClusters(sce, resolution = 0.5)
table(sce@meta.data$RNA_snn_res.0.5)  
load(file = 'phe-of-immune-or-not.Rdata')
table(phe$immune_annotation)
cells.use <- row.names(sce@meta.data)[which(phe$immune_annotation=='immune')]
length(cells.use)
sce <-subset(sce, cells=cells.use)  
sce

继续走Seurat标准流程之聚类分群

这样就拿到了免疫细胞单独进行后续分析,仍然是走Seurat标准流程之聚类分群,代码是:

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
sce
sce <- NormalizeData(sce, normalization.method =  "LogNormalize", 
                     scale.factor = 10000)
GetAssay(sce,assay = "RNA")
sce <- FindVariableFeatures(sce, 
                            selection.method = "vst", nfeatures = 2000) 
sce <- ScaleData(sce) 
sce <- RunPCA(object = sce, pc.genes = VariableFeatures(sce)) 
res.used <- 0.7
sce <- FindClusters(object = sce, verbose = T, resolution = res.used)
set.seed(123)
sce <- RunTSNE(object = sce, dims = 1:15, do.fast = TRUE)
DimPlot(sce,reduction = "tsne",label=T)
DimPlot(sce,reduction = "tsne",label=T, group.by = "patient_id")
table(sce@meta.data$seurat_clusters) 

然后走singleR的自动注释步骤:

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
sce_for_SingleR <- GetAssayData(sce, slot="data")
sce_for_SingleR
library(SingleR)
hpca.se <- HumanPrimaryCellAtlasData()
hpca.se
clusters=sce@meta.data$seurat_clusters
pred.hesc <- SingleR(test = sce_for_SingleR, ref = hpca.se, labels = hpca.se$label.main,
                     method = "cluster", clusters = clusters, 
                     assay.type.test = "logcounts", assay.type.ref = "logcounts")
table(pred.hesc$labels)
celltype = data.frame(ClusterID=rownames(pred.hesc), celltype=pred.hesc$labels, stringsAsFactors = F) 
sce@meta.data$singleR=celltype[match(clusters,celltype$ClusterID),'celltype']
DimPlot(sce, reduction = "tsne", group.by = "singleR")
phe=sce@meta.data
table(phe$singleR)
save(phe,file = 'phe-of-subtypes-Immune-by-singleR.Rdata')

分群结果的对比

但是我们的这个距离分群结果,跟文章是有一定差异的,如下:

分群是:

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
> as.data.frame(sort(table(phe$singleR)))
         Var1 Freq
1 Fibroblasts   91
2          DC  170
3     NK_cell  496
4         CMP  551
5 Neutrophils  645
6      B_cell 1957
7  Macrophage 2104
8    Monocyte 2803
9     T_cells 3698

但是它们仍然是可以细分的,如下:

文章是:

简单总结一下,文章把Macrophage和Monocyte合并起来,把B_cell区分成为两群。而我们的T_cells和NK_cell也需要被合并起来。而且 common myeloid progenitor (CMP) 这群细胞呢,在文章并没有出现,不知道是不是会被命名为 Mast-cells,需要相应的生物学背景。

文章给出的代码是;

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
# Annotate each of the clusters 
free_annotation <- c("T-cells","MF-Monocytes", "MF-Monocytes", "B-cells-PB", "MF-Monocytes", "T-cells", "T-cells", "Neutrophils", "Dendritic", "Mast-cells", "MF-Monocytes", "T-cells", "B-cells-M", "Unknown", "T-cells", "pDCs", "B-cells-M", "MF-Monocytes")
# free_annotation <- c("0","1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18")

但是我们中间很多步骤都没有精准的模仿文章的各个流程和参数,所以这个时候的差异就没办法解决了,我们是22个群,文章是19个群。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2020-09-23,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 单细胞天地 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
暂无评论
推荐阅读
编辑精选文章
换一批
CNS图表复现16—inferCNV结果解读及利用
我们的CNS图表复现之旅已经开始,你可以点击图表复现话题回顾。如果你感兴趣也想加入交流群,可以去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/11/02
13.1K0
CNS图表复现16—inferCNV结果解读及利用
CNS图表复现09—上皮细胞可以区分为恶性与否
如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/10/19
1.6K0
CNS图表复现09—上皮细胞可以区分为恶性与否
CNS图表复现13—使用inferCNV来区分肿瘤细胞的恶性与否
CNS图表复现之旅前面我们已经进行了12讲,你可以点击图表复现话题回顾。如果你感兴趣也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/11/02
10.6K1
CNS图表复现13—使用inferCNV来区分肿瘤细胞的恶性与否
​单细胞专题 | 8.单细胞类型注释之SingleR包详解
单细胞专题 | 1.单细胞测序(10×genomics技术)的原理 单细胞专题 | 2.如何开始单细胞RNASeq数据分析 单细胞专题 | 3.单细胞转录组的上游分析-从BCL到FASTQ 单细胞专题 | 4.单细胞转录组的上游分析-从SRA到FASTQ 单细胞专题 | 5.单细胞转录组的上游分析-从FASTQ到count矩阵 单细胞专题 | 6.单细胞下游分析——不同类型的数据读入 单细胞专题 | 7.单细胞下游分析——常规分析流程案例一
DoubleHelix
2022/12/16
10.1K0
​单细胞专题 | 8.单细胞类型注释之SingleR包详解
CNS图表复现15—inferCNV流程输入数据差异大揭秘
CNS图表复现之旅前面我们已经进行了14讲,你可以点击图表复现话题回顾。如果你感兴趣也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/11/02
2.7K0
CNS图表复现15—inferCNV流程输入数据差异大揭秘
CNS图表复现08—肿瘤单细胞数据第一次分群通用规则
如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/10/19
1.3K0
CNS图表复现08—肿瘤单细胞数据第一次分群通用规则
CNS图表复现03—单细胞区分免疫细胞和肿瘤细胞
交流群里大家讨论的热火朝天,而且也都开始了图表复现之旅,在这里我还是带大家一步步学习CNS图表吧。如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/09/23
2K0
CNS图表复现03—单细胞区分免疫细胞和肿瘤细胞
celldex数据库及singleR自动注释使用
不过一直没有使用singleR自己有的celldex数据包里面的参考数据库,所以这期就一起来看看下载需要的celldex参考数据库,以及使用singleR注释的方法。
生信漫漫学
2024/07/26
1.2K0
celldex数据库及singleR自动注释使用
CNS图表复现18—细胞亚群的比例展示
在前面的教程 CNS图表复现05—免疫细胞亚群再分类 ,我提到到免疫细胞通常是以CD45阳性为标志,第一次分群规则是 :
生信技能树jimmy
2020/12/11
2.3K0
CNS图表复现18—细胞亚群的比例展示
singleR的7个数据库文件下载失败的解决方案
singleR自带7个数据库文件,需要联网才能下载,其中5个是人类数据,2个是小鼠的数据:
生信技能树jimmy
2020/09/04
15.1K6
singleR的7个数据库文件下载失败的解决方案
​harmony、不harmony,这是个问题
harmony应用主成分分析,将转录组表达谱嵌入到低维空间中,然后应用迭代过程去除数据集特有的影响
生信技能树jimmy
2023/08/31
1.5K0
​harmony、不harmony,这是个问题
比较copykat和infercnv这两个从单细胞转录组数据推断肿瘤拷贝数变异技术差异
如果是现在才看到这个系列的小伙伴建议自己去先读一下CELL杂志的文章:Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing ,因为作者提供了全套代码,在:https://github.com/czbiohub/scell_lung_adenocarcinoma ,研究者们共收集到30位患者的49份活检样本(biopsy),分为三种类型:治疗前(TKI naive [TN]),靶向治疗后肿瘤消退或稳定(RD, residual disease state)以及靶向治疗后肿瘤仍然增长(PD, upon subsequent progressive disease),这样单细胞转录组数据就非常丰富!
生信技能树
2021/12/04
3.8K0
比较copykat和infercnv这两个从单细胞转录组数据推断肿瘤拷贝数变异技术差异
这个肝癌单细胞数据集居然真的有这么多NK细胞亚群啊!!!
根据生信技能树发布的学徒作业:SingleR说是NK细胞你就相信了吗, 验证一下看真的是有这么多NK细胞
生信技能树
2023/09/04
4720
这个肝癌单细胞数据集居然真的有这么多NK细胞亚群啊!!!
使用SingleR对小鼠免疫单细胞自动注释并不可靠
让我想起来了被免疫细胞里面的淋巴系和髓系的细胞细分群及注释支配的恐惧,虽然在前面的例子:人人都能学会的单细胞聚类分群注释 ,我们演示了第一层次的分群,但是第二层次的分群就很麻烦。绝大部分文献基本上没有参考价值,比如发表于2020年8月的文章是:《Single-cell RNA sequencing uncovers heterogenous transcriptional signatures in macrophages during efferocytosis》,链接在:https://www.nature.com/articles/s41598-020-70353-y ,做了6只 C57BL/6J品系的小鼠的单细胞,每个小鼠平均1400个单细胞,合起来还不到一万个单细胞,让我们来来看看作者是如何对细胞亚群进行注释:
生信菜鸟团
2021/04/13
2K0
使用SingleR对小鼠免疫单细胞自动注释并不可靠
CNS图表复现20—第三次分群,以T细胞为例
前面我们展现了 CNS图表复现08—肿瘤单细胞数据第一次分群通用规则,然后呢,第二次分群的上皮细胞可以细分恶性与否,免疫细胞呢,细分可以成为: B细胞,T细胞,巨噬细胞,树突细胞等等。实际上每个免疫细胞亚群仍然可以继续精细的划分,以文章为例:
生信技能树jimmy
2020/12/11
1K0
CNS图表复现20—第三次分群,以T细胞为例
CNS图表复现12—检查原文的细胞亚群的标记基因
CNS图表复现之旅前面我们已经进行了11讲,你可以点击图表复现话题回顾。如果你感兴趣也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/11/02
1.2K0
CNS图表复现12—检查原文的细胞亚群的标记基因
百万细胞舍我其谁(一晚上解决战斗)
而且,在为了,百万级别的单细胞转录组数据集会越来越多,首先是现在的10x单细胞转录组一个样品普遍只需要不到1万的人民币,而且国产是半价,只需要区区100万人民币就可以拿到100万甚至更多的细胞数量。这个经费对于张泽民院士来说,那就是毛毛雨了,他老人家的绝大部分cns文章都是标配百万级别单细胞数量。而且现在越来越多的多个单细胞数据整合文章,反正是公共数据集,多整合一些就很容易达到百万的数据量!
生信菜鸟团
2024/12/30
1610
百万细胞舍我其谁(一晚上解决战斗)
CNS图表复现06—根据CellMarker网站进行人工校验免疫细胞亚群
如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
生信技能树jimmy
2020/09/28
9620
CNS图表复现06—根据CellMarker网站进行人工校验免疫细胞亚群
很多时候你就是不知道如何提问
其中有一个环节是需要比较seurat分群以及singleR的分群,这样就可以合理的命名啦。
生信技能树
2020/09/14
4750
很多时候你就是不知道如何提问
使用singleR基于自建数据库来自动化注释单细胞转录组亚群
但是因为singleR的数据库资源陈旧而且很有限,满足不了日益增长的单细胞应用,后面我们都是主推第一层次降维聚类分群后的人工命名,通常我们拿到了肿瘤相关的单细胞转录组的表达量矩阵后的第一层次降维聚类分群通常是:
生信技能树
2024/04/13
1.3K0
使用singleR基于自建数据库来自动化注释单细胞转录组亚群
推荐阅读
相关推荐
CNS图表复现16—inferCNV结果解读及利用
更多 >
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档