理解差异表达与GO分析

记录跟差异基因分析相关的几个概念,主要摘自《R与Bioconductor》一书。

基因表达差异的显著性分析

简称表达差异分析,其目的是比较两个条件(包括种属、表型等)下的基因表达差异,通过一定的统计学方法,从中识别出与条件相关的特异性基因,然后进一步分析这些特异性基因的生物学意义。

基因表达差异分析的第一步是要识别在两个条件下有显著性表达差异的基因,简称差异表达基因。那么怎样才能称得上显著性表达差异?通常是指一个基因在两个条件中表达水平的检测值差异具有统计学意义,往往基于一定的统计假设。

常用的算法有三类:

  • 倍数分析 - 计算每一个基因在两个条件下的比值,若大于给定阈值,则为差异表达基因;
  • 经典统计模型(如t检验)方法 - 计算表达差异的置信度,选取一定P值以下的作为差异表达基因;
  • 机器学习 - 进行特征(基因)选择,包括贝叶斯模型、支持向量机或者随机森林等。

得到差异表达分析后,通常会进行基因本体论和通路分析(通常称为下游分析)。

基因本体论分析

基因本体论(Gene Onotology, GO)分析包括GO注释和富集分析。

GO是基因本体联合会(Gene Ontology Consortium)所建立的数据库,它由一组预先定义好的术语(GO term)组成,这组术语对基因和蛋白质功能进行限定和描述,适合于各种物种,并能伴随研究不断深入而更新。GO中的每个术语都有唯一的一个GO ID,GO的所有术语由有向无环图来相互联系,术语之间通常有三种关系:"is_a","part_of"和"regulates"。

Gene Ontology可分为分子功能(Molecular Function),生物过程(biological process)和细胞组成(cellular component)三个部分。蛋白质或者基因可以通过ID对应或者序列注释的方法找到与之对应的GO号,而GO号可对应到Term,即功能类别或者细胞定位。

GO注释,就是将表示基因或其产物的ID映射到一组GO的ID上,用这组GO term来描述这个基因。实际应用中,人们更关心差异表达基因(或一组基因)的共同点,分析它们所对应GO的分布情况。

通过加入一些统计模型,使分析更加深入,可能会产生更有生物意义的发现。GO富集分析,就是基于这么一个思想的分析方法。GO富集分析的统计学基础是超几何分布,简单而言就是根据Fisher精确检验对每个GO term计算一个P值。
$$
P=\frac{(\frac{M}{k})(\frac{N-M}{n-k})}{(\frac{N}{n})}
$$
公式显示查看博客

  • N表示此次研究中(比如一个芯片上)所有基因总数;
  • n表示N中差异表达基因的总数;
  • M表示N中术语某个GO term的基因个数
  • k表示n中属于某个GO term的基因个数

通路分析

通路分析包括通路注释和通路富集分析。通路富集分析的基本思路、统计模型等和GO富集分析如出一辙。常用的公共通路数据库主要有KEGG(Kyoto encyclopedia of genes and genomes)、BioCarta和GenMAPP(Gene map annotator and pathway profiler),最为著名的是KEGG库中的代谢通路,它几乎成了通路注释的代名词。

附上一个比较详细的介绍链接:基因芯片数据GO和KEGG功能分析

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

推荐阅读更多精彩内容

  • 以下介绍引用自 表达谱芯片数据的基因功能富集分析 刘 明1 王米渠2 丁维俊2 综述 毕 锋1 审校 基因富...
    赵晨西西西西阅读 3,467评论 0 8
  • 谈谈本体论 为什么需要本体论 作为一位大学统计棉花表皮毛的苦逼生物狗,深刻体会什么叫做经验,也就是人类模式识别能力...
    xuzhougeng阅读 2,976评论 0 10
  • 回顾 首先对基因表达分析(上)做一个简单的回顾 研究基因表达的有如下工具:RNA-Seq,microarray, ...
    xuzhougeng阅读 50,475评论 8 115
  • 定义为创建一组相关或相互依赖的对象提供一个接口,而且无需指定它们的具体类 特点抽象工厂是工厂方法模式的升级版本,用...
    ting723阅读 331评论 0 1
  • 公司项目紧,好久没有写文章。最近几天打开简书发现写的两篇关于maven的文章尽然还有读者在阅读,那必定maven在...
    zlwen阅读 3,657评论 4 23