《Biochemistry Research International》：Identifying Prognostic Biomarkers and Key Pathways in Renal Clear Cell Carcinoma: A Pilot Study Using Integrated miRNA and Gene Expression Analysis
編輯推薦：
本綜述通過整合miRNA與mRNA表達譜分析，系統鑒定了腎透明細胞癌（RCCC）中關鍵的預后生物標志物（如miR-26a-1-3p， CDK1， RUNX2）和失調的信號通路（如PI3K–Akt， MAPK）。研究揭示了miRNA-hub基因調控軸（如miR-26a-1-3p/RUNX2）在RCCC發生發展中的潛在作用，為理解其分子機制、開發預后工具及尋找新治療靶點提供了重要見解。
引言
腎透明細胞癌（RCCC）是腎細胞癌（RCC）中最常見且侵襲性強的亞型，約占所有RCC病例的85%，其高發病率、高死亡率及對常規治療反應不佳構成了重大的健康負擔。RCCC起源于近曲小管上皮細胞，其分子特征復雜，包括3p缺失和VHL腫瘤抑制基因失活等拷貝數改變。為了改善患者預后，識別準確的預后生物標志物和治療靶點至關重要。MicroRNAs（miRNAs）作為短鏈非編碼RNA，在基因轉錄后調控中扮演關鍵角色，其失調與腫瘤發生發展密切相關，因此成為癌癥生物標志物研究的 promising candidate。
材料與方法
本研究利用公共基因表達數據，對9例RCCC組織和11例正常腎組織進行了miRNA譜比較分析。通過GEO2R工具，以錯誤發現率（FDR）< 0.01和絕對log2倍數變化（|log2FC|）> 1.585為閾值，鑒定出差異表達miRNA（DEMs）。利用miRWalk 2.0和miRTarBase數據庫預測并驗證DEMs的靶基因，同時從DisGeNET數據庫獲取已知的RCCC相關基因。將兩組基因取交集，得到RCCC相關基因集。隨后，使用STRING數據庫構建蛋白互作網絡（PIM），并導入Cytoscape進行可視化與分析。通過計算度中心性和中介中心性，篩選出樞紐基因（hub genes）。利用MCODE插件識別網絡中的關鍵模塊（cluster）。功能富集分析通過g:Profiler工具進行，涵蓋基因本體（GO）術語（包括生物過程BP、分子功能MF、細胞組分CC）和KEGG等信號通路。預后分析方面，使用GEPIA2數據庫評估樞紐基因與RCCC患者總生存期（OS）的關系，并使用Kaplan–Meier plotter評估DEMs的預后價值。通過UALCAN數據庫驗證DEMs在腫瘤與正常組織中的表達差異。此外，利用GEO數據集GSE76351和TCGA–KIRC隊列（通過Kaplan–Meier plotter）獨立驗證關鍵樞紐基因（CDK1， RUNX2）的表達及預后意義。通過starBase數據庫分析關鍵miRNA與樞紐基因之間的表達相關性。最后，研究收集了5例腎癌和5例正常腎組織樣本，通過實時定量PCR（RT-PCR）驗證RUNX2的表達。
結果
1. RCCC中的差異表達miRNA
分析共鑒定出15個DEMs，其中9個上調，6個下調。上調的miRNA包括hsa-miR-542-5p（FC=20.535）、hsa-miR-144（FC=4.469）等；下調的miRNA包括hsa-miR-26a-1（FC=0.301）、hsa-miR-377（FC=0.089）等。
2. 模塊、樞紐基因、通路與GO術語
通過數據庫交集獲得了424個RCCC相關基因。基于此構建的PIM網絡包含396個蛋白和3704個互作。從中篩選出74個樞紐基因，其中度中心性最高的包括TP53、AKT1、MYC、EGFR、SRC等。MCODE分析識別出三個重要的基因子網絡（Cluster 1-3），其中Cluster 1包含56個基因，互作邊數最多（522條）。
通路富集分析顯示，“癌癥通路”（Pathways in cancer）、“PI3K–Akt信號通路”（PI3K–Akt signaling pathway）和“MAPK信號通路”（MAPK signaling pathway）是RCCC惡性轉化中最顯著失調的信號通路。在生物過程（BP）層面，“程序性細胞死亡調控”（Regulation of programmed cell death）、“信號轉導調控”（Regulation of signal transduction）和“凋亡過程調控”（Regulation of apoptotic process）最為富集。細胞組分（CC）中，“細胞質”（Cytoplasm）和“細胞外圍”（Cell periphery）受影響最大。分子功能（MF）方面，“蛋白結合”（Protein binding）和“酶結合”（Enzyme binding）最為顯著。
3. 預后標志物
生存分析表明，CRP、CDK1和RUNX2的過表達是RCCC的不良預后標志物。其中CRP的風險比（HR）最高，為1.8。這三個基因的組合（CRP+CDK1+RUNX2）也顯示出預后價值（HR=1.7）。相反，ERBB2、IGF1R等38個基因被鑒定為有利的預后標志物。對于miRNA，hsa-miR-26a-1-3p、hsa-miR-144-3p和hsa-miR-144-5p的低表達與腎透明細胞癌（KIRC）患者的不良預后顯著相關。UALCAN數據庫分析獨立驗證了miR-26a-1在腎癌組織中顯著下調，而hsa-miR-144顯著上調。
4. 基因表達模式與分期圖分析
在523例KIRC組織與100例正常腎組織的比較中，CDK1和RUNX2的mRNA在腫瘤組織中顯著上調，而CRP的mRNA水平未見顯著升高。分期圖分析進一步顯示，CDK1和RUNX2的表達在不同腫瘤分期（I-IV期）間存在顯著差異，提示它們可能與疾病進展相關；而CRP的表達在各分期間無顯著差異。
5. RT-PCR驗證
對5例腎癌組織和5例正常腎組織的RT-PCR分析證實，RUNX2在腎癌組織中的表達水平顯著高于正常組織。
6. RUNX2和CDK1表達與預后意義的獨立驗證
在獨立的GEO數據集GSE76351中，RUNX2（FC=1.89）和CDK1（FC=1.88）在RCCC組織中均顯著過表達。基于TCGA–KIRC隊列RNA-Seq數據的生存分析進一步證實，RUNX2（HR=1.93）和CDK1（HR=2.26）的高表達與患者總生存期縮短顯著相關。
7. miRNA-靶基因表達相關性分析
在TCGA–KIRC隊列中，hsa-miR-26a-1-3p與RUNX2表達呈微弱但顯著的負相關（r = -0.153）。hsa-miR-144-3p與RUNX2和CDK1也顯示出微弱的負相關。這些結果表明了潛在的miRNA-mRNA調控關系，但效應值較小，提示調控可能是多因素和情境依賴的。
討論
本研究鑒定的三個miRNA（miR-26a-1-3p， miR-144-3p， miR-144-5p）和三個樞紐基因（CRP， CDK1， RUNX2）在RCCC的預后中具有重要價值。已有研究支持這些發現：例如，miR-26a-1-3p在轉移性RCC中顯示出預后價值；miR-144-3p可通過與長鏈非編碼RNA NORAD相互作用影響MYCN表達；血清和腫瘤內CRP水平與RCCC不良預后相關；CDK1作為細胞周期關鍵激酶，其活性與RCCC增殖和遷移密切相關；RUNX2則通過與SCD1、Zic2、MAPK11等分子相互作用，促進RCCC的惡性進展。此外，本研究富集出的PI3K–Akt通路在RCCC中被廣泛報道，靶向該通路可能具有治療潛力。當然，本研究也存在一些局限性，如初始分析樣本量較小、部分標志物未進行實驗驗證、鑒定的miRNA-mRNA相關性需要進一步的功能實驗確認等。
結論
本研究通過整合生物信息學分析和多層面驗證，系統揭示了RCCC中15個DEMs和75個樞紐基因的失調。特別指出，miR-26a-1-3p、miR-144-3p和miR-144-5p的低表達，以及CDK1和RUNX2的高表達，是RCCC強有力的預后標志物。這些分子通過影響PI3K–Akt、MAPK等關鍵信號通路，參與調控細胞凋亡和細胞周期進程，共同驅動RCCC的發生發展。這些發現為深入理解RCCC的分子機制、開發新的預后評估工具以及探索個性化治療策略提供了重要的理論依據和潛在的干預靶點。
