GO钻蔑,KEGG啥刻,GSEA富集分析筆記

what is Gene Ontology ?——基因本體論

1 , 什么是本體論咪笑?
簡(jiǎn)單來(lái)說(shuō)可帽,就是我們對(duì)一個(gè)具體事物進(jìn)行分類(lèi)并予以描述。
例如:貓是一種哺乳動(dòng)物/貓是貓科動(dòng)物/貓是一種生活在陸地的生物等等窗怒。對(duì)于一事物映跟,我們可以用不同的分類(lèi)加以描述。因此扬虚,對(duì)于gene的本體論努隙,是對(duì)gene的一種描述。而對(duì)gene的描述大概分三種:
① Cellular component 簡(jiǎn)稱CC
② Biological process 簡(jiǎn)稱BP
③ Molecular Function 簡(jiǎn)稱MF
下面一一進(jìn)行介紹
① 主要是看gene存在在哪辜昵,比如存在在細(xì)胞質(zhì)還是細(xì)胞核荸镊,假如存在在細(xì)胞質(zhì),又存在在哪個(gè)細(xì)胞器呢堪置?亦或者存在在線粒體中躬存,是在線粒體膜上,還是線粒體機(jī)制中舀锨。
②指它能參與哪些生物學(xué)過(guò)程岭洲,比如說(shuō)能夠參與RNA的加工,DNA復(fù)制等
③指從分子層面的功能是什么坎匿,比如說(shuō)它具有什么催化功能分子功能等
所以盾剩, 我們可以從這三個(gè)方面得到基因的注釋信息雷激。

2, 讓我們來(lái)捋一捋具體的思路:
① RNA-Seq 一般分為處理組和對(duì)照組 彪腔, 即contral 和 treatment , 然后根據(jù)ctrl 和 treatment的表達(dá)譜進(jìn)行比較找出來(lái)DEG(differential expression gene : 差異表達(dá)基因侥锦,后續(xù)會(huì)更新怎么找出差異基因)。
② 我們得到DEG后德挣,對(duì)其進(jìn)行GO annotation (GO 注釋恭垦,即上面介紹的對(duì)gene 的描述), 我們期望可以在這三方面找到我們想要的以及未知的發(fā)現(xiàn)。
③那么下面一步就略微難理解格嗅, 比如說(shuō)番挺,我有200個(gè)DEG ,100個(gè)在細(xì)胞核中,那么我們可以說(shuō)基因富集在了細(xì)胞核嗎屯掖?所以引出另外一個(gè)問(wèn)題玄柏,How to test the GO is enriched ? 所以就有了GO 富集分析。對(duì)于模式生物(比如人)贴铜,GO注釋已經(jīng)有了完備的數(shù)據(jù)庫(kù)可以使用粪摘。

④KEGG富集分析對(duì)于人而言,也有了完備的數(shù)據(jù)庫(kù)可以使用绍坝,我們只要將DEG整理到一個(gè)數(shù)據(jù)集徘意,使用R便可得到富集結(jié)果(后續(xù)會(huì)更新怎么進(jìn)行GO ,KEGG富集分析)。個(gè)人理解的原理是轩褐,請(qǐng)思考一個(gè)問(wèn)題椎咧,人體有20000多個(gè)基因,假如由1000個(gè)基因能被注釋成通路A把介,而有500個(gè)差異基因勤讽,其中300個(gè)注釋到了通路A,那么這500個(gè)基因是否富集到了通路A拗踢?
計(jì)算公式

  個(gè)人通俗理解:從2000個(gè)基因碰到注釋為通路A的概率為1/20 , 而500個(gè)差異基因里就有300個(gè)脚牍,概率為3/5,所以說(shuō)明跟這條通路可能有關(guān)系巢墅,當(dāng)然具體還是要算P值诸狭。
⑤GSEA分析:Gene Set Enrichment Analysis (基因集富集分析)評(píng)估一個(gè)預(yù)先定義的基因集的基因在與表型相關(guān)度排序的基因表中的分布趨勢(shì),從而判斷其對(duì)表型的貢獻(xiàn)砂缩。需要輸入DEG基因集 作谚,還有其表達(dá)矩陣三娩,R會(huì)基因根據(jù)其于表型的關(guān)聯(lián)度(可以理解為表達(dá)值的變化)從大到小排序庵芭,然后判斷基因集內(nèi)每條注釋下的基因是否富集于表型相關(guān)度排序后基因表的上部或下部,從而判斷此基因集內(nèi)基因的協(xié)同變化對(duì)表型變化的影響雀监。與前面的GO和KEGG富集分析不同,GO和KEGG富集分析是設(shè)置閾值篩選差異基因双吆,再判斷差異基因在哪些注釋的通路存在富集眨唬;這涉及到閾值的設(shè)定,存在一定主觀性并且只能用于表達(dá)變化較大的基因好乐。而GSEA則不局限于差異基因匾竿,從基因集的富集角度出發(fā),可以包括我們篩選掉的表達(dá)變化小卻對(duì)通路細(xì)微作用和協(xié)調(diào)功能的基因蔚万。
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末岭妖,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子反璃,更是在濱河造成了極大的恐慌昵慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,635評(píng)論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件淮蜈,死亡現(xiàn)場(chǎng)離奇詭異斋攀,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)梧田,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,543評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門(mén)淳蔼,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人裁眯,你說(shuō)我怎么就攤上這事鹉梨。” “怎么了未状?”我有些...
    開(kāi)封第一講書(shū)人閱讀 168,083評(píng)論 0 360
  • 文/不壞的土叔 我叫張陵俯画,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我司草,道長(zhǎng)艰垂,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 59,640評(píng)論 1 296
  • 正文 為了忘掉前任埋虹,我火速辦了婚禮猜憎,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘搔课。我一直安慰自己胰柑,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,640評(píng)論 6 397
  • 文/花漫 我一把揭開(kāi)白布爬泥。 她就那樣靜靜地躺著柬讨,像睡著了一般。 火紅的嫁衣襯著肌膚如雪袍啡。 梳的紋絲不亂的頭發(fā)上踩官,一...
    開(kāi)封第一講書(shū)人閱讀 52,262評(píng)論 1 308
  • 那天,我揣著相機(jī)與錄音境输,去河邊找鬼蔗牡。 笑死颖系,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的辩越。 我是一名探鬼主播嘁扼,決...
    沈念sama閱讀 40,833評(píng)論 3 421
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼黔攒!你這毒婦竟也來(lái)了趁啸?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 39,736評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤督惰,失蹤者是張志新(化名)和其女友劉穎莲绰,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體姑丑,經(jīng)...
    沈念sama閱讀 46,280評(píng)論 1 319
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡蛤签,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,369評(píng)論 3 340
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了栅哀。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片震肮。...
    茶點(diǎn)故事閱讀 40,503評(píng)論 1 352
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖留拾,靈堂內(nèi)的尸體忽然破棺而出戳晌,到底是詐尸還是另有隱情,我是刑警寧澤痴柔,帶...
    沈念sama閱讀 36,185評(píng)論 5 350
  • 正文 年R本政府宣布沦偎,位于F島的核電站,受9級(jí)特大地震影響咳蔚,放射性物質(zhì)發(fā)生泄漏豪嚎。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,870評(píng)論 3 333
  • 文/蒙蒙 一谈火、第九天 我趴在偏房一處隱蔽的房頂上張望侈询。 院中可真熱鬧,春花似錦糯耍、人聲如沸扔字。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 32,340評(píng)論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)革为。三九已至,卻和暖如春舵鳞,著一層夾襖步出監(jiān)牢的瞬間震檩,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 33,460評(píng)論 1 272
  • 我被黑心中介騙來(lái)泰國(guó)打工系任, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留恳蹲,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,909評(píng)論 3 376
  • 正文 我出身青樓俩滥,卻偏偏與公主長(zhǎng)得像嘉蕾,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子霜旧,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,512評(píng)論 2 359