术前MRI在预测乳腺癌预后中的应用及展望

边小倩，杜思瑶，张立娜

2022 年美国癌症协会的统计数据显示，乳腺癌是女性发病率排名首位且致死率排名前5 的肿瘤[1]，其前5 年的复发风险为10.4%，后5年的复发风险为4.5%[2]。目前，临床已证实腋窝淋巴结转移、脉管侵犯、激素受体表达、组织学分级、分子亚型等是乳腺癌复发风险及预后的重要影响因素[3-6]。但以上指标只能通过活检或术后病理获得，存在侵入性和不可重复性。MRI 具有无创、易获得和可重复的特点，近年来已被广泛应用于乳腺癌预后相关研究，其部分影像特征及参数包括形态学、血流动力学、功能成像、影像组学等被证实为预测乳腺癌患者预后的成像标志物[7-8]，为患者制定更优化的治疗策略提供重要支持。

1 基于MRI的形态学特征

MRI 形态学特征是反映病灶基本属性的常见指标。其中肿瘤大小、边缘、非肿块样强化、边缘强化、瘤周水肿、背景强化等是目前已知的预后相关影像参数。

临床已经证实大于2 cm的肿瘤更易复发，基于MRI测量的肿瘤大小与术后病理大小一致性良好，复发组明显大于非复发组(平均大小3.9 cm vs. 2.7 cm,P=0.02)[9]。肿块不规则伴毛刺边缘更多见于三阴性亚型、存在脉管侵犯的侵袭性病例中，往往预后较差[10-12]。相对肿块，非肿块样强化表现出更差的预后[13-14]，这可能与非肿块样病变难以早期诊断且阳性切缘有关[15]。边缘强化是肿瘤外围对比于中心的增强[13]，被认为是由肿瘤内微血管密度降低引起的，而非肿瘤周围微血管密度的增加[16]。在三阴性亚组的多变量分析中，发现仅边缘增强与复发之间存在显著相关性，结果显示边缘增强的肿瘤复发率是非边缘增强的14 倍(HR=14.019；95%CI：1.773～110.864)[13,17]。瘤周水肿的存在是乳腺癌复发的独立危险因素(HR=2.77，P=0.022)[18]，在T2 加权成像(T2-weighted images，T2WI)中，不同的水肿类型(包括瘤周水肿、胸前水肿、皮下水肿和弥漫性水肿)与肿瘤侵袭性呈现不同程度的相关性[19]，而三阴性乳腺癌出现瘤周水肿的几率更大[20]，说明瘤周水肿与高侵袭性肿瘤相关，是瘤周脉管侵犯和炎症的间接指标。在目前的指南中，瘤周水肿并不是乳腺癌的关键预后因素，大数据高质量的相关研究可能会推动指南的未来更新。

除了原发灶及瘤周的特征，乳腺背景实质强化(background parenchymal enhancement，BPE)对乳腺癌预后也具有一定的预测价值。同侧乳腺的BPE 可能因乳腺癌的存在而受到血管化增加的影响，导致BPE的MRI评估假性抬高[21-22]，因此已提出对侧BPE的MRI评估作为改进乳腺癌决策的工具。根据乳腺成像报告和数据系统的BPE模式分类，在绝经后乳腺癌患者中，高BPE等级表现出更差的无复发生存结果(HR=3.086，P=0.003)[23]。但BPE的评估受主观因素、年龄和激素水平影响较大，因此，需要更多影响因素控制完善的相关数据以及量化技术的开发进一步证实。

综上，传统MRI 形态学特征受临床医生的主观判断及临床经验的影响较大，且大多为定性特征，因此对乳腺癌预后的价值相对受限，通过人工智能，如深度学习(deep learning，DL)全自动病灶分割，提高其可重复性，将传统的形态学指标进行定量化将是乳腺MRI的主流和趋势。

2 基于MRI的功能学参数

MRI 功能学参数是乳腺影像研究的热点之一。目前常用的功能学指标包括动态增强、弥散及衍生序列所产生的相关参数。尽管还存在着一些争议，系列研究已经证实了功能学参数在预测乳腺癌预后中的重要价值。

2.1 血流动力学

动态对比增强磁共振成像(dynamic contrast-enhanced magnetic resonance imaging，DCE-MRI)是评估血管功能最常用的方法，包括基于时间—强度曲线(time-intensity curve，TIC)的半定量参数，如：达峰时间(time-to-peak enhancement，Tpeak)、最大强化斜率(maximum rise slope，Slopemax)等和基于药代动力学计算模型的定量参数，如：转运常数(transfer constant，Ktrans)、回流常数(rate constant，Kep)和细胞外血管外体积分数(extracellular extravascular volume fraction，Ve)等。

基于TIC 计算得到的半定量参数中，不同研究的结果各异。相对预后良好的患者，肿瘤内较快的初始增强[24]、更高的峰值增强[25]，TIC 的Slopemax[26]、更短的Tpeak[27]以及较多的流出成分[28]是较差生存结局的重要预测参数。Niukkanen 等[29]利用3D 分割技术分析得到瘤内、瘤周的信号增强比(signal enhancement ratio，SER)与肿瘤大小、组织学分级、Ki-67 表达相关，进一步证实了瘤周区域TIC 相关的半定量参数与乳腺癌不良预后的关联性。

研究表明，Ktrans、Kep与雌激素受体(estrogen receptor，ER)、孕激素受体(progesterone receptor，PR)表达呈负相关；Ktrans与组织学分级、人表皮生长因子受体2 (the human epidermal growth factor receptor 2，HER-2)表达、淋巴结转移及Ki-67表达呈正相关；Kep与HER-2、Ki-67 表达呈正相关[30-31]，这可能与新生毛细血管增加导致对比剂流出及回流的速度加快有关。基于DCE-MRI 定量参数的直方图分析结果显示，高侵袭性乳腺癌(三阴性乳腺癌)、高Ki-67表达(＞20%)和高核级癌症(2级或3 级)表现出更高的Ve值变异系数和偏度，且具有显著差异性[32]，而这些癌症被认为具有高侵袭及高增殖性，往往预后不良。

虽然TIC形态与复杂的定量灌注参数显著相关[33]，但定量参数显然具有更好的客观性和可比性。基于药代动力学模型定量参数的精准计算需要依赖高时间分辨率的DCE-MRI序列扫描，并受不同设备、参数设置和所选计算模型的影响。但随着乳腺癌定量研究的不断深入，基于DCE-MRI影像定量参数的乳腺癌侵袭性或预后仍然是未来研究的重要趋势。随着磁共振波谱成像(magnetic resonance spectroscopic imaging，MRSI)、化学交换饱和转移成像技术(chemical exchange saturation transfer，CEST)、血氧水平依赖(blood oxygenation level-dependent，BOLD)、弛豫时间等更多定量扫描序列在乳腺领域的逐步应用，功能学定量参数不断丰富，我们对多参数定量模型的开发充满期待。

2.2 弥散加权成像及衍生技术

弥散加权成像(diffusion-weighted imaging，DWI)已在临床实践中常规使用，先进的DWI技术如体素内不相干运动(intravoxel incoherent motion，IVIM)、弥散张量成像(diffusion tensor imaging，DTI)、弥散峰度成像(diffusion kurtosis imaging，DKI)等在乳腺癌预后评估中也表现出潜在价值[34-36]。

传统的DWI 因出现时间早、扫描序列和后处理相对简便易行，已被广泛用于乳腺癌预后的研究中。肿瘤内表观弥散系数(apparent diffusion coefficient，ADC)差异值越大(＞0.698×10-3mm2/s)，发生远处转移的风险越高(HR=4.5，P＜0.001)[37]。经过7.2 年的平均随访时间，瘤周与瘤内ADC 比率较高的患者相对比率较低的患者具有更差的总生存率[38-39]，这可能与肿瘤内异质性及肿瘤侵袭性诱导的瘤周液体渗出相关[40]。IVIM 实现体素内水分子弥散和灌注成分的分离，反映组织细胞的数量和血管分布[41]。但它与乳腺癌预后之间的关系存在争议[34,42-43]。有研究[42-43]表明，其参数伪扩散系数与激素受体、Ki-67表达、组织学分级之间具有相关性，Song等[34]却认为IVIM 参数与预后因素之间并无显著关联，原因可能是定量参数的计算及拟合b 值的选取存在差异。DTI 获得体素内的各向异性程度及弥散情况，以表征乳腺癌微观结构的差异。先前的研究认为较大的肿瘤(＞2 cm)、高组织学分级(3级)、腋窝淋巴结转移与较低的平均弥散率(mean diffusivity，MD)显著相关，而淋巴结状态与各向异性分数(fractional anisotropy，FA)差异并无统计学意义[35,44]；ER 阴性和Ki-67 高表达与FA 值呈显著负相关[45]，这与肿瘤细胞密度增加、微观结构破坏，使得水分子弥散的幅度和方向降低相关[45]。DKI 用于研究b 值超过1000 s/mm2时水分子的非高斯弥散。组织学分级和Ki-67 表达被认为与平均峰度(mean kurtosis，MK)呈正相关，与MD 呈负相关，但与肿瘤大小无明显相关性[46-47]。高侵袭性肿瘤的MK值升高，反映了组织微观结构较复杂，偏离高斯分布较大，它增加了单体素中的细胞密度，进而影响水分子弥散[48]。但也有关于组织学分级的阴性研究[36]，其结果的差异可能是由b值选取、扫描参数不同等多种因素造成。

DWI 及其衍生技术在乳腺癌预后方面显示出巨大的预测潜力，但DWI 序列在乳腺扫描中也存在一些不足：信噪比和空间分辨率较低、空间失真较大等，会导致图像损失一些重要信息。目前DWI 序列仍作为DCE-MRI 的补充序列。DTI 和IVIM技术更详细地表征了组织微观结构，似乎更具优越性，但是其扫描和计算模型的复杂性更容易产生不一致的结果[34-35]。近年来，弥散相关序列扫描技术不断进步，有望弥补现有序列的缺陷，并将DWI 及衍生技术与人工智能相结合，以挖掘更多有价值的预后信息。

3 基于MRI的影像组学和人工智能

影像组学和人工智能是近几年来在各个领域科学研究的重点和热点，同样也应用在乳腺癌预后预测的相关研究中，并已经取得了很多乐观、积极的结果。

影像组学特征已被证明是预测乳腺癌预后的独立生物标志物[49]。一项基于肿瘤异质性纹理分析的研究发现，T2WI 的高熵值(≥6.013，HR=9.84)和动态增强T1 加权减影图像的低熵值(＜5.057，HR=4.55)与较差的无复发生存率显著相关[50]。虽然纹理分析可以在一定程度上测量整个肿瘤内的异质性，但这种测量依赖于肿瘤内异质性的良好混合，忽略了肿瘤内的区域表型变异[51]。而亚区分割技术更加关注肿瘤内的灌注异质性，相对于临床、影像组学等其他4 种模型，利用组学特征构建联合多亚区的复发风险评估模型表现出最好的预测性能(C-Index=0.760)[52]，因此，灌注异质性的量化是预测乳腺癌预后的一种潜在方法。然而，Yoon 等[53]从ADC 图中提取的纹理参数均未显示出对疾病无进展生存期的显著预测价值。造成结果不同的原因，可能是不同的研究在参数选择、量化标准、亚型组成及计算方法等方面存在差异。

目前，针对乳腺癌预后进行预测的机器学习(machine learning，ML)模型主要使用MRI 组学特征融合临床特征进行训练，并能达到满意的预测效能[3]。当使用包含多参数MRI 特征的8 种ML 算法(如：随机森林、决策树、K-最邻近等)对乳腺癌的重要预后因素Ki-67 和组织学分级进行评估时，6 个模型显示出对Ki-67 的等效性能，受试者工作特征(receiver operating characteristic，ROC)曲线下面积(area under the curve，AUC)为0.70；对组织学分级的预测，贝叶斯算法表现最佳(AUC=0.79)[54]。ML 允许提取内容丰富的成像信息，且能够量化人眼无法察觉的组织之间的差异，因此基于MRI 影像组学的ML在乳腺癌预后方面有望取得更多进展。

卷积神经网络(convolutional neural network，CNN)是目前最流行的用于图像分析的DL 结构[55]。近年来，乳腺癌预后预测相关的DL方法不断深入：Liu等[56]在DCE-MRI对比前后图像生成的肿瘤掩模基础上改进了3D-CNN 框架，认为由肿瘤掩模引导的3D-CNN 预测模型比掩模引导的完整图像或仅掩模体素的图像具有更高的5 年无复发生存预测精度。相比传统的3D-CNN 方法[57]，由掩模引导的预测模型能将注意力集中在与掩模不同的乳腺肿瘤内未知区域，更强调了与肿瘤预后相关的图像特征[56]。4D-CNN (随时间变化的3D-CNN)模型增加了时间维度，补充了2D 和3D 所缺乏的时空背景，在融合临床信息后，对腋窝淋巴结转移的预测效能及稳定性均有较大提升，也进一步证实了肿瘤和瘤周信息对预测乳腺癌转移的重要价值[58]。DL 已经在乳腺癌预后预测中显示出重要价值，其无创、可预知的特点有望将患者利益实现最大化。

由于缺乏明确的生物学原理来解释组学特征和预后之间的联系，影像组学的临床应用价值受到限制。同时，由于图像采集方法的不一致和计算方法选择的差异，影像组学研究结果的难以重复性也是其在各领域普遍处于瓶颈期的原因之一。尽管ML和DL对乳腺肿瘤的分割、诊断及预后等方面产生了深远影响，但其预测结果的低可解释性同样阻碍了他们的临床应用[59-60]，当前影像组学研究和人工智能相关研究正在致力于准确性、科学性和泛化能力，探索模型对预测结果的解释能力。另外在给予人工智能更大自主权的同时，对不能提供理由的黑箱决策尝试进行监督，尊重公众对人工智能决策的看法[60]。

4 结论和展望

综上所述，随着MRI 技术的不断开发与应用，术前对乳腺癌患者的预后进行无创预测的设想成为可能，并及时为预后不佳的患者提供恰当的辅助治疗，对提高患者生存率具有重要意义。除了已建立的乳腺癌预后相关的MRI 形态学、功能学特征外，基于MRI 的影像组学进一步揭示了更多与预后相关的高维度参数，计算机引导的人工智能正在兴起。然而，由于乳腺癌的异质性和影像表现的复杂性，未来在乳腺癌筛查、诊断、治疗和预后方面的研究还有很长的路要走，大队列、前瞻性、多中心研究值得进一步探索。

作者利益冲突声明：全体作者均声明无利益冲突。