前段時間,推過兩個教程,分別是:
1.《零基礎(chǔ)快速完成基因功能注釋 / GO / KEGG / PFAM...》
2.《「GO富集分析」從原理到實踐 ~ 零基礎(chǔ)掌握》
可以說,如果這系列教程還不完整,那么缺少的就是一個 KEGG 通路富集分析教程。而這個功能,也是 TBtools 的基本功能之一。很多時候,我更多撰寫或講演 GO 富集分析相關(guān)內(nèi)容,而忽略了 KEGG 通路富集相關(guān)。主要原因是,前者懂了,后者往往也就懂了,畢竟富集分析原理一樣,只是注釋信息有不同。
當然,在具體軟件使用上也有點點區(qū)別,甚至有不少時候,KEGG Pathway 富集分析能提供我們更為具體的分析結(jié)果。
使用這一功能,首先是打開對應功能
可以看到界面
注意到不同選項:
設(shè)置一個KEGG數(shù)據(jù)庫文件,一般在 TBtools 中,叫做 .KEGG.BackEnd 文件,類似 Gene Ontology 的 go-basic.obo 文件;這一文件可以用《TBtools Cookbook》中找到下載位置
https://www.yuque.com/cjchen/hirv8i/ghaxyl
也可以通過在選擇 6. 中選項,并點擊 7. Mack One Backend File From Web 來準備。不過這依賴于當前網(wǎng)絡與KEGG官網(wǎng)網(wǎng)絡的通暢程度。一般建議直接從Cookbook 給出的鏈接下載。
設(shè)置KEGG 注釋信息文件,基因ID和KEGG Knum各一列,使用制表符分隔。如果不知道如何準備,建議參考前述教程《零基礎(chǔ)快速完成基因功能注釋 / GO / KEGG / PFAM...》。
設(shè)置基因選擇集 ID 列表,一行一個ID,可以是差異表達基因集合,或者 GWAS 定位到的區(qū)間基因ID,或者是各類合理方式獲得的基因ID集合
設(shè)置一個輸出目錄,或可以給出輸出文件名字前綴
點擊 Start 即可開始
輸出兩個結(jié)果
一般只看.final.xls文件
具體含義,請參考《「GO富集分析」從原理到實踐 ~ 零基礎(chǔ)掌握》并自行理解。
當然,如果你想可視化
與 GO 富集分析可視化同個功能
輸出結(jié)果如下
具體用戶感興趣一些可視化參數(shù),請自行調(diào)整。
Emmm,果然,我還是不太想寫。不過我仍然相信,通過這個教程,大伙可以點點點,完成 KEGG Pathway通路富集分析。