GO,KEGG,GSEA富集分析笔记

what is Gene Ontology ?——基因本体论

1 , 什么是本体论?
简单来说,就是我们对一个具体事物进行分类并予以描述。
例如:猫是一种哺乳动物/猫是猫科动物/猫是一种生活在陆地的生物等等。对于一事物,我们可以用不同的分类加以描述。因此,对于gene的本体论,是对gene的一种描述。而对gene的描述大概分三种:
① Cellular component 简称CC
② Biological process 简称BP
③ Molecular Function 简称MF
下面一一进行介绍
① 主要是看gene存在在哪,比如存在在细胞质还是细胞核,假如存在在细胞质,又存在在哪个细胞器呢?亦或者存在在线粒体中,是在线粒体膜上,还是线粒体机制中。
②指它能参与哪些生物学过程,比如说能够参与RNA的加工,DNA复制等
③指从分子层面的功能是什么,比如说它具有什么催化功能分子功能等
所以, 我们可以从这三个方面得到基因的注释信息。

2, 让我们来捋一捋具体的思路:
① RNA-Seq 一般分为处理组和对照组 , 即contral 和 treatment , 然后根据ctrl 和 treatment的表达谱进行比较找出来DEG(differential expression gene : 差异表达基因,后续会更新怎么找出差异基因)。
② 我们得到DEG后,对其进行GO annotation (GO 注释,即上面介绍的对gene 的描述), 我们期望可以在这三方面找到我们想要的以及未知的发现。
③那么下面一步就略微难理解, 比如说,我有200个DEG ,100个在细胞核中,那么我们可以说基因富集在了细胞核吗?所以引出另外一个问题,How to test the GO is enriched ? 所以就有了GO 富集分析。对于模式生物(比如人),GO注释已经有了完备的数据库可以使用。

④KEGG富集分析对于人而言,也有了完备的数据库可以使用,我们只要将DEG整理到一个数据集,使用R便可得到富集结果(后续会更新怎么进行GO ,KEGG富集分析)。个人理解的原理是,请思考一个问题,人体有20000多个基因,假如由1000个基因能被注释成通路A,而有500个差异基因,其中300个注释到了通路A,那么这500个基因是否富集到了通路A?
计算公式

  个人通俗理解:从2000个基因碰到注释为通路A的概率为1/20 , 而500个差异基因里就有300个,概率为3/5,所以说明跟这条通路可能有关系,当然具体还是要算P值。
⑤GSEA分析:Gene Set Enrichment Analysis (基因集富集分析)评估一个预先定义的基因集的基因在与表型相关度排序的基因表中的分布趋势,从而判断其对表型的贡献。需要输入DEG基因集 ,还有其表达矩阵,R会基因根据其于表型的关联度(可以理解为表达值的变化)从大到小排序,然后判断基因集内每条注释下的基因是否富集于表型相关度排序后基因表的上部或下部,从而判断此基因集内基因的协同变化对表型变化的影响。与前面的GO和KEGG富集分析不同,GO和KEGG富集分析是设置阈值筛选差异基因,再判断差异基因在哪些注释的通路存在富集;这涉及到阈值的设定,存在一定主观性并且只能用于表达变化较大的基因。而GSEA则不局限于差异基因,从基因集的富集角度出发,可以包括我们筛选掉的表达变化小却对通路细微作用和协调功能的基因。
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,271评论 5 466
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,725评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,252评论 0 328
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,634评论 1 270
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,549评论 5 359
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 47,985评论 1 275
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,471评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,128评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,257评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,233评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,235评论 1 328
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,940评论 3 316
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,528评论 3 302
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,623评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,858评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,245评论 2 344
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,790评论 2 339