整合多组学识别PANX2为透明细胞肾细胞癌的铁死亡抑制因子和潜在治疗靶点

《Frontiers in Immunology》：Integrative multi-omics identifies PANX2 as a ferroptosis suppressor and potential therapeutic target in clear cell renal cell carcinoma

　　

背景：透明细胞肾细胞癌（ccRCC）在代谢上具有铁死亡（ferroptosis）易感性，然而铁死亡相关基因的预后意义及机制贡献仍不完全明确。本研究旨在识别具有预后价值的铁死亡相关生物标志物，并阐明其在ccRCC进展中的功能相关性。方法：研究人员整合了单细胞RNA测序（scRNA-seq）、批量RNA测序（bulk RNA-seq）、空间转录组学（spatial transcriptomics）和基于机器学习的特征选择，以识别ccRCC中与铁死亡相关的预后基因。构建了一个四基因风险模型和一个整合列线图（nomogram），并在独立队列中进行了评估。PANX2被优先用于实验验证，使用了稳定敲低模型、RNA测序、脂质过氧化和铁离子探针、氧化还原测定、蛋白质印迹（Western blot）、荧光素酶报告基因实验、异种移植瘤模型以及免疫健全小鼠肾癌模型。结果：识别出一个四基因预后特征（CA9、PVT1、RRM2和PANX2），并用于构建风险模型，在训练集和外部验证集中均具有一致的预测性能。在这些基因中，PANX2主要富集于上皮肿瘤区域，且此前在ccRCC中未被功能表征。敲低PANX2抑制了ccRCC细胞增殖，降低了抗氧化能力，增加了细胞内Fe²⁺积累和脂质过氧化，并使细胞对erastin诱导的铁死亡更敏感。机制上，PANX2缺失与Akt/mTOR通路活性降低及SLC7A11表达减少相关；使用Akt激活剂或过表达SLC7A11可减轻铁死亡相关表型。在免疫健全小鼠肾癌模型中，PANX2敲低伴随着CD45⁺白细胞、CD3⁺T细胞和CD8⁺T细胞浸润增加，支持PANX2依赖性铁死亡抵抗与肿瘤免疫微环境之间的潜在关联。结论：本研究将PANX2鉴定为ccRCC中一种铁死亡抑制因子，并支持PANX2-Akt/mTOR-SLC7A11相关信号轴参与氧化还原稳态和肿瘤进展。该铁死亡相关预后模型和列线图可辅助风险分层，而PANX2代表一个候选治疗脆弱点，有待进一步的机制和转化验证。

透明细胞肾细胞癌（ccRCC）是肾细胞癌最常见的组织学亚型，尽管手术、酪氨酸激酶抑制剂和免疫检查点抑制剂等治疗手段改善了临床管理，但复发和治疗耐药仍是主要挑战，尤其在晚期疾病中。当前预后评估仍主要依赖临床病理参数，而近期研究强调非侵入性成像和炎症指标可作为补充工具，但仍有必要建立整合临床风险分层与生物学可解释分子机制的预后框架。铁死亡（ferroptosis）是一种铁依赖性脂质过氧化物积累驱动的调节性细胞死亡（RCD）方式，区别于凋亡等其他RCD形式。由于恶性细胞通常维持较高的铁池和活性氧水平以支持代谢需求，诱导铁死亡可能成为针对常规化疗或分子靶向药物耐药癌症的有前景策略。然而，决定ccRCC预后、治疗反应和免疫背景的铁死亡相关基因尚不充分明确。单细胞RNA测序（scRNA-seq）能够高分辨率解析异质性肿瘤生态系统，空间转录组学保留组织架构并支持肿瘤和基质区室的空间解读，结合批量转录组队列和透明的特征选择策略，有助于识别既统计可重复又生物学可解释的预后生物标志物。本研究整合了scRNA-seq、批量RNA-seq、空间转录组学和机器学习，以识别ccRCC中铁死亡相关预后基因，构建了四基因预后特征和列线图，并优先对PANX2进行实验验证，因为其在ccRCC中的作用此前尚不清楚。研究人员证明PANX2缺失破坏氧化还原稳态，增加脂质过氧化和细胞内铁积累，并使ccRCC细胞对铁死亡敏感，进一步提供了PANX2与Akt/mTOR-SLC7A11信号通路相关并可能影响肿瘤免疫微环境的证据。该研究发表在《Frontiers in Immunology》。

主要关键技术方法包括：从癌症基因组图谱（TCGA）和国际癌症基因组联盟（ICGC）等公共数据库获取ccRCC转录组数据；使用Seurat软件包处理scRNA-seq数据集（GSE159115，包含8个肿瘤样本和6个正常样本），通过Harmony算法校正批次效应，并利用DoubletFinder去除双细胞；利用空间转录组数据集（GSE250163，10x Visium平台）进行细胞类型分解（RCTD）和空间聚类；采用LASSO回归、SVM-RFE、梯度提升机（GBM）和随机森林四种机器学习算法筛选铁死亡相关预后基因，最终通过多变量Cox回归确定四基因特征；构建风险评分模型和列线图，并在多个外部验证队列中评估其性能。实验验证方面，使用慢病毒介导的shRNA在ccRCC细胞系（786-O和Caki-1）中建立PANX2稳定敲低模型，进行RNA测序、脂质过氧化探针（BODIPY 581/591 C11、Liperfluo）和铁探针（FerroOrange）检测、氧化还原指标（NADPH/NADP⁺、GSH/GSSG比值）测定、蛋白质印迹、荧光素酶报告基因实验、以及异种移植瘤模型（裸鼠）和免疫健全小鼠（BALB/c）肾癌模型（Renca细胞）。

研究结果部分：
**Identification of ferroptosis-related genes in ccRCC patients**：通过差异表达分析获取scRNA-seq和TCGA-KIRC队列中肿瘤与正常组织之间的差异表达基因（DEGs），与GeneCards来源的铁死亡相关基因集取交集，得到23个候选基因。经四种机器学习算法筛选及多变量Cox回归，最终确定四个铁死亡相关预后基因：CA9、PVT1、RRM2和PANX2。
**Development of a prognostic signature of ccRCC**：基于四基因构建风险评分模型，计算公式为风险评分 = (?0.190 × CA9表达) + (0.415 × PVT1表达) + (0.242 × RRM2表达) + (0.193 × PANX2表达)。在训练队列中，模型C-index为0.68，3年和5年AUC分别为0.683和0.724；在外部ICGC队列中，3年和5年AUC分别为0.699和0.754，其他两个验证队列也表现良好。高评分组患者预后显著更差，且与晚期肿瘤分期正相关。
**Construction of a nomogram**：整合风险评分与年龄、肿瘤分期构建列线图，其预测准确性优于单个指标，校准曲线显示预测与观察生存率一致，决策曲线分析（DCA）证实其临床效用。
**Gene functional enrichment analysis of DEGs**：GSEA分析显示高危组富集磷脂外排通路，低危组富集脂肪酸氧化、脂质分解、铁硫簇组装等代谢通路，提示脂质代谢和氧化还原状态的差异。
**Ferroptosis-related prognostic signature is associated with distinct immune landscapes**：ESTIMATE和ssGSEA分析显示高危组免疫评分和基质评分较高，低危组富集NK细胞、中央记忆CD8⁺T细胞、未成熟树突状细胞和NKT细胞。四基因与髓系来源抑制细胞（MDSC）浸润正相关。在免疫健全小鼠模型中，PANX2敲低导致CD45⁺白细胞、CD3⁺T细胞和CD8⁺T细胞浸润增加。
**The prognostic signature correlates with immunotherapy response**：TIDE分析显示高危组T细胞功能障碍和排斥评分更高，多个免疫检查点基因（如CTLA4、LAG3）表达更高，且对阿西替尼敏感性降低。
**Single-cell and spatial transcriptomic characterization of ferroptosis-related features**：scRNA-seq分析表明PANX2主要表达于上皮肿瘤细胞，空间转录组学进一步证实PANX2信号主要分布于上皮肿瘤区域。
**Establishment of PANX2-knockdown models and transcriptomic profiling**：RNA-seq显示PANX2敲低后与氧化还原稳态和铁死亡相关的通路显著富集。
**PANX2 knockdown impairs antioxidant capacity and elevates lipid peroxidation in ccRCC cells**：PANX2敲低增加ccRCC细胞对erastin诱导的细胞死亡敏感性，且被铁死亡抑制剂ferrostatin-1特异性挽救，而非凋亡、坏死性凋亡或焦亡抑制剂。NADPH/NADP⁺和GSH/GSSG比值降低，脂质过氧化和Fe²⁺水平升高，回补PANX2可逆转这些表型。
**PANX2 regulates ferroptosis resistance via Akt/mTOR signaling and promotes tumor growth**：PANX2敲低导致SLC7A11表达下调，Akt/mTOR磷酸化降低，荧光素酶报告基因实验显示PANX2过表达增强SLC7A11启动子活性。CCK-8和EdU实验表明PANX2敲低抑制细胞增殖，异种移植瘤模型中肿瘤生长减缓，4-HNE染色增强。Akt激活剂SC79可恢复SLC7A11表达并减轻铁死亡，而SLC7A11过表达同样能逆转脂质过氧化和Fe²⁺积累。临床组织芯片（116例ccRCC标本）显示PANX2与SLC7A11在肿瘤组织中协同上调并正相关，高PANX2表达与晚期TNM分期（III–IV）、男性、包膜侵犯和更高Fuhrman核分级显著相关。

讨论部分总结：本研究通过整合scRNA-seq、批量RNA-seq、空间转录组学和机器学习，建立了铁死亡相关预后模型，并优先验证了PANX2在ccRCC中的铁死亡抑制功能。四基因特征捕捉了预后和生物学异质性，但应主要视为统计优化模型，而非单一线性调控通路。高危组磷脂外排富集，低危组脂肪酸氧化和铁硫簇组装富集，反映代谢差异。PANX2通过维持Akt/mTOR-SLC7A11轴促进铁死亡抵抗，但PANX2作为大孔通道蛋白对SLC7A11的调控可能通过上游信号间接实现，具体中间环节尚需研究。免疫分析提示PANX2介导的铁死亡抵抗可能影响肿瘤免疫微环境，但需在ccRCC特异性免疫治疗队列中验证。研究存在局限性，包括队列规模有限、基因间直接共调控未建立、免疫微环境实验使用鼠模型而非人类ccRCC等。结论部分翻译：总之，本研究整合了单细胞、批量、空间和机器学习分析，识别了ccRCC中铁死亡相关预后特征，并优先对PANX2进行功能验证。研究人员表明，PANX2通过维持氧化还原稳态和调节Akt/mTOR-SLC7A11相关通路，支持ccRCC细胞增殖和铁死亡抵抗。这些发现细化了PANX2在ccRCC特异性背景下的作用，提供了一个生物学可解释的预后框架，并将PANX2提名为进一步临床前和临床研究的潜在治疗脆弱点。