基于生物信息學預測分析C/EBPα靶基因及其在急性髓系白血病中的潛在機制

章瑞琪，郭 虹

（麗水學院醫學院，浙江 麗水 323000）

急性髓系白血?。╝cute myeloid leukemia，AML）是一種骨髓和外周血髓系細胞發育受阻的惡性腫瘤，主要表現為髓系原始細胞增殖不受控制，骨髓分化受抑制，從而抑制正常血細胞的發育。該病占白血病總發病數的50%以上，其中在AML 亞型之中M2與M5的發病率較高［1］。而AML的發生、發展、治療和預后與轉錄因子CCAAT 增強子結合蛋白α（CCAAT/enhancer binding protein alpha，C/EBPα）的異常表達密切相關。在中國AML成年患者中，約有11.7%患者的CEBPA基因突變陽性［2］，這干擾了CEBPA基因的正常表達，阻礙了其功能的發揮。而正常表達的轉錄因子C/EBPα因含有堿性亮氨酸拉鏈結構域，可以通過形成二聚體、與目標基因的啟動子區結合、調控轉錄激活等手段發揮其對維持血液系統穩定的功能：維持造血干細胞的自我更新、支持粒細胞的分化，并參與調節細胞周期與新陳代謝過程。

為進一步闡釋C/EBPα 的異常表達與AML 間的作用機制，筆者利用生物信息學技術探究C/EBPα的靶基因，并通過分析靶基因的功能與信號轉導通路，進一步了解C/EBPα 的轉錄調控功能，且重點揭示靶基因信號轉導通路（JAK-STAT信號通路、MAPK通路、PI3K-Akt通路），總結它們對維持正常血液系統的功能和在AML 發生、發展及治療過程中的作用與影響，揭開C/EBPα通過靶基因信號轉導通路在AML 中的潛在機制，并進而為AML臨床治療與預后提供依據。

1 材料與方法

1.1 數據獲取

通過GeneCards 數據庫（http：//www.genecards.org）檢索得到107個與AML相關的基因，將其作為候選預測靶基因。在基因表達數據庫（https：//www.ncbi.nlm.nih.gov/geo/）中檢索得到包含敲除Cebpa基因的實驗組與對照組的GSE146288、GSE153622、GSE71687數據集。

1.2 篩選差異表達基因

運用GEO 數據庫提供的GEO2R 分析工具對上述107 個候選預測靶基因的表達量進行綜合分析，篩選Cebpa敲除前后的差異表達基因，得出其在Cebpa基因敲除前后的比較結果。

1.3 差異基因的GO功能與KEGG富集分析

為進一步解釋差異表達基因在AML發病過程中發揮的作用，對C/EBPα的42個差異表達靶基因進行分析。借助DAVID 數據庫綜合分析工具（http：//david.ncifcrf.gov）進行基因本體論（gene ontology，GO）富集分析、京都基因與基因組百科全書（kyoto encyclopedia of genes and genomes，KEGG）富集分析，篩選顯著差異的數據（P＜0.05），探索差異基因潛在影響的功能與通路。

1.4 C/EBPα 與信號通路中差異表達基因的蛋白質互作網絡

利用STRING 數據庫（https：//cn.string-db.org/）分析調控網絡中的蛋白質-蛋白質相互作用（protein-protein interaction，PPI），對KEGG 富集分析結果中的JAK-STAT信號通路、MAPK信號通路與PI3K-Akt 信號通路中所包含的基因進行分析，構建其與C/EBPα之間的蛋白質互作網絡。

2 結果

2.1 篩選差異表達基因結果

運用GEO2R分析工具對上述107個候選預測靶基因的表達量進行綜合分析，得出其在Cebpa基因敲除前后的比較結果（見表1），將表達量隨Cebpa基因敲除而下降的基因作為預測靶基因，共42個，并對這些差異表達基因開展深入研究。

表1 Cebpa基因敲除后AML相關基因表達

2.2 差異基因的GO功能富集

借助DAVID 綜合分析工具對C/EBPα的42個差異表達靶基因進行GO富集分析，篩選顯著差異的數據（P＜0.05），得到GO 富集分析結果：主要富集于26個基因功能，按生物學過程、細胞組分和分子功能3類基因功能分類，分類結果如圖1所示。

圖1 C/EBPα預測靶基因的GO功能富集分析

由觀察分析結果可知：在生物學過程分類中，差異基因主要富集于RNA聚合酶Ⅱ啟動子轉錄的正調控、造血、基因表達的正向調控等。在細胞組分分類中，差異基因主要富集于細胞核、轉錄因子復合物、染色質等。在分子功能分類中，差異基因主要富集于蛋白質結合、染色質結合、轉錄因子結合、DNA結合等。

2.3 差異基因的KEGG通路分析

在GO功能富集分析的基礎上，利用DAVID綜合分析工具進行KEGG富集分析，篩選差異顯著的數據（P＜0.05），得到KEGG 富集分析結果，結果如圖2所示：差異基因主要富集于癌癥的通路、急性髓性白血病、癌癥中的轉錄失調、癌癥中的微小RNA、慢性粒細胞白血病、PI3K-Akt 信號通路、JAKSTAT信號通路、乙型肝炎、造血細胞譜系、MAPK信號通路。差異基因在癌癥的通路中富集最多，共有15 個基因富集于此通路。借助STRING 在線數據庫，繪制PI3K-Akt 信號通路、JAK-STAT信號通路與MAPK信號通路網絡圖，結果如圖3所示。

圖2 C/EBPα預測靶基因的KEGG信號轉導富集分析

圖3 KEGG信號通路的蛋白質互作網絡

2.4 C/EBPα與信號通路的蛋白質互作網絡

分別將JAK-STAT信號通路、MAPK信號通路與PI3K-Akt 信號通路所包含的蛋白基因導入STRING在線數據庫，獲取上述通路中所包含差異基因與C/EBPα的相互作用網絡，結果如圖4所示：C/EBPα與JAK-STAT信號通路所包含的差異基因間共有6 條相互作用關系；C/EBPα 與MAPK 信號通路所包含的差異基因間共有6條相互作用關系；C/EBPα 與PI3K-Akt 信號通路所包含的差異基因間共有5條相互作用關系。

圖4 C/EBPα與MAPK、JAK-STAT及PI3K-Akt信號通路的蛋白質互作網絡

3 討論

轉錄因子C/EBPα 的異常表達對AML 的發病及預后產生重要影響，研究表明C/EBPα能夠調節造血干細胞的細胞周期，促進中性粒細胞的分化，對骨髓發育起著重要作用。大部分AML患者的C/EBPα 表達量下降，但C/EBPα 不同的突變形式影響著AML 患者的預后，C/EBPα 單突變的AML 患者往往預后較差［3］，而C/EBPα雙突變的AML患者的死亡率較低，往往預后較好［4］。這說明C/EBPα的不同表達會對AML的發病與預后產生不同的影響，借助生物信息學技術，研究C/EBPα 靶基因功能及其信號通路，有利于進一步揭示AML 的發病機理，同時為臨床開發針對C/EBPα靶點的免疫治療藥物提供理論支持。

本研究主要運用生物信息學技術，在GEO 數據庫中篩選敲除Cebpa基因的實驗組及正常造血干細胞間的差異表達基因，發現差異表達基因共86 個，其中共有42 個基因隨Cebpa基因的敲除而表達下降，將其作為預測靶基因進行了GO 功能分析和KEGG 通路富集分析，結果發現這些基因在RNA聚合酶Ⅱ啟動子轉錄的正調控、造血、基因表達的正向調控等功能上富集，并且很可能通過JAK-STAT信號通路、MAPK信號通路及PI3K-Akt信號通路發揮作用。不僅如此，我們還證明了C/EBPα 與JAK-STAT 信號通路、MAPK 信號通路、PI3K-Akt信號通路中所包含的差異表達基因有互相作用，提示C/EBPα 功能的發揮與上述3 條通路的信號傳導具有密切聯系，推測轉錄因子C/EBPα可能通過JAK-STAT 信號通路、MAPK 信號通路、PI3K-Akt信號通路的傳導對急性髓系白血病的發病與預后產生一定影響。

已有研究已論證：JAK-STAT信號通路在非典型C/EBPα雙突變的AML患者中會被異常激活，破壞其原有抗腫瘤作用，使得穩態失衡，進而加劇AML患者的病情發展［5］；MAPK信號通路主要包括MEK1/2-ERK1/2 通路、JNKs 通路、p38 通路及MEK5-ERK5-MEF2C 通路4 類信號通路，在調控細胞增殖、凋亡及細胞周期方面發揮著重要作用［6-8］，其中MEK1/2-ERK1/2 通路活性的增強，會下調C/EBPα 的表達量，從而抑制粒細胞分化，進而誘導AML 的發生［9-10］，而MEK5-ERK5-MEF2C通路也能夠通過調控轉錄因子C/EBPα 與C/EBPβ來引導髓系分化［11］；PI3K-Akt 信號通路能夠控制C/EBPα 的磷酸化，以此介導骨髓細胞正常發育［12］，而該通路的過度激活往往會導致AML 細胞不受控制的增殖［13］。這些論據均說明C/EBPα 的異常表達對AML患者的影響與JAK-STAT信號通路、MAPK 信號通路、PI3K-Akt 信號通路的傳導密不可分。

綜上所述，本研究篩選出與C/EBPα表達相關的差異基因，并分析所述基因的功能及其相關信號通路，研究結果有助于探索尋求臨床治療AML的免疫靶點，為揭示AML 發生發展過程提供理論基礎。