临清网站建设费用,网络推广简短广告语,网站开发技术实验总结,emlog换wordpress导读
目的#xff1a;fMRI领域的一些研究使用瞬时相位(IP)表征(源自BOLD时间序列的解析表征)考察了脑区之间的同步性。本研究假设来自不同脑区的瞬时振幅(IA)表征可以为脑功能网络提供额外的信息。为此#xff0c;本研究探索了静息态BOLD fMRI信号的这种表征#xff0c;用于…导读
目的fMRI领域的一些研究使用瞬时相位(IP)表征(源自BOLD时间序列的解析表征)考察了脑区之间的同步性。本研究假设来自不同脑区的瞬时振幅(IA)表征可以为脑功能网络提供额外的信息。为此本研究探索了静息态BOLD fMRI信号的这种表征用于获得静息态网络(RSNs)并将其与基于IP表征的RSNs进行比较。
方法从HCP数据集的500名被试中选取100名健康成人(年龄为20-35岁54名女性)的静息态fMRI数据进行研究。使用3T扫描仪进行4次扫描(每次15min)相位编码方向为从左到右(LR)从右到左(RL)。这4次扫描是在两个session中完成的要求被试睁眼注视一个白色十字架。IA和IP表征来自使用希尔伯特变换的窄带滤波BOLD时间序列并使用基于种子的方法来计算大脑中的RSNs。
结果实验结果表明在0.01-0.1Hz频率范围内基于IA表征的RSNs在运动网络的两个session之间具有最高的相似性得分。而对于额顶网络基于IP的激活图在所有频段上的相似性得分最高。对于更高的频段(0.198-0.25Hz)IA和IP表征在两个session中获得的RSNs的一致性都降低了。与基于IP表征的RSNs相比融合IA和IP表征的RSNs使两个session中的默认模式网络之间的相似性得分提高了3-10%。此外同样的比较表明运动网络在频段[0.01-0.04Hz、0.04-0.07Hz、slow-5(0.01-0.027Hz)和slow-4(0.027-0.073Hz)]上提高了15-20%。在功能连接(FC)网络中使用瞬时频率(IF)表征的两个session之间的相似性得分与使用IP表征获得的相似性得分相当。
结论本研究结果表明基于IA表征的方法可以估计RSNs其session之间的可重复性与基于IP表征的方法相当。研究表明IA和IP表征包含BOLD信号的补充信息二者的融合改善了FC的结果。
前言
近十年来功能性磁共振成像(fMRI)作为一种强大的非侵入性成像技术利用血氧水平依赖(BOLD)对比来测量大脑活动和研究大脑功能。使用脑区对之间的BOLD时间序列的时间相关性来估计稳态功能连接(FC)。它通常通过计算跨脑区的fMRI时间序列之间的时间相关性来测量。FC主要是对被试在执行预定义任务(例如手指敲击或视觉任务)时获得的BOLD时间序列进行测量。然而静息态下的脑区在半球内和跨半球之间存在显著的相关性这引起了人们对定义静息态FC的兴趣。静息态下功能连接的区域称为静息态网络(RSNs)并且人们已经探索了各种统计方法来识别它们。基于种子的相关和数据驱动的方法如独立成分分析(ICA)、聚类和主成分分析(PCA)技术是FC分析中常用的方法。在基于种子的相关方法中首先根据假设或基于任务的激活图选择一个预定义的脑区然后将该区域的时间序列与大脑中所有其他区域的时间序列进行相关性分析以找到与之功能连接的区域网络。然而像PCA或ICA这样的技术则试图根据fMRI数据的内在结构来推断大脑活动。
人脑是一个复杂的非线性动态系统为了研究人脑各区域之间的非线性功能关系BOLD信号的解析表征一直是研究的热点。在这种方法中实值BOLD信号时间序列通过窄带滤波器(带通滤波器通过小范围的频率)所得信号被用来获得一个复解析信号。这种复解析信号有两个组成部分瞬时振幅(IA)和瞬时相位(IP)。IA始终为正并包含BOLD信号中呈现的时变幅值信息。另一方面IP信号提供了信号中存在的时变相位信息。在解析信号的这两种表征方法中大多数研究都集中探索以IP表征为主的FC。估计的IP时间序列用于测量两个空间上不同的脑区之间的同步性。这些研究发现在原始时间序列可能在时间上不相关的情况下大脑区域之间存在特定的依赖关系(主要是相位相关)。
Liard及其同事(2002)开发了一种应用于fMRI数据的相位同步(PS)技术以研究体素时间序列与基于事件的手指敲击任务的参考函数之间的相位锁定情况。分析显示在原始时间序列可能不相关的情况下它们的IP之间可能存在相互关系。Zhou等人(2010)探索了与电生理频率相关的不同频段的BOLD信号slow-5(0.01-0.027Hz)slow-4(0.027-0.073Hz)slow-3(0.073-0.198Hz)和slow-2(0.198-0.25Hz)。研究结果表明slow-5和slow-4频段主要与灰质有关而slow-3和slow-2频段主要与白质有关并伴有呼吸和心脏信号混叠。Huotari等人(2019)研究了采样率对静息态FC指标的影响量化了时域、频域、空间域的稳态和动态分析。Glerean及其同事(2012)使用0.04-0.07Hz范围的窄频带进行PS分析以避免被试在自然刺激下出现混叠伪影。该研究定义了几个基于IP的PS指标即基于种子的相位同步(SBPS)、个体间相位同步(IPS)和基于个体间种子点的相位同步(ISBPS)作为动态FC的测量指标以研究复杂的自然刺激和event/blocked设计范式。不同脑区的BOLD时间序列之间的全局相位同步反映了功能同步簇的动态关联和解离。Pedersen等人(2018)发现IP同步分析和基于相关性的滑动窗口方法在时间动态连接方面存在较强的关联。选择适当的窗长以及对fMRI数据进行窄带滤波(0.03-0.07Hz)IP方法和基于相关性的方法都提供了可比较的FC结果。基于小波变换的时频相干分析揭示了后扣带皮层与默认模式网络其他节点之间的相干性和相位变异性。
本研究分析了静息态下BOLD信号的IA和IP表征在检测功能连接脑区中的意义。在基于种子的FC框架下对静息态数据进行了IA分析并观察了RSNs。研究发现使用基于IA表征方法估计的大脑网络在整个session中是一致的并且与基于IP表征的脑网络估计结果相当。使用Jaccard相似性对两个session的RSNs空间图结果的一致性进行了量化。本文还探讨了相位展开的一阶导数即IF作为IP表征的替代方法因为它可以避免与相位计算相关的任何歧义并且可以进一步用于探索静息态的FC模式。此外本研究还探索了用于BOLD信号窄带滤波的不同频段将感兴趣的频段(0.01-0.1Hz)分割成更小的等宽频段0.01-0.04Hz、0.04-0.07Hz和0.07-0.1Hz以观察基于IA、IF和IP表征方法在每个子频段上的RSNs。除了这些频段之外本文还考察了Zhou等人(2010)中提到的低频段其频率范围为0.01-0.25Hz。在整个session中观察到使用不同频段产生的激活图的一致性。对于本研究中探讨的所有频段两个session中基于IA表征的RSNs获得的相似性得分与基于IP表征的方法获得的相似性得分相当。当考虑到所得RSNs在不同频段上的可重复性时对于部分RSNs基于IA的方法表现较好(如运动网络)而对于某些RSNs基于IP表征的方法表现较好(如额顶网络)。将基于IA和IP表征方法的RSNs融合后两个session之间的相似性评分比它们各自的评分都有所提高。
材料和方法
通过希尔伯特变换得到真实信号的解析表征(由Gabor提出)已被广泛用于语音信号处理、脑电图(EEG)和脑磁图(MEG)。然而近二十年来人们对fMRI信号的处理方法进行了大量的研究。解析表征是一种常用的调幅和调频(AM-FM)信号分解方法用于分析信号中呈现的时变振幅和频率内容。对于实值窄带通滤波BOLD信号x(n)利用希尔伯特变换可以得到解析信号xa(n) 其中H[.]为x(n)的希尔伯特变换j为虚数单位a(n)是IAφ(n)是IP。解析信号xa(n)与x(n)具有相同的傅里叶变换且负频率分量为零。根据Bedrosian定理如果带通滤波器的带宽足够窄则x(n)可以根据其IA和IP定义为 振幅a(n)和相位φ(n)分量对于信号的精确重建都很重要。信号参数φ(n)采用arc-tan估计取值范围为-π到π并且存在相位环绕问题。然而相位展开的一阶导数称为IF可用于避免相位环绕问题。窄带解析信号的IF可以解释为局部拟合信号的正弦波频率从而描绘信号中的时变频率内容。IF表征不需要相位信号计算可计算如下 其中Re为表达式的实部xa(k)为xa(n)的离散傅里叶变换F−1表示离散傅里叶逆变换N为信号长度。图1显示了高斯调制模拟信号x(n)的IA、IP和IF表征其频率内容先从0.2到0.03Hz下降(到信号长度的一半)然后再从0.03到0.2Hz增加(至信号末端)。从图中可以看出IA(a(n))给出了模拟信号的包络x(n)。在IP信号中可以看到频率的变化(φ(n))。还应该注意的是相位位于-π至π范围内因此是一个包裹相位。IF信号φ(n)表示信号的频率内容。 图1.模拟信号x(n)的IA、IP和IF表征x(n)是一个高斯调制信号其频率沿时间轴在0.03-0.2Hz范围内变化。IA(a(n))、IP(φ(n))和IF(φ(n))表征分别显示了解析信号(xa(n))中的时变振幅、相位和频率内容。
fMRI数据描述
本研究从人类连接组项目(HCP)数据集的500名被试中选取100名健康成人(年龄为20-35岁54名女性)的静息态fMRI数据进行研究。使用3T扫描仪进行4次扫描(每次15min)相位编码方向为从左到右(LR)从右到左(RL)。这4次扫描是在两个session中完成的要求被试睁眼注视一个白色十字架。采用重复时间(TR)720ms回波时间(TE)33.1ms层厚2.0mm回波平面成像(EPI)加速度因子8采集静息态图像。
分析框架
图2显示了进行FC分析的分析框架。整个流程分为以下三个阶段阶段1信号预处理阶段2信号转换和表征提取阶段3FC分析。阶段1包括fMRI数据的标准预处理然后在所需范围内对BOLD信号进行滤波(频带为fC1-fC2 Hz的Butterworth四阶带通滤波器其中fC1和fC2分别表示带通滤波器的低截止频率和高截止频率)。对于基于BOLD表征的方法期望的频率范围为0.01-0.1Hz这是rs-fMRI FC分析中广泛使用的频率范围。然而对于基于IA、IF和IP表征的方法本研究探索了fC1和fC2的多种选择。窄带滤波后的真实BOLD信号被转换为复解析信号。然后在第2阶段从该解析信号中提取IA、IF和IP表征。 图2.使用BOLD、IA、IP和IF表征对rs-fMRI数据进行FC分析的分析框架。
本研究使用“最小预处理”HCP的静息态数据即优化的空间预处理数据。预处理步骤包括去除各种空间伪影和头动校正。然后将时间序列数据与结构数据和2mm标准(MNI)空间对齐。对整个数据集进行全局强度归一化并去除非脑体素。第3阶段包括FC分析。这里使用基于种子的方法从各自的网络中选择种子点并计算种子表征时间序列与大脑中其余体素表征时间序列之间的线性相关性。然而对于基于IP的连接图相位锁定值(PLV)指标是根据种子的IP时间序列和其余体素的IP时间序列之间的相位差进行计算的。每种方法的激活图是每个session中被试特异性激活图的平均值。
结果
对于IA、IF和IP表征本文将BOLD信号在以下频率范围内进行窄带滤波0.01-0.04Hz0.04-0.07Hz0.07-0.1Hz0.01-0.027Hz0.027-0.073Hz0.073-0.198Hz和0.198-0.25Hz。图3展示了本研究中从横向视觉网络区域提取的两个BOLD时间序列的所有三种表征。种子和目标体素分别用黑色和蓝色突出显示(种子位置用绿色十字准线突出显示)。假设如果两个区域同步则它们的IP时间序列具有恒定的相位差。PLV指标在EEG/MEG数据中被广泛用于评估两个IP时间序列之间的PS如下所示 其中为n个样本上的平均值Δφ(n)为相位差时间序列。PLV0(表示两个时间序列不同步)PLV1(表示两个时间序列高度同步)。 图3.从外侧视觉区域的两个体素中提取滤波后的BOLD(0.04-0.07Hz)信号得到BOLD、IA、IF和IP表征。
Session 1中的前默认模式网络(aDMN)、后默认模式网络(pDMN)、背侧注意网络(DAN)、运动网络、外侧视觉网络(Lat.Visual)和额顶网络(FP)的激活图如(图4-5)所示。采用基于种子的FC方法获得基于BOLD、IA和IF表征的RSNs其中种子位置(MNI坐标)如下aDMN(244−10)pDMN(6−5426)DAN(60−2412)FP(−46−5046)Lat.Visual(31−942)和运动网络(30−3868)。对于IP表征使用PLV指标作为空间上不同脑区IP时间序列之间同步的度量。对于使用Fisher’s z变换FC值的所有表征都显示了显著的功能连接空间图(p0.05)。对于每个session包含两次run从左到右(LR)、从右到左(RL)对于每个被试分别计算每次run的FC矩阵然后计算两次run的平均FC矩阵以便进一步分析。生成的激活图(如图4和图5所示)是每个session中个体被试激活图的平均值。可以观察到基于BOLD信号的转换表征基于IAIP和IF表征的方法能够复现与基于BOLD表征方法获得的RSNs相似的静息态大脑FC模式但有些差异不显著。将这些基于表征方法的激活图(图6-7)成对地融合说明了此处所示的所有RSNs空间图的相似性。 图4.在BOLD、IA、IF和IP表征上使用基于种子的FC方法对session 1的外侧视觉网络、aDMN和FP网络的矢状面、冠状面和轴位面视图进行分析(此处BOLD信号用fC10.04Hz和fC20.07Hz进行窄带滤波并使用此滤波信号计算解析信号)。十字准线位于每个激活图的种子位置颜色条表示z分数值。 图5.在BOLD、IA、IF和IP表征上使用基于种子的FC方法对session 1的DAN、pDMN和运动网络的矢状面、冠状面和轴位面视图进行分析(此处BOLD信号用fC10.04Hz和fC20.07Hz进行窄带滤波并使用此滤波信号计算解析信号)。十字准线位于每个激活图的种子位置颜色条表示z分数值。 图6.基于BOLD-IP、BOLD-IA和BOLD-IF表征方法的RSNs激活图重叠。这里显示了两种激活图的加性融合为了可视化蓝色用于第一种表征方法红色用于第二种表征方法黄色用于两种表征方法的结合。例如在BOLD-IP重叠图中仅在BOLD激活图中激活的体素用蓝色表示仅在基于IP的激活图中激活的体素用红色表示而在BOLD和基于IP的方法中都激活的体素用黄色表示。 图7.基于IP-IA和IA-IF表征方法的RSNs激活图重叠。
本研究探索了截止频率的多种选择(fC1fC2)来对BOLD信号进行滤波(用于提取IA、IF和IP表征)。为了探索频率对估计FC的影响本研究还分析了比文献建议的PS分析频段(0.04-0.07Hz)更低和更高的频段(0.01-0.04Hz和0.07-0.1Hz)以涵盖传统静息态研究使用的整个低频范围(0.01-0.1Hz)。除了这些频段本文还探索了Zhou等人(2010)所描述的低频子频段这些子频段涵盖了与大脑灰质和白质组织相关的频率。图8显示了本研究中使用的基于IA、IF和IP表征的方法在不同频段获得的两个RSN(外侧视觉和aDMN)的加性融合。可以观察到所有三种表征的种子共激活区域在所有频段中占主导地位因此表明这些表征获得的RSNs具有相似性。 图8.基于IA、IF和IP表征的外侧视觉网络和aDMN网络在不同频段上的重叠。
本研究观察了每个频段的RSNs并且在两个静息态fMRI session之间计算这些表征用于估计FC的可靠性。利用Jaccard相似性(JS)得分计算两个session空间图之间的相似性作为可靠性的衡量指标。图9显示了基于IA、IP和IF表征的所有频段(0.01-0.04Hz、0.04-0.07Hz、0.07-0.1Hz、0.01-0.027Hz、0.027-0.073Hz、0.073-0.198Hz和0.198-0.25Hz)的RSN在两个session之间的JS分数。 图9.对不同频段的fMRI信号进行(IA、IP、IF和BOLD)表征得到两个session的RSNs之间的JS分数。
对于本研究中的所有RSNs基于IA表征的方法与基于IP表征的方法(PLV)具有相当的相似性得分。因此它可以被认为是计算FC的可靠度量。基于IF表征得到的估计FC是有噪声的(由一些在感兴趣网络区域之外活动的体素组成)因此是所有表征方法中最不可靠的。
BOLD信号的IA、IF和IP表征中的补充信息通过融合每个表征在特定频段(0.04-0.07Hz)所获得的RSNs(平均)来说明。融合IA和IP表征的RSNs使得默认模式网络和运动网络在两个session之间的JS得分比基于IP表征的得分分别提高了3-10%和15-20%。值得注意的是IF和IP表征的RSNs没有融合因为两者提供的都是相位信息(相位和相位的导数)。对于基于IA-IP融合的RSNs两个session在频段(0.01-0.04Hz0.04-0.07Hzslow-5和slow-4)上的相似性得分有所提高(图10)。但对于其余频段任何RSNs均未见显著提高。 图10.通过组合不同频段的fMRI信号(IA-IP和IA-IF)表征获得的两个session RSNs之间的Jaccard相似性得分。
结论
本研究使用解析信号的振幅(IA)、相位(IP)和频率(IF)表征对100名被试的静息态fMRI数据进行了FC图分析。使用基于种子的方法来演示RSNs该方法使用最为广泛可用于寻找与所选种子点存在功能连接的区域。使用这些基于表征的方法得到了aDMN、pDMN、DAN、FP网络、外侧视觉网络和运动网络的FC激活图。这些基于表征的方法(IP-IA和IP-IF)在空间图上的相似性表明所有表征方法都可以很好地探索静息态下的FC。通过JS分数计算的基于IA、IF和IP表征方法(对于在两个session中获得的RSNs)的可靠性表明基于IA表征的FC图在session之间具有良好的相似性得分与基于IP表征的方法相当。本研究发现使用IF表征获得的RSNs低于使用IA和IP表征方法获得的RSNs但并不差。这可能是因为BOLD信号是一个缓慢变化的信号因此由IF表征的展开相位梯度不能增强BOLD信号对共激活体素的判别因子。
fMRI FC的一致性已使用组内相关系数(ICC)指标进行了广泛的评估。ICC值越高表示可靠性越高而ICC值越低表示可靠性越差。一些研究表明高连接强度的连接具有高可靠性并且session之间的ICC分数通常在0.1-0.75之间。与上述研究相比本研究的相似性度量即Jaccard系数值更高范围从0.35-0.93。其原因可能是本研究考虑了各种窄带频率用于从BOLD时间序列中提取不同的表征(IA、IP和IF)并计算使用相同方法获得的两个session之间的功能连接图的相似性这可能会导致比之前的方法使用宽带BOLD时间序列获得的相似性分数更高。使用rs-fMRI信号的实验结果表明基于IA和IP表征的成功取决于所采用的滤波器的范围和带宽。理想情况下滤波后的信号(解析信号)应该具有单频分量以获得有意义的表征。但应该注意的是由于BOLD信号是缓慢变化的信号(0.01-0.25Hz)因此在低频范围(0.04-0.07Hz)中使用可变窄带滤波器将不会捕获IA或IP信号中的重要细节而IA或IP信号负责测量两个区域之间的同步。此外在高频范围内使用滤波频段也可能无法从缓慢变化的BOLD信号中提供有用的表征。这可能就是使用未滤波的BOLD信号获得的JS分数比使用IA、IF和IP表征获得的分数更好或略差的原因。
原文P. Mittal, A.K. Sao and B. Biswal, Impact of amplitude and phase of fMRI time series for functional connectivity analysis, Magnetic Resonance Imaging (2023), https://doi.org/10.1016/j.mri.2023.04.002
