基于贝叶斯优化增强的机器学习进行骨肉瘤风险分层的研究:基于鞘脂代谢

Bayesian Optimization-Enhanced Machine Learning for Osteosarcoma Risk Stratification Based on Sphingolipid Metabolism.

作者信息Yujian Zhong, Ruyuan He, Zewen Jiang, Queran Lin, Fei Peng, Wenyi Jin
PMID40688111
期刊Hum Mutat
发布时间2025-07-11
DOI10.1155/humu/2904964

实验完整度

包括基于TARGET队列的机器学习建模与验证、临床样本免疫荧光验证、体外细胞功能实验及体内异种移植瘤模型实验,涵盖多层级验证。

主要模型

TARGET骨肉瘤队列(85例) 143b骨肉瘤细胞系 U2OS骨肉瘤细胞系 143b和U2OS异种移植瘤裸鼠模型

重点核对

TARGET队列85例患者,随机按8:2分为训练集和测试集(固定随机种子) 特征选择:单变量Cox回归(p<0.01)后Elastic Net(最佳alpha=0.000391)筛选22个基因 SNEX超参数:lambda=0.000125, estimators=1022, max_depth=14, eta=0.848, gamma=0.000146, grow_policy=depthwise siRNA转染浓度:20 nM,作用48小时 细胞系:143b和U2OS,培养于高糖DMEM+10%FBS+1%青链霉素

摘要

背景:异质化的鞘脂代谢驱动骨肉瘤的发生及其促肿瘤微环境。最先进的生物信息学工具,如机器学习,对于通过研究鞘脂代谢的分子和细胞机制来剖析其预后价值至关重要。方法:通过整合Cox回归、5折交叉验证、Elastic Net、极致梯度提升(XGBoost)和贝叶斯优化(用于超参数调整),建立了一套定制的机器学习流程,以构建SM Elastic Net-XGBoost(SNEX)预后模型,并通过Shapley加性解释(SHAP)算法进行解释。进一步研究了SNEX识别的高风险骨肉瘤中分子通路和免疫微环境驱动的预后不良的改变。SNEX的预测结果也得到了临床和实验验证。结果:我们确定了22个关键的SM预后基因用于贝叶斯优化的SNEX。该模型对骨肉瘤患者的预后提供了出色的估计(C指数为1.000)。其稳健性在独立测试集中得到证实,1年、3年和5年的曲线下面积(AUC)分别高达0.875、0.930和0.930。SNEX的决策曲线净收益显著更高,AUC也更高,明显优于以往所有基因预后特征。ACTA2是SNEX阴性预测最关键基因,而BNIP3是阳性预测最关键基因。机制上,SNEX识别的高风险骨肉瘤由于许多关键代谢/炎症/免疫生物学过程失调和免疫抑制微环境而预后不良,14种免疫细胞(巨噬细胞、CD8+ T细胞、NK细胞等)浸润减少。值得注意的是,SNEX突出TERT是最显著的SM预后基因。高表达TERT的临床骨肉瘤比其他肿瘤表现出更显著的恶性特征,如其更高的增殖效率。此外,所有体外和体内实验均验证,抑制TERT丰度会降低骨肉瘤细胞的增殖、侵袭和迁移能力。结论:本研究是首次使用定制的机器学习流程来剖析鞘脂代谢在骨肉瘤中的预后价值和作用的报道。本研究构建了SNEX用于风险分层,具有出色的准确性,并提供了对骨肉瘤中SM介导的通路和微环境失调的深入见解。

实验结论

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

研究问题
鞘脂代谢相关基因能否用于骨肉瘤的风险分层,以及其影响预后的潜在机制。
核心机制
SNEX识别的高风险骨肉瘤中,鞘脂代谢相关基因(如TERT)表达异常,导致代谢/炎症/免疫通路失调,形成免疫抑制微环境,免疫细胞浸润减少,从而促进肿瘤恶性进展。
主要证据
基于TARGET队列(85例)建立SNEX模型,在训练集C-index为1.000,测试集AUC为0.875(1年)、0.930(3年)、0.930(5年);临床样本免疫荧光显示TERT高表达与增殖相关;体外实验(CCK-8、EdU、集落形成、Transwell、划痕实验)和体内异种移植瘤模型均证明敲低TERT抑制骨肉瘤细胞增殖、侵袭和迁移及肿瘤生长。
研究意义
本研究首次使用定制机器学习流程剖析鞘脂代谢在骨肉瘤中的预后价值,SNEX模型可辅助临床医生进行风险分层,并为骨肉瘤的个体化治疗提供潜在靶点(如TERT)。

研究路径

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

1

数据收集与预处理

获取骨肉瘤患者的基因表达数据和生存信息,并划分训练集和测试集。

从TARGET项目回顾性纳入85例骨肉瘤患者,使用scikit-learn的train_test_split函数按8:2比例随机分为训练集和测试集,固定随机种子以确保可重复性。

2

预后基因筛选

筛选与骨肉瘤预后显著相关的鞘脂代谢基因。

首先通过单变量Cox回归筛选预后相关基因(p<0.01),共得到28个显著基因;随后使用Elastic Net模型(结合Grid Search和5折交叉验证)进行二次筛选,移除共线性或影响较小的基因,最终得到22个关键基因。

3

SNEX模型训练与超参数优化

构建并优化基于鞘脂代谢的XGBoost预后模型。

使用XGBoost算法,结合贝叶斯优化(TPE算法)进行100次超参数调优,以C指数为评估指标、负对数似然为损失函数,最终得到最佳超参数组合并生成SNEX模型。

4

模型性能评估

评估SNEX模型的预测准确性和稳健性。

通过时间依赖性ROC曲线的AUC和C指数评估模型判别能力,用Kaplan-Meier曲线分析风险分层效果,并在测试集中独立验证。

5

模型解释

解释SNEX模型的预测机制,提高临床可理解性。

使用SHAP算法分析每个特征对预测的贡献,识别关键基因(如ACTA2、BNIP3、TERT)。

6

机制探索

探究SNEX识别的高风险骨肉瘤的分子通路和免疫微环境异常。

进行GO富集分析和GSEA分析,比较高低风险组的通路差异;使用ESTIMATE算法和ssGSEA分析免疫细胞浸润情况。

7

临床样本验证

验证SNEX识别的关键基因TERT在临床样本中的表达与恶性特征的关系。

使用免疫荧光检测临床骨肉瘤样本中TERT和Ki67的表达,比较TERT高表达和低表达样本的增殖能力。

8

体外功能实验

验证敲低TERT对骨肉瘤细胞增殖、侵袭和迁移的影响。

在143b和U2OS细胞中敲低TERT(siRNA),通过qRT-PCR和Western blot验证敲低效率,并通过CCK-8、EdU、集落形成、Transwell和划痕实验检测细胞功能。

9

体内实验

验证敲低TERT对骨肉瘤肿瘤生长的影响。

将转染si-TERT或NC的143b和U2OS细胞皮下注射到裸鼠右侧胁腹,定期测量肿瘤体积,最终取肿瘤进行免疫组化检测Ki67和TERT表达。

研究方法

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

产品清单

实验环节名称品牌货号
高糖DMEM培养基----
胎牛血清----
青链霉素----
0.25%胰蛋白酶----
Lipofectamine 3000转染试剂Thermo Fisher--
TERT特异性小干扰RNA----
RNeasy Plus Mini试剂盒Qiagen--
NanoDrop分光光度计Thermo Fisher--
PrimeScript RT试剂盒Takara--
SYBR Green PCR预混液Thermo Fisher--
QuantStudio 3实时PCR系统Applied Biosystems--
PVDF膜----
增强化学发光检测试剂----
CCK-8检测溶液----
EdU试剂----
Triton X-100----
结晶紫----
Transwell小室----
磷酸盐缓冲液----
雄性裸鼠Beijing HuaFukang Experimental Animal Center--
Ki67抗体Abcamab16667
TERT抗体Abcamab230527
Alexa Fluor 488偶联二抗ServicebioGB25303
Cy5偶联二抗ServicebioGB27303
TERT和GAPDH一抗----

关键环节

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

环节核对要点
患者队列
TARGET骨肉瘤队列样本量85例,随机划分8:2,固定随机种子
阅读提示:Methods 2.1节
特征选择
单变量Cox回归p<0.01筛选28个基因;Elastic Net最佳alpha=0.000391,5折交叉验证
阅读提示:Methods 2.2节;Results 3.2节
SNEX模型训练
XGBoost超参数:lambda=0.000125, estimators=1022, max_depth=14, eta=0.848, gamma=0.000146, grow_policy=depthwise
阅读提示:Methods 2.2节;Results 3.2节
体外细胞实验
细胞系143b和U2OS(Procell);si-TERT浓度20 nM,转染48小时;培养条件:高糖DMEM+10%FBS+1%青链霉素
阅读提示:Methods 2.8-2.16节
体内动物实验
雄性裸鼠(北京华阜康实验动物中心);每只注射5×10^6细胞于100μL PBS,皮下接种于右胁腹
阅读提示:Methods 2.17节