基于双模态MRI影像组学术前预测浸润性乳腺癌腋窝淋巴结转移

徐 敏，马宜传，张书海，王祥芝，汤晓敏，杨 丽，刘 浩，谢宗玉

乳腺癌是女性最常见恶性肿瘤之一，其发病率在女性恶性肿瘤中一直高居首位；在我国，乳腺癌死亡率高居女性恶性肿瘤致死率的第2位[1]。腋窝淋巴结转移(axillary lymph node metastasis，ALNM)不仅是乳腺癌预后的重要影响因素之一，而且与病人的手术方案制定有密切联系[2]。因此，术前准确判断乳腺癌病人有无ALNM具有重要临床意义。然而，目前常用的方式均为有创方式，包括腋窝淋巴结清扫术(axillary lymph node dissection，ALND)和前哨淋巴结清扫术(sentinel lymph node dissection，SLND)，这些方式可能具有一定的并发症风险，如神经和血管损伤、感染、血肿等[3-4]。由于MRI的无创性、无电离辐射、多序列成像、高软组织分辨率等优势，其在术前无创评估ALNM中具有巨大的应用潜力。然而既往研究[5-6]显示单一的MRI序列在评估ALNM中的能力有限。影像组学可以将数字医学图像转换为可挖掘的数据，分析数据，提取特征，并利用这些特征来直观、定量地描述病灶的生物学特点，全面评估肿瘤的异质性，从而提高诊断、预后和预测准确性，近年来被广泛应用于辅助疾病诊断、预测肿瘤病人的治疗疗效、预后等方面[7]。因此，本文旨在将一种基于双模态MRI的影像组学用于术前预测浸润性乳腺癌ALNM，以协助乳腺癌病人治疗方式的制定，改善病人预后。

1 对象与方法

1.1 研究对象 收集蚌埠医学院第一附属医院2018年11月至2020年11月期间，经手术病理证实的209例浸润性乳腺癌病人。纳入标准：(1)经病理证实为浸润性乳腺癌，且明确ALNM情况；(2)病人在肿瘤切除或活检前均行MRI检查。排除标准：(1)临床病理资料或MRI检查资料不完整；(2)图像质量差无法完成有效图像分割。最终纳入168例浸润性乳腺癌病人进行回顾性分析研究，64例有ALNM(ALN+)，104例无ALNM(ALN-)；病人均为女性，年龄23～68(48.67 ± 9.12)岁。

1.2 设备与扫描参数 采用Philips Achieva 3.0 T磁共振及SENSE 7通道矩阵乳腺专用线圈；主要序列及扫描参数如下：横断位T1WI序列(TR 400 ms，TE 10 ms，层厚5 mm，层间距1 mm，FOV 350×350，矩阵216×194，激励次数2，FA 90°)；横断位T2WI_SPAIR序列(TR 5 000 ms，TE 60 ms，层厚5 mm，层间距1 mm，FOV 350×350，矩阵260×274，激励次数1，FA 90°)；动态增强扫描(dynamic contrast enhanced magnetic resonance imaging，DCE-MRI)序列(TI 90 ms，TR 4.8 ms，TE 2.4 ms，层厚5 mm，层间距1 mm，FOV 350×350，矩阵252×339，激励次数1，FA 10°)，第1次扫描即预扫后注射钆对比剂(Gd-DTPA，0.2 mmol/kg，2 mL/s)，接着注射20 mL 0.9%氯化钠溶液(1.5 mL/s)，完成注射同时开始连续无间断扫描5次，每次扫描时间60 s，总共6期。

1.3 影像组学分析

1.3.1 图像分割 将所有病人的图像导入医准-达尔文科研平台(http://www.yizhun-ai.com)。分别在T2WI和DCE两个序列图像上，采用手动勾画感兴趣区(region of interest，ROI)的方式进行图像分割，在T2WI序列图像上选取肿瘤最大层面进行勾画，DCE上，根据时间-信号强度(TIC)曲线选取最高峰值信号的期相进行勾画，避开血管和坏死(见图1、2)。所有图像均由2名有丰富经验的影像医师共同进行勾画分割图像。根据分割图像后提取的组学特征采用相关系数(intraclass correlation coefficient，ICC)进行组间一致性分析。

1.3.2 特征提取、选择与模型构建 T2WI和DCE两个序列图像ROI各提取出939个特征，包括一阶特征、形状特征、纹理特征及小波特征，共1 878个影像组学特征，随机将病人按8∶2分为训练组(n=134)和验证组(n=34)。在训练组中，运用最小最大值归一化将每一维度特征线性拉伸至(0，1)之间，再使用Select K best、最小绝对收缩和选择算子(least absolute shrinkage and selection operator，LASSO )回归算法及迭代筛选特征对这些高维组学特征中进行降维，筛选出与浸润性乳腺ALNM高度相关的影像组学特征。利用logistic回归构建浸润性乳腺癌ALNM的影像组学预测模型，分别为T2WI、DCE和T2WI联合DCE预测模型。采用受试者工作特征(receiver operating characteristic，ROC)曲线分析每个模型在训练组和验证组中的表现，通过曲线下面积(area under curve，AUC)评价模型的预测效能，以最优模型生成影像组学列线图(nomogram)，使预测模型能更直观，有效地应用于临床。

1.4 统计学方法 采用t检验和χ2检验。

2 结果

2.1 各组病人临床资料的比较 淋巴结转移组病人的肿瘤长径大于无淋巴结转移组(P<0.01)，年龄和病理分级差异无统计学意义(P>0.05)；年龄、肿瘤长径及病理分级在训练组和验证组病人之间差异均无统计学意义(P>0.05)(见表1)。

表1 不同组间临床病理资料的比较

2.2 特征提取、选择和预测模型构建 MRI图像中提取的影像组学特征共1 878个。用最小最大值归一化、Select K best和LASSO回归算法及迭代筛选特征进行特征筛选，筛选出与浸润性乳腺癌ALNM相关的影像组学特征：T2WI序列4个特征(1个一阶特征、1个形状特征、1个纹理和1个小波特征)；DCE-MRI序列4个特征(1个形状特征、1个纹理特征和2个小波特征)；T2WI联合DCE序列5个特征(1个一阶特、1个形状特征、1个纹理特征和2个小波特征)。运用logistic回归构建的T2WI、DCE和T2WI联合DCE影像组学预测模型在训练组中的AUC和95%CI分别为0.75(0.63，0.88)、0.75(0.63，0.86)、0.80(0.68，0.92)；在验证组中分别为0.75(0.51，1.00)、0.73(0.44，0.99)、0.79(0.58，1.00)；T2WI联合DCE的AUC在训练组及验证组中均最高，高于单独T2WI和单独DCE预测模型，T2WI联合DCE预测模型在预测浸润性乳腺癌ALNM中预测效能最佳(见图3)。依据最优模型生成列线图(见图4)，该列线图可以实现浸润性乳腺癌病人发生ALNM的个体化预测，计算出的数值越高，病人发生ALNM的可能性越高。

3 讨论

乳腺癌病人在治疗前获得ALNM状态，对临床分期、治疗方案选择和预后评估至关重要[8]。乳腺MR现虽已广泛用于术前评估乳腺癌病人的ALNM状态，但乳腺表面线圈区域大小和腋窝空气界面的磁敏感伪影均限制了MRI对腋窝淋巴结的直接评估，而且单一序列图像的评估效能有一定局限性[9]。影像组学将影像图像资料数据化并深度挖掘、分析，可以显示出肉眼难以识别的肿瘤组织解剖学的微小变化，从而达到精准诊断[10]。目前，影像组学方法已被报道用于肺癌、结直肠癌、膀胱癌等恶性肿瘤术前淋巴结转移的预测[11-13]，然而对乳腺癌术前ALNM的预测尚属于初步探索阶段。因此，本研究聚焦于在MR图像上分割乳腺癌原发病灶提取影像组学特征预测ALNM状态。

本研究中ALNM组64例，无ALNM组104例，2组临床病理资料的病灶长径差异有统计学意义 (P<0.01)，与既往研究[14-15]结果一致；目前，病理分级与ALNM状态的关系仍有很多争论，本研究中病理分级与ALNM状态的差异无统计学意义(P>0.05)，与CHAI等[16]研究结果一致，而与LIU等[17]研究结果不一致，这可能与样本中组织学分型及浸润性乳腺癌的临床分期有关系。本研究临床资料在训练组和验证组间差异无统计学意义，这也反映了本研究随机分组的合理性。

logistic回归简便易行，一般用于二分类预测模型构建，也可用于多分类预测模型，是目前影像组学研究中较为常用的建模方法[18]，故本研究使用logistics回归基于T2WI、DCE-MRI提取的高度相关的影像组学特征，构建T2WI、DCE和T2WI联合DCE三个影像组学预测模型。T2WI加权图像可较好地反映正常组织结构和病灶形态，可以更准确地显示乳腺肿块。本研究中T2WI 影像组学预测模型的AUC(训练组0.75、验证组0.75)，对ALNM具有一定的预测效能，但低于TAN等[4](训练组0.880、验证组0.819)利用T2WI序列的放射组学特征预测ALNM中的研究(n=329)，这可能与本研究的样本量较少有关。DCE-MRI通过评价肿瘤的形态和血流动力学是诊断乳腺癌的有效手段，该方法能够提供高时间分辨率、高空间分辨率和高信噪比的图像。DCE-MRI有多期扫描期相，而目前选取特定的期相提取病灶特征能提供最好的预测效能还没有共识。本研究中DCE影像组学预测模型在训练组的AUC(0.75)低于LIU等[19]应用从原发肿瘤第一强化期提取的放射学特征预测前哨淋巴结转移，AUC值为0.806，该研究与本研究选取ROI方式相同，均选取最高强化峰值信号图像勾画ROI；然而也有研究[20]采用第二期图像进行勾画ROI，而AUC(0.86)高于上述研究，因此关于DCE-MRI图像勾画ROI的最佳期相选择有待于更多的研究来证实。本研究以T2WI联合DCE作为最优模型建立影像组学列线图，训练组AUC为0.80，验证组为0.79，在预测ALNM方面体现了较好的预测价值；MAO等[21]基于对比度增强光谱乳腺X线摄影建立影像组学列线图预测乳腺癌ALNM，验证组AUC为0.79；QIU等[22]使用了基于超声成像提取的21个纹理特征来预测乳腺癌淋巴结转移，验证组的AUC为0.759。本研究中,基于双模态MRI影像组学特征构建的列线图显示相当或更好的性能，并且可在术前通过无创方式获得，在预测淋巴结转移方面具有可接受的价值。

本研究的不足之处：(1)本研究为样本量尚小的回顾性研究，可能存在不可避免的选择性偏倚，且为单中心实验，后续仍需大样本、多中心的实验进一步验证。(2)本研究只分析了肿瘤图像的最大层面，特征提取不算全面；(3)放射科医师手工勾画乳腺癌病灶的ROI轮廓，主观经验影响较大。(4)本研究因DWI图像空间分辨率较低未纳入研究，旨在未来优化图像质量后，将DWI纳入研究。

综上，基于T2WI和DCE-MRI提取原发病灶的影像组学特征对浸润性乳腺癌ALNM具有较好的预测效能，可更好地协助临床、指导治疗。