張偉然 ，林雪峰，李 鑫，張 浩，王 猛，孫 偉，韓興鵬，孫大強(qiáng)

1. 天津醫(yī)科大學(xué)研究生院，天津 300070；2. 天津市胸科醫(yī)院胸外科，天津 300222；3. 天津醫(yī)學(xué)高等專(zhuān)科學(xué)校護(hù)理系，天津300222；4. 天津市南開(kāi)醫(yī)院胸外科，天津 300100

肺癌是世界上發(fā)病率和死亡率均排第一位的惡性腫瘤。據(jù)最新統(tǒng)計(jì)數(shù)據(jù)，2020 年全世界肺癌新發(fā)患者228 820 例，死亡病例135 720 例[1]。非小細(xì)胞肺癌（nonsmall-cell lung cancer，NSCLC）約占所有肺癌病例的85%[2]。肺腺癌是一種常見(jiàn)的NSCLC，約占NSCLC 的50%[3]。近年來(lái)，肺腺癌機(jī)制的研究得到了廣泛的開(kāi)展，研究人員發(fā)現(xiàn)遺傳因素影響了肺腺癌的發(fā)生和發(fā)展。其中有研究[4]發(fā)現(xiàn)，肺腺癌中最常出現(xiàn)的是ALK、ROS1 和RET 重排，另外還有80 多種基因易位。此外，還有很多常見(jiàn)的基因突變影響肺腺癌，包括TP53（46%）、EGFR（20%）、HER2（2%）、KRAS（1%）、ALK（1%）、BRAF（1%）和PIK3CA（1%）等[5]。這些發(fā)現(xiàn)對(duì)肺腺癌的靶向治療具有重要意義。例如，克唑替尼對(duì)ALK 或ROS1 融合的腫瘤有效，而吉非替尼、厄洛替尼和阿法替尼可用于EGFR 突變的腫瘤患者[6]。但目前臨床應(yīng)用的藥物很少，肺腺癌患者的死亡率很高，患者很難獲得長(zhǎng)期生存[7]。因此，肺腺癌的分子機(jī)制亟待進(jìn)一步探索。

目前，很多研究通過(guò)基因芯片數(shù)據(jù)分析揭示了一些與肺腺癌相關(guān)的基因。Tian 等[8]鑒定出9 個(gè)lncRNA 和6 個(gè)mRNA，它們可能是診斷或預(yù)測(cè)人類(lèi)肺腺癌細(xì)胞對(duì)紫杉醇耐藥性的有效標(biāo)志物。Liu 等[9]通過(guò)分析miRNA 基因芯片表達(dá)譜，發(fā)現(xiàn)miR-136、miR-376a 和miR-31 在小鼠肺癌中均顯著過(guò)表達(dá)，同時(shí)miR-31 作為致癌miRNA 在肺癌中靶向抑制特異性的腫瘤抑制因子。Zhu 等[10]報(bào)道了hsa_circ_0013958（一種環(huán)狀RNA）可以作為一種潛在的非侵入性生物標(biāo)志物用于肺腺癌的早期檢測(cè)和篩查。另一項(xiàng)微陣列圖譜研究[11]表明，hsa_circ_0007385 在NSCLC 腫瘤的形成中起重要作用，為NSCLC 的治療提供了一個(gè)潛在的靶點(diǎn)。本研究綜合分析了肺腺癌組織的基因、miRNA和circRNA 表達(dá)譜，以探索肺腺癌的潛在生物標(biāo)志物，并通過(guò)癌癥基因圖譜（the Cancer Genome Atlas，TCGA）對(duì)關(guān)鍵基因進(jìn)行驗(yàn)證。

1 資料與方法

1.1 基因芯片數(shù)據(jù)

芯 片 數(shù) 據(jù)GSE104854[12]、GSE102511、GSE85716、GSE83527 和GSE101586[13]從 基因表達(dá)數(shù)據(jù)庫(kù)（GEO，http://www.ncbi.nlm.nih.gov/geo/） 中 下 載。GSE104854、GSE102511 和GSE85716 是基因表達(dá)譜，分別包括3、16和6 個(gè)肺腺癌組織及其癌旁組織，測(cè)序平臺(tái)分別是Ion Torrent Proton（Homo sapiens）、Ion Torrent Proton（Homo sapiens）和Agilent-062918 OE Human lncRNA Microarray V4.0 028004。GSE83527 是一個(gè)基于Illumina Hiseq 2000（Homo sapiens）平臺(tái)的miRNA 表達(dá)譜，包含36 個(gè)肺腺癌樣本及其癌旁樣本。circRNA 表達(dá)譜芯片編號(hào)為GSE101586，有10 例樣本，包括5 例正常樣本及5 例癌癥樣本，芯片平臺(tái)為Agilent-069978 Arraystar Human CircRNA microarray V1。

1.2 差異表達(dá)分析

對(duì)所有的GEO 數(shù)據(jù)進(jìn)行處理后，通過(guò)Limma V3.34.8函數(shù)包（http://www.bioconductor.org/packages/3.6/bioc/html/limma.html）進(jìn)行差異表達(dá)基因（differentially expressed genes，DEGs）、差異表達(dá)miRNAs（differentially expressed miRNAs，DEMs） 和 差 異 表 達(dá)circRNAs（differentially expressed circRNAs，DECs）分析。GSE104854、GSE102511和GSE85716 的分析參數(shù)分別為|Log_FC|>3、|Log_FC|>2、|Log_FC|>2，均P<0.05。進(jìn)一步篩選出上述3 組表達(dá)譜的重疊DEGs。此外，DEMs 和DECs 的納入標(biāo)準(zhǔn)均為|Log_FC|>1 和P<0.05。采用circBase（http://www.circbase.org/）進(jìn)一步對(duì)DEGs 與DECs 的關(guān)系進(jìn)行分析。

1.3 功能富集分析

用注釋、可視化和綜合發(fā)現(xiàn)數(shù)據(jù)庫(kù)（David，https://david.ncifcrf.gov/summary.jsp，6.7 版）、京 都 基 因 和 基因組百科全書(shū)（KEGG PATHWAY，http://www.genome.jp/kegg）和Reactome（http://www.reactome.org）對(duì) 重 疊DEGs 進(jìn)行功能富集分析，得到P<0.05 的基因本體（gene ontology，GO）和重疊DEGs 所涉及的重要生物過(guò)程和 通路。

1.4 miRNA 靶向基因的篩選

為了探索肺腺癌樣本中DEMs 的調(diào)控機(jī)制，利用miRWall v3.0 數(shù) 據(jù) 庫(kù)（http://129.206.7.150/）獲 得DEMs的分子靶標(biāo)和miRNA-基因關(guān)系對(duì)，進(jìn)一步篩選出負(fù)調(diào)控的miRNA-基因關(guān)系對(duì)。

1.5 綜合調(diào)控網(wǎng)絡(luò)構(gòu)建

利用String V10.5（https://string-db.org）數(shù)據(jù)庫(kù)獲得重疊DEGs 的高置信（>400 分）的蛋白-蛋白對(duì)，基于蛋白-蛋白對(duì)和上述負(fù)調(diào)控的miRNA-基因關(guān)系對(duì)建立調(diào)控網(wǎng)絡(luò)，并用CytoCope V3.5.1 軟件（http://wwwyCytoCAPE.org/Sudio.php）進(jìn)行可視化。

1.6 關(guān)鍵基因的驗(yàn)證

從TCGA（https://cancer genome.nih.gov/）中下載15例肺腺癌患者的癌組織樣本及其癌旁組織進(jìn)行DEGs 篩選。TCGA 數(shù)據(jù)篩選參照中國(guó)臨床腫瘤學(xué)會(huì)腫瘤標(biāo)志物專(zhuān)家委員會(huì)和中國(guó)腫瘤驅(qū)動(dòng)基因協(xié)會(huì)的專(zhuān)家共識(shí)及相關(guān)文獻(xiàn)，由2 位病理學(xué)教授對(duì)每例樣本進(jìn)行病理定量和質(zhì)量控制分析。DEGs 以|Log_FC|>2、P<0.05 為閾值，記為DEGs-TCGA。將DEGs-TCGA 與上述基因表達(dá)譜的3 組DEGs 進(jìn)行比較。此外，采用聚類(lèi)算法對(duì)癌組織樣本和癌旁組織進(jìn)行分離，比較關(guān)鍵基因表達(dá)與患者總生存期和腫瘤分期間的相關(guān)性。

收集8 例2018 年1 月—10 月在天津市胸科醫(yī)院就診的肺腺癌患者的癌組織和癌旁組織，提取總RNA。使用SYBR?Premix Ex Taq ? Kit 試劑盒（TaKaRa，日本）進(jìn) 行SFTPC、ITLN2、SLC6A4、ADRA1A 以 及EGFR 的實(shí)時(shí)熒光定量PCR（qRT-PCR）檢測(cè)，hsa-mir-141-5p 和hsa-mir-191-3p 采 用All-in-One ? miRNA qPCR Kit 試 劑盒（GeneCopoeia Inc.，美國(guó)）檢測(cè)。所有引物均由寶日醫(yī)生物技術(shù)（北京）有限公司設(shè)計(jì)合成，如表1 所示。β-actin 和U6 為內(nèi)參。收集患者的外周血樣本，采用酶聯(lián)免疫吸附實(shí)驗(yàn)（ELISA）檢測(cè)血清中腫瘤標(biāo)志物癌胚抗原（CEA）和神經(jīng)元特異性烯醇化酶（NSE）的表達(dá)水平。

表1 SFTPC、ITLN2、SLC6A4、ADRA1A、hsa-mir-141-5p 和hsa-mir-191-3p 的引物序列Tab 1 Primer sequences of SFTPC, ITLN2, SLC6A4, ADRA1A, hsa-mir-141-5p and hsa-mir-191-3p

1.7 統(tǒng)計(jì)學(xué)分析

2 結(jié)果

2.1 DEGs、DEMs 和DECs 篩選

在GSE104854、GSE102511 和GSE85716 中， 與 癌旁組織相比，肺腺癌組織分別檢測(cè)到412 個(gè)（180 個(gè)上調(diào)和232 個(gè)下調(diào)）、201 個(gè)（71 個(gè)上調(diào)和130 個(gè)下調(diào)）和709 個(gè)（166 個(gè)上調(diào)和543 個(gè)下調(diào)）DEGs，圖1 顯示了3組DEGs 的Venn 圖。此外，3 組DEG 中包含61 個(gè)重疊DEGs（表2），其中P 值較小的重疊DEGs 包括SFTPC、SLC6A4、CA4、TMEM100 和ITLn2。這些重疊DEGs 具有相同的調(diào)節(jié)趨勢(shì)，即同時(shí)上調(diào)或下調(diào)。

圖1 GSE104854、GSE102511 和GSE85716 基 因 表 達(dá) 譜 中3 組DEGs 的Venn 圖Fig 1 Venn diagram of 3 sets of DEGs in gene expression profiles of GSE104854, GSE102511 and GSE85716

表2 GSE104854、GSE102511 和GSE85716 基因表達(dá)譜中的重疊DEGsTab 2 Overlapped DEGs in gene expression profiles of GSE104854, GSE102511 and GSE85716

Continued Tab

Continued Tab

在GSE83527 中，與癌旁組織相比，肺腺癌組織僅發(fā)現(xiàn)24 個(gè)DEMs，其中11 個(gè)上調(diào)，13 個(gè)下調(diào)（表3）。所有DEMs 在2 組間差異均有統(tǒng)計(jì)學(xué)意義（P<0.05）。

表3 肺腺癌組織相對(duì)于癌旁組織的DEMsTab 3 DEMs in lung adenocarcinoma tissues compared with their paracancerous tissues

Continued Tab

在GSE101586 中，肺腺癌組織和癌旁組織發(fā)現(xiàn)92 個(gè)DECs，包括10 個(gè)上調(diào)DECs 和82 個(gè)下調(diào)DECs。進(jìn)一步確定這些DECs 的位置，其中70 個(gè)位于編碼區(qū)的外顯子區(qū)，18 個(gè)位于編碼區(qū)的內(nèi)含子區(qū)，3 個(gè)位于內(nèi)部轉(zhuǎn)錄間隔區(qū)，1 個(gè)位于非編碼區(qū)的5'非翻譯區(qū)（5'UTR）。此外，篩選出與DECs 有關(guān)系的基因中，在上述3 組DEGs 中出現(xiàn)過(guò)的基因只有LDB2 和PREX1，而重疊DEGs 中沒(méi)有篩選出與DECs 相關(guān)的基因。前30 個(gè)差異最顯著的DECs 的位置和相應(yīng)的基因見(jiàn)表4。

表4 根據(jù)P 值排列的前30 個(gè)最顯著的DECs 及其位置Tab 4 Top 30 most significant DECs according to P value and their position and corresponding genes

2.2 功能通路富集

對(duì)DEGs 進(jìn)行功能和通路富集分析，共得到32 個(gè)通路，包括23 個(gè)生物過(guò)程、6 個(gè)細(xì)胞組分、1 個(gè)分子功能、1 個(gè)KEGG 通路和1 個(gè)RACTOME 通路（圖2）。從結(jié)果中可以看出大部分DEGs 顯著富集于系統(tǒng)發(fā)育、遷移調(diào)控、心率調(diào)控、血壓調(diào)控等生物學(xué)過(guò)程，主要富集于GPCR 及刺激神經(jīng)組織的受體配體相互作用信號(hào)通路。此外，富集到通路中最多的基因分別為S1PR1、EDN1 和ACVRL1，分別富集在24、24 和23 個(gè)通路中。

圖2 重疊DEGs 富集的GO 和通路Fig 2 Enriched GO and pathway terms of the overlapped DEGs

2.3 DEMs 的靶標(biāo)預(yù)測(cè)

根據(jù)芯片數(shù)據(jù)得到了24 個(gè)DEMs，并篩選出DEMs的靶基因，僅保留為DEGs 關(guān)系對(duì)的靶向基因。由此篩選出612 對(duì)miRNA-基因關(guān)系對(duì)，其中303 對(duì)為負(fù)調(diào)控。

2.4 重疊DEGs 與DEMs 的調(diào)控網(wǎng)絡(luò)

經(jīng)過(guò)String 數(shù)據(jù)庫(kù)分析，共鑒定出得分在400 分以上的19 對(duì)蛋白-蛋白相互作用對(duì)，涉及19 個(gè)重疊DEGs，EDNRB-EDN1、RAMP2-VIPR1 和RAMP3-VIPR1 是 打 分最高的3 個(gè)關(guān)系對(duì)（分別為997、947 和947 分）。基于19 個(gè)蛋白-蛋白對(duì)和以上303 個(gè)負(fù)調(diào)控的miRNA-基因?qū)?gòu)建了調(diào)節(jié)網(wǎng)絡(luò)（圖3），其中，hsa-mir-141-5p、hsa-mir-191-3p 和ADRA1A 的節(jié)點(diǎn)度最高，被認(rèn)為是潛在的生物標(biāo)志物。

2.5 DEGs 的驗(yàn)證

在TCGA 數(shù)據(jù)庫(kù)中共篩選出1 175 個(gè)DEGs，并與癌旁組織進(jìn)行比較，其中412 例上調(diào)，763 例下調(diào)。這些基因中有55 個(gè)基因（10 個(gè)上調(diào)，45 個(gè)下調(diào)）與上述的重疊DEGs（表2）相同并具有同樣的表達(dá)趨勢(shì)。其余6 個(gè)未鑒定的基因分別為T(mén)OX3、EDN1、SPDEF、FAM83A、FGFBP2 和COMP，說(shuō)明本研究所鑒定的DEGs 是可靠的。此外，SFTPC、ITLN2 和SLC6A4 的表達(dá)量較高，且在2 組DEGs 中均具有顯著差異。因此，選擇SFTPC、ITLN2 和SLC6A4 作為潛在的標(biāo)志物。

關(guān)鍵基因SFTPC、ITLN2、SLC6A4、hsa-mir-141-5p、hsa-mir-191-3p 和ADRA1A 的表達(dá)水平與患者總生存期和腫瘤分期間的相關(guān)性見(jiàn)表5。相關(guān)系數(shù)用來(lái)衡量2 個(gè)數(shù)據(jù)集合的線性關(guān)系，相關(guān)系數(shù)的絕對(duì)值越大，相關(guān)性越強(qiáng)。

圖3 重疊DEGs 和DEMs 的調(diào)控網(wǎng)絡(luò)Fig 3 Regulated network of overlapped DEGs and DEMs

經(jīng)過(guò)統(tǒng)計(jì)學(xué)分析，上述6 個(gè)基因中hsa-mir-191-3p 的表達(dá)與患者的腫瘤分期顯著相關(guān)（P=0.005），其相關(guān)系數(shù)為正，說(shuō)明基因的表達(dá)水平與腫瘤分期呈正相關(guān)；其他基因的表達(dá)與患者的總生存期和腫瘤分期間均無(wú)相關(guān)性（P>0.05）。

表5 SFTPC、ITLN2、SLC6A4、ADRA1A、hsa-mir-141-5p 和hsa-mir-191-3p 表達(dá)與總生存期和腫瘤分期的相關(guān)性Tab 5 Correlation of SFTPC, ITLN2, SLC6A4, ADRA1A, hsa-mir-141-5p and hsamir-191-3p expression with overall survival and tumor stage

RT-qPCR 結(jié) 果 如 圖4 所 示。癌 組 織 中SFTPC 和SLC6A4 的表達(dá)水平明顯低于癌旁組織，而hsa-mir-141-5p和hsa-mir-191-3p 的表達(dá)水平明顯升高（P<0.05）。ITLN2和ADRA1A 在癌組織中的表達(dá)水平雖然與癌旁組織存在差異，但兩者間差異無(wú)統(tǒng)計(jì)學(xué)意義（P>0.05）。EGFR 在癌組織中的表達(dá)水平高于癌旁組織，血清中CEA 和NSE 的表達(dá)水平分別為5.36 μg/L 和14.69 U/mL。

圖4 肺腺癌組織和癌旁組織中SFTPC、ITLN2、SLC6A4、ADRA1A、hsamir-141-5p、hsa-mir-191-3p 和EGFR 的相對(duì)表達(dá)水平Fig 4 Relative expression levels of SFTPC, ITLN2, SLC6A4, ADRA1A, hsa-mir-141-5p, hsa-mir-191-3p and EGFR in lung adenocarcinoma tissues and paracancer tissues

3 討論

GEO 數(shù)據(jù)庫(kù)是由美國(guó)國(guó)立生物技術(shù)信息中心（National Center for Biotechnology Information，NCBI）創(chuàng)建并維護(hù)的基因表達(dá)數(shù)據(jù)庫(kù)，其中收錄了世界各國(guó)研究機(jī)構(gòu)提交的高通量基因表達(dá)數(shù)據(jù)。GEO 數(shù)據(jù)庫(kù)的設(shè)計(jì)初衷是為了收集整理各種基因表達(dá)芯片數(shù)據(jù)，后來(lái)陸續(xù)加入了甲基化芯片、lncRNA、miRNA、CNV 芯片等各種芯片，甚至高通量測(cè)序數(shù)據(jù)。在本研究中，從GEO 數(shù)據(jù)庫(kù)中下載了3 個(gè)基因表達(dá)譜、1 個(gè)miRNA 表達(dá)譜和1 個(gè)circRNA表達(dá)譜，通過(guò)差異表達(dá)分析、功能富集分析、miRNA-基因關(guān)系對(duì)篩選以及綜合調(diào)控網(wǎng)絡(luò)構(gòu)建，最終篩選出hsamiR-141-5p、hsa-miR-191-3p、ADRA1A、SFTPC、ITLN2和SLC6A4 為肺腺癌潛在的標(biāo)志物。

hsa-mir-141-5p 和hsa-mir-191-3p 在肺癌中的作用尚不清楚，只有少數(shù)研究揭示了它們?cè)谄渌┌Y中的潛在作 用。Wang 等[14]確 定 了10 個(gè)mRNA（hsa-mir-141-5p、hsa-mir-1271-5p、hsa-mir-574-3p 等），對(duì) 人 類(lèi) 卵 巢 癌 組織和正常組織的敏感度和特異度分別為97%和92%。Plieskatt 等[15]報(bào)道了hsa-mir-141-5p 和其他miRNA 結(jié)合與viverrini 感染的膽管癌亞型有關(guān)。Wu 等[16]發(fā)現(xiàn)與非侵襲性垂體腺瘤相比，hsa-mir-191-3p 在侵襲性垂體腺瘤中的表達(dá)上調(diào)。ADRA1A 在精神分裂癥患者中受到廣泛關(guān)注，它是導(dǎo)致更嚴(yán)重代謝異常的危險(xiǎn)因素[17]。在腫瘤中，ADRA1A 是一種抑癌基因；與老年肝癌患者相比，年輕肝癌患者中ADRA1A 含量較少[18]。此外，ADRA1A 在子宮內(nèi)膜癌患者外周血中高表達(dá)，ADRA1A 的表達(dá)與國(guó)際婦產(chǎn)科協(xié)會(huì)（FIGO）分期和淋巴結(jié)轉(zhuǎn)移有關(guān)。因此，血清ADRA1A 的表達(dá)可用于評(píng)估生存率，并可能參與子宮內(nèi)膜癌的發(fā)病和轉(zhuǎn)移。本研究發(fā)現(xiàn)hsa-mir-141-5p、hsa-mir-191-3p 和ADRA1A 可能在肺腺癌中發(fā)揮一定的作用，并且調(diào)控網(wǎng)絡(luò)進(jìn)一步提供了可能的調(diào)控對(duì)，例如hsa-mir-141-5p 和hsa-mir-191-3p 可能通過(guò)靶向ADRA1A 發(fā)揮作用。

表面活性蛋白C（surfactant protein C，SFTPC）在肺癌中的研究已經(jīng)取得了一定的成果。它是肺表面活性物質(zhì)蛋白之一，與肺細(xì)胞的分化有關(guān)[19]。Saintigny 等[20]認(rèn)為SFTPC 在肺腺癌組織中下調(diào)，與本研究結(jié)果一致。2014年的一項(xiàng)研究[21]表明，SFTPC 的表達(dá)與NSCLC 較短的總生存期顯著相關(guān)。SLC6A4 是5-羥色胺轉(zhuǎn)運(yùn)體基因，其表達(dá)或甲基化與抑郁癥的發(fā)生、表型及預(yù)后密切相關(guān)[22-23]。同時(shí)，SLC6A4 變異與結(jié)直腸癌患者的低生存率相關(guān)[24]。此外，SLC6A4 變異是慢性阻塞性肺病合并肺癌的風(fēng)險(xiǎn)因素[25]。然而，并沒(méi)有發(fā)現(xiàn)另一潛在標(biāo)志物ITLN2 在肺癌中的作用。在本研究中，與正常組織相比，SFTPC、ITLN2和SLC6A4 在肺組織中的表達(dá)有顯著變化，提示它們可能在肺腺癌中起一定作用，但尚需實(shí)驗(yàn)室和臨床證實(shí)，這也是我們下一步要開(kāi)展的工作。

我們也試圖從circRNA 表達(dá)水平揭示肺腺癌的相關(guān)機(jī)制，并獲得了一些DECs。然而，只鑒定了3 組DECs中2 個(gè)相應(yīng)的基因LDB2 和PREX1，而沒(méi)有發(fā)現(xiàn)重疊的DECs。因此，推測(cè)circRNAs 的表達(dá)可能具時(shí)空特異性，可能與肺腺癌無(wú)關(guān)，具體原因尚需進(jìn)一步研究。

在本研究中，對(duì)于從GEO 數(shù)據(jù)庫(kù)中發(fā)現(xiàn)的重疊DEGs，通過(guò)TCGA 數(shù)據(jù)庫(kù)進(jìn)行了驗(yàn)證。TCGA 數(shù)據(jù)庫(kù)是由美國(guó)國(guó)家癌癥研究所和美國(guó)人類(lèi)基因組研究所聯(lián)合啟動(dòng)的一個(gè)項(xiàng)目。作為目前最大的癌癥基因信息數(shù)據(jù)庫(kù)，TCGA 收錄了各種人類(lèi)癌癥的臨床數(shù)據(jù)、基因組變異、mRNA 和miRNA 表達(dá)以及甲基化等數(shù)據(jù)，是腫瘤研究中資源最豐富、數(shù)據(jù)最權(quán)威的數(shù)據(jù)庫(kù)[26]。本研究得到的61 個(gè)重疊DEGs 中，大多數(shù)都在TCGA 數(shù)據(jù)庫(kù)中得到了驗(yàn)證，因此被認(rèn)為是準(zhǔn)確和可信的。只有6 個(gè)基因，即TOX3、EDN1、SPDEF、FAM83A、FGFBP2 和COMP 未被證實(shí)，它們?cè)谡{(diào)控網(wǎng)絡(luò)中沒(méi)有起到關(guān)鍵作用。此外，我們 對(duì)TCGA 數(shù) 據(jù) 中SFTPC、ITLN2、SLC6A4、hsa-mir-141-5p、hsa-mir-191-3p 和ADRA1A 的表達(dá)水平與患者總生存期和腫瘤分期進(jìn)行了相關(guān)性研究，其中hsa-mir-191-3p與患者的分期間存在顯著相關(guān)性。這種相關(guān)性很可能反映了hsa-mir-191-3p 對(duì)腫瘤進(jìn)展和細(xì)胞遷移的影響，有必要進(jìn)行更多的研究，以進(jìn)一步闡釋hsa-mir-191-3p 作為腫瘤分期標(biāo)志物的有效性及其在腫瘤細(xì)胞行為中的作用。其他基因的表達(dá)水平與患者的總生存期和腫瘤分期間均無(wú)相關(guān)性，但在癌組織和癌旁組織中的表達(dá)存在顯著差異，說(shuō)明這些基因在肺腺癌中的差異表達(dá)與患者的發(fā)病情況相關(guān)，與生存和預(yù)后無(wú)關(guān)。因此，這些基因在肺腺癌的臨床診斷和治療中可能具有潛在的應(yīng)用價(jià)值。

在對(duì)選定的分子標(biāo)志物SFTPC、ITLN2、SLC6A4、hsa-mir-141-5p、hsa-mir-191-3p 和ADRA1A 進(jìn)行驗(yàn)證時(shí)，癌組織中SFTPC 和SLC6A4 的表達(dá)水平明顯高于癌旁組織，hsa-mir-141-5p 和hsa-mir-191-3p 的表達(dá)水平明顯降低，這些結(jié)果與上述生物信息學(xué)分析結(jié)果一致。而ITLN2和ADRA1A 在癌組織中的表達(dá)水平雖然低于癌旁組織，但兩者間差異無(wú)統(tǒng)計(jì)學(xué)意義。芯片中的樣本來(lái)自各個(gè)國(guó)家和地區(qū)，包括白種人、亞洲人和非裔人群等，而本研究收集的病例大部分為本市患者，這種個(gè)體差異可能導(dǎo)致不同人群中基因表達(dá)水平的差異。因此，后續(xù)需要進(jìn)一步進(jìn)行臨床試驗(yàn)，對(duì)所選的標(biāo)志物進(jìn)行驗(yàn)證。此外，癌組織中EGFR 的表達(dá)水平明顯高于癌旁組織，這與以往的研究結(jié)果一致[27]。EGFR 是表皮生長(zhǎng)因子受體家族成員之一，EGFR 信號(hào)通路對(duì)細(xì)胞的生長(zhǎng)、增殖和分化等生理過(guò)程發(fā)揮重要的作用，EGFR 陽(yáng)性肺癌患者的生存期明顯短于EGFR 陰性患者[28]。CEA 和NSE 是肺癌中常用的腫瘤標(biāo)志物，在正常人體中的水平分別低于5.0 μg/L 和12.5 U/mL，在臨床中經(jīng)常可以發(fā)現(xiàn)肺腺癌患者血清CEA 及NSE水平升高[29]。本研究中檢測(cè)結(jié)果顯示CEA 和NSE 在肺腺癌患者中的表達(dá)水平分別為5.36 μg/L 和14.69 U/mL，這與以往的研究結(jié)果一致。

在本研究中，轉(zhuǎn)錄組測(cè)序和分析應(yīng)用于識(shí)別重要的miRNA 和與轉(zhuǎn)錄組測(cè)序相關(guān)的基因。ADRA1A、hsa-miR-141-5p、hsa-miR-191-3p、SFTPC、ITLN2 和SLC6A4 可能是肺腺癌的潛在標(biāo)志物。此外，以上6 個(gè)基因組合在臨床檢測(cè)中可能具有更高的準(zhǔn)確度、敏感度及特異度，為肺腺癌的早期診斷和靶向治療提供了理論依據(jù)和潛在的研究方向。然而，這些生物標(biāo)志物的篩選僅僅是通過(guò)生物信息學(xué)分析，還需要更多的臨床試驗(yàn)和更直接的證據(jù)進(jìn)行驗(yàn)證。

參·考·文·獻(xiàn)

[1] Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020[J]. CA Cancer J Clin, 2020, 70(1): 7-30.

[2] Pitroda SP, Stack ME, Liu GF, et al. JAK2 inhibitor SAR302503 abrogates PDL1 expression and targets therapy-resistant non-small cell lung cancers[J]. Mol Cancer Ther, 2018, 17(4): 732-739.

[3] Denisenko TV, Budkevich IN, Zhivotovsky B. Cell death-based treatment of lung adenocarcinoma[J]. Cell Death Dis, 2018, 9(2): 117.

[4] Huret JL. Lung: translocations in adenocarcinoma[J]. Atlas Genet Cytogenet Oncol Haematol，2015. DOI: 10.4267/2042/56443.

[5] McGuire S. World cancer report 2014. Geneva, Switzerland: world health organization, international agency for research on cancer, WHO press, 2015[J]. Adv Nutr, 2016, 7(2): 418-419.

[6] Ettinger DS, Wood DE, Aisner DL, et al. Non-small cell lung cancer, version 5.2017, NCCN clinical practice guidelines in oncology[J]. J Natl Compr Canc Netw, 2017, 15(4): 504-535.

[7] Huang WW, Zhu WZ, Mu DL, et al. Perioperative management may improve long-term survival in patients after lung cancer surgery: a retrospective cohort study[J]. Anesth Analg, 2018, 126(5): 1666-1674.

[8] Tian X， Zhang HY，Zhang BJ, et al. Microarray expression profile of long non-coding RNAs in paclitaxel-resistant human lung adenocarcinoma cells[J]. Oncol Rep, 2017, 38(1): 293-300.

[9] Liu X, Sempere LF, Ouyang HX, et al. MicroRNA-31 functions as an oncogenic microRNA in mouse and human lung cancer cells by repressing specific tumor suppressors[J]. J Clin Invest, 2010, 120(4): 1298-1309.

[10] Zhu XL, Wang XY, Wei SZ, et al. Hsa_circ_0013958: a circular RNA and potential novel biomarker for lung adenocarcinoma[J]. FEBS J, 2017, 284(14): 2170-2182.

[11] Jiang MM, Mai ZT, Wan SZ, et al. Microarray profiles reveal that circular RNA hsa_circ_0007385 functions as an oncogene in non-small cell lung cancer tumorigenesis[J]. J Cancer Res Clin Oncol, 2018, 144(4): 667-674.

[12] Ding XX, Zhang S, Li X, et al. Profiling expression of coding genes, long noncoding RNA, and circular RNA in lung adenocarcinoma by ribosomal RNAdepleted RNA sequencing[J]. FEBS Open Bio, 2018, 8(4): 544-555.

[13] Zong L, Sun QC, Zhang HP, et al. Increased expression of circRNA_102231 in lung cancer and its clinical significance[J]. Biomedecine Pharmacother, 2018, 102: 639-644.

[14] Wang L, Zhu MJ, Ren AM, et al. A ten-microRNA signature identified from a genome-wide microRNA expression profiling in human epithelial ovarian cancer[J]. PLoS One, 2014, 9(5): e96472. DOI: 10.1371/journal.pone.0096472.

[15] Plieskatt JL, Rinaldi G, Feng YJ, et al. Distinct miRNA signatures associate with subtypes of cholangiocarcinoma from infection with the tumourigenic liver fluke Opisthorchis viverrini[J]. J Hepatol, 2014, 61(4): 850-858.

[16] Wu SL, Gu Y, Huang YY, et al. Novel biomarkers for non-functioning invasive pituitary adenomas were identified by using analysis of microRNAs expression profile[J]. Biochem Genet, 2017, 55(3): 253-267.

[17] Cheng C, Chiu HJ, Loh El-W, et al. Association of the ADRA1A gene and the severity of metabolic abnormalities in patients with schizophrenia[J]. Prog Neuropsychopharmacol Biol Psychiatry, 2012, 36(1): 205-210.

[18] Wang HW, Hsieh TH, Huang SY, et al. Forfeited hepatogenesis program and increased embryonic stem cell traits in young hepatocellular carcinoma (HCC) comparing to elderly HCC[J]. BMC Genom, 2013, 14(1): 736.

[19] Setoguchi Y, Tamaoki F, Nakai M, et al. The frequency of mutations in the surfactant protein C gene (SFTPC) is much higher in pediatric interstitial lung disease (PILD) than either adult-onset familial or sporadic interstitial pneumonia (IP)[C]. San Diego: American Thoracic Society 2009 International Conference, 2009.

[20] Saintigny P, Massarelli E, Lin S, et al. CXCR2 expression in tumor cells is a poor prognostic factor and promotes invasion and metastasis in lung adenocarcinoma[J]. Cancer Res, 2013, 73(2): 571-582.

[21] Lindskog C, Fagerberg L, Hallstr?m B, et al. The lung-specific proteome defined by integration of transcriptomics and antibody-based profiling[J]. FASEB J, 2014, 28(12): 5184-5196.

[22] Dong C, Wong ML, Licinio J. Sequence variations of ABCB1, SLC6A2, SLC6A3, SLC6A4, CREB1, CRHR1 and NTRK2: association with major depression and antidepressant response in Mexican-Americans[J]. Mol Psychiatry, 2009, 14(12): 1105-1118.

[23] Philibert RA, Sandhu H, Hollenbeck N, et al. The relationship of 5HTT (SLC6A4) methylation and genotype on mRNA expression and liability to major depression and alcohol dependence in subjects from the Iowa Adoption Studies[J]. Am J Med Genet B Neuropsychiatr Genet, 2008, 147B(5): 543-549.

[24] Savas S, Hyde A, Stuckless SN, et al. Serotonin transporter gene (SLC6A4) variations are associated with poor survival in colorectal cancer patients[J]. PLoS One, 2012, 7(7): e38953. DOI: 10.1371/journal.pone.0038953.

[25] de Andrade M, Li Y, Marks RS, et al. Genetic variants associated with the risk of chronic obstructive pulmonary disease with and without lung cancer[J]. Cancer Prev Res (Phila), 2012, 5(3): 365-373.

[26] 謝龍祥, 閆中義, 黨藝方, 等. TCGA 數(shù)據(jù)庫(kù): 海量癌癥數(shù)據(jù)的源泉[J]. 河南大學(xué)學(xué)報(bào)(醫(yī)學(xué)版), 2018, 37(3): 223-228.

[27] 趙敏, 田暢, 邸鑫, 等. RBM5、EGFR 和KRAS 在非小細(xì)胞肺癌中的表達(dá)及其意義[J]. 中國(guó)實(shí)驗(yàn)診斷學(xué), 2019, 23(2): 303-308.

[28] Arfaoui A, Kriaa L, Znaidi N, et al. Over-expression of EGFR is closely correlated to poor prognosis in Tunisian patients with non-small cell lung adenocarcinoma[J]. J Immunoassay Immunochem, 2014, 35(3): 256-268.

[29] Peng Y, Wang Y, Hao XZ, et al. Utility of multiple increased lung cancer tumor markers in treatment of patients with advanced lung adenocarcinoma[J]. Chin J Lung Cancer, 2017, 20(10): 690-694.