苯并［a］芘惡性轉化細胞THBEc1中miR-584-5p/miR-211-5p/叉頭框蛋白A1軸參與IDH1和HSPB1轉錄調控

馬雪，劉芷毓，邢云昆，姚碧云，傅娟玲，趙鵬

（1.浙江省疾病預防控制中心，浙江 杭州 310051；2.北京大學公共衛(wèi)生學院毒理學系，食品安全毒理學研究與評價北京市重點實驗室，北京 100191）

苯并［a］芘（benzo［a］pyrene，BaP）屬于多環(huán)芳烴類化合物，廣泛存在于煙草煙霧及空氣、水和土壤等環(huán)境中，長期暴露會導致肺癌等多種癌癥的發(fā)生［1-2］。BaP 經(jīng)代謝酶代謝生成二氫二醇、酚類、環(huán)氧化物和醌類等中間產(chǎn)物及終致癌物7，8-二氫二醇-9，10-環(huán)氧BaP 發(fā)揮毒作用，可能的機制涉及與芳烴受體結合，與DNA 結合導致遺傳損傷，引起微RNA（microRNA，miRNA）失調等表觀遺傳改變及干擾信號通路等［3-10］。Jin 等［11］提出了BaP 始于芳香受體活化，繼而經(jīng)不同關鍵事件誘發(fā)肺功能減退、肺纖維化和肺癌等肺損傷有害結局路徑。盡管BaP 致癌機制研究取得了一定進展，但BaP 致癌的遺傳和表觀遺傳機制尚未闡明。

研究表明，miRNA表達失調是BaP誘發(fā)肺癌過程中的重要特征，基于二代測序（next-generation sequencing，NGS）技術的miRNA 表達譜分析可為進一步闡明BaP 致癌的表觀遺傳學機制提供線索。BaP可誘導miR-137，miR-150，miR-142-5p，miR-122，miR-34c，miR-34b-5p，miR-29b，let-7a，miR-21 和miR-29等表達失調，影響免疫應答、炎癥反應、細胞增殖、遷移、侵襲和細胞周期等生物學過程［12-14］。前期開展的miRNA組測序研究發(fā)現(xiàn)，BaP惡性轉化細胞THBEc1與對照細胞16HBE的miRNA存在明顯差異。功能研究發(fā)現(xiàn)，THBEc1細胞中miR-152-3p，miR-142-5p和miR-211-5p表達下調可導致活化白細胞黏附分子（activated leukocyte cell adhesion molecule，ALCAM）上調，影響THBEc1 細胞增殖、遷移和裸鼠體內成瘤等，提示miR-152-3p/miR-142-5p/miR-211-5p-ALCAM 軸可能在BaP 致癌過程中發(fā)揮作用［6］。本課題組前期對BaP 惡性轉化人支氣管上皮細胞T-16HBE-C1（THBEc1）和16HBE細胞的蛋白質組差異進行分析發(fā)現(xiàn)，THBEc1細胞中叉頭框蛋白A1（forkhead-box A1，F(xiàn)OXA1）蛋白表達水平上調。通過建立FOXA1敲除單克隆細胞THBEc1-ΔFOXA1-c34 和敲除對照單克隆細胞THBEc1-ctrl 開展表型錨定研究發(fā)現(xiàn)，F(xiàn)OXA1在BaP 致癌中具有癌基因作用［15］。但FOXA1在THBEc1細胞中表達上調的機制和發(fā)揮癌基因作用的信號通路或網(wǎng)絡尚未闡明。鑒于miRNA 在基因轉錄后調控中的重要作用，推測BaP 可能通過誘導miRNA 表達改變導致FOXA1 蛋白表達上調，繼而FOXA1 通過調控靶基因的表達在腫瘤發(fā)生發(fā)展中發(fā)揮作用。本研究采用人支氣管上皮細胞16HBE和本課題組構建的細胞THBEc1，THBEc1-ctrl 和THBEc1-ΔFOXA1-c34，綜合利用生物信息學和分子生物學技術，篩選靶向FOXA1的miRNA 以及FOXA1 的潛在靶基因，識別以FOXA1 為軸的上下游分子，以期為進一步探究FOXA1 在BaP 致癌中的作用和可能機制提供線索。

1 材料與方法

1.1 試劑和主要儀器

MEM 培養(yǎng)基和胎牛血清（美國Gibco公司）；胰蛋白酶、青霉素和鏈霉素（美國Sigma公司）；TRIzol和Lipofectamine?2000試劑（美國Invitrogen公司）；FastKing cDNA 第一鏈合成試劑盒、SuperReal-PreMix Plus（SYBR Green）和Pro-Light HRP 化學發(fā)光檢測試劑盒（北京天根生化科技有限公司）；人FOXA1基因cDNA 表達質粒和對照質粒pCMV-未標記的陰性對照質粒（北京義翹神州科技股份有限公司）；潮霉素B（hygromycin B，中國上海阿拉丁公司）；一抗，兔抗人FOXA1單克隆抗體（美國Abcam公司）；兔抗人異檸檬酸脫氫酶1（isocitrate dehydrogenase 1，IDH1）多克隆抗體和兔抗人熱休克蛋白B1（heat shock protein B1，HSPB1）多克隆抗體（北京義翹神州科技股份有限公司）；鼠抗人β 肌動蛋白單克隆抗體、辣根過氧化物酶標記山羊抗兔IgG 抗體和辣根過氧化物酶標記山羊抗小鼠IgG 抗體（二抗）（北京中杉金橋生物技術有限公司）；雙鏈miRNA 模擬物和雙鏈miRNA 模擬物陰性對照（上海吉瑪制藥技術有限公司）。SCI-I65D CO2恒溫培養(yǎng)箱（日本ASTEC 公司）；CFX96TM實時熒光定量PCR儀（美國Bio-Rad公司）。

1.2 細胞構建和培養(yǎng)

人支氣管上皮細胞16HBE，由中山大學公共衛(wèi)生學院陳雯教授惠贈，本實驗室保存。BaP 惡性轉化細胞THBEc1 由本實驗室構建并保存［16］。敲除FOXA1基因的THBEc1 單克隆細胞THBEc1-ΔFOXA1-c34和質粒對照細胞THBEc1-ctrl分別為利用CRISPR/Cas9 技術將FOXA1 雙切口酶質粒和對照雙切口酶質粒轉染到THBEc1 細胞后篩選所獲的單克隆株，由本實驗室構建并保存［15］。細胞培養(yǎng)均使用含10% 胎牛血清的MEM 培養(yǎng)液，培養(yǎng)于37 ℃、5% CO2、飽和濕度培養(yǎng)箱中。

1.3 靶向FOXA1的miRNA的預測和篩選

利用TargetScan7.2 數(shù)據(jù)庫（http：//www.targetscan.org/vert_72/）和ENCORI 數(shù)據(jù)庫（https：//starbase.sysu.edu.cn/）預測靶向FOXA1 的miRNA。同時結合本課題組前期miRNA 測序結果：與16HBE細胞比較，在THBEc1細胞中351個miRNA表達下調（差異倍數(shù)＜0.5 且錯誤發(fā)現(xiàn)率＜0.05）［6］。采用Venny2.1.0（https：//bioinfogp.cnb.csic.es/tools/venny/）對數(shù)據(jù)進行匯總分析，使用Venn圖對數(shù)據(jù)庫預測靶向FOXA1的miRNA與THBEc1細胞中表達下調的miRNA取交集，綜合預測可能靶向調控FOXA1的miRNA。

1.4 實時熒光定量PCR（real-time quantitative PCR，RT-qPCR）檢測16HBE和THBEc1細胞中hsamiR-584-5p（miR-584-5p），hsa-miR-142-5p（miR-142-5p）和hsa-miR-211-5p（miR-211-5p）表達水平

收集對數(shù)生長期的16HBE和THBEc1細胞，利用TRIzol 試劑提取細胞總RNA。以總RNA 為模板用FastKing cDNA 第一鏈合成試劑盒合成cDNA，RT-qPCR 檢測miRNA 水平。以U6 作為內參，miRNA相對表達水平以2-ΔΔCt表示。引物序列見表1。

Tab.1 MicroRNA（miRNA）primer sequences for realtime quantitative PCR（RT-qPCR）

1.5 靶向FOXA1的miRNA鑒定

1.5.1 miRNA轉染

將THBEc1 細胞以每孔2.5×105密度接種于6 孔板，當細胞處于對數(shù)生長期且細胞密度達到60%時進行轉染。轉染前棄掉原培養(yǎng)基，用無血清和無雙抗的MEM 培養(yǎng)基潤洗細胞3 次，按照Lipofectamine?2000 轉染試劑操作說明書，THBEc1 細胞分別轉染miR-584-5p 模擬物，miR-142-5p 模擬物和miR-211-5p 模擬物，同時設空白對照組（blank control，BC）、無模擬物含Lipofectamine?2000 試劑的轉染試劑對照組（mock control，MC）和轉染陰性對照模擬物的陰性對照組（negative control，NC）。于轉染后48 和72 h 收集細胞，Western印跡法檢測FOXA1表達水平。

1.5.2 Western 印跡法檢測miRNA 轉染后THBEc1細胞中FOXA1表達水平

取1.5.1 中轉染的THBEc1 細胞，使用RIPA 裂解液提取細胞總蛋白，Bradford 法測定蛋白濃度。將30 μg蛋白樣品經(jīng)12% SDS-PAGE分離，再經(jīng)濕式電轉法轉印至硝酸纖維素膜。電轉后的硝酸纖維素膜經(jīng)5%脫脂奶粉室溫封閉2 h 后，分別與FOXA1 單克隆抗體（1∶3000）和β 肌動蛋白抗體（1∶3000）4 ℃孵育過夜，繼而與辣根過氧化物酶標記的二抗（1∶20 000）室溫孵育2 h。經(jīng)TBST 清洗后，使用Pro-Light HRP 化學發(fā)光檢測試劑盒檢測蛋白條帶，利用天能化學發(fā)光成像分析系統(tǒng)獲取圖像。利用Image J軟件分析每個條帶積分吸收光度值，以目的蛋白與內參蛋白的積分吸光度比值表示目的蛋白表達水平。

1.6 FOXA1潛在靶基因的預測和篩選

利用hTF（Database of Human Transcription Factor Targets）數(shù)據(jù)庫（http：//bioinfo.life.hust.edu.cn/hTFtarget#！/），并結合Harmonizonme 數(shù)據(jù)庫（https：//maayanlab.cloud/Harmonizome/）進入ENCODE 和JASPAR 轉錄因子靶基因數(shù)據(jù)庫綜合預測FOXA1 的潛在靶基因。采用Venny2.1.0（https：//bioinfogp.cnb.csic.es/tools/venny/）對數(shù)據(jù)進行匯總分析，使用Venn圖對hTF，ENCODE和JASPAR 數(shù)據(jù)庫預測的FOXA1 潛在靶基因與本課題組前期對THBEc1-ctrl 和THBEc1-ΔFOXA1-c34細胞轉錄組測序檢出的差異表達基因取交集，綜合預測FOXA1的潛在靶基因［17］。

1.7 RT-qPCR 檢測THBEc1-ctrl 和THBEc1 -ΔFOXA1-c34 細胞中FOXA1 潛在靶基因的表達水平

收集對數(shù)生長期THBEc1-ctrl 和THBEc1-FOXA1-c34 細胞，利用TRIzol 試劑提取細胞總RNA。以總RNA為模板，使用FastKing cDNA第一鏈合成試劑盒合成cDNA。通過RT-qPCR 檢測mRNA 水平，包括跨膜蛋白98（transmembrane protein 98，TMEM98）、IKAROS 家族鋅指2（IKAROS family zinc finger 2，IKZF2）、IDH1、腫瘤壞死因子受體相關因子5（TNF receptor associated factor 5，TRAF5）、骨形態(tài)發(fā)生蛋白2 型受體（bone morphogenetic protein receptor type 2，BMPR2）、HSPB1、Runt 相關轉錄因子2（RUNX family transcription factor 2，RUNX2）、golgin A7 家族成員B（golgin A7 family member B，GOLGA7B）和維甲酸相關孤核受體A（RAR related orphan receptor A，RORA）共9 個基因表達水平。以GAPDH 作為內參，mRNA 相對表達水平以2-ΔΔCt表示。引物序列見表2。

Tab.2 mRNA primer sequences for RT-qPCR

1.8 FOXA1過表達質粒穩(wěn)定轉染細胞的構建和鑒定

1.8.1 構建

將THBEc1-ΔFOXA1-c34 細胞以每孔1.5×105密度接種于24 孔板中培養(yǎng)24 h，細胞密度60%～70%時進行轉染。按照Lipofectamine?2000 轉染試劑操作說明分別轉染FOXA1過表達質粒和陰性對照質粒，同時設轉染試劑組。于轉染后48 h，使用含潮霉素100 mg·L-1的MEM 培養(yǎng)基進行抗性篩選，每2 d更換1次含潮霉素的培養(yǎng)基。潮霉素篩選10 d，轉染試劑組細胞全部死亡后，使用不含潮霉素的MEM 完全培養(yǎng)基繼續(xù)培養(yǎng)細胞。待細胞生長成克隆后，將全部克隆相繼轉到6 孔板和25 cm2培養(yǎng)瓶擴增培養(yǎng)，從而得到穩(wěn)定轉染FOXA1 過表達質粒和對照質粒的細胞。以THBEc1-ΔFOXA1-c34細胞作為親本細胞構建的穩(wěn)定表達FOXA1的細胞命名為THBEc1-ΔFOXA1-c34-oe 細胞，穩(wěn)定轉染對照質粒的細胞命名為THBEc1-ΔFOXA1-c34-ctrl細胞。

1.8.2 鑒定

收集對數(shù)生長期的THBEc1-ctrl，THBEc1-ΔFOXA1-c34，THBEc1-ΔFOXA1-c34-ctrl和THBEc1-ΔFOXA1-c34-oe 細胞，用Western 印跡法檢測細胞中FOXA1表達水平，方法同1.5.2。

1.9 RT-qPCR 檢測THBEc1-ΔFOXA1-c34-ctrl 和THBEc1-ΔFOXA1-c34-oe 細胞中FOXA1 潛在靶基因的表達水平

收集對數(shù)生長期THBEc1-ΔFOXA1-c34-ctrl和THBEc1-ΔFOXA1-c34-oe 細胞，用RT-qPCR 檢測FOXA1潛在靶基因TMEM98，IKZF2，IDH1，TRAF5，BMPR2，HSPB1，RUNX2，GOLGA7B和RORA的表達水平，方法同1.7。

1.10 Western 印跡法檢測THBEc1-ΔFOXA1-c34-ctrl 和THBEc1-ΔFOXA1-c34-oe 細胞中l(wèi)DH1和HSPB1的表達水平

收集對數(shù)生長期THBEc1-ctrl，THBEc1-ΔFOXA1-c34，THBEc1-ΔFOXA1-c34-ctrl和THBEc1-ΔFOXA1-c34-oe細胞，用Western印跡法檢測IDH1和HSPB1的表達水平。其中電轉后的硝酸纖維素膜經(jīng)5%脫脂奶粉室溫封閉2 h 后，分別與IDH1 多克隆抗體（1∶500）、HSPB1 多克隆抗體（1∶10 000）和β 肌動蛋白抗體（1∶3000）4 ℃孵育過夜，繼而與辣根過氧化物酶標記的二抗（1∶20 000）室溫孵育2 h。其余方法均同1.5.2。

1.11 統(tǒng)計學分析

數(shù)據(jù)均為3 次重復實驗的結果，以±s表示。利用SPSS 26.0 軟件對數(shù)據(jù)進行統(tǒng)計學分析，Graphpad Prism 8.0 軟件進行圖表繪制。兩組間比較采用Studentt檢驗，多組間比較采用單因素方差分析（ANOVA），組內兩兩比較采用Dunnettt檢驗。P＜0.05表示差異具有統(tǒng)計學意義。

2 結果

2.1 預測和篩選靶向調控FOXA1的miRNA

利用TargetScan 7.2和ENCORI數(shù)據(jù)庫分別預測到213 和145 個可能靶向調控FOXA1的miRNA。對數(shù)據(jù)庫檢索到的miRNA 和THBEc1 細胞（與16HBE比較）中下調的miRNA取交集，發(fā)現(xiàn)THBEc1細胞中下調的miR-584-5p，miR-142-5p 和miR-211-5p 潛在靶向調控FOXA1（圖1A）。RT-qPCR結果顯示，THBEc1細胞中上述3個miRNA的表達水平均低于16HBE細胞（P＜0.01）（圖1B）。

2.2 鑒定靶向調控FOXA1的miRNA

Western 印跡結果如圖2 所示。在轉染后48和72 h，BC，MC 和NC 組間FOXA1 蛋白表達水平無明顯差異；miR-584-5p 模擬物和miR-211-5p 模擬物組THBEc1 細胞中FOXA1 蛋白表達水平均低于NC 組（P＜0.01）；而miR-142-5p 模擬物組THBEc1 細胞中FOXA1 蛋白表達水平與NC 組間無明顯差異。提示miR-584-5p和miR-211-5p參與FOXA1轉錄后調控，兩者均可下調FOXA1 蛋白的表達，miR-584-5p 和miR-211-5p 下調可能是BaP惡性轉化細胞THBEc1 中FOXA1 蛋白表達水平上調的機制之一。而miR-142-5p 不參與FOXA1轉錄后調控。

Fig.2 Expression levels of FOXA1 after transfection of miRNA mimics in THBEc1 cells for 48 h（A）and 72 h（B）by Western blotting.See Fig.1 for the cell treatment.BC：blank control cells；MC：mock control cells treated with Lipofection? 2000 reagent；NC：negative control cells transfected with negative mimics miR-584-5p：THBEc1 cells transfected with miR-584-5p mimics；miR-142-5p:THBEc1 cells transfected with miR-142-5p mimics;miR-211-5p:THBEc1 cells transfected with miR-211-5p mimics.A2 and B2 were the semi-quantitative results of A1 and B1，respectively.IA：integrated absorbance.±s，n=3.＊＊P＜0.01，compared with NC group.

2.3 預測和篩選FOXA1的潛在靶基因

通過檢索hTF，JASPAR 和ENCODE 數(shù)據(jù)庫，分別檢索到17 450，1 442 和11 209 個基因含F(xiàn)OXA1 結合位點。將3 個數(shù)據(jù)庫預測到的基因與THBEc1-ΔFOXA1-c34 和THBEc1-ctrl 細胞間的轉錄組測序后發(fā)現(xiàn)的1 195個差異表達基因（different expression genes，DEG）取交集，共發(fā)現(xiàn)20個基因在3個數(shù)據(jù)庫中均被預測為FOXA1的潛在靶基因，并且在THBEc1-ΔFOXA1-c34和THBEc1-ctrl細胞間差異表達，結果如圖3A所示。

Fig.3 lntersection set of potential target genes of FOXA1 predicted by hTF，JASPAR，ENCODE databases and DEGs between THBEc1-ΔFOXA1-c34 and THBEc1-ctrl cells (A) and mRNA expression levels of indicated genes in THBEc1-ΔFOXA1-c34 and THBEc1-ctrl cells（B）.THBEc1-ΔFOXA1-c34：monoclonal cell lines obtained after THBEc1 cells knocked out FOXA1；THBEc1-ctrl：THBEc1 cells transfected with negative control plasmid.±s，n=3.＊P＜0.05，＊＊P＜0.01，compared with THBEc1-ctrl cells.

本課題組前期工作證實了FOXA1在BaP惡性轉化細胞THBEc1中發(fā)揮促進裸鼠體內成瘤及THBEc1細胞體外增殖和遷移的作用。在上述20個FOXA1的潛在靶基因中，TMEM98，IKZF2，IDH1，TRAF5，BMPR2，HSPB1，RUNX2，GOLGA7B和RORA參與腫瘤的發(fā)生發(fā)展有文獻支持，故將這9 個基因作為候選基因，進一步研究FOXA1敲除對其mRNA表達水平的影響。結果如圖3B 所示，F(xiàn)OXA1 蛋白表達缺失后，TMEM98，IKZF2，IDH1，TRAF5，BMPR2，HSPB1和RUNX2mRNA 表達水平均下調（P＜0.05，P＜0.01），GOLGA7B和RORAmRNA表達水平均上調（P＜0.01，P＜0.05）。

2.4 構建FOXA1 過表達細胞THBEc1-ΔFOXA1-c34-oe

Western 印跡結果（圖4）顯示，與THBEc1-ΔFOXA1-c34-ctrl 細胞相比，THBEc1-ΔFOXA1-c34-oe 細胞中FOXA1 蛋白表達恢復（P＜0.01）。但THBEc1-ΔFOXA1-c34-oe細胞中FOXA1表達水平低于THBEc1-ctrl細胞（P＜0.01）；THBEc1-ΔFOXA1-c34 和THBEc1-ΔFOXA1-c34-ctrl 細胞中均未檢出FOXA1 蛋白表達。由此表明，穩(wěn)定表達FOXA1 的THBEc1-ΔFOXA1-c34-oe細胞模型構建成功。

Fig.4 Expression levels of FOXA1 protein in THBEc1-ctrl，THBEc1-ΔFOXA1-c34，THBEc1-ΔFOXA1-c34-ctrl and THBEc1-ΔFOXA1-c34-oe cells by Western blotting（A）.THBEc1-ΔFOXA1-c34-ctrl：THBEc1-ΔFOXA1-c34 transfected with negative control plasmid；THBEc1-ΔFOXA1-c34-oe：THBEc1-ΔFOXA1-c34 transfected with human FOXA1 gene ORF cDNA clone expressional plasmid.B was the semi-quantitative result of A.±s，n=3.＊＊P＜0.01，compared with THBEc1-ctrl cells；##P＜0.01，compared with THBEc1-ΔFOXA1-c34-ctrl cells.

2.5 FOXA1過表達對THBEc1-ΔFOXA1-c34細胞中FOXA1潛在靶基因mRNA水平的影響

RT-qPCR 結果（圖5）所示，THBEc1-ΔFOXA1-c34-oe 細胞中IDH1和HSPB1mRNA 表達水平較對照細胞THBEc1-ΔFOXA1-c34-ctrl上調（P＜0.01），余7個基因的mRNA表達水平無明顯差異。

Fig.5 mRNA expression levels of indicated genes in THBEc1-ΔFOXA1-c34-ctrl and THBEc1-ΔFOXA1-c34-oe cells by RT-qPCR.±s，n=3.＊＊P＜0.01，compared with THBEc1-ΔFOXA1-c34-ctrl cells.

2.6 敲除FOXA1 對THBEc1 細胞HSPB1 和lDH1蛋白表達的影響

Western 印跡結果（圖6）所示，與THBEc1-ctrl細胞相比，THBEc1-ΔFOXA1-c34 細胞中IDH1 和HSPB1 表達水平顯著下調（P＜0.01）；與THBEc1-ΔFOXA1-c34-ctrl 細胞相比，THBEc1-ΔFOXA1-c34-oe 細胞中IDH1 和HSPB1 表達水平顯著上調（P＜0.01）。表明FOXA1 表達水平的改變可影響IDH1和HSPB1的表達。

Fig.6 Expressions of HSPB1 and lDH1 in THBEc1-ctrl，THBEc1-ΔFOXA1-c34，THBEc1-ΔFOXA1-c34-ctrl and THBEc1-ΔFOXA1-c34-oe cells by western blotting(A).See Fig.4 for the cell treatment.B and C were the semi-quantitative results of A.±s，n=3.＊＊P＜0.01，compared with THBEc1-ctrl cells；##P＜0.01，compared with THBEc1-ΔFOXA1-c34-ctrl cells.

3 討論

本研究結果表明，在THBEc1 細胞中表達下調的miR-584-5p 和miR-211-5p 對FOXA1基因發(fā)揮轉錄后調控作用。通過對FOXA1敲除細胞THBEc1-ΔFOXA1-c34 與對照細胞THBEc1-ctrl 進行mRNA表達差異分析，以及對THBEc1-ΔFOXA1-c34 細胞進行FOXA1 功能回補，發(fā)現(xiàn)了可能受FOXA1 調控的2 個潛在靶基因，IDH1和HSPB1，提示miR-584-5p / miR - 211 - 5p-FOXA1-IDH1 /HSPB1軸可能在BaP致癌中發(fā)揮作用。

目前已發(fā)現(xiàn)多種miRNA 在多種腫瘤中表達失調，且影響腫瘤細胞的增殖、遷移和凋亡等生物學事件。miRNA 表達失調可進一步導致靶基因表達失調，繼而影響靶基因參與的信號通路，這可能是miRNA 參與腫瘤發(fā)生發(fā)展過程的重要方式。鑒于miRNA 在腫瘤發(fā)生發(fā)展中的重要作用，本研究重點關注了FOXA1的miRNA 調控，以探究BaP 惡性轉化細胞THBEc1 中FOXA1 蛋白表達水平上調的可能機制。已有文獻報道，miR-132，miR-200a 和miR-100 等多種miRNA 可通過靶向負調控FOXA1參與甲狀腺癌、膠質瘤和乳腺癌等多種腫瘤細胞的增殖、遷移和凋亡等，從而發(fā)揮腫瘤抑制因子的作用［18-20］。本研究綜合數(shù)據(jù)庫檢索、miRNA 測序、RT-qPCR 以及miRNA 模擬物轉染實驗的證據(jù)，首次證實了miR-584-5p 和miR-211-5p 對FOXA1的負調控作用，說明miR-584-5p 和miR-211-5p 下調是導致THBEc1 細胞FOXA1 蛋白表達水平上調的原因之一。miR-584-5p和miR-211-5p已被證實在非小細胞肺癌（non-small-cell lung cancer，NSCLC）、骨肉瘤、胃癌、肝癌等腫瘤中表達下調，并通過靶向相關基因，發(fā)揮腫瘤抑制因子作用，從而影響腫瘤的發(fā)生發(fā)展［21-23］。如miR-584-5p 在NSCLC 中表達下調，并通過靶向負調控YKT6 和基質金屬蛋白酶14 抑制NSCLC 的遷移、侵襲和血管生成［21-22］。下調miR-584-5p 可導致胃癌細胞中WW 結構域E3 泛素蛋白連接酶1 過表達，從而激活腫瘤生長因子-β 信號通路，促進增殖并抑制凋亡［23］。另外，在肝癌細胞中miR-584-5p 亦參與對細胞周期蛋白依賴性激酶16 的調控［24］。有關miR-211-5p，已有文獻報道其參與調控腫瘤相關的增殖、遷移、侵襲、血管生成和凋亡等，在腎細胞癌、膠質瘤、肝癌、骨肉瘤和三陰性乳腺癌等腫瘤中常表達下調［25-29］。在腎細胞癌中，過表達miR-211-5p 可導致上皮間充質轉化（epithelial-mesenchymal transition，EMT）關鍵轉錄因子蝸牛家族轉錄抑制因子1（encoding snail family zinc finger 1，SNAI1）表達下調，從而抑制腫瘤細胞體外遷移和侵襲及體內轉移［25］。在膠質瘤細胞中，miR-211-5p可以負調控同源盒基因C8 的表達，影響遷移和增殖［26］。體外研究表明，miR-211-5p 可通過負調控?；o酶長鏈家族成員4抑制肝癌細胞增殖、遷移和侵襲［27］。在骨肉瘤中，miR-211-5p 可通過負調控富含脯氨酸蛋白11，抑制腫瘤細胞遷移、侵襲，同時促進腫瘤細胞凋亡［28］。在三陰性乳腺癌中，miR-211-5p可通過靶向負調控沉默信息調節(jié)因子1，抑制腫瘤細胞增殖、侵襲、遷移和轉移［30］。亦有證據(jù)顯示，miR-211-5p 在前列腺癌、宮頸癌、甲狀腺癌和膠質母細胞瘤中可分別對殼多糖酶3 樣蛋白1、性別決定區(qū)Y 框蛋白4、丙酮酸脫氫酶激酶5 和陰陽-1 發(fā)揮轉錄后調控作用［30-33］。在上述文獻報道的miR-211-5p 的靶基因中，SNAI1 在BaP 惡性轉化細胞THBEc1 中mRNA表達水平較對照細胞16HBE 上調1.48 倍（P＜0.01），提示miR-211-5p/SNAI1 軸可能是THBEc1細胞呈高遷移和侵襲表型的機制之一。值得注意的是，課題組前期工作已證實ALCAM 受miR-211-5p 調控，敲除ALCAM可抑制THBEc1 細胞體外增殖和裸鼠體內成瘤，但促進THBEc1 細胞的體外遷移和體內轉移［6］。本研究證實，F(xiàn)OXA1亦是miR-211-5p 的下游靶基因，而敲除FOXA1可抑制THBEc1 細胞體外增殖、遷移以及裸鼠體內成瘤和轉移［15］。由此推斷，以miR-211-5p 為靶點，通過上調其表達將可能同時對FOXA1和ALCAM發(fā)揮負調控作用，抑制腫瘤細胞增殖而不增加轉移風險，相關研究工作值得開展。

研究發(fā)現(xiàn)，敲除FOXA1可導致BaP 惡性轉化細胞THBEc1 轉錄組顯著改變［17］。上述發(fā)現(xiàn)雖然提示FOXA1 直接或間接參與多種生物學事件，但欲闡明其直接參與的信號通路或網(wǎng)絡，則有賴于對FOXA1靶基因的鑒定。本研究在對FOXA1敲除細胞THBEc1-ΔFOXA1-c34 與對照細胞THBEc1-ctrl轉錄組差異分析的基礎上，通過對THBEc1-ΔFOXA1-c34 細胞進行FOXA1 功能回補，發(fā)現(xiàn)IDH1和HSPB1可能是受FOXA1 調控的2 個潛在靶基因。

IDH1 主要催化異檸檬酸氧化脫羧生成α-酮戊二酸。IDH1基因編碼的蛋白質是在細胞質和過氧化物酶體中發(fā)現(xiàn)的NADP+依賴性IDH，在細胞質NADPH 產(chǎn)生中起著重要作用。IDH1異常表達見于多種腫瘤。研究發(fā)現(xiàn)，IDH1在原發(fā)性膠質瘤中表達上調，并且其過表達可促進腫瘤細胞遷移［34］。Chen 等［35］通過對139 例NSCLC 患者血清IDH1 水平檢測發(fā)現(xiàn)，NSCLC 患者血清中IDH1 水平顯著高于健康對照和良性肺部疾病患者，提示IDH1 可能作為NSCLC 早期診斷的潛在標志物。Yang 等［36］研究發(fā)現(xiàn)，IDH1在結直腸癌中高表達并指示高級別腫瘤，而且IDH1 的高表達與結直腸癌對5-氟尿嘧啶的耐藥性有關。Zheng 等［37］研究發(fā)現(xiàn)，高IDH1水平與高骨肉瘤患者預后不良有關，HSP90 ATP酶活性激活因子1 可通過正向調節(jié)IDH1 促進骨肉瘤細胞生長和轉移。Charitou 等［38］研究發(fā)現(xiàn)，在宮頸癌細胞中敲低叉頭框O蛋白后過表達IDH1，可逆轉α-酮戊二酸和NADPH/NADP+水平下調。本研究發(fā)現(xiàn)IDH1受FOXA1 正調控，且IDH1在BaP 惡性轉化細胞THBEc1中mRNA表達水平是對照細胞16HBE 的1.69 倍（P＜0.01），提示FOXA1 可能通過上調IDH1發(fā)揮其癌基因作用。

HSPB1是本研究發(fā)現(xiàn)的另一個受FOXA1調控的潛在靶基因。該基因編碼小熱休克蛋白家族的一個成員，亦稱HSP27。在應對環(huán)境應激時，HSPB1蛋白從細胞質轉移到細胞核，作為分子伴侶促進其他蛋白質的正確折疊，從而參與蛋白的質量控制［39］。HSPB1基因的表達與多種人類癌癥的不良臨床預后相關，HSPB1蛋白可能促進癌細胞增殖和轉移，同時保護癌細胞免于凋亡［40-41］。多種腫瘤中HSPB1表達失調已有文獻報道。Huang 等［42］研究發(fā)現(xiàn)，HSPB1 在NSCLC 患者癌組織和血清中均上調，而且HSPB1 上調與腫瘤分化不良、淋巴轉移和不良預后有關。Sun 等［43］研究發(fā)現(xiàn)，在腫瘤細胞中敲低HSBP1可促進鐵死亡誘導劑Erastin引起的鐵死亡，提示HSPB1 可能是腫瘤細胞鐵死亡的負調節(jié)因子，在鐵代謝中發(fā)揮作用。此外，研究亦發(fā)現(xiàn)，HSPB1 與胃癌、肝癌、前列腺癌、肺癌和結直腸癌等多種腫瘤的耐藥性有關［44］。如在卵巢癌細胞中，HSPB1的過表達通過抑制P21從細胞核轉位至細胞質，從而誘導細胞對順鉑的耐藥性［45］。在前列腺癌中，HSPB1通過調節(jié)STAT3/Twist參與白細胞介素6 介導的EMT 從而驅動EMT，而其減弱逆轉EMT 并降低細胞遷移、侵襲和基質金屬蛋白酶活性［46］。在胃癌細胞中，F(xiàn)YN 過表達且可以通過FYN/殺傷性T 細胞來源的蛋白激酶HSPB1 軸介導細胞增殖和侵襲［47］。本研究結果表明，敲除FOXA1后HSPB1表達下調，而FOXA1回補可使HSPB1表達上調，提示FOXA1對HSPB1具有調控作用。值得注意的是，與對照細胞16HBE 比較，在FOXA1 表達上調的THBEc1 細胞中，HSPB1mRNA 表達水平無明顯改變（差異倍數(shù)為0.87，P=0.283），但HSPB1 蛋白表達水平顯著降低（差異倍數(shù)為0.14，P=0.001），提示除FOXA1 外，可能其他機制如miRNA、泛素降解系統(tǒng)等對HSPB1 的調控在THBEc1 細胞中發(fā)揮主導作用，相關機制值得深入研究。

綜上所述，BaP 惡性轉化細胞THBEc1 中miR-584-5p/miR -211-5p-FOXA1 軸參與IDH1和HSPB1轉錄調控。BaP 可能在誘發(fā)細胞惡性轉化過程中引起miR-584-5p 和miR-211-5p 表達下調，繼而導致FOXA1 蛋白表達上調。上調IDH1可能是FOXA1 在THBEc1 細胞中發(fā)揮癌基因作用的機制之一。BaP 誘導miR-584-5p 和miR-211-5p 表達下調的機制及IDH1 在THBEc1 細胞惡性表型中的作用仍需進一步研究。