**原文公众号：一起试验网
对下述分析方法感兴趣或者没有钻研思路的小伙伴，欢送踊跃探讨！**

早！明天小编和大家剖析一篇2022年10月发表在《Frontiers in Immunology》（IF:5.0）期刊上的文章《A necroptosis related prognostic model of pancreatic cancer based on single cell sequencing analysis and transcriptome analysis》。作者应用单细胞测序剖析和转录组剖析构建了胰腺癌坏死性凋亡相干特色。本钻研可为胰腺癌的诊断、医治及预后评估提供参考。该思路同样实用于其余肿瘤的钻研，有相干需要的老师欢送分割咱们。

背景&办法
背景：
胰腺癌（Pancreatic cancer, PCa）是影响人类衰弱最具侵略性的癌症之一。约80-85%的PCa患者在诊断时呈现不可切除状态或转移。承受手术切除的PCa患者的5年生存率仅为10-25%。目前全身化疗是次要医治办法，但生存终局往往不现实。坏死性凋亡（Necroptosis）是一种新定义的程序性细胞死亡类型。一方面，坏死性凋亡是一种癌症克制因子，诱导癌细胞坏死性凋亡能够逆转其对细胞死亡的抵抗力。另一方面，坏死性凋亡被认为是癌症的促进剂。在此过程中开释的免疫活性物质和活性氧促成许多癌症通路的激活，参加免疫微环境的调节，促成肿瘤细胞的增殖和侵袭。因而，坏死性凋亡及其调控机制可能是PCa肿瘤免疫治疗的重要医治靶点。

办法：
1.GEO和TCGA公共数据库获取PCa数据；
2.WGCNA用于辨认与坏死性凋亡评分相干的模块；
3.Seurat包用于单细胞数据分析；
4.单因素和LASSO用于构建预后模型；
5.ESTIMATE用于计算免疫和基质评分；ImmuneSubtypeClassifier包用于计算
每个样本的免疫亚型；TIMER 2.0用于剖析免疫细胞浸润；TIDE数据库用于计算每个样本的TIDE评分，以评估免疫治疗的疗效；
6.MAfTools用于渐变景观剖析；
7.Rms用于构建列线图；
8.基因敲低、CCK-8测定、菌落形成剖析、迁徙和侵袭测定以及伤口愈合细胞试验用于摸索要害基因EPS8的作用；
PCR测定要害基因EPS8的表白。
钻研后果
（一）WGCNA辨认与坏死性凋亡评分相干的模块

作者先通过ssGSEA计算每个样本的坏死性凋亡评分。并依据中位数分为高坏死性凋亡组和低坏死性凋亡组。而且，高坏死性凋亡组中PCa患者的坏死性凋亡评分较高，预后较差。将“坏死性凋亡评分”作为性状，并通过WGCNA辨认了与该评分显著相干的ME绿黄模块。

（二）单细胞测序剖析

本钻研中所纳入的样本间无显著批量效应。将所有细胞汇集成14个簇，依据每个簇的遗传特色，正文到8种细胞类型。而后确定每个细胞坏死性凋亡基因的百分比。依据中值，将细胞分为高NCPTs组和低NCPTs组。进一步鉴定了两组样本之间的2518个差别表白基因，而且，这些基因次要与mRNA代谢和蛋白质定位无关。

（三）TCGA队列中预后模型的构建和评估

作者进一步将ME绿黄模块与2518个差别表白基因穿插，取得805个与预后相干的坏死性凋亡基因（Necrotic apoptosis gene）。通过单因素Cox剖析和LASSO构建预后模型。依据中位数把TCGA-PCa患者分为高NCPTs_high组和低NCPTs_low组，NCPTs_high组患者预后较差。模型预测患者预后的AUC值显著优于年龄、性别分期等临床特色。决策曲线结果显示，患者从基于NCPTS的临床干涉中获益最多。而且，多因素Cox分析表明，NCPTS是PCa患者的独立预后因素。此外，该模型能够很好地区分NCPTS_high组和NCPTS_low组患者。

（四）预后模型的内部验证和评估

为了验证预后模型，作者别离下载了PA-AU队列、PA-CA队列和GSE85916作为内部验证。验证后果均与TCGA-PCa队列中的后果保持一致。

（五）免疫浸润剖析和PCa患者对免疫治疗反馈的应答

为了说明PCa预后与免疫微环境的关系，作者进一步剖析了NCPTS_low组和NCPTS_high组之间的免疫浸润景观。结果显示，与NCPTS_high组相比，NCPTS_low组有较高的预计评分、免疫评分和基质评分，较低的肿瘤纯度。而且，NCPTS_low组中有更高水平的免疫细胞浸润，如T细胞和B细胞等；另外，NCPTS_low组中ICPs和HLA的表白也较高。

为了摸索模型中的基因是否可能领导PCa患者的分子分型，进一步进行无监督聚类分析，结果显示，与簇1相比，簇2的预后较差。随后，NCPTS分组与簇以及免疫亚型之间的关系显示，簇1次要与NCPTS_low组无关，免疫亚型C3和C6次要散布在NCPTS_low组。簇2次要与NCPTS_high组相干，免疫亚型C1和C2次要散布在NCPTS_ high组。免疫亚型之间的生存剖析显示，免疫亚型C3组的预后最好，而C3次要散布在NCPTS_low组。总之，钻研结果表明，NCPTS_low组的免疫细胞浸润程度较高，ICPs和HLA的表白也较高，免疫亚型C3次要散布在该组。

另外，NCPTS_low组的TIDE较低，表明肿瘤免疫逃逸的可能性较低，而患者免疫治疗获益的可能性较大。
（六）渐变景观剖析

鉴于基因突变在肿瘤的倒退和患者的预后中起着重要作用，进一步对NCPTS_low组和NCPTS_high组进行渐变景观剖析。两组患者最常产生渐变的前5个基因别离为KRAS、TP53、SMAD4、CDKN2A和TTN。NCPTS_high组TMB较高，而且TMB与NCPTS呈正相干（P<0.001）。TMB通常与患者预后不良无关，这可能是NCPTS_high组预后不良的起因。另外，NCPTS_high组激活了与细胞周期相干的G2M_CHECKPOINT和E2F_TARGETS通路。

（七）NCPTS模型的临床意义

为了进一步钻研NCPTS在不同队列中的预后影响，随后进行了meta剖析。结果显示，NCPTS依然是导致患者预后不良的一个因素。为了进一步促成临床胰腺癌患者的危险因素分类和治理，于是对NCPTS和临床特色进行了决策树剖析，患者被分为RIS1、RIS2、RIS3和RIS4个危险组。其中，RIS4组预后较差（P < 0.001）。列线图中TCGA-2J-AABK患者的1年、3年和5年死亡率别离为18.5%、62.5%和75.2%，列线图预测的患者预后的AUC放弃在0.8左右，显著高于其余临床特色。基于列线图的患者及时承受临床医治的获益率高于其余临床特色。

（八）在PCa细胞系中验证要害基因EPS8的作用

最初，在细胞系中验证要害基因的作用。siEPS8-3在三种siRNA中体现出最高的敲低效率，因而被选作要害基因做进一步剖析。结果显示，在CAPAN-1和PANC-1细胞系中，EPS8敲低后表白和细胞活性显著升高，细胞系产生集落的能力也显著升高，PCa细胞的迁徙和侵袭能力也大大降低。伤口愈合试验显示，EPS8组的迁徙能力在si-EPS8组中较弱。

（九）在临床PCa组织样本中验证EPS8的表白

与癌旁组织相比，要害基因EPS8在PCa肿瘤组织中显著上调。