B2M通过焦亡调控ELANE影响胰腺癌进展

《Translational Oncology》：B2M regulates ELANE in pyroptosis to affect pancreatic cancer progression

　　

焦亡在胰腺癌中的作用仍存在争议。研究人员利用双样本孟德尔随机化（MR）整合FinnGen数据库（314,193例对照，731例病例）的全基因组关联研究（GWAS）数据、冰岛人群的蛋白质数量性状位点（pQTL）数据以及英国生物银行数据，系统研究了焦亡基因与胰腺癌之间的因果关联。研究发现，Beta-2-微球蛋白（B2M）通过上调中性粒细胞弹性蛋白酶（ELANE）间接增加胰腺癌风险——据研究人员所知，这是首次在胰腺癌背景下建立B2M与ELANE之间基于中介分析的因果遗传关联。中介分析显示，ELANE解释了该效应的20.572% [15.32%–25.81%]。敏感性分析确认了结果的稳健性，未发现显著的多效性，生物信息学验证支持了研究人员的MR结果。药物敏感性分析进一步识别出潜在的治疗药物。研究结果支持B2M作为胰腺癌的诊断生物标志物，因其在肿瘤中显著过表达且具有高诊断准确性（AUC = 0.976，95% CI：0.958–0.993），并突出了通过数据驱动的MR中介框架识别的B2M–ELANE轴作为有前景的治疗靶点。

**论文解读：B2M通过焦亡调控ELANE影响胰腺癌进展**

**研究背景与问题**
胰腺癌是恶性程度极高的消化系统肿瘤，具有高发病率、高复发转移率、高死亡率以及低早期诊断率、低有效治疗率、低5年生存率的“三高三低”特征。尽管手术和辅助化疗能使少数患者获益，但超过80%的患者因早期症状不典型而错过最佳治疗窗口，确诊时已出现远处转移或局部浸润。现有生物标志物如CA19-9，受良性病变影响且无法反映肿瘤微环境中的关键分子事件，在监测疗效和预测耐药方面存在局限。焦亡作为一种依赖gasdermin蛋白家族（GSDMs）的程序性细胞死亡，在肿瘤发生发展中具有双重作用：既可触发抗肿瘤免疫，也可能通过持续炎症促进肿瘤微环境形成。在胰腺癌这种免疫抑制且纤维化的基质中，焦亡的作用尤为复杂。已有线索表明，持续的KRAS-MAPK通路激活抑制caspase-4活性从而拮抗焦亡，吉西他滨耐药细胞系出现GSDME表达下调及焦亡抵抗表型，但尚无系统研究建立焦亡相关基因与胰腺癌风险之间的因果联系，或阐明关键焦亡效应蛋白的上游遗传调控。因此，研究人员假设关键焦亡执行蛋白的遗传决定表达与胰腺癌风险存在因果关联，且该关联由上游调控蛋白介导，并采用生物信息学和孟德尔随机化方法进行探索。该研究发表在《Translational Oncology》。

**主要关键技术方法**
研究人员采用双样本孟德尔随机化（MR）框架，整合来自FinnGen数据库（314,193例对照，731例胰腺导管腺癌病例）的GWAS数据、冰岛人群的pQTL数据（来自deCODE genetics）以及英国生物银行制药蛋白质组学项目（UKB-PPP）的pQTL数据。通过严格的工具变量筛选策略：选取目标基因转录起始位点（TSS）上下游1 Mb内的顺式pQTL（P < 5×10^-8，连锁不平衡修剪r2 < 0.1），获得1615个独立SNPs作为工具变量。采用逆方差加权法（IVW）作为主要分析方法，辅以MR-Egger回归、加权中位数法、MR-PRESSO等方法进行敏感性分析。中介分析采用三步法线性回归模型，以B2M为暴露、ELANE为中介、胰腺癌为结局，通过1000次bootstrap估计间接效应及其置信区间。外部验证使用UCSC XENA数据库的TCGA胰腺癌样本（n=178）和GTEx正常胰腺组织（n=171）的RNA-seq数据，以及GEO单细胞数据集GSE111672，进行差异表达、ROC曲线、生存分析、免疫浸润分析和药物敏感性分析。

**研究结果**

**MR分析焦亡与胰腺癌的关系**
采用IVW法评估焦亡相关基因对胰腺癌的因果效应，共使用199个SNPs作为工具变量。结果显示，中性粒细胞弹性蛋白酶（ELANE）和白细胞介素1β（IL1B）与胰腺癌显著相关。IL1B呈负相关（潜在保护作用），ELANE呈正相关（潜在促癌作用）。

**敏感性分析**
Cochran’s Q检验显示各SNPs间无异质性（P > 0.05），MR-Egger回归截距接近0（P > 0.05），MR-PRESSO检验无异常SNPs（P > 0.05），留一法分析未发现单个SNP显著改变因果估计，表明结果稳健。

**中介分析**
通过两步MR分析发现，B2M基因表达显著正向调控ELANE表达（路径a β=0.133，95% CI: 0.041–0.225），ELANE表达显著增加胰腺癌风险（路径b β=0.345，95% CI: 0.112–0.578），B2M对胰腺癌的总效应为0.223（95% CI: 0.151–0.295），其中通过ELANE介导的间接效应为0.046（95% CI: 0.015–0.077），中介比例为20.572%（95% CI: 15.32%–25.81%）。

**外部验证**
基于UCSC XENA数据库，Mann-Whitney U检验显示B2M在胰腺癌组织中显著高表达（P < 0.001）。ROC曲线分析显示B2M诊断胰腺癌的AUC为0.976（95% CI: 0.958–0.993）。Kaplan-Meier生存分析和多变量Cox回归（校正年龄、性别、肿瘤分期）显示B2M表达与总生存期（OS，HR=1.22，P=0.340）、无进展生存期（PFI，HR=1.13，P=0.546）和疾病特异性生存期（DSS，HR=1.11，P=0.660）无显著独立关联。ssGSEA免疫浸润分析显示B2M与多种免疫细胞亚型显著正相关，如活化树突状细胞（aDC，R=0.661，P<0.001）和Th1细胞（R=0.529，P<0.001）。单细胞分析（GSE111672）显示B2M在免疫细胞（CD8T细胞、单核/巨噬细胞）、恶性细胞和成纤维细胞中高表达，在腺泡细胞和导管细胞中低表达。药物敏感性分析（CellMiner数据库）发现B2M表达与磺胺替尼（Sulfatinib，Cor=0.402，P=0.001）和卡莫司汀（Carmustine，Cor=0.349，P=0.006）的敏感性最强正相关，与比沙可啶（Bisacodyl，Cor=-0.280，P=0.030）呈负相关。

**讨论与结论**
讨论部分指出，本研究的核心创新在于通过多层面遗传推断首次建立B2M与ELANE在胰腺癌中的因果中介遗传关联，填补了焦亡-胰腺癌研究中因果遗传证据的空白。B2M作为HLA-I类分子轻链，可能通过以下机制调控ELANE：①HLA-I类信号在肿瘤浸润中性粒细胞中通过NF-κB激活诱导ELANE转录；②B2M通过EGFR-VEGF通路上调促炎细胞因子（IL-1β、TNF-α），从而诱导ELANE表达；③B2M重塑肿瘤免疫微环境，通过成纤维细胞分泌HDACs和aDC分泌IL-6/STAT3通路间接增强ELANE表达。ELANE作为中性粒细胞弹性蛋白酶，可直接切割GSDMD的C268位点，释放具有膜穿孔活性的N端片段，触发焦亡。研究还讨论了代谢重编程与焦亡的关联，以及B2M-ELANE轴在胰腺癌中的生物学合理性。

研究结论部分翻译如下：
在本研究中，通过创新性整合孟德尔随机化框架和中介效应分析，研究人员系统识别出B2M是唯一通过孟德尔随机化和中介分析因果影响胰腺癌风险的关键焦亡执行蛋白ELANE的上游调控因子。研究结果表明，B2M基因通过调控ELANE的表达水平精确影响焦亡通路的激活强度，从而在胰腺癌的发生发展中发挥重要作用。这一发现为深入解析胰腺癌的发病机制开辟了新视角，不仅加深了研究人员对肿瘤表观遗传调控复杂性的科学认识，也凸显了焦亡在胰腺癌表观遗传调控机制中的关键作用。研究结果强有力地强调了表观遗传调控网络在肿瘤恶性转化中的核心作用，为基于表观遗传干预的胰腺癌精准治疗策略的开发奠定了坚实的理论基础。