在耐药性颞叶癫痫患者卵圆孔电极记录的网络分析

Abstract

癫痫患者的约30％是难治抗癫痫药。在这些情况下，外科手术是消除/控制癫痫发作的唯一选择。但是，患者的一个显著少数继续显示手术后癫痫，即使在癫痫发作的可疑源已被正确本地化和切除这些情况。这里介绍的协议相结合的颞叶癫痫（TLE）的患者的术前评价与网络分析的新技术中常规使用的临床过程。该方法允许对近中的网络参数的时间演化的评估。卵圆孔电极（FOE）的双边插入环池在颞叶数近中区域同时记录electrocortical活性。此外，施加到记录时间系列网络方法跟踪近中网络的时间演化既interictally并在癫痫发作。以这种方式，所提出的协议提供了可视化和量化考虑几个近中的区域，而不是一个单一的区域之间的关系的措施的独特方式。

Introduction

癫痫是一种致残性疾病，影响1 - 世界人口的2％。在大多数情况下，癫痫发作 - 癫痫的标志 - 可以完全控制或者与抗癫痫药物废除。然而，癫痫专利约30％的难治性药物治疗。在最常见的类型癫痫，颞叶癫痫^{（TLE）1，}幸运的是手术是改善病人的状况的有效替代。从荟萃分析的结果表明，几乎三分之二耐药TLE患者是无癫痫发作在切除手术^2,3后的第一个两到三年，虽然这一比例在几个因素而变化，最值得注意的是，海马的类型硬化^2。对于一个成功的结果的一个关键步骤是所谓的癫痫病灶的准确定位，皮层区负责发作的产生，其通常位于梅西亚颞叶l的面积。然而，即使在癫痫病灶已在手术过程中被正确识别和切除的情况下，患者的一个显著少数或者与手术后癫痫保持或必须在严格的抗癫痫药治疗控制发作被放置。因此，一个新的观点却出现了注意力不再仅仅关注偏远地区，而不是皮质的相互作用构成了目前的根本问题。这种“网络”方法在连接组概念^4，重点关注在不同区域之间的神经连接，而不是突出条块结构的作用接地。这一新的模式是在图论，专门用于图形拓扑和统计特性研究的一个数学框架中发现，适当的工具来表达它的基本结论。根据这种观点，大脑被认为是通过链路互连的一组节点

在世界上大多数的癫痫中心常规采用的众多微创神经生理学技术，卵圆孔电极（FOE）尤为显着。 FOE是半侵入性技术，因为没有必要执行开颅，这减少了手术相关的并发症^10。此外，FOE在环池¹¹的位置，使它们从参与发作产生和传播几个皮质结构，如内嗅皮层记录近中活性特别方便。因此，它因为使用它的外观是耐药TLE患者的术前评估普遍。传统上，这种技术用于定位在发作癫痫尖峰和尖锐波的形式刺激活性，更重要的是，要准确地识别近中发作起始的区域。

从委员会提议的分类和术语从对癫痫国际联盟（ILAE）的新定义表明，癫痫发作起源于特定网络¹²内的一些点。此外，几个研究已经证明，癫痫发作是由异常的网络活动，而不是由一个孤立的病变区域^13-16引起的。显然，这个新的角度使用新的数值的方法，如复杂的网络方法需要的先前获得的信息再分析。虽然这些分析的实际使用是在临床实践中仍然初期，几项研究已经证实其值^13-17。

下面描述的协议是对耐药TLE癫痫病人常规进行与网络分析的新技术在临床实践的结合。该方法允许对近中的网络参数的时间演化的评估。 FOE的双边插入环池在颞叶的几个领域近中同时记录electrocortical活动。施加到记录时间序列的网络的做法既interictally和癫痫发作期间跟踪近中网络的时间演化。以这种方式，所提出的协议提供了可视化和量化考虑几个近中区域之间的关系的措施的独特方式。

Protocol

在下面描述的协议中，步骤1,2和3属于两个研究和临床方案，这两者都是严格从仅通过临床标准选择的每个近中TLE候选切除手术，随后。步骤4和5专属于研究协议。这两个程序都是按照医院德拉公主伦理委员会的指导方针。

1.预植入程序

1. 解释了实验过程的参与者，指定哪个点对应的研究，哪些是适用于临床实践，陈述了研究过程中没有办法修改临床操作。要特别注意解释的电极植入手术的潜在风险。获取签署知情同意书的参与者。
2. 对于所有候选人切除手术，术前进行神经学和神经心理学EXAminations ^18。
   1. 20个国际系统与莫兹利的-根据10评估通过发作间期单光子发射计算机断层显像（SPECT）与⁹⁹锝HMPAO，磁共振成像使用25头皮电极（MRI）1.5 T和视频脑电图（V-EEG）病人协议^18。
   2. 在手术前V形脑电图记录期间，逐步锥度从第二天的抗癫痫药物到第四天（每天的剂量大约三分之一）。

2.植入程序（手术）

1. 施用抗癫痫药预可操作，并在全身麻醉下进行手术（3毫克/公斤异丙酚丸剂，接着0.2 - 0.3毫克/千克芬太尼和0.5毫克/千克罗库溴铵）。
2. 将两个六敌人接触用1厘米的中心到中心的距离为双侧使用克氏针技术¹⁹周围的蓄水池。
   1. 将病人Øn个手术台仰卧位，与颈部在15度轻轻地延长。与碘溶液制备患者的面颊，开始在切口部位和向外盘旋，并悬垂立即切口部位周围的区域。
   2. 根据哈特尔的地标²⁰与20号脊椎针头穿刺皮肤：一个入口点约3cm侧向朝向一个点口服合缝立即劣于前-后平面同侧瞳孔和点大致的同侧2.5厘米前在横向平面外耳道。
   3. 推进朝针透视下卵圆孔的区域。使用由透视图像提供横向的观点来确定针尖端的位置。当针通过卵圆孔，取出针芯，用电极取代它，并提前入环池（ 图1A
3. 评估通过在手术室²¹荧光透视成像正确注入;排除渗透到颅底，中椎间孔这一点很重要，如眶下裂（位于前部的卵圆孔）和颈静脉孔（位于后部的话）。这样的错位插管可能会导致严重的神经血管损伤^22。
4. 一旦电极被正确地定位在环境蓄水池，它们固定到与窗帘皮肤。醒来后病人，并导致他或她的恢复室。

3.录音FOE收购

1. 返回病人到V-脑电图室逗留约5.2±2.4天（平均值±SD）。
2. 根据国际10-20体系放置19电极。
   1. 用卷尺测量鼻根（鼻梁）和INION（枕骨粗隆）之间的距离，d标记与标记在中间点（锆石电极的位置）。测量并标记点的鼻根（所述FPZ电极的位置）以上的距离的10％。
      1. 重复同样的步骤为INION（盎司电极的位置），标志着距离从锆石20％，在鼻根和INION方向（分别为Fz的和PZ电极的位置）两者。
   2. 同时测量耳前点之间的距离，并标记的距离的左和右耳前点（T3和T4的电极，分别）在10％以上。然后，标记的距离的20％都T3和T4在直拉方向以上以获得C3和C4的位置。
   3. 使用卷尺到FPZ和盎司在两个电极以上在前面和在01（左）和O2（右）在后面的距离在FP1（左）和FP2（右）的5％链接创建一个圆周上。
   4. 在同一圆周上，在INION方向添加的距离的10％上升至OBT泉F7的位置，加10％，达到T3（它应位于耳前点之间的线的上方），并添加另外10％，得到的T5（O1电极）。标记各电极位置，并重复同样的程序为右（偶数）的电极。
   5. 测量和标记中途F7和Fz的，并在F3方向从FP1向上距离的20％之间的交集（F3电极位置）。重复此过程中的头部的每个象限获得F4（前右位置），P3（左后位置）和P4（背左位置）。
   6. 清洁和干燥皮肤。放置火棉胶的温和量与在各电极杯导电凝胶，并在坦然区域电极定位。干后用吹风机的火棉胶。
3. 全部由电​​线电极（头皮和敌人）的连接电极中，这已经是连接到一个electroencephalographer。确保电极信号都不错，并验证头皮电极IMPE舞蹈是在使用electroencephalographer为10kΩ。
4. 获取数字头皮脑电图（EEG）数据和FOE-数据1,024赫兹使用视频同步electroencephalographer（ⅴ脑电图），和利用带通滤波器在0.5过滤数据 - 100赫兹和一个陷波滤波器（50赫兹）与electroencephalographer。
5. 逐步消除来自第二抗癫痫药物第四天（每天的剂量约三分之一），以增加癫痫发作的可能性。此步骤取决于每个病人的特定的药物处方。
6. 同时使用发作阵发性发作活动约通过识别电极/频道里癫痫元素出现^23，包括慢复合波，多峰，快速棘波，尖波，尖和慢复合波的运行定位ictogenic区慢尖波，尖峰和棘慢波。记录发作起始和结束的时间，以及一个Y以外的临床症状或出现相关的研究。有在患者的头部的电极的位置和在其允许解剖识别所述致癫痫活性的显示位置脑电图软件的头部模型之间的一对一的映射。
7. 当研究完成后，轻轻地向外拉出来，而病人的口腔保持半开除去在V-EEG单位的敌人。 FOE删除后没有系统地进行拍摄，当出现神经系统症状时除外。在这样的情况下，执行紧急计算机断层扫描（CT）扫描。

4. FOE信号预处理

1. 导出在200Hz上存储electroencephalographer ASCII格式以适合于大约30分钟癫痫发作活动的（已经由专家神经生理学家标识）（ 图1C）的数值分析历元中的数据。避免含工件历元，如饱和电活动，肌肉活动，以及electrode位移。
2. 打开使用任何UNIX流编辑器导出的文件，并从导出的数据文件中删除所有的非数字字符，只留下时间戳和通道电压。除进一步的数值分析，修改后的文件。
   注：从现在起，执行使用R封装从R库或自制码（ 表1）所有计算。
3. 基于R软件，安装所需的R程序包，以及修改后的数据文件加载到第r环境。订购的所有通道，指定每一个到包含的所有数据，消除空渠道，他们引用的平均中线的参考法（FZ +锆+ Pz的）/ 3阵列的特定列。
   1. 使用快速傅立叶变换算法（R函数：FFT）并绘制得到的变量，以检查有效去除线频率（大约50Hz）的。使用频域来过滤其他杂散FRequencies可能污染信号。
4. 加载的数据转换为一个多元时间序列对象28列（MTS） - 16头皮和12的敌人-使用R功能的TS。划分MTS对象成每次5秒（在200赫兹1000个数据点）的非重叠时间窗口以减少文件大小和优化计算时间。

5.后处理计算（复杂网络分析）

注：计算在每个时间窗口下面所述的措施，开始痉挛发作（60窗口）前5分钟，在痉挛发作后5分钟（60窗口）结束，用可视化的时空演化的目的。

1. 计算单变量的措施，光谱功率，兴奋性和谱熵为每个列/通道不考虑不同时间序列之间的关系。
   1. 计算兴奋（S）为各V使用根据迅达²⁴提出的公式自制代码电压活动时间系列（见补充文件）。 S> 2.5被认为是癫痫，经验确定的阈值^17,25,26。
   2. （ - 7赫兹4），阿尔法（7 - 14赫兹），β（14 - 30赫兹为每个活动的时间序列，使用自制的代码为德尔塔（> 0.5赫兹和<4赫兹），西塔计算功率谱密度）和伽玛（> 30）。
   3. 计算香农熵使用每个时间序列代替的相应概率的时间序列的功率谱密度的自制代码。平均超过一组电极每一个信道得到的各个频谱熵（SE）的值。香农熵在补充文件中说明。
      注：在SE的减少应该被解释为在光谱的频率的数量的减少，因为SE是光谱的熵。
2. 网络措施
   不E：本节评估不同的时间序列的电极之间的相互作用。
   1. 利用在零延迟（：CCF R函数）计算出的线性交叉相关系数的绝对值计算在每个时间窗口中的每对电压的时间序列之间的功能连接。
      注：要消除同步的非代表值，建立了基于以往的研究^17,25,26的阈值。在这种特殊情况下使用的0.5的阈值。
   2. 安装的igraph - [R包^27。创建一个从邻接矩阵的IGRAPH对象（R功能：graph.adjacency）。使用在先前步骤中获得的相关矩阵，指定的图形被加权和无向。
   3. 在每个时间窗口中计算整个网络（头皮+ FOE）平均路径长度（APL）（R函数average.path.length），并为每个所述四个子网络：左头皮，右头皮，左FOE和正确的敌人。在电子Xactly公司以同样的方式，计算链接密度（DOL）（R功能：graph.density），模块化（MOD）（R功能：模块化），平均聚集系数（ACC）（R功能：传递）。
   4. 重复使用相位同步（自制R代码里面 ），为功能连接，而不是互相关函数的估计的先前步骤5.2.1至5.2.3。
3. 来表示变量的大小变化的影响，计算出标准化平均差（SMD）（R函数从包MBESS：SMD），所述preictal和发作阶段之间以及preictal和发作后阶段之间。
   1. 服用preictal为基线，选择30秒钟（6个值），五分钟前发作起始标志，作为preictal值。可以在癫痫发作过程中以量化的变化，对于preictal阶段，通过使用SMD选择的30个类似的时间窗口。

Representative Results

的FOE的最终位置是在环池，如在轴向和矢状MRI（ 图1A上图 ）看出。仇敌记录电活动从颞叶（ 图1A下图 ）的几个近中结构的接触。手术（ 图1B左图 ）后，患者被调度到视频脑电图室，头皮电极与10根据置于- 20系统（ 图1B右侧 ）。期间留在视频脑电图室，病人被连续监视，保存用于进一步分析头皮和FOE-录音，以及视频和重要常数。一个典型的原始头皮和FOE信号（ 图1C）显示了扣押左侧FOE的外观和其散布于头皮和右FOE联系人。

通过使用致痫活动的表示兴奋性的（S）从preictal到发作和发作后时期的过渡过程中对应于从图1C，原始脑电图记录（ 图2）。发作起始标记用实垂直线和时间（x轴）都提到这一点。 S（兴奋）>的值2.5表示刺激性或致癫痫活性^17,25,26。兴奋性较高（偏红色调）与左侧FOE触点（LFOE）较高的强度首先出现。作为由专家神经生理学家获悉这一结果是一致的左颞叶癫痫。

的几个网络措施时间动态以及从preictal到发作和发作后的阶段，对应于图1C和图2中显示的相同的发作在过渡期间的频谱熵（ 图3）。发作起始标有一个坚实的VERTI校准线和时间（x轴）都提到这一点。在这种情况下，该网络在全组电极，包括头皮和FOE-建造。 DOL和ACC值癫痫发作时均较高，随着APL和MOD的下降，表明整体连通的增加。在此期间还进行了观察和持续的兴奋性（垂直虚线）消失后SE水平较低。

的网络措施ACC，DOLS和APL和SE每个FOE（右和左）（ 图4）从preictal到发作和发作后阶段的过渡期间，分析。发作起始标记用实垂直线和时间（x轴）都提到这一点。这种措施的演进对应于相同的发作图1，图2和3。同侧（左）近中ACC，DOLS和APL提出较早，较contrala高变化或双侧的值，这可能是由在发作起始区的位置在左颞叶进行说明。在这种情况下，国防部不能计算，因为没有细分可用。

功能连接相同扣押图1，2，3时（ 图5），4代表视频呈现只是痉挛发作（时间0）后的一个关键的变化。在这一点，所有电极之间的连接显着增加，如通过增加的链路的数目和厚度（强度）的可见的那个边缘。这种增加左FOE之间开始于时间0.1和0.2，并在到达整个网络之前扩散到对侧。

FFT	4.3	（统计包）计算快速傅立叶TR的信号的ansform。
TS	4.4	（统计包）创建一个多元时间序列对象（MTS）。应提供的采样频率。
刺激性	5.1.1	基于差异 R函数（自制）的功能。计算的信号的斜率的绝对值，然后将它规范化为标准偏差短基线时期。应提供的阈值。
功率谱密度和谱熵	5.1.2	基于频谱和熵 R里面的函数（自制）的功能。计算归一化的功率谱和归一化的功率谱的香农熵
CCF	5.2.1	（基本软件包）通过使用Pearson相关零滞后，产生一个相关矩阵计算山对象的线性互相关。绝对值应该演算ated。
graph.adjacency	5.2.2	（IGRAPH包）创建的igraph图，通过以下的igraph功能使用的基本对象
average.path.length	5.2.3	（的igraph包）确定的曲线图的平均路径长度，通过所有网络节点的计算的沿最短路径的步骤的平均数目。
graph.density	5.2.3	（的igraph封装）通过计算链接的实际数量和网络的所有可能的链路之间的比率计算的曲线图的链接的密度。
模块化	5.2.3	（的igraph包）确定图形的模块化，其中节点组更它们之间连接比与网络的其他节点计算
传递	5。2.3	（的igraph包）确定的曲线图的平均聚类系数，通过计算也彼此的邻居的相邻节点的比例
相位同步	5.2.4	（自制）函数基于在计算平均相位相干性来0和1之间获得的值的fft R函数
SMD	5.3	（MBESS包）确定标准平均差-size effects-通过计算相对集中的差异群体之间平均的差异

表1：R功能，用于数据处理。

图1：卵圆孔电极。 （A）的最终位置FOE进入环池。上图显示的轴向（左）和矢状（右图）MRI图像显示FOE接触位置（白色箭头）。一个人体标本（尸体）与插入的FOE（下图，标有白色箭头触点）。 （B）FOE和头皮电极设置。患者只是FOE插入手术（左图）之后和视频脑电图住宿（右图）中的头。 （C）FOE和头皮录音。从左侧颞叶癫痫患者复杂部分性发作（5分钟后痉挛发作之前）。 RFOE1-RFOE6手段右FOE＃1至＃6和LFOE1-LFOE6手段左FOE＃1至＃6。发作起始由垂直红线和白色箭头标注的头。 请点击此处查看该图的放大版本。

图2：从兴奋量化为左颞叶癫痫病人的一个复杂部分发作的表示。色阶量化每个电极的兴奋性电平（S）。右卵圆孔电极（RFOE）和左卵圆孔电极（LFOE）分别代表的右触点和左卵圆孔电极（y轴）。 x轴标志着时间（分钟）相对于痉挛发作（粗垂直线）由专家确定的神经生理学家。 请点击此处查看该图的放大版本。

图3： 整个网络（头皮+ FOE）从图2的平均聚集系数（ACC），平均路径长度同一患者同一和扣押措施 （APL），链接（DOLS），模块化（MOD）和整个网络（头皮+ FOE）谱熵（SE）密度表示。垂直虚线表示的兴奋性（S）。 x轴标志着时间相对痉挛发作（垂直粗实线）。在连续十窗户移动平均线是由一个厚厚的黑色实线表示。 请点击此处查看该图的放大版本。

图4：从图2和3的平均聚集系数（ACC），平均路径长度（APL） 同一患者的近中的措施 ，为左和右孔两条链路（DOLS）和频谱熵（SE）密度未闭的电极（敌人）。垂直虚线标记的兴奋性。 x轴相对于与s的时间标记eizure发作（垂直粗实线）。在连续十窗户移动平均线是由一个厚厚的黑色实线表示。 请点击此处查看该图的放大版本。

图5： 动态连接模式的复杂部分发作期间。链接强度由边缘的厚度来表示。时报讯（数字越小），是相对于痉挛发作（时间0）。每一帧是5秒长。左，右卵圆孔电极（L1-L6和R1-R6）是由珊瑚和蓝色圆圈，分别代表。左，右头皮电极由橙色和青色圆圈，分别代表。 请点击此处到Download这部电影。

Discussion

传统上，癫痫是一种面向区域的方法，其中分离出特定领域，基本上是发作起始区的重要性，因为癫痫发作的病因独特涉猎。最近，一种强调皮层区之间相互作用的重要性一个真正的网络方法已被青睐在经典的面向区域的透视^13-17,28。然而，癫痫作为网络疾病的证据目前身体仍然高度分散，而且需要更多的研究。目前的工作旨在通过重新分析传统方法作为FOE提供的数据，复杂的网络方式下。这里介绍的协议描述一个循序渐进的方法步骤，在TLE患者进行半创记录的复杂网络和频谱分析。

上面描述的技术的应用已表明了网络的方法的有用性是与更传统的禄相比alized或面向区域的观点。在最近的作品^17,29已表明，使用同样的步骤，这里所描述的，在近中连接在耐火TLE患者的不平衡是显而易见的。近中连接在同侧无论是在发作间期²⁹和发作^17,29阶段减少。这个结果不能仅仅通过在那里的致癫痫活性发生的地区寻找可以预见。这在某种程度上令人惊讶的结果也利用功能磁共振成像上信号^30,31网络理论描述。此外，FOE +网络理论的组合技术的应用已经显示出近中活动的等价期间发作和下致癫痫活性的启动子的影响，因为它是依托咪酯³²的药理学施用。

此处所描述的技术能够在短间期录音检测近中网络不平衡至多一个O持久- [R2小时^29。以这种方式，在分析中的时间和病人住院时间大幅减少，可以实现的。另外，从治疗的角度看，在TLE患者现有的不平衡可以通过使用长期植入（由神经外科医生）设备，多，因为它是在深部脑刺激做的方式“解决”。

要使用在本协议中提供的信息，取得最佳的效果，有些问题应提前考虑。首先，电极植入应由有经验的神经外科医生来进行，因为他们的不正确的位置可能会产生严重的神经系统的后果和误导性的记录。此外，为进一步分析相应时期的选择完全依赖于原始脑电图的神经生理学家的解释;因此，在临床脑电图分析方面的经验是强制性的。从脑电图导出文件的数据格式取决于汉邦丘拉尔品牌;因此，需要良好的编程技能相适应的脚本以不同的数据格式。最后，为了确保数据的可靠性，质量控制应该被应用到的结果。过度估计和误报可能具有大量的相关性的工作时出现。在这种情况下，统计方法，以提高灵敏度，应使用。在这方面，必须建立在相关性的阈值丢弃并不代表真正的底层同步的值是重要的。因此，在该协议中，节点i和j之间的边缘将仅视为存在如果这些节点之间的相关性的绝对值大于0.5，以前使用^17,26的标准。在0.2至0.8的范围内的其他的阈值应该被用来验证相似的结果，并确保从一个阈值到以下阈值的平滑过渡。除了阈值，其他methodologie之间可以用于获得可靠的结果，如Bonferroni校正或替代数据的测试。此外，脑电图数据时，必须记住，大脑网络与非线性动力学复杂的系统是非常重要的;因此，除了线性相关，应使用其它非线性同步措施，以确保结果，诸如交互信息或相位同步³³的质量。

直接从头皮上的电极连接的计算，因为它是在这项工作完成了一部分可能会带来一定的风险。主要的问题，休息，由于体积传导，总是存在用头皮记录污染的效果。克服这个问题的一种方式是通过在源的空间，由许多研究中采用的吸引人的选择工作。另一种方法要求使用同步措施的影响幅度污染降至最低。通过使用相位同步（也称为相大号ocking值）我们最小化体积传导的影响，因为它是在几个作品³⁴证实。

如在其他侵入性神经生理技术，从FOE记录不能从对照受试者得到，这一事实严重地限制使用某些研究协议。从FOE录音数据提供有关内侧颞叶活动^17,29,35有价值的信息，尤其是不对称的致痫侧颞叶癫痫病人³³时。与侵入技术相比，FOE技术是非外伤性的脑并涉及相对简单的操作，其记录的高品质在一段¹¹长时间。相较于MRI，FOE录音提供electrocortical活动更好的时间分辨率。此外，存在许多可能性，探索比在这项工作中使用的其他措施。这些事实也增加了一些分析生物医学记录的可能性同时。 FOE录音用复杂的网络和频谱分析相结合的这些优点使得该技术对癫痫研究与临床实践中的应用潜力的强大工具。

Materials

Name	Company	Catalog Number	Comments
Foramen Ovale Electrodes	AD-Tech, Racine, USA	FO06K-SP10X-000	Six-contact platinum
Electroencephalograph	XLTEK, Canada	XLT-EEG32T	Natus XLTEK
MRI machine	General Electric
SPEC machine	General Electric