脊髓由大量神经元、胶质细胞及血管等组织构成，其结构在发育过程中会随年龄增长出现不同程度的退化^[1,2]。扩散张量成像（diffusion tensor imaging, DTI）虽然可以从微观上评估组织结构特性^[3]，但对神经纤维较多且结构复杂的区域反映不够准确（如灰质)^[4]。所以需要一种既能显示纤维束空间方向性又能反映其结构复杂性的影像学手段，来如实反映脊髓微观结构改变。

       近年来出现的扩散峰度成像（diffusion kurtosis imaging, DKI)技术是以研究水分子受限状态下扩散的特性为基础的成像方法^[5]，比DTI更真实地反映组织结构中水分子扩散受限的不均质性^[6,7]。从而能更准确地评估脊髓灰质神经元、白质纤维素等微观结构及其损伤和再生修复的动态过程^[8]。DKI在脊髓的应用研究国内外报道较少，本研究旨在探讨正常人颈髓DKI序列参数值变化规律及其与年龄相关性。

1　材料与方法

1.1　研究对象

       选取2015年8月至2016年3月中国医科大学附属盛京医院36例成年志愿者(≥18岁），对其进行常规颈髓磁共振检查和DKI序列检查。其中男20例，女16例，年龄18~ 65岁，平均46岁。按年龄分为3组：30岁及以下为青年组（12例），31~ 50岁为中年组（12例），51岁及以上为老年组（12例）。入选标准：以往无中枢神经系统疾患和其他影响中枢神经系统的全身性疾患病史。常规MRI扫描图像质量满意，且颈髓形态和信号无异常者。检查前均签署知情同意书。

1.2　仪器设备及参数

       使用GE 3.0 T超导型磁共振成像仪（GE Medical System, Healthcare Discovery MR750 MRI, USA)，选用头颈联合线圈进行磁共振扫描。采取快速扰相梯度回波序列(fast spoiled gradient echo, FSPGR)进行定位扫描，最上层面达中脑水平，最下达胸1椎体水平。使用快速自旋回波技术序列（fast spin echo, FSE）行矢状面及横轴面T2WI扫描，反转恢复序列（fluid attenuated inversion recovery, FLAIR)行矢状面T1WI扫描，平面回波序列（echo-planar imaging, EPI）进行DKI扫描。

       常规扫描：矢状位T1WI：TR=2100 ms，TE=2.8 ms；T2WI：TR=2500 ms，TE=124.8 ms；层厚6 mm，层间距0，FOV 240 mm，矩阵为256 × 256。横轴位T2WI：TR=3861 ms，TE=127.7 ms；层厚3 mm，FOV 200 mm，矩阵为512×512. DKI扫描：扩散加权系数（b值）分别为0、1000、2000 s/mm²，扩散敏感梯度取25个不同方向。TR=3600 ms，TE=108 ms，层厚3 mm，层间距0，FOV 200 mm，矩阵为512×512，共6个间盘层；采集时间约7 min。

1.3　数据处理及图像分析

       数据全部存储在GE工作站内，ADW 4.6 (Advantage Workstation；Sun Microsystems, Santa Clara, Calif)，使用软件（Functool 9.4.05a software, GE Healthcare）进行DKI数据处理。依据常规正中矢状位T2WI来定位颈髓轴位C1~ 7间盘层面，每一个层面分别获得弥散各向异性分数（fractional anisotropy, FA)、平均弥散系数（mean diffusivity, MD)、平均弥散峰度（mean kurtosis, MK）参数图，并测量各参数值，同时生成伪彩图。图像评价由两位高年资的影像专家分别完成，并达成一致意见。采用设置感兴趣区（region of interest, ROI）法进行DKI参数值测定，分别将4个相同大小的ROI（面积3.4 mm²，约3个像素大小），放置于脊髓的双侧前角灰质(gray matter, GW)和侧索白质（white matter, WM)，每个ROI重复测量3次，取平均值。

1.4　统计分析

       采用SPSS 11.0软件包，计量资料以均数±标准差(±s）表示，选用单因素方差分析、Pearson相关分析。P<0.05为差异有统计学意义。

2　结果

       所选志愿者均获得满意的常规颈髓MRI图像及DKI图像。颈髓各间盘层面FA、MD及MK参数图及伪彩图见图1。伪彩图显示：FA为均匀的红色信号，MD为较均匀的淡绿色信号，MK为较均匀的淡红色信号。颈髓各层面DKI不同参数测量值见表1，不同层面参数值变化趋势见图2。从C1~ 2至C6~ 7间盘水平，GM和WM的DKI总体趋势为：MD逐渐增高，FA和MK逐渐减低。C2~ 3层面FA-GM和MD-WM与其他层面存在显著性差异(F=10.01，P<0.05；F=16.94，P<0.05)，C6~ 7层面MK-GM和MD-GM亦与其他层面存在显著性差异(F=10.11，P<0.05；F=10.01，P<0.05)，其他各参数在各节段层面无显著性差异(P<0.05)。

       颈髓FA、MD及MK各年龄组参数值与年龄相关性见图3、表2，FA-WM年龄呈负相关(r=-0.507，P<0.05)，MK-GM和MK-WM年龄呈正相关(r=0.436，P<0.01；r=0.249，P<0.05)。

       FA-WM和MK-GM在青年组、中年组、老年组间均有统计学意义（表2，P<0.05)，MK-WM在青年组与老年组、青年组与中年组间有统计学意义（表2，P<0.05)，而在中年组与老年组间无统计学意义。FA-GM相关性无统计学意义(P<0.05)，但均有随年龄增长而减少的趋势；MD-GW及MD-WM值与年龄的相关性无统计学意义(P<0.05)，但均有随年龄增长而增加的趋势。

点击查看大图

图1  男，35岁，颈2-3间盘水平脊髓。A为轴位T2WI图，B~ D为FA、MD、MK参数图，E~ G为图B~ D的伪彩图
Fig. 1  Male, 35 years old, Spinal on C2- 3. A is T2WI axis image, B—D is the FA, MD, MK parameters image after scanning, E—G is pseudo color map of figure B—D.

点击查看大图

图2  颈段脊髓DKI不同参数值在各间盘层面变化趋势图
Fig. 2  Cervical spinal cord DKI different parameter values in every dish of level changes between the trends of average figure

点击查看大图

图3  A~ F分别为DKI不同参数值随年龄变化的散点图。FA-GM和FA-WM随年龄有下降趋势，且有相关性；MK-GM和MK-WM随年龄有增加趋势，且有相关性
Fig. 3  A—F are scatter diagram for different DKI parameter values change along with the age, respectively. FA-GM and FA-WM decline with age, MK-GM and MK-WM tendency to increase with age.

点击查看表格

表1  正常人C1-C6椎间隙水平DKI参数测量值(±s)
Tab. 1  The volunteers in different segments of DKI parameter mean value(±s)

点击查看表格

表2  健康志愿者各年龄组间颈部脊髓不同部位DKI参数值比较(±s)
Tab. 2  Each age group of the volunteers in different of the cervical spinal DKI parameter values(±s)

3　讨论

3.1　DKI脊髓应用基本原理及技术

       Jensen^[8]于2005年首次提出DKI技术理论，该技术在非高斯模式（non-Gaussian diffusion model)测定水分子扩散运动的程度，比扩散加权成像(diffusion weighted imaging, DWI）和DTI更能准确反映水分子扩散的受限程度^[6,7,8]，同时还能得到DTI相关参数，如表观扩散系数（apparent diffusion coefficient, ADC)、FA等^[8,9]。而脊髓DKI技术能更好地反映更丰富的组织内水分子非正态分布特征，尤其是在灰质比DTI成像能显示更多的微观改变^[7,8]。有研究显示^[9]，使用3.0 T MR进行脊髓DKI扫描时，至少需要3个不同的b值和15个梯度方向（motion probing gradient, MPG)(30个MPG最佳）。本研究DKI序列采用了3个b值25个方向，得到了满意的颈髓DKI图像，与Hori等人^[6]研究结构相一致。本研究扫描时间约7 min，低于Wheeler- Kingshott等^[9]30 min的扫描时间，更易在临床得到推广。

3.2　DKI在正常人颈髓的应用

       脊髓白质主要结构为纵行的神经纤维，其粗细差异很大，多数是有髓神经纤维，主要起传导作用。灰质主要由神经元的胞体构成，前角内大多是躯体运动神经元，胞体大小不等，其轴突较粗细不均，且走行方向比较复杂^{[1, 10]}。Raz等^[11]认为DKI比DTI能够更全面地描述脊髓内正常结构及病变特征。因为DKI能敏感地反映脊髓轴突退变和脱髓鞘以及再髓鞘化过程中水分子扩散各向异性的改变，尤其是对灰质不均质性方面有明显优势^{[5, 12]}。

       本研究结果显示C2~ 3层面FA-GM和MD-WM与其他层面存在显著性差异(F=10.01 ，P<0.05； F=16.94，P<0.05)，C6-7层面MK-GM和MD-GM亦与其他层面存在显著性差异(F=10.11，P <0.05；F=10.01，P<0.05)，FA最大值在C2~ 3，MK最大值在C6~ 7，其他灰白质各参数在不同层面未见显著差异性。Wheeler-Kingshott等^[9]提出FA值随着脊髓位置而变化与本研究结果相一致，在C2~ 3水平达到最大值（图1，表1)且越向尾侧FA值越低（图1A、图1B)。由于灰质密度在不同部位下降速率也不尽相同，有研究显示白质和灰质的MK峰值会逐步减少^[13]。本研究结果MK最大值在C1~ 2(MK-GM)、C2~ 3 (MK-WM)，且不同层面变化趋势性不明显，可能是受测量误差及年龄因素的影响^[11]。

       对于脊髓微观结构的改变与年龄的相关性，有研究^[11]认为FA值与年龄负相关，与本研究结果相一致。可能与脊髓轴突、髓鞘变性退化时含水量减少，纤维方向性发生改变和结构复杂程度也随之减低有关。另有文献表明^[7]，脊髓不同间盘水平的扩散值不完全一致，并且随着年龄的增加，但并未把脊髓白质和灰质分开评估。而本研究将脊髓灰质和白质分别进行测量，结果显示FA-GM和FA-WM年龄呈负相关(r=-0.507，P<0.05)，MK-GM和MK-WM年龄呈正相关(r=0.436， P<0.01；r=0.249，P<0.05)；虽然FA-GM也有随年龄增长而减少的趋势，但没有统计学意义。除原来可能的测量误差外，还可能是由于脊髓灰质本身结构相对复杂，纤维交叉成分较多造成。

       MK是目前DKI技术最常用的参数，代表组织沿空间各方向弥散峰度的平均值，值越大表明弥散受限越严重^[8]。Falangola等^[10]认为不同的年龄阶段表现的MK值也不同，MK值与年龄之间有着明显的关联。可能是神经组织发育过程中，不同层面的白质和灰质比例与数量变化有关。本研究结果显示MK-GM和MK-WM年龄呈正相关(r=0.436，P<0.01；r=0.249，P<0.05)，可能是由于随年龄增长，脊髓内水分子弥散受限程度会逐渐增加。MK-WM在青年组与老年组、青年组与中年组间有统计学意义(P<0.05)，而在中年组与老年组间无统计学意义。可能由于中年组变化幅度小于老年组，或样本量较少及测量误差所致，下一步将深入研究。

       MK-GM年龄呈正相关强于MK-WM（图3)，可能是由于MK对灰质复杂结构的敏感性高于纵行纤维的白质，但尚需加大样本量进一步研究。所以在对于脊髓白质和灰质的分析，需充分考虑年龄因素和脊髓不同节段DKI参数差异。虽然DKI在脊髓灰质方面的研究还比较少，但至少说明对于脊髓灰质损伤后的改变，其参数值会发生变化。

4　局限性

       DKI技术虽然是在非高斯基础上研究水分子的扩散情况^[3,4,5,6]，但还并没达到理想状态下测量水分子所有方向扩散的程度，描述水分子的扩散还不够全面^[8]。

       本研究所选取的样本量相对较少，均为成年人，且年龄跨度相对较窄，需增加样本量进一步研究。ROI选取方法也可能存在误差，拟进一步优化。

       年龄及性别因素：随着社会发展，现在人们生活习惯发生改变，颈椎病的发生趋向年轻化，DKI参数与年龄相关趋势变化有待于进一步研究。有研究显示^[14]，性别和扩散指标之间显著相关性，但有文献报道，DKI参数值与性别未见明显相关性。下一步研究将充分考虑性别因素^[7]，并期待衡量DKI的客观标准出现^{[7, 11]}。

       综上所述，DKI技术能反映出健康成人脊髓灰、白质微观结构的变化规律^[11]，其参数值变化与年龄存在相关性，MK可能是反映神经系统复杂微观结构和损伤改变的高度敏感的指标，将为判断脊髓损伤程度及评估预后方面提供很大帮助。