转移性结直肠癌免疫微环境及预后相关的免疫相关差异表达基因分析

作     者：刘艳红 丁杰 葛美玲 高玒 王小巍 胡月 Liu Yanhong;Ding Jie;Ge Meiling;Gao Hong;Wang Xiaowei;Hu Yue

作者机构：南京鼓楼医院生物样本库南京210008 

出 版 物：《肿瘤研究与临床》 (Cancer Research and Clinic)

年 卷 期：2022年第34卷第9期

页      面：641-647页

学科分类：1002[医学-临床医学] 100214[医学-肿瘤学] 10[医学] 

基　　金：江苏省科技资源(重大疾病生物样本)统筹服务平台南京鼓楼医院生物样本库(BM2015004) 江苏省重大疾病生物资源样本库开放课题(SBK202006002、SBK202006003)。 

主　　题：结直肠肿瘤 肿瘤微环境 预后 免疫细胞 

摘      要：目的通过生物信息学方法, 分析转移性结直肠癌原发肿瘤组织及转移灶组织的免疫微环境差异, 筛选与预后相关的转移性结直肠癌特异性免疫相关差异表达基因。方法从基因表达综合(GEO)数据库中下载结直肠癌及转移灶基因芯片数据集GSE131418, 包括莫菲特癌症中心MCC队列的样本数据517例和Consortium队列样本数据618例;从免疫学数据库和分析门户IMMPORT中下载免疫相关基因集, 包含免疫相关的基因2 483个。从癌症基因组图谱(TCGA)数据库中下载结直肠癌组织及癌旁组织的RNA测序数据共695例, 临床信息627例。采用ESTIMATE算法计算转移灶组织和原发肿瘤组织的基质细胞评分、免疫细胞评分、基质免疫总评分, 并用CIBERSORT反卷积算法对结直肠癌原发肿瘤组织及转移灶组织的22种免疫细胞浸润情况进行比较分析。筛选免疫相关差异表达基因并进行京都基因和基因组百科全书(KEGG)信号通路富集分析。根据各免疫相关差异表达基因表达水平的中位数将患者分为高、低表达组, 分别使用Kaplan-Meier法和Cox比例风险模型对免疫相关差异表达基因进行分析, 筛选两种方法结果中与预后显著相关的基因(均P0.01), 用Cox回归进行多因素分析。比较TCGA数据库和GEO数据库GSE131418数据集中各基因在肿瘤组织和癌旁组织、原发肿瘤组织和转移灶组织中的表达差异, 并进行生存分析。结果转移性结直肠癌转移灶组织基质细胞评分、免疫细胞评分、基质免疫总评分均低于原发肿瘤组织(均P0.001)。与原发肿瘤组织相比, 转移灶组织中活化的自然杀伤(NK)细胞、单核细胞、CD8+ T细胞、T细胞、活化的树突细胞比例增高, 而未活化的肥大细胞、未活化的树突细胞、未活化的NK细胞、活化的记忆CD4+ T细胞、M1型巨噬细胞以及中性粒细胞比例降低。转移性结直肠癌转移灶组织与原发肿瘤组织免疫相关差异表达基因289个, 其中上调基因101个, 下调基因188个。KEGG信号通路富集结果显示, 在转移性结直肠癌转移灶组织的免疫微环境中, 发生了血管内皮生长因子(VEGF)信号通路、程序性死亡受体配体1(PD-L1)表达与程序性死亡受体1(PD-1)检查点通路、辅助性T细胞(Th)1和Th2细胞分化、Th17细胞分化、NF-κB信号通路、白细胞介素17(IL-17)信号通路、趋化因子信号通路、T细胞受体信号通路、MAPK信号通路、自然杀伤细胞介导的细胞毒性等信号通路的富集。筛选出7个预后相关的免疫相关差异表达基因, 分别为ANGPTL5、FPR1、HSPA8、NR2E3、PSMD2、PSMD8、SBDS。Cox回归多因素分析结果显示, 转移性结直肠癌免疫相关差异表达基因ANGPTL5(HR=2.69, 95%CI 1.22~5.92, P0.05)、HSPA8(HR=0.57, 95%CI 0.33~0.97, P0.05)、SBDS(HR=2.23, 95%CI 1.18~4.21, P0.05)是转移性结直肠癌独立的预后影响因素。ANGPTL5在肿瘤组织中表达低于正常组织, 在转移灶组织中的表达高于原发肿瘤组织, 而肿瘤组织高表达ANGPTL5的患者预后更差。HSPA8在肿瘤组织中表达高于正常组织, 在转移灶组织中的表达低于原发肿瘤组织, 而肿瘤组织高表达HSPA8的患者预后更好。SBDS在肿瘤组织中表达低于正常组织, 在转移灶组织中的表达低于原发肿瘤组织, 而肿瘤组织高表达SBDS的患者预后更差。结论转移性结直肠癌的免疫微环境与原发肿瘤相比存在较大差异, 其免疫细胞浸润程度降低, 整体处于免疫抑制状态, 筛选出预后相关的转移性结直肠癌特异性免疫相关差异表达基因可能是新的转移性结直肠癌治疗靶点。