馬幫靖，郝翔麟，韓飛，陳洪強，劉晉祎，曹佳， *，張愛華，*

（1．貴州醫(yī)科大學(xué)公共衛(wèi)生學(xué)院毒理教研室，貴州貴陽550004；2．陸軍軍醫(yī)大學(xué)軍事預(yù)防醫(yī)學(xué)系毒理學(xué)研究所，重慶400038）

肺癌是嚴(yán)重威脅人類健康和生命的惡性腫瘤，高居全球男性、發(fā)達國家女性癌癥死亡率之首[1]。近年來，我國肺癌的發(fā)病率與死亡率已超越其他惡性腫瘤躍升至第1位[2-3]。目前，肺癌發(fā)生發(fā)展的分子機制尚不十分清楚，對其深入研究有利于改善肺癌的診療及預(yù)后。

SOX30基因?qū)儆赟OX基因家族的H亞族，在雄性生殖系統(tǒng)發(fā)生發(fā)育及功能維持中起重要作用[4-5]。本課題組研究[6]發(fā)現(xiàn)， SOX30是肺癌中的新抑癌基因，通過轉(zhuǎn)錄激活p53抑制肺癌細胞增殖、誘導(dǎo)肺癌細胞凋亡，在肺腺癌中的高表達有利于患者預(yù)后[7]。最近發(fā)現(xiàn)SOX30能夠抑制肺癌轉(zhuǎn)移，但分子機制尚不清楚。因此，本研究利用基因表達譜數(shù)據(jù)功能富集分析、轉(zhuǎn)錄結(jié)合位點預(yù)測、TCGA數(shù)據(jù)庫分析篩選SOX30調(diào)控的關(guān)鍵基因；通過差異表達細胞模型和基因敲除動物模型進一步驗證相關(guān)基因的表達；采用平板克隆形成實驗及劃痕愈合實驗分析相關(guān)基因在SOX30抑癌中的主要作用，為進一步研究 SOX30基因的作用及臨床應(yīng)用提供資料。

1 材料與方法

1.1 細胞株及其培養(yǎng)

A549細胞購自中國科學(xué)院上海生命科學(xué)研究院細胞資源中心，本實驗室前期[6]已將SOX30表達載體轉(zhuǎn)入A549細胞成功構(gòu)建過表達SOX30的A549-SOX30穩(wěn)定轉(zhuǎn)染細胞系，這些細胞株均使用含10%胎牛血清的F-12K培養(yǎng)液，培養(yǎng)于37℃、CO2體積分?jǐn)?shù)為5%的常規(guī)培養(yǎng)箱中，采用0.25%胰蛋白酶消化、傳代，實驗用細胞為對數(shù)生長期細胞。

1.2 SOX30基因敲除小鼠

SOX30基因敲除小鼠品系為C57BL6，委托南京大學(xué)模式動物研究所構(gòu)建，于陸軍軍醫(yī)大學(xué)實驗動物中心(SPF級)飼養(yǎng)，由本實驗室負責(zé)繁殖維護并利用PCR對新生小鼠基因型進行鑒定。

1.3 主要試劑及儀器

F-12K培養(yǎng)基、胰酶、胎牛血清均購自Gibco公司，總RNA提取試劑Trizol Reagent以及pcDNA6.2-GW/EmGFP-miR-DSC3miRNA干擾質(zhì)粒購自Invitrogen公司， Via FectTMTransfection Reagent(E4982)轉(zhuǎn)染試劑盒、GoScriptTM逆轉(zhuǎn)錄試劑盒(A5001)及2×Taq PCRMix購自Promega公司，殺稻瘟菌素(Blasticidin)購自Sigma公司，普通RT-PCR儀購自美國Bio-Rad 公 司，化學(xué)發(fā)光凝膠成像系統(tǒng)購自英國Syngene公司。

1.4 實驗方法

1.4.1 SOX30蛋白調(diào)控基因功能富集分析基于我們前期獲得的SOX30差異表達細胞系表達譜數(shù)據(jù)[6]，使用Gene Spring Software11.0(Agilent Technologies，Santa Clara，CA，US)軟件進行數(shù)據(jù)處理，采用兩樣本t檢驗篩選差異表達基因，顯著差異表達基因 的 閾值設(shè)定為差異倍數(shù)≥2且 P <0.05，隨后于網(wǎng)站http://www.geneontology.org/在線對差異基因進行基因功能富集分析(GO enrichment analysis)，獲得關(guān)聯(lián)的主要通路及相關(guān)基因的功能。

1.4.2 DSC3啟動子區(qū)SOX30蛋白結(jié)合位點預(yù)測采用JASPAR數(shù)據(jù)庫(http://jaspar.binf.ku.dk/)在線預(yù)測DSC3(desmocollin3)啟動子區(qū)SOX30蛋白結(jié)合位點。

1.4.3 臨床標(biāo)本數(shù)據(jù)來源及熱圖繪制肺癌組織中SOX30與 DSC3基因的表達數(shù)據(jù)于UCSC Xena(https://xenabrowser.net)下載，二者數(shù)據(jù)集ID均為TCGA.LUNG.sampleMap/HiSeqV2；通過提取并進行相應(yīng)的數(shù)據(jù)匹配獲得具有 SOX30、DSC3基因表達信息的肺癌及癌旁組織共1129例，其中肺腺癌及癌旁組織576例，肺鱗癌及癌旁組織553例。肺癌組織及癌旁組織中SOX30與DSC3基因的表達熱圖在網(wǎng)站https: //genome-cancer.ucsc.edu/proj/site/hgHeatmap/在線繪制。

1.4.4 DSC3miRNA干擾細胞模型的構(gòu)建SOX30差異表達細胞模型來源于前期構(gòu)建[6]的穩(wěn)定表達pIRES2-EGFP-SOX30質(zhì)粒(A549-SOX30)和對照質(zhì)粒(A549-Vector)。按每孔3.0×105個細胞接種A549-SOX30于6孔板，37℃、CO2體積分?jǐn)?shù)為5%條件下培養(yǎng)過夜；參照Via FectTMTransfection Reagent(E4982)轉(zhuǎn)染試劑盒說明書配制pcDNA6.2-GW/EmGFP-miR-DSC3miRNA轉(zhuǎn)染復(fù)合物，按每孔200μL分別加至6孔板中，搖勻后于37℃、CO2體積分?jǐn)?shù)為5%的常規(guī)培養(yǎng)箱中繼續(xù)培養(yǎng)，48 h后換用含3μg/mL殺稻瘟菌素的F-12K完全培養(yǎng)液篩選細胞，直至出現(xiàn)熒光陽性克隆；挑取熒光陽性克隆，擴大培養(yǎng)，獲得A549-SOX30+DSC3miRNA穩(wěn)定表達細胞系并用于后續(xù)分析。

1.4.5 細胞與組織的總RNA提取及cDNA合成細胞株總RNA的抽提參照Trizol(Invitrogen)試劑使用說明書進行，并采用Nano Drop分光光度計測定所提總RNA的濃度和純度，然后參照GoScriptTM逆轉(zhuǎn)錄試劑盒(A5001)使用說明書合成cDNA。 SOX30基因敲除小鼠肺組織首先進行液氮研磨，之后參照上述方法提取總RNA并合成cDNA。

1.4.6 逆轉(zhuǎn)錄PCRSOX30引物，上游序列5′-GATG TCCCGCTCACCGTGTTGC-3′，下游序列5′-GACAGGG CTTGGGCTCTGGACT-3′；DSC3引物，上游序列5′-GACCCTCGTGATCTTCAGTCG-3′，下游序列5′-TCAC TTGACCGGATGAGGTCT-3′；內(nèi)參基因 β -actin引物，上游序列5′-GAGCTACGAGCTGCCTGACGG-3′，下游序列5′-CCTAGAAGCATTTGCGGTGG-3′。逆轉(zhuǎn)錄PCR反應(yīng)體系(20μL)包括：上、下游引物各1μL，2×Taq PCRMix10μL，cDNA模板1μL，Nuclease-free水7 μL。擴增程序為：95℃、3min；94℃、30s；58℃、30s；72℃、30s；72℃、10min；33個循環(huán)。擴增結(jié)束后，各取10μL PCR產(chǎn)物進行瓊脂糖凝膠電泳，電泳結(jié)束后使用凝膠成像系統(tǒng)進行曝光分析，拍照并保存圖片。實驗重復(fù)3次。

1.4.7 平板克隆形成實驗按每孔200個細胞接種A549-Vector、A549-SOX30、A549-SOX30+DSC3miRNA穩(wěn)定轉(zhuǎn)染細胞于6孔板，搖勻后于37℃、CO2體積分?jǐn)?shù)為5%的常規(guī)培養(yǎng)箱中，用含3μg/mL殺稻瘟菌素的F-12K完全培養(yǎng)液培養(yǎng)2周，直至出現(xiàn)肉眼可見的克隆時，終止培養(yǎng)，取出細胞，PBS清洗細胞2次，4%多聚甲醛室溫固定30min；結(jié)晶紫染色5min；PBS清洗細胞2次，空氣中干燥；拍照保存，分析結(jié)果。

1.4.8 細胞劃痕實驗分別將A549-Vector、A549-SOX30、A549-SOX30+DSC3miRNA穩(wěn)定轉(zhuǎn)染細胞系按每孔5×105個細胞接種于6孔板中，搖勻后于37℃、CO2體積分?jǐn)?shù)為5%的常規(guī)培養(yǎng)箱中培養(yǎng)；待細胞鋪滿皿底(融合度>90%)，用10μL槍頭在6孔板中劃痕，PBS清洗細胞2次；然后加入無血清的F-12K培養(yǎng)液繼續(xù)培養(yǎng)，于劃痕后0、48h在同一位置利用倒置顯微鏡進行拍照；統(tǒng)計分析細胞遷移情況。

1.5 統(tǒng)計分析

采用SPSS17.0軟件進行統(tǒng)計分析，計量資料均采用x±s 表 示，兩樣本均數(shù)比較采用成組t檢驗，肺癌中SOX30與DSC3基因表達的相關(guān)性分析采用Spearman秩相關(guān)分析。檢驗水準(zhǔn) α=0.05。

2 結(jié)果

2.1 SOX30轉(zhuǎn)錄調(diào)控基因全基因組表達譜數(shù)據(jù)功能富集分析

為了尋找 SOX30基因抑制腫瘤的可能下游通路，我們首先對 SOX30基因差異表達細胞系基因芯片數(shù)據(jù)進行了功能富集分析(GO分析)，結(jié)果發(fā)現(xiàn)SOX30廣泛與黏著斑(focal a dhesion)、錨定連接(anchoring junction)、黏著連接(adherens junction)、細胞連接(cell juncion)、細胞骨架(cytoskeleton)等在腫瘤侵襲和轉(zhuǎn)移中具有關(guān)鍵調(diào)控作用的細胞連接相關(guān)基因表達顯著相關(guān)，提示SOX30抑癌功能發(fā)揮可能與細胞黏附相關(guān)通路有關(guān)(見圖1)。

圖1 基于SOX30轉(zhuǎn)錄調(diào)控基因表達譜經(jīng)功能富集分析潛在的下游通路

2.2 生物信息學(xué)預(yù)測發(fā)現(xiàn) SOX30基因可能調(diào)控DSC3基因的表達

為了進一步明確哪些細胞間連接相關(guān)基因可能受SOX30蛋白的調(diào)控，我們利用JASPAR數(shù)據(jù)庫(http://jaspar.binf.ku.dk)對SOX30蛋白潛在的轉(zhuǎn)錄調(diào)控基因進行了在線預(yù)測，結(jié)果發(fā)現(xiàn)在 DSC3基因啟動子區(qū)存在多個SOX30蛋白結(jié)合位點(見表1)，提示SOX30很可能調(diào)控DSC3的表達。

表1 DSC3啟動子區(qū) SOX30蛋白結(jié)合位點預(yù)測

2.3 TCGA數(shù)據(jù)庫分析 SOX30與 DSC3基因在肺腺癌中表達的相關(guān)性

DSC3是一種鈣依賴性糖蛋白，與腫瘤細胞增殖、集落形成、侵襲和遷移等調(diào)控有關(guān)。由此我們推測SOX30基因的抑癌作用可能與DSC3有關(guān)。為了進一步明確SOX30與DSC3在肺癌中的潛在調(diào)控關(guān)系，我們利用TCGA數(shù)據(jù)庫分析了 SOX30基因與 DSC3基因在非小細胞肺癌中表達 的 相關(guān)性，結(jié)果發(fā)現(xiàn)二者在非小細胞肺癌中表達無明顯相關(guān)(見圖2A)。由于前期研究[7]表明SOX30主要在肺腺癌中發(fā)揮作用，為此我們進一步分析了二者在肺腺癌和肺鱗癌中的表達相關(guān)性，結(jié)果顯示二者在肺腺癌中表達正相關(guān)(見圖2B)，在肺鱗癌中表達無明顯相關(guān)(見圖2C)，提示二者在臨床肺腺癌中存在調(diào)控關(guān)系。

圖2 SOX30 與 DSC3基因在 TCGA數(shù)據(jù)庫肺癌中的表達及相關(guān)性分析

2.4 SOX30蛋白是 DSC3基因轉(zhuǎn)錄調(diào)控的關(guān)鍵分子

為進一步確證 SOX30 與 DSC3基因之間的調(diào)控關(guān)系，我們分析了SOX30差異表達細胞系中DSC3的表達， SOX30基因高表達的A549-SOX30細胞中DSC3 mRNA表達也顯著上調(diào)(見圖3)，說明 SOX30基因能夠調(diào)控 DSC3基因的表達，但是這種調(diào)控是不是必須的呢？為此，我們利用 SOX30基因敲除動物模型，使用逆轉(zhuǎn)錄PCR檢測了野生型(SOX30+/+) 、雜合子(SOX30+/-)以及純合子(SOX30-/-)小鼠肺組織中DSC3mRNA 的 表達情況，結(jié)果顯示純合子(SOX30-/-)小鼠肺組織中DSC3 mRNA表達水平顯著低于雜合子(SOX30+/-)及野生型(SOX30+/+)小鼠(見圖4)，表明SOX30蛋白是調(diào)控 DSC3基因表達的關(guān)鍵分子。

圖3 各細胞中SOX30與 DSC3mRNA 的表達

2.5 SOX30蛋白對A549細胞的集落形成抑制依賴于DSC3基因的表達

圖4 SOX30基因敲除小鼠肺組織中SOX30與DSC3m RNA 的表達

為確定 DSC3基因是否參與SOX30蛋白抑制肺癌細胞增殖的過程，我們進行了克隆形成實驗。結(jié)果顯示，過表達SOX30的A549-SOX30細胞集落形成數(shù)較對照A549-Vector細胞顯著降低，表明SOX30高表達抑制腫瘤細胞集落形成，而在A549-SOX30+DSC3miRNA細胞中， DSC3基因的表達被干擾后細胞集落形成數(shù)較A549-SOX30細胞顯著增加(見圖5)。這提示 SOX30基因抑制腫瘤細胞增殖的作用依賴于 DSC3基因的表達。

圖5 干擾DSC3基因表達減弱SOX30基因?qū)?A 549細胞集落形成的抑制

2.6 SOX30基因?qū)549細胞的遷移抑制依賴于DSC3基因的表達

由于DSC3是橋粒關(guān)鍵構(gòu)成蛋白，為探索 DSC3基因是否參與了 SOX30基因?qū)δ[瘤轉(zhuǎn)移的抑制，我們利用劃痕愈合實驗分析表明，過表達SOX30的A549-SOX30細胞相對遷移率較對照A549-Vector細胞顯著降低，表明SOX30高表達抑制腫瘤細胞遷移；而在A549-SOX30+DSC3miRNA細胞中，即 DSC3基因的表達被干擾后顯著減弱 SOX30基因?qū)549細胞遷移的抑制作用，細胞遷移率較A549-SOX30細胞顯著升高(見圖6)。這提示SOX30基因抑制腫瘤轉(zhuǎn)移與DSC3高度相關(guān)。

圖6 干擾DSC3基因表達減弱 SOX30基因?qū)?A 549細胞遷移的抑制

3 討論

脊椎動物細胞間連接主要有黏著連接(adherens junction，AJ)和橋粒(半橋粒)兩種類型。橋粒(desmosome)在細胞間信息交流、細胞定位、細胞形態(tài)、細胞結(jié)構(gòu)和組織機械張力等調(diào)控中起著重要作用[8-9]。腫瘤細胞局部浸潤必要步驟就是原發(fā)腫瘤細胞失去細胞間連接，包括黏著連接和橋粒結(jié)構(gòu)的分解和調(diào)控[10]。多項研究表明，橋粒相關(guān)蛋白參與腫瘤細胞的增殖、侵襲及轉(zhuǎn)移等過程，在腫瘤的發(fā)生發(fā)展中起重要作用[11]。近年研究發(fā)現(xiàn)，轉(zhuǎn)錄因子對橋?；虻霓D(zhuǎn)錄激活是橋粒基因表達調(diào)控的一種重要途徑[12-16]。SOX30是典型的轉(zhuǎn)錄因子，作為一個新的抑癌基因，在肺癌的發(fā)生中起著重要作用，但是相關(guān)的機制還不清楚。本研究對SOX30差異表達基因進行了功能富集分析，發(fā)現(xiàn)SOX30基因的功能主要與黏著斑、錨定連接、黏著連接等細胞連接信號通路顯著相關(guān)，同時通過轉(zhuǎn)錄結(jié)合位點分析發(fā)現(xiàn)橋粒關(guān)鍵基因 DSC3可能受SOX30蛋白的轉(zhuǎn)錄調(diào)控。 DSC3基因定位于人類染色體18q12.1，編碼產(chǎn)物含896個氨基酸多肽。DSC3蛋白是一種鈣依賴性糖蛋白，且是鈣黏蛋白超家族的橋粒子家族的成員，可與其他家族成員相互作用而促成細胞橋粒介導(dǎo)的細胞黏附[17]。多項研究表明，DSC3在多種腫瘤中表達下調(diào)，包括結(jié)直腸癌[18-20]、前列腺癌[21]和乳腺癌[22]。Cui等[14]研究發(fā)現(xiàn)DSC3是肺癌中一個新的抑癌基因，它可以抑制肺癌的增殖、侵襲和遷移。我們推測SOX30基因可能通過轉(zhuǎn)錄調(diào)控DSC3而發(fā)揮抑癌功能。為此，我們進一步利用差異表達細胞模型和基因敲除動物模型驗證，結(jié)果發(fā)現(xiàn)，過表達SOX30后顯著上調(diào)A549細胞中DSC3mRNA 的 表達，敲除 SOX30基因后，小鼠肺組織中DSC3mRNA 的 表達顯著下調(diào)，表明SOX30對DSC3 的表達有著關(guān)鍵的調(diào)控作用，提示SOX30可能是一個新的橋?；蜣D(zhuǎn)錄調(diào)控因子。

由于DSC3受SOX30表達調(diào)控，因此有可能在SOX30抑癌中發(fā)揮重要作用。我們前期的研究[7]也表明，SOX30在肺腺癌中是一個典型的抑癌基因，其表達與肺腺癌患者預(yù)后正相關(guān)，而在肺鱗癌中無明顯的抑癌作用，因此，SOX30與DSC3的表達僅在肺腺癌中相關(guān)，與 SOX30基因的功能一致。為明確DSC3是否是SOX30基因在肺腺癌中發(fā)揮抑癌功能的可能下游分子，我們利用集落形成實驗和細胞劃痕愈合實驗檢測干擾 DSC3基因的表達后是否會影響SOX30對A549細胞的增殖和遷移抑制，結(jié)果顯示干擾 DSC3基因的表達顯著減弱SOX30對A549細胞的增殖和遷移抑制，表明DSC3是SOX30在肺腺癌中發(fā)揮抑癌功能的重要分子。但是，二者之間的這種調(diào)控作用，如何在腫瘤發(fā)生過程中影響橋粒的形態(tài)結(jié)構(gòu)？是否二者進一步調(diào)控其他下游信號通路，都值得進一步研究。

綜上所述，本研究首次發(fā)現(xiàn)SOX30是橋粒基因DSC3的轉(zhuǎn)錄調(diào)控因子，且DSC3是SOX30發(fā)揮抑癌功能的重要下游靶分子。本研究為進一步探索 SOX30基因功能以及SOX30-DSC3在肺腺癌發(fā)生發(fā)展中的作用和臨床意義提供了重要依據(jù)。

[1]TORRE L A，BRAY F，SIEGEL R L，et al.Global cancer statistic，2012[J].CA CancerJClin，2015，65(2)：87-108.

[2]CHEN W，ZHENG R，ZUO T，et al.National cancer incidence and mortality in China，2012[J].ChinJCancer Res，2016，28(1)：1-11.

[3]CHEN W，ZHENG R，ZHANG S，et al.Cancer incidence and mortality in China in2013：an analysis based on urbanization level[J].ChinJCancer Res，2017，29(1)：1-10.

[4]HAN F，DONGY，LIU W，et al.Epigenetic regulation of Sox30is associated with testis development in mice[J].PLoSOne，2014，9(5)：e97203.

[5]FENG C A，SPILLER C，MERRINER DJ，et al.SOX30is required for male fertility in mice[J].Sci Rep，2017，7(1)：17619.

[6]HAN F，LIU W，JIANG X，et al.SOX30，a novel epigenetic silenced tumor suppressor，promotes tumor cell apoptosis by transcriptional activating p53in lung cancer[J].Oncogene，2015，34(33)：4391-4402.

[7]HAN F，LIU W，XIAO H，et al.High expression of SOX30is associated with favorable survival in human lung adenocarcinoma[J].Sci Rep，2015，5：13630.

[8] GARROD D， CHIDGEYM.Desmosome structure，composition and function[J].Biochim Biophys Acta，2008，1778(3)：572-587.

[9]HUBER O，PETERSEN I.150thanniversary series：desmosomes and the hallmarks of cancer[J].Cell Commun Adhes，2015，22(1)：15-28.

[10]DUSEK R L，ATTARDI L D.Desmosomes：new perpetrators in tumour suppression[J].Nat Rev Cancer，2011，11(5)：317-323.

[11]ZHOU G，YANG L，GRAY A，et al.Therole of desmosomes in carcinogenesis[J].Onco Targets Ther，2017，10：4059-4063.

[12]VANDEWALLE C，COMIJNJ，DE CRAENE B，et al.SIP1/ZEB2induces EMTby repressing genes of different epithelial cell-cell junctions[J].Nucleic Acids Res，2005，33(20)：6566-6578.

[13]BAZZI H，DEMEHRI S，POTTER C.Desmoglein4is regulated by transcription factors implicated in hair shaft differentiation[J].Differentiation，2009，78(5)：292-300.

[14]CUI T，CHENY，YANG L，et al.The p53target gene desmocollin3acts as a novel tumor suppressor through inhibiting EGFR/ERK pathway in human lung cancer[J].Carcinogenesis，2012，33(12)：2326-2333.

[15]MAOX，CHOMJT，ELLEBRECHTC T，et al.Stat3 regulates desmoglein3transcription in epithelial keratinocytes[J].JCI Insight，2017，2(9)：1- 1 0.

[16]LIUY，CHIDGEYM，YANG V W，et al.Krüppel-like factor 5is essential for maintenance of barrier function in mouse colon[J].AmJPhysiol Gastrointest Liver Physiol， 2017，313(5)：478-491.

[17]EZZATE S，TAHOUN N.Therole of Napsin-A and Desmocollin-3in classifying poorly differentiating non-small cell lung carcinoma[J].J Egypt Natl Canc Inst，2015，28(1)：13-22.

[18]CUI T，CHENY，YANG L，et al.DSC3expression is regulated by p53，and methylation of DSC3DNA is a prognostic marker in human colorectal cancer[J].BrJCancer，2011，104(6)：1013-1019.

[19]BUJKOM，KOBER P，MIKULAM，et al.Expression changes of cell-cell adhesion-related genes in colorectal tumors[J].Oncol Lett，2015，9(6)：1-14.

[20]KN?SEL T， CHENY， HOTOVY S， et al.Loss of desmocollin1-3and homeobox genes PITX1and CDX2are associated with tumor progression anDSurvival in colorectal carcinoma[J].IntJColorectal Dis，2012，27(11)：1391-1399.

[21]PANJ，CHENY，MO C，et al.Association of DSC3mRNA down-regulation in prostate cancer with promoter hypermethylation and poor prognosis[J].PLoSOne，2014，9(3)：e92815.

[22]OSHIROMM，KIM CJ，WOZNIAK RJ，et al.Epigenetic silencing of DSC3is a common event in human breast cancer[J].Breast Cancer Res，2005，7(5)：669-680.