机器学习增强的单细胞和空间转录组学发现胶质母细胞瘤基底膜相关基因签名

Machine learning-enhanced discovery of a basement membrane-related gene signature in glioblastoma via single-cell and Spatial transcriptomics.

作者信息Ziyuan Liu, Yufan Yang, Haiting Fang, Miaoxi Han, Panpan Qi, Yiqi Fan, Jianlong Li, Shuai He, Lijie Wang
PMID41267067
发布时间2025-11
DOI10.1186/s12967-025-06918-0

实验完整度

高

包含多队列转录组分析、单细胞和空间转录组学、独立队列IHC/mIF验证及siRNA细胞功能实验,多个独立证据层级。

主要模型

U138-MG细胞系 U251-MG细胞系 GBM患者样本(17例独立队列) 四个GBM患者队列(Rembrandt、Gravendeel、Murat、TCGA-GBM)

重点核对

模型使用RSF+Lasso组合,包含FOS、ADM、PLOD3、AEBP1四个基因 Rembrandt队列为训练集,Gravendeel、TCGA-GBM、Murat队列为验证集 siRNA转染浓度及时间:48小时,使用si-FOS1 细胞活力检测时间点:0h、24h、48h、72h 细胞接种密度:4.5×10^3细胞/孔(96孔板)

摘要

背景:胶质母细胞瘤(GBM)的复杂侵袭性和异质性阻碍了肿瘤的完全根除。基底膜(BM)的侵袭发生在肿瘤细胞向脑膜扩散和转移之前,增加了疾病复发率,导致患者预后不良。方法:本研究汇总了四个GBM患者队列,通过二分类机器学习(ML)确认BM候选基因。我们描述了它们的特征,并通过ML组合开发了一种新的BM相关基因(BMRGs)风险模型,并与之前的模型进行比较。同时,使用包含CIBERSORT、xCell和ssGSEA的免疫浸润算法分析高风险和低风险组的免疫特征。基于免疫治疗评估工具TIDE和IRnet的评估以及药物敏感性分析,提出了新的免疫治疗和化疗策略。最后,利用单细胞和空间转录组学揭示高BMRGs表达肿瘤细胞的特征,并通过包含17例患者的42个样本的独立队列进行验证。结果:根据ML组合,四个BMRGs为预测模型提供了基础,用于将患者分为高风险和低风险组。观察到风险评分增加与总生存率降低相关。免疫浸润分析显示,高风险组以M2样巨噬细胞浸润增加和CD8 T细胞减少为特征。同时,注意到高风险组的免疫治疗获益评分降低。药物敏感性分析显示,高风险组对五种药物表现出更高的敏感性。此外,单细胞RNA测序(scRNA-seq)分析发现,BMRGs评分高的肿瘤细胞表现出更有效的细胞间通讯。最终观察结果表明,BMRGs在GBM样本中表达增强并显著共定位。结论:本研究阐明了GBM与BM之间的直接分子相关性,建立了BMRGs签名风险模型,并探讨了它们之间的相互联系。研究结果为GBM的临床结果、治疗反应和潜在治疗策略提供了关键见解。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
基底膜(BM)相关基因能否作为胶质母细胞瘤(GBM)的预后和免疫治疗反应预测标志物,并揭示其潜在机制?
核心机制
FOS作为转录因子可能调节ADM、PLOD3和AEBP1的表达,共同参与基底膜重塑,并通过PI3K、VEGF、缺氧等通路影响肿瘤进展。
主要证据
通过多队列转录组分析构建BMRGs风险模型,结合单细胞和空间转录组学发现高表达BMRGs的肿瘤细胞具有更强的细胞通讯,独立队列IHC和mIF验证共表达,siRNA敲低FOS下调ADM、PLOD3、AEBP1并抑制细胞活力。
研究意义
BMRGs风险模型有望成为GBM患者预后预测和个性化治疗(免疫治疗、化疗)的工具,为BM靶向治疗提供新方向。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

候选基因筛选

鉴定GBM中关键的BM相关基因。

整合四个GBM队列,进行差异表达分析,并结合14种二分类机器学习算法筛选候选基因。

2

BMRGs风险模型构建与验证

构建并验证基于BMRGs的预后风险模型。

基于候选基因进行单变量Cox回归,结合10种ML算法构建117个模型,选择C-index最高的RSF+Lasso模型,并进行生存分析和ROC曲线验证。

3

免疫微环境分析

比较高风险和低风险组的免疫细胞浸润和免疫治疗反应。

使用ssGSEA、xCell、CIBERSORT、TimiGP分析免疫细胞浸润;使用TIDE和IRnet预测免疫治疗反应;分析免疫检查点基因表达差异。

4

药物敏感性分析

预测不同风险组对化疗药物的敏感性。

使用oncoPredict和GDSC数据库预测药物IC50,比较高低风险组的敏感性差异。

5

单细胞转录组分析

阐明BMRGs在单细胞水平的表达特征及细胞通讯能力。

整合4个scRNA-seq数据集,进行细胞聚类、BMRGs评分、细胞通讯分析(CellChat),并比较高低BM肿瘤细胞的通讯强度。

6

空间转录组分析

验证BMRGs在空间转录组中的共定位关系。

使用空间转录组数据(GSE235676),通过semla包分析FOS、ADM、PLOD3、AEBP1的共定位。

7

独立队列验证

通过独立样本验证BMRGs的表达和共定位。

使用17例患者样本进行IHC和mIF实验,验证BMRGs在GBM中的表达和共定位。

8

FOS功能验证

验证FOS对ADM、AEBP1、PLOD3表达的调控作用及其对细胞活力的影响。

使用siRNA敲低FOS,通过qRT-PCR检测目标基因表达,通过CCK-8检测细胞活力。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
杜氏改良Eagle培养基Gibco--
胎牛血清ExCell Bio--
青霉素-链霉素Gibco--
靶向FOS的siRNAHanyi Bio--
CALNPTM RNAi转染试剂D-Nano Therapeutics--
Opti-MEM培养基Gibco--
RNA纯化试剂盒EZBioscience--
SYBR Green预混液EZBioscience--
CCK-8溶液Dojindo--
酶标仪BioTek--
ADM抗体Proteintech10778-1-AP
AEBP1抗体Biossbs-0322R
PLOD3抗体Proteintech11027-1-AP
FOS抗体HuabioHA722666
酪胺信号放大试剂盒----
尼康共聚焦显微镜NikonAX
ImageJ软件----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
细胞培养
细胞系:U138-MG和U251-MG;培养基:DMEM添加10% FBS和1%青霉素-链霉素;培养条件:37°C、5% CO2。
阅读提示:见Methods: Cell culture
siRNA转染
转染试剂:CALNPTM RNAi;使用Opti-MEM培养基;转染时间:48小时;需验证siFOS-1的效果。
阅读提示:见Methods: SiRNA treatment and qRT-PCR
qRT-PCR
使用RNA纯化试剂盒提取RNA,SYBR Green Master Mix进行扩增,ΔΔCt法计算相对表达量。
阅读提示:见Methods: SiRNA treatment and qRT-PCR
细胞活力检测
接种密度:4.5×10^3细胞/孔(96孔板);检测时间:0、24、48、72小时;CCK-8孵育时间:2小时;测量波长450nm。
阅读提示:见Methods: Cell viability assay
免疫组化和多重免疫荧光
抗体稀释度:ADM 1:100,AEBP1 1:500,PLOD3 1:500,FOS 1:2000;使用TSA试剂盒;图像采集使用Nikon共聚焦显微镜。
阅读提示:见Methods: Immunohistochemistry and multiplex Immunofluorescence staining