基于脑电图的聚类分析显示脊髓电刺激手术前慢性疼痛患者的显著分离

EEG-Based Clustering Shows Distinct Separation of Chronic Pain Patients Before Spinal Cord Stimulation Surgery

作者信息Ayden Dunn, Jay Gopal, Marisa DiMarzio, Julie G Pilitsis, Ilknur Telkes
PMID42134463
发布时间2026-05-13
DOI10.1016/j.neuroimage.2026.122002

实验完整度

研究包含明确的实验流程(EEG记录、预处理、聚类分析、统计比较),但缺乏体外/动物等独立证据层级,且样本量小,仅进行探索性分析。

主要模型

慢性疼痛患者(计划接受SCS植入,n=16)

重点核对

患者群体:16名慢性疼痛患者,诊断包括PSPS I型和II型 EEG记录条件:60通道,采样率2400Hz,参考耳垂,接地AFz,静息闭眼10分钟 预处理:带通滤波1-150Hz,陷波60Hz和120Hz,ICA去伪迹(MARA),平均参考 聚类方法:K-means (K=3),基于PCA降维,采用肘部法则和轮廓分析确定K值 统计比较:Kruskal-Wallis检验,Tukey-Kramer事后检验,显著性水平α=0.05

摘要

慢性疼痛与皮层活动紊乱有关,但个体在这些神经模式上的差异仍知之甚少。脑电图(EEG)提供了一种非侵入性地表征这些动力学的方法,可能有助于识别与脊髓电刺激(SCS)预后相关的患者亚型。本研究将无监督机器学习应用于术前脑电图,以确定慢性疼痛患者是否表现出不同的神经生理聚类,以及这些组在临床特征或患者报告结局(PROs)上是否存在差异。PROs包括疼痛强度及相关领域(如残疾、灾难化和情绪)的测量,采用数字疼痛评定量表(NRS)、Oswestry残疾指数(ODI)、疼痛灾难化量表(PCS)、麦吉尔疼痛问卷(MPQ)和Beck抑郁量表(BDI)进行评估。对16名计划接受SCS植入的患者记录了静息态头皮脑电图。在标准预处理后,提取频谱特征,并应用K均值聚类(K=3)识别脑电图特征空间内的结构。出现了三个聚类。聚类1的特征是全局α和低β功率降低,以及θ功率和熵降低,提示感觉运动相关区域内振荡活动和频谱复杂性降低的模式。相反,聚类2显示α和低β功率升高,与不同的振荡特征一致,可能反映网络动力学的改变。聚类3表现出中等α和低β功率,以及最高的θ熵,表明神经活动模式具有更大的频谱复杂性和变异性,将此类与其他聚类区分开来,并可能对疼痛处理相关。人口统计学变量在各组间相似,但NRS“最严重”评分存在显著差异(p=0.046)。特征重要性分析确定初级躯体感觉(S1)、次级躯体感觉(S2)、初级运动(M1)和顶枕(PO)区域的峰值低β和α功率是聚类分离的最强贡献者。S1、S2、M1和PO的峰值低β功率显示出最稳健的聚类间差异(所有p≤0.004),S1和M1的峰值α功率也存在显著差异(p≤0.008)。在聚类2和3中,多个脑电图特征与术后PROs改善显著相关,提示症状变化的潜在神经标志物。聚类特异性脑电图特征与术后改善相关,支持脑电图在识别慢性疼痛表型和指导个体化神经调控方法中的潜在效用。

实验结论

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

研究问题
术前静息态脑电图能否通过无监督聚类识别出慢性疼痛患者中神经生理学上不同的亚组,并且这些亚组在临床特征和患者报告结局上是否存在差异?
核心机制
文中未明确说明
主要证据
基于16名计划接受SCS植入的慢性疼痛患者的术前静息态脑电图,提取频谱特征后进行K-means聚类,识别出三个神经生理学上不同的聚类,其主要差异体现在α、低β和θ频段活动上,且这些聚类与疼痛强度(NRS最严重评分)及术后结局的改善存在关联。
研究意义
研究提示基于静息态脑电图的聚类可能捕获传统临床评估未反映的慢性疼痛异质性,为未来个性化脊髓电刺激策略提供基础。

研究路径

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

1

采集静息态脑电图与临床数据

记录慢性疼痛患者术前静息态脑电图,并收集临床、人口统计学和患者报告结局数据。

对计划接受SCS植入的16名慢性疼痛患者记录10分钟闭眼静息态脑电图,同时收集年龄、性别、诊断、病程、阿片类药物使用量(MME)以及多种量表评分(NRS、BDI、ODI、MPQ、PCS)。

2

脑电图预处理与特征提取

去除伪迹并提取反映疼痛神经生理学的频域特征。

使用EEGLAB和自定义MATLAB例程进行预处理,包括带通滤波、陷波滤波、ICA去伪迹(MARA)、平均参考。采用改进的Welch周期图估计中位数频谱,提取θ、α、低β、β、γ频段的功率(相对于高γ频段归一化)、峰值功率、峰值频率和α-θ比值等特征。

3

降维与聚类分析

识别脑电特征空间内的潜在患者亚组。

对标准化脑电特征矩阵应用主成分分析(PCA)降维,基于肘部法则确定主成分数。对降维后的特征集进行K均值聚类,使用肘部法则、轮廓分析和组内平方和选择最优聚类数(K=3)。

4

聚类特征比较与重要性分析

确定区分聚类的关键脑电特征。

计算聚类质心之间的平均欧氏距离作为特征重要性,选择重要性最高的十个特征,使用Kruskal-Wallis检验和Tukey-Kramer事后检验评估聚类间差异。生成地形图可视化关键频段的空间模式。

5

聚类与临床及结局的关联分析

探讨聚类与临床特征及术后改善之间的关系。

比较各聚类的人口统计学和临床特征(使用Kruskal-Wallis检验)。对聚类2和聚类3(因聚类1样本量过小)计算聚类特异性脑电特征与从基线到3个月随访PROs改善之间的相关性(使用MATLAB corrcoef函数)。

研究方法

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

产品清单

实验环节名称品牌货号
g.USBampg.tec medical engineering GmbH--
60通道有源电极蒙太奇----
EEGLAB----
MARA插件----
MATLAB R2024bMathWorks--

关键环节

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

环节核对要点
患者选择与数据采集
患者人群:16名计划接受SCS植入的慢性疼痛患者,诊断包括PSPS I型和II型;年龄43-83岁;病程4-40年。
阅读提示:Methods 2.1节,Table 1
EEG记录
记录条件:60通道,采样率2400Hz,参考耳垂,接地AFz,电极阻抗低于5kOhm,静息闭眼10分钟。
阅读提示:Methods 2.1节
信号预处理
预处理步骤:带通滤波1-150Hz,陷波滤波60Hz和120Hz,ICA(扩展Infomax算法),MARA去伪迹,平均参考。
阅读提示:Methods 2.1节
特征提取
频谱特征:θ、α、低β、β、γ频段相对高γ频段的功率比,峰值功率、峰值频率、α-θ比值;使用1秒汉宁窗、50%重叠的中位数频谱。
阅读提示:Methods 2.1节
聚类分析
聚类方法:PCA降维,K-means聚类(K=3),确定K值的方法(肘部法则、轮廓分析)。
阅读提示:Methods 2.2节
统计分析
统计检验:Kruskal-Wallis检验用于组间比较,Tukey-Kramer事后检验用于两两比较,显著性水平α=0.05;相关性分析使用MATLAB corrcoef函数。
阅读提示:Methods 2.2节,Results 3.3节