《International Journal of Molecular Sciences》：Integrated Transcriptomic and Single-Cell Analyses Identify HILPDA as a Hypoxia-Mediated Regulator of Ferroptotic Signaling in Glioblastoma
Nelin Hacioglu
编辑推荐：
本文针对胶质母细胞瘤（GBM）缺氧微环境与铁死亡（一种铁依赖性脂质过氧化相关的细胞死亡过程）信号之间联系不清的问题，整合转录组学与单细胞RNA测序分析，揭示了缺氧诱导脂滴相关蛋白（HILPDA）是缺氧-铁死亡-免疫微环境轴的关键关联基因。研究发现，HILPDA高表达与铁死亡相关基因特征、铁-ROS轴分子（如HMOX1、NOX4、STEAP3）及促炎免疫微环境显著相关，并可能通过缺氧通路间接介导铁死亡信号。该研究为理解GBM中缺氧驱动的代谢适应与治疗抵抗提供了新的转录组学关联框架。
胶质母细胞瘤（Glioblastoma, GBM）是成人中最常见且最具侵袭性的原发性恶性脑肿瘤，其治疗一直充满挑战。标准疗法虽包含手术、放疗和化疗，但患者预后依然很差，这主要归因于肿瘤的高度异质性、对治疗的抵抗以及一个关键特征——广泛存在的缺氧微环境。这个氧气不足的区域就像是肿瘤内部的“黑暗森林”，它不仅滋养了最具攻击性的肿瘤细胞，还重塑了肿瘤的代谢与免疫状态，使其能躲避治疗的“追杀”。近年来，一种名为铁死亡（ferroptosis）的新型细胞死亡方式为癌症治疗带来了新希望。铁死亡不同于细胞凋亡，它依赖于铁离子和致命的脂质过氧化。理论上，诱导铁死亡或许是消灭那些对传统疗法“刀枪不入”的肿瘤细胞的“秘密武器”。然而，在GBM的复杂“战场”上，铁死亡究竟是“盟友”还是“内鬼”？它与缺氧信号这个强大的“环境塑造者”之间，存在何种隐秘的关联？这正是研究者们渴望破解的谜题。因为只有理清缺氧如何调控铁死亡及相关基因表达网络，才有可能找到GBM的新治疗靶点，打破其治疗抵抗的僵局。
为此，研究人员开展了一项整合性研究。他们从公共数据库（如TCGA）获取了GBM和低级别胶质瘤（LGG）的大量RNA测序数据，并利用单细胞RNA测序数据进行验证。核心的生物信息学分析方法包括：差异表达基因分析、基因集富集分析（GSEA/GSVA）、生存分析、相关性与中介分析，以及单细胞转录组数据的降维与细胞聚类分析。这些技术手段帮助研究者在不依赖新实验样本的情况下，系统地探索缺氧信号、铁死亡相关基因特征和免疫微环境之间的转录组学关联。
1. 鉴定胶质瘤中与缺氧相关的铁死亡基因特征
研究者首先将缺氧条件下差异表达的基因与已知的铁死亡相关基因集进行比较，发现了29个位于二者交集的核心候选基因，这提示缺氧信号与铁死亡相关基因表达程序在转录组层面存在重叠。热图分析显示，在根据缺氧评分划分的高危组中，这些特征基因的表达呈现一致性的上调。这些基因的表达水平在不同级别胶质瘤（GBM与LGG）间存在显著差异，且多个基因的高表达与患者较差的总生存期相关，表明这个缺氧相关的铁死亡基因特征可能作为胶质瘤进展的转录组标志物。
2. 铁死亡活性在胶质瘤中的预后意义
通过基因集变异分析计算每个患者的铁死亡评分（FerroScore），研究者发现，在TCGA-GBM队列中，FerroScore高的患者总生存期显著更短，但在LGG队列中无显著差异。这表明铁死亡相关基因表达模式可能与高级别胶质瘤的不良预后相关。多变量Cox回归分析显示，年龄、性别、IDH野生型状态和未甲基化的MGMT启动子状态是独立的预后因素，而HILPDA表达本身并非独立的预后因子。对FerroScore高表型相关基因的通路富集分析发现，它们显著富集于炎症和应激反应通路，如细胞因子-细胞因子受体相互作用、TNF信号通路、HIF-1信号通路等，这支持了铁死亡相关转录程序与缺氧及炎症微环境存在关联的观点。
3. HILPDA与铁死亡和缺氧特征的关联
在核心候选基因中，缺氧诱导脂滴相关蛋白（HILPDA）因其与缺氧反应直接相关且在脂质代谢和应激适应中的作用而被选为重点候选基因。分析发现，HILPDA表达与FerroScore呈显著正相关。功能富集分析表明，与HILPDA正相关的基因显著富集于炎症和免疫应答程序，而与HILPDA负相关的基因则富集于神经元信号相关过程。这暗示HILPDA高表达的胶质瘤细胞在转录层面可能向分化程度更低、代谢重编程和炎症表型转变。中介分析进一步揭示，HILPDA与铁死亡相关基因特征之间的关联，可能是通过缺氧相关的转录程序介导的。
4. HILPDA与铁/ROS通路及缺氧相关铁死亡特征的关联
更深入的分析显示，HILPDA高表达的样本中，与铁稳态、活性氧生成和缺氧反应相关的基因（如HMOX1、STEAP3、NOX4）也同步上调。相关性分析证实HILPDA与这些铁-ROS轴分子表达呈强正相关。功能富集评分也显示，HILPDA高表达组的铁死亡和缺氧评分均显著增加。另一项中介分析再次支持了HILPDA→缺氧→铁死亡的间接关联模型，表明HILPDA可能通过影响缺氧相关的转录活动，进而与铁死亡相关基因表达模式产生关联。
5. HILPDA与免疫微环境的关联
对肿瘤免疫微环境的分析发现，HILPDA高表达组具有更高的CD8+T细胞评分、细胞溶解活性、IFN-γ反应、自然杀伤细胞和巨噬细胞特征评分，表明存在免疫细胞浸润增加。然而，HILPDA表达也与更高的PD-L1（CD274）表达和T细胞耗竭评分相关。这表明HILPDA高表型与一个“发炎但功能失调”的肿瘤免疫微环境相关，即虽然存在免疫细胞浸润和炎症反应，但可能伴随免疫检查点分子的上调和T细胞功能耗竭。
6. HILPDA-缺氧-铁死亡关系的单细胞转录组学分析
单细胞RNA测序分析在细胞水平验证了这一关联。HILPDA表达富集于特定的肿瘤/干细胞样细胞亚群，并与缺氧评分区域显著重叠。在细胞水平，HILPDA与缺氧评分显著相关，但与铁死亡评分的直接相关性很弱。结构方程模型（SEM）分析证实，HILPDA对铁死亡评分的直接影响不显著，但其通过缺氧评分产生的间接效应显著。这进一步支持了在单细胞水平，HILPDA与铁死亡相关基因特征的关系也是由缺氧相关的转录活动介导的。
综合以上研究结果，本文的结论是：本研究通过整合大量和单细胞转录组学分析，系统地揭示了胶质母细胞瘤中缺氧信号、铁死亡相关基因特征和免疫微环境重塑之间的转录组学关联。研究发现，铁死亡相关基因特征，尤其在GBM中，与不良预后相关，并可能与缺氧驱动的代谢和炎症重编程有关。缺氧诱导脂滴相关蛋白（HILPDA）被鉴定为一个关键的缺氧相关基因，其表达与铁死亡相关基因特征、铁-ROS轴分子（如HMOX1、NOX4、STEAP3）以及一个促炎但功能失调的免疫微环境显著相关。重要的是，中介分析和单细胞数据均表明，HILPDA与铁死亡相关基因表达模式之间的关联，很可能是通过缺氧相关的转录程序间接介导的，而非HILPDA直接调控铁死亡。因此，HILPDA应被视为一个反映肿瘤缺氧微环境压力的生物标志物，而非独立的临床预后因子。
这项研究的核心意义在于，它提供了一个整合的转录组学框架，将缺氧信号、铁死亡相关基因表达和免疫微环境重塑联系起来，深化了我们对GBM这种致命肿瘤内部复杂相互作用网络的理解。它强调了在GBM的背景下，铁死亡相关基因特征可能更多地与肿瘤的适应性和侵袭性表型相关，而非单纯的细胞死亡易感性。研究将HILPDA定位为该调控轴中的一个关键关联节点，为未来探索缺氧环境下脂质代谢如何影响铁死亡和免疫应答的功能性研究奠定了基础，有望为开发针对GBM缺氧微环境的新治疗策略提供新的思路和潜在靶点。本研究已发表在《International Journal of Molecular Sciences》期刊上。
