预测癫痫治疗中丙戊酸盐反应的多模态机器学习原型:一项转化研究

A translational multimodal machine-learning prototype predicting valproate response in epilepsy treatment

作者信息Simeon Platte, Afsheen Kumar, Giorgia Guerini, Massimo Pandolfo, Denise Haslinger, Colin B Josephson, Guillermo Delgado-García, Navprabhjot Kaur, Michaela-Pauline Lux, Heiko Stempfle, Chantal Depondt, Reetta Kälviäinen, Felix Rosenow, Sophie von Brauchitsch, Karl Martin Klein, Andreas G Chiocchetti
PMID42159090
期刊Epilepsia
发布时间2026-08
DOI10.1002/epi.70300

实验完整度

研究包含体外iPSC衍生神经元模型、临床队列(发现与验证队列)及多模态机器学习模型开发与验证,但缺少体内动物模型或功能实验,主要依赖临床队列与计算分析。

主要模型

iPSC衍生皮质神经元 Epi25发现队列(比利时、芬兰、德国) 加拿大独立验证队列

重点核对

VPA暴露浓度:0.33 mmol·L⁻¹和1 mmol·L⁻¹(iPSC神经元) 响应定义:≥12个月癫痫无发作(应答者)与>50%复发或无应答停止VPA(无应答者) 模型验证:独立验证队列(n=156)中平衡准确度63%,AUC=0.73 特征选择:Boruta算法选出的22个特征(16个PBSNP、2个PBSEQ、3个临床)

摘要

目的:癫痫影响约1%的全球人口,通常需要终身抗癫痫药物治疗。丙戊酸(VPA)是常用的一线抗癫痫药物,但只有约一半的患者能实现持续的癫痫无发作。治疗选择在很大程度仍凭经验。我们旨在开发并独立验证一个多模态预测模型,以评估对VPA的反应并支持更具个体化的治疗策略。方法:这项横断面治疗反应建模研究使用了国际Epi25队列(比利时、芬兰、德国)的一个子集。纳入标准为接受VPA单药治疗且有可用遗传或临床数据的癫痫患者。发现队列(1965-2021年,58%女性)分为训练集(n=196)和测试集(n=133)。独立验证在加拿大队列(2021-2022年,n=156,40%女性)中进行。主要结局是二分类的VPA反应。应答者定义为归因于VPA达到≥12个月无发作;无应答者定义为>50%的癫痫复发或由于无效、不良反应或不明原因停止VPA。预测算法整合了VPA药代动力学和药效学基因的常见及罕见变异特征、体外神经元VPA反应测量以及临床特征。模型性能通过准确度、预测值(阴性预测值[NPV]/阳性预测值[PPV])和曲线下面积(AUC)评估。结果:在独立验证队列中,多模态分类器实现了63%的平衡准确度(95%置信区间[CI]=52%-73%)、70%的NPV(95%CI=51%-85%)、60%的PPV(95%CI=46%-72%)和0.73的AUC(95%CI=0.63-0.83)。仅使用单一或双数据模态的模型始终表现出较低的预测性能。意义:这项概念验证研究表明,整合遗传、细胞和临床数据能够以具有临床意义的准确度预测VPA治疗反应。尽管尚未准备应用于临床,但这种方法支持基于生物标志物的抗癫痫药物选择的可行性,并有望最终缩短达到有效癫痫控制的时间。

实验结论

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

研究问题
能否整合遗传、细胞和临床数据构建预测VPA治疗反应的模型,以支持更个体化的抗癫痫药物选择?
核心机制
VPA反应与VPA药代动力学/药效学相关基因的常见和罕见变异,以及VPA在iPSC衍生神经元中诱导的差异表达基因网络相关。
主要证据
多模态模型在独立验证队列中达到63%的平衡准确度和0.73的AUC,优于单一模态模型;体外iPSC神经元实验识别出22个VPA响应差异表达基因,包括ATP1A2和KCNA1。
研究意义
该概念验证支持基于生物标志物的抗癫痫药物选择,可能减少试错延迟约60%,但需在更多样化的队列中进一步验证。

研究路径

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

1

VPA响应基因集的体外定义

识别VPA暴露下人iPSC衍生神经元中差异表达基因及共表达网络,以定义基因集1。

将两名健康供体的iPSC分化40天为皮质神经元,在D35暴露于VPA(0.0625-1.0 mmol·L⁻¹),D40提取RNA进行3' RNA测序,通过DESeq2(线性混合模型)和WGCNA分析差异表达基因及模块。

2

药代动力学基因集定义

收集VPA代谢相关基因以构建基因集2。

从DrugBank(ID: DB00313)检索12个VPA代谢相关基因(细胞色素P450和UDP-葡萄糖醛酸转移酶)。

3

遗传特征工程

从常见变异和罕见变异中提取压缩特征以用于预测模型。

基于1000 Genomes数据训练自编码器(AE)生成PBSNP特征,包括基因±2kb的SNP和GTEx脑eQTL,剪枝后;从WES数据中计算PBSEQ图特征,基于VPA响应基因模块移除LGD变异后的网络度量。

4

多模态模型训练与验证

整合遗传、细胞和临床特征预测VPA反应并评估性能。

发现队列(n=329)分为训练集和测试集(60/40),使用SMOTE平衡类别,特征经过Boruta选择,XGBoost分类器训练,在测试集和独立验证队列中评估性能。

研究方法

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

产品清单

实验环节名称品牌货号
StemMACS™Miltenyi Biotec--
神经元基础培养基Thermo Fisher Scientific--
B27添加剂Thermo Fisher Scientific--
层粘连蛋白----
Hoechst 33342Invitrogen--
RNeasy试剂盒Qiagen--
安捷伦2100生物分析仪Agilent--

关键环节

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

环节核对要点
iPSC衍生神经元模型
iPSC来源、分化方案、VPA暴露浓度(0.0625-1.0 mmol·L⁻¹)、暴露时间、RNA提取和测序方法
阅读提示:补充材料Methods S1,图S1和S2
遗传特征生成
SNP筛选标准、自编码器架构、LGD变异识别方法、基因模块定义
阅读提示:材料与方法2.1、2.3、补充表S3
模型训练与验证
训练集/测试集划分、特征选择算法、XGBoost超参数、验证队列预处理
阅读提示:材料与方法2.2、2.5、补充表S4
临床队列
患者入选标准、VPA响应定义、无应答定义、临床变量
阅读提示:材料与方法2.2.3、表1和表2