首先憋飞,一般我們分析水稻的數(shù)據(jù)時候敏弃,選用的是MSU和RAPDB這兩個數(shù)據(jù)庫的genome和gtf文件笛辟,這里只介紹MSU的ID羽利,RAPDB的同理宫患。
水稻GO分析的網(wǎng)站很多(這里不涉及R包,只說一些典型的在線的)
1 AGRIGO2 http://systemsbiology.cau.edu.cn/agriGOv2/ 支持多種ID这弧,包括MSU
2 RIGW http://rice.hzau.edu.cn/cgi-bin/rice2/enrichment 只支持水稻的轉(zhuǎn)錄本ID娃闲,可做KEGG,個人覺得結(jié)果還可以
3 PlantGSEA http://structuralbiology.cau.edu.cn/PlantGSEA/analysis.php 只支持MSU ID
4 PANTHER http://www.pantherdb.org/ 可視化比較漂亮的当宴,也是本文所涉及的
5 CARMO:http://bioinfo.sibs.ac.cn/carmo/result.php?job_id=1625924324108758969
只更新到 2015年畜吊,支持 LOC ID
1、我們通過差異分析户矢,得到一串像這樣的基因
但是PATHER支持的ID為
所以我們需要將MSU ID轉(zhuǎn)換為 Uniprot ID玲献,這里介紹一個在線網(wǎng)站PlantGSEA
轉(zhuǎn)換完成后,為以下這樣的結(jié)果梯浪,在EXCEL中將Uniprot ID復(fù)制出來即可
-2 進入正片捌年,將Uniprot ID粘貼到PANTHER中
-3 出圖, 五個參數(shù)對應(yīng)的功能如字面意思,包括GO分析挂洛,蛋白功能注釋礼预,Pathway分析
我特別喜歡的一個小功能,就是可以將你感興趣的那個通路虏劲,“露出尖尖角”
-4 當然我們要選擇顯著性分析
顯著性分析結(jié)果