三指数模型在胶质瘤术前分级中的价值:3D兴趣区直方图分析

2019-03-15 09:17李建红许乙凯王春红
中国实用神经疾病杂志 2019年1期
关键词:扩散系数直方图胶质瘤

李建红 杜 勇 肖 翔 张 静 许乙凯 王春红△

1)信阳市中心医院医学影像科,河南 信阳 464000 2)南方医科大学附属南方医院影像中心,广东 广州 510515

脑胶质瘤是是颅内原发肿瘤中发病率最高的肿瘤[1]。世卫组织将脑神经胶质瘤分为低级别(WHOⅠ、Ⅱ)和高级别(WHOⅢ、Ⅳ)胶质瘤[2-4]。对胶质瘤进行术前分级在临床实践中至关重要,对治疗方案的选择和预后判断有重大意义。

常规序列的MRI检查在神经胶质瘤术前分级中有一定价值,主要依据病变周围的水肿及强化程度等特征判断,但有时仅依靠病变强化判定胶质瘤分级可能存在困惑[5]。目前,一些新的技术,如波谱、扩散加权成像和灌注成像成为临床评价脑肿瘤的成像方法[6-9]。多b值扩散加权成像不依赖对比剂,表现出明显的优势,并可以同时反映病变内的扩散和灌注相关信息[10]。多b值序列可以用单指数模型、拉伸指数模型和双指数模型进行拟合。部分研究探讨了这些拟合模型的价值[11-15]。最近,三指数模型被提出[16-17],通过拟合可以计算出慢扩散系数(Ds)、快速扩散系数(Df)、灌注相关扩散系数(Dp)及其对应的分数Fs、Ff和Fp。有研究探索了该模型在脑内应用的可行性,认为新的拟合模型很有价值,其可以使用非侵入方式收集大脑灌注和扩散信息[18]。本研究尝试使用3D兴趣区直方图分析法探讨三指数模型相关参数在鉴别低级别与高级别胶质瘤中的价值,目前直方图分析方法还未在三指数模型中应用,有部分文献研究了其他参数的直方图应用[7,19-20]。

1 资料与方法

1.1一般资料收集2016-08—2018-02信阳市中心医院有病理结果的胶质瘤患者46例,男28例,女18例,其中低级别组17例(46.9%),高级别组29例(53.1%)。

入组标准:(1)病理结果证实为胶质瘤;(2)初次发病,行MRI 检查前未经过手术、放疗、化疗及抗血管生成等治疗。排除标准:(1)病理结果不是胶质瘤或未能得到病理结果;(2)资料不全或图像质量不能满足诊断和研究要求;(3)手术与MRI检查间距时间>1个月。

1.2磁共振检查与图像分析

1.2.1 磁共振检查:所有患者磁共振检查全部在信阳市中心医院医学影像科3.0T磁共振扫描系统上完成(Achieva;Philips Healthcare,Netherlands),使用8通道正交头部线圈。扫描方案包括轴位T1WI、轴位T2WI、轴位FLAIR、矢状位T2WI和轴位多b值DWI序列(14b值),所有轴位图像采用统一FOV(24 cm)、层厚(6 mm)和层数(20层);DWI序列(14b值):TR=6 600 ms,TE=Shortest,flip angle=90°,matrix=128×128,NEX=1,Directions=3,b Value=0、10、20、40、60、100、140、160、200、300、500、1 000、1 500、2 000 s/mm2。

1.2.2 多b值扩散加权序列数据处理:多b值扩散加权序列数据处理在第三方软件包上进行(基于Matlab自主开发软件,郭义昊,陈茂东,南方医科大学生物医学工程学院)。

三指数模型拟合参考文献[18]分两步进行。首先,根据公式(1)拟合,其中S代表给定b值的信号强度,Sint代表b值为0时的的信号强度,Ds为慢扩散相关系数,此次拟合使用b值>1 000 s/mm2数据。第二步拟合,其中Df值根据文献在3.0T场强下固定为3.0×10-3mm2/s,然后根据公式(2)拟合,可以得出其他参数,Dp代表灌注相关系数,Fs、Ff和Fp分别代表各扩散系数的相应分数。

S=Sint·exp(-bDs)

(1)

S=S0·[Fp·exp(-bDp)+Ff·exp(-bDf)+Fs·exp(-bDs)]

(2)

1.2.3 ROI的定义:选择包括肿瘤的所有层面计算出各参数图(兴趣区除去肿瘤坏死、出血、钙化及可见的大血管等区域),最后将所勾画ROI拷贝至各参数图上,计算每种参数图上ROI区域直方图相关参数,包括第5、10、90、95百分位数和均数、中位数。

1.3统计学分析使用SPSS 20.0软件进行数据分析,低级别组与高级别组比较,采用独立样本t检验(变量符合正态分布)或Mann-WhitneyU检验(变量不符合正态分布),性别比较采用χ2检验。P<0.05为差异有统计学意义。

2 结果

2.1临床资料低级别组年龄(38.94±14.64)岁,中位年龄40岁(14~63岁);高级别组年龄(51.45±13.88)岁,中位年龄53岁(10~72岁),高级别组年龄略大(t=-2.891,P=0.006)。性别分布比较,低级别组男9例,女8例;高级别组男19例,女10例,2组间性别分布相近,差异无统计学意义(χ2=0.712,P=0.399)。患者临床症状无相对特异性,主要以头痛、癫痫、记忆力低下、反应迟钝和肢体感觉、运动障碍为主。

2.2病理结果低级别组17例,均为WHOⅡ级,高级别组29例,其中WHOⅢ级11例,Ⅳ级18例。组织学类型结果如下:低级别组包括星形细胞瘤12例(图1D),少突胶质星形细胞瘤3例,纤维型星形细胞瘤2例;高级别组包括多形性胶质母细胞瘤17例,间变星形细胞瘤4例,间变少突胶质星形细胞瘤4例,间变少突胶质细胞瘤3例(图2D),小细胞胶质母细胞瘤1例。

2.3 2组三指数模型拟合相关参数比较传统序列有时不能很好地鉴别低级别和高级别胶质瘤(图1A、2A)。低级别组与高级别组比较,Ds-P5、Ds-P10、Ff-mean和Ff-median差异均有统计学意义(P<0.05)。Ds-P5和Ds-P10分别为(0.76±0.08)×10-3mm2/s vs(0.73,0.06)×10-3mm2/s(IQR,P=0.010)、(0.84±0.10)×10-3mm2/s vs(0.80,0.09)×10-3mm2/s (IQR,P=0.013),低级别组高于高级别组(图1B、2B、表1);Ff-mean和Ff-median分别为0.20±0.04 vs 0.22±0.04(P=0.033)、0.15±0.06 vs 0.20±0.04(P=0.003),低级别组低于高级别组(图1C、2C、表1)。2组其他相关参数差异无统计学意义(P>0.05)。

图1 患者男,26岁,左侧额叶肿块影,FLAIR序病灶瘤周水肿较轻,显示3D ROI勾勒方法(A,单层示意图);Ds-p5值为0.71×10-3mm2/s,Ds-p10值为0.81×10-3mm2/s(B);Ff-mean值为0.19,Ff-median值为0.16(C);病理结果为星形细胞瘤(D,WHOⅡ级)Figure 1 Male,26 years old,with amass in the left frontal lobe,with a relatively light edema on FLAIR image,which showed the ROI delineation method (A,single slice);Ds-p5 was 0.71×10-3mm2/s,and Ds-p10 was 0.81 × 10-3mm2/s (B);The Ff-mean value was 0.19,and the Ff-median value was 0.16(C);The pathological results were astrocytoma (D,WHOⅡ)

参数低级别组高级别组P值DsP50.76±0.080.73,0.060.010&DsP100.84±0.100.80,0.090.013&Ffmean0.20±0.040.22±0.040.033#Ffmedian0.15±0.060.20±0.040.003#

注:#t检验;&Mann-WhitneyU检验

2.4诊断指标效能评价区分低级别组和高级别组胶质瘤,当Ds-P5的截断值取0.77×10-3mm2/s时,其敏感性和特异性分别为52.9%和93.1%;当Ds-P10的截断值取0.86×10-3mm2/s时,其敏感性和特异性分别为47.1%和96.6%;Ff-mean的截断值取0.23时,其敏感性和特异性分别为41.4%和89.2%;Ff-median的截断值取0.17时,其敏感性和特异性分别为75.9%和64.7%(表2)。

Ds-P5的ROC曲线下面积(AUC)为0.728(P=0.010),Ds-P10的AUC为0.722(P=0.013),Ff-mean的AUC为0.665(P=0.064),Ff-median的AUC为0.721(P=0.013)。这些参数能为鉴别低级别与高级别胶质瘤提供有价值信息(表2)。

图2 患者男,55岁,右侧顶叶肿块影,FLAIR序病灶瘤周水肿不明显,显示3D ROI勾勒方法(A,单层示意图);Ds-p5值为0.75×10-3mm2/s,Ds-p10值为0.84×10-3mm2/s(B);Ff-mean值为0.19,Ff-median值为0.14(C);病理结果为间变少突胶质细胞瘤(D,WHOⅢ级)Figure 2 Male,55 years old,with amass in the right parietal lobe,with a relatively light edema on FLAIR image,which showed the ROI delineation method (A,single slice);Ds-p5 was 0.75×10-3mm2/s,and Ds-p10 was 0.84×10-3mm2/s (B).The Ff-mean value was 0.19,and the Ff-median value was 0.14(C).The pathological results were anaplastic oligodendrocytoma (D,WHO Ⅲ)

参数AUCP值截断值敏感性/%特异性/%DsP50.7280.0100.7752.993.1DsP100.7220.0130.8647.196.6Ffmean0.6650.0640.2341.489.2Ffmedian0.7210.0130.1775.964.7

3 讨论

胶质瘤术前分级对治疗方法的选择和预后评价很重要。磁共振传统序列在术前分级中有时会有局限,因不同级别胶质瘤的影像表现会有重叠。部分研究已表明,血管渗透性的信息有助于胶质瘤分级[7]。灌注成像技术由于能够评价肿瘤内毛细血管通透性,成为一种最有效的胶质瘤分级方法[11]。

随着多b值扩散加权成像不同拟合模型的提出[9-10],可同时评价肿瘤组织中的扩散和灌注相关信息,为胶质瘤分级提供了新的模式,而三指数模型被认为是适合描述组织的扩散和灌注信息的拟合模型。

细胞密度和有丝分裂活动是神经胶质瘤组织病理学重要的一个分级标准,有文献认为DWI序列相关参数ADC值和D值与细胞密度相关性较高,对分级有较高价值[14,21]。本研究显示,Ds-P5和Ds-P10在低级别和高级别组间有明显差异,能显示不同级别肿瘤内细胞密度及水分子扩散的差异,对胶质瘤分级有潜在的价值。

肿瘤新生血管形成是恶性神经胶质瘤的特点,因此,可通过定量评价肿瘤内的微血管为分级提供依据[22-23]。本研究显示,Ff-mean和Ff-median在低级别和高级别组间有明显差异,与以往有关灌注相关参数的研究结果一致。

由于参数拟合的鲁棒性较差,本次研究未评估Dp值,与以往研究结果有出入[16-17]。部分研究中给出的Dp值结果的标准差很大[21],可能由于Dp值与双指数模型中的D*值一样信噪比较低所致[11]。

本研究仍存在局限性。首先,三指数模型相关参数值会受到b值分布的影响[11,16,21,24],最优b值方案,仍需进一步探讨。更高的b值能更好地消除灌注相关的扩散,然后得到一个更接近实际的扩散系数值和灌注相关扩散系数值,本研究使用的b值最高为2 000 s/mm2,这比其他作者使用的b值要小[16,18],可能影响相关参数的准确性。其次,本研究中ROI的定义方法使用的是包括全瘤的异常信号区域,虽可见的坏死区、出血、大血管可有排除,但部分较小的坏死、出血等区域可能会影响到结果的准确性。最后,本研究未探讨三指数模型直方图分析相关参数与肿瘤的分子标志物的关联,有文献证实分子标志物与肿瘤的诊断和预后判断有明显关系[25-29],希望在下一步的研究中探讨相关问题。

三指数模型作为一种新的拟合模型,可在一个多b值的扩散加权序中分离出扩散相关信息和灌注相关信息,通过3D直方图分析能为胶质瘤的术前分级提供有价值信息。

猜你喜欢
扩散系数直方图胶质瘤
成人高级别脑胶质瘤术后复发相关因素分析
符合差分隐私的流数据统计直方图发布
青藏工程走廊冻融土热扩散系数特性与预测模型研究
BFAR 在胶质瘤中的表达及其与胶质瘤预后的关系
恐惧应激对胶质瘤影响机制及干预研究进展
基于FPGA的直方图均衡图像增强算法设计及实现
基于相关性分析方法建立混凝土RCM扩散系数多因素模型
磁共振体素内不相干运动在肺良恶性肿瘤的诊断效能及肺癌化疗疗效评估的应用
用直方图控制画面影调
中考频数分布直方图题型展示