2019-12-18 Day9學(xué)習(xí)- 3分0代碼生信套路

看的第二遍了 第一遍看的忘記了 所以整理一下 以防又忘
雖然套路很老 圖也沒R做的好 但借鑒的地方還是很多


image.png

IF 3分左右
PMID: 28901457 PMCID: PMC5780015 DOI: 10.3892/or.2017.5946
在肝癌中的生信分析-確定關(guān)鍵基因
GEO上下載基因,然后對基因進行篩選,篩選差異基因丧靡,然后做PPI網(wǎng)絡(luò)眨攘,蛋白互作網(wǎng)絡(luò)可視化,上下調(diào)基因標(biāo)注,差異基因功能富集,關(guān)鍵基因富集,生存分析看是否有關(guān)系继低。

GEO數(shù)據(jù)下載和篩選差異基因

GEO網(wǎng)站 GSEXXX


image.png

image.png

注意

比較時注意 cancer在下 Nomal在上



按住Ctrl選擇 Shifit連續(xù)選 比較 SAVE ALL
保存TXT 復(fù)制到EXCEL
然后篩選差異基因 條件logFC(fold change)的絕對值≥1 , adj.p-value<0.01
image.png

對多個芯片取交集差異基因(Venn圖)

對篩選好的基因做Venn圖


image.png

保存EXCEL和TXT數(shù)據(jù)

做PPI及上下調(diào)基因標(biāo)注Cytoscape

蛋白互作網(wǎng)路 百度STRING


PPI

image.png

PPI互作網(wǎng)絡(luò),將所篩選的全部差異基因輸入稍走,然后download TSV格式文件袁翁,EXCEL或Cytoscape打開,Cytoscape婿脸。按File+Import+Network+File.
Style美化

上下調(diào)基因篩選

在之前EXCEL文件中對基因列和logFC列進行上下調(diào)標(biāo)注粱胜,VLOOK函數(shù)應(yīng)用,上調(diào)UP盖淡,下調(diào)DOWN,logFC≥1 up年柠,logFC≤-1 down。然后復(fù)制到TXT文檔中,導(dǎo)入Cytoscape冗恨,F(xiàn)ile+Import+table+File答憔。File collar ,按group分組掀抹,選擇顏色形狀虐拓。

顯著Module和關(guān)鍵基因篩選

上下調(diào)基因篩選完成之后,APP+Manger+install MCODE 然后analysis傲武,得出關(guān)鍵基因
記下篩選參數(shù)


image.png

image.png

顯著Module


image.png

image.png

每個seed download

差異基因和最顯著Module的GO及KEGG富集分析

image.png


打開DAVID網(wǎng)站 Functional Annotation,導(dǎo)入差異基因蓉驹,選擇OFFICIAL-GENE-SYMBOL,GENE LIST-SUBMIT LIST.


image.png

image.png

FDR值 矯正后p值

關(guān)鍵基因相互作用網(wǎng)絡(luò)及功能分析

將之前SED保存的關(guān)鍵基因到EXCEL中星澳,復(fù)制到cbioportal中做NETWORK感混。


image.png

image.png

image.png

黑色為關(guān)鍵基因,其余為共表達(dá)基因区转。

功能分析
image.png

Cytoscape疟位,APPS+BINGO然后導(dǎo)入關(guān)鍵基因瞻润,選擇GO種類,物種甜刻,start BINGO


image.png
熱圖

使用在線網(wǎng)站UCSC绍撞,選擇癌癥類型,Genomic得院,Gene Expression


image.png

image.png

image.png

image.png

然后選擇亞型 Grade Symple type傻铣,最后拼接一起


image.png

hub gene生存分析


也是使用在線網(wǎng)站cbioportal,選擇TCGA數(shù)據(jù)祥绞,然后survival

關(guān)鍵基因 OncomineMeta分析

image.png

Search 關(guān)鍵基因非洲,選擇cancer type,選擇研究對象勾上就谜,compare比較怪蔑,記錄p值和fold change,然后拼接

關(guān)鍵基因 Oncomine分析

選擇關(guān)鍵基因 cancer typ 分級


image.png
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末丧荐,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子喧枷,更是在濱河造成了極大的恐慌虹统,老刑警劉巖,帶你破解...
    沈念sama閱讀 218,525評論 6 507
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件隧甚,死亡現(xiàn)場離奇詭異车荔,居然都是意外死亡,警方通過查閱死者的電腦和手機戚扳,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,203評論 3 395
  • 文/潘曉璐 我一進店門忧便,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人帽借,你說我怎么就攤上這事珠增〕瑁” “怎么了?”我有些...
    開封第一講書人閱讀 164,862評論 0 354
  • 文/不壞的土叔 我叫張陵蒂教,是天一觀的道長巍举。 經(jīng)常有香客問我,道長凝垛,這世上最難降的妖魔是什么懊悯? 我笑而不...
    開封第一講書人閱讀 58,728評論 1 294
  • 正文 為了忘掉前任,我火速辦了婚禮梦皮,結(jié)果婚禮上炭分,老公的妹妹穿的比我還像新娘。我一直安慰自己剑肯,他們只是感情好欠窒,可當(dāng)我...
    茶點故事閱讀 67,743評論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著退子,像睡著了一般岖妄。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上寂祥,一...
    開封第一講書人閱讀 51,590評論 1 305
  • 那天荐虐,我揣著相機與錄音,去河邊找鬼丸凭。 笑死福扬,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的惜犀。 我是一名探鬼主播铛碑,決...
    沈念sama閱讀 40,330評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼虽界!你這毒婦竟也來了汽烦?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,244評論 0 276
  • 序言:老撾萬榮一對情侶失蹤莉御,失蹤者是張志新(化名)和其女友劉穎撇吞,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體礁叔,經(jīng)...
    沈念sama閱讀 45,693評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡牍颈,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,885評論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了琅关。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片煮岁。...
    茶點故事閱讀 40,001評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出画机,到底是詐尸還是另有隱情冶伞,我是刑警寧澤,帶...
    沈念sama閱讀 35,723評論 5 346
  • 正文 年R本政府宣布色罚,位于F島的核電站碰缔,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏戳护。R本人自食惡果不足惜金抡,卻給世界環(huán)境...
    茶點故事閱讀 41,343評論 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望腌且。 院中可真熱鬧梗肝,春花似錦、人聲如沸铺董。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,919評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽精续。三九已至坝锰,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間重付,已是汗流浹背顷级。 一陣腳步聲響...
    開封第一講書人閱讀 33,042評論 1 270
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留确垫,地道東北人弓颈。 一個月前我還...
    沈念sama閱讀 48,191評論 3 370
  • 正文 我出身青樓,卻偏偏與公主長得像删掀,于是被迫代替她去往敵國和親翔冀。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,955評論 2 355

推薦閱讀更多精彩內(nèi)容