摘要

背景与目的:胰腺癌具有高度迁移性和侵袭性的特点,预后不良,但癌基因的突变及分子调控关系等的改变会影响胰腺癌的恶性生物学行为。泛素/ISG15结合酶E2L6 (UBE2L6)在泛癌中的表达及胰腺癌中的生物学功能依然不清楚。本研究拟通过生物信息学分析及实验验证,探讨UBE2L6对胰腺癌增殖、迁移侵袭能力的影响及可能的分子机制。方法:从UCSC Xena下载TCGA和GTEx的RNA-seq数据,用R (4.0.2)软件的“limma”包分析UBE2L6在泛癌中的差异表达。利用从GEO数据集下载数据及qRT-PCR分别验证UBE2L6在胰腺癌组织及细胞中的差异表达。胰腺癌细胞转染靶向UBE2L6的siRNA后用CCK-8细胞增殖、克隆形成、Transwell迁移和侵袭实验分析UBE2L6对胰腺癌细胞生物学功能的影响。利用UBE2L6的分子相关性、蛋白相互作用及基因富集分析(GSEA)探讨其可能的作用机制。进一步利用TCGA数据库的胰腺癌临床数据,分析UBE2L6的表达与其临床病理特征的关系及临床应用价值。结果:表达差异分析提示UBE2L6在泛癌中表达增高,在胰腺癌组织及BxPC-3 (t=33.82,P<0.000 1)、PANC-1 (t=7.36,P=0.001 8)、AsPC-1 (t=9.61,P=0.000 7)、SW1990 (t=10.26,P=0.000 5)及MIA PaCa-2(t=12.65,P=0.000 2)等胰腺癌细胞系中亦高表达。细胞功能实验显示,UBE2L6可促进胰腺癌PANC-1和AsPC-1细胞的增殖、迁移及侵袭。Spearman相关性分析提示BCL2A1蛋白与UBE2L6的相关性最大(r=0.442)。蛋白相互作用分析显示,其可能与HERC6、 RPS27A、 HERC5、 UBA52、 ISG15、 UBA7、DDX58、UBC、UBB、ARIH1等蛋白有相互作用。同时,GSEA富集分析提示UBE2L6与肿瘤免疫微环境等密切相关。而临床病理相关性分析显示,UBE2L6高表达与胰腺癌组织病理学分级密切相关(χ2=6.966,P=0.031)。单因素及多因素Cox回归分析表明,年龄及淋巴结转移是胰腺癌患者的独立危险因素(P<0.05),UBE2L6的表达水平与胰腺癌患者的预后明显相关(P<0.05),但不是其独立预后因素(P>0.05)。ROC曲线及Kaplan-Meier生存分析提示,UBE2L6对胰腺癌具有较高的临床诊断价值(AUC=0.970,95%CI=0.950~0.989)且与胰腺癌患者的预后相关(P<0.05)。结论:UBE2L6在胰腺癌中表达增高并促进胰腺癌细胞的增殖、迁移及侵袭,其分子机制可能与细胞凋亡、泛素化修饰、MAPK信号通路及肿瘤免疫微环境等有关,且UBE2L6表达水平与胰腺癌的诊断及预后密切相关。