单样本的gsea分析图表

最近看到文章 Cancer Cell . 2020 Dec ,标题是《Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation》

链接是:https://pubmed.ncbi.nlm.nih.gov/33450198/  有一个图表很适合布置给学徒作为作业。

研究者对荷兰癌症研究所(NKI)发表的乳腺癌基因表达数据进行分析,绘制出来了如下所示的单样本的gsea分析图表:

单样本的gsea分析图表

图例是:(H) Single-sample gene set enrichment analysis of CTSC correlation to neutrophil-related gene sets in the NKI dataset (n = 295). Scale bar, 50 mm. p values were obtained by two-tailed unpaired t test (B, E, and F) or Pearson’s correlation analysis (others)

我看了看文章:

  • The NKI microarray data and clinical information were obtained from previous publications   (Minn et al., 2007; van de Vijver et al., 2002).

  • The enrichment of neutrophilrelated gene sets in these tumors were analyzed by ssGESA in GenePattern platform with the previously published neutrophil gene sets (Bindea et al., 2013; Ebert et al., 2008; Kim et al., 2019).  These gene sets are listed in Table S6.

  • Download spreadsheet (13KB) ,Table S6. Gene sets used for single-sample gene set enrichment analyses, related to the STAR methods.

数据分析大纲

  • 下载表达量矩阵
  • 读文献拿到3个基因集的基因列表
  • 针对基因列表去表达量矩阵里面走ssGSE分析
  • 对ssGSE分析值与CTSC表达量进行相关性计算
  • 绘制热图

老规矩,学徒作业需要以markdown形式,描述清楚你的分析过程,并且发邮件到我邮箱(jmzeng1314@163.com)

(0)

相关推荐

  • clusterProfiler|GSEA富集分析及可视化

    GSEA(Gene Set EnrichmentAnalysis),即基因集富集分析,无需设定阈值来区分上调下调基因,使用所有的基因进行分析. GO 和 KEGG 可参考:R|clusterProfi ...

  • 关于功能富集分析的基础知识

    富集分析基因富集分析(gene set enrichment analysis)是在一组基因或蛋白中找到一类过表达的基因或蛋白.研究方法可分为三种:Over-Repressentation Analy ...

  • R语言GSEA分析(三)

    GAEA df_all_sort <- df_all[order(df_all$logFC, decreasing = T),]#先按照logFC降序排序 gene_fc = df_all_so ...

  • R语言GSEA分析(一)

    安装并导入要用到的R包 BiocManager::install("clusterProfiler") #感谢Y叔的clusterprofiler包 BiocManager::in ...

  • R语言GSEA分析(二)

    转换基因ID 如基因名是symbol,需要将基因ID转换为Entrez ID格式.Entrez ID实际上是指的Entrez gene ID,是对应于染色体上一个gene location的.每一个发 ...

  • 对于数据可视化分析图表,你了解多少?

    随着可视化技术的快速发展,数据也在逐渐进步.借助于可视化大屏,数据可视化的数据分析图表,将数据以图表的形式进行分析,缩短了用户对于数据的理解时间,同时对于整体的布局也进行了一些美化. 那么,数据可视化 ...

  • GSEA分析一文就够(单机版+R语言版)

    通过前面的讲解,我们顺利的了解了GEO数据库以及如何下载其数据,得到我们想要的表达矩阵,但,这只是分析的开始,最经典的分析就是GSEA了,看看基因全局表达量的变化是否有某些特定的基因集合的倾向性. 历 ...

  • GSEA分析合理性讨论

    g写在前面 昨天发了一个差点捅了马蜂窝的推文:费九牛二虎之力也无法重现的GSEA图 ,当时的确花了两三个小时也没能成功重复文章的GSEA图,觉得很有趣,就跟粉丝分享了,纯粹是交流数据处理技巧. 后来跟 ...

  • 做GSEA分析你的基因到底该如何排序

    大家都知道,GSEA最重要的就是数据集的所有背景基因按照某种指标排序好,这样才能说明你感兴趣的基因集是否在背景基因集里面出现了统计学显著的富集情况,如下: 很明显,我们感兴趣的这个基因集,也就是 hu ...

  • 单基因GSEA分析策略(数据分析免费做活动继续)

    我写了部分常见的TCGA数据库用法: TCGA的28篇教程-免疫全景图 TCGA的28篇教程-指定癌症查看感兴趣基因的表达量 TCGA的28篇教程-对TCGA数据库的任意癌症中任意基因做生存分析 TC ...

  • 生存分析图表美化

    2019中秋节我们的全国巡讲本来应该是去贵州或者兰州,但无奈报名人数太少,临时调到广州.中秋节专场吸引到了几位中山大学临床统计方向博士研究生,其中一位小伙伴把她以前的统计学技能全部R语言化了,会出一系 ...