根据国际癌症研究机构发布的统计数据显示
[1],2022年乳腺癌的全球新发约230万例,在女性恶性肿瘤中位居前列,病死率占全球女性所有恶性肿瘤死亡病例的16.7%。腋窝淋巴结转移(axillary lymph node metastasis,ALNM)对预测乳腺癌患者总体复发率和生存率具有显著意义
[2-3]。目前,腋窝淋巴结(axillary lymph node,ALN)清扫或前哨淋巴结活检(sentinel lymph node biopsy,SLNB)是评估乳腺癌患者ALNM的金标准,但仍可能存在发生血管神经损伤、切口感染、淋巴水肿等并发症的风险,且因SLNB切除前哨淋巴结数量的差异,其假阴性率不容忽视
[4-5]。临床上评估ALN状态的影像学方法主要包括超声(ultrasound,US)、磁共振成像(magnetic resonance imaging,MRI)、钼靶X线(mammography,MG)等
[6]。US因其操作简便、可重复性强已成为乳腺癌及乳腺疾病最常用的检查方法
[7]。MG检查对乳腺微小钙化的检测具有高敏感度,对筛查早期乳腺癌具有重要意义
[8]。MRI可多方位显示肿块的位置及与深部组织的关系,其对软组织的分辨率更高,对操作者的依赖性更小
[9]。3种影像学手段各有优势、互为补充。当前,已有预测模型多依赖单一模态
[10]或2种模态的组合
[11],预测乳腺癌ALNM时均存在一定不足。而融合US、MRI、MG影像特征能够整合多模态影像的优势、克服单一模态的局限性,且进行多模态系统性整合分析的相关报道较少。多模态影像通过整合不同模态的图像特征,可为各系统疾病提供更全面、准确的信息,弥补单一模态数据的局限性,并且已经被广泛地用于不同疾病的诊断中
[12-13]。基于此,本研究拟通过分析多模态影像特征及临床病理参数,构建预测ALNM的列线图模型,并对其诊断性能进行验证,旨在为乳腺癌患者术前ALN评估提供个性化量化工具。
1 资料与方法
1.1 一般资料
选取2018年9月—2025年2月在哈尔滨医科大学肿瘤医院、山东省临沂市人民医院、哈尔滨医科大学附属第二医院就诊的女性乳腺癌患者332例,年龄24~81岁,平均(49.2±10.9)岁,分别收集患者术前US、MRI、MG图像及临床病理资料,读取多模态影像特征。将收集到的数据按约7∶3比例随机分成训练集(n=232)与验证集(n=100)。
纳入标准:①经病理学检查确诊为原发性乳腺癌,接受SLNB或ALN清扫,并获得了研究相关的临床及病理结果;②术前均接受乳腺MRI、MG及US检查;③影像学检查前均未进行乳腺手术、放疗及化疗等相关治疗。排除标准:①图像不清晰者;②已发生远处转移的患者。
本研究的伦理审查经哈尔滨医科大学附属第二医院审批通过,由于为回顾性及非侵入性研究,经过医学伦理委员会批准通过,免除受试者知情同意。
1.2 仪器与方法
1.2.1 检查仪器
US检查采用Philips Epiq 7、Philips Epiq 5超声诊断设备,配备5~15 MHz高频线阵探头,对患者双侧乳腺及腋窝区域进行系统扫查。MRI检查使用联影uMR 780 3.0T高强场MRI扫描仪进行动态增强磁共振成像检测,使用乳腺专用八通道相控阵线圈进行扫描。嘱患者俯卧位,双乳自然垂放于线圈中央。扫描方案包括平扫及动态增强扫描。先行平扫,获取T1WI及 脂肪抑制T2WI序列;随后进行动态增强扫描,经肘静脉注射对比剂钆喷酸葡胺(gadopentetate dimeglumine,Gd-DTPA)0.2 mmol/kg,速率2.5 mL/s,注射后立即追加生理盐水冲洗。对所获影像资料进行评估,在绘制病灶时间―信号强度曲线(time-signal intensity curve,TIC),选择病灶中最能体现恶性特征的区域作为感兴趣区,尽量避开病灶内坏死或囊变区。MG检查采用德国Senographe Pristina数字化乳腺X线摄影系统,在自动曝光控制条件下,对双侧乳腺执行内外斜位与头尾位标准投照。
1.2.2 多模态图像分析
依据BI-RADS诊断标准
[14]进行图像分析并记录多模态影像特征,包括①US检查记录肿块大小、形状、纵横比、边界、钙化、CDFI血流分级、长短径比值、淋巴门、皮质厚度、淋巴结形状;②MRI检查记录纤维腺体组织量、背景实质强化程度、肿块大小、形状、边界、强化特点、强化特征、TIC曲线类型、淋巴结短径、淋巴门结构;③MG检查记录病变类型、肿块形状、肿块边界、毛刺征、钙化。分别选US、MRI、MG专业各2名高年资医师(具有5年以上乳腺疾病诊断经验)在不知晓病理结果的情况下独立对乳腺肿块及ALN的图像进行分析,遇到分歧时共同协商解决。
1.2.3 临床病理资料
从医院电子信息系统中提取患者的临床资料(包括年龄、月经情况、淋巴结触诊结果、肿瘤所在象限)以及病理参数(包括病理类型、雌激素受体与孕激素受体状态、人类表皮生长因子受体2表达水平、Ki-67指数、分子分型以及是否存在脉管浸润)。
1.3 统计学方法
应用R 4.4.2统计软件分析数据。不符合正态分布的计量资料以[M(QR)]表示,比较采用Mann-Whitney U检验。计数资料以例数(%)表示,比较采用χ²检验或Fisher精确概率法。在R软件中设定随机种子(set.seed=999),采用随机分层抽样方法将总样本按7∶3比例分为训练集和验证集。在训练集中将单因素分析中P<0.05的变量纳入多因素Logistic回归分析,采用双向逐步回归法进行变量筛选,确定独立预测因素,构建预测ALNM的列线图模型,并在验证集中进行验证。在训练集与验证集中,通过绘制受试者工作特征(receiver operating characteristic,ROC)曲线计算曲线下面积(area under the curve,AUC),以评估模型的区分能力,Delong检验比较模型间AUC值差异的显著性;模型的校准度通过校准曲线评价,临床实用性通过决策曲线进行分析。P<0.05为差异有统计学意义。
2 结 果
2.1 训练集与验证集比较
332例患者中,训练集232例,验证集100例。2组US、MRI、MG特征及临床病理特征等比较差异均无统计学意义(
P均>0.05),说明2组所有变量具有良好的均衡性,有利于模型的稳定性验证,使模型的预测结果更具稳定性,见表
1~
4。病例展示见
图1。
2.2 构建列线图模型
2.2.1 基于多模态影像联合临床病理特征预测ALNM风险的列线图模型
对训练集中乳腺癌患者的多模态影像特征联合临床病理特征进行单因素及多因素逻辑回归分析,结果显示US-肿块边界不清晰、US-淋巴门消失、MG-存在毛刺征、MRI-无淋巴门及脉管浸润阳性与ALNM显著相关(
P<0.05),见
表5,是乳腺癌ALNM的独立影响因素。基于上述5个独立危险因素构建多模态影像―临床病理模型(multimodal imaging-clinicopathologic model,MICP),见
图2A。
2.2.2 基于多模态影像特征及临床资料预测ALNM风险的列线图模型
对训练集中乳腺癌患者的多模态影像特征及临床资料进行单因素及多因素逻辑回归分析,结果显示,US-肿块边界不清晰、US-淋巴门消失、MG-存在毛刺征、MG-有钙化及MRI-无淋巴门是乳腺癌ALNM的独立影响因素(
P<0.05),见
表6。基于上述5个独立危险因素构建模型多模态影像―临床模型(multimodal imaging-clinical,MIC),
图2B。
基于模型中各个预测变量对结局变量的贡献程度分别进行赋分,再将各个评分相加计算总评分,从而估测出该个体结局事件的预测概率。
2.3 评估列线图模型
MICP在训练集与验证集中的AUC值分别为0.870(95%
CI:0.824~0.917)和0.828(95%
CI:0.744~0.913);MIC在训练集与验证集中的AUC值分别为0.846(95%
CI:0.797~0.895)和0.807(95%
CI:0.715~0.899)(
图3A~B),列线图模型具有较好的区分度。训练集与验证集中,模型的校准曲线提示其对ALNM的预测概率与实际概率的一致性较好,均具有较好的拟合精度(
图3C~F)。临床决策曲线显示列线图模型在训练集与验证集中均获得了最优的临床净收益,具有较强的临床适用性(
图3G~H)。Delong检验结果显示MICP与MIC的AUC值差异无统计学意义(
Z=1.517,
P=0.129),见
图4。
3 讨 论
在我国,乳腺癌是严重威胁女性健康的主要恶性肿瘤之一
[15]。早期乳腺癌通常缺乏明显症状,难以被察觉
[16]。ALN状态始终是乳腺癌最为重要的预后因素,对于明确治疗决策有着极为重要的意义
[17]。乳腺影像学检查方法如MG、US、MRI等,在乳腺疾病的诊断与治疗中发挥至关重要的作用。本研究分析了术前US、MRI及MG的多模态影像特征与ALNM的相关性,并联合临床病理参数,构建列线图预测模型。通过评估模型诊断效能显示,能够较好地预测ALNM的风险,并能够直观的向临床医生展示,提高临床工作效率。
本研究基于US、MRI及MG的多模态影像及临床病理特征构建了预测ALNM的列线图模型MICP,显示原发灶边界不清、毛刺征、淋巴门结构异常以及脉管浸润是ALNM的独立预测因素,究其原因,本研究认为原发灶边缘区域聚集了大量的癌细胞,当癌细胞向周围正常组织浸润时,可导致影像学上出现特征性表现,US图像上多表现为肿块边界不清晰,而MG图像上可能呈现为向周围组织呈现锐利的放射状条带,即毛刺征,这些征象在一定程度上反映了癌细胞的侵袭能力
[18-20]。影像中出现边界不清和毛刺征时,表明癌细胞侵袭能力越强,并刺激病灶周围组织生成大量新生血管,为癌细胞经淋巴管扩散、转移至ALN提供条件
[21]。而在淋巴转移过程中,癌细胞在ALN内浸润、增殖,逐步破坏淋巴结内部正常结构,致使在US和MRI中表现为皮质不均匀增厚,严重者出现淋巴门结构不清晰甚至消失
[22-24]。既往研究
[25-26]表明,当出现淋巴门消失时常提示淋巴结受累程度较重,本研究结果与之一致。脉管浸润是指淋巴管、血管内能够检测到癌细胞,从病理机制角度分析,这为肿瘤细胞淋巴结转移提供了重要途径
[27-28]。已有研究
[29]表明脉管浸润对患者的总生存期、无病生存期及无转移生存期有着显著影响,因此其在ALNM预测中具有重要意义。本研究未显示Ki-67、Her-2与ALNM具有相关性,可能与患者的入组标准不同有关。
此外,由于脉管浸润需要通过粗针穿刺或手术切除标本获得
[30],属有创参数,临床应用具有一定的局限性。因此,本研究使用多模态影像特征及临床资料构建了预测乳腺癌ALNM的列线图模型MIC,最终证实US影像的边界特征及淋巴门消失、MG毛刺征及钙化、MRI淋巴门结构不清是ALNM的危险因素。训练集和验证集的AUC值分别为0.846和0.807,说明模型的预测准确性良好,能够较好的区分高风险人群与低风险人群。在临床实际应用中,部分患者可能因肿瘤位置不佳、身体状况差、凝血功能异常或穿刺取材受限等原因无法获得脉管浸润结果。此时,仅依据术前影像特征构建的MIC仍具有较好的预测效能,可作为有效的补充工具,辅助临床决策。近年来,诸多研究基于单独US、单独MRI、单独MG或两种方法联合应用构建乳腺癌ALNM预测模型。Wang等
[31]分析了162例三阴性乳腺癌的ALN的US特征,评估淋巴结良恶性,并获得了0.690的AUC值。陈洁文等
[32]基于507例早期乳腺癌患者的MRI特征及临床病理信息构建了预测模型,该模型区分有无ALNM的AUC值为0.762。方静等
[33]研究显示,结合术前超声造影、动态增强MRI定量参数及临床病理指标构建的多参数模型,对浸润性乳腺癌ALNM的诊断AUC为0.812。相较于以往研究,本研究纳入乳腺癌患者的术前US、MRI、MG多模态影像特征联合临床病理特征绘制成列线图模型MICP,结合三者的检查优势,可以弥补不同影像学方法的不足,从而提高模型的准确性,本研究结果显示MICP在预测ALNM风险方面表现出良好的判别效能,较现有临床模型明显提高(AUC=0.870)。此外,本研究使用多模态影像特征构建了预测乳腺癌ALNM的列线图模型MIC,其诊断效能较单一检查方法联合病理特征或两者方法联合病理特征的预测模型稍高,AUC为0.846。Delong检验结果表明两模型差异无统计学意义,2个模型的95%
CI高度重叠,进一步证实了MICP与MIC的预测效能相近。尽管本研究构建的多模态影像模型诊断性能均优于以往的临床模型,但较影像组学预测模型诊断效能较低。Yu等
[34]基于MRI构建影像组学模型用于识别早期乳腺癌ALNM的风险,训练组和验证组的AUC值分别为0.92和0.85。然而,由于影像组学技术仍需要大规模高质量影像数据训练模型,图像采集、预处理等环节仍需人工操作,且尚无统一标准,其临床转化仍面临一系列挑战。因此,与操作复杂的影像组学模型相比,基于多模态影像特征的列线图模型因其高效、直观、易于操作等特点在实际应用中更具优势。校准曲线及决策曲线显示,本研究构建的多模态影像模型具有良好的校准度和临床适用性,其预测乳腺癌ALNM能够为患者带来最大的临床净获益。综上,本模型通过整合多模态影像及临床病理特征,该模型能够辅助临床医生实现个性化综合评估,有望为ALN术前精准评估提供新思路。
然而本研究仍存在一定的局限性:①本研究为回顾性研究,样本量有限,存在选择偏倚;②人工阅片获取影像特征,可能存在一些主观误差;③本研究纳入分析的MRI、MG影像中淋巴结相关特征较少;④本研究在纳入数据时采用统一的纳排标准,划分训练集与验证集时未按中心进行分层随机分配,但随机分组后检验显示2组各变量无统计学差异,由于未采用中心分层随机化,中心间潜在的不可控差异仍可能存在。未来本研究将结合前瞻性、大样本数据以及多中心外部验证进一步完善和验证模型的稳健性。
综上所述,基于术前多模态影像特征联合临床病理参数构建的乳腺癌ALNM列线图预测模型,可以有效地评估乳腺癌ALN状态,协助临床医生辅助制定个体化治疗方案,避免乳腺癌患者的过度治疗及相关术后并发症,具有一定的临床应用潜力。
国家自然科学基金项目(82371984)
黑龙江省教育厅项目(2024-KYYWF-0221)