復(fù)旦大學(xué)附屬腫瘤醫(yī)院麻醉科，復(fù)旦大學(xué)上海醫(yī)學(xué)院腫瘤學(xué)系，上海 200032

肺癌是全球人群中發(fā)病率最高的腫瘤（占所有腫瘤的11.6%），其死亡率也位居各類(lèi)腫瘤之首（占18.4%）［1］。目前，手術(shù)治療仍是中早期肺癌治療的首選方法，但仍有許多患者喪失手術(shù)機(jī)會(huì)或在術(shù)后出現(xiàn)復(fù)發(fā)和轉(zhuǎn)移，因而早發(fā)現(xiàn)、早診斷、早治療在肺癌的診治中至關(guān)重要。對(duì)于部分無(wú)法手術(shù)或術(shù)后出現(xiàn)復(fù)發(fā)和轉(zhuǎn)移的患者，分子靶向治療作為新的治療手段取得了良好的效果，正越來(lái)越受到重視。與傳統(tǒng)化療相比，靶向藥物具有特異性強(qiáng)、療效明顯、不良反應(yīng)較少等優(yōu)點(diǎn)［2］。然而，靶向藥物僅適用于目標(biāo)基因有突變的患者，且治療一段時(shí)間后，大多數(shù)患者對(duì)靶向藥物產(chǎn)生耐受。這表明，肺癌復(fù)發(fā)轉(zhuǎn)移的機(jī)制異常復(fù)雜，有必要進(jìn)一步探索肺癌細(xì)胞的遷移和侵襲機(jī)制，以尋找更好的生物標(biāo)志物，為腫瘤治療、靶向藥物的研發(fā)提供理論支持。

近年來(lái)，越來(lái)越多的研究發(fā)現(xiàn)，上皮-間質(zhì)轉(zhuǎn)化（epithelial-mesenchymal transition，EMT）介導(dǎo)了腫瘤的轉(zhuǎn)移、復(fù)發(fā)和腫瘤耐藥［3］。在EMT過(guò)程中，表皮表型的腫瘤細(xì)胞失去細(xì)胞極性和細(xì)胞間的緊密連接，轉(zhuǎn)化為沒(méi)有細(xì)胞極性且遷移侵襲能力較強(qiáng)的間質(zhì)表型。此外，EMT還受到Twist家族BHLH轉(zhuǎn)錄因子1（Twist family BHLH transcription factor 1，Twist1）、Twist2、Snail超家族鋅指轉(zhuǎn)錄因子（Snail family zinc finger，SNAI）1、SNAI2、鋅指E盒同源結(jié)合蛋白1（zinc finger E-box binding homeobox 1，ZEB1）、ZEB2等轉(zhuǎn)錄因子的調(diào)控，它們常被視作EMT的生物標(biāo)志物，可統(tǒng)稱(chēng)為EMT誘導(dǎo)轉(zhuǎn)錄因子（EMT-inducing transcription factor，EMT-TF）［3］。研究發(fā)現(xiàn)，EMT還可激活腫瘤干細(xì)胞，導(dǎo)致耐藥［4］。因此研究EMT-TF在肺癌中的作用有助于肺癌的病理診斷與預(yù)后分析。

Twist1和Twist2是堿性的螺旋-環(huán)-螺旋蛋白（basic helix-loop-helix protein，bHLH）家族中的兩個(gè)重要轉(zhuǎn)錄因子。在多種腫瘤中，Twist1呈現(xiàn)高度甲基化，且表達(dá)明顯增加［5］，它可通過(guò)調(diào)控EMT促進(jìn)腫瘤的增殖和轉(zhuǎn)移［6］。此外，Twist1還可通過(guò)誘導(dǎo)EMT、激活腫瘤干細(xì)胞，致使多種腫瘤對(duì)靶向藥物發(fā)生耐受［7-9］。然而，與Twist1相反，Twist2對(duì)腫瘤有明顯的抑制作用，它在多種腫瘤組織中，無(wú)論是mRNA水平還是蛋白水平都有明顯的下降［10］。ZEB家族在激活腫瘤干細(xì)胞、調(diào)控腫瘤細(xì)胞凋亡、促進(jìn)腫瘤血管生成和化療藥耐藥等方面發(fā)揮著重要作用，因而有可能成為腫瘤診斷、預(yù)后分析以及治療的靶點(diǎn)［11］。ZEB家族有兩個(gè)重要轉(zhuǎn)錄因子ZEB1和ZEB2，它們主要在細(xì)胞核內(nèi)發(fā)揮作用。SNAI是保守的鋅指家族的轉(zhuǎn)錄因子之一，它對(duì)果蠅中胚層的形成至關(guān)重要［12］。SNAI在轉(zhuǎn)錄水平受到很多生長(zhǎng)因子和信號(hào)分子的調(diào)節(jié)，這些因子通過(guò)結(jié)合SNAI基因的啟動(dòng)子區(qū)域發(fā)揮調(diào)控作用［13-14］。SNAI1能夠與組蛋白去乙?；窰DAC1/2相互作用，通過(guò)改變?nèi)旧|(zhì)的局部結(jié)構(gòu)，抑制E-cadherin的表達(dá)，進(jìn)而調(diào)節(jié)EMT［15］。隨著研究的深入，人們?cè)絹?lái)越關(guān)注EMT-TF在腫瘤轉(zhuǎn)移復(fù)發(fā)中的作用，但目前有關(guān)EMT-TF與腫瘤分期、預(yù)后預(yù)測(cè)方面的研究很少，因此我們期望利用腫瘤在線數(shù)據(jù)庫(kù)，篩選對(duì)肺癌診斷、分期和預(yù)后預(yù)測(cè)方面有一定價(jià)值的轉(zhuǎn)錄因子。

1 材料和方法

1.1 Oncomine Database分析、GEPIA和Lung Cancer Explore

Oncomine Database數(shù)據(jù)庫(kù)中共有肺癌數(shù)據(jù)集74個(gè)，包括6 289個(gè)樣本，我們通過(guò)Oncomine Database分析EMT TF在不同腫瘤中的表達(dá)水平，進(jìn)而分析EMT TF在肺癌組織和正常組織中的表達(dá)差異（cut-off值：P=0.01，fold change of 2）。GEPIA公共在線數(shù)據(jù)庫(kù)包含TCGA數(shù)據(jù)庫(kù)中9 736例腫瘤標(biāo)本以及GTEx數(shù)據(jù)庫(kù)8 587例正常組織樣本［16］。我們應(yīng)用GEPIA比較了EMT-TF基因在肺癌（肺腺癌和鱗狀細(xì)胞癌）和正常組織中的表達(dá)差異，同時(shí)還將分析EMT-TF mRNA表達(dá)水平與腫瘤臨床病理學(xué)分期的相關(guān)性。

1.2 Kaplan-Meier Plotter

Kaplan-Meier Plotter能夠檢測(cè)近54 000個(gè)基因?qū)?1種腫瘤生存預(yù)后的影響［17］。我們應(yīng)用Kaplan-Meier Plotter Lung Cancer數(shù)據(jù)庫(kù)分析EMT TF mRNA水平對(duì)肺癌患者5年生存率的影響，例如總生存期（overall survival，OS）、無(wú)進(jìn)展生存期（progression-free survival，PFS）、進(jìn)展后生存期（post-progression survival，PPS）。

1.3 STRING和Mentha數(shù)據(jù)庫(kù)

STRING數(shù)據(jù)庫(kù)收集了目前研究的蛋白質(zhì)間的相互作用，我們通過(guò)多基因分析篩選了EMTTF的蛋白質(zhì)-蛋白質(zhì)相互作用（protein-protein interaction，PPI）網(wǎng)絡(luò)。Mentha數(shù)據(jù)庫(kù)的數(shù)據(jù)來(lái)自于依附IMEx聯(lián)盟的人工管理的PPI數(shù)據(jù)庫(kù)，聚集的數(shù)據(jù)為科研工作提供了一系列工具，在PPI網(wǎng)絡(luò)的背景下分析選定的蛋白質(zhì)。通過(guò)Mentha數(shù)據(jù)庫(kù)，分析EMT-TF相互作用的系數(shù)以及與之作用的分子。

2 結(jié)果

2.1 EMT-TF在肺癌組織中的表達(dá)水平

應(yīng)用ONcomine數(shù)據(jù)庫(kù)，我們對(duì)比了多種腫瘤組織和正常組織中EMT-TF的轉(zhuǎn)錄水平（圖1），并將肺癌有關(guān)的結(jié)果單獨(dú)歸納于表1。在Hou肺癌數(shù)據(jù)庫(kù)中，鱗癌、腺癌組織中Twist1 mRNA水平分別是正常組織的7.7倍和2.54倍，其差異有統(tǒng)計(jì)學(xué)意義（P＜0.01）［18］。在包括Hou等［18］報(bào)道的數(shù)據(jù)庫(kù)在內(nèi)的4個(gè)肺癌數(shù)據(jù)庫(kù)中，鱗癌SNAI2轉(zhuǎn)錄水平均顯著高于正常組織：Yamagata數(shù)據(jù)庫(kù)中，SNAI2轉(zhuǎn)錄水平增加2倍［19］；Wachi等［20］和Hou等［18］報(bào)道，SNAI2轉(zhuǎn)錄水平增加3倍多；Garber數(shù)據(jù)庫(kù)中［21］，SNAI2轉(zhuǎn)錄水平增加4倍多，差異均有統(tǒng)計(jì)學(xué)意義（P＜0.05）。

圖1 EMT-TF在不同腫瘤中的轉(zhuǎn)錄水平Fig.1 Transcription level of EMT-TF in different types of cancer

表1 EMT-TF轉(zhuǎn)錄水平在肺癌和肺組織中的改變Tab.1 The significant changes in transcription level of EMT-TF between different types of lung cancer and normal tissues

2.2 EMT-TF的mRNA表達(dá)與肺癌臨床病理學(xué)分期的關(guān)系

通過(guò)搜索基因表達(dá)譜數(shù)據(jù)動(dòng)態(tài)分析（Gene Expression Profiling Interactive Analysis，GEPIA）數(shù)據(jù)庫(kù)比較了EMT-TF在肺癌組織和正常組織中的表達(dá)差異。結(jié)果顯示，與正常組織相比，在肺鱗癌中Twist1和SNAI2 mRNA表達(dá)增加，Twist2、ZEB1、ZEB2和SNAI1表達(dá)降低；而在腺癌中，Twist2、ZEB1、ZEB2和SNAI1 mRNA表達(dá)水平降低，Twist1、SNAI2 mRNA表達(dá)無(wú)差異（圖2A）。

為探究EMT-TF的表達(dá)水平與腫瘤分期的相關(guān)性，我們采用stage plot方法對(duì)GEPIA數(shù)據(jù)進(jìn)行了分析，結(jié)果顯示，Twist1、SNAI1和SNAI2表達(dá)水平與肺癌臨床病理分期具有明顯相關(guān)性，而Twist2、ZEB1和ZEB2的表達(dá)水平與臨床病理學(xué)分期無(wú)相關(guān)性（圖2B）。

2.3 Twist1與肺癌患者預(yù)后相關(guān)

我們采用Kaplan-Meier Plotter工具分析了EMT-TF mRNA水平與肺癌5年生存率之間的關(guān)系，結(jié)果發(fā)現(xiàn)，Twist1和Twist2的表達(dá)增加提示肺癌患者的不良預(yù)后（P＜0.05）：例如Twist1、Twist2高表達(dá)組患者OS、PFS均較低（P＜0.05），SNAI1高表達(dá)組患者PFS、OS和PPS均較低（P＜0.05），ZEB1、ZEB2低表達(dá)與肺癌的不良預(yù)后相關(guān)，其中ZEB1的表達(dá)降低與肺癌OS降低相關(guān)，ZEB2的低表達(dá)與肺癌的PFS和OS預(yù)后不良顯著相關(guān)（圖3）。

2.4 EMT-TF的相互作用及通路

我們首先采用STRING分析了EMT-TF之間的PPI，結(jié)果顯示，EMT-TF之間存在直接或間接的相互作用，其中Twist1和Twist2在信號(hào)網(wǎng)絡(luò)中起關(guān)鍵作用（圖4A）。隨后，我們采用Mentha數(shù)據(jù)庫(kù)進(jìn)一步展示了TF之間的相互作用以及與它們相互作用的分子，并且顯示了分子作用的部位（圖4B）。KEGG通路分析發(fā)現(xiàn)，EMT-TF受到細(xì)胞外界微環(huán)境的變化，如細(xì)胞因子等的調(diào)控，發(fā)揮促進(jìn)腫瘤遷移和侵襲的作用，進(jìn)而導(dǎo)致轉(zhuǎn)移的發(fā)生。

圖2 EMT-TF的mRNA水平以及與肺癌臨床病理分期的關(guān)系Fig.2 The association of EMT-TF mRNA transcription level and tumor stage in lung cancer patients

圖3 EMT-TF與肺癌預(yù)后的相關(guān)性Fig.3 The correlation of EMT-TF expression level and prognosis of lung cancer patients

圖4 PPI和EMT-TF相互作用的分子Fig.4 PPI network of EMT-TF and the cofunction molecules

3 討 論

EMT是肺癌侵襲和轉(zhuǎn)移的重要環(huán)節(jié)，非小細(xì)胞肺癌患者可通過(guò)EMT機(jī)制對(duì)腫瘤治療藥物產(chǎn)生耐藥。研究發(fā)現(xiàn)，TGF-β可通過(guò)下列途徑介導(dǎo)EMT：①通過(guò)Smad途徑介導(dǎo)EMT；② 通過(guò)誘導(dǎo)Smad復(fù)合體浸入細(xì)胞核進(jìn)而調(diào)節(jié)EMT相關(guān)基因和因子的表達(dá)；③可增加SNAI1、SNAI2和ZEB1/2等EMT-TF的表達(dá)［22］。LncRNA也能通過(guò)影響EMT的進(jìn)程和腫瘤干細(xì)胞特性，發(fā)揮調(diào)節(jié)非小細(xì)胞肺癌轉(zhuǎn)移的作用［23-24］。研究表明，miRNA也能調(diào)控EMT，如抑癌因子miR-451a的過(guò)表達(dá)可破壞肺癌細(xì)胞的EMT，顯著增強(qiáng)肺癌細(xì)胞對(duì)多柔比星治療的敏感性［25］。在肺癌中，很多促癌基因也能促進(jìn)EMT，比如FAM83D［26］、FOXK1［27］、KLK4［28］等。因此識(shí)別并阻斷EMT進(jìn)程，對(duì)降低腫瘤的侵襲和轉(zhuǎn)移非常重要，同時(shí)也可降低進(jìn)展期肺癌患者的負(fù)荷，進(jìn)而提高對(duì)手術(shù)的耐受和預(yù)后，這對(duì)肺癌乃至所有腫瘤患者都很重要。

研究發(fā)現(xiàn)Twist1在肺癌中受很多因素的調(diào)節(jié)，如LncRNA JPX上調(diào)Twist1后通過(guò)Wnt/β-catenin信號(hào)通路誘導(dǎo)EMT和肺癌細(xì)胞侵襲［29］；HR23A與Twist1相互作用促進(jìn)了Twist1的泛素介導(dǎo)的蛋白酶體降解，從而抑制了EMT、癌細(xì)胞遷移和干細(xì)胞特性［30］；此外還包括表觀遺傳相關(guān)分子和miRNA。目前對(duì)Twist2的研究寥寥無(wú)幾，并且對(duì)其在肺癌中的作用并無(wú)明確一致的觀點(diǎn)［31-32］。本研究發(fā)現(xiàn)Twist1在肺癌中的表達(dá)增加，并且與肺癌的臨床病理學(xué)分期相關(guān)，Twist1表達(dá)增加提示肺癌患者的不良預(yù)后，以上皆提示Twist1在肺癌診斷和研究中的重要性和作為肺癌潛在標(biāo)志物的可能。

SNAI1是一種重要的EMT-TF，可以促進(jìn)腫瘤的轉(zhuǎn)移，從機(jī)制上講，SNAI1向其目標(biāo)啟動(dòng)子招募多個(gè)參與組蛋白去乙?；?、甲基化和泛素化的抑制蛋白復(fù)合物，并發(fā)揮其抑制功能［33］。臨床上發(fā)現(xiàn)，SNAI1的表達(dá)和化療耐藥、生存率降低、復(fù)發(fā)率較高和預(yù)后較差相關(guān)［15，33-34］。然而SNAI1是不穩(wěn)定的，通過(guò)泛素化介導(dǎo)的蛋白酶體途徑很快被降解，通過(guò)穩(wěn)定SNAI1能夠促進(jìn)肺癌細(xì)胞的遷移［33］。腫瘤在線數(shù)據(jù)庫(kù)分析發(fā)現(xiàn)SNAI1在肺癌組織中的表達(dá)低于正常組織，并且與分期密切相關(guān)，但其預(yù)后和在腫瘤中的表達(dá)存在分歧。目前對(duì)于SNAI1的研究甚少，SNAI1在肺癌中降低的機(jī)制，以及它調(diào)控的基因仍有待進(jìn)一步探討；同時(shí)我們分析不同的研究存在樣本的差異，以及臨床數(shù)據(jù)收集樣本的特異性等，使得無(wú)法對(duì)其原因進(jìn)行深入分析。

ZEB1/2的作用相似，在肺癌中，各種調(diào)節(jié)因子通過(guò)作用于ZEB1/2的啟動(dòng)子，調(diào)節(jié)它們的表達(dá)，進(jìn)而調(diào)節(jié)EMT和肺癌的侵襲和轉(zhuǎn)移［35-37］。ZEB1/2通過(guò)數(shù)據(jù)庫(kù)分析發(fā)現(xiàn)其表達(dá)水平在肺癌中降低，但和臨床病理學(xué)分期無(wú)相關(guān)性，ZEB1/2表達(dá)降低與患者的不良預(yù)后相關(guān)。同時(shí)EMT-TF在肺癌中的作用逐漸明了，本研究也發(fā)現(xiàn)了其與肺癌分期和預(yù)后的相關(guān)性，這可能輔助于臨床診斷、分期和預(yù)后預(yù)測(cè)，同時(shí)未來(lái)的研究需要進(jìn)一步探索外界及腫瘤微環(huán)境中調(diào)節(jié)EMT-TF的因子，為腫瘤的預(yù)防和治療提供理論依據(jù)。

Twist1、SNAI2、ZEB1和ZEB2在肺癌中的表達(dá)量和正常組織不同，Twist1和SNAI2的表達(dá)量與腫瘤的分期與預(yù)后相關(guān)，ZEB1和ZEB2在肺癌中表達(dá)量降低與肺癌患者的預(yù)后不良相關(guān)。因此，Twist1和SNAI2可能是肺癌診斷和研究的潛在的生物標(biāo)志物，ZEB1和ZEB2可作為肺癌預(yù)后的預(yù)測(cè)分子。