整合性批量分析和单细胞分析揭示ASL在胶质母细胞瘤中MIF相关髓系程序中的作用

《Functional & Integrative Genomics》：Integrative bulk and single-cell analyses implicate ASL in MIF-associated myeloid programs in glioblastoma

　　

精氨酸代谢与胶质母细胞瘤（GBM）的进展和免疫调控相关，但具有明确细胞背景的预后生物标志物仍未被充分表征。研究人员分析了TCGA GBM中18个精氨酸代谢相关基因，采用LASSO Cox回归结合十折交叉验证及后续单变量/多变量Cox模型筛选预后基因，依据最优截断值对患者进行分层并进行Kaplan-Meier分析，在CGGA和GEO数据集中验证，通过ROC曲线和决策曲线分析（DCA）评估预测性能。对5例GBM样本的单细胞RNA测序（scRNA-seq）数据使用Seurat处理，在单细胞分辨率下量化精氨酸代谢活性，对髓系细胞进行亚群细分，通过伪时间分析、RNA velocity和CellChat分析状态动态。ASL、FAH和NAGS被确定为候选预后基因，其中ASL在TCGA队列中被表征为与高危分子亚型（如IDH野生型和1p/19q非共缺失）明确相关的预后指标。ASL高表达在TCGA中预测较短总生存期，并在CGGA和GEO数据集中得到验证，ASL在DCA中显示出最佳临床预后区分度和净获益。在scRNA-seq中，38,263个高质量细胞用于下游分析，ASL表达和精氨酸代谢活性在单核细胞/巨噬细胞群体中富集。轨迹分析将MKI67⁺单核细胞置于上游，并显示ASL沿伪时间递减，RNA velocity支持定向转变。CellChat凸显了髓系网络中巨噬细胞迁移抑制因子（MIF）、MHC-Ⅱ和SPP1相关信号通路，与ASL和MIF轴基因之间的Spearman相关模式一致。ASL被确定为精氨酸代谢相关预后信号，尤其与GBM高危分子亚群相关，并与肿瘤微环境中以髓系为中心的免疫程序和MIF相关通讯特征显示出强相关性。

**论文解读文章**

**研究背景与问题**
胶质母细胞瘤（GBM）是成人中最常见且最具侵袭性的原发性恶性脑肿瘤，尽管采用最大安全手术切除联合替莫唑胺同步放化疗（肿瘤电场治疗为近年补充），其中位总生存期仍仅约15个月。GBM的高复发率和治疗耐药性主要源于其高度侵袭性及肿瘤微环境（TME）的动态重塑。TME中，肿瘤相关髓系细胞（尤其是单核细胞和单核细胞来源的巨噬细胞）占肿瘤内活细胞的50%，它们常被极化为抗炎表型，促进肿瘤进展、侵袭和耐药。然而，针对这些细胞表型可塑性的治疗策略仍处于探索阶段。精氨酸代谢在肿瘤生物学中发挥关键调控作用，2019年一项针对精氨酸琥珀酸合成酶1（ASS1）缺陷型复发GBM患者的Ⅰ期临床试验显示，聚乙二醇化精氨酸脱亚氨酶（PEG-ADI）治疗安全性良好，80%患者达到疾病稳定。精氨酸琥珀酸裂解酶（ASL）作为精氨酸生物合成中的关键酶，催化精氨酸琥珀酸裂解为精氨酸和延胡索酸，其在乳腺癌、肝细胞癌、结直肠癌及GBM中均被报道与癌细胞稳态维持相关。然而，ASL在GBM中的预后价值及其与免疫微环境的关系尚不明确，且缺乏基于单细胞转录组学的细胞特异性和动态分析。为此，研究人员整合批量RNA测序和单细胞RNA测序数据，系统探究精氨酸代谢相关基因在GBM中的预后意义及髓系免疫程序中的作用。

**研究内容与结论**
研究人员通过分析MSigDB数据库中的18个精氨酸代谢相关基因，在TCGA GBM队列中采用LASSO Cox回归、单变量及多变量Cox模型筛选出ASL、FAH和NAGS三个候选预后基因，其中ASL在多变量分析中保留独立预后意义（HR=1.904，95%CI: 1.241–2.922，p=0.0032）。Kaplan-Meier生存分析显示，ASL高表达患者总生存期显著缩短，并在CGGA和三个独立GEO数据集（GSE4412、GSE13041、GSE7696）中得到验证。ROC曲线和决策曲线分析（DCA）表明ASL具有最佳临床预后区分度和净获益。ASL表达与PRS分类（原发性vs复发性）、IDH突变状态（野生型vs突变型）及1p/19q共缺失状态显著相关，但在多变量分析中并非独立于这些分子标志物，提示其预后价值主要源于富集于高危分子亚群（IDH野生型、1p/19q非共缺失）。Gene Ontology（GO）及KEGG富集分析显示精氨酸代谢相关通路显著富集，基因集富集分析（GSEA）进一步揭示ASL高表达组富集干扰素γ应答、炎症反应、同种异体移植排斥、TNF-α信号通路（经NF-κB）及干扰素α应答等免疫相关通路。
单细胞RNA测序分析（5例GBM患者，共38,263个高质量细胞）鉴定出5个主要细胞群体：T淋巴细胞、星形胶质细胞、内皮细胞、单核细胞/巨噬细胞及组织干细胞。ASL表达及精氨酸代谢活性评分（AUCell）在单核细胞/巨噬细胞群体中显著富集。对髓系细胞亚群分析后，鉴定出12个亚群，单核细胞中ASL表达显著高于巨噬细胞。轨迹分析（Monocle3和Slingshot）将MKI67⁺单核细胞确定为发育起点，ASL表达在伪时间早期（0-5）逐渐升高后下降，RNA velocity支持该定向分化过程。细胞通讯分析（CellChat）揭示髓系网络中巨噬细胞迁移抑制因子（MIF）、MHC-Ⅱ及SPP1相关信号通路显著富集，且ASL与MIF轴基因（MIF、CD44、CD74、CXCR4）呈显著正相关，与MKI67呈负相关（Spearman相关性）。
该研究将ASL鉴定为GBM中精氨酸代谢相关预后信号，尤其与高危分子亚群相关，并揭示其与髓系中心免疫程序及MIF相关通讯特征的强关联。论文发表在《Functional》。

**主要关键技术与方法**
研究人员利用TCGA（美国癌症基因组图谱）、CGGA（中国脑胶质瘤基因组图谱）及GEO（GSE4412、GSE13041、GSE7696）的批量RNA测序数据，结合5例GBM患者（GSE271379）的单细胞RNA测序数据。关键方法包括：LASSO Cox回归（十折交叉验证）用于筛选预后基因；单变量/多变量Cox比例风险模型确定独立预后因素；通过最大化选择秩统计确定最优截断值进行Kaplan-Meier生存分析；时间依赖性接收者操作特征曲线（ROC）和决策曲线分析（DCA）评估预测性能；使用Seurat包处理单细胞RNA测序数据，AUCell算法量化精氨酸代谢活性；Monocle3和Slingshot进行伪时间轨迹分析；scVelo进行RNA velocity分析推断细胞状态转换方向；CellChat包推断细胞间配体-受体介导的通讯网络。

**研究结果**
**1. 鉴定具有预后意义的精氨酸相关基因**
通过LASSO Cox回归和单变量/多变量Cox分析，从18个精氨酸代谢相关基因中筛选出ASL、FAH和NAGS，其中仅ASL在多变量分析中保持独立显著性（HR=1.904，p=0.0032），被确定为后续核心基因。

**2. 验证GBM预后标志物**
在TCGA队列中，ASL高表达组总生存期显著缩短（p=0.00026），并在CGGA（p=0.048）及三个GEO数据集（GSE4412、GSE13041、GSE7696）中均得到验证（p<0.05）。ROC分析显示ASL在1年、2年、3年总生存期预测中AUC分别为0.643、0.592、0.578；DCA显示ASL在所有阈值概率下均提供最佳净获益。

**3. 整合精氨酸代谢特征与已确立的分子预后因素**
ASL表达与PRS分类（p.adj=1.94e-06）、IDH突变状态（p.adj=1.20e-11）及1p/19q共缺失状态（p.adj=1.65e-06）显著相关，但与MGMT启动子甲基化无显著关联（p.adj=0.532）。多变量分析提示ASL的预后价值主要源于富集于IDH野生型和1p/19q非共缺失等高危分子亚群。

**4. 精氨酸相关基因富集分析**
GO分析显示生物过程富集于精氨酸代谢过程、氨基酸代谢过程；细胞组分富集于线粒体基质、膜筏；分子功能富集于有机酸结合、羧酸结合。KEGG通路富集于精氨酸生物合成、氨基酸生物合成、精氨酸和脯氨酸代谢。GSEA揭示ASL高表达组显著富集干扰素γ应答（NES=3.32）、炎症反应（NES=3.22）、TNF-α信号经NF-κB（NES=3.18）等免疫相关通路。

**5. 描绘GBM肿瘤微环境异质性**
单细胞转录组分析鉴定出5个主要细胞类型，ASL表达和精氨酸代谢活性评分在单核细胞/巨噬细胞群体中显著富集，提示其细胞类型特异性作用。

**6. 髓系细胞代表肿瘤微环境代谢异质性**
对髓系细胞亚群分析后，鉴定出12个亚群，包括单核细胞和巨噬细胞相关状态。单核细胞中ASL表达显著高于巨噬细胞，精氨酸代谢活性评分在髓系群体中广泛富集。

**7. 单核细胞免疫动态中的精氨酸代谢活性与互作**
伪时间分析将MKI67⁺单核细胞置于发育起点，ASL表达沿伪时间早期升高（约0-5）后下降，RNA velocity支持定向分化。Slingshot分析验证了该轨迹，并确认MKI67⁺单核细胞为起点。

**8. 多方法轨迹分析鉴定MKI67⁺单核细胞为GBM相关单核细胞的前体群体**
Monocle3和Slingshot一致性表明MKI67⁺单核细胞为分化起点，ASL表达在该群体中较高，随分化进程递减。

**9. 细胞通讯分析提示髓系群体中潜在的MIF相关信号网络**
CellChat分析显示CD44-CD74、CD44-CXCR4、SPP1-CD44等配体-受体对在MKI67⁺和S100A12⁺单核细胞间显著富集。MIF、MHC-Ⅱ及SPP1通路在髓系网络中占主导，且ASL与MIF轴基因（MIF、CD44、CD74、CXCR4）呈显著正相关（Spearman r=0.4-0.54），与MKI67负相关（r=-0.37）。

**讨论与结论**
讨论部分重点总结了现有证据：代谢重编程是恶性肿瘤的标志，精氨酸代谢在GBM中高度活跃，靶向精氨酸可用性可能成为治疗突破口。研究人员通过整合分析，将ASL确立为GBM中与精氨酸代谢相关的预后信号，尤其与高危分子亚群（IDH野生型、1p/19q非共缺失）密切相关。单细胞证据表明ASL表达特征化早期MKI67⁺单核细胞群体，并与其分化路径一致。ASL与MIF轴成分的相关性提示精氨酸代谢与髓系中心通讯特征之间的潜在关联。这些发现推动了以精氨酸靶向药物（如ADI-PEG20）为基础的个体化治疗策略的发展。研究局限性包括：scRNA-seq样本量小、缺乏蛋白水平验证和功能实验、未考虑空间复杂性。未来需结合更大队列、空间转录组学和体内/外模型验证ASL的因果作用。

**结论**
通过综合单细胞转录组谱和批量RNA测序分析，本研究系统描绘了GBM肿瘤微环境的复杂异质性，揭示了ASL在GBM发病机制和患者预后中的先前未被认识的作用。关键发现是，ASL作为与高危分子亚型（包括IDH野生型和1p/19q非共缺失）密切相关的预后指标，其与MIF信号轴呈强正相关，该轴与MKI67⁺单核细胞群体的分化程序紧密相连。这种ASL相关的单核细胞分化程序与GBM进展密切相关。这些分子见解不仅增进了对GBM生物学的理解，还揭示了ASL可作为针对高危分子谱患者的有前景的治疗靶点，可能为GBM的靶向干预策略开辟新途径。