血管迂曲结构对伪连续动脉自旋标记的影响及流驱动绝热翻转的修正作用
2021-10-16袁伟文江桂华罗泽斌
袁伟文,江桂华,罗泽斌
1.广东医科大学第一临床医学院,广东湛江524023;2.广东省第二人民医院影像科,广东广州510020
前言
动脉自旋标记(Arterial Spin Labeling,ASL)使用绝热脉冲对动脉血水分子的氢(H)原子核净磁场(Net Magnetization,M)进行方向上的翻转[1]。理想状态中,净磁场M翻转180°后,标记后的水分子携带异于人体其他组织的负性磁共振信号-M,成为血液的内源性示踪剂,当标记水分子随着血液流入相应组织,将特异性降低该组织的整体信号强度。2015年Alsop 等[2]起草的ASL 临床使用指南将伪连续动脉自选标记(Pseudo-Continuous Arterial Spin Labeling, PCASL)列为众方案的推荐首选。既往文献对PCASL可靠性、适用性及可重复性进行了验证,但仅停留在受检者脑血流量(Cerebral Blood Flow,CBF)测量结果的比较层面,并没有总结其技术原理层面的原因[3-4]。在PCASL 的实际操作中,多种因素可能影响其测量结果,标记血管的迂曲结构被认为是其中之一[5]。本实验通过比较129 名受检者标记层面与血管迂曲结构不同几何关系下PCASL的测量差异性来探讨血管迂曲结构对PCASL 的实际影响,并通过流驱动绝热翻转的原理特征来解释产生该现象的原因,为PCASL的临床操作提供参考。
1 资料与方法
1.1 一般资料
随机选取2018~2019年进行头部磁共振灌注成像的受检者140 名,排除无法完整进行扫描卡片所有序列及因幽闭恐惧症、躁动等原因导致图像出现较明显伪影的个体后,最终纳入受检者129 名,平均年龄(44±7)岁,根据默认方案PCASL 定位时标记层面与血管迂曲结构的几何关系对其进行分组。控制组:标记层面交于双侧颈内动脉、椎动脉走行笔直处,共34名(男15名,女19名)。实验组:标记层面交于双侧颈内动脉或(与)椎动脉血管迂曲处,共95 名(男46名,女49名)。
1.2 仪器设备
采用荷兰Philips公司Ingenia 3.0T磁共振扫描仪及配套32通道阵控式头部线圈。
1.3 方法
1.3.1 扫描卡片设计磁共振扫描卡片包括快速血管成像序列(PCA)、血流速度检测序列(Qflow)、两次默认方案PCASL 序列、两次修改方案PCASL 序列、弥散加权成像(DWI)序列。PCASL 两次默认方案扫描用于验证序列本身的稳健性与可重复性;两次修改方案分别将标记层面与成像区域中心距离d改为95 与100 mm,以避让血管迂曲结构;PCA 用于确定受检者血管迂曲结构与标记层面的几何关系,对受检者进行分组;DWI 用于检查受检者是否存在可能导致PCASL 的CBF 测量结果异常的颅内病变;Qflow 用于监控受检者扫描过程中血流速度有无明显异常改变。
1.3.2 PCASL 数据处理PCASL 的测量结果以CBF图呈现,由磁共振机器配套系统生成。受检者个体CBF 值取成像区域中间层面脑灰质随机10个圆形感兴趣区域的平均值。感兴趣区域面积均为10 mm2,尽量避开血管、边缘部位。
1.4 统计学方法
控制组与实验组均以第一次PCASL默认方案为基准,比较第二次默认方案、第一次修改方案、第二次修改方案与第一次默认方案测量结果的差异性。测量结果的统计学差异由Bland-Altman检验与配对t检验进行判定,P<0.05 为差异有统计学意义;通过山形图与盒形图对数据的总体分布状况与趋势进行补充说明。
2 结果
PCASL 序列默认方案定位时标记层面交于血管笔直部位者纳入控制组,共34 名。控制组第一次默认方案与第二次默认方案、第一次修改方案的CBF测量结果均无统计学差异(P>0.05),第一次默认方案与第二次修改方案的CBF 测量结果具有统计学差异(P<0.05)。
PCASL序列默认方案定位时标记层面交于血管迂曲结构者纳入实验组,共95名,其中17名血管迂曲结构多次穿行标记层面。实验组第一次默认方案与第二次默认方案、第一次修改方案的CBF测量结果均无统计学差异(P>0.05),第一次默认方案与第二次修改方案的CBF测量结果具有统计学差异(P<0.05)。
控制组与实验组的受检者均未在DWI序列中发现可导致颅内血流发生明显改变的病变,Qflow 序列作为佐证同样未发现有受检者在数次PCASL扫描中出现显著血液流速差异。
控制组与实验组的组内比较以山形图的形式表示(图1)。山形图的横坐标为默认方案一与后3 次PCASL 扫描方案的CBF 测量结果差值(CBF默认方案一-CBF后续方案),纵坐标为差值百分位数。两次扫描方案的数据越接近,山形图越趋于以零点为中心的对称图形。山形图向负轴偏移,则表示后续扫描方案所测CBF 高于作为参照标准的默认方案一;山形图向正轴偏移,则表示后续扫描方案所测CBF 低于默认方案一。
图1 控制组与实验组后3次扫描与默认方案一的CBF测量结果山形图Fig.1 Mountain plots of CBF measurement results between the last 3 protocols with the first default protocol in control group and experimental group
无论是控制组或实验组,默认方案一与默认方案二的山形图基本以零点对称分布,显示两次CBF测量结果无明显区别。两组的默认方案一与修改方案一的山形图均稍向正轴偏移,显示修改方案一CBF 测量结果整体稍低于默认方案一,但区别不大。两组的默认方案一与修改方案二的山形图同样继续向正轴偏移,以实验组偏移更明显,显示修改方案二CBF 扫描结果较默认方案一明显下降,且实验组下降尤甚。
盒形图显示修改方案一所测CBF 数据整体上与两次默认方案的数据分布区别不大,修改方案二所测结果整体呈下降趋势(图2a、2b),实验组的下降幅度大于控制组(图2c)。
图2 控制组与实验组各次PCASL测量结果的盒形图Fig.2 Box plots of PCASL measurement results in control group and experimental group
3 讨论
3.1 流驱动绝热翻转的实现原理
PCASL 的标记通过流驱动绝热翻转来进行[6]。在标记层面中设置一个双极梯度场,其中心与标记层面中心重叠。梯度场在标记层面最下沿有最大正值Gmax,此后线性衰减,在标记层面中心衰减为0,越过标记层面中心后负性增长,在标记层面最上沿有最大负值-Gmax。血液在流入标记层面前,仅受到z 轴方向的B0作用,其净磁场M的进动频率ω为:
式(1)为经典的拉莫尔进动公式,γ为常数磁旋比。当血液流入标记层面后,除了外磁场B0,还受到梯度场的作用。结合血液流动性与梯度场的极性,血液水分子所受梯度场强度G(t)为时间t的因变量:
其中,Gave为梯度场平均强度;z为标记层面宽度;v为血液垂直流速;t为血液在标记层面中移动时间。结合式(1)与式(2),血液水分子标记层面内垂直流动时,其净磁场M的进动频率Ω为:
用选择坐标系替换传统空间坐标系,在x轴施加射频脉冲B1,并设定xy平面绕z轴旋转的频率fxy与射频脉冲B1频率f、血液未流入标记层面净磁场M进动频率ω相等,即:
净磁场M在旋转坐标系中的进动频率变为Ω与f的差值,根据式(1)与式(3),有:
将拉莫尔进动公式代入式(5),旋转坐标系中净磁场M受到的外磁场从B0改为Bz,且:
PCASL 的脉冲间隔非常短,旋转坐标系中的净磁场M可以视为持续受到Bz与B1两个磁场的作用,遵循矢量法则合成综合磁场Beff
射频脉冲B1为固定值,因此Beff大小与方向变化取决于Bz(G(t))的变化,并分为以下5 步:(1)血流进入标记层面最下沿,G(t)=Gmax≫B1,Beff的合成由G(t)主导,其方向基本与G(t)一致,指向坐标系z轴;(2)血流从标记层面最下沿流向标记层面中心,G(t)强度不断衰减,B1逐渐占据主导地位,Beff方向向B1所在x 轴倾斜;(3)血流到达标记层面中心时,G(t)衰减为0。Beff=B1,方向指向x 轴;(4)血流从标记层面中心流向标记层面最上沿,进入梯度场负极。G(t)方向指向-z 轴并负性增长,|G(t)|逐渐大于B1,G(t)重新占据主导地位。Beff方向从x 轴逐渐向G(t)所在-z 轴偏转;(5)血流到达标记层面最上沿,|G(t)|=|-Gmax|≫B1,Beff的方向基本与|G(t)|主导一致,指向-z轴。
上述5 步即综合磁场Beff的翻转过程,净磁场M绕Beff进动,因此随其一并翻转,Beff与z轴夹角等同于净磁场M的翻转角α。翻转角α以血液在标记层面中移动时间t为自变量的函数表达式为:
式(8)与De Graaf 等[7]描述的绝热翻转一致,区别在于传统绝热翻转需要不断调整射频脉冲频率以达成Bz的变化[8],而流驱动绝热翻转则利用了血液流动性与梯度场极性的结合[9]。
3.2 流驱动绝热翻转的物理特性
3.2.1 流驱动绝热翻转的位移依赖性在式(8)中,血液在标记层面内的运动时间t与血液垂直流速v以乘积的方式出现,表示流驱动绝热翻转效果并非直接由血液流速或射频脉冲作用时间决定,而由血液垂直位移决定,即位移依赖性。
绝热翻转的位移依赖性质最大程度削减了血流差异对净磁场M翻转的影响,但并非流速对绝热翻转毫无相干。血液流速通过结合其它因素影响绝热翻转的最终效果。净磁场M会在绝热翻转进行的同时受到纵向弛豫与横向弛豫的影响,当血液流速过慢时,弛豫影响增大,翻转效率下降。磁场M在xz平面的翻转相当于绕y 轴进行的进动,因此y 轴将产生Bz与B1以外的第三个磁场By,当血液流速过快时By较大,其参与Beff的合成,导致翻转效率下降。因此,绝热翻转对血液流速具有一定的适用条件,Maccotta等[10]将其命名为流驱动绝热翻转的绝热条件,数学表达式如下:
式(9)表明了流驱动绝热翻转的最适流速范围。将T2blood(3.0T 场强下为147 ms)、Gave(0.36 mT/m)与B1(1.814 3 μT)代入式(8)与式(9),可知血流速度在3.6~38.3 cm/s 范围内变动时,PCASL 的平均标记效率在89%左右浮动,显示出良好的稳定性。Václavů等[11]在2018年ISMRM 会议所作报告中以同样的参数获得了相同的结果(图3)。
图3 流驱动绝热翻转在适合流速范围内可稳定保持高标记效率Fig.3 Flow-driven adiabatic inversion maintains high labeling efficiency within a suitable flow rate range
3.2.2 流驱动绝热翻转的绝热性绝热在量子物理学中定义为运动系统能量状态的缓慢改变,但该改变只能逼近而无法达到另一能级(图4a)[12]。式(8)属于反正切函数,其在定义域内单调递增,然而值域为开口状态。该数学意义的实践表达为:无论如何调节PCASL 序列参数,受检者标记血管内的动脉血流经标记层面时,血液水分子净磁场M的翻转角α最多只能接近180°,无法像共振翻转一样,达到180°甚至超过180°。梯度场的负极设置在标记层面最上沿,确保血液水分子在任何位置进入标记层面,最终到达标记层面最上沿时,其在旋转坐标系中所受综合磁场Beff均指向-z方向。以上两点完美契合了绝热的物理学概念。Hernandez-Garcia 等[13]在不同序列参数条件下测试了流驱动绝热翻转的绝热性(图4b),净磁场逼近但无法达到-M,因此不存在重复、过度翻转,值得注意的是时间轴越过零点后曲线上扬代表净磁场的纵向弛豫。
图4 绝热理论的量子物理学表述及其在PCASL中的应用Fig.4 Adiabatic theorem in quantum mechanics and its application in PCASL
3.3 血管迂曲结构对PCASL 的影响及流驱动绝热翻转对其的修正
血管迂曲结构可分为先天性与后天性,其在颈内动脉颅外段的发生率约为4%~66%[14],考虑到枢椎层面的先天性弯曲,其在椎动脉颅外段的发生率为100%[15]。ASL临床使用指南将血管迂曲结构列为可能影响PCASL 的CBF 测量结果的因素,但没有展开进一步的讨论。当血管迂曲结构交于标记层面时,有完全位于标记层面内与多次穿行于标记层面两种情况。
位于标记层面内的血管迂曲结构最直接的影响为改变血流的路径,相当于降低血液的垂直流速,增加血液水分子在标记层面中的运动时间及其所受射频脉冲数量。在共振翻转中,翻转角的大小为磁旋比、射频脉冲场强与时间三者的乘积[16],因此血液流速的改变直接影响翻转效果。在流驱动绝热翻转中,其位移依赖性可修正血液流速变化对翻转效果的影响。实验组默认方案与第一次修改方案的测量结果判定无统计学差异,验证了流驱动绝热翻转的位移依赖特性。
多次穿行于标记层面的血管迂曲结构可将标记后的血液重新导入标记层面,使血液水分子再次接受标记脉冲的作用。Amukotuwa 等[17]曾报道某受检者PCASL 序列与其它磁共振序列结果相悖的案例,并认为出现CBF 测量异常的原因在于该受检者左侧血管迂曲结构导致的重复翻转。根据前文数学模型论证,流驱动绝热翻转的绝热性确保了针对血液水分子净磁场M的翻转不会真正达到或超越180°,因此PCASL不存在所谓的重复翻转。实验组中共有17名受检者的血管迂曲结构重复穿行于标记层面,无一出现异常CBF 测量结果,进一步验证绝热特性对PCASL标记稳健性的保障作用。
标记后的血液水分子净磁场为-M,与脑组织水分子净磁场方向相反,因此具有降低灌流区域脑组织磁共振信号强度的示踪剂效应,但纵向弛豫会逐渐将其恢复为与脑组织水分子净磁场方向相同的M,降低其对脑组织信号的干扰能力,Zuo等[18]将纵向弛豫对内源性示踪剂的损耗作用命名为T1效应。位于标记层面与成像区域间的血管迂曲结构最直接的影响是延长示踪剂的运输路程,增加示踪剂的运输时间,其将放大T1 效应,导致PCASL 的测量结果下降。T1效应对标记血液的损耗发生在其完成绝热翻转离开标记层面以后,因此流驱动绝热翻转对其无修正作用。本研究中,实验组与控制组的第二次修改方案的测量结果均较第一次默认方案下降,且实验组下降趋势较控制组明显,验证了该影响。Alsaedi等[19]在PASL序列的实验中观察到相似的表现。
4 结论
流驱动绝热翻转通过血液流动性、梯度场极性与射频脉冲链的巧妙结合而使PCASL的标记具有独特的血液位移依赖性及绝热性,使PCASL 序列能广泛、稳健地运行在不同血流速度条件下的受检者身上,并且可以防止重复、过度翻转的出现。位于标记层面内部与多次穿行标记层面的血管迂曲结构主要改变受检者血流速度与导致标记后血液重复进入标记层面,以上两点均可为流驱动绝热翻转的数理特性所修正,一般情况下对PCASL 的测量结果无实际影响。在实际操作中,PCASL 与常规序列出现较明显相悖状况时,不应简单归咎为血管迂曲结构所致失误,宜重复进行PCASL 或进一步检查。位于标记层面与成像区域间的血管迂曲结构可因T1效应而导致PCASL 的CBF 测量结果下降,且流驱动绝热翻转对其无修正作用。在实际操作中,标记层面距离不宜过度增加,以防止血管迂曲结构进一步放大T1效应。