Richard C D, Tanenbaum A, Audit B, Arneodo A, Khalil A, Frankel W N
The Jackson Laboratory, Bar Harbor, ME 04609 USA; Graduate School for Biomedical Sciences and Engineering, University of Maine, Orono, ME 04469 USA.
Department of Neurology, School of Medicine, Washington University, St. Louis, MO 63130 USA; CompuMAINE Lab, Department of Mathematics, University of Maine, Orono, ME 04469 USA.
J Neurosci Methods. 2015 Mar 15;242:127-40. doi: 10.1016/j.jneumeth.2014.12.016. Epub 2014 Dec 27.
Spike-wave discharges (SWD) found in neuroelectrical recordings are pathognomonic to absence epilepsy. The characteristic spike-wave morphology of the spike-wave complex (SWC) constituents of SWDs can be mathematically described by a subset of possible spectral power and phase values. Morlet wavelet transform (MWT) generates time-frequency representations well-suited to identifying this SWC-associated subset.
MWT decompositions of SWDs reveal spectral power concentrated at harmonic frequencies. The phase relationships underlying SWC morphology were identified by calculating the differences between phase values at SWD fundamental frequency from the 2nd, 3rd, and 4th harmonics, then using the three phase differences as coordinates to generate a density distribution in a {360°×360°×360°} phase difference space. Strain-specific density distributions were generated from SWDs of mice carrying the Gria4, Gabrg2, or Scn8a mutations to determine whether SWC morphological variants reliably mapped to the same regions of the distribution, and if distribution values could be used to detect SWD.
To the best of our knowledge, this algorithm is the first to employ spectral phase to quantify SWC morphology, making it possible to computationally distinguish SWC morphological subtypes and detect SWDs.
RESULTS/CONCLUSIONS: Proof-of-concept testing of the SWDfinder algorithm shows: (1) a major pattern of variation in SWC morphology maps to one axis of the phase difference distribution, (2) variability between the strain-specific distributions reflects differences in the proportions of SWC subtypes generated during SWD, and (3) regularities in the spectral power and phase profiles of SWCs can be used to detect waveforms possessing SWC-like morphology.
神经电记录中发现的棘波放电(SWD)是失神癫痫的特征性表现。SWD的棘波复合波(SWC)成分的特征性棘波形态可以通过可能的频谱功率和相位值的一个子集进行数学描述。莫雷小波变换(MWT)生成的时频表示非常适合识别与该SWC相关的子集。
SWD的MWT分解显示频谱功率集中在谐波频率处。通过计算SWD基频与第二、第三和第四谐波的相位值之间的差异,确定SWC形态背后的相位关系,然后使用这三个相位差作为坐标,在{360°×360°×360°}相位差空间中生成密度分布。从携带Gria4、Gabrg2或Scn8a突变的小鼠的SWD中生成特定品系的密度分布,以确定SWC形态变异是否可靠地映射到分布的相同区域,以及分布值是否可用于检测SWD。
据我们所知,该算法是第一个采用频谱相位来量化SWC形态的算法,使得在计算上区分SWC形态亚型和检测SWD成为可能。
结果/结论:SWDfinder算法的概念验证测试表明:(1)SWC形态的主要变化模式映射到相位差分布的一个轴上;(2)特定品系分布之间的变异性反映了SWD期间产生的SWC亚型比例的差异;(3)SWC的频谱功率和相位分布规律可用于检测具有SWC样形态的波形。