本方案描述了一种用于追踪颞叶癫痫(TLE)患者内侧网络测量指标演变过程的方法。该方法基于颅内记录与一种新颖的数值数据分析技术的结合。具体而言,我们提出了一种针对卵圆孔记录数据的网络分析方案。
本方案描述了一种用于追踪颞叶癫痫(TLE)患者内侧网络测量指标演变过程的方法。该方法基于颅内记录与一种新颖的数值数据分析技术的结合。具体而言,我们提出了一种针对卵圆孔记录数据的网络分析方案。
约30%的癫痫患者对抗癫痫药物耐药。对于这些患者,手术是消除或控制癫痫发作的唯一替代方案。然而,即使在癫痫发作疑似起源区域已被准确定位并切除的情况下,仍有相当一部分患者术后继续出现癫痫发作。本文介绍的方案将一种在颞叶癫痫(TLE)患者术前评估中常规使用的临床操作与一种新颖的网络分析技术相结合。该方法可用于评估内侧网络参数的时间演化过程。通过双侧卵圆孔电极(FOE)插入环池,可同时记录颞叶多个内侧区域的皮层电活动。此外,将网络分析方法应用于所记录的时间序列数据,可追踪间歇期及癫痫发作期间内侧网络的时间动态变化。因此,本方案提供了一种独特的方法,用于可视化和量化多个内侧区域之间的相互关系,而非仅关注单一区域的指标。
癫痫是一种致残性疾病,影响全球1%至2%的人口。在大多数情况下,癫痫的标志性症状——癫痫发作——可以通过抗癫痫药物得到完全控制或消除。然而,约30%的癫痫患者对药物治疗具有耐药性。在最常见的癫痫类型——颞叶癫痫(TLE)1中,幸运的是手术是一种有效的替代治疗手段,可改善患者病情。来自荟萃分析的结果显示,在接受切除手术后的最初两到三年内,近三分之二的药物难治性TLE患者可实现无发作2,3,尽管这一比例因多种因素而异,其中最重要的是海马硬化类型2。实现成功治疗的关键步骤是准确定位所谓的癫痫灶,即负责产生癫痫发作的大脑皮层区域,该区域通常位于颞叶内侧。然而,即使在那些癫痫灶已被准确定位并在手术中成功切除的病例中,仍有相当一部分患者术后仍存在癫痫发作,或必须接受严格的抗癫痫药物治疗以控制发作。因此,一种新的研究视角应运而生:研究重点不再局限于孤立的脑区,而是转向皮层之间的相互作用,这已成为核心问题。这种"网络"研究方法基于“连接组”(connectome)概念4,强调不同脑区之间的神经连接,而非孤立结构的作用。这一新范式借助图论得以实现,图论是一种用于研究图的拓扑与统计特性的数学框架,是表达该领域基本发现的合适工具。在此视角下,大脑被视为由节点通过连接构成的网络5-9,其中节点代表电极所覆盖的皮层区域,节点之间的连接则由其同步化程度决定。因此,这种网络方法已被用于分析侵入性电极记录数据,并为理解癫痫发作产生与传播的潜在机制提供了新的信息。
在全球大多数癫痫中心常规使用的多种侵入性神经生理学技术中,卵圆孔电极(FOE)尤为突出。FOE是一种半侵入性技术,因其无需进行开颅手术,从而降低了与手术相关的并发症风险10。此外,FOE位于环池11内的位置使其特别适合记录来自多个参与癫痫发作产生与传播的皮层结构的内侧活动,例如内嗅皮层。因此,自该技术问世以来,已广泛应用于药物难治性颞叶癫痫(TLE)患者的术前评估中。传统上,该技术用于定位表现为发作间期癫痫样尖波和棘波的激惹性放电活动,更重要的是,可精确定位内侧癫痫发作起始区域。
国际抗癫痫联盟(ILAE)分类与术语委员会提出的新定义认为,癫痫发作起源于特定神经网络中的某个部位12。此外,多项研究已表明,癫痫发作是由异常的网络活动引起,而非由孤立的病理区域所致13-16。显然,这一新观点要求采用新的数值分析方法(如复杂网络方法)对先前获取的信息进行重新分析。尽管这些分析方法在临床实践中的应用仍处于初步阶段,但已有若干研究证实了其应用价值13-17。
下述方案结合了临床上对耐药性颞叶内侧癫痫患者常规实施的医疗操作与一种新型的网络分析技术。该方法可用于评估颞叶内侧网络参数随时间的演变过程。通过将 FOE 双侧插入环池,可同时记录颞叶多个内侧区域的皮层电活动。对记录所得时间序列应用网络分析方法,可追踪癫痫间期及发作期间颞叶内侧网络随时间变化的动态特征。因此,本方案提供了一种独特的方法,用于可视化并量化多个颞叶内侧区域之间相互关系的相关指标。
在下述方案中,步骤1、2和3属于研究和临床方案共有的部分,所有仅根据临床标准筛选出的前颞叶癫痫手术切除候选患者均需严格遵循这些步骤。步骤4和5仅属于研究方案。两种程序均符合Princesa医院伦理委员会的指导方针。
1. 植入前操作步骤
2. 植入操作(手术)
3. FOE 记录的获取
4. FOE 信号预处理
5. 后处理计算(复杂网络分析)
注意:计算每个时间窗内以下所述的指标,时间窗从癫痫发作前5分钟(60个时间窗)开始,至发作后5分钟(60个时间窗)结束,以可视化其时间演变过程。
FOE 的最终位置位于环境脑池,轴向和矢状位 MRI 可见图1A 上图)。FOE 记录可采集颞叶多个内侧结构的电活动图1A 下层面板)。手术后(图1B 左图),患者被送往视频脑电图室,按照10-20系统放置头皮电极(图1B 右侧)。在视频脑电图室停留期间,患者会接受持续监测,以保存用于进一步分析的头皮和颅内电极脑电图记录,以及视频和生命体征数据。典型的原始头皮和颅内电极信号(图1C显示左侧额眶区(FOE)癫痫发作的出现及其向头皮和右侧额眶区(FOE)电极接触点的扩散。
使用兴奋性(S)表征致痫活动图2)对应于原始脑电图记录的 图1C, 从先兆期向发作期及发作后期过渡的过程中。发作起始以一条实心垂直线标记,时间(x轴)以此点为参考。S(兴奋性)的值 >2.5 代表激惹性或致癫痫样活动17,25,26较高的兴奋性(偏红色区域)首先在左侧眶额叶皮层接触点(LFOE)上表现出更强的强度。该结果与神经生理学专家报告的左侧内侧颞叶癫痫一致。
在从发作前期到发作期及发作后期的过渡过程中,多个网络指标以及时频谱熵(图3)的时间动态变化,对应于 图1C 和 2 中所示的同一发作事件。发作起始时间以一条实心垂直线标记,时间(横轴)以此为参考点。本例中,网络基于全部电极构建,包括头皮电极和颅内电极。在发作期间,DoL 和 ACC 值升高,同时 APL 和 Mod 值降低,提示整体连接性增强。在此期间,SE 水平也降低,并在兴奋性消失后(垂直虚线)持续保持较低水平。
在从先兆期向发作期及发作后期过渡过程中,对每个半球功能区(右侧和左侧)(图4)的皮层网络指标——ACC、DoLs 和 APL 以及 SE 的分析。癫痫发作起始时间以一条实心垂直线标记,时间(x 轴)以此时间点为参考。这些指标的变化过程对应于图1、图2和图3中同一发作事件。同侧(左侧)内侧 ACC、DoLs 和 APL 的变化早于对侧且幅度更大,这可由癫痫发作起始区位于左侧颞叶加以解释。本例中,由于缺乏细分区域,无法计算 Mod。
一段代表性视频展示了与图1、2、3和4同一发作期间的功能连接性(图5),在发作起始(时间0)后立即出现显著变化。此时,所有电极之间的连接性急剧增强,表现为连接数量增多以及边的粗细(强度)增加。这种增强最初出现在左侧FOE区域,时间为0.1至0.2,随后扩散至对侧,最终波及整个网络。
| fft | 4.3 | (stats 包)计算信号的快速傅里叶变换。 |
| ts | 4.4 | (stats 包)创建多变量时间序列对象(mts),需提供采样频率。 |
| 兴奋性 | 5.1.1 | (自编函数)基于 R 的 diff 函数,计算信号斜率的绝对值,并将其归一化至短基线期的标准差,需提供阈值。 |
| 功率谱密度与谱熵 | 5.1.2 | (自编函数)基于 R 的 spectrum 和 entropy 函数,计算归一化的功率谱及其香农熵。 |
| ccf | 5.2.1 | (base 包)利用皮尔逊相关系数在零滞后时计算 mts 对象的线性互相关,生成相关矩阵,应计算其绝对值。 |
| graph.adjacency | 5.2.2 | (igraph 包)创建一个 igraph 图,作为后续 igraph 函数使用的基本对象。 |
| average.path.length | 5.2.3 | (igraph 包)通过计算网络中所有节点间最短路径的平均步数,确定图的平均路径长度。 |
| graph.density | 5.2.3 | (igraph 包)通过计算实际连接数与网络中所有可能连接数的比值,计算图的连接密度。 |
| modularity | 5.2.3 | (igraph 包)通过计算哪些节点组之间的连接比与其他网络节点的连接更紧密,确定图的模块性。 |
| transitivity | 5.2.3 | (igraph 包)通过计算彼此相邻的邻接节点所占比例,确定图的平均聚类系数。 |
| 相位同步 | 5.2.4 | (自编函数)基于 R 的 fft 函数,计算平均相位一致性,以获得介于 0 和 1 之间的数值。 |
| smd | 5.3 | (MBESS 包)通过计算组间均值差异相对于合并标准差的比值,确定标准化均值差(效应量)。 |
表1:用于数据处理的R函数。

图 1:卵圆孔电极。(A) 卵圆孔电极(FOE)置入脑池后的最终位置。上排图像显示轴向(左)和矢状面(右)MRI 图像,标示出 FOE 电极接触点的位置(白色箭头)。下排为插入 FOE 的人类标本(尸体),电极接触点以白色箭头标示。(B) FOE 与头皮电极的布设。左图为患者在 FOE 插入术后头部情况,右图为视频脑电图(video-EEG)监测期间的头部电极布置。(C) FOE 与头皮记录图。来自一名左侧颞叶癫痫(TLE)患者的复杂部分性发作(发作前后各 5 分钟)。RFOE1–RFOE6 表示右侧 FOE#1 至 #6,LFOE1–LFOE6 表示左侧 FOE#1 至 #6。发作起始时刻以垂直红色线条和白色箭头头标示。请点击此处查看该图的放大版本。

图 2:通过兴奋性量化的一例左侧颞叶癫痫患者复杂部分性发作的表征。 颜色标尺量化了每个电极的兴奋性水平(S)。右侧卵圆孔电极(RFOE)和左侧卵圆孔电极(LFOE)分别代表右侧和左侧卵圆孔电极的记录触点(y轴)。x轴表示相对于癫痫发作起始时间(由专业神经生理学家确定,以垂直粗线表示)的时间(单位为分钟)。请点击此处查看该图的放大版本。

图 3: 同一患者、同一发作事件(如图2所示)的全网络(头皮+FOE)指标。 图中展示了全网络(头皮+FOE)的平均聚类系数(ACC)、平均路径长度(APL)、连接密度(DoLs)、模块度(Mod)和谱熵(SE)。垂直虚线表示兴奋性(S)。横轴表示相对于发作起始的时间(粗实垂直线)。连续十个时间窗的移动平均值以粗实黑线表示。请点击此处查看该图的放大版本。

图 4:来自图 2 和图 3 同一患者的中线测量指标。 左侧和右侧卵圆孔电极(FOEs)的平均聚类系数(ACC)、平均路径长度(APL)、连接密度(DoLs)和谱熵(SE)。垂直虚线标记了兴奋性时间点。x 轴表示相对于癫痫发作起始的时间(粗的垂直实线)。连续十个时间窗的移动平均值由粗的黑色实线表示。请点击此处查看该图的放大版本。

图 5: 复杂部分性癫痫发作期间连接模式的动态变化。 连接强度通过边的粗细表示。时间(下方数字)相对于癫痫发作起始时间(时间 0)计算。每一帧持续 5 秒。左侧和右侧卵圆孔电极(L1–L6 和 R1–R6)分别用珊瑚色和蓝色圆圈表示。左侧和右侧头皮电极分别用橙色和青色圆圈表示。请点击此处下载该视频。
传统上,癫痫的研究采用基于脑区的方法,将特定脑区——尤其是癫痫发作起始区——孤立出来,视其为癫痫发作的唯一病因。然而最近,强调整皮层区域之间相互作用重要性的真正网络化方法,已逐渐取代经典的基于脑区的观点13-17,28。尽管如此,目前将癫痫视为一种网络性疾病的证据仍高度分散,尚需更多研究。本研究旨在以复杂网络的视角,重新分析传统方法(如致痫区定位)所提供的数据。本文所提出的方案,详细描述了一种对颞叶癫痫患者半侵入性脑电记录进行复杂网络分析和谱分析的逐步方法学流程。
上述技术的应用表明,与传统的局部化或区域导向的视角相比,网络方法具有更高的实用性。在近期的研究中17,29,采用与本文所述完全相同的方法,发现难治性颞叶癫痫(TLE)患者存在内侧连通性的失衡。在发作间期29和发作期17,29,同侧半球的内侧连通性均有所降低。仅通过观察癫痫活动起源的脑区,无法预测这一结果。这一略显意外的发现也通过基于fMRI信号的网络理论得到了描述30,31。此外,FOE与网络理论相结合的技术应用表明,癫痫发作期间的内侧活动与使用促癫痫活动药物(如依托咪酯的药理学给药)后的内侧活动具有等效性32。
本文所述技术能够在持续最多一至两小时的短暂发作间期记录中检测出内侧网络的失衡29。通过这种方法,可显著缩短分析时间并减少患者住院时长。此外,从治疗角度来看,可借鉴深部脑刺激的方法,利用由神经外科医生长期植入的设备,对颞叶癫痫(TLE)患者存在的网络失衡进行"干预调节"。
为获得本方案中所述信息的最佳结果,需预先考虑若干问题。首先,电极的植入应由经验丰富的神经外科医生操作,因为电极位置不当可能导致严重的神经功能后果及误导性的记录结果。此外,用于后续分析的合适时间段的选择完全依赖于神经生理学家对原始脑电图(EEG)的解读,因此必须具备临床脑电图分析的经验。脑电图仪导出文件的数据格式取决于具体品牌,因此需要具备良好的编程技能,以使脚本适应不同的数据格式。最后,为确保数据的可靠性,应对结果实施质量控制。当进行大量相关性分析时,可能出现高估和假阳性结果。在此情况下,应采用统计学方法以提高敏感性。为此,重要的是设定相关性阈值,以剔除不能代表真实潜在同步性的数值。因此,在本方案中,节点间的一条边 i 和 j 仅当这些节点间相关性的绝对值大于0.5时,才认为其存在,该阈值标准此前已被采用17,26. 应采用0.2至0.8范围内的其他阈值以验证结果的一致性,并确保从一个阈值到下一个阈值的平滑过渡。除阈值法外,还可采用其他方法以获得可靠结果,例如Bonferroni校正或替代数据检验。此外,在处理脑电图(EEG)数据时,需注意脑网络是具有非线性动力学特性的复杂系统;因此,除线性相关外,还应使用其他非线性同步性度量方法以确保结果质量,例如互信息或相位同步。33.
直接根据头皮电极计算连接性(如本研究中部分采用的方法)存在一些风险。主要问题在于容积传导带来的污染效应,这种效应在头皮记录中始终存在。解决此问题的一种方法是在源空间进行分析,这是许多研究人员采用的一种可行替代方案。另一种方法是使用能够最小化振幅效应污染的同步性度量指标。通过采用相位同步(也称为相位锁定值),我们可以最小化容积传导的影响,这一点已在多项研究中得到证实34。
与其他侵入性神经生理学技术一样,无法从对照受试者获取颅内电极脑电图(FOE)记录,这一事实严重限制了某些研究方案的应用。FOE记录所提供的数据能够揭示内侧颞叶活动的宝贵信息17,29,35,特别是在颞叶癫痫(TLE)患者中对致痫侧进行定侧分析时尤为有用33。与侵入性技术相比,FOE技术对大脑无创伤,操作相对简便,并且可在较长时间内获得高质量的记录数据11。与磁共振成像(MRI)相比,FOE记录具有更高的时间分辨率,能够更精确地反映皮层电活动。此外,除本研究中所采用的指标外,还存在多种其他可探索的分析指标。这些特点也增加了同时分析多种生物医学信号的可能性。FOE记录技术结合复杂的网络分析与谱分析方法,使其成为癫痫研究中一种强有力的工具,并在临床实践中具有广阔的应用前景。
作者声明无利益冲突。
本工作由卡洛斯三世健康研究所通过项目PI10/00160和PI12/02839资助,部分资金来自欧洲区域发展基金(F.E.D.E.R)和马德里互助会。A.S-G. 获得了马德里互助会的博士后奖学金。3D模拟使用BioDigital Human软件(www.biodigital.com)和ZygoteBody Professional软件(www.zygotebody.com)创建
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 卵圆孔电极 | AD-Tech, Racine, USA | FO06K-SP10X-000 | 六触点铂金 |
| 脑电图仪 | XLTEK, Canada | XLT-EEG32T | Natus XLTEK |
| 磁共振成像设备 | 通用电气 | ||
| 单光子发射计算机断层成像设备 | 通用电气 |