基于生物信息學研究細胞周期蛋白依賴性激酶抑制因子2A基因與宮頸癌的關系及調控機制▲

劉 珍 羅永金,2 胡曉霞

(1 廣西壯族自治區人民醫院婦科，廣西南寧市 530000； 2 廣西醫科大學臨床醫學院，廣西南寧市 530021)

宮頸癌是一種常見的婦科惡性腫瘤，高危型人乳頭瘤病毒感染已明確是其病因之一。盡管在過去幾十年里，宮頸癌的診斷和治療技術有很大進步，但其發病率和死亡率仍高居女性癌癥發病率的第四位，僅2020年全世界就有新發宮頸癌604 127例，有341 831人因宮頸癌死亡[1]。手術和放化療是治療宮頸癌的主要手段，分期較晚、腫瘤復發和轉移是影響宮頸癌患者預后的主要因素[2]，早發現、早治療對于改善宮頸癌患者的預后至關重要。因此，進一步探索宮頸癌診斷和治療的靶點有重要意義。細胞周期蛋白依賴性激酶抑制因子2A(cyclin-dependent kinase inhibitor 2A，CDKN2A)基因的缺失可通過負反饋導致細胞周期蛋白依賴性激酶(cyclin-dependent kinase,CDK)4、CDK6的活性增加，從而促進細胞周期進程[3]。CDKN2A基因異常表達或啟動子區過度甲基化在腫瘤的發生、發展中可能起著至關重要的作用[4-6]。本研究采用生物信息學技術分析CDKN2A基因在宮頸癌組織中的表達情況及其與宮頸癌患者臨床病理特征、預后、腫瘤免疫細胞浸潤的關系，并通過繪制蛋白-蛋白相互作用(protein-protein interaction，PPI)網絡、京都基因與基因組百科全書(Kyoto Encyclopedia of Genes and Genomes，KEGG)通路富集分析、分析CDKN2A基因啟動子的甲基化水平和上游miRNA，來探討CDKN2A基因在宮頸癌中的可能作用機制，為進一步研究宮頸癌的診斷和治療方法提供新思路。

1 資料與方法

1.1 數據來源 從GEO數據庫(https：//www.ncbi.nlm.nih.gov/geo/)中下載同時含有宮頸癌組織與正常宮頸組織基因表達的芯片GSE7803、GSE64217和GSE63514用于篩選差異性表達基因，然后分別采用GPL96、GPL10558和GPL570平臺采集芯片數據。

1.2 基因篩選及驗證 利用GEO數據庫自帶軟件包GEO2R對芯片(GSE7803、GSE64217和GSE63514)數據進行差異性表達基因的篩選，基因篩選標準為|log2FC|>1且adj.P值<0.01；再利用Venn Diagram在線工具(http://bioinformatics.psb.ugent.be/webtools/Venn/)繪制韋恩圖，對篩選出的差異性表達基因取交集，獲得共同的差異性表達基因。利用GEPIA數據庫(http://gepia.cancer-pku.cn/)驗證CDKN2A基因在宮頸鱗狀細胞癌組織中的表達水平。

1.3 差異性表達基因的通路富集分析 針對共同的差異性表達基因，利用DAVID 在線工具( https：//david.ncifcrf.gov/)分析差異性表達基因的KEGG的通路富集情況。

1.4 構建PPI網絡 將3個芯片的差異性表達基因導入STRING數據庫 (www.string-db.org/)，提取含有差異性表達基因的文件；再將獲得的文件導入Cytoscape軟件(Vision 3.8.2)，通過軟件中的Cytohubba插件預測和查找基因的重要節點和子網絡，利用其中的MCC算法，在差異性表達基因編碼蛋白中篩選與CDKN2A基因編碼蛋白相互作用密切的蛋白，得到PPI網絡。

1.5 CDKN2A基因表達與宮頸癌患者臨床病理特征、預后的關系及宮頸癌組織中 CDKN2A基因啟動子甲基化水平 利用 UALCAN在線工具( http://ualcan.path.uab.edu/cgi-bin/ualcan-res.pl)下載TCGA數據庫中宮頸癌鱗狀細胞癌患者基因的表達數據， 分析CDKN2A基因在不同臨床病理特征宮頸癌患者中的表達情況，其中臨床病理特征包括患者體重、年齡、腫瘤分期、腫瘤分化程度、淋巴結轉移情況；分析不同CDKN2A基因表達水平宮頸癌患者的生存期差異；分析宮頸癌組織中CDKN2A基因啟動子區域甲基化水平。以P<0.05為差異具有統計學意義。

1.6 CDKN2A基因表達水平與免疫細胞浸潤的關系 利用TIMER數據庫(https://cistrome.shinyapps.io/timer/)分析宮頸癌中CDKN2A基因表達水平與各類型免疫細胞浸潤水平的相關性，包括B淋巴細胞、CD8+T淋巴細胞、CD4+T淋巴細胞、巨噬細胞、中性粒細胞和樹突狀細胞。

1.7 上游miRNA的篩選 利用在線軟件ENCORI(https://starbase.sysu.edu.cn)和mirDIP(http://ophid.utoronto.ca/mirDIP)分析與CDKN2A基因有結合靶點的miRNA，利用Venn Diagram在線工具繪制韋恩圖，獲取兩個軟件均能預測到的miRNA交集。

2 結 果

2.1 CDKN2A基因在宮頸癌組織中的表達情況 從數據集GSE7803、GSE64217和GSE63514中分別得到差異性表達基因919個、4 324個和3 845個，經過Venn Diagram在線工具取交集后獲得347個共同差異性表達基因 ，其中CDKN2A基因在3個基因芯片數據中均顯示為上調基因，見圖1。經GEPIA數據庫驗證，CDKN2A基因在宮頸鱗狀細胞癌組織組織中呈高表達水平(P<0.05)，見圖2。

圖1 3個數據集篩選出的差異性表達基因的交集

圖2 GEPIA數據庫驗證結果

2.2 KEGG通路富集分析結果 針對347個差異性表達基因進行KEGG通路富集分析，發現CDKN2A基因參與了細胞周期信號通路、p53信號通路及膀胱癌信號通路，見表1。

表1 CDKN2A基因參與的KEGG信號通路及其他參與基因

2.3 PPI網絡分析結果 利用Cytoscape軟件中Cytohubba程序進行分析，發現共有22個蛋白與CDKN2A基因編碼蛋白存在密切的相互作用，見圖3。

圖3 CDKN2A基因編碼蛋白與其他差異性表達基因編碼蛋白的相互作用

2.4 CDKN2A基因表達情況與宮頸癌患者臨床病理特征、預后的關系 與淋巴結轉移陰性宮頸癌患者相比，淋巴結轉移陽性宮頸癌患者的CDKN2A基因表達水平更高(P=0.035)，但在不同腫瘤分期、體重、年齡及腫瘤分化程度患者之間，CDKN2A基因表達水平差異均無統計學意義(均P>0.05)，見圖4。CDKN2A基因高表達水平的宮頸癌患者的生存期短于低表達者(P=0.024)，見圖5。

圖4 圖4 不同臨床病理特征宮頸癌患者之間 CDKN2A基因的表達情況

圖5 CDKN2A基因高表達水平與低表達水平宮頸癌患者的生存曲線比較

2.5 宮頸癌組織的CDKN2A基因啟動子區域甲基化水平 宮頸癌組織中CDKN2A基因啟動子區域甲基化水平較正常宮頸組織升高(P=0.008)，見圖6。

圖6 宮頸癌組織的CDKN2A基因啟動子區域甲基化水平

2.6 CDKN2A基因與腫瘤免疫細胞浸潤的相關性 CDKN2A基因的表達水平與巨噬細胞浸潤水平呈負相關，與中性粒細胞、樹突狀細胞浸潤水平均呈正相關(均P<0.05)，見圖7。

圖7 CDKN2A基因表達水平與巨噬細胞、中性粒細胞、樹突狀細胞浸潤水平的相關性注：橫坐標代表浸潤水平，縱坐標代表CDKN2A表達水平

2.7 上游miRNA的篩選結果 通過ENCORI和mirDIP軟件分析，分別得到90個和318個可能與CDKN2A基因有結合靶點的miRNA，取交集后獲得9個miRNA，分別為hsa-miR-4436a、hsa-miR-449a、hsa-miR-125a-5p、hsa-miR-125b-5p、hsa-miR-124-3p、hsa-miR-1286、hsa-miR-331-3p、hsa-miR-4319、hsa-miR-506-3p。

3 討 論

雖然近幾年宮頸癌的篩查和治療手段不斷發展，但宮頸癌仍是發病率最高的女性生殖道惡性腫瘤[7]，嚴重威脅女性的生命健康，宮頸癌新的診斷和治療方法仍有待進一步研究。與根治性手術、放射治療、化學治療等傳統的治療方法相比，宮頸癌基因靶向治療是一種新的治療方法，通過對宮頸癌患者實施基因治療能夠有效抑制宮頸癌細胞的生長[8]。因此，尋找與宮頸癌發生和發展相關的關鍵基因具有重要意義。我們根據GEO數據庫的3個宮頸癌基因芯片信息整合結果及GEPIA軟件的驗證結果，發現CDKN2A基因在宮頸癌組織中呈高表達；通過UALCAN在線工具對TCGA數據庫提供的數據進行分析后發現，宮頸癌組織中CDKN2A基因表達水平升高時宮頸癌患者淋巴結轉移可能性大，且患者的預后更差，這表明CDKN2A基因可能在宮頸癌的發生及發展過程中起著重要作用，是宮頸癌患者預后評估的潛在分子標志物。

研究表明，免疫細胞的浸潤與腫瘤的發生、發展息息相關，許多腫瘤不良預后相關基因的表達水平與腫瘤免疫細胞的浸潤有關[9-10]。免疫治療已成為宮頸癌治療的新手段。因此，我們利用TIMER數據庫進行分析，發現CDKN2A基因表達水平與腫瘤免疫細胞浸潤相關，其中與巨噬細胞浸潤水平呈負相關，與中性粒細胞、樹突狀細胞浸潤水平均呈正相關(均P<0.05)，這或可為研究宮頸癌的免疫治療提供新的線索。

本研究中，我們從多方面進一步探討CDKN2A基因的作用機制。蛋白質的功能主要是通過與其他蛋白質相互作用來實現的，因此我們構建了CDKN2A基因編碼蛋白的PPI網絡，發現CDKN2A基因編碼蛋白與多個基因編碼蛋白之間存在相互作用，其中已有研究報告CCNE1可以促進宮頸癌的進展[11]。通過研究此網絡中與CDKN2A關系密切的關鍵蛋白，或可為探索CDKN2A在宮頸癌中的作用機制提供新線索。同時，我們進行了KEGG信號通路分析，結果顯示CDKN2A基因參與了細胞周期信號通路、p53信號通路及膀胱癌信號通路。研究表明，CDKN2A基因是一個細胞周期依賴的蛋白，編碼細胞周期抑制劑蛋白p16和p14，可與 CDK4、CDK6 結合，形成具有激酶活性的復合物，從而阻斷該復合物對Rb蛋白磷酸化，導致細胞周期停止在G期，從而抑制細胞的增殖[12]。此外，DNA甲基化是癌癥發生和發展過程中常見的早期分子遺傳學變化。已有研究表明在多種腫瘤細胞中CDKN2A基因啟動子過度甲基化導致CDKN2A基因突變，從而影響腫瘤的進展，如胰腺癌、胃癌、食管癌等[3，5，13-15]，但在宮頸癌中暫無此類研究。因此，我們利用UALCAN在線工具提取TCGA數據庫中的數據進行分析，發現宮頸癌組織中CDKN2A基因的啟動子區域甲基化水平升高。這提示宮頸癌的發展可能受CDKN2A基因的啟動子區域甲基化水平的影響，但具體調控機制仍需進一步研究。

miRNA廣泛參與基因轉錄后水平的調控，是腫瘤領域研究的熱點，其主要通過靶向調控靶基因mRNA的3′末端非翻譯區參與細胞增殖、侵襲、遷移和凋亡等重要生物學過程[16]。本研究結果顯示，CDKN2A基因可能與9個miRNA存在結合靶點，推測CDKN2A基因可能受上游miRNA的調控發揮生物學功能。但這些miRNA與CDKN2A基因表達的相關性及結合靶點仍需進一步實驗研究驗證。

綜上所述，宮頸癌組織中 CDKN2A基因呈高表達水平，且與患者淋巴結轉移、疾病預后、腫瘤免疫細胞浸潤有關，其可作為宮頸癌診斷及預后評估的潛在分子標記物，并可為研究宮頸癌的免疫治療提供新線索；此外，CDKN2A基因可能通過啟動子區域甲基化和上游miRNA調控發揮生物學功能，這為進一步研究宮頸癌的診斷和治療方法提供了新思路。