《Frontiers in Oncology》：Tumor–neuronal hybrid programs reveal a spatially adaptive malignant state in glioblastoma
编辑推荐：
目的：胶质母细胞瘤（GBM）表现出显著细胞可塑性，包括部分获得神经元相关转录特征。研究人员测试了公共GBM空间和单核数据集是否揭示具有空间组织且区别于经典前神经元样身份的肿瘤-神经元程序。方法：研究人员重新分析了来自Greenwald等人的公共多区域GBM空间
目的：胶质母细胞瘤（GBM）表现出显著细胞可塑性，包括部分获得神经元相关转录特征。研究人员测试了公共GBM空间和单核数据集是否揭示具有空间组织且区别于经典前神经元样身份的肿瘤-神经元程序。方法：研究人员重新分析了来自Greenwald等人的公共多区域GBM空间转录组数据，并使用来自Ravi等人的健康皮层数据作为非恶性参考。肿瘤-神经元（TN）分类通过使用置换法导出的经验阈值同时富集恶性相关和神经元相关基因程序来定义。通过自举重采样、重复阈值估计和Jaccard重叠评估分类稳定性。使用Greenwald衍生的CNA/恶性评分评估TN分类空间单元中的肿瘤相关信号。在独立的Wang恶性GBM细胞/单核数据集中进行外部验证，包括在排除原始评分程序中使用的基因后进行留出TN特异性评分。在TCGA GBM中进行了探索性整体队列分析。结果：在肿瘤区域中检测到TN分类空间单元，其频率高于健康皮层，并在空间发现标本中显示区域组织。它们表达突触和囊泡相关基因，包括SYT1和STMN2，同时保留恶性/胶质相关特征，如GFAP、VIM和CLU。成熟神经元亚型标记物富集较少，支持部分神经元模拟而非谱系转换。TN相关基因形成一个以突触、神经突相关和细胞骨架蛋白为中心的连接蛋白相互作用网络。CNA/恶性评分分析显示TN分类空间单元中的肿瘤相关信号。在Wang恶性GBM细胞验证数据集中，相对于非TN和前神经元高非TN恶性细胞/核，TN分类细胞/核保留了留出TN特异性富集。结论：公共空间和单核GBM数据集识别了一个肿瘤相关神经元模拟程序，该程序在GBM发现标本中空间组织，并在独立恶性细胞验证数据集中可检测。该程序邻近前神经元状态，但不可简化为经典前神经元身份。
论文解读文章
**研究背景与问题**
胶质母细胞瘤（GBM）是成人最常见的恶性原发性脑肿瘤，尽管采用最大安全切除、放射治疗和替莫唑胺化疗等积极多模式治疗，几乎总是复发，患者生存率极低。这种恶性行为归因于GBM细胞显著的细胞可塑性，即细胞能够在不同表型状态间转换，支持侵袭、治疗逃逸和适应局部微环境。其中一种表现是神经元模拟，即恶性细胞获得部分神经元相关特征，而不一定经历完全神经元谱系转换。近年来，肿瘤-神经元通信成为胶质瘤生物学的重要组分，神经元活动通过活性依赖性信号（如神经连接蛋白-3释放）促进胶质瘤生长，并已证实存在直接神经元-胶质瘤突触和电化学相互作用。然而，单细胞和空间转录组研究虽揭示了恶性细胞异质性，但空间转录组学保留了解剖背景，使得原位评估肿瘤-神经元程序成为可能。当前研究的问题是：是否存在一种空间组织化的肿瘤-神经元混合程序，且该程序是否可区别于经典前神经元样身份？研究人员假设，如果神经元模拟是一种适应性肿瘤相关状态，TN分类空间单元应满足：同时富集恶性和神经元程序、非随机空间分布、携带肿瘤相关CNA/恶性信号、显示可重复基因表达特征。本研究旨在通过公共空间和单核数据集验证这一假设。
**研究内容与结论**
研究人员重新分析了Greenwald等人的多区域GBM空间转录组数据（包括一个GBM标本的三个解剖区域：肿瘤核心、非核心肿瘤和浸润边缘），并使用Ravi等人的健康皮层数据作为非恶性参考。通过置换法定义经验阈值，将空间单元或细胞/核分类为肿瘤程序、神经元程序、肿瘤-神经元（TN）或其他。在发现标本中，TN分类空间单元在肿瘤区域中检测到，频率高于健康皮层，并显示区域组织。这些单元表达突触和囊泡相关基因（如SYT1、STMN2），同时保留恶性/胶质特征（如GFAP、VIM、CLU），成熟神经元亚型标记物富集较少，支持部分神经元模拟而非谱系转换。蛋白质-蛋白质相互作用（PPI）分析显示TN相关基因形成一个以突触、神经突相关和细胞骨架蛋白为中心的网络。CNA/恶性评分分析显示TN分类空间单元携带肿瘤相关信号。在Wang恶性GBM细胞/单核验证数据集中，TN分类细胞/核在排除原始评分基因后仍保留留出TN特异性富集，且相对于非TN和前神经元高非TN细胞/核更为显著。探索性TCGA GBM整体队列分析显示，TN富集评分高的患者总生存期更短（log-rank检验p=0.014459），但连续Cox模型未达显著（p=0.0633），提示潜在临床关联但需进一步验证。结论：公共空间和单核GBM数据集发现一个空间组织化的肿瘤相关神经元模拟程序，该程序邻近前神经元状态但不可简化为经典前神经元身份。
**主要关键技术方法**
本研究使用的主要关键技术方法包括：公共空间转录组数据（10x Genomics Visium平台，来自Greenwald等人的多区域GBM标本，包含肿瘤核心、非核心肿瘤和浸润边缘，以及Ravi等人的健康皮层参考）和单核/单细胞RNA-seq数据（Wang等人的恶性GBM细胞/核数据集，来自4例患者的配对肿瘤核心和瘤周脑组织，共34,023个细胞/核）。通过置换法（基于1000个随机对照基因集生成经验阈值，取第99百分位）定义恶性程序和神经元程序的显著富集，从而进行TN分类。分类稳定性通过自举重采样、重复阈值估计和Jaccard重叠评估。使用Greenwald衍生的CNA/恶性评分（基于拷贝数改变）评估空间单元肿瘤信号。外部验证中，在排除原始评分基因后计算留出TN特异性评分。探索性整体队列分析使用TCGA GBM批量转录组数据（518个样本），通过TN富集评分进行生存分析。
**研究结果**
**3.1 空间识别GBM发现标本中的肿瘤-神经元程序富集**
通过整合空间转录组学，研究人员在GBM发现标本中绘制了肿瘤、神经元和双程序富集空间单元。在肿瘤核心，大多数空间单元被分类为肿瘤程序，但TN分类单元仍被检测到；非核心肿瘤区域包含混合类别；浸润边缘显示TN单元与神经元程序和其他单元交错。健康皮层以神经元程序和其他为主，TN比例低于肿瘤区域。这支持TN双程序富集具有区域组织而非随机分布。
**3.2 TN分类空间单元表现出部分神经元模拟程序**
差异表达分析显示，TN分类单元相对于非TN单元上调的基因包括SNAP25、SYT1、STMN2和MAP2，指示神经元相关特征（囊泡处理、神经突组织和细胞骨架调控）。但成熟神经元亚型标记物（如SLC17A7、CCK）富集较少，同时保留恶性/胶质标记物（如CLU、ACTB、GFAP、VIM），支持部分神经元模拟而非完全转换。
**3.3 TN相关基因形成生物学相关互作网络**
PPI分析显示，TN上调基因形成43个节点、84条观察边的连接网络（预期8条边），平均节点度3.91，平均局部聚类系数0.407，PPI富集p值<1.0×10^-16。网络中心为突触、神经突相关和细胞骨架蛋白（如SNAP25、SYT1、NEFL、MAP2、STMN2），表明生物学一致性。
**3.4 互补类富集分析揭示额外TN特征**
TN分类观测还显示细胞外基质重塑（BGN、MCAM）、可塑性调节因子（ID1、ID3）和神经信号相关基因（APP、SCRG1）的富集，提示TN特征超越双程序阳性本身。
**3.5 TN邻近前神经元但不可还原为经典前神经元程序**
尽管TN分类观测在前神经元相关背景下最频繁，但与高前神经元非TN观测相比，TN单元仍转录上独特，上调神经元和突触基因（CNTNAP2、SCN1A、GRIA1、NRXN1、SATB2）同时保留恶性/反应性基因（GFAP、EGFR、CD44、VIM、FN1、TNC、F3、PTN、VEGFA）。前神经元评分与TN特异性评分弱相关（Spearman rho=0.051），支持TN不可简化为前神经元身份。
**3.6 CNA/恶性评分验证支持TN空间单元中的肿瘤相关信号**
在5430个肿瘤区域空间单元中，TN分类单元的中位CNA/恶性评分为0.493，高于神经元程序（0.305）和其他（0.275），低于肿瘤程序（0.646），且所有配对比较显著，提示TN单元携带中间水平的肿瘤相关信号。
**3.7 外部恶性细胞验证支持留出TN特异性富集**
在Wang数据集的34,023个恶性GBM细胞/核中，TN分类细胞/核在排除原始评分基因后仍显示更高留出TN特异性评分，相对于非TN和前神经元高非TN细胞/核，且与肿瘤程序细胞/核无显著差异，支持TN框架转移至独立数据集。
**3.8 TN富集评分的探索性整体队列分析**
TCGA GBM中，高TN评分（前20%）患者总生存期短于低评分（后20%），log-rank检验显著（p=0.014459），但连续Cox模型未达显著（HR=1.17，p=0.0633），提示探索性关联。
**讨论与结论**
讨论部分强调，本研究识别了GBM中一种肿瘤相关神经元模拟程序，该程序具有空间组织、CNA/恶性信号支持，并在独立恶性细胞数据集中验证。TN程序表现为部分获得神经元特征而非完全转换，保留恶性状态并选择性激活神经元模块。该程序邻近前神经元但不等同，可能与肿瘤-脑边界处的组织适应、侵袭和神经整合相关。探索性TCGA分析提示潜在临床关联，但需更大规模空间队列和功能实验验证。研究局限包括：主要发现基于一个GBM标本的空间单元分类，而非单细胞；评分依赖预设基因列表，但通过置换阈值和留出验证增强稳定性；CODEX数据仅用于标记可用性；TCGA分析为探索性，未调整所有协变量。结论：公共GBM空间和单核数据集识别了一个空间组织化的肿瘤-神经元模拟程序，该程序保留肿瘤身份同时选择性激活神经元程序，支持部分恶性可塑性状态，不可简化为经典前神经元身份，其与粘附、细胞外基质和界面特征关联提示在肿瘤-脑边界局部适应中的潜在作用。
