在生命科學(xué)與生物信息學(xué)領(lǐng)域，KEGG（京都基因與基因組百科全書）是一個至關(guān)重要且廣泛使用的數(shù)據(jù)庫資源。它系統(tǒng)地整合了基因組、化學(xué)和系統(tǒng)功能信息，尤其以其直觀、標準化的通路圖而聞名。本文旨在快速解讀KEGG的核心價值，并概述通路圖相關(guān)數(shù)據(jù)處理的基本思路。
一、KEGG數(shù)據(jù)庫是什么？
KEGG不僅僅是一個數(shù)據(jù)庫，而是一個集成化的知識庫系統(tǒng)。它主要包含以下幾個核心模塊：
- KEGG PATHWAY：核心所在。以圖形方式展示細胞過程、代謝、信號轉(zhuǎn)導(dǎo)等分子相互作用網(wǎng)絡(luò)。每個通路圖都是一個標準化的參考圖譜。
- KEGG GENES： 收集了來自眾多已完成測序生物體的基因/蛋白質(zhì)信息。
- KEGG ORTHOLOGY (KO)： 直系同源基因分類系統(tǒng)，是連接基因與通路的功能單元。將基因注釋到KO編號是進行通路富集分析的關(guān)鍵步驟。
- KEGG COMPOUND/GLYCAN/REACTION： 包含化學(xué)物質(zhì)、糖類、生化反應(yīng)等信息，是代謝通路圖的化學(xué)基礎(chǔ)。
- KEGG BRITE： 提供基因、化合物等實體的功能層級分類體系。
KEGG的核心哲學(xué)是將生物系統(tǒng)視為分子網(wǎng)絡(luò)，通路圖則是這些網(wǎng)絡(luò)的“地圖”。
二、如何“讀懂”KEGG通路圖？
一張KEGG通路圖是一個高度標準化的信息可視化工具：
- 圖形元素：矩形通常代表基因、蛋白質(zhì)或酶；圓形代表化學(xué)化合物；線條箭頭代表相互作用（如激活、抑制、轉(zhuǎn)化）。
- 顏色編碼：這是將自身實驗數(shù)據(jù)（如差異表達基因）映射到通路圖的關(guān)鍵。在官網(wǎng)工具中，用戶可以上傳基因列表，系統(tǒng)會根據(jù)預(yù)設(shè)規(guī)則（如上調(diào)/下調(diào)）將圖中的對應(yīng)基因框著色（如紅/綠），從而直觀看到哪些通路區(qū)域在特定條件下發(fā)生了顯著變化。
- KO標識：圖中的每個基因框都關(guān)聯(lián)著KO編號（如K04503），這是功能統(tǒng)一標識符，不依賴于具體物種。
三、通路相關(guān)數(shù)據(jù)處理的核心流程
處理KEGG通路數(shù)據(jù)通常是為了回答“我的基因列表在哪些通路上有顯著富集？”或“我的代謝物變化影響了哪些通路？”這類問題。基本流程如下：
- 基因標識符轉(zhuǎn)換與KO注釋：
- 輸入：通常是一組有意義的基因列表（如差異表達基因）。
- 關(guān)鍵步驟：將基因標識符（如Gene Symbol, Ensembl ID）映射到KEGG的KO編號。這可以通過KEGG官方API、DAVID、KOBAS或clusterProfiler等R/Bioconductor包完成。KO編號是實現(xiàn)跨物種比較和功能解讀的橋梁。
- 通路富集分析：
- 原理：基于超幾何分布檢驗或類似方法，判斷目標基因集在某個KEGG通路中是否顯著過出現(xiàn)（即“富集”）。
- 工具：常用工具包括R語言的
clusterProfiler、enrichKEGG函數(shù)，在線平臺如DAVID、KOBAS-i等。
- 輸出：得到一個包含通路名稱、富集基因數(shù)、P值、錯誤發(fā)現(xiàn)率（FDR）的表格，用于識別最相關(guān)的生物學(xué)通路。
- 結(jié)果可視化：
- 富集氣泡圖/條形圖：展示排名靠前的富集通路及其顯著性，是展示分析結(jié)果的經(jīng)典方式。
- 通路映射圖：使用KEGG官網(wǎng)的“Color”工具或R包
pathview，將基因表達量、代謝物濃度等定量數(shù)據(jù)映射到通路圖上，生成著色的自定義通路圖，直觀展示變化發(fā)生在通路的具體位置。
- 通路網(wǎng)絡(luò)圖：展示不同通路之間的基因重疊關(guān)系，揭示更上層的功能模塊。
四、實用工具與注意事項
- 常用工具：對于生信分析，R包
clusterProfiler和pathview是進行KEGG富集分析和可視化的一站式利器。在線用戶可直接訪問KEGG官網(wǎng)使用其搜索、瀏覽和著色工具。 - 注意物種差異：KEGG通路是參考圖譜，分析時需確保使用正確的物種數(shù)據(jù)庫進行映射（如hsa代表人類，mmu代表小鼠）。
- 理解局限性：通路富集分析結(jié)果是統(tǒng)計關(guān)聯(lián)，需結(jié)合生物學(xué)知識進行解讀。KEGG通路圖是靜態(tài)的參考模型，不一定反映所有細胞類型或狀態(tài)下的具體情況。
###
掌握KEGG數(shù)據(jù)庫與通路圖，本質(zhì)上是獲得了一套解讀基因組數(shù)據(jù)功能意義的“詞典”和“地圖”。從基因列表到KO注釋，再到通路富集分析和可視化映射，這一標準化流程已成為功能基因組學(xué)、代謝組學(xué)等研究的常規(guī)分析手段，能有效將海量的分子列表轉(zhuǎn)化為可理解的生物學(xué)故事。
