摘 要:化学交换饱和转移磁共振成像通过施加特定频率饱和脉冲来标记包含某类可交换质子的分子,而被标记质子通过化学交换引起水信号的变化实现间接测量该分子及其微环境.与传统氢谱(1H- MRS)相比,化学交换饱和转移磁共振成像可达到更高的灵敏度和空间分辨率.然而,对于活体成像,水信号的降低不光来自 化学交换饱和转移效应,还来自直接水饱和、半固体磁化转移效应和核奥氏效应等竞争效应.要实现精确可靠 的化学交换饱和转移磁共振成像量化,就必须减弱饱和脉冲对其他竞争效应的影响.介绍了化学交换饱和转 移磁共振成像定量技术的基本原理和最新进展,并阐述各种量化方法的优缺点,供对本领域感兴趣者参考及在此基础上进一步地研究开发.
关键词:化学交换饱和转移;磁共振成像;量化 CEST;数学模型
1 化学交换饱和转移磁共振成像(Chemi- calExchangeSaturationTransfer,CEST- MRI)量化的意义
CEST-MRI成像是目前备受关注的一种分子影像技术,其使用特定频率的饱和脉冲来标记待检测的溶质分子,该分子上可交换质子与水发生多次化学交换引起水信号的降低从而被检测[1].尽管生物组织中这些小溶质池的浓度通常仅在毫摩尔浓度范围内,只要选择合适的实验参数,可交换质子与水的化学交换的累积效应就能实现信号放大, 因此 CEST- MRI具有较高的灵敏度(与传统的MRI相比,可达到102~106)[2].自从 Ward和 Bal-aba 等人 2000 年提出 CEST - MRI 的 概 念 以来[1],由于其具有可灵敏检测特定(类型)分子的独特优势,且对温度和 pH 值[3]等微环境的变化敏感,被认为有临床转换潜能的可包含生化信息的对比机制,迅速引起广泛的兴趣[4].与正常组织相比,肿瘤本身显示共振频率距离水~3.5 的高信号,被认为主要是来自多肽和蛋白质的酰胺质子,称之为酰胺质子转移 (amideprotontransfer,APT)成像[5-6].CEST- MRI 技术也可检测蛋白质和多 肽[7-8]、肌酸[9]、葡萄糖[10]、谷氨酸[11]、糖原[12]和糖胺聚糖[13]等,有潜力在临床中用于各种疾病和代谢性紊乱,包括:乳腺癌[14-15]、前列腺癌[16]和中风[17]等的检测.该领域已经从最初简单的 CEST加权图像转变为更量化的 CEST-MRI分析.为了进一步揭示与疾病相关的病理生理特征,并对同一 组织不同时间点进行监测,需要更加精准的量化. 因此,笔者简要地概述 CEST- MRI的不同方法和最新进展,供对该领域有兴趣的人员参考.
CEST-MRI的基本原理是通过水信号的减少来间接实现对特定分子的检测.然而,对于活体CEST-MRI成像,水信号的降低不光来自 CEST效应,还来自直接水饱和度(directwatersatura- tion,DS),半 固体磁化转移效应 (magnetization transfercontrast,MTC)和 核奥氏效应 (nuclear overhauserenhancement,NOE)等 竞争效应[18].同时,由于这些竞争效应受到静磁场 B0、饱和功率(saturationpower,B 1-sat)和其他实验参数的 影响,致使 CEST- MRI成为一个复杂的技术.如 何从水信号中提取出 CEST 效应,尽量减弱饱和脉冲对其他效应的影响,实现对特定分子的精确量化,进 而更加准确地 对疾病进行诊 断,一 直 是CEST-MRI研究的热点.为了证明特定的 CEST效应,通常通过获取大量的水信号强度作为射频(RF)饱和频率偏移的函数来产生Z 谱.Z 谱是通过施加脉冲后水信号的强度Ssat 与施加脉冲前水信号的强度S0 的比值得到 .理论上,Z 谱中显示的CEST 增强取决于池的大小、交换率和可交换质子的弛豫时间.这个简单的概念与其他相互竞争的影响导致对 CEST 量化变得复杂.因此,笔者拟定从量化的角度对 CEST-MRI作一个回顾.
2 CEST-MRI理论模型
CEST- MRI 成 像 机 制 复 杂,不 仅 取 决 于CEST-MRI试剂的浓度、交换和弛豫性质,而且随B0 和B 1-sat等实验条件的变化而变化.因此,对于研究这些最优条件,数值模拟是有用和有效的. CEST 效应的量化是一个复杂的过程,除了与质子浓度和交换率有关之外,还受到如 RF 辐射强度、持续时间、化学位移、主磁场强度以及纯水的横向、纵向弛豫时间等因素的影响[20].为了更精确地描述 CEST 成像机制,通常使用包含可交换质子池的布洛赫方程(Bloch- McConnellequation).2004年Zhou等人将包含可交换质子池的最简化的 2池模型的Bloch方程引入到 CEST-MRI中,并推导出解析解[21].常见的Bloch方程描述的模型包括 2池模型[20](水池和可交换质子池)和3 池模式[22](水池、可交换质子池和大分子固体池).下面以2池模型为例来介绍 Bloch 方程的量化方法,同时,给出3池模型的框图和仿真结果,并说明磁化转移(magnetizationtransfer,MT)对 CEST 量化的影响.2池模型和3 池模型的 CEST 交换示意图如图1所示.
在2 池模型中,描述2 池交换过程用12 个参数描述.其中:T1i(i=a,b)表示i 池的纵向弛豫时间;T2i(i=a,b)表示i 池的横向弛豫时间;Ma 表示水池核自旋的化学位移;Mb 表示可交换质子池核自旋的化学位移;Moa 表示水池的可交换氢原子核在初始条件下的平衡浓度;γM0b 表示可交换质子池的可交换氢原子核在初始条件下的平衡浓度与水池之比;kba 表示可交换质子池核自旋交换到水池的速率;kab 表示水池的核自旋交换到可交换质子池的速率;ω1 表示描述预饱和射频脉冲的强度;t 表示脉冲施加时间(下面公式中kab和kba 的意义相同).
在核磁矩和磁化强度矢量的概念基础上,通过在旋转坐标系下,MRI系统中B1-sat 是从x 方向上加入的和2池模型的交换,得到描述图1A 所示2池模型的Bloch方程如公式1所示.
其中,描述b和a 分别表示可交换质子池和水池.例如,Mb (t)表示在时间t 时,可交换质子池的磁化在x 轴上的分量;Mb 和 Ma 分别表示在可交换质子池和水池在z 轴上的热平衡磁化;ω 表示 RF 照射的频率;ωb 和ωa 分别表示可交换质子池和水池的拉莫尔频率;Δωb 和Δωa 分别表示ωb -ω 和ωa-ω .
图1模型示意图中,仿真的各参数如表1所示.序列的参数设置如下:Z 的初始值为1;B0 =7 T;2池模型的饱和时间tp =5s,3池模型的饱和时间tp =80s;B1-sat 分别为:1μT,2μT,3μT,4 μT,5μT 和6μT.模型的仿真结果如图2 所示.其中,上面大的光滑的曲线为饱和交换后水的 Z 谱, 左下方比较矮的曲线为 Diff_Z 谱,即饱和交换之前水的Z谱和饱和交换之后水的Z谱的差.
由2池和3池的对比结果可以看出:在 CEST发生化学交换(~2ppm)的地方,2 池模型在仿真 B1_sat 范围内,Z 谱的下降都比较明显,而3 池模型在B 1_sat4μT 时Z 谱的下降已经看不到.因此, MT 成分大大降低了 CEST 的特异性,尤其是在B1_sat 大的时候.但是 B1_sat 小的时候,病变部位和正常组织的差异较小,即:区别病变部位和正常组织需要较大的 B1_sat.因 此,对 于量化 CEST -MRI,如何降低 MT 成分的影响是需要考虑的很重要的一个因素.
Bloch 方程不仅可以对 CEST 实验进行模拟仿 真,而 且 可 以 对 CEST 测 量 值 进 行 数 值 拟合[19,23].此外,用扩展的 Bloch 方程描述的多池模型,考虑到 RF 照射 CEST 的 MT 和 NOE[24]伴随效应.Bloch方程的一个最大的不足就是计算复杂.因此,有各种假设对 Bloch 方程进行简化,如强饱和脉冲的假设和弱饱和脉冲的假设[25].将 CEST数据与用非线性最小二乘法优化的 Bloch 方程拟合可以提供相关 MRI参数的值,这有助于设计特定应用的最优 PARACEST 试剂[19].但是,在实际应用中,由于 Bloch 方程可变参数繁多,直接拟合难以标准化,因此研究者们提出了其他较为简化的量化方法.
3 CEST-MRI量化方法介绍
3.1 简 单 磁 化 转 移 不 对 称 性 (magnetization transfer ratio based on asymmetry analysis, MTRasym )分析及其改进
由于Bloch方程比较复杂,以水频率为中心将Z 谱 两 侧 相 减 的 磁 化 转 移 率 不 对 称 分 析 法(MTRasym )成为一种常用的 CEST- MRI量化方法[21,26],MTR 的定义如下:
其中,S0 是 RF 照射前的水信号的强度,Ssat(Δω)和Ssat(-Δω)分别是经过 RF 照射后在标记频率和参考频率处的信号(后面公式 S0、Ssat(Δω)和 Ssat(-Δω)表示的意义相同).
该量化方法有以下缺点:(1)容易受到 B0 场不均匀性的影响[27].(2)MTR
方法容易受到各种混淆参数的影响.包括组织松弛、MT、MTC[28].
更重要的是 MTRasym 不能从 CEST 对比中区别上场的 NOE 效应.(3)MTRasym 方法没有校正水的纵向弛豫效应,这是影响 CEST 信号幅度的主要因素[29].尽管如此,MTR 仍是目前 CEST- MRI领域一 种 简 单 有 效 的 量 化 CEST - MRI 的 方法[30].
针对 MTRasym 量化 CEST-MRI存在的不足,提出了一系列的改进措施.例如:针对 MTRasym 容易受B0 场不均匀性影响,提出了水饱和移位参考(watersaturationshiftreference,WASSR)[31]和同时校正B0 场和B1 场的不均匀性(Simultaneous mappingofwatershiftandB ,WASABI)[32]两种方法.针对减弱 CEST 的竞争效应,提出了三频偏方法[2,33]和洛伦兹拟合[34-37]等量化方法.针对没有校正水的纵向弛豫效应的问题,提出了表观交换依赖弛豫(apparentexchange-dependentrelaxa- tion,AREX)[22,29,38-39]的量化方法.
3.2 三频偏的方法
三频偏的方法通过用特定共振频率照射如酰胺、胺基质子等,用和共振频率两个相近频率的平均值作为参考信号,用特定频率处的信号作为标记信号的一种量化方法[2,33].三频偏方法被证明可以有效减弱 MT 不对称性和 NOE 效应[29].三频偏及其量化 APT 和 NOE 的定义如下[33]:
其中,APT* 表示使用三频偏方法得到的酰胺质子转移.NOE* 表示使用三频偏方法得到的 NOE效应.
三频偏方法可以减弱 CEST 的一些竞争效应,例如:MT 不对称性和 NOEs.该量化方法已经成功的用于肿瘤[2]和中风[33]的检测.但是三频偏量化的方法存在如下的缺点:(1)三频偏方法线性假设过于简单,明显低估了 APT 和rNOE 在~3.5(NOE(~3.5))的 饱和转移[35,40].B _较大时,APT 和 NOE 的峰值变宽,导 致低估 APT* 和NOE* [33];同时,随着B _增大,饱和效率达到最大值,饱和溢出和 MT 效应一直增加[25,41],从而导致 APT* 和 NOE* 对比度的下降;(2)APT* 使用三频偏的方法量化时,会受到相邻共振质子的干扰,如胺基质子在2ppm 到3ppm[29].
3.3 洛伦兹拟合
洛伦兹拟合是一个较为简单的最小二乘法Z 谱拟合的量化方法[36].洛伦兹拟合的定义如下[40]:
其中,Ai 、ωi 和δi 分别表示第i个池的幅度、频偏和线宽,N 表示拟合池的个数,一般的洛伦兹拟合是指N =1的情况,如果N 大于1就是多池洛伦兹拟合.
洛伦兹差(Lorentziandifference,LD)是洛伦兹拟合的曲线与Z 谱数据的差[35,42].在该方法中,用洛伦兹函数去拟合Z谱中的-10ppm 到-6.25 ppm、-0.5ppm 到 0.5ppm 和 6.25ppm 到 10 ppm,接着用样条插值完成整个拟合过程,并用拟合的频谱结果作为代表 DS和 MT 效应参考信号;最后,用拟合的结果和 Z 谱数据的差作为 CEST信号进行量化[35,43].
3.4 倒Z谱分析法
最近的一项研究表明 CEST 信号、DS 信号和MT 信号并不是线性叠加在一起,而是反向加在一起的[38].AREX 和 MTR 是两种常见的倒Z谱分析法.其中,AREX 是一种从稳态下得到的标记信号的倒数减去参考信号的倒数的一种量化方 法[29].AREX
可以校正CEST竞争效应中半固体MT 效应、T1、DS效应和水的弛豫效应.AREX 和MTRRex的定义分别为[44]:
图3是不同 B1_sat 时Z 谱和 AREX 的仿真结果,仿真Z 谱的参数设置和图2一样.其中,上面光滑的为Z 谱,左下方为 AREX 的图.由图2 和图3的结果可以看出:相 同功率下,AREX 的峰值比Diff_Z 谱的峰值大(尤其是3池模型),说明 AREX
量化方法中 CEST 的特异性更强.同时也说明:相比饱和交换之前水的Z 谱和饱和交换之后水的Z 谱的差值,AREX 可以降低 MT 效应对 CEST 效应的影响.其中,在0ppm 附近,由于分母为0,所以 AREX 非常的大,这和文献[45]的研究结果一致.
3.5 CEST 比率的方法(CESTratio,CESTR)和参考值归一化后的 CESTR 的方法(CESTRnr)
MRI信号的强度取决于多种参数,包括质子浓度、交换质子的数目、质子交换率、T1、T2、饱和时间和饱和效率.最常用的 CEST 图像的量化是磁化转移率 (magnetizationtransferratio,MTR), MTR 的定义如下:
CEST 比率的定义如下[45]:
在最初对 APT 量化的过程中,选 取 -3.5 ppm 作为参考信号,就是前面3.1 节的 MTRasym 的量化方法[46].最近,Heo等人通过使用插值半固体参考 信 号 (extrapolated semi-solid MT refer- ence,EMR)量化 CESTR[6,47].
参考 值 归 一 化 后 的 CEST 比 率 的 定 义 如下[45]:
CESTR 和 CESTRnr 中的参考信号 (也 就是ZEMR)和 标记信号分别是拟合 2 池 MT 模型和CEST 度量计算中的经过 B0 校正过的Z 谱数据. 仿真证明,在临床的3T 和4.7 T 的时候,CESTR和 CESTRnr 量化更可靠[45].CESTR 的 EMR 谱(ZEMR)可以 通过相应的 ω1和频谱范围 Δ 获得.
CESTR 的Zlab是通过5池Bloch方程获得.具体获得方法只需对前面介绍的2 池模型的 Bloch 方程进行扩展即可.图4和图5分别为B0 =3T和9.4T 时5池模型在B1_sat 分别等于0.5μT,1μT,1.5 μT,2μT,2.5μT 和3μT 六种饱和功率下的基于 Bloch方程的Z 谱和几种常见量化方法:CESTR、 CESTRnr、MTRRex和 AREX 的仿真结果.仿真Z 谱中使用的参数和文献[45]的一样.从图5到图6可以看出:相同 B1_sat 时,3 T 时 Z 谱在3.5ppm 和 2 ppm 的下降均没有9.4 T 时的明显;特别地,在高场(B0 =9.4 T)、B1_sat 较低的时候 (尤 其是在 B1_sat <=2μT),APT 在3.5ppm 和 Amine在2 ppm 的下降看得很清楚;对于 CESTR、CESTRnr、 MTRRex和 AREX 四种量化方法,随着 B0 场强的增大,2ppm 和3.5ppm 的峰值都变得更加明显,而且3.5ppm 的峰值比较窄,2ppm 的峰值比较宽,因此,低场时的 CEST 信息被隐藏.这也是德国神经退行性疾病中心(GermanCenterforNeuro- degenerativeDiseases,DZNE)Zaiss 等人在 2019年2月最新发表深度低场预测高场的信息这篇论文的目的[48].
4 小结与展望
各种量化方法会有具体的适合的实验参数和对应的应用范围,因此可以根据实际需要选取合适的 CEST-MRI的量化方法.同时最近也出现了用比较流行的机器学习中的深度学习进行 CEST-MRI量化[48].在低场的时候采集的数据 SNR 低、硬件容易实现,而且现在的采集工作多数在低场环境下进行.但是,低场强下存在采集的数据 SNR 低、关键信息被干扰隐藏等缺陷.高场中可以得到隐藏在低场中的信息(从图5 和图6 可以看出).但是,高场中采集数据的时候,受到 B0 场不均匀性的影响较大,SNR 也较大,最重要的是由于产生高磁场的设备较为昂贵,因而难以普及,受硬件条件限制较大.Zaiss等人提出同时用低场3 T 和高场 9.4T 采集同一个人的脑部(包括正常人和脑部肿瘤的患者),通过输入3T 时Z谱和输出9.4T的洛伦兹拟合参数,训练深度神经网络的参数.然后利 用训练好的深度神经网络参数,只需要采集低场3 T 的Z 谱,即可预测高场9.4 T 的洛伦兹参数,从而实现预测高场9.4T 的不同 CEST 效应[48].深度 CEST-MRI可以将少数超高场的优点和洞察力带到更广泛的临床应用中.未来深度神经网络的方