安罗替尼在肝癌中耐药的差异mRNA筛选
DOI: 10.3969/j.issn.1001-5256.2021.02.022
Differentially expressed mRNA involved in the resistance of liver cancer to anlotinib
-
摘要:
目的 运用ceRNA芯片筛选可能参与肝癌细胞对安罗替尼耐药过程的mRNA。 方法 利用大剂量冲击联合低剂量诱导的方法建立对安罗替尼耐药的肝癌细胞,用CCK8实验进行验证耐药细胞在安罗替尼作用下的细胞增殖差异;运用ceRNA芯片检测耐药肝癌细胞与正常肝癌细胞的基因表达差异;运用实时荧光定量PCR(real-time PCR)对部分芯片测出的部分基因差异进行验证。计量资料两组间比较采用独立样本t检验,Kaplan-Meier法对肝癌样本的总生存期进行生存分析,log-rank检验比较生存率差异。芯片筛选结果使用Fisher精确检验。 结果 耐药肝癌细胞与正常肝癌细胞基因表达差异较大,通过缩减范围筛选出差异最大的10个基因进行分析。与耐药和肿瘤生长相关的基因有4个,分别为BIRC2、BIRC7、ABCC2、MAPK8。其中BIRC2、ABCC2、MAPK8表达水平下降(P值分别为0.001 4、0.001 2、0.011 8),BIRC7的表达水平增多(P<0.001)。real-time PCR的验证结论与芯片一致(t值分别为10.74、32.65、18.34、2.80,P值分别为0.000 4、0.000 1、0.000 1、0.044 8)。BIRC7的高表达与MAPK8的低表达对应显著减少的生存期(P值分别0.022 0、0.005 6)。 结论 BIRC2、BIRC7、ABCC2、MAPK8在对安罗替尼耐药的肝癌细胞中差异表达,可能参与了肝癌细胞对安罗替尼耐药的过程。 Abstract:Objective To screen out the mRNAs involved in the resistance of hepatoma cells to anlotinib using ceRNA microarray. Methods High-dose shock combined with low-dose induction was used to culture hepatoma cells resistant to anlotinib, and CCK8 assay was used to verify the difference in the proliferation of drug-resistant hepatoma cells treated by anlotinib. The ceRNA microarray was used to screen out the differentially expressed genes between drug-resistant hepatoma cells and normal hepatoma cells, and real-time PCR was used to verify the differentially expressed genes detected by some microarrays. the independent samples t-test was used for comparison of continuous data between two groups, and the Kaplan-Meier method was used to analyze the overall survival of hepatoma cells samples, and the log-rank test was used to compare survival rates. Fisher's exact test was used for chip screening. Results There was a significant difference in gene expression between drug-resistant hepatoma cells and normal hepatoma cells, and 10 genes with the greatest difference were screened out for analysis by reducing the range. There were 4 genes associated with drug resistance and tumor growth, i.e., BIRC2, BIRC7, ABCC2, and MAPK8. There were significant reductions in the expression levels of BIRC2, ABCC2, and MAPK8 (P=0.001 4, 0.001 2, and 0.011 8), and there was a significant increase in the expression of BIRC7 (P < 0.001). The results of real-time PCR were consistent with those of microarray (t=10.74, 32.65, 18.34, and 2.80; P=0.000 4, 0.000 1, 0.000 1, and 0.044 8). The high expression of BIRC7 and the low expression of MAPK8 were associated with the significant reduction in survival time (P=0.022 0 and 0.005 6). Conclusion BIRC2, BIRC7, ABCC2, and MAPK8 are differentially expressed between anlotinib-resistant hepatoma cells and normal hepatoma cells and may be involved in the resistance of hepatoma cells to anlotinib. -
Key words:
- Liver Neoplasms /
- Anlotinib /
- Drug Resistance, Neoplasm /
- RNA, Messenger
-
原发性肝癌是一种常见的恶性肿瘤,具有起病隐匿、进展迅速、恶性程度高、治疗难度大等特点[1]。索拉非尼作为晚期肝癌的一线治疗药物,能有效延长大多数患者的生存期[2]。但是,对于一些耐药的患者,索拉非尼很难取得良好的成效。安罗替尼是一种新型小分子多靶点酪氨酸激酶抑制剂,能有效抑制血管内皮细胞生长因子受体(vascular endothelial growth factor receptor,VEGFR)、Met等激酶,可以抑制肿瘤及其血管的生成和生长[3]。目前包括非小细胞肺癌、结直肠癌等多种癌症的临床试验正在开展[4]。通过与其他酪氨酸激酶进行比较发现, 安罗替尼对VEGFR2的抑制作用显示出极高的选择性[5]。但是安罗替尼在肝癌中的研究目前较少,尤其是关于其耐药方面的研究更是少有报道。本研究主要通过培养建立对安罗替尼耐药的肝癌细胞并通过ceRNA芯片进行比较,最后筛选出可能存在的与耐药相关的mRNA。为今后的临床和科研工作提供基础。
1. 材料与方法
1.1 材料
人肝癌细胞Huh7购自美国模式菌种收集中心;安罗替尼试剂购自正大天晴药业集团;CCK8试剂盒购自日本同仁公司;DMEM培养基、胎牛血清(FBS)、0.25%胰蛋白酶均购自GIBCO公司;引物由擎科生物公司合成;培养细胞总RNA提取试剂盒购自百序生物科技有限公司;细胞裂解液购自碧云天生物科技有限公司;SBC Human (4*180K) ceRNA芯片购自上海伯豪生物有限公司;PrimeScriptTM RT Reagent Kit购自日本TAKARA公司;2X PowerUp SYBR Green Master Mix购自武汉塞维尔科技有限公司。
1.2 方法
1.2.1 培养细胞
以45 ml DMEM培养液中添加5 ml FBS和溶于50 μl DMSO中不同浓度的安罗替尼为培养液;细胞置于37 ℃、5%CO2、85%~90%湿度的恒温培养箱中。每天至少观察细胞1次,当细胞生长状态良好、密度达到90%时进行传代。
1.2.2 Huh7耐药细胞的建立
采用大剂量冲击联合低剂量诱导的方法筛选耐药细胞,取处于对数期的人肝癌细胞Huh7置于终浓度为0.5 μmol/L安罗替尼的DMEM完全培养液中进行培养。密切观察1 d,若细胞状态不好,及时更换成不含安罗替尼的DMEM完全培养液来调整细胞的状态。若细胞状态好则常规更换培养液,正常消化传代,当细胞完全适应后,将细胞置于2 μmol/L安罗替尼的DMEM完全培养液中培养24 h来进行大剂量冲击。逐步提升药物浓度,直至Huh7细胞可以在药物浓度为6 μmol/L的DMEM完全培养液中可以正常增殖传代,并且在不添加药物培养2周后,仍可以在药物浓度为6 μmol/L的DMEM完全培养液中正常增殖传代,经过6个月的培养,获得对安罗替尼耐药的细胞Huh7/R,并将其维持在含2 μmol/L的安罗替尼培养液中。
1.2.3 CCK8细胞增殖实验
取消化好的处于对数期的肝癌细胞Huh7和Huh7/R,调整细胞浓度为3×104~5×104/ ml,向96孔板内每孔加入100 μl细胞悬液。设置空白组、对照组与实验组,每组设置5个复孔;边缘用PBS填充;于空白组与对照组中加入100 μl不含安罗替尼的培养液,实验组分别加入100 μl不同浓度安罗替尼的培养液;1~2 h后,取出96孔板置于全自动酶标仪中,测定450 nm光吸收值;测出的对照组吸收值以接近1最为理想;根据每孔所得的光密度值来计算增殖抑制率;增殖抑制率=(对照组OD-实验组OD)/(对照组OD-空白组OD)×100%。
1.2.4 ceRNA芯片筛选差异基因
采用试剂盒的操作步骤对正常肝癌细胞Huh7和耐药肝癌细胞Huh7/R进行RNA的初步提取;对提取的RNA进行纯化;纯化后的RNA进行反转录和荧光标记;对杂交后的芯片进行洗涤和扫描;对差异的mRNA进行Go和KEGG pathway富集分析。
1.2.5 实时荧光定量(real-time PCR)反应
用RNA提取试剂盒对处于对数期的正常肝癌细胞株Huh7和耐药肝癌细胞株Huh7/R进行RNA的提取。按照PrimeScriptTM RT Reagent Kit说明书将RNA进行反转录成cDNA。以GADPH作为内参,按照2X PowerUp SYBR Green Master Mix说明书进行real-time PCR反应,条件为:95 ℃预变性30 s,而后进行40个循环(95 ℃ 15 s,60 ℃ 10 s,72 ℃ 30 s)。引物序列5′-3′:BIRC2(上游引物-CATAGTAGCTTGTTCAGTGGT,下游引物-CTCTAGAATTAAGAGGGTTTGGAG),ABCC2(上游引物-GGACACTCTTACAGGGTGAC,下游引物-GGATCAGGATCTGGAATCCG),MAPK8(上游引物-GATGTTTACATAGTCATGGAGCTC,下游引物-TACTGGGCTTTAAGTCCCG),BIRC7(上游引物-CTCCTTCTATGACTGGCCG,下游引物-CCTCACCTTGTCCTGATGG),GADPH(上游引物-TCAAGATCATCAGCAATGCC,下游引物-CGATACCAAAGTTGTCATGGA)。
1.3 统计学方法
使用SPSS 21.0软件对数据进行统计分析。计量资料以x±s表示,两组间比较采用独立样本t检验。采用Kaplan-Meier法对肝癌样本的总生存期进行生存分析,log-rank检验比较生存率差异。芯片筛选结果使用Fisher精确检验。P<0.05为差异有统计学意义。
2. 结果
2.1 肝癌细胞的抑制率
正常肝癌细胞24 h的IC50为(2.99±0.93)μmol/L,耐药肝癌细胞的24 h IC50为(27.85±1.42)μmol/L。通过公式:耐药指数=耐药细胞IC50/亲本细胞的IC50,计算出耐药肝癌细胞相对于正常肝癌细胞对安罗替尼的耐药指数为9.31(图 1)。
2.2 差异mRNA的筛选
芯片的数据来自R/bioconductor的数据包clusterProfiler,使用Fisher精确检验,筛选标准为term/GO上差异的基因数目不小于2,P<0.05。研究发现共有2790个差异mRNA,其中上调mRNA有611个,下调mRNA有2179个(图 2)。
2.3 Go和KEGG pathway分析
通过分析term/GO发现差异基因分布最多的3个基因集为cell、cell part和cellular progress(图 3)。KEGG中基因改变最多的通路为信号传导通路(图 4)。
2.4 前10个差异mRNA
通过缩减P值筛选出10个差异最大的基因分别为BIRC2、BIRC7、ABCC2、MAPK8、MAP2K6、SERPING1、MX1、RPP30、UTP14A、KLHL13。通过文献检索,发现有4个基因可能与耐药有关,分别是:BIRC2、BIRC7、ABCC2、MAPK8。其中BIRC2、ABCC2、MAPK8表达水平下降(P值分别为0.0014、0.0012、0.011 8),BIRC7的表达水平增加(P<0.001)。
2.5 real-time PCR验证差异表达mRNA
利用传统的real-time PCR对这BIRC2、BIRC7、ABCC2、MAPK8 4个mRNA进行了再次验证,并通过GraphPad Prism 8软件进行处理,结论与芯片一致(t值分别为10.74、32.65、18.34、2.80,P值分别为0.000 4、0.000 1、0.000 1、0.044 8)(图 5)。
2.6 生存分析
对4个差异最大的mRNA进行生存分析,其中BIRC7的高表达与MAPK8的低表达对应显著减少的累积生存率(P值分别为0.022 0、0.005 6)(图 6)。
3. 讨论
肝癌治疗方法主要有肝切除术、肝移植、局部射频消融术、经肝动脉化疗栓塞术(TACE)、放疗和系统治疗等[6]。但是由于肝癌起病隐匿,很多患者发现时已处于肝癌晚期。安罗替尼作为一种新型的小分子酪氨酸激酶抑制剂,目前有研究[7-9]表明原发性肝癌患者在TACE术后行安罗替尼治疗安全有效,并且不良反应可以耐受,这为患者带来了新的希望。然而,一旦出现安罗替尼耐药,晚期患者的治疗将面临更大的困难。
本研究通过大剂量冲击联合低剂量诱导的方法成功培养出对安罗替尼耐药的肝癌细胞,并通过ceRNA芯片对耐药肝癌细胞与正常肝癌细胞进行测序比较发现共有差异mRNA 2790个,其中上调mRNA有611个,下调mRNA有2179个。通过Go和KEGG pathway分析发现,差异mRNA分布最多的3个基因集为cell、cell part和cellular progres,KEGG中基因改变最多的通路为信号传导通路。通过缩小P值筛选出差异性较大的10个mRNA中有4个mRNA与耐药和细胞生长与死亡有关,分别为BIRC7、BIRC2、ABCC2、MAPK8。
杆状病毒IAP重复序列7(BIRC7)是最近发现的一个凋亡抑制蛋白家族成员。BIRC7的过表达与凋亡、化疗和放疗抵抗力增加、肿瘤复发和患者生存率下降有关[10-11]。杆状病毒IAP重复序列2(BIRC2)与BIRC7同属凋亡抑制蛋白家族。在NF-κB的调控中发挥关键作用[12-13],同时,抗凋亡基因BIRC2和BIRC3是潜在的致癌基因染色体9A1扩增子的驱动因子,可以增强肿瘤细胞的活性并抑制其凋亡[14]。有研究[15]表明BIRC2不仅可以抑制HIV1的转录,还可被Smac模拟物靶向作用以促进病毒潜伏期的逆转。ABCC2包含32个外显子,其中4个主要在肝脏中表达[16-17],与Dubin-Johnson综合征和遗传性胆汁淤积性疾病有密切关系。研究[18]表明ABCC2与苯妥英钠、卡马西平、丙戊酸钠、拉莫三嗪、托吡酯和左乙拉西坦等抗癫痫药物的耐药机制相关。MAPK8,也称为c-Jun N末端激酶,是MAPK8家族中的一员[19],可以通过MAPK信号通路介导对替莫唑胺的抗性和胶质母细胞瘤细胞凋亡[20]。
综上所述,通过分析安罗替尼耐药的肝癌细胞与正常肝癌细胞的基因表达差异,从中筛选出了10个差异最大的mRNA,进一步研究发现,在耐药细胞中BIRC2、ABCC2、MAPK8表达水平下降,而BIRC7的表达水平增加;其中BIRC7高表达以及MAPK8低表达均与肝癌患者的生存率降低相关,而BIRC2与ABCC2的表达水平与肝癌患者生存率没有显著相关性,这对以后的实验和临床工作存在一定的借鉴意义。
-
[1] Bureau of Medical Administration, National Health Commission of the People's Republic of China. Guidelines for diagnosis and treatment of primary liver cancer in China (2019 edition)[J]. J Clin Hepatol, 2020, 36(2): 277-292. (in Chinese) DOI: 10.3969/j.issn.1001-5256.2020.02.007中华人民共和国国家卫生健康委员会医政医管局.原发性肝癌诊疗规范(2019年版)[J].临床肝胆病杂志, 2020, 36(2): 277-292. DOI: 10.3969/j.issn.1001-5256.2020.02.007 [2] WANG S, WU QW, LI XK, et al. Interpretation of guidelines for diagnosis and treatment of primary liver cancer in China (2019 edition)[J]. J Clin Hepatol, 2020, 36(5): 996-999. (in Chinese) DOI: 10.3969/j.issn.1001-5256.2020.05.009王姗, 吴庆旺, 李小科, 等. 《原发性肝癌诊疗规范(2019年版)》解读[J].临床肝胆病杂志, 2020, 36(5): 996-999. DOI: 10.3969/j.issn.1001-5256.2020.05.009 [3] SUN B, YANG XZ, XIE F, et al. Effect of sorafenib combined with transcatheter arterial chemoembolization and radiofrequency ablation on the survival of patients with hepatocellular carcinoma[J]. J Clin Hepatol, 2019, 35(1): 104-108. (in Chinese) DOI: 10.3969/j.issn.1001-5256.2019.01.019孙斌, 杨晓珍, 谢放, 等.索拉非尼联合经肝动脉化疗栓塞术和射频消融对肝细胞癌患者生存期的影响[J].临床肝胆病杂志, 2019, 35(1): 104-108. DOI: 10.3969/j.issn.1001-5256.2019.01.019 [4] HE Z, YU Y, NONG Y, et al. Hepatitis B virus X protein promotes hepatocellular carcinoma invasion and metastasis via upregulating thioredoxin interacting protein[J]. Oncol Lett, 2017, 14(2): 1323-1332. DOI: 10.3892/ol.2017.6296 [5] LIN B, SONG X, YANG D, et al. Anlotinib inhibits angiogenesis via suppressing the activation of VEGFR2, PDGFRβ and FGFR1[J]. Gene, 2018, 654: 77-86. DOI: 10.1016/j.gene.2018.02.026 [6] YE SL. Downstaging and conversion treatment of intermediate and advanced liver cancer should be taken seri-ously[J]. J Clin Hepatol, 2020, 36(2): 249-251. (in Chinese) DOI: 10.3969/j.issn.1001-5256.2020.02.001叶胜龙.重视中晚期肝癌的降期转化治疗[J].临床肝胆病杂志, 2020, 36(2): 249 -251. DOI: 10.3969/j.issn.1001-5256.2020.02.001 [7] WANG YM, MENG Y, SHI XN, et al. Analysis of clinical efficacy, safety and prognosis of anlotinib hydrochloride in the treatment of advanced primary liver cancer[J]. Chin J Hepatol, 2020, 28(7): 619-624. (in Chinese) DOI: 10.3760/cma.j.cn501113-20191020-00385王亚梦, 孟宇, 师霄楠, 等.盐酸安罗替尼治疗晚期原发性肝癌的临床疗效、安全性及预后分析[J].中华肝脏病杂志, 2020, 28(7): 619-624. DOI: 10.3760/cma.j.cn501113-20191020-00385 [8] PAN YF, FENG XL, YU X, et al. The clinical efficacy of transarterial chemoembolization combined with anlotinib in the treatment of advanced hepatocellular carcinoma[J]. Chin J Exp Surg, 2020, 37(3): 531-533. (in Chinese) DOI: 10.3760/cma.j.issn.1001-9030.2020.03.040潘延凤, 冯晓丽, 余雪, 等.肝动脉化疗栓塞联合安罗替尼治疗晚期肝癌的临床研究[J].中华实验外科杂志, 2020, 37(3): 531-533. DOI: 10.3760/cma.j.issn.1001-9030.2020.03.040 [9] YANG B, XIE H, WANG CP, et al. Study on Inhibition of arotinib on human intrahepatic cholangiocarcinoma cell line HCC-9810[J]. Chin J Med Guide, 2017, 19(12): 1389-1391. (in Chinese) DOI: 10.3969/j.issn.1009-0959.2017.12.037杨斌, 谢辉, 王春平, 等.安罗替尼对人肝内胆管细胞癌细胞系HCC-9810作用研究[J].中国医药导刊, 2017, 19(12): 1389-1391. DOI: 10.3969/j.issn.1009-0959.2017.12.037 [10] LIU K, YU Q, LI H, et al. BIRC7 promotes epithelial-mesenchymal transition and metastasis in papillary thyroid carcinoma through restraining autophagy[J]. Am J Cancer Res, 2020, 10(1): 78-94. [11] DOU H, BUETOW L, SIBBET GJ, et al. BIRC7-E2 ubiquitin conjugate structure reveals the mechanism of ubiquitin transfer by a RING dimer[J]. Nat Struct Mol Biol, 2012, 19(9): 876-883. DOI: 10.1038/nsmb.2379 [12] SAMANTA D, HUANG TY, SHAH R, et al. BIRC2 expression impairs anti-cancer immunity and immunotherapy efficacy[J]. Cell Rep, 2020, 32(8): 108073. DOI: 10.1016/j.celrep.2020.108073 [13] GYRD-HANSEN M, MEIER P. IAPs: From caspase inhibitors to modulators of NF-kappaB, inflammation and cancer[J]. Nat Rev Cancer, 2010, 10(8): 561-574. DOI: 10.1038/nrc2889 [14] MA O, CAI WW, ZENDER L, et al. MMP13, Birc2 (cIAP1), and Birc3 (cIAP2), amplified on chromosome 9, collaborate with p53 deficiency in mouse osteosarcoma progression[J]. Cancer Res, 2009, 69(6): 2559-2567. DOI: 10.1158/0008-5472.CAN-08-2929 [15] CAMPBELL GR, BRUCKMAN RS, CHU YL, et al. SMAC mimetics induce autophagy-dependent apoptosis of HIV-1-infected resting memory CD4+ T cells[J]. Cell Host Microbe, 2018, 24(5): 689-702. DOI: 10.1016/j.chom.2018.09.007 [16] CORPECHOT C, BARBU V, CHAZOUILLÈRES O, et al. Genetic contribution of ABCC2 to Dubin-Johnson syndrome and inherited cholestatic disorders[J]. Liver Int, 2020, 40(1): 163-174. DOI: 10.1111/liv.14260 [17] RAZALI RH, NOORIZHAB M, JAMARI H, et al. Association of ABCC2 with levels and toxicity of methotrexate in malaysian childhood acute lymphoblastic leukemia (ALL)[J]. Pediatr Hematol Oncol, 2020, 37(3): 185-197. DOI: 10.1080/08880018.2019.1705949 [18] WU L, LI Y, SONG Y, et al. A recurrent ABCC2 p.G693R mutation resulting in loss of function of MRP2 and hyperbilirubinemia in Dubin-Johnson syndrome in China[J]. Orphanet J Rare Dis, 2020, 15(1): 74. DOI: 10.1186/s13023-020-1346-4 [19] NGAN HL, LIU Y, FONG AY, et al. MAPK pathway mutations in head and neck cancer affect immune microenvironments and ErbB3 signaling[J]. Life Sci Alliance, 2020, 3(6): e201900545. DOI: 10.26508/lsa.201900545 [20] XU P, ZHANG G, HOU S, et al. MAPK8 mediates resistance to temozolomide and apoptosis of glioblastoma cells through MAPK signaling pathway[J]. Biomed Pharmacother, 2018, 106: 1419-1427. DOI: 10.1016/j.biopha.2018.06.084 期刊类型引用(0)
其他类型引用(1)
-