basedongraylevelco-occurrencematrix[J].JClinRadiol,2019,38(3):386-389.
[13] 谭红娜,武明辉,顾建钦,等.乳腺X线图像纹理特征预测乳asaspossible
腺癌腋窝淋巴结转移[J].中国医学影像技术,2017,33(12):1774-1778.
TANHN,WUMH,GUJQ,etal.Mammogramtextureanalysisinpredictionofaxillarylymphnodemetastasisforbreastcarcinoma[J].ChinJMedImagingTechnol,2017,33(12):1774-1778.
[14] HOLM,SAMEIE,MAZUROWSKIMA,etal.Cantextureanaly
美国2016大选sisbeusedtodistinguishbenignfrommalignantadrenalnodulesonunenhancedCT,contrast-enhancedCT,orin-phaseandop posed-phaseMRI[J].AJRAmJRoentgenol,2019,212(3):
崇文门新世界电话554-561.
[15] SONGG,XUEF,ZHANGC.Amodelusingtexturefeaturesto
differentiatethenatureofthyroidnodulesonsonography[J].JUl trasoundMed,2015,34(10):1753-1760.
[16] ZHANGGM,SUNH,SHIB,etal.QuantitativeCTtextureanaly
sisforevaluatinghistologicgradeofurothelialcarcinoma[J].Ab domRadiol
(NY),2016,42(2):561-568.[17] GANESHANB,GOHV,MANDEVILLEH,etal.Non-smallcell
lungcancer:histopathologiccorrelatesfortextureparametersatCT[J].Radiology,2013,266(1):32
6-336.
(编校:蔺癑)
临床影像组学列线图鉴别前列腺良恶性结节的应用价值
邢朋毅1,彭雯佳1,付 贝1,段绍峰2,阳青松1,陆建平1
Applicationvalueofclinicalradiomicsnomogramindifferentialdiagnosisofbenignandmalignantprostaticnodules
XINGPengyi1,PENGWenjia1,FUBei1,DUANShaofeng2,YANGQingsong1,LUJianping
英语单词学习11
DepartmentofRadiology,theFirstAffiliatedHospital,NavalMedicalUniversity,Shanghai200433,China;2GEHealthcare,Shanghai
210000,China.
【Abstract】 Objective:ToestablishaclinicalradiomicsnomogrambasedondualparameterMRandevaluateitsapplicationvalueindifferentiatingbenignandmalignantprostaticnodules.Methods:FromApril2014toJuly2017,227patients(99benignprostatichyperplasiaand128prostatecancer)withprostateMRIandconfirmedbypathologywereselected.RadiomicsfeatureswereextractedfromaxialT2WIandADCimagesofeachpatient.Theminimumab solutecontractionandLASSOregressionwereusedtoreducethedimension,andtheeffectivefeatureswereselectedandtheradiomicsfeaturemodelwasconstructed.Combinedwithindependentclinicalriskfactors,multivariateLogisticregressionmodelwasusedtoestablishradiomicsnomogram.Thecalibrationandd
iagnosisefficiencyofnomogramwasevaluatedinthetrainingsetandverifiedintheverificationset.Finally,theclinicalapplicationvalueofnomogramwasevaluatedbydecisioncurveanalysis.Results:Theageofthepatient,PSA,theareaofthelesionandtheselectedima gingfeaturesweresignificantlycorrelatedwiththedifferentialdiagnosisofbenignandmalignantprostaticnodules.Theradiomicsmodelhadgooddiscriminationefficiencyforbenignandmalignantprostaticnodules.Thetrainingset(AUC=0.91,95%CI:0.87~0.96)andvalidationset(AUC=0.90,95%CI:0.81~1.00)wereeffective.Theradiomicsnomogram,combinedwiththeradiomicscharacteristicsofdualparameterMRimagesandindependentclinicalriskfactors,showedbetterdifferentiationefficiencyandpathologicalcoincidencerateintrainingset(AUC=0.95,95%CI:0.91~0.98)and
validationset(AUC=0.92,95%CI:0.83~1.00).Thedecisioncurveshowedtheclinicalap plicationvalueoftheclinicalradiomicscombinationmodel.Conclusion:TheradiomicsnomogramofdualparameterMRIcombinedwithclinicalriskfactorscaneasilyidentifybenignandmalignantprostatenodules,andprovidesugges tionsforfurtherclinicaldecision-making.
【收稿日期】 2020-08-29【修回日期】 2020-09-23
【基金项目】 国家临床重点专科军队建设项目(原总后卫生部)【作者单位】
1海军军医大学第一附属医院影像医学科,上海 2004332
progress可数吗
(GE)通用医疗,上海 210000
【作者简介】 邢朋毅(1979-),男,河南洛阳人,副主任医师,博士,主要从事腹部影像学研究。E-mail:drxingpengyi@163.com
【通讯作者】 陆建平(1959-),男,江苏太仓人,主任医师,教授,博士生导师,主要从事腹部影像学研究。E-mail:cjr.lujianping@vip.163.com
·
9041·现代肿瘤医学 2021年04月 第29卷第08期 MODERNONCOLOGY,Apr 2021,VOL 29,No 08
【Keywords】prostatecancer,magneticresonanceimaging,radiomics,machinelearning
ModernOncology2021,29(08):1409-1414
【摘要】 目的:建立基于双参数MR的临床影像组学列线图,评价其鉴别前列腺良恶性结节的应用价值。方法:选择2014年4月至2017年7月间接受前列腺MRI检查并经病理证实的前列腺结节患者227例(前列腺结节增生99例,前列腺癌128例)。从每个患者的轴位T2WI和ADC图像中提取影像组学特征。使用最小绝对收缩和LASSO回归进行降维,选择有效特征并构建影像组学特征模型。结合独立的临床危险因素,采用多元Logistic回归模型建立影像组学列线图。列线图的校准和诊断效能在训练集中进行评估,并在验证集中进行验证。最后通过决策曲线分
析评价列线图的临床应用价值。结果:患者年龄、PSA、病灶所在区域以及筛选出的影像组学特征与前列腺结节的良恶性鉴别显著相关。影像组学模型对前列腺良恶性结节有较好的鉴别效能,训练集(AUC=0.91,95%CI:0.87~0.96),验证集(AUC=0.90,95%CI:0.81~1.00)。影像组学列线图结合了双参数MR图像的影像组学特征和独立的临床危险因素,在训练集(AUC=0.95,95%CI:0.91~0.98)和验证集(AUC=0.92,95%CI:0.83~1.00)中表现出更好的鉴别效能和病理符合率,决策曲线展示了临床影像组学结合模型的临床应用价值。结论:结合双参数MRI影像组学特征和临床危险因素的影像组学列线图可以简便地鉴别出前列腺结节的良恶性,为下一步临床决策提供建议。
【关键词】前列腺癌;磁共振成像;影像组学;机器学习
【中图分类号】R737.25 【文献标识码】A DOI:10.3969/j.issn.1672-4992.2021.08.028
【文章编号】1672-4992-(2021)08-1409-06
前列腺癌(PCa)是男性常见的恶性肿瘤之一,我国前列腺癌的发生率亦逐年上升[1]。前列腺癌的过度诊断和过度治疗是现代前列腺癌管理中备受关注的问题[2]。伴随着多参数磁共振(mpMRI)的推广应用,大量的高级别前列腺癌被检出,降低了不必要的前列腺活检手术的发生率及
其并发症,当前影像组学在肿瘤异质性的研究受到越来越多的关注[3],影像组学利用计算机技术对医学图像中的大量数据进行定量分析,其结果有助于临床决策,本研究旨在建立并验证基于双参数MRI的影像组学特征和临床危险因素相结合的临床影像组学结合模型,用于鉴别前列腺结节的良恶性。1 资料与方法
1.1 临床资料
回顾性纳入海军军医大学第一附属医院2014年4月至2017年7月间接受3TMR检查的前列腺结节病人,其中前列腺结节增生99例,前列腺癌128例,年龄为37~90岁,平均年龄为(66.70±8.74)岁,血清前列腺特异抗原(prostatespe cificantigen,PSA)为4.11~234.00ng/mL,平均为(29.19±21.4)ng/mL。所有病人皆行MR引导下前列腺穿刺活检或前列腺全切根治术,并获得病理及完整临床资料。临床表现为有排尿困难、尿频、尿急等症状或PSA升高(>4ng/mL)。纳入标准:①前列腺穿刺活检或全切术后,病理确诊者;②术前行3TMR检查且检查设备相同,序列一致;③MR检查后,两周内获得病理结果。排除标准:①行MR检查前,前列腺曾进行过活检、手术、放疗或内分泌治疗;②结节体积过小(最大直径<3mm),病灶边界难以勾画;③MR图像存在伪影,并影响病灶分割。其他临床指标有前列腺病灶体积(prostatelesionvolume,PLV)为0.35~90.30mm3,平均为(7.37±14.07)mm3,病灶位置(外周带70例
、移行带126例、弥漫分布31例)。肿瘤病例的Gleason评分:3+3分36例、3+4分24例、4+3分18例、4+4分23例、4+5分15例、5+4分7例、5+5分5例。本研究经海军军医大学第一附属医院伦理委员会审批。
1.2 设备与方法
采用西门子3T超导MR系统(MAGNETOMSkyra,Sie mensHealthcare,Erlangen,Germany),采用标准18通道相控阵体线圈和32通道集成脊柱线圈采集信号。横断面T2WI参数为重复时间/回波时间(TR/TE)=5460/104ms,视野(FOV)为180mm2×180mm2,矩阵384×384,层数24,层厚4mm,层间距为0,回波链长度18,采集时间为3min49s;横断面DWI的主要参数为TR/TE=5100/89ms,FOV为224mm2×280mm2,矩阵120×150,层数20,层厚4mm,层间距为0,加速因子2,b值=0,1500s/mm2,采集三个正交方向的信号,采集时间为7min59s。
1.3 影像组学分析
病灶分割和特征提取:预处理后的T2WI和ADC图像中成像层面应保持匹配,MRI上所描绘的感兴趣区(ROI)要与病理描述的位置相一致,外周带(PZ)病变勾画以ADC图为主,移行带
(TZ)病变勾画以T2WI为主。选择病灶最大截面,沿着病灶边界勾画感兴趣区,如图1所示,以获得二维数据,对于多灶前列腺癌病灶,选择Gleason评分(GS)最高的病变勾画感兴趣区,如果GS相同,则选择直径最大的病变勾画病灶,对于活检病理证实的非癌前列腺区域的患者,勾画感兴趣区的直径不小于3mm。采用组内和组间相关系数(ICC)评价影像组学特征提取的重复性。两位在前列腺癌诊断方面有6年(医生A)和10年(医生B)经验的影像科医生在T2WI和ADC图像上手动将病灶感兴趣区分割出来,病理结果对他们盲化。首先,两位医生随机分析了30张图片,以评估组间的重复性。一周后,医生A重复了同样的程序,ICC大于0.75表示特征提取的一致性很好,其余的图像分割由医生A执行。图像预处理及特征计算使用AK软件(Artifi cialIntelligenceKit;version3.2.0,GEHealthcare)。
·
0
1
4
新东方英语培训学校1
·邢朋毅,等 临床影像组学列线图鉴别前列腺良恶性结节的应用价值
图1 病灶分割
患者,72岁,移行带前列腺癌,PSA为24.33ng/mL,Gleason评分为3+4分。以ADC图和DWI为参照,在T2WI上勾画肿瘤感兴趣区。A-B:轴位T2加权成像(T2WI);C:扩散加权成像(DWI),b值=1500s/mm2;D:表观扩散系数图(ADC)。
Fig.1 Lesionsegmentation
Thepatient,72yearsold,hadtransitionalzoneprostatecancerwithaPSAof24.33ng/mL,Gleasonscoreof3+4.WithADCmapandDWIasreference,theregionofinterestoftumorwasdelineatedonT2WI.A-B:AxialT2-weightedimaging(T2WI).C:Diffusion-weightedimaging(DWI),bvalue=1500s/mm2.D:Apparentdiffusioncoefficient(ADC).
1.3.1 影像组学特征选择与标签构建 采用适用于高维数据回归的最小绝对收缩选择算子(LASSO)Logistic回归算法[4],从训练集中选取与前列腺良恶性结节鉴别最相关的非零系数特征,并对所选特征各自的系数加权进行线性组合,为每个患者计算影像组学评分(Radscore),以量化影像组学特征。在训练集和验证集中,影像组学特征的预测准确性通过受试者操作特征曲线下的面积(AUC)来量化。1.3.2 影像组学列线图的构建 多元Logistic回归分析用于筛选与前列腺结节良恶性鉴别的独立预测因素,其中包括潜在的预测因素,如影像组学特征和临床危险因素。为了给临床医生提供个体化鉴别良恶性结节的定量工具,我们在多元Logistic回归模型的基础上构建了兼具影像组学和临床特征的影像组学列线图。
1.3.3 影像组学列线图的性能和验证 量化影像组学列线图的鉴别性能,使用从训练数据导出的公式计算验证数据中每个患者的Radscore,同时除计算AUC外,亦进行了校准和Hosmer-Lemeshow检验。通过量化不同阈值概率下联合模型训练集和验证集的净效益,进行决策曲线分析(DCA)。1.4 统计学分析
采用R软件(3.5.1版,www.Rproject.org)进行统计分析。采用连续变量的t检验和分类变量的Fisher精确检验来检测非PCa组和PCa组的临床特征差异。经统计学处理,差异有统计学意义(P<0.05)。计算组内相关系数(ICC)以确定观察者间和观察者内的一致性,并保留ICC值>0.75的特征。LASSO回归是使用“glmnet”包执行的,
ROC分析用“PROC”进行,校准曲线和Hosmer-Lemeshow检验分别用“ModelGood”和“DescTools”进行。最后,使用“rmda”软件包进行决策曲线分析,以确定影像组学特征及列线图的临床实用性。
2 结果
2.1 临床模型
在临床因素方面,单因素Logistic分析显示,年龄、PSA、病灶所在区域是预测有临床意义的PCa的显著因素(表1)。多因素Logistic回归分析显示,年龄、PSA和病灶所在区域有统计学差异(P<0.05),这3个临床因素可以作为独立的预测因素(表1)。最后,根据选定的临床特征建立Logistic回归分类器。临床模型鉴别效能见表2。
表1 临床因素单变量和多变量Logistic分析结果Tab.1 UnivariateandmultivariateLogisticanalysisresults
Clinical
factors
Univariate
OR(95%CI)P
Multivariate
OR(95%CI)P
Age1.07(1.03~1.12)0.0011.06(1.01~1.12)0.040
Zone0.18(0.15~0.49)<0.0010.13(0.06~0.28)<0.001PSA1.06(1.03~1.11)0.0031.10(1.05~1.17)<0.001PLV1.01(0.99~1.05)0.340--注:OR:优势比;PSA:前列腺特异抗原;Zone:前列腺病灶所在区域;PLV:前列腺病灶体积。
Note:OR:Oddsratio.PSA:Prostatespecificantigen.Zone:Zoneofpros tatelesion.PLV:Prostatelesionvolume.
2.2 影像组学模型
共从感兴趣区中提取了396个影像性特征。在观察者内和观察者间一致性检验中,所有特征的IC
C均大于0.75。使用LASSO算法和10倍交叉验证将所有特征精简为5个潜在的预测因子,这些预测因子基于训练队列中的159名患者,运用这些预测因子来建立LASSO回归模型(图2)。模型中使用的特征及影像学标签如图3所示。
这些选中的特征被用来计算影像组学评分,Radscore是由构建模型的特征乘以各自对应的系数之和所得到,公式为“Radscore=βi×Xi+常数”。加权特征如下所示:Radscore=1.501×ADC_log_sigma_3_0_mm_3D_glrlm_RunLengthNon Uniformity+-0.146×T2_log_sigma_3_0_mm_3D_glszm_Size ZoneNonUniformityNormalized+0.226×T2_log_sigma_3_0_mm_3D_firstorder_Maximum+0.054×T2_log_sigma_3_0_mm_3D_glrlm_GrayLevelVariance+-0.751×ADC_wavelet_LLL_fir storder_10Percentile+0.682。
影像组学模型对前列腺结节良恶性鉴别的预测效率较高,见表2。
2.3 临床影像组学联合模型
结合年龄、PSA和病灶所在区域和影像组学评分建立临床影像组学列线图如图4所示。临床-影像
组学联合模型显示了对前列腺良恶性结节良好的鉴别能力,训练组和验证组的AUC值分别为0.95(95%CI:0.91~0.98)和0.92(95%CI:0.83~1.00),较高的病理符合率在校准曲线(图5)中体现,表2中列出了三种模型的AUC、准确性、敏感性和特异性。我们绘制了这三种模型ROC曲线及决策曲线进行比较
·
1
1
4
1
·
现代肿瘤医学 2021年04月 第29卷第08期 MODERNONCOLOGY,Apr 2021,VOL 29,No 08
(图6、7)。结果表明:临床影像组学模型对前列腺良恶性结
节的鉴别能力高于单纯的临床模型及影像组学模型。
表2 临床、影像组学及联合模型鉴别前列腺良恶性结节的AUC结果
Tab.2 TheAUCresultsofclinical,radiomicsandcombinedmodelstoidentifybenignandmalignantnodules
IndexClinicalmodel
Trainingset
Validationset
Radiomicsmodel
Trainingset
Validationset
商务英语培训班
Combinedmodel
Trainingset
Validationset
Cut-off0.02
英语在线-0.21
-0.47
AUC(95%CI)0.84(0.78~0.90)
0.87(0.78~0.96)
0.91(0.87~0.96)
0.90(0.81~1.00)
0.95(0.91~0.98)
0.92(0.83~1.00)
Accuracy0.800.810.860.850.890.90Sensitivity0.790.820.900.950.880.88Specificity0.810.800.810.730.910.93PPV0.840.840.860.820.930.95NPV
0.750.770.860.920.850.83P(vscombinedmodel)
<0.001
<0.001
<0.001
<0.001
-
opus
-
注:95%CI:95%置信区间;PPV:阳性预测值;NPV:阴性预测值。
Note:95%CI:95%confidenceinterval.PPV:Positivepredictvalue.NPV:Negativepredictivevalue
.
图2 LASSO回归模型的建立
A:双参数MR影像组学模型二项式偏差随参数λ变化曲线图。纵轴为二项式偏差,横轴表示log(λ)值,上方数字表示筛选出特征的数量,模型的二项式偏差最小处的λ为最优值(垂直虚线)。B:双参数MR模型影像组学特征系数随λ变化的曲线图。上方数字表示筛选出特征的数量。
Fig.2 TheestablishmentofLASSOregressionmodel
A:CurveofbinomialdeviationofbiparameterMRradiomicsmodelvaryingwithparameterλ.Theverticalaxisisbinomialdeviation.Thehorizontalaxisre presentsthelog(λ)value.Thenumberaboverepresentsthenumberofselectedfeatures,andtheλattheminimumbinomialdeviationofthemodelistheoptimalvalue(thecurveoftheimagegroupcharacteristiccoefficientoftheverticaldottedline).B:BiparameterMRmodelchangingwithλ.Thenumberaboveindicatesthenumberoffeaturesfilteredout
.
图3 组学模型中使用的特征以及影像组学标签
A:双参数MR模型筛选出的影像组学特征;B:双参数MR模型训练集(左)与验证集(右)影像组学评分的比较,蓝色标签为前列腺结节增生,黄色标签为前列腺癌。
Fig.3 Featuresandradiomicslabelsusedinthegroupmodel
A:Imagingfeaturesscreenedbytwo-parameterMRmodel.B:Comparisonofimagingscorebetweentwo-parameterMRmodeltrainingset(left)andvalidationset(right).Thebluelabelisprostatenodularhyperplasiaandtheyellowlabelisprostatecancer.
·2141·邢朋毅,等 临床影像组学列线图鉴别前列腺良恶性结节的应用价值
图4 临床影像组学列线图
Zone中0表示外周带,1表示移行带,2表示弥漫分布。
Fig.4 Clinicalradiomicsnomogram
Inzone,0indicatestheperipheralzone,1indicatesthetransitionzone,and2indicatesthediffusedistribution.
3 讨论
影像组学近年来在肿瘤研究中得到了广泛的应用,其优势在于它可以自动从图像中过滤获得全面的数据[5],在宏观图像和分子水平的基因或蛋白质变化中识别肿瘤的异质性[3],目前已被广泛用于对头颈部肿瘤、肺癌和结直肠癌等病变的检测、分类、治疗反应评估及预测[6],影像组学具有无创、快速、可重复性强和低成本的特点,以磁共振成像为基础的影像组学作为一个发展趋势和前景广阔的研究领域,已被广泛应用于肿瘤异质性的评价,并取得了较高的诊断效率[7]。
在本研究中,我们建立并验证了一个基于MR双参数影像组学特征的列线图来鉴别前列腺结节的良恶性。此临床影像组学列线图提供了一种易于使用、可量化和个体化的筛查前列腺癌的工具,有助于避免对前列腺癌患者进行不必要的治疗和侵入性检查,并防止及推迟了低级别前列腺癌的进展。影像组学特征在前列腺非癌组织与癌组织之间以及移行带、外周带肿瘤之间均有显著差异存在[8]。WIBMER等[9]对147例患者的ADC、T2WI图像进行全病灶纹理分析,结果表明Haralick纹理特征可区分前列腺组织的癌性和非癌性,ADC
图上的肿瘤能量和熵与格里森得分相关。
图5 临床影像组学列线图预测前列腺结节良恶性的结果与病理结果一致性的校准曲线(左侧为训练集,右侧为验证集)
Fig.5 Thecalibrationcurveofclinicalradiomicsnomogramforpredictingtheconsistencyofbenignandmalignantprostatenoduleswithpathologicalre sults(thetrainingsetisontheleft,andthevalidationsetisontheright
)
图6 三种模型预测前列腺良恶性结节的受试者操作特征(ROC)曲线(左侧为训练集,右侧为验证集)
Fig.6 Predictionofreceiveroperatingcharacteristic(ROC)curveofbenignandmalignantprostaticnodulesbythreemodels(thetrainingsetisontheleft,andthevalidationsetisontheright
重庆心理咨询师培训)
图7 三种模型的临床决策曲线
Fig.7 Clinicaldecisioncurveofthethreemodels LITJENS等[10]对70例前列腺全切患者术前mpMRI的影像研究发现,高b值DWI图像在鉴别良性前列腺增生(BPH)和前列腺癌方面优于ADC,而ADC最适合鉴别前列腺上皮内瘤变(PIN)和前列腺癌。用T2加权和DWI训练的U-Net实现了与临床前列腺成像报告和数据系统评估类似的性能[11]。影像组学特征的诊断效率高于临床独立危险因素,表明影像组学方法比传统的临床特征更能反映前列腺恶性结节存在的风险。因此,影像组学方法可以通过我们现有的MRI图像结合简易的临床指标无创地量化肿瘤的异质性。
·
3
1
4
1
·
现代肿瘤医学 2021年04月 第29卷第08期 MODERNONCOLOGY,Apr 2021,VOL 29,No 08