前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >单细胞数据(scRNAseq)可以做GSEA吗?

单细胞数据(scRNAseq)可以做GSEA吗?

作者头像
生信编程日常
发布2020-09-07 09:58:55
3.8K0
发布2020-09-07 09:58:55
举报
文章被收录于专栏:生物信息学、python、R、linux

单细胞测序数据也可以做gsea,步骤跟用RNAseq的数据差不多,主要是要用到差异基因并且根据Fold change来排序。

代码语言:javascript
复制
library(msigdbr)
library(fgsea)
library(dplyr)
library(ggplot2)

选择自己数据的物种以及要做的GSEA的数据库类型

代码语言:javascript
复制
##查看物种的数据 msigdbr_show_species()
m_df<- msigdbr(species = "Homo sapiens", category = "H")
代码语言:javascript
复制
#将m_df的基因与通路取出并改成一个通路对应相应基因的格式
fgsea_sets<- m_df %>% split(x = .$gene_symbol, f = .$gs_name) 

这里选用seurat求差异基因,并将差异基因按显著性排序

代码语言:javascript
复制
#每一个细胞类型的GSEA按显著性进行降序排序
gesa_TvsC_allgenes<-FindMarkers(seurat.combined, ident.1 = "TCF7KO", ident.2 = "NC", verbose = FALSE,test.use ="roc",logfc.threshold = 0.01,only.pos =F)
gesa_TvsC_allgenes$gene<-rownames(gesa_TvsC_allgenes)

gsea_genes<-gesa_TvsC_allgenes %>%
  arrange(desc(myAUC), desc(avg_diff)) %>%
dplyr::select(gene,avg_diff)
代码语言:javascript
复制
library(tibble)
ranks <- deframe(gsea_genes)
代码语言:javascript
复制
fgseaRes <- fgsea(pathways = fgsea_sets,
                  stats = ranks ,
                  minSize=5,
                  maxSize=500,
                  nperm=10000)

对结果进行整理

代码语言:javascript
复制
library(dplyr)
fgseaResTidy <- fgseaRes %>%
  as_tibble() %>%
  arrange(desc(NES))

fgseaResTidy %>%
  dplyr::select(-leadingEdge, -ES, -nMoreExtreme) %>%
  arrange(padj) %>%
  head()
代码语言:javascript
复制
ggplot(fgseaResTidy %>% filter(padj < 0.5) %>% head(n= 15), aes(reorder(pathway, NES), NES)) +
  geom_col(aes(fill= NES < 0)) +
  coord_flip() +
  labs(x="Pathway", y="Normalized Enrichment Score",
       title="Hallmark pathways NES from GSEA") +
  theme_minimal() 
代码语言:javascript
复制
plotEnrichment(fgsea_sets[["HALLMARK_APICAL_JUNCTION"]],
               ranks) + labs(title="HALLMARK_APICAL_JUNCTION")

参考:https://jishuin.proginn.com/p/763bfbd24004

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

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

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

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
相关产品与服务
数据库
云数据库为企业提供了完善的关系型数据库、非关系型数据库、分析型数据库和数据库生态工具。您可以通过产品选择和组合搭建,轻松实现高可靠、高可用性、高性能等数据库需求。云数据库服务也可大幅减少您的运维工作量,更专注于业务发展,让企业一站式享受数据上云及分布式架构的技术红利!
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档