胃癌免疫预后模型

今天给大家带来一个老套路文章,文章分析方法其实都是大家已经司空见惯的方法,但是它为啥能发到6分多,让我们一起去探个究竟~

论文题目:TP53相关胃癌免疫预后模型的构建和验证

研究背景

目前已经有研究表明,TP53突变可预测胃癌患者整体生存和无复发生存期。虽然TP53突变影响胃癌患者预后已得到证实,但其对STAD患者免疫功能的具体影响尚未得到充分研究。

方法及结果

➀ TP53在胃癌样本中的突变情况

从TCGA数据库中下载胃癌(STAD)突变数据、临床数据、基因表达数据。对STAD的基因的突变频率进行统计,结果发现TP53排名第二。

➁ TP53突变对免疫功能影响

将STAD样本分为TP53突变组和TP53野生组。两组样本的基因表达数据进行基因富集分析(GSEA)。结果发现有7个免疫相关的通路富集在TP53野生组,TP53突变组没有免疫相关通路富集。

➂ 建立免疫相关预后模型

上述GSEA分析中7个免疫相关通路共包含766个基因。单因素Cox回归分析,结果显示有45个基因与预后相关。LASSO回归分析、多因素Cox回归分析、最终筛选出9个基因。用这9个基因建立预后模型,公式:Risk score = ECAMK4*0.456 + ECYP19A1*0.823 + EZFP36*0.195 + ECLCF1 *0.226 + ECACTIN *(−0.653) + ENRARP*(−0.170) + ECFHR4*1.050 + EAMBP* 0.181 + EFBN1*0.149。模型测试结果(TCGA数据测试,图ABC)和验证结果(GEO数据验证,图DEF)如下图所示。该预后模型对1、3、5年的整体生存的预测AUC值均能达到0.5以上,说明该模型的预测效果良好。

➃ 根据TP53突变状态进行预后特征的生存分析

TCGA数据库中STAD样本的总体生存和TP53的图片没有显著相关性(图A)。为了验证预后模型是否独立于TP53突变,作者将样本STAD分为了两种亚型:TP53突变型和野生型,发现预后模型能够很好的预测这两个亚型的患者的预后(图BC)。由于不同类型的TP53突变通常会导致不同的生物学后果,作者分析了TP53各类突变与STAD患者预后之间的关系,然而,结果表明TP53突变并不影响STAD患者的预后。作者将TP53突变频率最高的类型(missense mutation subgroup)的样本提出来,发现预后模型对该亚型的预后也有很好的预测效果(图E)。

➄ STAD患者的免疫状况

为了探究预后模型的免疫特性,作者应用ESTIMATE计算STAD样本的基质评分、肿瘤纯度、免疫评分和ESTIMATE评分。将预后模型计算的风险分数和上述评分进行相关性分析,结果表结果发现风险分数与基质评分、免疫评分和ESTIMATE评分呈正相关,与肿瘤纯度呈负相关(图ABCD)。此外,作者还研究了低、高风险STAD患者基质评分、肿瘤纯度、免疫评分和ESTIMATE评分的差异。高风险患者的免疫、基质和ESTIMATE评分均高于低风险患者,肿瘤纯度低于低风险患者(P < 0.05,图EFGH)。

➅ 患者的浸润免疫细胞情况分析

肿瘤浸润免疫细胞比例和水平的变化可能代表不同个体的内在特征。CIBERSORT估算高风险组和低风险组患者的22种免疫细胞组成(图A)。22个肿瘤免疫细胞呈弱或中度相关(图B)。高风险组和低风险组患者的免疫细胞组成进行比较,发现高风险组STAD患者的单核细胞、M0巨噬细胞和滤泡辅助T细胞比例显著高于低风险组,静息树突状细胞和CD4记忆激活T细胞比例显著低于低风险组(P < 0.05,图C)。

➆ 预后模型与免疫检查点关系分析

免疫检查点药物已被证明具有逆转肿瘤免疫抑制作用的抗肿瘤作用,因此作者分析了风险分数与不同免疫检查点蛋白(PD-L1, PD-1, CTLA-4, TIGIT, LAG-3 and TIM-3)表达水平的相关性。结果发现风险分数与与TIM3和BTLA显著相关(图AB)。近一步的分析表明高风险患者TIM3和BTLA的表达明显高于低风险患者(P < 0.05)(图6C和D)。以上结果提示高风险患者不良预后可能与免疫抑制微环境有关。

➇ STAD富集通路分析

以logFC > 1和FDR < 0.05为阈值选择低风险和高风险患者差异表达的免疫基因。筛选与风险分数有显著相关性的的基因(Pearson相关系数> 0.2和P < 0.05)。一共筛选到39个基因(图E)。对l这39个基因进行GO和KEGG分析,结果发现STAD中与风险分数相关的基因主要富集于免疫效应过程的调控、体液免疫反应的调控、体液免疫反应的调控、ECM受体相互作用等(图6F和G)。

⑨ 基于预后模型的列线图构建

单因素、多因素Cox回归分析,发现预后模型可以作为独立的预后因素。nomogram生存预测模型具有良好的预后。

相信大家看到这里也能明白本文为啥能发6分多了。和别的文章不同的是,本文特征基因的筛选、预后模型的建立只是文章的一小部分,文章大量的分析都在探究预后模型与免疫之间的关系,分析细致、紧贴主题,是文章的一大亮点。另外由一个单基因引发的一系列研究,也是文中的一大亮点,老师们如果有想要研究的单基因,也可以尝试这种思路~

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,189评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,577评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,857评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,703评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,705评论 5 366
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,620评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,995评论 3 396
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,656评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,898评论 1 298
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,639评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,720评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,395评论 4 319
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,982评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,953评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,195评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 44,907评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,472评论 2 342

推荐阅读更多精彩内容