《Frontiers in Computational Neuroscience》：Quantifying cortex-wide traveling brain waves of complex patterns with a graph-based algorithm
编辑推荐：
行波(TW)是沿神经组织传播的神经振荡。近期，行波检测算法已成功将TW与感知、运动等多种脑功能相联系。然而，现有多数方法并不适用于跨高度弯曲曲面多脑区的大规模记录系统。研究人员开发了一种融合广义相位分析(generalized phase analysis)与
行波(TW)是沿神经组织传播的神经振荡。近期，行波检测算法已成功将TW与感知、运动等多种脑功能相联系。然而，现有多数方法并不适用于跨高度弯曲曲面多脑区的大规模记录系统。研究人员开发了一种融合广义相位分析(generalized phase analysis)与图算法的框架，用于皮层行波检测，可将大规模复杂相位模式分解为多个同时传播的TW，并在单试次数据中量化其属性。该算法被应用于两只执行视觉引导眼跳任务绒猴的半球皮层电图(ECoG，82与122通道)记录。研究发现，眼跳偏移后最初100 ms内枕叶皮层20–50 Hz活动构成宏观TW，波起源于初级视觉皮层(V1)并向喙侧传播至颞叶与顶叶区域；且眼跳后TW起源于V1内与眼跳目标视网膜拓扑表征一致的具体位置，波延迟、时长、空间覆盖与振幅等参数均依赖眼跳方向。结合大规模神经记录，所述图框架为阐释由TW介导的高度动态脑区间通信与协调提供了通用方法。
研究背景与问题提出：局部场电位(LFP)相位在空间上可去同步并沿皮层表面产生时间延迟，从而形成行波(TW)。既往基于解析信号的方法（如相位速度场、相位梯度场、光流、局部圆-线拟合）多面向二维平面数据（钙成像、微电极阵列、局部ECoG），难以直接推广至非平面的脑电图(EEG)、脑磁图(MEG)及跨越高度弯曲曲面的大规模ECoG阵列。非平面数据面临三重困难：传感器分布不均使数值梯度法受限；传播模式超越简单平面/径向/旋转形式，预定义模板分类失效；复杂模式常由多个同时发生的波叠加，需先分解再分类。为此，研究人员在《Frontiers in Computational Neuroscience》提出一种结合广义相位分析与图表示的算法，约束图连接依据皮层表面物理邻近性，以检测距离依赖相位延迟传播的第一类TW，并在模拟与绒猴半球ECoG数据中验证。
主要关键技术方法：研究人员对两只普通绒猴（Marmoset S右半球96通道阵列取82有效通道，Marmoset L左半球128通道阵列取122有效通道）硬膜外植入半球ECoG阵列，采集视觉引导眼跳任务（注视中心环250 ms后跳至周边4°–9°偏心、8方向，眼跳后注视250 ms）下LFP（降采样至1000 Hz）与眼动数据；采用Hilbert变换对20–50 Hz带通滤波单试次LFP求瞬时相位，去除递减段做广义相位线性化；构建以电极为节点、测地距离为权重的邻近图，用最小平均绝对相位差(MAPD)判断边波流方向与对称性，生成有向相位梯度图，从仅出度根节点做子图搜索得最大有向生成树(MDST)，以累积测地距离-相位差线性相关剔除伪子图，跨时拼接同邻根子图形成TW时空表征；模拟数据采用概率深度优先搜索(pDFS)生成各向异性多波并加高斯白噪声评测。
研究结果：
2.1 图算法：研究人员将TW建模为邻域电极间渐变单调相位变化，用周期相似性与相位梯度上限（对应最小传播速度500 mm/s）剪枝邻边；以MAPD平移对齐判波流，建有向图取根节点，子图搜索结合MDST与根到节点累积距离-相位相关，实现单时刻多波分解与跨时追踪。
2.2 模拟基准：在1–4源、噪声标准差10%–50%条件下各生成10000个模拟波。中等噪声（≤20%）约80%被准确判为单波，漏检<10%；时间真阳性率(TPR)近100%，空间TPR>90%、假发现率(FDR)<10%；振幅估计中位偏差约?15%（受带通滤波截止衰减），频率偏差<1%，速度低估随噪声线性增加但30%噪声时仍<10%；额外模拟颞顶交界螺旋波亦被成功检测。
2.3 眼跳后枕叶宏观TW：眼跳偏移后100 ms内枕叶ERP呈三相，V1→V2→V3/V4波谷延迟沿尾-喙轴梯度出现；20–50 Hz相位锁定功率具同向梯度。广义相位图显示单试次中从V1喙向颞顶传播的清晰时空梯度，图算法成功将该TW与同时发生其他活动（绿框）及边界相位重合区（橙框）分离。
2.4 眼跳后TW属性：V1起源最大TW集中于偏移前10至后40 ms；反侧眼跳（Marmoset S 180°、Marmoset L 0°）源位于V1尾内侧电极，向下眼跳（270°）源位于背内侧，符合V1视网膜拓扑；反侧/向下眼跳TW时长（中位~30 ms）、空间覆盖（超视区电极数：S 32、L 45）、根振幅更大，50–100 ms后减弱；中位频率~30 Hz，传播速度1000–2500 mm/s。
讨论部分总结：算法相较既有方法优势为适用任意排布曲面、少模式假设、保广义相位宽频波形、模块化可替换边度量（MAPD可换为滞后相关、Granger因果、传递熵等）、可分解重叠多波；局限为假设根植源/汇（纯平面/旋转波不适用）、恒速假设、噪声偶发伪波需靠试次一致性筛选。生理上该20–50 Hz宏波不同于Kaneko等8 Hz θ波（更慢、平均相位图分析），其V1起源、尾-喙传播匹配眼跳中视网膜扫掠目标外周→中央表征轨迹，也兼容眼跳 corollary discharge 重置视觉区兴奋性的解释。
结论翻译：本研究主要贡献为提供图框架，在非平面记录系统中高效将复杂大规模相位模式分解为多个同时行波。研究人员将该框架实现为算法，成功识别并量化绒猴半球ECoG中眼跳后起源于V1的宏观TW。该框架灵活模块化，可据记录模态与计算需求替换局部传播度量，结合大规模电生理或神经影像数据，可作为研究TW介导的动态脑区间通信与协调的通用方法，进而塑造感知与行为。
