KEGG富集分析

　　分析模块，输入差异基因列表文件、全部基因列表文件和KEGG注释信息文件。差异基因作为foreground，所有基因作为background，进行差异基因KEGG富集分析，输出富集结果。

分析模块，使用软件KOBAS（http://kobas.cbi.pku.edu.cn）进行富集分析，默认使用方法为超几何检验/Fisher精确检验。通常情况下，当经过校正后的p值≤0.05时，认为该ko通路显著富集。

！！对于主要物种，软件团队从KEGG网站上，下载并整理了对应物种的KEGG注释信息。访问，VG软件官方网站：（http://www.vgenomics.cn/），进行下载。

　　输入：

       1、差异基因列表文件，可由分析模块“Batch Mode: Retrieve Diff Genes list”获得。

　　示例：

BM590_B0190

BM590_A0615

BM590_B0191

BM590_B0407

BM590_B0166

BM590_A0618

……

       2、全部基因列表文件，可由分析模块“Fetch all genes(row names) list from matrix”通过原始的FPKM矩阵或Count矩阵获得。

　　示例：

BM590_A0001

BM590_A0002

BM590_A0003

BM590_A0004

BM590_A0005

BM590_A0006

BM590_A0007

BM590_A0008

BM590_A0009

……

       3、对应的KEGG注释信息文件，其中，第一列为基因名，第二列为对应的K号。Kddddd表示，在所有同源物种中具有相似结构和功能的一类同源蛋白。如K04456表示丝氨酸/苏氨酸蛋白激酶。

　　示例：

BM590_A0001         K02313

BM590_A0003         K03629

BM590_A0005         K12972

BM590_A0006         K13896

BM590_A0007         K13895

……

　　输出：

　　差异基因KEGG富集分析结果文件。

　　示例：

#Term       Database ID     Input number   Background number         P-Value     Corrected P-Value     Input         Hyperlink

Thiamine metabolism       KEGG PATHWAY     ko00730   4       7       0.00601059318131 0.140625236043         BM590_A0226|BM590_A0220|BM590_A0222|BM590_A0223    http://www.genome.jp/kegg-bin/show_pathway?ko00730/K00788%09red/K03149%09red/K00941%09red/K03707%09red

Nitrogen metabolism         KEGG PATHWAY     ko00910   5       13     0.00760136411042 0.140625236043         BM590_B0293|BM590_B0239|BM590_B0240|BM590_B0290|BM590_B0268   http://www.genome.jp/kegg-bin/show_pathway?ko00910/K00374%09red/K00376%09red/K04561%09red/K02305%09red/K00370%09red

…….

注：其中，每一列的含义如下所示：

l   #Term，代谢通路描述。

l   Database ，数据库类型。

l   ID，ko号，代谢通路名称，表示一个特定的生物路径。

l   Input number，差异基因中，注释到该ko通路的差异基因个数。

l   Background number，全部基因中，注释到该ko的基因个数。

l   P-Value，统计检验的p值。

l    Corrected P-Value，校正后的p值。

l   Input，具体富集的差异基因。

l   Hyperlink，超链接，链接到KEGG官网上对应的代谢通路图。

分析模块引用了KOBAS（v2.0-20150126）软件（http://kobas.cbi.pku.edu.cn）。

　　相关文献如下所示：

       Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., Kong, L., Gao, G., Li, C. and Wei, L. (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res, 39, W316-322.

       KEGG 库(Kyoto encyclopedia of genes and genomes 数据库)：京都基因和基因组百科全书，是系统分析基因功能、联系基因组信息和功能信息的知识库。利用KEGG数据库，可将基因按照参与的pathway通路或行使的功能分类。

      ！！更新KOBAS程序后台的KEGG数据库，通过KOBAS官网链接（http://kobas.cbi.pku.edu.cn/），下载最新的ko.db.gz文件，解压后覆盖软件根目录database文件夹下的同名文件。