人脑膜瘤空间分辨单细胞分析鉴定影响肿瘤微环境与进展的新型细胞状态

《Nature Genetics》：Spatially resolved single-cell analyses of human meningioma identify novel cell states influencing tumor microenvironment and progression

　　

摘要：近期关于脑膜瘤分子图谱认知的进展，深化了对该疾病生物学特性及异质性的理解，并展现出临床价值。然而，仍需在单细胞分辨率下解析肿瘤内异质性，以为潜在治疗路径提供依据。本研究利用大样本队列（单核RNA测序(snRNA-seq)与高分辨率空间转录组，及批量DNA甲基化与RNA测序(n=712)、批量蛋白质组(n=88)和血浆甲基化(n=59)），系统考察了脑膜瘤的细胞类型与状态广度。研究人员证实，髓系细胞状态在不同分子分组脑膜瘤中存在差异，并从硬脑膜到肿瘤发生有意义的演化。髓系细胞状态还与独特的髓系-肿瘤（myeloid–neoplastic）互作及肿瘤基因程序相关联，提示其在塑造微环境中发挥作用。最后，多种非肿瘤细胞状态能在分子分组之外细化结局预测，提示其在脑膜瘤进展中的作用。

论文解读：《人脑膜瘤空间分辨单细胞分析鉴定影响肿瘤微环境与进展的新型细胞状态》

一、研究背景与立项依据

脑膜瘤是成人最常见的脑肿瘤，其临床异质性极为显著，而现行的世界卫生组织（WHO）分级并未完全涵盖这种异质性。此前，研究人员已建立起四种脑膜瘤分子分组（Molecular Groups, MGs）：免疫原性（MG1）、NF2野生型（MG2）、高代谢（MG3）以及增殖型（MG4），并证实细胞遗传学的改变（特别是1p缺失和22q缺失）可作为WHO分级的辅助指标。然而，肿瘤内异质性与微环境对疾病进展的塑造作用，尚需在单细胞及空间分辨率下解析。受其他脑肿瘤（如胶质母细胞瘤）中微环境重要性的启发，本研究旨在通过单核RNA测序（snRNA-seq）与空间转录组学，深度表征脑膜瘤微环境的细胞构成与状态转换，并探索其临床预后价值。该论文发表于《Nature Genetics》。

二、主要关键技术方法概览

研究人员构建了多组学、多区域的大样本队列：单核RNA测序（snRNA-seq，n=102样本，580,185个高质量细胞核）与空间转录组（Visium HD，n=43样本，7,146,538个分节细胞）为核心；配对批量DNA甲基化与RNA-seq（n=712）用于分子分组与验证；辅以批量蛋白质组（n=88）与血浆甲基化测序（n=59）。此外，设置了“内部区域”队列（同一肿瘤多区域采样）与“空间硬脑膜”队列（肿瘤及周围硬脑膜多部位采样）。生物信息学上，采用非负矩阵分解（NMF）识别重复转录程序（元程序，Metaprograms, MPs），利用CellChat与CellNEST推断细胞-细胞互作，结合贝叶斯反卷积（Bayesprism）与生存分析（Cox回归、Kaplan–Meier）关联临床结局。

三、研究结果解读

1. 脑膜瘤单细胞图谱揭示重复性转录元程序

研究人员定义“图谱”队列（n=68样本；392,818细胞核），发现大部分为肿瘤细胞（74.9%），髓系细胞主导非肿瘤区（22.3%）。通过NMF识别出肿瘤细胞11种元程序（MPs），其中MP1（细胞周期）与MP11（代谢/核糖体）跨样本共享；其余与干扰素（IFN）、肿瘤坏死因子（TNF）、上皮-间质转化（EMT）等相关，显示肿瘤间转录多样性。非肿瘤细胞中，鉴定出髓系7种MPs（如细胞周期、IFN/监视、EMT混合、TNF/炎症、组织修复、代谢等）、淋巴与内皮细胞特异性MPs，并通过公共数据正交验证。

2. 髓系细胞状态在不同现有生物学分类中差异富集

聚焦髓系细胞，通过伪 bulk与单样本基因集富集（ssGSEA）量化MPs：细胞周期MP在增殖型（MG3/4）及高WHO级别中上调；EMT MP在免疫原性（MG1）中上调；TNF MP在MG1中下调；组织修复MP在NF2野生型（MG2）中上调，且在WHO 1级更常见。细胞组成上，组织修复髓系细胞占比最高（57.6%），在MG2中尤其富集；EMT髓系在MG1富集；代谢髓系在MG2较少；细胞周期髓系在MG3/4略增。结果经批量反卷积验证，提示髓系状态与肿瘤整体生物学及临床分级紧密相关。

3. 部分髓系元程序与胶质母细胞瘤髓系程序保守

与胶质母细胞瘤（GBM）类比，细胞周期、IFN、TNF相关髓系MPs强相关，组织修复与代谢中度相关，EMT无对应GBM程序，说明部分髓系转录程序跨脑肿瘤保守。

4. 批量蛋白质组验证髓系元程序的生物学相关性

蛋白质通路与各MP的RNA特征大体对应，并筛选出候选蛋白标记：LMNB1（细胞周期）、STAT1（IFN）、NFIA（EMT）、SLC2A3（TNF）、MRC1（组织修复）、PHPT1（代谢），强调多基因特征优于单标记。

5. 部分脑膜瘤在同一肿瘤空间上存在显著肿瘤转录异质性

利用“内部区域”队列（4例肿瘤，每例3区域），通过ANOVA的F统计量评估区域转录差异。髓系细胞区域差异不显著；而某例（P3）肿瘤细胞显示较高异质性。说明肿瘤转录可空间异质，但髓系组成相对稳态。

6. 肿瘤染色体1p与22q缺失在整体脑膜瘤中大致稳定

尽管少数肿瘤转录异质，染色体22q与1p缺失比例在同一肿瘤不同区域相似（P3的1p缺失比例分别为0.90、0.86、0.87；22q为0.74、0.84、0.86），批量DNA甲基化也验证1p/22q拷贝数及甲基化风险评分在区域内稳定，支持单一临床样本通常足够；但图谱队列中偶见少量不一致病例，提示罕见克隆异质性可能。

7. 染色体22q与1p缺失伴随更大空间肿瘤转录异质性

空间转录组用Moran指数刻画局部自相关：TNF样MP（MP5）中位数最高；存在22q或1p缺失的肿瘤，多个肿瘤MP（细胞周期、IFN样、炎症、EMT、代谢）Moran指数更高，即局部转录异质性更强，或部分解释此类肿瘤预后差及耐药。

8. 空间解析采样揭示脑膜瘤免疫重编程从硬脑膜到肿瘤

“空间硬脑膜”队列（肿瘤+多部位硬脑膜）显示：肿瘤细胞比例从肿瘤、覆盖硬膜向远处（1 cm、2 cm）锐减；成纤维与内皮在远侧硬膜为主；髓系在肿瘤样本激增至42.7%，远/覆盖硬膜<5%，提示髓系招募参与发病。批次校正后，肿瘤与硬膜来源髓系细胞转录分离：硬膜髓系以TNF样程序为主，肿瘤髓系以代谢样为主；EMT与代谢MP在肿瘤髓系富集，TNF MP在肿瘤被抑制。空间转录组跨脑膜瘤-硬膜边界以100 μm分箱分析，证实TNF MP在入肿瘤前2 mm内表达渐降，提示脑膜瘤进展需抑制该髓系程序。

9. 髓系-肿瘤互作具状态特异性并塑造微环境

CellChat推断髓系（配体）-肿瘤（受体）及反向的州特异性互作（如IGF1(髓系)-IGF1R(肿瘤)关联组织修复MP；MERTK受体相关互作）。空间邻域分析（半径50 μm）发现：特定髓系MP邻近的肿瘤细胞呈现独特基因相关性；有趣的是，与髓系IFN/TNF MP相关的肿瘤基因常重叠于髓系自身MP定义基因，暗示髓系可“传染”自身转录特征给邻近肿瘤细胞。CellNEST深度学习进一步在空间上复现部分互作。髓系与肿瘤MP表达还关联局部细胞密度：TNF MP与两细胞类型密度正相关，EMT/组织修复/代谢MP弱正关联；部分肿瘤MP（EMT1、代谢）仅与肿瘤而非髓系密度相关，暗示“免疫活跃区”与“肿瘤富集区”的全局空间依赖。

10. 脑膜瘤细胞状态具预后价值，独立于分子分组与细胞丰度

批量RNA反卷积+ssGSEA：髓系细胞周期与代谢MP关联较短无进展生存（PFS），调整MG后仍为危险因子（HR分别为4.04与194.84）；EMT与组织修复MP具保护效应。按中位分割可进一步分层WHO 2/3级及血浆甲基化推断的循环髓系MP同样预后。整合7个预后MP（髓系EMT/TNF/组织修复/代谢；肿瘤炎症；淋巴细胞周期/代谢）与MG构建综合细胞状态（ICS）评分，高低风险组差异显著（调整HR=2.23），且细化MG、甲基化风险、RNA风险及放疗后复发的预测。提示微环境细胞状态可补充现有分子分型。

四、讨论与结论总结

研究人员呈现了超50万高质量细胞核的脑膜瘤单细胞分辨率图谱，结合~700万空间分割细胞的转录组与712例批量多组学，系统描绘了髓系细胞状态的空间演变、状态依赖性互作及对临床行为的影响。核心结论如下：

• •
  脑膜瘤微环境由多样且临床重要的细胞状态构成，髓系细胞尤为突出；
• •
  髓系状态依分子分组与WHO级别分布不同，并从硬脑膜向肿瘤发生转录重编程（TNF MP抑制，代谢/EMT MP增强）；
• •
  染色体1p/22q缺失在肿瘤内空间稳定，但其存在伴随更高局部肿瘤转录异质性；
• •
  髓系-肿瘤互作状态特异（如IGF1-IGF1R、MERTK相关），邻近肿瘤基因谱可“镜像”髓系MP；
• •
  多种非肿瘤细胞状态（尤其髓系EMT、组织修复、代谢、细胞周期）独立预测进展，整合为ICS可 refinement 现有分子分组与放疗反应预测；
• •
  血浆甲基化可非侵入性推断部分预后MP。

研究强调，脑膜瘤微环境是具空间结构的生态系统，其细胞状态可超越肿瘤固有分子分组影响结局，为未来靶向髓系重编程的疗法及精细风险分层提供理论与资源基础。