回顾性分析2021年10月至2022年9月于苏北人民医院就诊的宫颈病变病例。纳入标准：（1）年龄18～80岁；（2）原发性宫颈鳞状细胞癌；（3）患者术前行T2WI序列及DR-CSI序列盆腔扫描；（4）由妇产科临床医师根据国际妇产科联盟（International Federation of Gynecology and Obstetrics, FIGO）2018临床分期标准进行准确分期^[1]；（5）出现性交后出血、阴道不规则排血、排液、腰骶部疼痛等临床症状；（6）病灶均为首发初诊。排除标准：（1）图像存在严重伪影；（2）临床资料不完整。最终纳入病例60例，年龄24～78岁，平均56岁。根据中国抗癌协会妇科肿瘤专业委员会发布的《宫颈癌诊断与治疗指南（2021年版）》：ⅠA期、ⅠB1期、ⅠB2期、ⅡA1期患者推荐手术治疗；ⅠB3期、ⅡA2及ⅡB期及以上患者首选同步放化疗^{[14, 15]}。研究对象分组：早期组（ⅠB1、ⅠB2期、ⅡA1期）26例，晚期组（ⅠB3期、ⅡA2期、ⅡB期、Ⅲ期、Ⅳ期）34例。本研究遵循《赫尔辛基宣言》，研究方案经苏北人民医院伦理委员会批准，免除受试者知情同意，批准文号：2022ky206。

患者术前采用3.0 T磁共振仪（MR 790，中国上海联影医疗）和12通道腹部体线圈进行盆腔MRI检查，所有患者行MRI检查前均禁食4～6 h。检查前0.5 h肌内注射盐酸消旋山莨菪碱5 mg（国药集团容生制药有限公司），以减少肠道蠕动，适当排空膀胱以减少尿液电解质伪影。患者取仰卧位，脚先进，放松平静呼吸，扫描范围为双侧髂骨上缘至股骨颈水平。

脂肪抑制T2WI序列行矢状位、轴位扫描。矢状位扫描参数：TR 6200 ms，TE 94 ms，层数29，层厚4 mm，层间距10 mm，视野350 mm×280 mm，翻转角90°，矩阵352×352；轴位扫描参数：TR 11 817 ms，TE 94 ms，层数46，层厚4 mm，层间距10 mm，视野350 mm×280 mm，翻转角90°，矩阵352×352。DR-CSI行单次激发自旋回波平面成像序列，获得5个回波时间（47.6、60.0、70.0、80.0、90.0 ms）和6个b值（0、200、400、800、1200、1600 s/mm²）组合的轴位图像（在子宫T2WI矢状位定位并与宫颈管垂直），扫描参数：TR 2000 ms，层数5，层厚4 mm，层间距10 mm，FOV 260 mm×200 mm，翻转角90°，矩阵128×128，带宽2070 Hz，扫描时间10 min。

将DR-CSI扫描的5个TE和6个b值组合的30组图像导入Matlab软件（2021a, MathWorks, Natick, MA, USA），由1名具有5年盆腔MRI诊断经验的放射科主治医师，以T2轴位图像为参照，在TE=60 ms和b=800 s/mm²组合序列选取肿瘤横断面经线最大的层面进行ROI勾画，避开出血、坏死、囊变、水肿区域。再由1名具有10年盆腔MRI诊断经验的放射科主任医师对勾画结果进行检查并确认最终ROI。

DR-CSI将每个体素的MRI信号建模为以T2值和D值（ADC值）为特征的连续指数衰减函数求和：Sx,y,TE,b=∬wx,y,T2,Dexp-TET2exp-bDdT2dD=Lw，其中Sx,y,TE,b表示每个TE-b编码的体素测量信号，wx,y,T2,D表示要重建的体素T2-D谱（T2和ADC值构建扩散-弛豫相关谱），L表示拉普拉斯变换^[12]。通过应用非负性约束和空间总变化约束：w=argmin||S-Lw|||F2+TVxyw,subjecttow≥0，使用40个D值（0.001～5.000）×10^-3 mm²/s和40个T2值（0～300）ms的任意组合来构建每个体素的T2-D谱。体素T2-D谱根据T2值和D值范围的不同划分为六个区室，即为六室信号模型（图1），六室信号模型的分界值：T2值为50 ms，D值分别为1×10^-3 mm²/s、3×10^-3 mm²/s。微观区室的信号占比分数（例如f_A对应区室A），通过对区室平均的T2-D谱上的每个区室信号进行积分，然后进行归一化来计算测量：fAx,y=∫AreaundercompartmentAwx,y,T2,DdT2dD∫entireT2-Dspacewx,y,T2,DdT2dD。计算ROI的定量参数：ADC、T2值、各区室信号占比分数（f_A～f_F）。

点击查看大图

图1

扩散-弛豫相关谱成像（DR-CSI）六区室信号模型图。D（表观扩散系数）值范围为（0.001～5.000）×10^-3 mm²/s，T2值范围为0～300 ms。T2（50 ms）、D（1×10^-3、3×10^-3 mm²/s）分界值将扩散-弛豫相关谱（T2-D谱）划分为六个区室（A～F）。

Fig. 1

Diffusion-relaxation correlation spectrum imaging (DR-CSI) six-compartment signal model. D (apparent diffusion coefficient) values range from (0.001-5.000) ×10^-3 mm²/s, and T2 values range from 0-300 ms. The divided six areas (ie, six compartments A-F) on diffusion-relaxation correlation spectrum (T2-D spectrum) space according to the cut-off values of T2 value 50 ms, D value 1×10^-3 mm²/s and 3×10^-3 mm²/s, respectively.

点击查看大图

图1

扩散-弛豫相关谱成像（DR-CSI）六区室信号模型图。D（表观扩散系数）值范围为（0.001～5.000）×10^-3 mm²/s，T2值范围为0～300 ms。T2（50 ms）、D（1×10^-3、3×10^-3 mm²/s）分界值将扩散-弛豫相关谱（T2-D谱）划分为六个区室（A～F）。

Fig. 1

Diffusion-relaxation correlation spectrum imaging (DR-CSI) six-compartment signal model. D (apparent diffusion coefficient) values range from (0.001-5.000) ×10^-3 mm²/s, and T2 values range from 0-300 ms. The divided six areas (ie, six compartments A-F) on diffusion-relaxation correlation spectrum (T2-D spectrum) space according to the cut-off values of T2 value 50 ms, D value 1×10^-3 mm²/s and 3×10^-3 mm²/s, respectively.

使用SPSS 26.0软件进行统计学分析。采用Kolmogorov-Smirnov检验方法分析ADC、T2和f值是否符合正态分布。符合正态分布的参数以均数±标准差（x¯±s）格式表示，两组间比较采用两独立样本t检验，非正态分布的参数以中位数（上下四分位数）[M（P₂₅，P₇₅）]表示，两组间比较采用Mann-Whitney U检验。采用二元logistic回归分析进行DR-CSI六区室联合和ADC联合T2参数分析。使用MedCalc 20.022软件进行受试者工作特征（receiver operator characteristic, ROC）曲线绘制并计算曲线下面积（area under the curve, AUC），根据约登指数来确定各参数的最佳诊断临界值、敏感度、特异度。采用DeLong检验比较各参数AUC的差异，P＜0.05为差异有统计学意义。

鉴别早晚期宫颈癌DR-CSI参数中，晚期宫颈癌以f_B最高（0.30±0.12），早期宫颈癌以f_C最高（0.31±0.11）。ADC（P=0.022）、f_C（P=0.006）在早期宫颈中明显高于晚期宫颈癌，而f_B（P=0.009）在晚期宫颈癌中明显高于早期宫颈癌（表1，图2）。

点击查看表格

表1

早晚期宫颈鳞癌ADC、T2及各微观区室信号占比分数比较结果

Tab. 1

Comparative results of ADC and T2 with the fractions of each compartment in early and advanced stage squamous cervical carcinoma

表1

早晚期宫颈鳞癌ADC、T2及各微观区室信号占比分数比较结果

Tab. 1

Comparative results of ADC and T2 with the fractions of each compartment in early and advanced stage squamous cervical carcinoma

参数	早期	晚期	t值/Z值	P值
ADC/（×10-3 mm2/s）	1.32±0.22	1.08±0.22	-2.502	0.022
T2/ms	84.97±19.60	100.65±20.56	1.784	0.090
fA	0.17（0.13，0.19）	0.18（0.15，0.21）	-0.845	0.398
fB	0.15±0.12	0.30±0.12	2.888	0.009
fC	0.31±0.11	0.20±0.05	-3.082	0.006
fD	0.11（0.07，0.27）	0.12（0.07，0.17）	-0.282	0.778
fE	0.20±0.07	0.16±0.07	-1.130	0.272
fF	0.02±0.02	0.02±0.02	0.153	0.880

注：ADC为表观扩散系数；f_A～f_F为扩散-弛豫相关谱划分的6个微观区室信号占比分数。

点击查看大图

图2

女，55岁，早期宫颈鳞癌，T2-D谱（2A）和DR-CSI参数ADC（2B）、f_B（2C）、f_C（2D）图，ADC=1.07×10^-3 mm²/s，f_B、f_C分别为0.21和0.25，临床分期为早期（ⅠB2期）。图3　女，26岁，宫颈鳞癌患者，T2-D谱（3A）和DR-CSI参数ADC（3B）、f_B（3C）、f_C（3D）图，ADC=0.99×10^-3 mm²/s，f_B、f_C分别为0.40和0.15，临床分期为晚期（ⅢC期）。T2-D：T2值与ADC值构建扩散-弛豫相关谱；DR-CSI：扩散-弛豫相关谱成像；ADC：表观扩散系数；f_B、f_C：T2-D谱划分的微观区室信号占比分数。

Fig. 2

A 55-year-old female early SCC patient. Corresponding T2-D spectrum (2A) and parameter maps [ADC (2B), f_B (2C), f_C (2D)] are shown, ADC=1.07×10^-3 mm²/s, f_B, and f_C are 0.21 and 0.25 respectively. The clinical stage is early (stage ⅠB2).

Fig. 3

A 26-year-old female SCC patient. Corresponding T2-D spectrum (3A) and parameter maps [ADC (2B), f_B (2C), f_C (2D)] are shown, ADC=0.99×10^-3 mm²/s, f_B and f_C are 0.40 and 0.15 respectively. The clinical stage is advanced (stage ⅢC). SCC: squamous cervical carcinoma; T2-D: T2 and ADC values to construct diffusion-relaxation correlation spectrum (T2-D spectrum); DR-CSI: diffusion-relaxation correlation spectrum imaging; ADC: apparent diffusion coefficient; f_B, f_C: percentage fraction of microscopic compartment signal by the T2 and ADC values to construct diffusion-relaxation correlation spectrum.

点击查看大图

图2

女，55岁，早期宫颈鳞癌，T2-D谱（2A）和DR-CSI参数ADC（2B）、f_B（2C）、f_C（2D）图，ADC=1.07×10^-3 mm²/s，f_B、f_C分别为0.21和0.25，临床分期为早期（ⅠB2期）。图3　女，26岁，宫颈鳞癌患者，T2-D谱（3A）和DR-CSI参数ADC（3B）、f_B（3C）、f_C（3D）图，ADC=0.99×10^-3 mm²/s，f_B、f_C分别为0.40和0.15，临床分期为晚期（ⅢC期）。T2-D：T2值与ADC值构建扩散-弛豫相关谱；DR-CSI：扩散-弛豫相关谱成像；ADC：表观扩散系数；f_B、f_C：T2-D谱划分的微观区室信号占比分数。

Fig. 2

A 55-year-old female early SCC patient. Corresponding T2-D spectrum (2A) and parameter maps [ADC (2B), f_B (2C), f_C (2D)] are shown, ADC=1.07×10^-3 mm²/s, f_B, and f_C are 0.21 and 0.25 respectively. The clinical stage is early (stage ⅠB2).

Fig. 3

A 26-year-old female SCC patient. Corresponding T2-D spectrum (3A) and parameter maps [ADC (2B), f_B (2C), f_C (2D)] are shown, ADC=0.99×10^-3 mm²/s, f_B and f_C are 0.40 and 0.15 respectively. The clinical stage is advanced (stage ⅢC). SCC: squamous cervical carcinoma; T2-D: T2 and ADC values to construct diffusion-relaxation correlation spectrum (T2-D spectrum); DR-CSI: diffusion-relaxation correlation spectrum imaging; ADC: apparent diffusion coefficient; f_B, f_C: percentage fraction of microscopic compartment signal by the T2 and ADC values to construct diffusion-relaxation correlation spectrum.

鉴别早晚期宫颈癌诊断效能中，单参数诊断效能以f_B最高（AUC=0.818），敏感度、特异度分别为63.6%和90.0%，其次是f_C（AUC=0.809）和ADC（AUC=0.773），三者最佳诊断临界值分别为0.276、0.205、1.043×10^-3 mm²/s。多参数联合以DR-CSI六区室联合（f_A～f_F）诊断效能最高，AUC为0.973，明显优于ADC（AUC=0.773，P=0.031）和T2（AUC=0.790，P=0.040）的诊断效能，较ADC联合T2的诊断效能（AUC=0.809，P=0.065）有所提高（表2，图3）。

点击查看表格

表2

ADC、T2及各微观区室信号占比分数鉴别早晚期宫颈癌的诊断效能

Tab. 2

Diagnostic efficacy of ADC and T2 with the fractions of each compartment in identifying early and advanced stage squamous cervical carcinoma

表2

ADC、T2及各微观区室信号占比分数鉴别早晚期宫颈癌的诊断效能

Tab. 2

Diagnostic efficacy of ADC and T2 with the fractions of each compartment in identifying early and advanced stage squamous cervical carcinoma

参数	AUC（95% CI）	敏感度/％	特异度/％	最佳诊断临界值
ADC/（×10-3 mm2/s）	0.773（0.540～0.924）	54.5	90.0	1.043
T2/ms	0.790（0.473～0.884）	72.7	70.0	91.056
fB	0.818（0.591～0.950）	63.6	90.0	0.276
fC	0.809（0.580～0.945）	100.0	60.0	0.205
ADC＋T2	0.809（0.580～0.945）	63.6	90.0	—
DR-CSI	0.973（0.793～1.000）	100.0	99.9	—

注：ADC为表观扩散系数；f_B、f_C为扩散-弛豫相关谱划分的微观区室信号占比分数；DR-CSI为扩散-弛豫相关谱成像；AUC为受试者工作特征曲线下面积。

本研究将DR-CSI扫描的序列进行区室划分，生成每个体素相关的六室信号模型。六室信号模型的分界值T2为50 ms、ADC分别为1×10^-3 mm²/s、3×10^-3 mm²/s。此分区的优势在于可以根据宫颈癌扩散和弛豫的特点划分出利于分期鉴别的区室。ADC=3×10^-3 mm²/s被认为可以区分扩散和灌注，本研究中宫颈癌组织均为扩散受限的肿瘤组织，因此ADC＞3×10^-3 mm²/s的灌注成分不具有鉴别价值，与本试验中f_E、f_F差异无统计学意义的结果相符。对于ADC临界值，韩萍等^[19]利用常规MRI和DWI序列对宫颈癌进行术前分期诊断的研究中，ⅠB、ⅡA、ⅡB、Ⅲ～Ⅳ期平均ADC分别为（1.16±0.25）×10^-3、（0.95±0.16）×10^-3、（0.82±0.12）×10^-3、（0.71±0.11）×10^-3 mm²/s，由此可见宫颈癌病灶ADC在1×10^-3 mm²/s上下浮动，且晚期通常小于1×10^-3 mm²/s，早期通常大于1×10^-3 mm²/s，与本研究中晚期组以f_B最大，而早期组以f_C最大，且f_B、f_C差异具有统计学意义的结果相符。对于T2临界值，LI等^[8]通过T2 mapping评估宫颈癌病理学在特征中发现，肿瘤组织的平均T2值为（88.45±12.01）ms，而本试验T2值定为50 ms，也可有利于宫颈癌组织中早晚期的鉴别。

本研究通过计算ADC、T2、f值，并在早晚期宫颈癌中进行对比，晚期宫颈癌以f_B最高，早期宫颈癌以f_C最高，ADC、f_B、f_C在两组间差异具有统计学意义（P＜0.05）。ADC可以反映水分子在组织中的扩散信息，本研究中晚期组ADC显著低于早期组，与既往研究结果一致^{[7,20, 21]}。T2反映局部物理和化学微环境的变化。随着宫颈癌分期的增加，肿瘤组织的异质性和复杂性增加，这可以从体素中每个组织隔室的f变化来体现。在我们的研究中，晚期宫颈鳞癌以f_B最高，表明晚期宫颈鳞癌以ADC值较低、T2值较高的区室信号为主。与早期宫颈鳞癌相比，晚期宫颈鳞癌通常具有肿瘤增殖速度较快、肿瘤细胞密度增加、细胞外间隙缩小的特点，明显限制水分子的扩散，使ADC值更低^{[22, 23]}。此外，随着肿瘤分期的增加，肿瘤体积逐渐增大，血供更加丰富，肿瘤血管更加杂乱迂曲，微坏死区逐渐增多^{[24, 25]}，以上原因可能导致的T2横向弛豫时间延长，因此，早期宫颈癌中以ADC值较高而T2值较低的区室的f_C为主。SLATOR等^[26]研究体内人类胎盘的T2^*-D光谱，在对照组（健康对照）中观察到三个具有相似T2^*值的不同ADC峰，慢扩散成分与胎盘内的区域相关，而中等和快速扩散在子宫壁中更为突出；而对于试验组（慢性高血压、先兆子痫和胎儿生长受限），存在向较短的T2^*值、较慢的扩散转变的趋势。由此可见DR-CSI具有反映组织内部病理生理变化能力。此部分研究结果显示，ADC、f_B、f_C可作为早晚期宫颈癌的鉴别指标，为宫颈癌的术前分期诊断提供价值。

本研究通过对ROC曲线的分析，在宫颈癌早晚期鉴别中单参数中f_B的AUC最大，对宫颈癌早晚期鉴别的诊断效能优于传统ADC和T2。多参数联合分析中，DR-CSI序列的六区室联合分析的诊断效能优于ADC联合T2。这是由于在分析来自包含多个组织隔室的体素的MRI信号时，ADC来自单指数模型，是所有隔室扩散系数的加权平均值，而不是直接反映各个组织隔室的比例。相对于ADC，f值可以反映各个组织隔室对体素的贡献，因此与微观组织隔室的相关性更强^[12]。而随着宫颈癌分期的变化，其内部微环境发生变化，体现在各个组织隔室贡献发生变化，这可能是DR-CSI诊断效能更高的原因。正如CHATTERJEE等^[13]将DR-CSI序列组织隔室与离体前列腺组织成分（上皮、基质、管腔）相对应测得的组织成分的体积分数，相较于ADC、T2，与Gleason分级的Spearman相关系数显著性更高。DR-CSI另一优势是可以包含扩散和弛豫两方面耦合信息，比单一的ADC或T2提供更多生物学信息。以往的研究通常使用T2WI和DWI的组合独立地测量T2值和ADC值，因此隐含地假设每个体素中的T2和ADC是独立的^{[27, 28]}。然而，已有多项研究证实ADC和T2经常是耦合的，如MAI等^[29]在研究定量T2值对前列腺癌的检测和分级中发现，T2值和ADC值有较强的相关性（r=0.772），短T2耦合低ADC值提示前列腺癌恶性程度增加。ADC和T2两者均可以显示肿瘤内部的异质性，对早晚期宫颈癌进行鉴别。此部分研究结果显示，DR-CSI较常规ADC和T2在肿瘤诊断和侵袭性评估方面具有更大的优势。

本研究存在一定的局限性：第一，本试验样本量较少，未来研究将持续扩展样本量；第二，参考既往文献，为避免不同病理类型宫颈癌之间的差异影响，本研究仅选用宫颈鳞癌作为研究对象，未来还需纳入更多的宫颈癌病理学类型进行研究；第三，由于受试者没有接受离体宫颈癌MRI扫描，影像学数据无法与具体的微观组织结构详细关联分析，未来将致力于离体试验将DR-CSI微观区室与组织微结构相结合，发挥其更大的优势。