依据铁死亡相关长链非编码RNA特征对结肠癌患者进行预后分析

摘要

目的探讨应用基于铁死亡相关长链非编码RNA(lncRNA)构建的预后模型预测结肠癌患者预后的价值。方法从FerrDb数据库下载铁死亡相关基因, 从癌症基因组图谱(TCGA)数据库下载建库至2021年11月结肠癌患者RNA测序基因数据和临床数据。应用R3.6.3软件对从TCGA数据库获取的结肠癌基因表达数据与FerrDb数据库铁死亡相关基因进行分析, 获取结肠癌与正常组织差异表达的铁死亡相关基因;应用R3.6.3软件计算结肠癌铁死亡相关基因和lncRNA之间表达的相关性, 确定结肠癌铁死亡相关lncRNA。筛选出与生存相关的差异表达铁死亡相关lncRNA, 将其纳入多因素Cox比例风险模型, 构建结肠癌预后模型;根据lncRNA表达量, 通过预后模型计算结肠癌患者预后的风险评分。根据中位风险评分, 将从TCGA数据库收集的临床病例分为高风险组和低风险组, 各223例, 对两组进行Kaplan-Meier生存分析。应用受试者工作特征(ROC)曲线分析应用预后模型风险评分及临床特征预测所有患者生存的效果。应用GSEA 4.1.0软件对高风险组和低风险组lncRNA进行基因集富集分析(GSEA), 采用R3.6.3软件的ggpubr程序包对高风险组和低风险组间差异表达lncRNA进行免疫细胞及免疫功能的单样本基因集富集分析(ssGSEA)。结果依据数据库获取的铁死亡相关基因和差异表达基因的交集获得了65个差异表达的铁死亡相关基因, 分析最终纳入24个结肠癌预后相关lncRNA, 基于这些lncRNA构建了预后模型。Kaplan-Meier生存分析显示低风险组的生存优于高风险组(P<0.001);ROC曲线分析显示, 应用预后模型风险评分预测所有患者1、2、3年生存的曲线下面积(AUC)均>0.750, 风险评分预测所有患者1年生存的AUC大于年龄、性别、美国综合癌症网络(NCCN)临床分期、T分期、N分期、M分期。GSEA结果显示, 高风险组和低风险组差异表达lncRNA集中在肿瘤和免疫相关通路;ssGSEA结果显示, 高风险组和低风险组间T细胞、巨噬细胞、肥大细胞、中性粒细胞、免疫刺激、人白细胞抗原(HLA)及Ⅰ型、Ⅱ型干扰素反应方面均存在差异(均P<0.05), 其中免疫检查位点的CD200、TNFRSF14表达量差异明显(均P<0.01)。结论铁死亡相关lncRNA可能在结肠癌的肿瘤免疫中发挥重要作用, 其可用于结肠癌患者的预后分析。

关键词