基于生物信息学分析核仁纺锤体相关蛋白1在肝细胞癌中的表达及对临床预后的影响

靳鑫旻，杨通旺 ，徐庆国，刘 欢，关 鸽，臧运金

1 青岛大学 医学部，山东 青岛 266000；2 青岛大学附属医院 器官移植中心，山东 青岛 266000

肝癌作为消化系统最常见的恶性肿瘤之一，其发病率与致死率在各类肿瘤中一直位居前列[1]。我国每年肝癌新发病例数与死亡病例数约占全球总病例数的一半[2]，其中肝细胞癌占绝大多数，因此本文主要讨论肝细胞癌。肝癌的发展，除因环境因素影响外，往往还与基因的突变或异常表达有关[3]。由于肝癌具有早期症状不明显、进展隐匿等特点，尽管在肝癌发现后可以进行手术为主的多种治疗方法，但总体预后不佳[4-6]。肿瘤发生与进展的关键在于正常的细胞周期被破坏，失去分化的能力，最终发生恶性增殖。因此，可以将抑制肝癌细胞的恶性增殖、恢复其正常的细胞进程作为治疗肝癌的潜在靶向。

核仁纺锤体相关蛋白1(nuclelar spindle-associated protein 1,NUSAP1)作为一种微管结合蛋白，在人体内高度保守，其作用在于调控纺锤体的正确装配以及染色体的形成，在有丝分裂过程中发挥不可或缺的作用[7-9]。在正常的组织中，NUSAP1的表达水平受到严格的调控，以防止细胞增殖异常。当NUSAP1表达异常增高时，会导致细胞有丝分裂进程无法终止，引起细胞的无限增殖以及遗传物质的不稳定。现有研究[10-16]已经证实，NUSAP1在多种肿瘤中表达异常增高，提示其密切参与肿瘤的发生与发展。

生物信息学分析在近几年越来越多的用于研究肿瘤进展过程中的基因表达变化[17]。通过对大量共享数据的挖掘、分析及解释，可以方便、快捷地发现具有表达差异的基因，找到疾病治疗的潜在靶标以及预后评估的可能标志物[18]。本研究通过对GEO及TCGA数据库公共数据的整理分析，对NUSAP1在肝癌组织及细胞亚群中的表达差异及其对预后的影响进行了初步的探索。

1 资料与方法

1.1 数据收集与筛选 从TCGA数据库下载肝癌转录组信息及临床数据，选择临床信息完善的数据进行基因表达差异分析，以剔除肿瘤复发类型为非肝癌的R0切除术后患者的临床信息评估其预后。在GEO官网下载GSE57957、GSE14520、GSE22058、GSE46444、GSE54236、GSE36376、GSE64041、GSE76297、GSE76427、GSE102079共10个种属为Homo sapiens的肝癌数据集。在GEO官网下载单细胞测序数据集GSE149614。选取2018年1月—2019年11月在青岛大学附属医院进行同种异体肝移植术的随访信息完善、病理诊断明确、生物样本取材规范的乙型肝炎相关肝癌患者42例，统计其临床资料，包括性别、血型、年龄、BMI、血清AFP水平、MELD评分、Child-Pugh分级、肿瘤直径、肿瘤数目、TNM分期及Okuda分期等。

1.2 实时荧光定量PCR 对选取的42例临床样本利用Trizol法提取总RNA，并通过TaKaRa反转录试剂盒合成cDNA。利用TaKaRa PCR试剂盒进行实时荧光定量-PCR反应，以GAPDH为内参照，测定NUSAP1相对表达水平。GAPDH引物序列为：F：5′-CGGGCTCTCCAGAACATC-3′，R：5′-ATGACCTTGCCCACAGCCT-3′；NUSAP1引物序列为F：5′-ACCCTGACTCACAGCAGAATC-3′，R：5′-CAGCTTCTTGGTGCTCGTCT-3′。

1.3 伦理学审查 本研究方案经由青岛大学附属医院伦理委员会审批，批号：QYFYWZLL26046。

1.4 统计学方法 R 4.0软件用于所有数据集基因表达差异分析及单细胞测序数据的聚类分析。采用SPSS 21.0软件对数据进行统计分析。计数资料2组间比较采用χ2检验。以Kaplan-Meier方法生存分析NUSAP1表达水平与总体生存率及无病间隔生存率之间的关系，运用log-rank进行检验，通过Graphpad prim 7使其结果可视化。P<0.05为差异有统计学意义。

2 结果

2.1 NUSAP1在肝癌组织及癌旁组织中的表达 为研究NUSAP1在肝癌组织中的表达，同时从GEO和TCGA中下载肝癌和癌旁组织基因表达数据。结果显示，NUSAP1在10个GEO肝细胞数据集中的肝癌组织中高表达，在TCGA数据库中也呈高表达(P值均<0.001)(附录1)。

2.2 NUSAP1在肝癌细胞中的表达 为避免肿瘤异质性对本研究的影响，从GEO中下载了肝细胞癌的单细胞测序数据集GSE149614并进行聚类分析。依据细胞标志物和组织来源的差异，将所有细胞分为来源于肿瘤组织的细胞与来源于癌旁组织的细胞。结果显示，来源于肿瘤组织的细胞其NUSAP1高表达细胞明显高于癌旁组织细胞(附录2)。

2.3 NUSAP1在临床肝癌样本中的表达 数据库的测序结果，通过实时荧光定量PCR检测42例患者肝癌和癌旁组织中NUSAP1的相对表达。与GEO和TCGA数据库分析结果一致，NUSAP1在29例患者的肝癌组织中高表达，在13例患者的肝癌组织中低表达(P=0.003)(图1)。

图1 NUSAP1在临床肝癌样本中高表达

2.4 NUSAP1表达水平与肝癌患者临床病理特征的关系

NUSAP1在肝细胞癌组织中高表达组及低表达组间，Child-Pugh分级、TNM分期、Okuda分期比较差异均有统计学意义(P值均<0.05)(表1),表明NUSAP1的表达水平与肝癌患者的病情严重程度有关。

表1 NUSAP1表达水平与肝癌患者临床病理特征的关系

2.5 NUSAP1在肝细胞癌中高表达与预后的关系 NUSAP1高表达的肝癌患者总体生存率与无病间隔生存率均显著低于低表达的肝癌患者(图2)，NUSAP1低表达患者1、3、5年生存率分别为92.96%、83.80%、76.76%，明显高于高表达的患者(76.76%、67.60%、64.78%)(χ2值分别为8.156、7.853、7.671,P值均<0.05)，NUSAP1高表达的患者中位生存时间明显低于低表达的患者(61.7 个月 vs 108.6个月)(χ2=7.955,P<0.05)。

图2 TCGA数据库中NUSAP1与总体生存率及无病间隔生存率之间的关系

3 讨论

肝癌的发生是一个复杂的过程，有多种细胞参与肝癌的发生与进展[19]。肝癌因其恶性程度高、易复发等特点，常见的治疗手段并不能显著改善肝癌患者的生存时间[20]。遗传因素作为限制肝癌患者预后的重要因素持续受到人们的关注。相比于传统的实验方法，随着生物信息学的发展以及肿瘤基因数据的共享，使得肿瘤相关基因的发现周期大为缩短，为研究肝癌发生发展机制，探寻潜在治疗靶点提供了重大便利。

NUSAP1作为有丝分裂中的一种关键蛋白，通过参与有丝分裂中纺锤体的形成得以调控细胞的增殖，其表达水平受到严格的调控，一旦NUSAP1表达异常升高，细胞

将无限增殖，甚至影响遗传物质[9]。已有研究[14-15]证实，NUSAP1在多种肿瘤(如胶质母细胞癌、肾癌、乳腺癌等)中表达均显著提高，且均与预后有关，可以作为该型肿瘤的潜在标志物。本研究通过对大量GEO芯片数据及TCGA数据集的分析，证实NUSAP1在肝癌组织中的表达明显升高。通过对单细胞测序数据的分析，使NUSAP1高表达的肝癌细胞来源于肝癌组织中的推论更为明确，进一步证实肝癌细胞的无限增殖与NUSAP1表达升高的关系。由此推断，在肝脏因遗传因素或外界理化因素刺激后，NUSAP1表达出现异常增高，使肝细胞持续增殖，最终诱导癌变的产生。通过对TCGA数据库中患者预后信息的分析发现，NUSAP1的表达水平明显影响患者的预后，高表达组的患者总体生存率、无病间隔生存率均显著低于低表达组。进一步临床资料分析证实了上述结论，并发现NUSAP1在肝癌中的表达水平与Child-Pugh分级、TNM分期、Okuda分期显著相关，表明NUSAP1的表达升高使得肿瘤细胞增殖旺盛，使肝癌进展迅速，患者病情恶化加快，从而影响预后。综上所述，NUSAP1可以作为肝癌不良预后的新标志物并可以成为肝癌治疗的潜在新靶点。

综上所述，本研究将高通量基因芯片数据、单细胞测序数据的分析与现有临床样本PCR验证相结合，避免了单一研究方式的不足，结论更有说服力。研究验证了NUSAP1在肝癌细胞中的异常表达升高且与患者预后相关，为肝癌的靶向用药提供了新的潜在靶点。

利益冲突声明：本研究不存在研究者、伦理委员会成员、受试者监护人以及与公开研究成果有关的利益冲突。

作者贡献声明：靳鑫旻、杨通旺、徐庆国负责课题设计，资料分析，撰写论文；靳鑫旻、刘欢、关鸽参与收集数据，修改论文；臧运金等负责拟定写作思路，指导撰写文章并最后定稿。

附录1见二维码

附录2见二维码