基于功能磁共振流入效应量化脑脊液运动的心脏、呼吸及低频成分

《MAGNETIC RESONANCE IN MEDICINE》：Quantifying Cardiac, Respiratory, and Low Frequency Components of CSF Motion From fMRI Inflow Effects

　　

目的：脑脊液(CSF)流动振荡被认为是反映脑清除功能的重要标志物，但其测量通常依赖专用磁共振成像序列。本研究旨在利用广泛使用的功能磁共振成像(fMRI)数据，定量评估与心脏、呼吸及低频周期相关的CSF流动。方法：研究人员开发了一种将fMRI衍生的CSF流入信号转化为定量流速的方法，该模型描述了振荡分子群的纵向磁化历史。通过在模拟心脏、呼吸及低频振荡的体模实验中进行验证，并将该方法应用于48名老年人(68–82岁，女性19人)的静息态数据集，以表征枕骨大孔水平的CSF流动。结果：体模实验显示，估计流速与真实流速在心脏及呼吸类似频率下相关性良好(r=0.94和0.97)，在低频类似振荡中呈中等相关(r=0.58)。在人群队列中，心脏周期的中位CSF每搏输出量为0.77[0.57,1.09]mL，呼吸周期为0.38[0.26,0.88]mL，低频周期为0.26[0.14,0.39]mL。结论：所提出的纵向磁化历史建模方法可利用常规fMRI数据集定量估计CSF流动成分，并显示心脏周期主导枕骨大孔处的CSF运动。

该研究由瑞典研究团队发表于《MAGNETIC RESONANCE IN MEDICINE》，旨在解决当前脑脊液(CSF)动力学定量评估对专用磁共振脉冲序列的依赖问题，并利用广泛采集的功能磁共振成像(fMRI)数据实现多生理周期的CSF流动分离与量化。CSF流动受心脏搏动、呼吸压力变化及低频血管振荡驱动，与脑清除功能（类淋巴系统）及神经退行性疾病密切相关。然而，传统相位对比磁共振成像(PC?MRI)虽可定量CSF流速，却难以在大样本回顾性分析中普及。fMRI的边缘层面可观测CSF流入效应，但因信号受纵向磁化历史影响且采样率常低于心脏频率的奈奎斯特极限，导致混叠与定量困难。研究人员提出一种基于Bloch方程的正向模型反演方法，将fMRI流入信号转化为定量流速，并通过体模及人群数据验证其可靠性。结果显示，心脏周期在枕骨大孔水平CSF运动中占主导地位，且该方法兼容标准fMRI协议，无需额外扫描。

关键技术方法方面，研究纳入瑞典乌默奥人口登记的61名老年人，最终分析48人(68–82岁，女性19人)，均完成静息态fMRI及生理记录。MRI采用3T GE Discovery MR750扫描仪，BOLD敏感梯度回波平面成像(EPI)序列(TR/TE=2000/30ms，翻转角80°)，同时采集脉搏血氧与呼吸波纹信号。CSF感兴趣区(ROI)手动绘制于fMRI最低层面，利用K?means聚类分离CSF与非CSF体素。模型假设CSF速度由心脏、呼吸及低频成分叠加构成，利用希尔伯特变换提取各生理信号的瞬时相位，结合纵向磁化历史模拟与粒子群优化反演流速参数。体模实验采用空心圆柱模拟振荡流，数字体模用于评估空间分辨率、噪声及几何假设的稳健性。

研究结果如下：

3.1 体模验证

物理体模实验中，心脏频率(0.93Hz)、呼吸频率(0.23Hz)及低频(0.077Hz)的估计与真实流速显著相关，其中低频在≤1.5mm/s范围内相关性达r=0.99。数字体模证实模型在不同空间分辨率、噪声水平及非理想几何条件下仍保持高精度，且塞状流假设拟合的流速与层流、环空流的平均流速高度一致(r=0.999)。

3.2 CSF流速与容量

在43名受试者中，模型拟合优度R²中位数为0.57，显著高于随机配对对照。心脏、呼吸及低频的速度幅值中位数分别为12.0[8.1,18.4]mm/s、1.21[0.76,2.14]mm/s、0.18[0.11,0.24]mm/s；对应每搏输出量为0.77[0.57,1.09]mL、0.38[0.26,0.88]mL、0.26[0.14,0.39]mL。固定组平均T₁值与个体T₁值估算结果差异小于3%，且无外部门控的低频估算与全模型高度相关(r=0.83)。

讨论部分指出，该方法首次在标准fMRI框架下实现多生理周期CSF流动的定量分解，避免了PC?MRI的专用序列限制，并与既往PC?MRI报告的枕骨大孔流速及每搏输出量一致。低频CSF流速缺乏参考值，但基于颅内压慢波推算的估计值与本结果同量级。局限性包括头部运动干扰、B₀漂移、塞状流假设及未纳入净CSF产生流速，但这些因素对振荡流速影响有限。结论强调，通过纵向磁化历史建模，可从常规fMRI数据可靠量化心脏、呼吸及低频驱动的CSF运动，心脏周期在枕骨大孔水平占主导，为大规模回顾性分析脑清除与神经退行疾病关联提供了可行工具。