胰腺神经内分泌肿瘤近年来发病率持续升高,约55%~85%为无功能性胰腺神经内分泌肿瘤(nonfunctioning pancreatic neuroendocrine tumours, NF-PanNETs)
[1],具有异质性生物学行为特性,从惰性到高度侵袭性不等
[2],预后差异显著,5年生存率约为30%~66%
[3]。手术是目前治疗指南推荐的首选治疗方案,欧洲神经内分泌肿瘤学会共识指南建议>2 cm的胰腺神经内分泌肿瘤(pancreatic neuroendocrine tumours, PanNETs),应常规进行胰腺切除术联合区域淋巴结清扫术
[4]。而美国国家综合癌症网络指南建议>2 cm的肿瘤或淋巴结阳性患者施行根治性手术联合区域淋巴结清扫术,而对于1~2 cm的肿瘤,应考虑淋巴结切除术
[5]。淋巴结转移是NF-PanNETs预后的独立危险因素
[6-7]。此外,淋巴结转移的存在会增加术后复发的风险,显著降低NF-PanNETs患者的生存率
[8]。因此,术前淋巴结转移(lymph node metastatic, LNM)的准确评估对手术决策和预后预测至关重要。但目前术前对淋巴结转移的评估具有挑战性,使得手术淋巴结区清扫决策变得困难
[9]。目前术前评估淋巴结转移主要依赖于影像学征象,包括增强CT、MRI、超声等,然而,影像征象评估存在主观性,且在早期发现微小转移、鉴别反应性增生与转移性淋巴结时仍存在局限性
[7,10]。PET/CT特别是
68Ga-DOTATOC PET和
68Ga-DOTATATE PET通过靶向神经内分泌肿瘤细胞上高表达的生长抑素受体,能够检测原发肿瘤和转移灶,但该检查价格昂贵、辐射有害且仍存在敏感性、特异性问题
[11]。影像组学和深度学习技术在术前无创预测淋巴结转移中具有较大的潜力,但目前不同研究之间模型的可重复性差,这限制了模型的泛化性和临床适用性
[10]。亟需无创、准确、简便的术前预测无功能性胰腺神经内分泌患者淋巴结转移的方法。能谱CT是一种新兴的成像技术,可以重建出多参数图像并提供更多定量信息,在检测和诊断胰腺类疾病具有巨大的潜力
[12-13]。蓝燕芬等
[14]的研究指出,与常规CT相比,低能单能量虚拟图像结合无水碘图可以提高病灶的检出率。此外,已有研究验证能谱CT参数在胃癌、直肠癌预测淋巴结转移的能力
[15-16]。目前能谱CT参数对无功能性胰腺神经内分泌患者淋巴结转移预测效能尚未明确。因此,本研究旨在探讨能谱CT定量参数联合临床影像特征在术前预测NF-PanNETs淋巴结转移中的价值。
1 材料与方法
1.1 研究对象
本研究已获得中山大学附属第一医院伦理委员会批准通过并豁免了知情同意书(批件号:伦审[2021]721号)。回顾性分析2019年9月至2024年7月术前行能谱CT增强检查且有病理报告的胰腺无功能性神经内分泌肿瘤患者82例,其中淋巴结转移组30例、无淋巴结转移组52例。纳入标准如下:①经病理报告证实为神经内分泌肿瘤患者;②原发性病灶位于胰腺;③术前两星期内行能谱CT腹部增强扫描。排除标准:①术前经过局部或全身治疗;②有功能性胰腺神经内分泌肿瘤患者;③临床数据缺失;④影像数据不全或图像质量欠佳影响感兴趣区(region of interest, ROI)勾画。病人的纳入排除流程图详情请扫右侧附图阅读。
临床数据包括年龄、性别、身高、体质量、体质量指数(body mass index, BMI)、黄疸、腹痛、肿瘤大小、肿瘤位置、有无钙化、血管侵犯状态、器官侵犯状态、胰管状态、胆总管状态、肝转移状态、肿瘤分级、淋巴结状态、碳水化合物抗原(carbohydrate antigen,CA)19-9、CA125、癌胚抗原(carcinoembryonic antigen, CEA)、神经元特异性烯醇化酶(neuron-specific enolase, NSE)。
1.2 图像采集
所有患者均行能谱CT(IQon spectral CT,飞利浦)平扫及双期增强扫描,扫描前禁食至少4 h。扫描参数设置:管电压120 kVp;自动调制管电流:(166~255) mAs;螺旋螺距0.798;管旋转时间0.5 s;激光器宽度(64×0.625) mm;层厚1 mm;层间距0.8 mm。静脉注射碘对比剂(Ultravust 370,优维显)2.5 mL/s (用量:1.5 mL/kg,最大容量100 mL)。分别于注射对比剂后40 s和65 s获得胰腺实质期和门静脉期图像
[16]。患者平均剂量长度乘积为1 606 mGy × cm。
图像采集后,将由迭代重建技术(iDose4,3级,Philips Healthcare)的常规120 kV混合能量图像、光谱重建算法(spectral,level3,Philips Healthcare)的光谱基础数据导入专用后处理工作站(IntelliSpace Portal 10.0,Philips Healthcare)分别重建出胰腺实质期与门脉期40、70、100 keV单能量图像、碘密度图(iodine concentration, IC)、有效原子序数图像(Z-effective atomic number, Zeff)、虚拟平扫图像(virtual monoenergetic images, VMIs),重建层厚为1 mm,层间距为0.8 mm,与常规120 kVp混合能量图像参数且层面完全相同。
1.3 图像分析
由2名具有10年、15年腹部放射诊断经验的医师对120 kVp混合能量图像进行定性分析。若存在征象的分歧,此时2人与具有22年腹部放射诊断经验的医师共同进行讨论分析决定。定性分析的特征包括:①胰腺肿瘤的位置(头颈/体尾);②肿瘤的边缘(清晰/模糊);③是否存在钙化;④胰管是否存在扩张;⑤胆总管是否存在扩张;⑥是否存在器官侵犯;⑦是否存在血管侵犯;⑧是否存在肝转移;⑨是否存在淋巴结转移。肿瘤边缘清晰的定义为边界平滑且边界清晰可见。肿瘤边缘模糊的定义为瘤周超过90%出现毛刺或浸润
[17]。当CT值大于160 Hu定义为钙化,其中钙化分为点状钙化与粗大钙化(直径<5 mm为点状钙化、直径>5 mm为粗大钙化)。主胰管扩张定义为直径>3 mm。胆管扩张的定义为肝内胆管>2 mm或肝外胆管>8 mm。器官侵犯的定义为肿瘤与十二指肠、胃、脾、肾上腺、胆总管、结肠间正常脂肪间隙消失,且分界模糊不清
[18]。血管侵犯的判断标准为肿瘤血栓形成、血管闭塞、狭窄或轮廓畸形、血管周径超过半数与肿瘤接触。肝转移定义为肝内出现与胰腺原发灶强化方式相似的病灶
[19]。CT评估淋巴结转移定义为胰周或腹膜后淋巴结短径>10 mm或淋巴结内出现坏死
[20-22]。此外,肿瘤直径是基于CT轴位图像中肿瘤最大层面测量的肿瘤最大直径。
定量分析是基于120 kVp混合能量图像、40、70、100 keV虚拟单能量图像、无水碘图、有效原子序数图像、虚拟平扫图像,包括胰腺实质期和门静脉期。2名放射科医师独立勾画肿瘤感兴趣区域,并记录每个病人的测量数值。3个月后再安排其中一位医师勾画所有病人的肿瘤感兴趣区域并记录测量数值,使用平均值用于后续统计分析。ROI勾画原则:尽可能接近肿瘤的边缘,同时避开钙化、明显的血管以及大面积的坏死区域。复制ROI在其他相对应的图像上测量肿瘤的最大层面及上下层面,取3个ROI的平均值作为病灶的最终测量值。同时测量相同层面的主动脉相对应的值为后续计算使用,图像ROI勾画如
图1所示。
常规CT的定量参数包括胰腺实质期和门脉期病灶的CT值和标准化CT值。能谱CT的定量参数包括胰腺实质期和门脉期病灶的无水碘图、归一化无水碘图,以及40、70、100 keV虚拟单能量图像的CT值。在这些参数后附加字母“a”和“v”,分别表示胰腺实质期和门脉期。计算标准化指标的算式如下所示
[23]:
nCT = CTlesion/CTaorta
nVMI = VMIlesion/VMIaorta
nIC = IClesion/ICaorta
nZeff = Zefflesion/Zeffaorta
1.4 统计学分析
使用SPSS软件25.0版(SPSS Inc.,Chicago,IL,USA)和R语言(R Development Core Team)对所有数据进行统计学分析。计算组间和组内的相关系数以确保观察者间和观察者内定量参数的可重复性。采用卡方检验或Fisher精确检验对定性参数进行比较。根据Shapiro-Wilk检验结果,采用t检验或Mann-Whitney U检验对各参数进行定量参数的组间比较分析。单因素logistic回归有意义的指标采用最小绝对收缩和选择算子(least absolute shrinkage and selection operator,LASSO)回归分析筛选出独立相关参数。采用多因素logistic回归分析确定构建模型的参数。绘制受试者工作特征(receiver operating characteristic, ROC)曲线,计算曲线下面积(area under the curve, AUC)及计算95%置信区间(confidence interval, CI)、敏感度、特异度、准确度。P<0.05为差异有统计学意义。
2 结 果
本研究共纳入了82例无功能性胰腺神经内分泌肿瘤患者,其中转移组患者30例(平均年龄49.2±14.5,男性14例,女性16例),无转移患者52例(平均年龄50.0±14.9岁,男性23例,女性29例)。两组的临床参数与影像征象详见
表1。
2.1 临床指标
临床特征包括年龄、BMI、性别、肿瘤标志物CA19-9、CA125、CEA和NSE。除了2组间患者BMI的存在统计学差异,其余参数均无显著统计学差异。
2.2 CT定性模型
定性参数包括肿瘤边缘、器官侵犯、血管侵犯、肝转移、肿瘤分级、CT评估淋巴结转移的差异存在统计学意义。而组间肿瘤位置、腹痛、黄疸、钙化、胰管扩张、胆总管扩张的没有统计学意义。相比无转移组,转移组的肿瘤边缘多不清晰(77%
vs. 42%,
P=0.005),且更容易发生器官侵犯(57%
vs. 15%,
P<0.001),血管侵犯(53%
vs. 17%,
P=0.002),以及CT评估的淋巴结转移(47%
vs. 21%,
P=0.030),转移组的肿瘤高分级的更多。在肿瘤位置、腹痛、黄疸、钙化、胰管扩张、胆总管扩张方面,两组之间没有明显差异(
P>0.05)。单因素logistic分析(
表2)具有显著统计学差异的指标(BMI、肿瘤边缘、器官侵犯、血管侵犯、CT评估淋巴结转移)纳入LASSO回归分析(
图2),确定BMI、器官侵犯及血管侵犯为相关因素。多因素logistic回归分析显示BMI和器官侵犯是发生淋巴结转移的独立预测因素。基于临床指标和影像征象的模型(BMI+器官侵犯)AUC为0.815(95%CI:0.719~0.910;
表3),敏感度为0.780,特异度为0.769,准确度为0.780。
2.3 常规CT定量模型
常规CT定量参数包括了平扫CT值、胰腺实质期标准化的CT值(normalized CT value in the pancreatic phase,nCTa)、门静脉期标准化的CT值(normalized CT value in the portal-venous phase,nCTv)。单因素分析显示转移组CT值均低于无转移组,其中nCTa、nCTv存在统计学差异(均
P<0.05),而组间平扫CT值没有显著统计学差异。单因素分析显示nCTa、nCTv具有显著统计学差异的指标,经过多因素logistic回归确定nCTv为相关独立因素。基于nCTv的定量模型AUC为0.750(95%CI:0.642~0.858;
表3),敏感度为0.600,特异度为0.827,准确度为0.744。
2.4 能谱CT定量参数模型
能谱CT定量参数包括胰腺实质期和门静脉期的n40、n70、n100、nIC、nZeff。单因素分析显示淋巴结转移组能谱CT定量参数值均低于无转移组,具有显著统计学差异(所有
P<0.05)。使用LASSO回归分析胰腺实质期和门静脉期的能谱CT参数(n40、n70、n100、nIC、nZeff),确定门静脉期标准化碘浓度(normalized iodine concentration in the portal-venous phase,nICv)为淋巴结转移的相关因素。基于nICv的定量模型AUC为0.829(95%CI:0.733~0.925;
表3),敏感度为0.633,特异度为0.942,准确度为0.829。
2.5 综合模型
各模型的4个参数(BMI、器官侵犯、nCTv和nICv)采用多因素Logistic回归分析,其中器官侵犯和nICv纳入建立综合模型进行分析。基于器官侵犯和nICv的综合模型AUC为0.879(95%CI:0.806~0.953;
表3),敏感度为0.867,特异度为0.788,准确度为0.817。
2.6 诊断性能比较
4个模型ROC曲线如
图3所示。根据DeLong测试结果显示能谱CT定量参数模型和综合模型诊断效率高于常规CT定量模型(
P=0.018及
P=0.004),但能谱CT定量参数模型和综合模型之间没有观察到显著统计学差异(
P>0.05)。
3 讨 论
本研究构建了术前无功能性胰腺神经内分泌肿瘤淋巴结转移的预测模型。对临床指标、影像征象、常规CT及能谱CT参数进行分析,筛选出2个能预测淋巴结转移的独立危险因素:器官侵犯、nICv。基于这2个因素构建综合预测模型AUC为0.879,具有较好的预测性能,优于常规CT参数构建的预测模型。
手术切除是无功能性神经内分泌肿瘤的首选治疗方法,术前预测淋巴结转移对术中是否进行淋巴结区清扫的决策以及预后预测至关重要。低风险患者淋巴结转移率仅9.5%,可考虑不清扫;高风险患者淋巴结转移率达到41.2%,需要积极清扫
[9]。因此,若能术前准确预测淋巴结转移可以减少低风险患者过度治疗。CT是胰腺神经内分泌肿瘤的常规的检查方法
[24]。双层光谱探测器CT是CT成像领域的先进技术,可直接生成多参数图像,为胰腺疾病诊断提供更丰富的定量信息
[12]。
本研究分别构建单独临床参数模型(器官侵犯+BMI,AUC=0.815)、常规CT参数模型(nCTv,AUC=0.750)、能谱CT参数模型(nICv,AUC=0.829)和临床参数与能谱CT参数构建的综合预测模型(器官侵犯+nICv,AUC=0.879),并证实能谱CT参数模型的预测效能优于常规CT参数模型。将来若能发掘更多有意义的能谱CT新参数并由此开发新模型可能进一步提升模型的AUC值。有研究报道采用常规CT影像组学参数与临床参数建立的预测模型AUC=0.80,但影像组学需要收集大通量纹理数据进行分析筛选,其过程相对复杂,花费时间较长
[25],本研究采用能谱CT参数与临床参数结合的综合模型不仅能取得更好的表现,而且数据处理更简便。值得注意的是,本研究中能谱CT参数模型与临床参数模型表现相当,区别在于临床特征(器官侵犯)来源于诊断医生对CT影像的定性描述,能谱CT参数为定量数据,能谱参数更具有客观、定量、泛用的优势。能谱CT参数模型的表现优于常规CT参数模型,且能谱CT一次扫描可以同时获得常规与能谱CT定量参数供诊断医生进行分析,由此可见能谱CT为临床评估淋巴结转移提供极大的便利,展现较高的临床应用价值。
本研究发现能谱CT中门静脉期nICv是淋巴结转移的独立影响因素,也是最重要的定量预测因子。其中淋巴结转移组各能谱CT参数值均低于无淋巴结转移组,可能与分组中神经内分泌肿瘤的分级构成有关。既往研究报道G3级PanNETs与G1/G2级具有不同的CT增强强化特征:G3级肿瘤增强强化较弱
[26],高分级肿瘤强化程度显著低于G1/G2级
[27-30],组织病理学分析结果显示高级别PanNETs肿瘤血管生成程度低于低级别病例
[31]。肿瘤的强化特征与新生血管形成及血管密度有关
[32]。此前有研究将G1、G2、G3分类成高分化NF性神经内分泌肿瘤
[33],表现为过度表达血管内皮生长因子受体和血小板衍生生长因子受体通路,导致新生血管形成增加及微血管密度升高
[34]。相反,低分化NF性神经内分泌肿瘤(神经内分泌癌)则表现出微血管密度降低
[34-35]。IC因对碘具有高敏感性,能有助于显示病灶血供情况
[36],与本研究观察的结果一致。本研究中转移组中高级别肿瘤占比较无转移组更高(无转移组:G1和G2占95%,G3占6%,转移组:G1和G2占67%,G3和NEC占33%),转移组nIC值均较无转移组低(PnICa=0.007,PnICv<0.001),侧面也验证了以往的研究观点
[26-30]。
此外,本研究中还发现器官侵犯是预测NF-PanNETs的独立影响因素,与Harisinghani等
[37]研究发现相同。低级别NETs通常细胞分化良好且形态均一,侵袭程度较低。相比之下,高级别NETs属于侵袭性肿瘤,细胞分化程度低,且通常无功能活性
[34]。转移组中高级别NF-PanNETs的占比较无转移组多(转移组:G3和NEC占33%,无转移组:G3和NEC占6%),与既往研究报道
[34,37]一致,表明高级别NF-PanNETs侵袭程度高,容易出现器官侵犯,是预测淋巴结转移的重要因素。
本研究仍存在一定局限性。首先,本研究是回顾性设计,不可避免会存在病例纳入偏倚。其次,本研究纳入数据量有限,有可能会存在统计误差,但本研究中采用LASSO回归分析,尤其适用于小样本研究筛选变量建模,一定程度上减少模型的不可靠性;同时,研究中在变量筛选阶段仅纳入具有明确病理生理意义的指标,减少冗余信息干扰,并采用多因素logistic回归进行独立危险因素验证,仅保留具有统计学意义指标,提升模型简约性。此外,本研究中受限于回顾性设计,未能分析淋巴结的能谱CT参数特征与病理结果相关性,后续将进一步优化研究方案,积极开展前瞻性研究,以期为临床提供更为精准的淋巴结转移的影像学依据。最后,本研究为单中心、单机器研究,缺乏多中心及多机器外部数据验证模型性能。后续将补充多中心、多机器数据,进一步扩大样本量,验证能谱CT参数预测NF-PanNETs淋巴结转移的效能。
综上所述,本研究证明了能谱CT定量参数能在术前无创、准确预测NF-PanNETs淋巴结转移,能谱CT定量参数结合临床影像特征可以提升诊断效能,且优于常规CT定量参数,为NF-PanNETs手术治疗决策及预后预测提供影像支持,改善患者预后。
国家自然科学基金(82472096)
广东省自然科学基金(2024A1515011968)