鄭艷 張茜 孫菲

甲狀腺癌是一種好發(fā)于女性的內(nèi)分泌系統(tǒng)常見惡性腫瘤，近年來甲狀腺癌的發(fā)病率逐年增加，據(jù)報道，目前我國新發(fā)甲狀腺癌患者約占全身惡性腫瘤患者的3%，約占頭頸部惡性腫瘤患者的5%[1，2]。根據(jù)組織病理學特征，甲狀腺癌可分為乳頭狀癌、濾泡狀癌、髓樣癌及未分化癌等病理類型，其中，以甲狀腺乳頭狀癌最為常見，約占甲狀腺癌總數(shù)的85%，而甲狀腺未分化癌最為少見，但惡性程度高，預(yù)后差[3，4]。目前，臨床上對于甲狀腺癌患者的主要治療手段包括手術(shù)切除、I131放射治療及內(nèi)分泌治療等[5]。早期患者能取得較好的治療效果，但容易出現(xiàn)腫瘤復(fù)發(fā)；而對于晚期患者，治療效果較差，5年生存率不容樂觀[6]。因此，探究甲狀腺癌的發(fā)病機制，尋找治療晚期甲狀腺癌患者的新途徑仍是當前臨床醫(yī)生所關(guān)注的問題。

本研究基于GEO（Gene Expression Omnibus）數(shù)據(jù)庫中人甲狀腺癌數(shù)據(jù)集GSE50901和GSE138198，使用R軟件進行差異基因的篩選并利用TCGA數(shù)據(jù)庫對差異基因進行驗證，篩選出靶基因S100A1，隨后通過細胞和分子實驗探究S100A1對人甲狀腺癌SW579細胞增殖的影響及其可能的作用機制，有望為探索甲狀腺癌藥物治療的新靶點提供依據(jù)。

材料與方法

1 實驗材料 人甲狀腺癌SW579細胞系購自上海通派生物科技有限公司；DMEM培養(yǎng)基、青鏈霉素雙抗購自武漢普諾賽生命科技有限公司；RIPA裂解液購自上海碧云天生物技術(shù)有限公司；S100A1特異性siRNA序列（siS100A1）及對照組序列（siNeg）購自蘇州泓迅生物科技股份有限公司；PC3.1-S100A1-HA過表達載體和對照組載體（PC3.1-EGFP-HA）購自廣州云舟生物科技有限公司；HighGene轉(zhuǎn)染試劑購自武漢愛博泰克生物科技有限公司；Cell Counting Kit-8 （CCK-8）試劑盒和胰酶購自武漢賽維爾生物科技有限公司；抗S100A1一抗購自賽默飛世爾科技公司；抗β-catenin，c-myc和cyclinD1一抗均購自美國CST（Cell Signaling Technology）公司；抗GAPDH一抗購自武漢三鷹生物技術(shù)有限公司。

2 方法

2.1 細胞培養(yǎng)：人甲狀腺癌SW579細胞使用DMEM完全培養(yǎng)基（含10%胎牛血清和1%雙抗）常規(guī)培養(yǎng)，取對數(shù)生長期細胞進行后續(xù)細胞實驗。

2.2 細胞轉(zhuǎn)染及分組：按照轉(zhuǎn)染試劑說明書，將甲狀腺癌SW579細胞接種在6孔板中，“8”字法搖勻，細胞轉(zhuǎn)染時棄去培養(yǎng)基，PBS洗滌2次后加入不含血清的DMEM培養(yǎng)基，將siRNA或質(zhì)粒按照比例（3 μL siRNA或3 μg質(zhì)粒對應(yīng)5 μL轉(zhuǎn)染試劑）混勻，室溫靜置15 min后加入6孔板中，搖勻，轉(zhuǎn)染48 h后提取總蛋白驗證敲低及過表達效率，共包括4組，即敲低對照組（siNeg組），S100A1敲低組（siS100A1組），過表達對照組（EGFP組）和S100A1過表達組（S100A1組）。

2.3 細胞增殖檢測：將質(zhì)?；騭iRNA轉(zhuǎn)染48 h后的SW579細胞用胰酶消化重懸，均勻接種在96孔板中，調(diào)整細胞密度為5×103個/孔。待接種8h、24 h、48 h和72 h后，加入CCK-8試劑，用全自動酶標儀檢測各孔細胞吸光度值（OD Value，波長450 nm）。

2.4 Western Blot：轉(zhuǎn)染48 h后，用胰酶消化，離心后收集細胞，加入200 μLRIPA裂解液提取總蛋白；BCA法定量后，取20 μg總蛋白進行電泳分離；濕法轉(zhuǎn)膜及5%脫脂奶粉封閉后，根據(jù)目的蛋白分子大小裁膜，置入相應(yīng)的一抗（β-catenin，c-myc，cyclinD1和GAPDH）中孵育過夜；TBST洗膜后，加入辣根過氧化物酶標記的二抗，室溫孵育1h后使用ECL試劑盒進行顯影及曝光，使用ImageJ軟件進行定量分析。

2.5 RNA-Seq測序：收集轉(zhuǎn)染siNeg和siS100A1的SW579細胞，TRizol法提取總RNA，取0.5 μg總RNA進行逆轉(zhuǎn)錄和實時熒光定量PCR以檢測siRNA的敲低效率，其余總RNA置于-80℃保存，由深圳華大基因股份有限公司完成RNA-Seq測序，使用GSEA4.0.3軟件進行基因集富集分析（Gene Set Enrichment Analysis）。

2.6 差異基因篩選及富集分析：數(shù)據(jù)集GSE50901和GSE138198下載自GEO數(shù)據(jù)庫，使用R軟件及Limma包等功能包篩選差異基因，設(shè)置篩選條件為|log2（Fold Change）| ＞1且P.Val＜0.05，使用ggplot2、pheatmap等功能包進行繪圖和數(shù)據(jù)可視化。

3 統(tǒng)計學處理 應(yīng)用SPSS20.0統(tǒng)計軟件進行數(shù)據(jù)分析，結(jié)果以均數(shù)±標準誤（ ±SEM）表示，使用配對t檢驗進行統(tǒng)計學分析，P＜0.05表示差異有統(tǒng)計學意義。

結(jié) 果

1 甲狀腺癌生物信息學分析 選取GEO數(shù)據(jù)庫中甲狀腺癌數(shù)據(jù)集GSE50901，共含有65個甲狀腺組織樣本，其中正常組織4個，甲狀腺癌組織61個，對該數(shù)據(jù)集進行分析發(fā)現(xiàn)，與正常甲狀腺組織相比，甲狀腺癌組織中有454個基因表達水平發(fā)生改變，其中208個基因表達下調(diào)，246個基因表達上調(diào)（圖 1A）；選取GEO數(shù)據(jù)庫中數(shù)據(jù)集GSE138198，共含有17個甲狀腺組織樣本，其中正常組織3個，甲狀腺癌組織14個，對該數(shù)據(jù)集進行分析發(fā)現(xiàn)，與正常甲狀腺組織相比，甲狀腺癌組織中有2 324個基因表達水平發(fā)生改變，其中953個基因表達下調(diào)，1 371個基因表達上調(diào)（圖1B）。

2 S100A1在甲狀腺癌中高表達 對數(shù)據(jù)集GSE50901和GSE138198進行聯(lián)合分析發(fā)現(xiàn)，有97個基因在兩個數(shù)據(jù)集中均發(fā)生上調(diào)，46個基因均發(fā)生下調(diào)（圖2A），其中，S100A1基因在甲狀腺癌中高表達（圖2B和2C）。此外，對TCGA數(shù)據(jù)庫中甲狀腺癌樣本數(shù)據(jù)進行分析發(fā)現(xiàn)，與正常組織相比，甲狀腺癌組織的S100A1基因表達水平明顯上調(diào)（圖 2D）。

3 S100A1對甲狀腺癌SW579細胞增殖的影響 對人來源的正常甲狀腺細胞系HT-ori3和甲狀腺癌細胞系（TPC-1，K1，SW579，BCPAP）進行S100A1蛋白檢測，發(fā)現(xiàn)甲狀腺癌SW579細胞中S100A1蛋白水平顯著升高（圖 3A和3B）。通過siRNA干擾或者質(zhì)粒轉(zhuǎn)染的方式敲低或過表達SW579細胞的S100A1蛋白（圖 3C和3D），發(fā)現(xiàn)敲低S100A1能抑制SW579細胞的增殖，過表達S100A1則能促進SW579細胞的增殖（圖 3E和3F）。

圖1 甲狀腺癌組織生物信息學分析

圖2 S100A1在甲狀腺癌組織中高表達

圖3 S100A1在甲狀腺癌細胞系中的表達及對SW579細胞增殖的影響

4 S100A1對WNT/β-catenin信號通路的影響 為進一步探究S100A1促進甲狀腺癌SW579細胞增殖的作用機制，本研究采用RNA-Seq測序分析，對基因表達矩陣進行基因集富集分析（GSEA）發(fā)現(xiàn)，敲低SW579細胞中S100A1后，WNT信號顯著富集，且表現(xiàn)為信號通路活性降低（圖4A和4B）。Western blot結(jié)果表明，siS100A1組細胞中β-catenin，c-myc和cyclinD1的蛋白表達水平顯著低于siNeg組（圖4C和4E），而S100A1組細胞中β-catenin，c-myc和cyclinD1的蛋白表達水平顯著高于EGFP組（圖4D和4F）。

討 論

S100A1屬于鈣結(jié)合蛋白S100家族中的一員，其基因定位于人1號染色體長臂21區(qū)，在體內(nèi)主要以同源二聚體的形式存在[7]。S100A1蛋白羧基末端含Ca2+結(jié)合位點，通過與Ca2+結(jié)合，控制Ca2+內(nèi)流，引起細胞內(nèi)Ca2+穩(wěn)態(tài)失衡，導致下游一系列改變[8]；同時，S100A1蛋白與Ca2+結(jié)合后會發(fā)生蛋白構(gòu)象改變，暴露出疏水裂隙，通過與不同的靶蛋白結(jié)合，參與調(diào)控各種生物學過程[9]。既往研究表明，S100A1廣泛參與調(diào)控細胞收縮、細胞凋亡、細胞分化及基因轉(zhuǎn)錄等生物學過程[10，11]，而越來越多的研究發(fā)現(xiàn)，S100A1在多種惡性腫瘤中異常表達，參與腫瘤的發(fā)生發(fā)展，調(diào)控腫瘤細胞的增殖、轉(zhuǎn)移、耐藥等[12，13]。GUO[14]等發(fā)現(xiàn)，S100A1在肝癌中高表達，并與腫瘤分級、分期及預(yù)后密切相關(guān)，通過與LATS1蛋白相互作用，調(diào)控下游Hippo信號通路從而影響肝癌細胞的增殖。TIAN[15]等發(fā)現(xiàn)，S100A1在卵巢癌腫瘤樣本中高表達，并與卵巢癌淋巴結(jié)轉(zhuǎn)移密切相關(guān)，進一步細胞實驗結(jié)果表明，過表達S100A1能顯著提高卵巢癌細胞的遷移和侵襲能力。這些研究提示，S100A1可能作為癌基因參與調(diào)控惡性腫瘤的發(fā)生發(fā)展，與腫瘤細胞的增殖、遷移、侵襲等生物學行為密切相關(guān)，而目前關(guān)于S100A1對甲狀腺癌調(diào)控作用的報道較少。

本研究通過分析人甲狀腺癌數(shù)據(jù)集GSE50901和GSE138198，發(fā)現(xiàn)有97個基因在此兩個數(shù)據(jù)集中均發(fā)生上調(diào)，46個基因均發(fā)生下調(diào)，其中S100A1基因在甲狀腺癌中高表達。進一步分析TCGA數(shù)據(jù)庫中甲狀腺癌組織樣本的轉(zhuǎn)錄組數(shù)據(jù)，發(fā)現(xiàn)S100A1基因在腫瘤組織中明顯上調(diào)，提示S100A1基因的高表達可能參與甲狀腺癌的發(fā)生和發(fā)展，因此本研究選擇S100A1基因進行后續(xù)研究。為進一步探究S100A1基因在甲狀腺癌中的作用，本研究檢測了敲低或過表達S100A1對人甲狀腺癌SW579細胞增殖的影響。結(jié)果表明，敲低S100A1能抑制SW579細胞的增殖，而過表達S100A1則能促進SW579細胞的增殖，說明S100A1在甲狀腺癌的細胞增殖中具有重要作用。

圖4 S100A1對WNT/β-catenin信號通路蛋白表達的影響

為探究S100A1促進甲狀腺癌細胞增殖的作用機制，本研究通過siRNA干擾的方式敲低SW579細胞的S100A1蛋白表達，隨后進行RNA-Seq測序分析，對基因表達矩陣進行基因集富集分析發(fā)現(xiàn)，敲低SW579細胞的S100A1后，WNT信號顯著富集，且表現(xiàn)為信號通路活性降低。WNT信號通路是哺乳動物細胞中一條高度保守的信號轉(zhuǎn)導通路，由WNT配體蛋白和膜蛋白受體及下游信號分子組成，其中最為經(jīng)典的是WNT/β-catenin信號，通過一系列級聯(lián)反應(yīng)，可引起β-catenin蛋白入核，與細胞核中不同的轉(zhuǎn)錄因子結(jié)合調(diào)控下游靶基因Myc，CCND1等的表達[16]。既往研究表明，WNT/β-catenin信號通路在惡性腫瘤中異常激活，并與腫瘤的發(fā)生、發(fā)展密切相關(guān)，參與調(diào)控腫瘤細胞的增殖、分化、耐藥、轉(zhuǎn)移等生物學過程[17]。邵宇鑫[18]等發(fā)現(xiàn)β-catenin在分化型甲狀腺癌中異常表達，并與其進展和轉(zhuǎn)移密切相關(guān)。張貴軍[19]等研究表明，β-catenin及WNT/β-catenin信號通路下游蛋白cyclin D1在甲狀腺癌組織中高表達，并與腫瘤病理類型、臨床分期及轉(zhuǎn)移相關(guān)，這些研究均提示W(wǎng)NT/β-catenin信號通路激活與甲狀腺癌的發(fā)生、發(fā)展密切相關(guān)，而抑制WNT/β-catenin信號通路活性能發(fā)揮抗甲狀腺癌作用。因此，為進一步探究S100A1促進甲狀腺癌SW579細胞增殖的作用機制，本研究檢測了S100A1敲低或過表達對WNT/β-catenin信號通路活性的影響。Western blot結(jié)果表明，敲低S100A1能降低β-catenin，c-myc和cyclinD1蛋白表達水平，而過表達S100A1則能提高β-catenin，c-myc和cyclinD1的蛋白表達水平，說明S100A1蛋白可能是通過調(diào)控WNT/β-catenin信號通路活性影響甲狀腺癌SW579細胞的增殖。為了深入探討WNT/β-catenin信號通路在S100A1蛋白介導的甲狀腺癌細胞增殖中的作用，后續(xù)研究需要在敲低S100A1蛋白的同時激活WNT/β-catenin信號通路，或過表達S100A1蛋白的同時使用抑制劑等方式抑制WNT/β-catenin信號通路以探究能否逆轉(zhuǎn)S100A1蛋白對甲狀腺癌細胞增殖的影響。此外，本研究尚未探究S100A1調(diào)控WNT/β-catenin信號通路活性的具體機制，這也是后續(xù)進一步研究的重點。

總之，本研究表明，S100A1在甲狀腺癌中高表達，敲低S100A1可能是通過抑制WNT/β-catenin信號通路活性從而抑制甲狀腺癌細胞增殖，S100A1可能成為潛在的甲狀腺癌藥物治療的新靶點。

利益沖突所有作者均聲明不存在利益沖突