富集分析

新知识

关于对于生信学习的反思,最近一直接触大量的新知识,说实话,很难,因为这些新知识没有办法和我现有的知识形成联系,所有的概念都是新的,如果我要记笔记,参考的文章为十,那我几乎要把十全部整理下来,这让我的心态很不好。

原理及目的

基因富集分析(GeneSet Enrichment Analysis,GSEA)一种对基因进行富集分析的方法,检验已知功能的基因集(即gene set,可以是从GO/KEGG/hallmark/MsigDB中拎出的某一特定类别的基因集合,也可以是自定义的功能基因集合),在一个依据与表型的相关度进行排序的基因列表(即两组样品的表达谱数据,依据基因在两种表型中的表达量的高低进行排序,因未对基因进行显著差异的筛选,因而可以将全部基因与不同表型的相关性均考虑进去)中是随机排列还是主要集中在列表的顶部或底部。若研究的已知功能基因集是非随机分布的,则说明该已知功能基因集与表型相关,根据其基因的集中情况,则可以推断出该已知的功能具体和哪种表型更为接近。

分类

GO富集分析

Gene Ontology: 描述基因的层级关系【基于ORA算法】可以算得上是高通量数据分析的标配,转录组、甲基化、ChIP-seq、重测序等,都会用到对一个或多个集合的基因进行功能富集分析,来找这个基因集的功能偏好性
这是一个基因本体联合会组织(Gene Ontology Consortium)建立的数据库,规范统一了对于不同物种的基因和蛋白描述。

条目标准定义

id:也就是GO编号,如:GO:0031985
name:全称,Golgi cisterna
ontology:命名分类cellular_component
definition:定义,Any of the thin, flattened membrane-bounded compartments that form the central portion of the Golgi complex.Source:GOC:mah

条目之间的关系

有向无环图

有向无环图

有向指的是term之间的单向指向性关系,比如termA是内质网,termB是细胞器,规定A是B,却不能说B是A;无环指的是从任何一点开始沿着规定的指向都不能回到原点
具体作用
基因分门别类放入细胞组分CC、分子功能MF和生物过程BP,分别对应基因产物在那里发挥功能怎样发挥功能发挥什么样的功能
GO分析可以帮我们找到提交的基因集中各个基因是否有共同的GO条目,或者有没有共同的上级GO条目,可以发现某些具有共同特点的基因(比如在某条共同的通路中起作用)
表现形式
气泡图表示显著性不同的GO条目对应的基因数量

富集气泡图

X轴:RichFactor,富集因子,是指前景基因集中属于这个term的基因的数量/背景基因集中富集在这个term中所有基因的数量;
Y轴:GO term名称;
气泡颜色:Q值(也可以用P值绘图),代表富集显著程度,在这个图形当中,颜色越红代表Q值越小,富集程度越高;
气泡大小:数量,前景基因集中属于这个term的基因数量。
气泡越红越大,富集的基因数目越多,富集越显著

条形图表示三个分类的前多少条目对应基因的数量
网络图可以表示显著性较强的条目之间的相互关系
map图表示显著性较强的GO条目的层级关系,以树状图的形式展现

KEGG

Kyoto Encyclopedia of Genes and Genomes: 系统分析基因产物和化合物在细胞中的代谢途径以及这些基因产物的功能的数据库【基于ORA算法】

包括

代谢通路(KEGG PATHWAY)、药物(KEGG DRUG)、疾病(KEGG DISEASE)、功能模型(KEGG MODULE)、基因序列(KEGG GENES)及基因组(KEGG GENOME)等等
有一套完整的KO注释系统(KO是蛋白质或酶的一个分类体系,将同一条通路上功能相似、序列相似的蛋白质归为一类),可完成新测序物种的基因组或转录组的功能注释
因此可以将基因一个个归置到代谢网络指定位置

命名规则

K(大写) +num基因ID号,表示所有同源物种中具有相似结构或功能的一类同源蛋白,如:K04456表示丝氨酸蛋白激酶;
ko+num代谢通路,表示特定的生物路径,如:ko04151表示PI3K-Akt信号通路【也是我们常用的代谢通路】;
M+num表示模块,如:M00676表示PI3K-Akt信号模块
C+num表示化合物,如:C00533表示一氧化碳
ECx.x.x.x表示酶,如:EC2.7.11.1表示丝氨酸
R+num表示反应名称
举个例子:三个字母表示物种,hsa表示Homo sapiens;具体的KO号,如K12407表示和葡糖激酶glucokinase序列和功能相近的蛋白质/酶类,当然一个KO号有可能会对应好几个数字(基因登录号),表示细胞中存在几种不同的葡糖激酶,分别由以上几种数字表示的基因编码


KEGGpathview

主要看颜色:红色表示上调的差异基因,绿色表示下调,白色的表示没有差异基因成功注释


KEGG条形图

挑选显著分析的前20左右的pathway进行展示
X轴:Gene Percent(%),柱子长短代表前景基因富集在该pathway上数目占所有前景基因的百分比。柱子上的数字为基因数量,和对应的q值;
Y轴:Pathway名称;
柱状图颜色:Q值(也可以用P值绘图),代表富集显著程度,颜色越深代表Q值越小,富集程度越高。

GSEA

Gene Set Enrichment Analysis 基因集富集分析,用于评估一个基因集的基因在表型相关度排序中的分布趋势,进而判断它们对表型的贡献
FCS算法
区别

  1. GO是先筛选基因(需要人为设定阈值),再判断差异基因在哪些通路有富集
  2. GSEA可以考虑那些表达差异不大却功能重要的基因对通路影响,相比GO和KEGG能保留更多信息

首先将我们要研究的基因集里的基因按照差异倍数或者差异显著性P值等进行排序,再和已知功能的基因集进行对比,看后者在前者的分布情况,然后根据富集得分(ES值)去判断我们要研究的基因集是不是和我们关注的通路相关。

GSEA富集图像

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

推荐阅读更多精彩内容