賈晨飛，陳恩立，宋惠玲，司華蕊，王娟

1河北省人民醫(yī)院，石家莊050051；2石家莊市人民醫(yī)院

食管癌屬于原發(fā)于食管黏膜上皮的惡性腫瘤，是臨床最常見的消化道惡性腫瘤之一，在我國的發(fā)病率和病死率均位居前列［1］，由于其早期臨床特征不明顯，很多患者確診時(shí)已是中晚期。即便進(jìn)行標(biāo)準(zhǔn)的同步放化療序貫手術(shù)切除后，仍有較高的復(fù)發(fā)和轉(zhuǎn)移率，5年生存率不足30%［2］。目前食管癌發(fā)生發(fā)展的分子機(jī)制還不清楚。miRNA是一種小分子非編碼RNA，作用于mRNA的降解或抑制從而影響基因表達(dá)［3］，多項(xiàng)研究表明其可能參與了腫瘤的發(fā)生發(fā)展，是一種具有廣闊前景的生物標(biāo)記物，但與食管癌發(fā)生發(fā)展相關(guān)miRNA的研究少見，2020年10—12月，我們通過下載GEO數(shù)據(jù)庫中食管鱗癌與正常組織的miRNA和mRNA表達(dá)譜芯片，構(gòu)建miR‐NA-mRNA調(diào)控網(wǎng)絡(luò)，尋找食管癌發(fā)生發(fā)展的相關(guān)分子標(biāo)記物，為其分子機(jī)制的闡明提供實(shí)驗(yàn)依據(jù)。

1 材料與方法

1.1 食管鱗癌相關(guān)基因芯片來源GEO數(shù)據(jù)庫（https：//www.ncbi.nlm.nih.gov/geo）。共下載兩張組織芯片（GSE75241、GSE114110）。GSE75241為mRNA表達(dá)數(shù)據(jù)，包括15個(gè)食管鱗癌樣本和15個(gè)配對(duì)正常食管黏膜組織樣本，檢測平臺(tái)為GPL5175。GSE114110為miRNA表達(dá)數(shù)據(jù)，包括30個(gè)食管鱗癌樣本和10個(gè)正常組織樣本，檢測平臺(tái)為GPL5175。

1.2 食管鱗癌組織芯片中差異表達(dá)miRNA及差異表達(dá)mRNA的篩選通過與正常食管黏膜組織標(biāo)本相比較，利用R語言limma函數(shù)包篩選出食管鱗癌組織樣本中差異表達(dá)的mRNAs和miRNAs。篩選閾值為P＜0.05和|log2FC|＞2（logFC＞0表示上調(diào)差異表達(dá)基因，logFC＜0表示下調(diào)差異表達(dá)基因）。

1.3 食管鱗癌組織樣本miRNA-mRNA調(diào)控網(wǎng)絡(luò)的構(gòu)建利用FunRich3.1.3軟件預(yù)測食管鱗癌組織樣本差異表達(dá)miRNA調(diào)控的mRNA，剔除非差異表達(dá)mRNA，之后利用Cytoscape軟件構(gòu)建miR‐NA-mRNA調(diào)控網(wǎng)絡(luò)圖，并且篩選出節(jié)點(diǎn)度值為4～8的miRNA（節(jié)點(diǎn)度大的RNA為調(diào)控網(wǎng)絡(luò)中參與miRNA-mRNA關(guān)系對(duì)較多的RNA）。

1.4 miRNA-mRNA調(diào)控網(wǎng)絡(luò)中mRNA生物學(xué)功能分析利用DAVID在線工具（https：//david.ncifcrf.gov）對(duì)調(diào)控網(wǎng)絡(luò)中的mRNA進(jìn)行基因本體論分析（GO富集）和京都基因與基因組百科全書分析（KEGG通路分析），篩選標(biāo)準(zhǔn)為P＜0.05。其中GO富集包括生物過程（BP）、細(xì)胞組分（CC）、分子功能（MF）三部分。

2 結(jié)果

2.1 食管鱗癌組織中差異表達(dá)的miRNAs與正常食管黏膜組織標(biāo)本相比，食管鱗癌組織樣本中差異表達(dá)miRNAs有114個(gè)，其中42個(gè)表達(dá)下調(diào)、72個(gè)表達(dá)上調(diào)；差異表達(dá)mRNAs 224個(gè)，其中124個(gè)表達(dá)上調(diào)、100個(gè)表達(dá)下調(diào)。（見表1-2）。

表1 食管鱗癌組織樣本中差異表達(dá)的前30個(gè)miRNAs

2.2 食管鱗癌組織樣本中差異表達(dá)miRNA-mRNA調(diào)控網(wǎng)絡(luò)食管鱗癌組織樣本中差異mRNA和差異miRNA共同參與的miRNA-mRNA關(guān)系對(duì)66個(gè)，其中具有較高節(jié)點(diǎn)度的miRNA為hsa-miR-92b-3p、hsa-miR-195-5p、hsa-miR-143-3p，mRNA為HAS2、ITGA6、COL1A1等（節(jié)點(diǎn)度范圍4～8），見圖1。

表2 食管鱗癌組織樣本中差異表達(dá)明顯的前30個(gè)mRNAs

2.3 miRNA-mRNA調(diào)控網(wǎng)絡(luò)中mRNA的生物學(xué)功能miRNA-mRNA調(diào)控網(wǎng)絡(luò)中的mRNAs富集于211個(gè)GO terms，包 括177個(gè)BP、19個(gè)CC和15個(gè)MF。主要富集在細(xì)胞外基質(zhì)、成骨作用、維生素應(yīng)答、血小板α顆粒等生物學(xué)功能。另外富集在ECM-受體相互作用、黏著斑、PI3K-Akt、TGF-β等13條通路中。

3 討論

圖1 食管鱗癌組織樣本miRNA-mRNA調(diào)控網(wǎng)絡(luò)圖（三角形代表miRNA，圓形代表mRNA）

miRNA是一種小分子非編碼RNA，具有多個(gè)生物學(xué)功能，近年研究表明，其可能在基因表達(dá)的全過程發(fā)揮作用，在腫瘤細(xì)胞侵襲轉(zhuǎn)移中也起關(guān)鍵作用，但不同的遺傳背景和細(xì)胞狀態(tài)都會(huì)影響miRNA對(duì)腫瘤細(xì)胞運(yùn)動(dòng)的調(diào)控。Hsa-miR-92b通過直接抑制p57蛋白的表達(dá)對(duì)胚胎干細(xì)胞G1期實(shí)現(xiàn)調(diào)控，已有Hsa-miR-92b在神經(jīng)系統(tǒng)腫瘤、淋巴瘤、非小細(xì)胞肺癌、肝癌、食管癌等惡性腫瘤中影響細(xì)胞增殖的報(bào)道［4］。MA等［5］發(fā)現(xiàn)，不伴有淋巴結(jié)轉(zhuǎn)移的食管癌樣本中Hsa-miR-92b高表達(dá)，且預(yù)后好。本研究中Hsa-miR-92b是miRNA-mRNA調(diào)控網(wǎng)絡(luò)中節(jié)點(diǎn)度最高的miRNA，與癌旁組織相比，其表達(dá)顯著下調(diào)，說明在抑制食管癌侵襲轉(zhuǎn)移過程中起著關(guān)鍵作用。

本研究結(jié)果還顯示hsa-miR-195-5p、hsa-miR-143-3p可能參與了食管鱗癌的發(fā)生發(fā)展。miR-195隸屬miR-15/107miRNA家族，目前已知在惡性腫瘤、心力衰竭和精神分裂癥等多種疾病中被誘導(dǎo)和激活。其通過靶向于多種蛋白的mRNA作為癌基因或腫瘤抑制因子，在包括消化系統(tǒng)腫瘤、生殖系統(tǒng)腫瘤、骨腫瘤等多種腫瘤中發(fā)揮作用。已有文獻(xiàn)表明，miR195在口腔鱗癌組織中表達(dá)低于正常組織，其因?qū)DC42（已被證明在多種人類癌癥中過表達(dá)）的抑制實(shí)現(xiàn)抗增殖和抗侵襲作用［6］。本研究結(jié)果顯示，食管癌樣本中miRNA-195表達(dá)下調(diào)，提示其在食管癌中也可能發(fā)揮相同作用。SUN等［7］認(rèn)為，miR-195/CDC42的聯(lián)合表達(dá)比兩者單獨(dú)表達(dá)更有利于作為預(yù)后標(biāo)記。miR-195的行為機(jī)制十分復(fù)雜，參與包括Rb-E2F、PI3K/AKT、MAPK/REK在內(nèi)的多種信號(hào)通路，雖然其在腫瘤發(fā)生發(fā)展復(fù)雜機(jī)制中作用的研究還處在早期階段，但作為食管鱗癌患者生物標(biāo)記物的潛力值得期待。hsa-miR-143是多種癌癥的抑瘤因子，包括miR143-3p和miR143-5p兩種亞型，其中miR-143-3p是最著名的抑癌miRNA之一，在乳腺癌、胃癌、結(jié)直腸癌、膀胱癌、腎癌等多種癌癥中均有下調(diào)。有研究認(rèn)為，HOTAIR通過抑制miR-143和miR-125來下調(diào)HK2［8］，從而增強(qiáng)了食管癌細(xì)胞的侵襲性和遷移，本研究結(jié)果與文獻(xiàn)報(bào)道一致。

透明質(zhì)酸（HA）是ECM的重要組成部分，在腫瘤發(fā)生的各個(gè)階段（包括腫瘤起始、轉(zhuǎn)移以及化療耐藥性的獲得）都是關(guān)鍵因素［9］。HAS2 mRNA的表達(dá)與HA的產(chǎn)生密切相關(guān)，已有研究證明，HAS2表達(dá)的增加與乳腺癌、子宮內(nèi)膜癌和星形細(xì)胞瘤的進(jìn)展和轉(zhuǎn)移相關(guān)，其不僅在HA代謝和信號(hào)通路中起協(xié)調(diào)作用，而且在miRNAs的作用中還起一系列其他通路的作用［10］。本研究結(jié)果顯示，食管鱗癌樣本中HAS2作為表達(dá)上調(diào)的差異mRNA，在miRNA-mRNA調(diào)控網(wǎng)絡(luò)中節(jié)點(diǎn)度較高，表示HAS2可能在食管鱗癌的增殖、轉(zhuǎn)移中起到了促進(jìn)作用。整合素α6（ITGA6）常被用于識(shí)別各種組織中的癌癥干細(xì)胞群體，它是一種主要的層粘連蛋白受體，研究證明其在胰臟癌、乳腺癌和食管癌中可增強(qiáng)細(xì)胞增殖和侵襲能力［11］。本研究結(jié)果顯示，食管鱗癌樣本中ITGA6表達(dá)顯著上調(diào)，這與既往研究一致。值得一提的是，ITGA6是與細(xì)胞外配體相互結(jié)合的細(xì)胞表面受體［12］，基于這種特性，其具有良好的藥理學(xué)前景，另外，Hsa-miR-92b可直接錨定于ITGA6［5］，說明調(diào)節(jié)Hsa-miR-92b有可能可以治療食管癌。

細(xì)胞外基質(zhì)由基底膜和細(xì)胞間基質(zhì)構(gòu)成，是復(fù)雜的大分子網(wǎng)絡(luò)，主要成分是糖胺聚糖、蛋白聚糖和和纖維蛋白等［13］。細(xì)胞外基質(zhì)的降解是惡性腫瘤侵襲與轉(zhuǎn)移不可或缺的環(huán)節(jié)，基質(zhì)金屬蛋白酶（MMPs）是目前發(fā)現(xiàn)這一環(huán)節(jié)最重要的酶類，研究表明有淋巴結(jié)轉(zhuǎn)移的食管癌患者中MMP-9蛋白表達(dá)高于對(duì)照組［14］。另外MMPs也具有促進(jìn)腫瘤血管生長的作用。本研究結(jié)果顯示，調(diào)控網(wǎng)絡(luò)中差異mRNA主要富集在細(xì)胞外基質(zhì)的成分構(gòu)成、參與受體間相互作用的通路、腫瘤相關(guān)蛋白聚糖通路等，證明細(xì)胞外基質(zhì)與食管癌發(fā)生發(fā)展密切相關(guān)。除此之外，差異表達(dá)的mRNA還富集在多條與腫瘤生長相關(guān)的信號(hào)通路中，如PI3K-AKT、TGF-β。PI3K-AKT通路已成為人類癌癥最常激活的通路之一，靶向于PI3K和AKT的特異性抑制劑治療腫瘤已成為現(xiàn)實(shí)［15］。提示PI3K-AKT通路的異常激活參與了食管鱗癌的生長過程。TGF-β是免疫耐受性的關(guān)鍵因素，它能抑制免疫系統(tǒng)許多成分的活性和功能，參與腫瘤微環(huán)境免疫抑制，并促進(jìn)腫瘤的生長，還有文獻(xiàn)認(rèn)為其與免疫治療不良反應(yīng)有關(guān)［16］，根據(jù)本研究結(jié)果，這條通路也與食管鱗癌的發(fā)生發(fā)展有關(guān)。

綜上所述，食管鱗癌發(fā)生發(fā)展相關(guān)miRNA有hsa-miR-195-5p、hsa-miR-143-3p、hsa-miR-31-5p、hsa-miR-127-3p，靶向mRNA主要有細(xì)胞外基質(zhì)、成骨作用等生物學(xué)功能，參與的信號(hào)通路有ECM-受體相互作用、黏著斑、PI3K-Akt、TGF-β等通路。