GSEA富集分析:从概念理解到界面实操

共 2246字,需浏览 5分钟

 ·

2020-05-21 23:22

GSEA定义

Gene Set Enrichment Analysis (基因集富集分析)用来评估一个预先定义的基因集的基因在与表型相关度排序的基因表中的分布趋势,从而判断其对表型的贡献。


输入数据包含两部分:

  • 一是已知功能的基因集 (可以是GO注释、MsigDB的注释或其它符合格式的基因集定义);

  • 一是表达矩阵,软件会对基因根据其于表型的关联度(可以理解为表达值的变化)从大到小排序,然后判断基因集内每条注释下的基因是否富集于表型相关度排序后基因表的上部或下部,从而判断此基因集内基因的协同变化对表型变化的影响。


(The gene sets are defined based on prior biological knowledge, e.g., published information about biochemical pathways or coexpression in previous experiments. The goal of GSEA is to determine whether members of a gene set S tend to occur toward the top (or bottom) of the listL, in which case the gene set is correlated with the phenotypic class distinction.)


这与之前讲述的GO富集分析不同。GO富集分析是先筛选差异基因,再判断差异基因在哪些注释的通路存在富集;这涉及到阈值的设定,存在一定主观性并且只能用于表达变化较大的基因,即我们定义的显著差异基因。


而GSEA则不局限于差异基因,从基因集的富集角度出发,理论上更容易囊括细微但协调性的变化对生物通路的影响。

GSEA原理

给定一个排序的基因表L和一个预先定义的基因集S (比如编码某个代谢通路的产物的基因, 基因组上物理位置相近的基因,或同一GO注释下的基因)。


GSEA的目的是判断S里面的成员sL里面是随机分布还是主要聚集在L的顶部或底部。


这些基因排序的依据是其在不同表型状态下的表达差异,若研究的基因集S的成员显著聚集在L的顶部或底部,则说明此基因集成员对表型的差异有贡献,也是我们关注的基因集。


aa89045eea424f127d595de9f1bffd6a.webp


GSEA计算中几个关键概念:


  1. 计算富集得分 (ES, enrichment score). ES反应基因集成员s在排序列表L的两端富集的程度。计算方式是,从基因集L的第一个基因开始,计算一个累计统计值。当遇到一个落在s里面的基因,则增加统计值。遇到一个不在s里面的基因,则降低统计值。每一步统计值增加或减少的幅度与基因的表达变化程度(更严格的是与基因和表型的关联度)是相关的。富集得分ES最后定义为最大的峰值。正值ES表示基因集在列表的顶部富集,负值ES表示基因集在列表的底部富集。

  2.  评估富集得分(ES)的显著性。通过基于表型而不改变基因之间关系的排列检验 (permutation test)计算观察到的富集得分(ES)出现的可能性。若样品量少,也可基于基因集做排列检验 (permutation test),计算p-value。

  3.  多重假设检验矫正。首先对每个基因子集s计算得到的ES根据基因集的大小进行标准化得到Normalized Enrichment Score (NES)。随后针对NES计算假阳性率。(计算NES也有另外一种方法,是计算出的ES除以排列检验得到的所有ES的平均值)

  4.  Leading-edge subset,对富集得分贡献最大的基因成员。

GSEA分析

软件和基因集下载:http://software.broadinstitute.org/gsea/downloads.jsp

181aae286eab87cc91e97942ec79f8ae.webp

0d9ab4ec02083e7ad642e77cbcb408aa.webp


输入数据准备


  1.  表达矩阵。常见表达矩阵格式,tab键分割,txt格式,第一列为基因名字(名字与注释数据库一致,同为GeneSymbol或EntrezID或其它自定义名字),第一行为标题行,含样品信息。也可为gct文件,具体见 http://blog.genesino.com/2014/08/GSEA-usages/


    cf4f95218e6e75d152441834753deb82.webp

  2. 样品分组信息

    decdc6f2983fc23399568f3928549886.webp

    分组信息示例


    b9e8242f9c6a57f07dace2fe5fb85588.webp

  3. 基因集信息

    6c856495f209b0b66969bf8cf0459238.webp

    基因集信息示例


    7e191564a6e78431e4f8301f17fc189d.webp

软件运行 (每一步的步骤如有不明确的参考文后第一个链接)


  1. 导入数据


    3bc3eb7e18af745ac3b9cd41a79c01d1.webp


  2. 运行GSEA (若每组样品都有多于7个样品,则Permutation type选择phenotype,结果理论上更好;否则选择gene_set)


    7c62bb52524a69ca777939d56eaadbdd.webp


  3. 设置好参数后,点击正下方的run,等待运行结束,左侧出现success


    04f5932d2c878ecee4e18113bebda871.webp


  4. 点击success,查看结果


    8fa838e8a635b8615b7701c19c302c61.webp

顺着网页的导航一步步去查看结果,有耐心就好。主要的解释下,最常见的这种图。

32bc55aef340e8085d2d74e32bc3385c.webp


  1. 图最上面部分展示的是ES的值计算过程,从左至右每到一个基因,计算出一个ES值,连成线。最高峰为富集得分(ES)。在最左侧或最右侧有一个特别明显的峰的基因集通常是感兴趣的基因集。

  2. 图中间部分每一条先代表基因集中的一个基因,及其在基因列表中的排序位置。

  3. 最下面部分展示的是基因与表型关联的矩阵,红色为与第一个表型(MUT)正相关,在MUT中表达高,蓝色与第二个表型(WT)正相关,在WT中表达高。

  4. Leading-edge subset对富集得分贡献最大的基因成员。若富集得分为正值,则是峰左侧的基因;若富集得分为负值,则是峰右侧的基因。

  5. FDRGSEA默认提供所有的分析结果,并且设定FDR<0.25为可信的富集,最可能获得有功能研究价值的结果。但如果样品数目少,而且选择了gene_set作为Permumation type则需要使用更为严格的标准,比如FDR<0.05

Leading-edge分析

主要对筛选感兴趣的基因有意义;选择一个或多个显著富集的基因集,查看其内Leading-edge基因的表达和重叠状态。


d67a9c53c96b1815bb8c970b12b869e9.webp

1f02dd3c449b13aa5396c7448d6c39ad.webp


MSigDB

GSEA团队整理好的基因集,可用于注释,也可下载下来搜寻自己感兴趣的方向的基因作为一个补充。每个注释都提供了基于Gene SymbolEntrez ID的索引表格。

1ca3278747dc843ea0d62d7b188b7a73.webp

参考

  1. 较早记录的一篇GSEA的使用,有脚本可以转换表达矩阵为gctcls文件作为GSEA的输入。文档为英文,但软件操作步骤还算详细,可配合着看:

    http://blog.genesino.com/2014/08/GSEA-usages/

  2. 最开始学习的教程,每一步操作都比较详细:

    http://www.baderlab.org/Software/EnrichmentMap/Tutorial

  3. GSEA软件和数据集下载:

    http://software.broadinstitute.org/gsea/downloads.jsp

  4. 原文对GSEA原理的讲解是很清晰的,可以读下,关键的内容也都摘录在第一个链接里:

    https://www.ncbi.nlm.nih.gov/pubmed/16199517


推荐阅读

往期精品(点击图片直达文字对应教程)

620cfc27a62010537af24b0754f9b7e3.webp

f5361c0e4aa1d87e76b40b623ee40133.webp

0b6cebdd21a6926117c9ffd0c1d2c862.webp

6b1aa65bbf6d96813eccad98ac5c4e14.webp

239a552ad2ef960d5c4efd3daec675a9.webp

cb20bdda6d078010f9d8c12936ddf097.webp

3e2dac826a7ec8550e00e68a77a0fb47.webp

4f2bda3c33ca2cc7c14ff71be46bef78.webp

c9eac0f23dccd79ec3f0e561bea54f07.webp

8ec09f8a38570212e28b4974a403ec8b.webp

593d26de6efacded43d83e724d04f1bf.webp

41ca5a099defe2c8a7336fa57eeb8477.webp

149548054d6035bfd932bfed011a6b44.webp

0a509223598476f82e5870f2f5383ae7.webp

f62a1352acc6592239957514c2a3711d.webp

7c8ded0a698c6020ee915ac423efbd9c.webp

eedd81d7398376ffcc7e2dd6a7a20828.webp

418df2af3b6eefcb9f5b6e42bd150217.webp

614459a26f9b822eec9d99f6f2b027ca.webp

09fc1e565ce64e2839029fe59bc7d1b1.webp

970157b3ea11ceff0f9067b9ee769413.webp

a48acbdea450dfde65358e10cac5c26d.webp

a17474de5d2e01bb89dcfedc4fc5acb2.webp

656cc9caec87a2bc65ed8bfd83cef077.webp

57c9a88d52e91103d70a60f04b9334c3.webp

73ed32260d7b801991a409ee8d8ff9b8.webp

ac5a8e28758087c4f7072c72772652be.webp

708ccee5f4b88708971565cda8b5fab2.webp


后台回复“生信宝典福利第一波”或点击阅读原文获取教程合集

642c8111ca588b2e9aae67d0e0098cbb.webp





浏览 97
点赞
评论
收藏
分享

手机扫一扫分享

分享
举报
评论
图片
表情
推荐
点赞
评论
收藏
分享

手机扫一扫分享

分享
举报