Name: 同时MEG / EEG和解剖约束的最小模估计皮质动态映射:听觉注意示例 (视频) | JoVE Journal
Uploaded: 2012-10-24T13:00:00.000+00:00
Duration: 8 min 45 s
Description: 学习使用MEG/EEG与MRI联合技术绘制听觉注意相关的皮层动态图谱。该方法通过最小模估计逆成像和共同表面坐标系统，实现高时空分辨率的脑活动映射，为认知神经科学研究提供有效工具。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

作者感谢Matti S. Hämäläinen、Lilla Zöllei以及三位匿名审稿人提出的有益意见。资助来源：R00DC010196（AKCL）；T32DC000018（EDL）；T32DC005361（RKM）。

1. 解剖数据的采集与处理
<ol>
 <li>获取受试者的一次磁化准备快速梯度回波（MPRAGE）MRI扫描。具体扫描时间取决于所采用的扫描方案，通常需要5-10分钟。</li>
 <li>如果使用脑电图（EEG）数据进行逆成像分析，则需额外获取两次快速低角度激发（FLASH）MRI扫描（翻转角分别为5°和30°），因为FLASH序列提供的组织对比度与标准MPRAGE序列不同1。</li>
 <li>使用FreeSurfer软件（见下表）2, 3重建皮层表面，并建立个体化的M/EEG偶极子源空间。
 <ol>
 <li>该源空间被限制在从MPRAGE扫描中分割出的灰质/白质边界上。每个半球包含约100,000个潜在顶点，间距约为1 mm。在进行偶极子振幅估计（见下文）时，使用7 mm的网格间距，每个半球可得到约3,000个偶极子。</li>
 </ol>
 </li>
 <li>利用MNE（见下表）和FreeSurfer，从MPRAGE和FLASH图像中重建皮肤、颅骨外层和颅骨内层表面。使用这些表面生成三层边界元模型（BEM）。</li>
</ol>
2. M/EEG 数据采集
<ol>
 <li>准备受试者以进行脑磁图/脑电图（M/EEG）记录。
 <ol>
 <li>有关眼电图、参考电极的准备，以及受试者解剖标志点、头部位置指示线圈（HPI）和脑电（EEG）电极的数字化等详细步骤，请参见 Liu et al 4。</li>
 </ol>
 </li>
 <li>当受试者在脑磁图（MEG）设备中就座后，使用 HPI 线圈测量其头部位置。</li>
 <li>开始记录，并启动听觉与视觉刺激的呈现。
 <ol>
 <li>目前有多种硬件和软件方案可用于刺激呈现（例如，Presentation、E-Prime）。我们采用 Tucker-Davis Technologies RZ6 系统进行听觉刺激呈现和触发标记，同时使用 Psychtoolbox 5 在 MATLAB 控制下实现视觉刺激呈现。实验前，可通过连接至屏幕的麦克风和光电二极管测试听觉与视觉刺激的延迟，并确认无明显时序抖动（必要时可将投影仪设置为其原生分辨率），以确保时间精度。</li>
 </ol>
 </li>
 <li>受试者在执行行为任务的同时，通过光学按钮盒对听觉和视觉刺激作出反应。</li>
 <li>保存所有刺激材料、实验参数及数据文件，用于离线分析。</li>
</ol>
3. M/EEG 与 MRI 扫描的共配准及数据处理
<ol>
 <li>使用 MNE 软件加载数字化仪数据和受试者的重建 MRI 头部模型。选取解剖标志点以启动配准过程，并执行自动对齐程序完成坐标变换（图 2）。</li>
 <li>为将源空间中每个偶极子的位置与每个传感器的位置关联起来，结合记录的 HPI 数据（见 2.2）并利用三层边界元模型（BEM）（见 1.4）计算正向解。</li>
 <li>检查所有记录的 M/EEG 数据，识别出方差异常高或完全平坦的通道，并将这些通道标记为坏通道。</li>
 <li>使用信号空间投影 6 或其他降噪技术（例如信号空间分离 7），以投影或分离出由环境背景磁场污染或其他不希望的生理信号（如眼动和心电伪迹）引起的时空场模式（图 3）。
 <ol>
 <li>应用时域伪迹去除方法（例如，移除因通道尖峰导致信号幅值异常高的时段）和频域伪迹去除方法（例如，在 50 或 60 Hz 工频处进行带阻滤波），进一步提高信噪比。</li>
 </ol>
 </li>
 <li>确定一个受试者未执行任何任务的基线期（例如，每次试验开始前 200 ms 的时间段）。对这些基线时段进行平均，以获得噪声估计值（也称为协方差矩阵）。</li>
 <li>确定感兴趣的时段（例如，仅收集行为反应正确的时段），并定义实验对比条件（例如，受试者将其听觉注意力从最初提示的半侧视野切换到对侧半侧视野的时段——“切换”条件——与受试者保持注意力在原始半侧视野的时段——“保持”条件）。为每个定义的条件生成平均响应。
 <ol>
 <li>这些平均值是否进行基线校正取决于实验参数（见 8）；此处展示的数据已进行基线校正。</li>
 </ol>
 </li>
 <li>结合协方差矩阵（3.5）和计算得到的正向解（3.2），获得一个分布式的皮层约束最小范数逆算子，用于将传感器测量值与源空间中偶极子电流估计值关联起来。
 <ol>
 <li>可近似约束或固定偶极子方向为皮层法线方向 9。</li>
 </ol>
 </li>
 <li>为每个实验条件生成分布式偶极子估计的“脑电影”（即，源空间中每个偶极子位置随时间变化的电流估计值）（图 4）。
 <ol>
 <li>根据实验设计的时间特性，可通过使用不重叠的时间窗对电流估计值进行平均，将数据按时间分段。</li>
 </ol>
 </li>
</ol>
4. 基于通用表面坐标的统计推断
<ol>
 <li>基于一种能够最优对齐个体脑沟-脑回模式的表面坐标系，将每位受试者的“脑电影”变形至一个共同（平均）的皮层空间3。这使得我们能够在所有受试者之间比较或平均皮层活动。（图5）。</li>
 <li>存在多种不同的统计推断方法。本文将重点介绍三种可能的方法。未在软件包中实现的方法可通过自定义软件编写——在我们的示例中，使用MATLAB执行非参数时空聚类置换检验。尽管这些数据具有高维度（空间 × 时间 × 受试者），但所有这些方法均可在标准的现代台式计算机硬件上在数秒（感兴趣区域方法；4.3节）到数小时（非参数聚类方法；4.5节）内完成。</li>
 <li>感兴趣区域（ROI）方法
 <ol>
 <li>可依据解剖结构（例如，由自动分区算法定义1）和/或功能特性（例如，通过记录功能定位任务，如“执行/不执行”眼跳任务，以识别眼动控制相关脑区）来定义感兴趣区域。</li>
 <li>可进一步将分析限制在与实验范式相适应的特定时间范围内（例如，声音刺激呈现前后的一段时间）。也可采用与时间序列分析相关的其他统计推断方法。</li>
 </ol>
 </li>
 <li>全脑Bonferroni校正或错误发现率（FDR）校正
 <ol>
 <li>若需进行全脑、全时间点分析，可采用Bonferroni校正或FDR校正。</li>
 <li>在每个偶极子位置和每个时间点，使用适当的检验统计量（如t检验或针对近似正态分布数据的被试内重复测量方差分析）生成统计图。例如，当结合对估计值间相关性的校正（如保守的Greenhouse-Geisser校正）时，可使用MNE估计值在固定偶极子源上的动态统计参数映射所得的z分数10。</li>
 <li>对于Bonferroni校正，通过将显著性水平设为0.05除以比较次数（偶极子数量乘以时间点数量）来确定显著的时空点。若采用较不保守的方法，可使用FDR p值校正11。</li>
 </ol>
 </li>
 <li>非参数时空聚类方法
 <ol>
 <li>使用该方法（基于对文献12的简单扩展），可在控制族系误差率的前提下，识别出具有大范围且一致的空间和时间激活的脑区，其保守性低于Bonferroni校正，且相较于FDR更不易产生I类统计错误。
 <ol>
 <li>由于该方法采用置换检验或蒙特卡洛重采样技术，不依赖于数据正态性假设，仅假设在零假设下条件标签可交换。尽管其计算量大于前两种方法，但在现代台式计算机硬件上，仍可在数小时内于单台机器上完成。</li>
 </ol>
 </li>
 <li>在每个偶极子位置和每个时间点，使用适当的检验统计量（如t检验）生成统计图。</li>
 <li>以初步显著性阈值对该统计图进行阈值化处理，例如 p &lt; 0.05。</li>
 <li>根据时空邻近性对这些初步显著的点进行聚类，例如，时间上相距5毫秒以内、空间上测地距离5毫米以内的显著点归为同一聚类。使用超体积或总显著性（例如，聚类内各点t分数之和）对每个聚类进行评分。</li>
 <li>执行标准的置换重采样检验（对于较大数据集，例如 受试者数量N &gt; 10，为节省计算资源可采用蒙特卡洛重采样），使用最大统计量（参见12中的置换检验示例）。简而言之，从一组随机子集的受试者中（选择0至N名受试者），在获得统计图前重新标记待比较的条件，对新的统计图进行聚类，并获取该次重标记下的最大聚类评分。对最多2N种可能的重标记重复此过程，以获得最大统计量的分布；执行全部2N种可能的重标记即构成置换检验，而仅使用少于2N的随机子集则构成蒙特卡洛（或随机）置换检验。</li>
 <li>通过计算最大聚类评分超过原始聚类评分的比例，确定给定原始聚类（来自原始标签）的显著性，例如，若某聚类大于95%的最大统计量聚类，则可判定其具有显著性。
 <ol>
 <li>关于MEG分布式源成像中统计推断的深入讨论，参见文献13。</li>
 </ol>
 </li>
 </ol>
 </li>
 <li>所得数据文件可通过多种方式可视化，包括使用MNE软件原生支持的用于存储时空皮层估计值的格式，即STC文件。这些文件以及对应显著区域生成的标签，均可通过MNE提供的MATLAB和Python工具箱生成。</li>
</ol>
5. 代表性结果
图6展示了使用图4中所述行为范式获得的一组代表性结果。通过非参数时空聚类分析方法（4.5），发现被试在执行重定向任务时，右侧额眼区（FEF）相较于标准任务具有显著性（图6 左）。采用感兴趣区域（ROI）分析方法（4.3），展示了右侧FEF的时间进程，以及这两种条件之间存在显著差异的时间段。
<img alt="静态平衡示意图；M/EEG 数据处理流程：MRI 输入、传感器位置、统计推断。" src="/files/ftp_upload/4262/4262fig1.jpg" data-altupdatedat="1787630408000" data-alt="图 1"> 
 图 1. 利用皮层约束的最小范数偶极子估计生成“脑电影”的工作流程（参见，Liu et al.，2010 年的图 1）。
<img alt="脑电图示意图；带有传感器位置的脑部图谱；神经活动分析。" src="/files/ftp_upload/4262/4262fig2.jpg" data-altupdatedat="1787630408000" data-alt="图 2"> 
 图 2. 使用 MNE 软件将 EEG 通道和 HPI 位置共配准至一名受试者的 MRI 坐标空间。
<img alt="脑电数据处理示意图，包括滤波和投影；头部模型显示电极放置位置。" src="/files/ftp_upload/4262/4262fig3.jpg" data-altupdatedat="1787630408000" data-alt="图 3"> 
 图 3. 使用SSP去除心电（以橙色标出）和眨眼（以蓝绿色标出）伪迹并进行低通滤波以消除工频干扰前后的脑磁图（MEG）数据。 <a href="https://www.jove.com/files/ftp_upload/4262/4262fig3large.jpg" target="_blank"> 点击此处查看大图</a>。
<img alt="神经激活、视听刺激示意图、dSPM结果、脑活动图谱" src="/files/ftp_upload/4262/4262fig4.jpg" data-altupdatedat="1787630408000" data-alt="图4"> 
 图4. 在受试者原生皮层空间上的“脑部动态图”以及一种实验范式中视听刺激呈现的时间（其中听觉刺激在600 ms时呈现，视觉刺激在-600 ms时呈现）（注：在最终的视频片段中将以动态视频形式展示）
<img alt="脑图谱绘制过程；示意图；受试者数据映射到球形空间以构建平均脑模型。" src="/files/ftp_upload/4262/4262fig5.jpg" data-altupdatedat="1787630408000" data-alt="图 5"> 
 图 5. 假设的感兴趣区域（ROI）在受试者原始皮层空间与变形至公共皮层空间后的对比。
<img alt="右侧FEF激活，随时间变化的电流图分析；重定向反应与标准反应比较。" src="/files/ftp_upload/4262/4262fig6.jpg" data-altupdatedat="1787630408000" data-alt="图6"> 
 图6. 与所测试的两种实验条件相关的代表性时空聚类及时间进程。

<ol><li>Fischl, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Automatically+parcellating+the+human+cerebral+cortex%5BTitle%5D)+AND+%22Cerebral+Cortex%22%5BJournal%5D)">Automatically parcellating the human cerebral cortex.</a> Cerebral Cortex. 14, 11-22 (2004).</li>
<li>Dale, A., Sereno, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Improved+localization+of+cortical+activity+by+combining+EEG+and+MEG+with+MRI+cortical+surface+reconstruction%3A+A+linear+approach%5BTitle%5D)+AND+%22Journal+of+Cognitive+Neuroscience%22%5BJournal%5D)">Improved localization of cortical activity by combining EEG and MEG with MRI cortical surface reconstruction: A linear approach.</a> Journal of Cognitive Neuroscience. 5, 162-176 (1993).</li>
<li>Fischl, B., Sereno, M. I., Dale, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cortical+surface-based+analysis.+II%3A+Inflation%2C+flattening%2C+and+a+surface-based+coordinate+system%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system.</a> NeuroImage. 9, 195-207 (1999).</li>
<li>Liu, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+Mapping+with+Simultaneous+MEG+and+EEG%5BTitle%5D)+AND+%22Journal+of+Visualized+Experiments%22%5BJournal%5D)">Functional Mapping with Simultaneous MEG and EEG.</a> Journal of Visualized Experiments. (40), e1668(2010).</li>
<li>Brainard, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Psychophysics+Toolbox%5BTitle%5D)+AND+%22Spatial+Vision%22%5BJournal%5D)">The Psychophysics Toolbox.</a> Spatial Vision. 10, 433-436 (1997).</li>
<li>Uusitalo, M. A., Ilmoniemi, R. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Signal-space+projection+method+for+separating+MEG+or+EEG+into+components%5BTitle%5D)+AND+%22Med.+Biol.+Eng.+Comput%22%5BJournal%5D)">Signal-space projection method for separating MEG or EEG into components.</a> Med. Biol. Eng. Comput. 35, 135-140 (1997).</li>
<li>Taulu, S., Simola, J., Kajola, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Applications+of+the+signal+space+separation+method%5BTitle%5D)+AND+%22IEEE+Transactions+on+Signal+Processing%22%5BJournal%5D)">Applications of the signal space separation method.</a> IEEE Transactions on Signal Processing. 53, 3359-3372 (2005).</li>
<li>Urbach, T. P., Kutas, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Interpreting+event-related+brain+potential+%28ERP%29+distributions%3A+Implications+of+baseline+potentials+and+variability+with+application+to+amplitude+normalization+by+vector+scaling%5BTitle%5D)+AND+%22Biological+Psychology%22%5BJournal%5D)">Interpreting event-related brain potential (ERP) distributions: Implications of baseline potentials and variability with application to amplitude normalization by vector scaling.</a> Biological Psychology. 72, 333-343 (2006).</li>
<li>Lin, F. -H., Belliveau, J. W., Dale, A. M., Hämäläinen, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distributed+current+estimates+using+cortical+orientation+constraints%5BTitle%5D)+AND+%22Human+Brain+Mapping%22%5BJournal%5D)">Distributed current estimates using cortical orientation constraints.</a> Human Brain Mapping. 27, 1-13 (2006).</li>
<li>Dale, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+statistical+parametric+mapping%3A+combining+fMRI+and+MEG+for+high-resolution+imaging+of+cortical+activity%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity.</a> Neuron. 26, 55-67 (2000).</li>
<li>Nichols, T., Hayasaka, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Controlling+the+familywise+error+rate+in+functional+neuroimaging%3A+a+comparative+review%5BTitle%5D)+AND+%22Statistical+Methods+in+Medical+Research%22%5BJournal%5D)">Controlling the familywise error rate in functional neuroimaging: a comparative review.</a> Statistical Methods in Medical Research. 12, 419-446 (2003).</li>
<li>Nichols, T. E., Holmes, A. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nonparametric+permutation+tests+for+functional+neuroimaging%3A+a+primer+with+examples%5BTitle%5D)+AND+%22Human+Brain+Mapping%22%5BJournal%5D)">Nonparametric permutation tests for functional neuroimaging: a primer with examples.</a> Human Brain Mapping. 15, 1-25 (2001).</li>
<li>Pantazis, D., Leahy, R. M. Statistical Inference in MEG Distributed Source Imaging. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPantazis%2C+D+%22MEG%3A+An+Introduction+to+Methods%22">MEG: An Introduction to Methods</a>. Hansen, P., Kringelbach, M., Salmelin, R. , Oxford University Press. New York. 245-272 (2010).</li>
<li>Ahveninen, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Attention-driven+auditory+cortex+short-term+plasticity+helps+segregate+relevant+sounds+from+noise%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences%22%5BJournal%5D)">Attention-driven auditory cortex short-term plasticity helps segregate relevant sounds from noise.</a> Proceedings of the National Academy of Sciences. , 1-6 (2011).</li>
<li>Sharon, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+advantage+of+combining+MEG+and+EEG%3A+comparison+to+fMRI+in+focally+stimulated+visual+cortex%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">The advantage of combining MEG and EEG: comparison to fMRI in focally stimulated visual cortex.</a> Neuroimage. 36, 1225-1235 (2007).</li>
<li>Herrmann, B., Maess, B., Hasting, A. S., Friederici, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Localization+of+the+syntactic+mismatch+negativity+in+the+temporal+cortex%3A+An+MEG+study%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Localization of the syntactic mismatch negativity in the temporal cortex: An MEG study.</a> NeuroImage. 48, 590-600 (2009).</li>
<li>Baillet, S. The Dowser in the Fields: Searching for MEG Sources. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBaillet%2C+S+%22MEG%3A+An+Introduction+to+Methods%22">MEG: An Introduction to Methods</a>. Hansen, P., Kringelbach, M., Salmelin, R. , Oxford University Press. New York. 83-123 (2010).</li>
<li>Gutschalk, A., Micheyl, C., Oxenham, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Neural+correlates+of+auditory+perceptual+awareness+under+informational+masking%5BTitle%5D)+AND+%22PLoS+Biology%22%5BJournal%5D)">Neural correlates of auditory perceptual awareness under informational masking.</a> PLoS Biology. 6, e138(2008).</li>
</ol>

未声明任何利益冲突。

为了从采集到的脑磁图（MEG）/脑电图（EEG）数据中估算大脑皮层上的偶极子激活，我们需要求解一个逆问题。除非施加适当的解剖学和生理学约束，否则该问题不存在唯一且稳定的解。通过使用个体受试者经磁共振成像（MRI）获得的解剖学约束，并采用最小范数作为估计准则，我们可以得到一个与传感器测量结果一致的皮层电流源逆向估计结果。该方法已被证明在听觉处理14以及其他领域（如视觉处理15和语言处理16）的研究中具有重要应用价值。
还有许多其他逆问题求解方法。然而，所有这些方法均可归纳为两类：定位方法（例如，等效电流偶极子建模）或成像方法（例如，最小范数估计（MNE）、波束成形技术）。此外，每种逆问题求解方法都有其权衡取舍（详见文献17的深入讨论）。例如，由于本方法采用最小范数约束，其电流估计结果在空间上必然是分布式的。这种最小范数估计方法特别适用于涉及广泛皮层网络参与的任务。而对于刺激诱发的局灶性源活动的早期反应映射（例如，听觉刺激诱发的、被认为定位于双...

<table><tbody><tr><td>306通道 Vectorview MEG 系统</td><td>Eleka-Neuromag Ltd,</td><td></td><td></td></tr><tr><td>1.5-T Avanto MRI 扫描仪</td><td>西门子医疗系统公司</td><td></td><td></td></tr><tr><td>FreeSurfer</td><td></td><td></td><td>&lt;a href=&quot;http://freesurfer.net/&quot; target=&quot;_blank">http://freesurfer.net/&lt;/a&gt;</td></tr><tr><td>MNE 软件</td><td></td><td></td><td>&lt;a href=&quot;http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php&quot; target=&quot;_blank">http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php&lt;/a&gt;</td></tr><tr><td>脑电图电极</td><td>Brain Products, Easycap GmbH</td><td></td><td></td></tr><tr><td>3Space Fastrak 系统</td><td>Polhemus</td><td></td><td></td></tr><tr><td>光学按钮盒（FIU-932）</td><td>Current Designs</td><td></td><td></td></tr></tbody></table>

利用同步脑磁图/脑电图与基于解剖约束的最小范数估计绘制皮层动态：以听觉注意为例

Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example

脑磁图与脑电图（MEG/EEG）是神经影像技术，具有高时间分辨率，特别适用于研究参与动态感知与认知任务的皮层网络，例如在鸡尾酒会情境中注意不同声音的过程。以往许多研究仅使用在传感器水平记录的数据， 即，头皮表面及外部记录的磁场或电位，通常关注与刺激呈现时间锁定的脑活动。此类事件相关场/电位分析在仅有少量可在时空上分离和识别的偶极子模式时尤为适用。 alternatively，通过利用解剖学信息，这些不同的场模式可被定位于大脑皮层上的电流源。然而，对于持续时间较长、可能不与特定刺激时间锁定的脑响应（例如，为基于提示的听觉特征从两个同时呈现的 spoken digits 中选择一个做准备）或可能分布在多个未知的空间位置上 先验的仅使用有限数量的局部源可能无法充分捕捉到分布式皮层网络的招募。
在此，我们描述一种利用个体解剖磁共振成像（MRI）数据，通过最小范数估计（minimum-norm estimates, MNE）方法建立传感器信号与大脑皮层偶极子激活之间关系的实验流程。这种逆成像方法为我们提供了进行分布式源分析的工具。为便于说明，我们将以免费提供的 FreeSurfer 和 MNE 软件为例，详细介绍全部操作步骤。我们将总结生成前向模型所需的 MRI 序列及分析步骤，该模型可用于将分布于皮层上的偶极子所产生的预期场模式与 M/EEG 传感器信号关联起来。接下来，我们将逐步介绍对传感器数据进行去噪处理的关键流程，以消除环境干扰和生理噪声的影响。随后，我们将概述将 MEG/EEG 传感器数据整合并映射到皮层空间的操作方法，从而生成与各实验条件相关的、反映大脑皮层表面偶极子激活动态变化的时间序列集合（即“脑部电影”）。最后，我们将重点介绍若干统计分析技术，这些方法可支持我们在被试群体水平上进行科学推断即基于一个共同的皮层坐标空间进行组水平分析。

我们采用脑磁图和脑电图（MEG/EEG），结合磁共振成像（MRI）获取的解剖学信息，来绘制与听觉注意相关的皮层网络动态图谱。

脑磁图与脑电图（MEG/EEG）是两种神经影像技术，具有高时间分辨率，特别适用于研究……

本实验的总体目标是绘制人类不同认知状态下潜在的皮层动态图谱。首先，在受试者执行行为任务时同步采集脑磁图与脑电图数据（简称 MEG 与 EEG）；第二步是采用相应的 MRI 序列获取解剖学数据。接下来进行共配准，以建立MEG和EEG传感器位置与解剖结构信息之间的空间对应关系。最后一步是采用逆成像方法，将个体的脑活动映射到皮层空间。最终，基于通用的基于皮层表面的坐标系进行统计推断，以发现能够区分不同认知状态的显著时空模式。本实验操作由 Eric Lawson 和 Ross Maddox 演示。对于我实验室的博士后研究人员，请首先获取受试者的结构磁共振（MR）图像。首先，使用磁化准备快速梯度回波序列（magnetization prepared rapid gradient echo，MP-RAGE）或类似序列获取结构磁共振扫描图像。该序列可能需要5到10分钟，具体时间取决于所使用的特定扫描分辨率和成像方案，特别是当脑电图（EEG）数据将用于逆成像分析时。此外，还需获取两次快速低角度激发（FLASH）磁共振成像扫描。这些FLASH序列可提供与标准MP-RAGE序列不同的组织对比度。成像完成后，使用 MNE 和 FreeSurfer 软件从 MP-RAGE 和 FLASH 图像中重建皮肤、外颅骨和内颅骨表面。然后利用这些表面在 MEG 实验前生成三层边界元模型（BEM）。首先测试听觉和视觉潜伏期，以确保时间同步的准确性。使用连接到屏幕的麦克风和光电 DDE，并确保无可见的抖动。这可能需要将投影仪设置为其原生分辨率。接下来，准备受试者的记录工作，具体步骤可参考 luital 之前发布的视频文章，了解心电图和参考电极的制备方法，以及受试者解剖标志点、头部位置指示线圈和脑电图电极的数字化过程。当受试者在MEG测量头部位置舒适就位后，使用头部位置指示器（HPI线圈）开始记录，并呈现听觉和视觉刺激。注意，HPI测量也可持续进行。受试者在执行视听行为任务时，应通过光学按钮盒对听觉和视觉刺激作出反应。在此过程中，受试者根据视觉线索提示，报告来自半侧视野的轮辐和数字。有时，在试验过程中会通过视觉提示要求受试者将注意力切换至对侧半视野。为了研究听觉注意力的转换，目前已有多种硬件和软件解决方案可用于刺激呈现。此处使用 Tucker Davis Technologies RZ6 系统呈现听觉刺激，并利用 Psychtoolbox 呈现视觉刺激并标记触发信号，两者均由 MATLAB 控制。开始数据处理时，使用 MNE 软件将 EEG 数据与结构磁共振图像进行共配准：首先将定位器数据载入受试者的重建 MRI 头部模型中；接着选择解剖标志点以启动共配准过程，然后执行自动对齐程序完成坐标变换。然后将源空间中每个偶极子的位置与每个传感器的位置相关联。结合记录的头部位置指示器数据，利用三层边界元模型计算正向解，以进一步提高数据的信噪比。应用时域伪迹去除方法，例如去除由于通道尖峰引起的异常高幅信号所包含的事件片段。还可应用频域伪迹去除方法，例如在50或60赫兹工频处进行带阻滤波，使用信号空间投影或其它降噪技术（如信号空间分离），以投影或分离出由环境背景磁场污染或其他不希望的生理信号（如与眨眼和心电伪迹相关的信号）引起的时空场模式。接下来，针对每种实验条件，生成一个分布偶极子估计的脑部动态图像，即在源空间中每个偶极子位置上的电流估计随时间变化的可视化结果。根据实验设计的时间特性，可通过使用无重叠的时间窗对电流估计值进行时间上的平均，以实现数据在时间维度上的整合，从而继续进行形态分析。先前为每位受试者创建的大脑电影被映射到一个基于表面坐标系统的公共皮层空间，该空间能够最优地对齐个体的皮层功能模式。这使得可以在不同受试者之间比较或平均皮层活动。采用感兴趣区域（ROI）方法时，可通过解剖学方式定义ROI，例如使用自动脑区分割算法，或通过功能定位任务（如“go no-go”选择任务）进行功能定义，以识别眼动相关脑区。分析还可进一步限定于与所用实验范式相适应的特定时间窗口，例如限制在声音刺激呈现前后的一段时间内。其他与时间序列分析相关的统计推断也可使用上述行为范式进行。此处展示的是采用非参数时空聚类方法得到的代表性结果。当个体受试者执行重定向任务与标准任务相比时，发现右侧额眼区具有显著性。采用感兴趣区域（ROI）方法，展示了右侧额眼区的时间进程，以及这两种条件存在显著差异的时间段。观看本视频后，您将能够充分了解如何利用脑磁图-脑电图（M-E-G-E-E-G）和磁共振成像（MRI）来绘制不同行为任务中的皮层动态图谱。通过采用适当的统计学方法，您可以发现能够区分不同认知状态的时空模式差异。感谢观看，祝您的实验顺利。

mapping-cortical-dynamics-using-simultaneous-megeeg-anatomically

Research

科研

JoVE Journal

Neuroscience

神经科学

同时MEG / EEG和解剖约束的最小模估计皮质动态映射：听觉注意示例

我们采用脑磁图（MEG）和脑电图（EEG）来绘制参与简单感觉刺激处理的脑区。...

同步 MEG 与 EEG 的功能定位映射

注意控制包括对目标信号的增强和对干扰信号的抑制。本文描述了一种方法，可在持续跨模态注意过程中，利用一种被动控制条件，分别但同时地测量注意和忽略的神经生理机制；在此被动控制条件下，上述两种过程均未持续参与。...

注意力控制中忽略与注意过程的神经生理信号测量

本文介绍了一种脑电图-功能磁共振成像（EEG-fMRI）多模态成像方法，即时空fMRI约束的脑电图源成像方法。该方法采用条件激活的功能磁共振成像子图（或先验信息），以提高脑电图源定位的空间特异性，并减少错误结果。...

同步脑电图与功能磁共振成像记录及整合分析用于动态皮层活动成像

本研究提出了一种利用神经影像学方法研究仅通过可穿戴设备实现左右反转听觉适应的实验方案，该方法可作为揭示人类在听觉领域对新环境适应能力的有效工具。...

一种研究适应左右反转听觉的方法

本文介绍如何设置一项实验，以利用脑磁图（MEG）检测物体感知前的源水平影响。内容涵盖刺激材料、实验设计、MEG记录以及数据分析。...

利用脑磁图检测物体感知中刺激前源水平效应

脑磁图（MEG）与高密度脑电图（HD-EEG）很少被同时记录，尽管二者可提供相互验证且互补的信息。本文展示了同步记录MEG与HD-EEG的实验设置，以及分析这些数据的方法，旨在定位药物难治性癫痫患儿的致痫区和功能重要脑区。...

药物难治性癫痫患儿术前评估中的电磁源成像

一种结合多种磁共振成像（MRI）数据模态指导下的个性化经颅磁刺激与脑电图（TMS-EEG）皮层网络定位的新方案。...

神经影像引导的经颅磁刺激–脑电图用于实时皮层网络定位

本文介绍了在大鼠和小鼠中进行高密度微电皮层记录（µEcoG）的两种方案，包括手术、植入和记录方法。µEcoG记录分别与大鼠听觉皮层的层状多通道电极记录或小鼠体感皮层的光遗传学神经活动操控相结合进行。...

通过整合层状多通道电极、光遗传学与微阵列皮层电图在啮齿类动物中进行皮层处理的多尺度研究

本方案展示了如何利用基于物理的神经模拟来解读神经治疗药物的电生理生物标志物，并揭示其对神经回路的影响，从而为神经治疗药物的开发提供一种具有机制基础的方法。...

利用人类新皮层神经求解器揭示神经治疗对脑电图影响的神经机制研究方案

鸡的听觉脑干由负责双耳声音处理的核团组成。单个冠状切片制备可保持完整的神经环路，而培养方法则为在分子、细胞和网络水平上研究神经元结构和听觉功能的发育提供了独特的制备手段。...

鸡听觉脑干切片的制备与培养

学习使用MEG/EEG与MRI联合技术绘制听觉注意相关的皮层动态图谱。该方法通过最小模估计逆成像和共同表面坐标系统，实现高时空分辨率的脑活动映射，为认知神经科学研究提供有效工具。

利用同步脑磁图/脑电图与基于解剖约束的最小范数估计绘制皮层动态：以听觉注意为例

探索更多视频