综合PANoptosis分析鉴定免疫抑制亚型及具有MLKL功能验证的胶质母细胞瘤预后标志

《Annals of Clinical and Translational Neurology》：Integrated PANoptosis Profiling Identifies Immunosuppressive Subtypes and a Prognostic Signature With Functional Validation of MLKL in Glioblastoma

　　

目的：胶质母细胞瘤（GBM）的预后仍然非常不利，主要由于高度的肿瘤异质性和免疫抑制微环境。然而，PANoptosis在此背景下的功能作用尚不清楚。方法：研究人员通过K-means聚类对患者进行分层。利用Cox回归确定的预后相关基因构建风险评分模型，并在独立队列中进行验证。使用CIBERSORT和ESTIMATE分析免疫浸润。单细胞RNA测序（scRNA-seq）描绘了肿瘤微环境。在MLKL敲低后进行功能实验。结果：基于PANoptosis相关基因鉴定出两种分子亚型，具有不同的生存和免疫特征。一个五基因（MLKL、YWHAG、GZMB、ELANE、CASP4）风险评分作为独立的预后因素。高风险组表现出炎症性但功能失调的肿瘤免疫微环境，以更高的PD-L1表达、T细胞功能障碍和Merck18评分为特征。scRNA-seq证实了GBM中PANoptosis活性升高。最后，MLKL敲低被证实可抑制恶性表型并诱导凋亡。结论：研究人员的发现将PANoptosis与GBM异质性联系起来，提供了一个预后模型，并提名MLKL作为关键功能介质，这可能有助于患者分层和靶向治疗的发展。

胶质母细胞瘤（GBM）是中枢神经系统高度侵袭性肿瘤，预后极差，主要归因于显著的肿瘤异质性和免疫抑制微环境。现有治疗手段（手术、放疗、化疗）效果有限，亟需解析新机制并发现新靶点。程序性细胞死亡（PCD）在组织稳态和异常细胞清除中至关重要，而PANoptosis作为一种整合焦亡、凋亡和坏死性凋亡的炎症性PCD，在肿瘤免疫调控中具有潜力，但其在GBM中的功能角色尚不明确。基于此，研究人员利用癌症基因组图谱（TCGA）、中国脑胶质瘤基因组图谱（CGGA）和基因表达综合数据库（GEO），对PANoptosis相关基因（PANRGs）进行生物信息学分析，构建并验证风险评分模型，并通过体外实验验证核心基因MLKL（混合谱系激酶结构域样蛋白）的功能。该研究发表于《Annals of Clinical and Translational Neurology》，建立了PANoptosis与GBM异质性的关联，提供了预后模型并提名MLKL为关键功能介质，为患者分层和靶向治疗策略开发奠定基础。

研究人员采用的主要技术方法包括：从TCGA（169例GBM样本）、CGGA（176例GBM样本）和GEO（GSE43378，50例样本）获取转录组和临床数据；通过K-means共识聚类进行分子亚型鉴定；利用Cox回归和LASSO回归（最小绝对收缩和选择算子）构建并验证五基因风险评分模型；使用CIBERSORT和ESTIMATE算法分析免疫浸润景观；应用单细胞RNA测序（scRNA-seq，GSE162631数据集，含3例GBM和3例正常脑组织）解析肿瘤微环境（TME）细胞组成；通过细胞功能实验（CCK-8、划痕愈合、克隆形成、流式细胞术）及分子检测（RT-qPCR、Western blot、ELISA）评估MLKL敲低后的恶性表型和炎症因子变化。

**3.1 基于PANRGs表达谱的分子亚型鉴定及其免疫微环境特征**
通过共识聚类将GBM患者分为两个亚型，Kaplan-Meier分析显示Cluster 2总生存期（OS）显著短于Cluster 1；CIBERSORT分析揭示两亚型免疫细胞浸润（如巨噬细胞、单核细胞）存在显著差异，差异表达基因（DEGs）富集于免疫和炎症通路。

**3.2 风险评分模型的构建与验证**
单变量Cox回归从66个PANRGs中筛选出5个预后相关基因（MLKL、YWHAG、GZMB、ELANE、CASP4），经LASSO-Cox回归确认后构建风险评分模型。在TCGA训练集中，ROC曲线显示模型对1、3、5年生存具有良好预测效能；高风险组（HRS）OS显著短于低风险组（LRS），并在CGGA和GSE43378独立队列中得到验证。多变量Cox回归确认风险评分为独立预后因素，独立于年龄、性别和IDH突变状态。

**3.3 临床列线图构建与基因组特征分析**
整合风险评分、替莫唑胺（TMZ）治疗和IDH突变状态的列线图可预测1、2、3年生存，校准曲线显示良好一致性。亚组分析显示，在接受放疗或TMZ治疗的患者中，HRS组OS仍显著更差。肿瘤突变负荷分析显示HRS组以PTEN（43%）、TP53（34%）、TTN（24%）突变为主，LRS组以TP53（36%）、EGFR（31%）、TTN（30%）为主。

**3.4 风险评分与TIME的相关性分析**
HRS组ESTIMATE、免疫、基质评分显著高于LRS组；CIBERSORT分析显示两组在活化记忆CD4⁺ T细胞和活化自然杀伤（NK）细胞比例上存在差异。风险评分与活化记忆CD4⁺ T细胞、活化肥大细胞正相关，与静息肥大细胞和活化NK细胞负相关。TIDE分析表明HRS组具有更高的癌症相关成纤维细胞（CAFs）浸润、CD8⁺ T细胞浸润、CD274（PD-L1）表达和IFN-γ信号，但T细胞功能障碍评分更高；LRS组则具有更高的髓系来源抑制细胞（MDSCs）浸润。Merck18评分在HRS组显著升高，提示更强的预测免疫治疗应答潜力。

**3.5 风险分层的功能富集与治疗反应预测**
免疫表型评分（IPS）分析显示，仅在CTLA4(?)/PD-1(?)亚组中HRS组评分显著低于LRS组，提示该类“冷”肿瘤可能从抗PD-1/CTLA4单药治疗中获益有限。基因集富集分析（GSEA）显示HRS组富集免疫炎症通路（如“细胞因子–细胞因子受体相互作用”、“IL-17信号”），并下调“细胞周期”。药物敏感性分析鉴定出HRS组对KU-5933和SB2-16763更敏感，LRS组对daporinad和pyridostatin更敏感；YWHAG和MLKL表达与多种药物（如crizotinib、TMZ）的IC50负相关，而GZMB、ELANE、CASP4表达与crizotinib和sorafenib的IC50正相关。

**3.6 单细胞转录组分析揭示风险模型基因表达与TME细胞景观**
scRNA-seq分析鉴定出10种主要细胞类型，GO富集显示GBM中小胶质细胞、巨噬细胞和T细胞中炎症免疫通路显著富集。细胞间通讯分析显示GBM中相互作用数量和强度增加，形成以髓系细胞为中心的信号网络。PANoptosis基因集在GBM中活性显著高于正常组织；五个风险基因在GBM的巨噬细胞、T细胞和内皮细胞中表达更高且分布更广。

**3.7 关键基因MLKL的临床意义与表达验证**
Kaplan-Meier分析显示MLKL、CASP4、GZMB、YWHAG高表达与不良预后显著相关。GEPIA数据库证实MLKL在GBM组织中表达高于正常脑组织，CCLE数据库显示其在多种GBM细胞系中表达。基于实验室可获取的细胞系，选择高表达MLKL的U343和U118MG细胞进行后续功能实验。

**3.8 MLKL敲低对GBM细胞恶性表型的影响**
设计三个特异性siRNA，筛选出敲低效率最高的si-MLKL#3，Western blot确认其有效降低MLKL及其磷酸化形式（p-MLKL）蛋白水平。功能实验显示，MLKL敲低显著抑制细胞活力（CCK-8）、减慢划痕愈合、减少克隆形成。

**3.9 MLKL敲低影响炎症因子表达与细胞凋亡**
基于MLKL表达的GSEA显示高表达组富集炎症通路。在U343和U118MG细胞中，MLKL敲低显著抑制IL-1α mRNA但上调HMGB1 mRNA，ELISA证实HMGB1分泌增加，Western blot显示细胞内IL-1α蛋白水平下降。流式细胞术表明MLKL敲低显著增加两种细胞系的早期凋亡。

讨论部分总结：本研究构建的PANoptosis相关分子分类和预后风险模型揭示了其在GBM发病机制中的核心作用，并与TME异质性和治疗反应相关。两个亚型具有不同生存和免疫特征，PANRGs不仅调节细胞死亡，还塑造免疫微环境。五基因风险评分模型在多个独立队列中验证为独立预后因素，优于传统因素。HRS组呈现“免疫富集但功能失调”表型，表现为高免疫/基质评分、CD8⁺ T细胞浸润、PD-L1表达和IFN-γ信号，但T细胞功能障碍评分升高，提示单药免疫检查点抑制剂获益有限，联合治疗可能更有效。scRNA-seq验证了PANoptosis在GBM中活性升高，风险基因在髓系、T细胞和内皮细胞中广泛分布，髓系细胞作为信号中枢。MLKL敲低抑制GBM细胞恶性表型并诱导凋亡，同时差异调控IL-1α和HMGB1等炎症因子，可能解释高风险组“高炎症但免疫抑制”表型，靶向MLKL或可重塑TME以增强免疫治疗。

研究结论翻译：本研究阐明了PANRGs在GBM预后、免疫微环境异质性和肿瘤细胞恶性中的关键作用。构建的五基因风险评分模型可作为可靠且独立的预后指标。机制研究确定MLKL为关键功能基因，通过调控细胞增殖、细胞死亡和炎症因子分泌驱动GBM进展。这些发现加深了对GBM中细胞死亡与免疫相互作用的理解，为预后评估提供了有价值工具，并为靶向PANoptosis以增强免疫治疗反应的新型联合疗法奠定了理论基础。