基于MRI DWI的影像组学模型对骨肉瘤新辅助化疗疗效的评估价值

于荭 张胜男 段丽莎 孔蒙娟 崔建岭

引用本文:
Citation:

基于MRI DWI的影像组学模型对骨肉瘤新辅助化疗疗效的评估价值

    通讯作者: 崔建岭, jianlingcui@sina.com

Evaluation of the efficacy of neoadjuvant chemotherapy in osteosarcoma based on MRI DWI radiomics model

    Corresponding author: Jianling Cui, jianlingcui@sina.com
  • 摘要: 目的 探讨基于MRI弥散加权成像(DWI)的影像组学模型对骨肉瘤新辅助化疗疗效的评估价值。 方法 回顾性分析河北医科大学第三医院2015年6月至2017年11月经术后组织病理学检查结果证实,且在接受新辅助化疗前、后均行MRI检查的41例骨肉瘤患者[男性26例、女性15例,年龄(22.0±11.0)岁,范围11~49岁]的病历及影像资料。根据术后组织病理学检查结果,将肿瘤组织坏死率≥90%者纳入疗效好组,<90%者纳入疗效差组。分别于新辅助化疗前、化疗一期结束后5 d内和完成整个化疗后测量所有患者的表观扩散系数(ADC,分别记为ADC0、ADC1、ADC2,比较疗效好组和疗效差组ADC间的差异。在化疗一期结束后的DWI(b=1000 s/mm2)和ADC图像上手动勾画病灶的感兴趣区,提取影像组学特征,用随机分组法将数据按6∶4的比例分为训练集和测试集,采用SMOTE算法对训练集上的数据进行扩充,采用方差阈值、SelectKBest、最小绝对收缩和选择算子(LASSO)法进行影像组学特征筛选,采用逻辑回归分类器构建出影像组学模型。采用独立样本t检验或Wilcoxon秩和检验进行2组间比较。采用受试者工作特征(ROC)曲线评估传统影像学(ADC)及影像组学模型对骨肉瘤新辅助化疗疗效的预测效能。 结果 疗效好组10例、疗效差组31例。2组患者ADC0的差异无统计学意义[(0.95±0.05)×10−3 mm2/s对(1.05±0.05)×10−3 mm2/s,t=1.14,P>0.05)];疗效好组的ADC1、ADC2高于疗效差组,且差异均有统计学意义[(1.44±0.10)×10−3 mm2/s对(1.10±0.06)×10−3 mm2/s,t=−2.92,P<0.05;1.68(1.55,1.85)×10−3 mm2/s对(1.33±0.06)×10−3 mm2/s,Z=−2.61,P<0.01]。ROC曲线分析结果显示,当ADC1≥1.34×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为80%,特异度为81%,曲线下面积(AUC)为0.797(95%CI:0.629~0.965);当ADC2≥1.51×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为90%,特异度为71%,AUC为0.777(95%CI:0.588~0.967)。从化疗一期结束后的DWI和ADC图像中共提取出1409个影像组学特征,按6∶4的比例随机分为训练集和测试集[24(疗效好:6,疗效坏:18)∶17(疗效好:4,疗效坏:13)],将训练集数据扩充为70(疗效好:20,疗效坏:50),经影像组学特征筛选后,最终得到5个最优影像组学特征,包括Interquartile Range、Skewness、Uniformity、Median、Maximum。采用逻辑回归分类器构建影像组学模型,训练集中该模型预测骨肉瘤新辅助化疗疗效的ROC曲线的AUC为0.881(95%CI:0.811~0.942),灵敏度为90%,特异度为74%;测试集中AUC为0.769(95%CI:0.515~0.933),灵敏度为75%,特异度为69%。 结论 基于MRI DWI的影像组学模型在评估骨肉瘤新辅助化疗疗效中的效能优于传统影像学(ADC),在临床应用中潜力较大。
  • 图 1  骨肉瘤患者弥散加权成像(A,b=1000 s/mm2)和表观扩散系数图像(B)上勾画的病灶感兴趣区

    Figure 1.  Regions of interest delineated on diffusion weighted imaging (A, b=1000 s/mm2) and apparent diffusion coefficient (B) images of osteosarcoma patient

    图 2  ADC预测骨肉瘤患者新辅助化疗疗效的受试者工作特征曲线

    Figure 2.  Receiver operating characteristic curves of apparent diffusion coefficient values in prediction of the efficacy of neoadjuvant chemotherapy for osteosarcoma

    图 3  影像组学特征的均方误差随α参数变化趋势图(A)及影像组学降维图(B)

    Figure 3.  Mean square error of radiomics features change trend according to α parameters (A) and dimension reduction diagram of radiomics features (B)

    图 4  影像组学模型预测骨肉瘤患者新辅助化疗疗效的受试者工作特征曲线

    Figure 4.  Receiver operating characteristic curves of radiomics features in evaluation of the treatment response in osteosarcoma after neoadjuvant chemotherapy

    表 1  新辅助化疗后不同疗效的2组骨肉瘤患者的ADC比 较(×10−3 mm2/s)

    Table 1.  Comparison of apparent diffusion coefficient values between two groups of osteosarcoma patients with different therapeutic effects after neoadjuvant chemotherapy (×10−3 mm2/s)

    组别ADC0
    $\bar x\pm s $
    ADC1
    $\bar x\pm s $
    ADC2
    [MQ1Q3)或$\bar x\pm s $]
    疗效好组(n=10)0.95±0.051.44±0.101.68(1.55, 1.85)
    疗效差组(n=31)1.05±0.051.10±0.061.33±0.06
    检验值1.14a−2.92a−2.61b
    P0.2600.0060.009
    注:a表示t值;b表示Z值。ADC为表观扩散系数;ADC0为新辅助化疗前的ADC;ADC1为化疗一期结束后5 d内的ADC;ADC2为完成整个化疗后的ADC
    下载: 导出CSV
  • [1] Sadykova LR, Ntekim AI, Muyangwa-Semenova M, et al. Epidemiology and risk factors of osteosarcoma[J]. Cancer Invest, 2020, 38(5): 259−269. DOI: 10.1080/07357907.2020.1768401.
    [2] Jundt G. Updates to the WHO classification of bone tumours[J]. Pathologe, 2018, 39(2): 107−116. DOI: 10.1007/s00292-017-0396-4.
    [3] Manfredi R, Maresca G, Smaniotto D, et al. Cervical cancer response to neoadjuvant therapy: MR imaging assessment[J]. Radiology, 1998, 209: 819−824. DOI: 10.1148/radiology.209.3.9844681.
    [4] Lim HJ, Ong CAJ, Tan JWS, et al. Utility of positron emission tomography/computed tomography (PET/CT) imaging in the evaluation of sarcomas: a systematic review[J]. Crit Rev Oncol Hematol, 2019, 143: 1−13. DOI: 10.1016/j.critrevonc.2019.07.002.
    [5] Hao YW, An R, Xue YS, et al. Prognostic value of tumoral and peritumoral magnetic resonance parameters in osteosarcoma patients for monitoring chemotherapy response[J]. Eur Radiol, 2021, 31(5): 3518−3529. DOI: 10.1007/s00330-020-07338-y.
    [6] Oh C, Bishop MW, Cho SY, et al. 18F-FDG PET/CT in the management of osteosarcoma[J]. J Nucl Med, 2023, 64(6): 842−851. DOI: 10.2967/jnumed.123.265592.
    [7] Wang JF, Sun ML, Liu DW, et al. Correlation between apparent diffusion coefficient and histopathology subtypes of osteosarcoma after neoadjuvant chemotherapy[J]. Acta Radiol, 2017, 58(8): 971−976. DOI: 10.1177/0284185116678276.
    [8] Yu H, Gao L, Shi RQ, et al. Monitoring early responses to neoadjuvant chemotherapy and the factors affecting neoadjuvant chemotherapy responses in primary osteosarcoma[J]. Quant Imaging Med Surg, 2023, 13(6): 3716−3725. DOI: 10.21037/qims-22-1095.
    [9] Degnan AJ, Chung CY, Shah AJ. Quantitative diffusion-weighted magnetic resonance imaging assessment of chemotherapy treatment response of pediatric osteosarcoma and Ewing sarcoma malignant bone tumors[J]. Clin Imaging, 2017, 47: 9−13. DOI: 10.1016/j.clinimag.2017.08.003.
    [10] Habre C, Dabadie A, Loundou AD, et al. Diffusion-weighted imaging in differentiating mid-course responders to chemotherapy for long-bone osteosarcoma compared to the histologic response: an update[J]. Pediatr Radiol, 2021, 51(9): 1714−1723. DOI: 10.1007/s00247-021-05037-4.
    [11] 王芳, 夏雨薇, 柴象飞, 等. 影像组学分析流程及临床应用的研究进展[J]. 中华解剖与临床杂志, 2021, 26(2): 236−241. DOI: 10.3760/cma.j.cn101202-20200701-00200.
    Wang F, Xia YW, Chai XF, et al. Analysis process and clinical application of radiomics[J]. Chin J Anat Clin, 2021, 26(2): 236−241. DOI: 10.3760/cma.j.cn101202-20200701-00200.
    [12] He LJ, Yang HN, Huang JS. The tumor immune microenvironment and immune-related signature predict the chemotherapy response in patients with osteosarcoma[J/OL]. BMC Cancer, 2021, 21(1): 581[2022-11-22]. https://bmccancer.biomedcentral.com/articles/10.1186/s12885-021-08328-z. DOI: 10.1186/s12885-021-08328-z.
    [13] Prabowo Y, Setiawan I, Kamal AF, et al. Correlation between prognostic factors and the histopathological response to neoadjuvant chemotherapy in osteosarcoma: a retrospective study[J]. Int J Surg Oncol, 2021, 2021: 8843325. DOI: 10.1155/2021/8843325.
    [14] Ritter J, Bielack SS. Osteosarcoma[J]. Ann Oncol, 2010, 21(Suppl 7): Svii320−vii325. DOI: 10.1093/annonc/mdq276.
    [15] Kerwin WS, O'Brien KD, Ferguson MS, et al. Inflammation in carotid atherosclerotic plaque: a dynamic contrast-enhanced MR imaging study[J]. Radiology, 2006, 241(2): 459−468. DOI: 10.1148/radiol.2412051336.
    [16] Laux CJ, Berzaczy G, Weber M, et al. Tumour response of osteosarcoma to neoadjuvant chemotherapy evaluated by magnetic resonance imaging as prognostic factor for outcome[J]. Int Orthop, 2015, 39(1): 97−104. DOI: 10.1007/s00264-014-2606-5.
    [17] Choi SH, Jung SC, Kim KW, et al. Perfusion MRI as the predictive/prognostic and pharmacodynamic biomarkers in recurrent malignant glioma treated with bevacizumab: a systematic review and a time-to-event meta-analysis[J]. J Neurooncol, 2016, 128(2): 185−194. DOI: 10.1007/s11060-016-2102-4.
    [18] Zeitoun R, Shokry AM, Khaleel SA, et al. Osteosarcoma subtypes: magnetic resonance and quantitative diffusion weighted imaging criteria[J]. J Egypt Natl Canc Inst, 2018, 30(1): 39−44. DOI: 10.1016/j.jnci.2018.01.006.
    [19] Liu HH, Zhang CY, Wang LJ, et al. MRI radiomics analysis for predicting preoperative synchronous distant metastasis in patients with rectal cancer[J]. Eur Radiol, 2019, 29(8): 4418−4426. DOI: 10.1007/s00330-018-5802-7.
    [20] Kumar V, Gu YH, Basu S, et al. Radiomics: the process and the challenges[J]. Magn Reson Imaging, 2012, 30(9): 1234−1248. DOI: 10.1016/j.mri.2012.06.010.
    [21] Kayal EB, Kandasamy D, Khare K, et al. Texture analysis for chemotherapy response evaluation in osteosarcoma using MR imaging[J]. NMR Biomed, 2021, 34(2): e4426. DOI: 10.1002/nbm.4426.
    [22] 苗文杰, 杨光杰, 聂佩, 等. 高分辨率CT影像组学联合传统影像学征象预测肺腺癌微血管浸润的价值[J]. 国际放射医学核医学杂志, 2020, 44(9): 541−547. DOI: 10.3760/cma.j.cn121381-201909006-00075.
    Miao WJ, Yang GJ, Nie P, et al. Value of HRCT radiomics combined with traditional imaging features in predicting microvascular invasion of lung adenocarcinoma[J]. Int J Radiat Med Nucl Med, 2020, 44(9): 541−547. DOI: 10.3760/cma.j.cn121381-201909006-00075.
    [23] Lee SK, Jee WH, Jung CK, et al. Prediction of poor responders to neoadjuvant chemotherapy in patients with osteosarcoma: additive value of diffusion-weighted MRI including volumetric analysis to standard MRI at 3T[J/OL]. PLoS One, 2020, 15(3): e0229983[2022-11-22]. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0229983. DOI: 10.1371/journal.pone.0229983.
  • [1] 石琴张文付志明卢晓娟贝金玲 . 肿块型与非肿块型乳腺癌MRI弥散加权成像表观扩散系数与分子生物学指标的相关性研究. 国际放射医学核医学杂志, 2020, 44(4): 225-230. doi: 10.3760/cma.j.cn121381-201903008-00019
    [2] 高平王茜赵赟赟岳明纲李河北李原99Tcm-MDP SPECT/CT显像用于骨肉瘤肺转移的诊断. 国际放射医学核医学杂志, 2014, 38(3): 148-151. doi: 10.3760/cma.j.issn.1673-4114.2014.03.002
    [3] 陈伟君孙达99Tcm-MIBI显像在乳腺癌新辅助化疗中的应用价值. 国际放射医学核医学杂志, 2015, 39(6): 487-492. doi: 10.3760/cma.j.issn.1673-4114.2015.06.011
    [4] 李雯冯彦林18F-FDG PET/CT评估三阴性乳腺癌新辅助化疗疗效的研究进展. 国际放射医学核医学杂志, 2014, 38(3): 197-201. doi: 10.3760/cma.j.issn.1673-4114.2014.03.013
    [5] 刘建井边海曼马文娟王子阳陈薇朱磊戴东徐文贵 . 基于18F-FDG PET/CT图像建立的多元影像组学模型对乳腺癌原发灶HER-2表达状态的预测价值. 国际放射医学核医学杂志, 2022, 46(7): 430-440. doi: 10.3760/cma.j.cn121381-202203015-00206
    [6] 曹胜男舒政陈艾琪卫军杜小萌马宜传 . 基于骰骨MRI的糖尿病足影像组学特征. 国际放射医学核医学杂志, 2021, 45(12): 759-766. doi: 10.3760/cma.j.cn121381-202011036-00122
    [7] 侯震李双双闫婧 . 影像组学与深度学习在肿瘤放疗中的研究进展. 国际放射医学核医学杂志, 2020, 44(12): 744-749. doi: 10.3760/cma.j.cn121381-202004004-00107
    [8] 苗文杰杨光杰聂佩颜蕾郭柄州龚爱迪赵钰鋆崔景景王振光 . 高分辨率CT影像组学联合传统影像学征象预测肺腺癌微血管浸润的价值. 国际放射医学核医学杂志, 2020, 44(9): 541-547. doi: 10.3760/cma.j.cn121381-201909006-00075
    [9] 陈旭唐光才王富林韩福刚 . DWI联合增强MRI对肝硬化背景下小肝癌诊断价值的Meta分析. 国际放射医学核医学杂志, 2019, 43(2): 152-159. doi: 10.3760/cma.j.issn.1673-4114.2019.02.010
    [10] 高爽王三刚 . 急性期脑梗死患者DWI与MRA诊断的一致性及影响因素分析. 国际放射医学核医学杂志, 2024, 48(2): 82-90. doi: 10.3760/cma.j.cn121381-202212011-00348
    [11] 潘登于丽娟18F-FDG PET/CT代谢参数与非小细胞肺癌病理类型相关性的研究进展. 国际放射医学核医学杂志, 2021, 45(8): 527-531. doi: 10.3760/cma.j.cn121381-202009004-00082
  • 加载中
图(4)表(1)
计量
  • 文章访问数:  2779
  • HTML全文浏览量:  2244
  • PDF下载量:  8
出版历程
  • 收稿日期:  2022-11-23
  • 网络出版日期:  2023-09-04
  • 刊出日期:  2023-09-25

基于MRI DWI的影像组学模型对骨肉瘤新辅助化疗疗效的评估价值

    通讯作者: 崔建岭, jianlingcui@sina.com
  • 1. 河北医科大学第三医院CT/MR室,河北省骨科生物力学重点实验室,石家庄 050011
  • 2. 天津市天津医院放射科,天津 300299

摘要:  目的 探讨基于MRI弥散加权成像(DWI)的影像组学模型对骨肉瘤新辅助化疗疗效的评估价值。 方法 回顾性分析河北医科大学第三医院2015年6月至2017年11月经术后组织病理学检查结果证实,且在接受新辅助化疗前、后均行MRI检查的41例骨肉瘤患者[男性26例、女性15例,年龄(22.0±11.0)岁,范围11~49岁]的病历及影像资料。根据术后组织病理学检查结果,将肿瘤组织坏死率≥90%者纳入疗效好组,<90%者纳入疗效差组。分别于新辅助化疗前、化疗一期结束后5 d内和完成整个化疗后测量所有患者的表观扩散系数(ADC,分别记为ADC0、ADC1、ADC2,比较疗效好组和疗效差组ADC间的差异。在化疗一期结束后的DWI(b=1000 s/mm2)和ADC图像上手动勾画病灶的感兴趣区,提取影像组学特征,用随机分组法将数据按6∶4的比例分为训练集和测试集,采用SMOTE算法对训练集上的数据进行扩充,采用方差阈值、SelectKBest、最小绝对收缩和选择算子(LASSO)法进行影像组学特征筛选,采用逻辑回归分类器构建出影像组学模型。采用独立样本t检验或Wilcoxon秩和检验进行2组间比较。采用受试者工作特征(ROC)曲线评估传统影像学(ADC)及影像组学模型对骨肉瘤新辅助化疗疗效的预测效能。 结果 疗效好组10例、疗效差组31例。2组患者ADC0的差异无统计学意义[(0.95±0.05)×10−3 mm2/s对(1.05±0.05)×10−3 mm2/s,t=1.14,P>0.05)];疗效好组的ADC1、ADC2高于疗效差组,且差异均有统计学意义[(1.44±0.10)×10−3 mm2/s对(1.10±0.06)×10−3 mm2/s,t=−2.92,P<0.05;1.68(1.55,1.85)×10−3 mm2/s对(1.33±0.06)×10−3 mm2/s,Z=−2.61,P<0.01]。ROC曲线分析结果显示,当ADC1≥1.34×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为80%,特异度为81%,曲线下面积(AUC)为0.797(95%CI:0.629~0.965);当ADC2≥1.51×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为90%,特异度为71%,AUC为0.777(95%CI:0.588~0.967)。从化疗一期结束后的DWI和ADC图像中共提取出1409个影像组学特征,按6∶4的比例随机分为训练集和测试集[24(疗效好:6,疗效坏:18)∶17(疗效好:4,疗效坏:13)],将训练集数据扩充为70(疗效好:20,疗效坏:50),经影像组学特征筛选后,最终得到5个最优影像组学特征,包括Interquartile Range、Skewness、Uniformity、Median、Maximum。采用逻辑回归分类器构建影像组学模型,训练集中该模型预测骨肉瘤新辅助化疗疗效的ROC曲线的AUC为0.881(95%CI:0.811~0.942),灵敏度为90%,特异度为74%;测试集中AUC为0.769(95%CI:0.515~0.933),灵敏度为75%,特异度为69%。 结论 基于MRI DWI的影像组学模型在评估骨肉瘤新辅助化疗疗效中的效能优于传统影像学(ADC),在临床应用中潜力较大。

English Abstract

  • 骨肉瘤是青少年最常见的原发性恶性骨肿瘤,其致死和致残率高[1]。目前新辅助化疗、手术加术后化疗使骨肉瘤患者的5年生存率提高到了70%~80%[2]。新辅助化疗后骨肉瘤的肿瘤组织坏死率被认为是评价预后的“金标准”,而预后的好坏是确定是否可以行局部切除保肢治疗的依据。然而,这个“金标准”只能在术后获得,无法用于判断保肢手术进行与否。目前,用于术前评价新辅助化疗疗效的方法主要有临床评价和传统影像学评价,前者依据肿瘤大小变化、疼痛评分等化疗前后的症状或体征进行评价;后者主要包括对肿瘤边界、大小等形态学变化进行评价。这些方法虽然有一定的参考价值,但缺乏量化、标准化的指标,不能准确反映化疗疗效[3]。影像学评价方法具有很大的价值,其包括形态学评估及功能学评估,例如动态对比增强(DCE)-MRI、18F-FDG PET/CT等[4-6],但这些检查目前缺乏特异性高的标准化参数。MRI弥散加权成像(diffusion weighted imaging, DWI)功能显像方法简单、方便、无创,近年来的很多研究结果显示,表观扩散系数(apparent diffusion coefficient,ADC)在反映骨肉瘤等恶性肿瘤新辅助化疗疗效中具有重要价值[7-10]。但是由于肿瘤的异质性很大,其结果并不稳定,无法满足精准医学的要求。影像组学是将影像数据定量分析与机器学习相结合,筛选出具有诊断价值的影像组学特征,建立诊断模型,为临床诊疗提供更准确的信息[11]。本研究探讨基于MRI DWI的影像组学模型对骨肉瘤新辅助化疗疗效评估的价值。

    • 回顾性分析2015年6月至2017年11月在河北医科大学第三医院接受治疗的41例骨肉瘤患者的病历及影像资料,其中男性26例、女性15例,年龄(22.0±11.0)岁,范围11~49岁。纳入标准:(1)经穿刺或手术活检的组织病理学检查结果证实为原发性骨肉瘤;(2)手术前接受新辅助化疗,且化疗前行常规MRI及DWI检查。排除标准:(1)未做手术,无法获得肿瘤组织坏死率;(2)MRI图像伪影重,影响诊断。手术切除肿瘤后,根据Huvos分级系统[12]对所有标本进行组织学分析,将肿瘤组织坏死率≥90%者纳入疗效好组,<90%者纳入疗效差组。本研究已申请并通过了河北医科大学第三医院医学伦理委员会批准(批准号:2022-059-1),且所有患者及其家属在进行MRI检查前均签署了安全调查表及知情同意书。

    • 所有患者采用相同的新辅助化疗方案。每组化疗药物包含:(1)顺铂20 mg/m2×5 d+吡柔比星20~30 mg/m2×3 d;(2)异环磷酰胺2 g/m2×5 d;(3)甲氨蝶呤8~12 g/m2×1 d。每组化疗药用药结束为一个化疗周期,每个化疗周期结束后间隔2~3周开始用下一组药,第3组化疗药用药结束后2~3周进行手术。分别于新辅助化疗前、化疗一期结束后5 d内和完成整个化疗后测量所有患者的ADC,分别记为ADC0、ADC1、ADC2。

    • 应用德国Siemens公司Verio 3.0T超导型MRI仪进行扫描。根据肿瘤位置的不同,患者可采取仰卧位头先进或仰卧位足先进。根据扫描部位选取不同的表面线圈,如八通道膝关节线圈或八通道体部柔软线圈。常规行T1加权轴位显像(重复时间 700 ms,回波时间 12 ms);质子加权成像脂肪抑制序列(PDWI-FS)横断面、矢状面、冠状面显像(重复时间4000 ms,回波时间 83 ms);DWI采用单次激发自旋回波平面回波序列(SE-EPI),扫描参数:重复时间5800 ms,回波时间100 ms,视野10 mm×210 mm,矩阵192×192,层厚3 mm,层间距0.75 mm,层数20,体素1.1 mm×1.1 mm×3.0 mm,b值为0、500、1000 s/mm2

    • 将入组患者各期的原始图像调入德国Siemens公司Verio Workstation 3.0后处理工作站,由具有3年和5年临床工作经验的2名放射科医师在未知临床和组织学信息的情况下,在b=1000 s/mm2的DWI图像和ADC图像上,沿肿瘤边缘勾画病灶ROI(图1),并在ADC图像上测量病灶ADC,测量时避开肿瘤内出血、钙化及液化坏死区域。

      图  1  骨肉瘤患者弥散加权成像(A,b=1000 s/mm2)和表观扩散系数图像(B)上勾画的病灶感兴趣区

      Figure 1.  Regions of interest delineated on diffusion weighted imaging (A, b=1000 s/mm2) and apparent diffusion coefficient (B) images of osteosarcoma patient

    • 在Python 3.8.5(http://www.rodiomics.io/pyradio-mics.html)中使用PyRadiomics 3.0.1工具软件,对化疗一期结束后DWI和ADC图像上勾画的病灶ROI进行影像组学特征的提取,提取的影像组学特征分为以下4组:(1)一阶统计量特征,由126个描述符组成,用来描述MRI图像中voxel强化度的分布,包括最大值、最小值、中值、均值等;(2)空间几何特征,包含14个反映区域形状和大小的三维特征,包括体积、表面积、球形度、最大直径等;(3)纹理特征,从灰色水平运行长度和灰色水平共发生纹理矩阵计算,有525个可以量化区域异质性差异的纹理特征,包括灰度共生矩阵、灰度梯度矩阵、领域灰度差分矩阵等;(4)变换特征,包括对数变换、指数变换、小波变换、平方根变换等。将提取的影像组学特征用随机分组法按6∶4的比例分为训练集和测试集,采用SMOTE算法对训练集上的数据进行扩充。

    • 首先,将上述获取的所有特征分别进行组内相关系数(intraclass correlation coefficient,ICC)分析,选择ICC>0.75的特征;然后,采用方差分析(ANOVA)进行特征的初步筛选,选择P<0.05的特征后,采用方差阈值、SelectKBest、最小绝对收缩和选择算子(LASSO)法进行降维,采用LASSO法得到偏差和方差的最佳平衡点,剔除冗杂、错误和无相关性的特征值,筛选出最有价值的特征,采用逻辑回归分类器构建影像组学模型。

    • 应用SPSS 22.0软件进行统计学分析。采用Shapiro-Wilk 检验法对计量资料进行正态性检验,符合正态分布的计量资料以$\bar x\pm s $表示,呈偏态分布的计量资料以MQ1Q3)表示。符合正态分布且方差齐的数据采用独立样本t检验进行2组间比较;不符合正态分布或方差不齐的数据采用两独立样本Wilcoxon秩和检验法。绘制传统影像学(ADC)预测骨肉瘤患者新辅助化疗疗效的ROC曲线,评估其诊断效能。绘制影像组学模型的ROC曲线,计算AUC、灵敏度、特异度,评估影像组学模型对骨肉瘤新辅助化疗疗效的预测效能。采用ICC评价2名医师提取的病灶ROI的影像组学特征的一致性:ICC<0.40表示一致性较差;ICC为0.40~0.59表示一致性一般;ICC为0.60~0.74表示一致性较好;ICC>0.75表示一致性良好。P<0.05为差异有统计学意义。

    • 41例骨肉瘤患者中,疗效好组10例、疗效差组31例。由表1可知,2组患者的ADC0的差异无统计学意义(t=1.14,P>0.05);疗效好组的ADC1、ADC2高于疗效差组,且差异均有统计学意义(t=−2.92,P<0.05;Z=−2.61,P<0.01)。ROC曲线分析结果显示,当ADC1≥1.34×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为80%,特异度为81%,AUC为0.797(95%CI:0.629~0.965);当ADC2≥1.51×10−3 mm2/s时,其对骨肉瘤新辅助化疗疗效评估的灵敏度为90%,特异度为71%,AUC为0.777(95%CI:0.588~0.967)(图2)。

      组别ADC0
      $\bar x\pm s $
      ADC1
      $\bar x\pm s $
      ADC2
      [MQ1Q3)或$\bar x\pm s $]
      疗效好组(n=10)0.95±0.051.44±0.101.68(1.55, 1.85)
      疗效差组(n=31)1.05±0.051.10±0.061.33±0.06
      检验值1.14a−2.92a−2.61b
      P0.2600.0060.009
      注:a表示t值;b表示Z值。ADC为表观扩散系数;ADC0为新辅助化疗前的ADC;ADC1为化疗一期结束后5 d内的ADC;ADC2为完成整个化疗后的ADC

      表 1  新辅助化疗后不同疗效的2组骨肉瘤患者的ADC比 较(×10−3 mm2/s)

      Table 1.  Comparison of apparent diffusion coefficient values between two groups of osteosarcoma patients with different therapeutic effects after neoadjuvant chemotherapy (×10−3 mm2/s)

      图  2  ADC预测骨肉瘤患者新辅助化疗疗效的受试者工作特征曲线

      Figure 2.  Receiver operating characteristic curves of apparent diffusion coefficient values in prediction of the efficacy of neoadjuvant chemotherapy for osteosarcoma

      从化疗一期结束后的DWI和ADC图像中共提取出1409个影像组学特征,按6∶4的比例随机分为训练集和测试集[24(疗效好:6,疗效坏:18)∶17(疗效好:4,疗效坏:13)],将训练集数据扩充为70(疗效好:20,疗效坏:50),经筛选降维后获得5个最优影像组学特征,分别为Interquartile Range、Skewness、Uniformity、Median、Maximum,均为纹理特征,如图3所示,确定了LASSO法的最优α参数为−log(α)=2.32。以最优α参数对应的5个特征及其系数分别为:original_firstorder_Interquartile Range(系数为−0.289894)、original_firstorder_Skewness(系数为−0.409912)、original_firstorder_Uniformity(系数为0.015519)、original_firstorder_Median(系数为−0.373091)、original_firstorder_Maximum(系数为0.811615)。

      图  3  影像组学特征的均方误差随α参数变化趋势图(A)及影像组学降维图(B)

      Figure 3.  Mean square error of radiomics features change trend according to α parameters (A) and dimension reduction diagram of radiomics features (B)

      使用此影像组学模型预测骨肉瘤患者新辅助化疗疗效,结果显示,训练集中,影像组学模型预测骨肉瘤患者新辅助化疗疗效的ROC曲线的AUC为0.881(95%CI:0.811~0.942) ,灵敏度为90%,特异度为74%,测试集中AUC为0.769(95%CI:0.515~0.933),灵敏度为75%,特异度为69%(图4)。

      图  4  影像组学模型预测骨肉瘤患者新辅助化疗疗效的受试者工作特征曲线

      Figure 4.  Receiver operating characteristic curves of radiomics features in evaluation of the treatment response in osteosarcoma after neoadjuvant chemotherapy

    • 骨肉瘤的恶性程度高,预后差,虽然近年来新辅助化疗和手术切除的综合治疗方法显著提高了骨肉瘤患者的5年生存率[2],但仍有30%~40%的患者出现局部复发或远处转移,导致5年生存率显著降低,仅为 23%~29%[13],这可能与部分骨肉瘤患者对化疗不敏感有关,另一个重要原因是术前化疗效果不理想。术前对骨肉瘤化疗疗效的准确评估至关重要,不仅对化疗方案的及时调整具有指导意义,还能在一定程度上预测患者预后[8]。目前,采用 Huvos 分级系统进行化疗疗效评估,可在一定程度上预测复发风险,但该分级系统是基于术后切除标本的肿瘤组织坏死率来进行评估的[12],具有滞后性,无法为是否可以进行保肢治疗提供有力依据,也无法很好地指导术后化疗。传统影像检查如X射线、CT和MRI等可根据肿瘤体积变化、肿瘤内部液化坏死及钙化情况、周围骨质破坏等形态学变化对化疗效果进行评估[14],但无法反映肿瘤细胞的存活情况,与肿瘤组织坏死率的相关性较差,在骨肉瘤的分级、疗效及预后评估中具有明显的局限性[15],且相关研究结果显示肿瘤体积变化与化疗疗效间并无明显相关性[16]。动态对比增强(DCE)-MRI、DWI等能够反映恶性肿瘤血管新生情况及肿瘤细胞活性等信息[5,17-19],具有非常好的应用价值和发展前景,但由于肿瘤的异质性,目前仍缺乏特异性较高的定量指标。因此,为了能够早期判断肿瘤的化疗敏感性、及时调整化疗及手术方案、预测术前化疗疗效及预后,临床迫切需要一种简便、准确、可重复的新辅助化疗疗效评估方法。

      Kumar等[20]完善了影像组学的概念,高通量地从MRI、PET及CT等常规医学影像中提取出大量高维的定量影像特征,并进行进一步分析,从而能够定量地描述影像中的时间和空间异质性,揭示出肉眼无法识别的图像特征。影像组学研究的主要流程包括图像获取、ROI分割、影像组学特征的提取、特征值选择和降维、预测模型的性能评估5个步骤[21]。通过提取整合常规医学影像中肉眼无法识别的海量数字化信息,并经过机器学习算法、建立模型,有效避免了观察者的主观差异性,具有量化分析肿瘤异质性的功能[22],对疾病的精准诊疗具有重大的价值。并且已有研究基于MRI DWI构建影像组学模型,对骨肉瘤化疗反应及预后进行预测,如Lee等[23]的研究选择DWI和ADC图像特征构建模型,结果显示构建的模型可提高对骨肉瘤患者接受新辅助化疗后反应不良的诊断准确性。影像组学的最终目的是使建立的模型不仅对现有的数据有很好的预测能力,而且对未知的数据也有很好的预测能力。作为一种新兴的技术,影像组学在临床中的应用价值已在越来越多的研究中得到证实。

      本研究利用影像组学的方法,初步探索了基于MRI DWI的影像组学模型对骨肉瘤新辅助化疗疗效的评估效能。对所有提取到的影像组学特征进行统计分析、降维后,最终得到5个具有预测价值的特征,构建影像组学模型。训练集中,AUC为0.881(95%CI:0.811~0.942) ,灵敏度为90%,特异度为74%,测试集中AUC为0.769(95%CI:0.515~0.933),灵敏度为75%,特异度为69%,这表明本研究所获得的影像组学模型对于骨肉瘤新辅助化疗疗效具有较好的预测效能。

      本研究存在一定的局限性:(1)收集的骨肉瘤数据样本量较小且是单中心研究,缺乏对该模型在不同条件下的稳定性验证、前瞻性验证,为提高模型的诊断效能,今后尚需进行前瞻性研究、多中心合作,扩大样本量;(2)本研究采用人为方法勾画ROI,不仅工作量大、耗时久,不适用于大样本量的数据处理,而且由于肿瘤病灶形态欠规则,对病灶边界的确定存在个体差异性,这可能会对纹理特征筛选的准确性、模型的可重复性和稳定性造成影响,今后仍需进一步研究进行ROI自动分割的方法;(3) 本研究仅单一地分析了影像组学特征在预测骨肉瘤新辅助化疗疗效中的效能,未与其他常规影像学特征、组织病理学特征以及临床特征相结合,有待下一步研究和验证。

      综上所述,本研究结果表明,基于MRI DWI的影像组学模型对于骨肉瘤新辅助化疗疗效具有良好的预测效能,能够为骨肉瘤临床诊疗决策的制定提供有价值的指导信息,从而实现骨肉瘤的个体化精准诊疗。

      利益冲突 所有作者声明无利益冲突

      作者贡献声明 于荭负责研究命题和方案的建立、论文的撰写和审阅;张胜男负责数据的收集和分析、论文的撰写、文献的查阅;段丽莎、孔蒙娟负责数据的采集和处理、论文的修订;崔建岭负责研究的设计和支持、论文的审阅

参考文献 (23)

目录

    /

    返回文章
    返回