FJX1基因在結(jié)腸癌中的臨床意義

薛 麗，黃 奔

南京醫(yī)科大學第一附屬醫(yī)院檢驗學部，江蘇南京 210009

人類四連接同源 (FJX1) 基因與果蠅高爾基體固有跨膜激酶(Fj) 基因高度同源，位于11號染色體(11p13)，可以編碼含有437個氨基酸的表面分泌蛋白[1]。FJX1基因可以在神經(jīng)系統(tǒng)、多個器官的上皮細胞和肢體發(fā)育過程中表達，對生物體的發(fā)育起著重要的作用[2-3]。有研究指出，F(xiàn)JX1 基因可以作為致癌基因，通過多種信號途徑影響腫瘤的發(fā)生、發(fā)展。FJX1 基因被報道在鼻咽癌[3]、黑色素瘤[4]、頭頸癌[5]等多種腫瘤中發(fā)揮了重要的作用。雖然有研究指出，F(xiàn)JX1 mRNA和蛋白在晚期結(jié)直腸癌上皮細胞中表達上調(diào)[6]，但研究者并未探討FJX1基因表達與結(jié)腸癌患者臨床病理參數(shù)、預后的關(guān)系。因此，本研究通過數(shù)據(jù)挖掘癌癥基因組圖譜(TCGA)數(shù)據(jù)庫，利用生物信息學的方法，綜合多種數(shù)據(jù)庫加以驗證，探討FJX1 基因在結(jié)腸癌中的臨床意義、預后價值以及FJX1 基因在結(jié)腸癌中的分子作用機制，旨在為結(jié)腸癌的診治提供新型分子標志物。

1 材料與方法

1.1數(shù)據(jù)下載與處理 登陸TCGA官方數(shù)據(jù)庫(https://portal.gdc.cancer.gov/)，下載結(jié)腸癌轉(zhuǎn)錄組測序數(shù)據(jù)(level 3,HTSeq-FPKM) 和患者的臨床資料?；颊叩呐R床數(shù)據(jù)包括年齡、性別、腫瘤TNM分期、生存狀態(tài)以及隨訪時間。共下載結(jié)腸癌組織標本473例、正常組織標本41例。在后續(xù)分析中，刪除臨床數(shù)據(jù)不完整或者缺失的標本86例，納入387例結(jié)腸癌組織標本，其中男性患者203例，女性患者184例；年齡>65歲患者235例，≤65歲患者 152例；Stage分期Ⅰ～Ⅱ期患者224例，Ⅲ～Ⅳ期患者163例；存活患者291例，死亡患者96例；最長隨訪時間為4 502 d。使用R 4.1.2軟件提取基因表達數(shù)據(jù)并進行歸一化處理。在后續(xù)分析中，根據(jù)FJX1基因在結(jié)腸癌患者中表達水平的中位值，將患者分為高表達(>中位值)、低表達(≤中位值)。

1.2FJX1基因在結(jié)腸癌中的表達以及與臨床特征關(guān)系的分析 根據(jù)TCGA結(jié)腸癌數(shù)據(jù)集的臨床資料，分析FJX1基因在結(jié)腸癌組織中的表達情況，探討FJX1基因表達與結(jié)腸癌患者臨床特征(年齡、性別、腫瘤分期)的關(guān)系。同時，使用UALCAN數(shù)據(jù)庫(http://ualcan.path.uab.edu/index.html) 驗證分析FJX1基因在結(jié)腸癌組織中的表達情況。在UALCAN數(shù)據(jù)庫中，共下載結(jié)腸癌組織標本286例、正常組織標本41例。在分析TP53基因突變的時候，刪除4例信息缺失的標本。在后續(xù)分析中，刪除臨床數(shù)據(jù)不完整或者缺失的標本7例，納入279例結(jié)腸癌組織標本，其中男性患者156例，女性患者123例；年齡>65歲患者177例，≤65歲患者102例；Stage分期Ⅰ～Ⅱ期患者155例，Ⅲ～Ⅳ期患者124例；存活患者183例，死亡患者96例；最長隨訪時間為4 000 d。

1.3FJX1基因表達與結(jié)腸癌患者預后關(guān)系的分析 根據(jù)結(jié)腸癌患者的生存資料，結(jié)合FJX1基因的表達情況，分析高、低表達患者總體生存率的差異。在本研究中，利用UALCAN數(shù)據(jù)庫、GEPIA2數(shù)據(jù)庫(http://gepia2.cancer-pku.cn/)和OncoLnc在線數(shù)據(jù)庫(http://www.oncolnc.org) 驗證分析FJX1基因表達與患者總體生存率的關(guān)系。其中，在GEPIA2數(shù)據(jù)庫中，共下載結(jié)腸癌組織標本275例、正常組織標本41例。275例結(jié)腸癌組織標本中男性患者192例，女性患者83例；年齡>65歲患者165例，≤65歲患者110例；Stage分期Ⅰ～Ⅱ期患者149例，Ⅲ～Ⅳ期患者126例；存活患者165例，死亡患者110例；最長隨訪時間為5 400 d。在OncoLnc數(shù)據(jù)庫中，共下載結(jié)腸癌組織標本440例，正常組織標本41例。440例結(jié)腸癌組織標本中男性患者239例，女性患者201例；年齡>65歲患者215例，≤65歲患者225例；Stage分期Ⅰ～Ⅱ期患者239例，Ⅲ～Ⅳ期患者201例；存活患者342例，死亡患者98例；最長隨訪時間為4 502 d。

1.4基因富集分析 使用R4.2.2軟件對FJX1基因進行功能富集分析，探討FJX1基因在結(jié)腸癌中的分子作用機制。采用加權(quán)富集法分析，分析的通路來自Molecular Signature Database (Msig-DB)數(shù)據(jù)庫中的 c2.cp.kegg.v7.0.symbols.gmt 數(shù)據(jù)集。參數(shù)設置隨機組合次數(shù)1 000次。錯誤發(fā)現(xiàn)率(FDR)小于0.05和P<0.05的基因集作為顯著富集的基因集。此外，為了充分研究FJX1基因的生物學功能，使用STRING數(shù)據(jù)庫 (https://cn.string-db.org/) 探索與FJX1基因相互作用的基因網(wǎng)絡，對網(wǎng)絡中的基因進行基因本體 (GO)富集分析，探討其參與的生物學過程 (BP)、細胞學組分(CC)和分子功能(MF)。

1.5統(tǒng)計學處理 使用SPSS Statistics 25.0統(tǒng)計學軟件及R 4.2.2軟件進行統(tǒng)計學分析。采用Mann-WhitneyU檢驗比較FJX1基因在癌組織和正常組織中的表達差異，采用χ2檢驗分析在患者不同臨床病理參數(shù)下FJX1基因的表達差異，采用Kaplan-Meier法和Cox回歸分析探討FJX1基因表達對預后的判斷價值。以P<0.05為差異有統(tǒng)計學意義。

2 結(jié) 果

2.1FJX1基因在結(jié)腸癌組織中的表達 在TCGA結(jié)腸癌數(shù)據(jù)集中，相比于正常組織(FJX1基因表達水平中位值為0.261)，F(xiàn)JX1基因在結(jié)腸癌組織中的表達(中位值為1.761)顯著升高，差異有統(tǒng)計學意義(P<0.001),見圖1A。與相應的癌旁組織(中位值為0.261)相比，癌組織中的FJX1基因表達水平(中位值為1.960)顯著升高，差異有統(tǒng)計學意義(P<0.001),見圖1B。此外，在UALCAN數(shù)據(jù)庫中，結(jié)果也顯示FJX1基因在結(jié)腸癌組織表達(中位值為2.150)顯著高于正常組織(中位值為0.350)，差異有統(tǒng)計學意義(P<0.001),見圖1C；與未發(fā)生TP53基因突變的組織標本(中位值為1.939)相比，發(fā)生TP53基因突變的組織標本FJX1基因表達水平(中位值為2.255)顯著升高，差異具有統(tǒng)計學意義(P=0.031)，見圖1D。

注：A為在TCGA結(jié)腸癌數(shù)據(jù)集中FJX1基因表達情況；B為在TCGA結(jié)腸癌數(shù)據(jù)集中癌組織和相應癌旁組織中FJX1基因表達情況;C為在UALCAN數(shù)據(jù)庫中FJX1基因表達情況；D為在TCGA結(jié)腸癌數(shù)據(jù)集中具有TP53基因突變標本與TP53基因未突變標本之間FJX1基因表達情況。

2.2FJX1基因表達與結(jié)腸癌患者臨床特征的關(guān)系 根據(jù)FJX1基因表達水平中位值1.761，將結(jié)腸癌患者分為高、低表達。結(jié)果顯示，F(xiàn)JX1基因表達水平在患者不同Stage分期、M分期中比較，差異均有統(tǒng)計學意義(均P=0.038)，見表1。

表1 FJX1基因表達與結(jié)腸癌患者臨床特征的關(guān)系

在TCGA結(jié)腸癌數(shù)據(jù)集中，F(xiàn)JX1基因在不同分期的患者中表達情況如圖2所示。其中，Stage分期Ⅳ期患者FJX1基因的表達量顯著高于Stage分期Ⅰ期的患者，差異具有統(tǒng)計學意義(P=0.023)；M1期患者的FJX1基因的表達量顯著高于M0期的患者，差異具有統(tǒng)計學意義 (P= 0.044)。

注：A表示FJX1基因在患者不同Stage分期中的表達情況；B表示FJX1基因在T分期中的表達情況；C表示FJX1基因在M分期中的表達情況;D表示FJX1基因在N分期中的表達情況。

2.3FJX1基因表達與結(jié)腸癌患者預后的關(guān)系 Kaplan-Meier生存分析結(jié)果顯示，高表達FJX1基因組的患者總體生存率顯著低于低表達的患者(P=0.007)，見圖3A。此外，在UALCAN數(shù)據(jù)庫、GEPIA2數(shù)據(jù)庫和OncoLnc數(shù)據(jù)庫中，均顯示了高表達FJX1基因組的患者預后不良 (P<0.05)，見圖3B～D。

注：A為在TCGA結(jié)腸癌數(shù)據(jù)集中FJX1基因表達與結(jié)腸癌患者總體生存率的關(guān)系；B為UALCAN數(shù)據(jù)庫中FJX1基因表達與結(jié)腸癌患者總體生存率的關(guān)系；C為GEPIA2數(shù)據(jù)庫中FJX1基因表達與結(jié)腸癌患者總體生存率的關(guān)系；D為在OncoLnc數(shù)據(jù)庫中FJX1基因表達與結(jié)腸癌患者總體生存率的關(guān)系。

2.4Cox回歸分析 單因素Cox回歸分析結(jié)果顯示，年齡、TNM分期和FJX1基因表達均在結(jié)腸癌不同預后患者中差異有統(tǒng)計學意義(P<0.05)，見表2。而多因素Cox回歸分析結(jié)果顯示，年齡(HR=1.043,95%CI:1.022～1.065,P<0.001)、T分期(HR=1.981,95%CI:1.177～3.336,P=0.010)、FJX1基因表達(HR=1.133,95%CI:1.025～1.253,P=0.015)可以作為患者預后的獨立影響因素。

表2 Cox回歸分析

2.5基因富集分析 基因富集分析結(jié)果顯示，F(xiàn)JX1基因主要參與了致癌相關(guān)信號途徑、細胞黏附分子、趨化因子信號通路、補體和凝血級聯(lián)、細胞因子受體相互作用、細胞外基質(zhì)(ECM)受體相互作用、Hedgehog(Hh)信號通路等分子作用通路。使用STRING數(shù)據(jù)庫探索與FJX1基因相互作用的基因網(wǎng)絡，結(jié)果顯示，共有10個基因與FJX1基因互相作用 (圖4A)。GO富集分析結(jié)果顯示，這些基因與FJX1基因共同參與了細胞黏附、細胞因子受體作用、癌癥信號途徑等生物學過程；參與了細胞頂端、細胞連接、細胞外泌體等細胞學組分；參與了鈣離子、蛋白激酶激活、受體配體作用等分子功能 (圖4B)。因此，F(xiàn)JX1基因可以通過細胞黏附、細胞因子受體作用、ECM受體作用以及癌癥相關(guān)信號途徑發(fā)揮其生物學功能。

注：A圖為STRING數(shù)據(jù)庫中與FJX1基因相互作用的基因網(wǎng)絡圖；FAT1～4即脂肪非典型鈣黏蛋白1～4;DCHS1～2即鈣依賴性細胞黏附蛋白;COMMD9即銅代謝結(jié)構(gòu)域蛋白-9;SYNGR3即神經(jīng)細胞囊泡循環(huán)蛋白3;TRIM44即含三聯(lián)基元蛋白44;FAM20A即序列相似家族20成員A;B圖為GO富集分析結(jié)果圖。

3 討 論

人類 FJX1基因與在果蠅和小鼠中的直系同源物具有 29% 和 88% 的同源性[6]。在果蠅體內(nèi)，F(xiàn)j是生長發(fā)育和平面細胞極性調(diào)節(jié)所必需的[7]，而在小鼠體內(nèi)，F(xiàn)JX1基因存在于發(fā)育中的大腦中，純合型FJX1基因小鼠在大腦海馬中會表現(xiàn)出異常形態(tài)的樹突結(jié)構(gòu)，在體外培養(yǎng)時，這些小鼠的海馬神經(jīng)元會表現(xiàn)出更高水平的樹突延伸和分支[8]。因此，這些研究強調(diào)了FJX1基因在不同生物體發(fā)育中的重要性。盡管FJX1基因在人類中的確切生物學功能尚不完全清楚，但來自發(fā)育研究的數(shù)據(jù)表明，F(xiàn)JX1基因參與調(diào)節(jié)細胞的生長和分化，在多種人類疾病中發(fā)揮了重要的作用[9]，F(xiàn)JX1基因在癌癥發(fā)展中的作用可能與其在促進生長和細胞運動中的作用有關(guān)[3]。目前許多研究已經(jīng)證實FJX1基因在人類癌癥中過表達，如乳腺癌[10]、肺癌[11]等，這表明了FJX1基因可能是癌癥的分子驅(qū)動因素。

在鼻咽癌組織中，F(xiàn)JX1基因和蛋白水平均表達上調(diào)，促進了癌細胞增殖和侵襲[3]。與正常卵巢內(nèi)皮細胞相比，F(xiàn)JX1 mRNA在卵巢腫瘤內(nèi)皮細胞中表達上調(diào)[12-13]。CHANG等[14]的研究發(fā)現(xiàn)，與無子宮內(nèi)膜異位癥的女性相比，F(xiàn)JX1基因表達水平在異位子宮內(nèi)膜分泌期顯著升高。在結(jié)直腸癌中，CUI等[15]研究者構(gòu)建了miRNA-mRNA調(diào)節(jié)網(wǎng)絡以揭示結(jié)直腸癌的發(fā)病機制，并指出了FJX1等基因可以作為結(jié)直腸癌的關(guān)鍵基因。然而，關(guān)于FJX1基因在結(jié)腸癌中的臨床意義和預后作用，目前鮮有研究報道。

本研究通過數(shù)據(jù)分析發(fā)現(xiàn)，F(xiàn)JX1基因在結(jié)腸癌組織中表達上調(diào)，并且與患者的Stage分期和遠處轉(zhuǎn)移相關(guān)，這表明FJX1基因表達參與了結(jié)腸癌患者的惡性進展。高表達FJX1基因患者總體生存率較低，這表明FJX1基因影響患者的生存，可以作為結(jié)腸癌潛在的預后分子標志物。有研究指出，在結(jié)直腸癌組織中，F(xiàn)JX1基因的表達與血管生成基因如低氧誘導因子-1α(HIF-1α)、血管內(nèi)皮生長因子受體、血管內(nèi)皮生長因子C的變化顯著相關(guān)，F(xiàn)JX1基因可以通過HIF-1α依賴性方式促進內(nèi)皮細胞毛細管的形成[6]。因此，研究者認為FJX1基因是一種新型的腫瘤進展調(diào)節(jié)因子，F(xiàn)JX1基因增強血管生成的能力可能解釋了FJX1基因高表達的患者總體生存率較低。

基因富集分析結(jié)果顯示，F(xiàn)JX1基因表達參與了多種信號途徑，如調(diào)控細胞分化與增殖、ECM受體途徑[16]，與腫瘤轉(zhuǎn)移和惡性進展相關(guān)的通路[17]，與腫瘤發(fā)生、侵襲和轉(zhuǎn)移相關(guān)的細胞因子途徑[18-19]等。分子機制研究指出，F(xiàn)JX1基因的過表達可以通過上述多種途徑參與腫瘤的發(fā)生與發(fā)展。Hippo通路在調(diào)節(jié)細胞增殖和凋亡方面發(fā)揮著關(guān)鍵作用，F(xiàn)JX1基因是該通路中YAP和TAZ癌基因的直接靶標[20]，F(xiàn)JX1蛋白過表達還可以進一步反饋調(diào)節(jié)YAP的激活，參與癌細胞的增殖途徑[21]。FJX1基因過表達后細胞周期蛋白D1(cyclin D1)和E1(cyclin E1)的表達水平也相應增加，這表明FJX1基因可以協(xié)同調(diào)節(jié)細胞周期的關(guān)鍵蛋白來增強癌細胞增殖[3]。FJX1基因還可以通過參與調(diào)節(jié)平面細胞極性(PCP)來增強癌細胞的侵襲，雖然PCP是傷口修復和發(fā)育的關(guān)鍵過程，它參與促進了原腸形成、神經(jīng)管閉合等過程，該生物學過程的失調(diào)可能將導致癌細胞產(chǎn)生侵襲和轉(zhuǎn)移的行為[3,22]。此外，F(xiàn)JX1基因表達還可以通過促進血管生成素等促血管生成基因的表達，誘導血管生成肽的產(chǎn)生來增強癌癥進展[23]。

由于FJX1基因的關(guān)鍵致癌作用，靶向FJX1基因被認為是一種潛在的結(jié)直腸癌治療策略。LIU等[24]的研究發(fā)現(xiàn)，抑制非編碼RNA PVT1的表達可以通過靶向調(diào)控miR-106b-5p/FJX1信號軸，抑制結(jié)直腸癌細胞的增殖、遷移和侵襲，為結(jié)直腸癌患者治療策略的發(fā)展提供了新思路。此外，miR-532-3p[25]、miR-1249[26]均可以靶向調(diào)控FJX1基因表達，抑制結(jié)腸腺癌細胞的增殖和侵襲。因此可以推斷，F(xiàn)JX1基因有希望成為結(jié)直腸癌潛在的分子標志物。

然而，本研究也存在一定的局限性。首先，本研究是通過生物信息學的方法，綜合多種數(shù)據(jù)庫進行驗證分析FJX1基因在結(jié)腸癌中的表達及對患者總體生存率的影響，初步揭示了FJX1基因在結(jié)腸癌中的分子意義。在分析患者的臨床數(shù)據(jù)時，雖然排除了臨床信息不完整、生存數(shù)據(jù)缺失的患者資料，但由于每位患者的情況不同，接受的臨床治療也存在差異，本研究是利用患者的隨訪數(shù)據(jù)初步顯示了FJX1基因?qū)Y(jié)直腸癌預后的判斷作用，結(jié)果可能存在偏倚。因此，在后續(xù)研究中仍需大標本的醫(yī)學實驗進一步驗證。其次，通過基因富集分析指出了FJX1基因在致癌過程中的潛在分子作用機制，仍需要進行生物體內(nèi)的醫(yī)學實驗進行驗證。

綜上所述，本研究顯示FJX1基因在結(jié)腸癌組織中異常表達上調(diào)，與患者的惡性進展特征和不良預后相關(guān)。FJX1基因通過參與調(diào)節(jié)細胞黏附、細胞因子受體作用、ECM受體作用以及癌癥相關(guān)信號途徑促進了腫瘤的進展。FJX1基因可作為結(jié)腸癌患者的分子預后標志物和潛在的治療靶點，這為進一步探討FJX1基因在結(jié)腸癌中的生物學功能提供了線索和依據(jù)。