单细胞分析鉴定出一个与髓系OSM信号和CEBPD活性相关的VIM+胶质瘤干细胞样血管-免疫状态

《Frontiers in Immunology》：Single-cell analysis identifies a VIM+ glioma stem-like vascular–immune state linked to myeloid OSM signaling and CEBPD activity

　　

背景：胶质瘤表现出肿瘤间和肿瘤内异质性，异柠檬酸脱氢酶（IDH）状态与生物学行为和临床侵袭性相关。证据表明，胶质瘤进展依赖于胶质瘤干细胞样肿瘤细胞（GSTCs）与肿瘤微环境之间的相互串扰，其中免疫和血管重塑促进可塑性、侵袭和免疫逃逸。然而，在不同IDH背景下的GSTC状态架构及其与血管-免疫信号的关联尚不清楚。研究人员旨在定义IDH相关的GSTC状态及其功能和调控特征。方法：使用Seurat整合、聚类、差异表达分析和inferCNV辅助恶性细胞鉴定，分析了来自33名胶质瘤患者（GSE278456）的单细胞RNA测序数据。重新聚类的GSTCs通过基因本体论（GO）/Hallmark富集、基于曲线下面积的细胞层面基因集富集分析（AUCell）、CytoTRACE、Monocle、Slingshot、CellChat和pySCENIC进行表征，以描述状态、发育轨迹、细胞-细胞通讯和调控网络。在U251和LN229细胞中沉默CEBPD后，进行了定量实时PCR（qRT-PCR）、细胞计数试剂盒-8（CCK-8）、集落形成、Transwell迁移和凋亡实验。结果：鉴定了6个细胞群和5个GSTC亚群。一个VIM⁺ GSTC状态（C0），在IDH野生型肿瘤中富集，高表达CHI3L1、VIM、LGALS3、VEGFA和SEC61G。与其他状态相比，C0表现出更强的血管生成、缺氧、IL6/JAK/STAT3信号和上皮-间充质转化相关活性，并富集氧化磷酸化、黏着斑、细胞-底物连接和肌动蛋白/钙黏蛋白结合程序。轨迹分析将C0置于一个早期到中间节点，连接不同的后期状态，表明其可塑性不同于纯增殖状态。CellChat提示髓系细胞到C0的显著通讯，以OSM为候选输入，LIFR/IL6ST为受体基础。pySCENIC鉴定了一个以CEBPD、RUNX2、POU2F3、TEAD3和TEAD4为中心的C0偏好调控子景观，以及一个与应激适应和细胞内稳态相关的M2调控子模块。CEBPD敲低减少了增殖、集落形成和迁移，同时增加了凋亡。结论：研究人员鉴定了一个IDH野生型富集的VIM⁺胶质瘤干细胞样状态，具有血管生成、缺氧、间充质样和炎症特征，并提名髓系来源的OSM信号和CEBPD作为这种血管-免疫相关恶性表型的候选介质。这些发现提供了一个基于状态的假设生成视角，用于理解不同IDH背景下胶质瘤的进展，并突出了候选生物标志物和用于验证及治疗分层的机制切入点。

**论文解读：单细胞分析揭示VIM⁺胶质瘤干细胞样血管-免疫状态及其与髓系OSM信号和CEBPD活性的关联**

**研究背景与问题**
胶质瘤，尤其是高级别胶质瘤和胶质母细胞瘤（GBM），是成人中枢神经系统中最具侵袭性的原发性恶性肿瘤。当前标准治疗（手术、放疗和替莫唑胺化疗）疗效有限，主要障碍包括肿瘤广泛浸润、血脑屏障限制、显著的瘤内异质性和高复发率。异质性不仅存在于肿瘤细胞本身，还体现在以髓系细胞富集和淋巴细胞浸润受限为特征的免疫抑制微环境，以及肿瘤细胞与周围微环境（如神经元、胶质细胞）的持续复杂相互作用中。同时，GBM常伴有明显的血管增生，提示恶性细胞状态改变与微环境重塑相互交织，共同驱动肿瘤进展。在分子分层中，IDH状态与肿瘤生物学行为密切相关：IDH野生型（IDH-WT）肿瘤的瘤内异质性更高，与更侵袭性表型相关。然而，不同IDH背景下胶质瘤干细胞样肿瘤细胞（GSTCs）的内部状态组成、状态转换位置及其与髓系微环境的关键通讯关系尚未系统阐明。

**研究目的与意义**
研究人员旨在通过分析公开的单细胞转录组数据，描绘不同IDH分子背景下胶质瘤的细胞组成和GSTCs的内部状态异质性，特别关注在IDH-WT肿瘤中相对富集的GSTC状态。通过整合功能分析、发育轨迹推断、细胞-细胞通讯分析和转录调控网络分析，研究人员表征了该状态相关的间充质样、缺氧适应、炎症和血管生成程序，并评估其与髓系来源信号的潜在联系。最后，通过体外实验对单细胞分析提名的候选调控因子CEBPD进行功能评估。该研究发表在《Frontiers in Immunology》上，为理解不同IDH背景下胶质瘤恶性进展提供了基于状态的框架，并突出了与微环境相关恶性进展的潜在调控节点。

**关键技术方法**
本研究主要基于公开数据集GSE278456（包含33例胶质瘤患者的单细胞RNA测序数据），采用Seurat整合与聚类、inferCNV辅助恶性细胞鉴定、AUCell评估干性基因集活性、CytoTRACE、Monocle和Slingshot进行发育轨迹推断、CellChat推断细胞-细胞通讯（基于配体-受体共表达），以及pySCENIC重建基因调控网络。样本队列来源为Gene Expression Omnibus。此外，在U251和LN229胶质瘤细胞系中进行CEBPD沉默（siRNA），并通过qRT-PCR、CCK-8、集落形成、Transwell迁移和流式细胞凋亡实验检测表型变化。

**研究结果**
1. **跨IDH分子背景胶质瘤单细胞图谱构建及GSTCs状态异质性揭示**
通过整合分析，研究人员鉴定出6个主要细胞群（髓系细胞、GSTCs、少突胶质细胞样干细胞肿瘤细胞OSTCs、成纤维细胞、T/NK细胞、内皮细胞EC）和5个GSTC亚群（C0 VIM⁺ GSTCs、C1 DLL3⁺ GSTCs、C2 XIST⁺ GSTCs、C3 CENPF⁺ GSTCs、C4 DCX⁺ GSTCs）。C0在IDH-WT中相对富集，高表达CHI3L1、VIM、LGALS3、VEGFA和SEC61G。

2. **功能分析显示C0 VIM⁺ GSTCs具有间充质样和微环境适应复合特征**
Hallmark通路分析表明，C0在血管生成、上皮-间充质转化（EMT）、IL6/JAK/STAT3信号和缺氧通路中活性显著高于其他亚群。GO富集分析显示C0在生物学过程层面富集氧化磷酸化和细胞呼吸，在分子功能层面富集钙黏蛋白结合和肌动蛋白结合，在细胞组分层面富集黏着斑和细胞-底物连接。

3. **发育轨迹和CellChat通讯推断提示C0位于状态转换节点并可能接受髓系来源OSM输入**
CytoTRACE显示C0具有较高发育可塑性。Monocle和Slingshot轨迹分析表明C0位于早期到中间节点，连接分叉的后期状态（C1和C2）。CellChat分析显示髓系细胞与C0之间通讯显著，其中OSM作为候选配体，C0表达受体LIFR/IL6ST，提示髓系来源OSM信号可能为C0的外部输入。

4. **pySCENIC分析揭示CEBPD相关调控子特征和C0状态的M2模块偏好**
pySCENIC鉴定出C0偏好性调控子（RUNX2、POU2F3、CEBPD、TEAD3、TEAD4）。调控子模块分析将M2模块强烈富集于C0，该模块与泛素介导的蛋白水解、内质网蛋白加工、自噬和线粒体自噬等应激适应通路相关。

5. **CEBPD敲低抑制胶质瘤细胞生长和迁移并增加凋亡**
在U251和LN229细胞中，CEBPD沉默后集落形成能力、CCK-8增殖活性和Transwell迁移能力均降低，而凋亡比例增加，支持CEBPD参与恶性表型维持。

**讨论与结论**
本研究系统描绘了不同IDH背景下胶质瘤的细胞组成，解析了GSTCs的内部状态异质性，并通过体外实验对候选调控因子进行了功能评估。C0 VIM⁺ GSTCs代表一个相对IDH-WT富集的恶性状态，整合了间充质样、缺氧适应、血管生成和炎症特征，将其与血管-免疫微环境背景联系起来，并提名CEBPD作为候选调控节点。研究结论指出：研究人员鉴定了一个IDH野生型富集的VIM⁺胶质瘤干细胞样状态，具有血管生成、缺氧、间充质样和炎症特征，并提名髓系来源的OSM信号和CEBPD作为这种血管-免疫相关恶性表型的候选介质。这些发现提供了一个基于状态的假设生成视角，用于理解不同IDH背景下胶质瘤的进展，并突出了候选生物标志物和用于验证及治疗分层的机制切入点。需要指出，OSM-STAT3-CEBPD框架目前基于计算推断和常规细胞系实验，仍需空间原位证据和患者来源模型验证。