聯(lián)合臨床、MR T2WI及表觀彌散系數(shù)圖影像組學(xué)特征列線圖預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移

姬健智，張 倩，牛 猛，郭順林*，翟亞楠

(1.蘭州大學(xué)第一臨床醫(yī)學(xué)院，甘肅 蘭州 730000；2.蘭州大學(xué)第一醫(yī)院放射科，甘肅 蘭州 730000)

前列腺癌易發(fā)生骨轉(zhuǎn)移。99Tcm亞甲基二磷酸鹽(99Tcm-methylene diphosphonate,99Tcm-MDP)全身骨顯像是檢查骨轉(zhuǎn)移癌的常用方法[1]，敏感度高而特異度低[2]；而常規(guī)骨顯像并非適用于所有初診前列腺癌患者[1]。2014版《中國前列腺癌診斷治療指南》建議對初診前列腺癌患者行99Tcm-MDP骨掃描，但對其是否應(yīng)作為所有初診前列腺癌患者的常規(guī)檢查仍存在爭議[3]。本研究基于臨床聯(lián)合T2WI及表觀彌散系數(shù)(apparent diffusion coefficient, ADC)圖影像組學(xué)特征建立聯(lián)合模型列線圖，觀察其預(yù)測初診前列腺癌骨轉(zhuǎn)移的價值。

1 資料與方法

1.1 一般資料 回顧性收集2017年1月—2022年1月110例于蘭州大學(xué)第一醫(yī)院接受前列腺M(fèi)R檢查、且經(jīng)病理證實的初發(fā)前列腺癌患者，年齡50～84歲、平均(68.7±7.2)歲；根據(jù)99Tcm-MDP骨顯像分為骨轉(zhuǎn)移組(n=50)和無骨轉(zhuǎn)移組(n=60)；MR檢查前均未接受組織穿刺活檢或手術(shù)；99Tcm-MDP骨顯像與MR檢查間隔≤14天[4]，且臨床及影像學(xué)資料完整。排除標(biāo)準(zhǔn)：①前列腺手術(shù)、放射化學(xué)治療、激素治療和靶向治療史；②無法確定骨轉(zhuǎn)移癌來源；③合并其他骨損傷或骨代謝??；④圖像質(zhì)量差；⑤病灶過小或邊界不清。

110例中，腺泡腺癌69例、腺癌41例；T分期結(jié)果為48例T2、24例T3及38例T4；N分期結(jié)果為64例N0、46例N1分期；前列腺癌Gleason評分5～10分，中位評分9分，其中4例Gleason評分≤6分、23例為7分、83例≥8分；總前列腺特異性抗原(total prostate specific antigen, tPSA) 3.36～100.00 ng/ml(參照文獻(xiàn)[5]，將tPSA>100 ng/ml視為100 ng/ml)，中位tPSA為55.26 ng/ml；堿性磷酸酶(alkaline phosphatase, ALP) 39～4 189 U/L，中位ALP為90 U/L，其中78例ALP<121 U/L、32例ALP≥121 U/L；血鈣1.45～2.83 mmol/L，中位血鈣2.22 mmol/L；血磷0.63～1.67 mmol/L，中位血磷1.09 mmol/L。

1.2 儀器與方法 采用Siemens Magnetom Skyra 3.0T MR儀、腹部相控陣線圈采集前列腺M(fèi)RI。囑患者仰臥，頭先進(jìn)，采集軸位T2WI，TR 6 580 ms，TE 99 ms，層厚3 mm，層間距0.8 mm，矩陣256×100，F(xiàn)OV 200 mm×100 mm；行軸位彌散加權(quán)成像(diffusion weighted imaging, DWI)，TR 4 700 ms，TE 88 ms，層厚3 mm，層間距0.8 mm，F(xiàn)OV 200 mm×90 mm，b=0、800、1 400、3 000 s/mm2。

1.3 提取及篩選影像組學(xué)特征 由2名具有10年以上影像學(xué)診斷經(jīng)驗、不知曉骨顯像結(jié)果的副主任醫(yī)師，參考病理所見腫瘤區(qū)域，結(jié)合DWI，以3D Slicer軟件(版本4.10.2，http://www.slicer.org)于軸位T2WI和軟件自動獲得的ADC圖上逐層手動勾畫病灶ROI(圖1)，存在多發(fā)病灶時，根據(jù)病理結(jié)果勾畫Gleason評分最高者；意見發(fā)生分歧時經(jīng)討論達(dá)成一致。

將病灶ROI導(dǎo)入3D Slicer軟件，基于T2WI和ADC圖分別提取851個影像組學(xué)特征，包括14個形態(tài)學(xué)(shape)特征、18個一階直方圖(first order)特征、14個灰度依賴矩陣(gray Level dependence matrix, Gldm)特征、24個灰度共生矩陣特征(gray level cooccurrence matrix, Glcm)特征、16個灰度運(yùn)行長度矩陣(gray level runlength matrix, Glrlm)特征、16個灰度區(qū)域大小矩陣(gray level size zone matrix, Glszm)特征、5個臨域灰度差矩陣(neighbouring gray tone difference matrix, Ngtdm)特征和744個小波(wavelet)濾波特征；將其標(biāo)準(zhǔn)化后，以5折交叉驗證、最小絕對收縮和選擇算子(least absolute shrinkage and selection operator, LASSO)算法篩選最佳影像組學(xué)特征，計算相應(yīng)系數(shù)，并以系數(shù)加權(quán)求和獲得影像組學(xué)評分(rad-score)。

1.4 統(tǒng)計學(xué)分析及構(gòu)建模型 采用R語言(版本4.1.2，https://www.r-project.org/)和SPSS 26統(tǒng)計分析軟件。以中位數(shù)(上下四分位數(shù))表示不符合正態(tài)分布的計量資料。因ALP結(jié)果差異較大，以受試者工作特征(receiver operating characteristic, ROC)曲線判定其診斷骨轉(zhuǎn)移的最佳截斷值，設(shè)定<截斷值定義為0，≥截斷值定義為1，轉(zhuǎn)換為二分類變量。以Wilcoxon秩和檢驗比較無骨轉(zhuǎn)移組與骨轉(zhuǎn)移組間rad-score的差異。以組內(nèi)相關(guān)系數(shù)(intra-class correlation coefficient, ICC)評估2名醫(yī)師勾畫病灶ROI的一致性,ICC>0.75為一致性較好。將最佳影像組學(xué)特征納入多因素logistic回歸分析，構(gòu)建影像組學(xué)模型。以單因素及多因素logistic回歸分析篩選初發(fā)前列腺癌骨轉(zhuǎn)移的臨床相關(guān)獨(dú)立危險因素，建立臨床模型，并以方差膨脹系數(shù)(variance inflation factor, VIF)判斷變量間是否存在共線性。采用多因素logistic回歸構(gòu)建rad-score聯(lián)合臨床獨(dú)立危險因素的聯(lián)合模型，并繪制列線圖將其可視化，以C指數(shù)評估其區(qū)分度，并通過Bootstrap重抽樣1 000次的方法進(jìn)行內(nèi)部驗證，計算校正C指數(shù)，應(yīng)用Hosmer-Lemeshow檢驗評估聯(lián)合模型的擬合優(yōu)度，繪制其校正曲線。繪制各模型預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移的ROC曲線，以DeLong檢驗比較其曲線下面積(area under the cure, AUC)的差異。以決策曲線分析(decision curve analysis, DCA)評價聯(lián)合模型的臨床凈獲益。P<0.05為差異有統(tǒng)計學(xué)意義。

2 結(jié)果

2.1 構(gòu)建模型 2名醫(yī)師勾畫各層面病灶ROI的一致性較好，ICC為0.87[95%CI(0.92,0.95)]～0.93[95%CI(0.90,0.93)]。經(jīng)LASSO分析，最終基于T2WI及ADC圖篩選出11個關(guān)鍵特征(表1)，構(gòu)建影像組學(xué)模型。骨轉(zhuǎn)移組rad-score[0.16 (-0.15，0.41)]高于無骨轉(zhuǎn)移組[-0.49(-0.95, -0.09),Z=-6.159，P<0.001)]。

表1 預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移的最佳影像組學(xué)特征及其系數(shù)

ROC曲線結(jié)果顯示，以ALP 121 U/L為預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移的截斷值，其敏感度為60.00%，特異度為96.70%，約登指數(shù)為0.567。單因素和多因素logistic分析顯示，tPSA、ALP及N分期是骨轉(zhuǎn)移的臨床相關(guān)獨(dú)立危險因素(表2)，三者的VIF值均為1.01，表明變量間不存在多重共線性；以之構(gòu)建臨床模型，并與rad-score構(gòu)建聯(lián)合模型。

表2 單因素和多因素logistic回歸分析初發(fā)前列腺癌骨轉(zhuǎn)移的臨床相關(guān)特征

2.2 各模型效能及繪制列線圖 影像組學(xué)模型預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移的AUC為0.82[95%CI(0.74，0.90)]，敏感度和特異度分別為0.98和0.53；臨床模型AUC為0.93[95%CI(0.88，0.98)]，敏感度和特異度分別為0.82和0.93；聯(lián)合模型AUC達(dá)0.96[95%CI(0.94，0.98)]，高于臨床模型(Z=-2.066，P=0.039)和影像組學(xué)模型(Z=-3.451，P<0.001)，其敏感度和特異度分別為0.84和0.97，見圖2。繪制聯(lián)合模型列線圖(圖3)，以Bootstrap重抽樣1 000次后，校正C指數(shù)為0.96，表明聯(lián)合模型的區(qū)分度良好；Hosmer-Lemeshow檢驗顯示其擬合優(yōu)度佳(χ2=1.13，P=0.99)；校正曲線(圖4)顯示聯(lián)合模型預(yù)測結(jié)果與實際結(jié)果的一致性良好。DCA(圖5)結(jié)果顯示，閾值概率為0～0.98時，聯(lián)合模型的臨床凈獲益大于臨床模型。

3 討論

前列腺癌發(fā)生骨轉(zhuǎn)移后患者死亡率較高[6]，及時診斷和干預(yù)改善預(yù)后尤為重要。目前，影像組學(xué)已廣泛應(yīng)用于醫(yī)學(xué)研究領(lǐng)域，WANG等[7]建立基于動態(tài)增強(qiáng)MRI(dynamic contrast-enhanced MRI, DCE-MRI)聯(lián)合T2WI影像組學(xué)模型，其預(yù)測初發(fā)前列腺癌發(fā)生骨轉(zhuǎn)移的AUC為0.90。ZHANG等[8]聯(lián)合DCE-MRI、T2WI和DWI建立預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移影像組學(xué)模型的AUC為0.86。本研究基于T2WI聯(lián)合ADC圖選出11個最佳影像組學(xué)特征，包括形態(tài)學(xué)、Glrlm及Glszm等特征，以之構(gòu)建影像組學(xué)模型，結(jié)果顯示骨轉(zhuǎn)移組rad-score高于無骨轉(zhuǎn)移組，該模型預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移的AUC為0.82。

本研究針對常見臨床因素行單因素及多因素logistic分析，發(fā)現(xiàn)tPSA、ALP和N分期是初發(fā)前列腺癌骨轉(zhuǎn)移的臨床相關(guān)獨(dú)立危險因素，與陳渤文等[5，8-10]的結(jié)果基本一致。前列腺癌可致血-上皮屏障破壞，PSA直接進(jìn)入血液，使tPSA增高；腫瘤惡性程度越高，組織破壞越嚴(yán)重，骨轉(zhuǎn)移風(fēng)險亦增高。ALP水平變化可評估骨骼變化和成骨細(xì)胞活性及前列腺癌骨轉(zhuǎn)移等[5]。本組ALP預(yù)測骨轉(zhuǎn)移的最佳截斷值為121 U/L，與既往研究[11-12]相似，但高于陳渤文[5]的結(jié)果(84.5 U/L)，可能與樣本不同及樣本量較小有關(guān)。另外，前列腺癌區(qū)域淋巴結(jié)轉(zhuǎn)移提示預(yù)后較差[13]，但對于N分期可否作為預(yù)測骨轉(zhuǎn)移的獨(dú)立危險因素尚存爭議。本研究部分前列腺癌患者僅接受穿刺活檢，其N分期是臨床醫(yī)師結(jié)合病理學(xué)及影像學(xué)所見進(jìn)行評估的結(jié)果，難以避免主觀因素的影響，可能存在偏差。

本研究基于臨床相關(guān)獨(dú)立危險因素建立的臨床模型預(yù)測骨轉(zhuǎn)移的AUC為0.93，高于ZHANG等[8]的臨床模型(0.85)；以臨床結(jié)合rad-score建立的聯(lián)合模型的AUC為0.96，高于文獻(xiàn)報道[7-8]的0.92，且預(yù)測效能優(yōu)于臨床模型和影像組學(xué)模型，其在閾值概率0～0.98時的臨床凈獲益大于臨床模型，并可通過列線圖可視化其中多個變量之間的函數(shù)關(guān)系，應(yīng)用更為便利。聯(lián)合模型列線圖可便捷、有效地預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移，尤其適用于基層醫(yī)療機(jī)構(gòu)。

本研究的主要局限性：①樣本量較小，且缺乏外部驗證；②以骨顯像評估骨轉(zhuǎn)移可能存在假陽性結(jié)果[8]；③未分析其他序列MR圖像及前列腺特異性膜抗原、PSA密度等[14-15]臨床因素。

綜上，聯(lián)合臨床和MR T2WI及ADC圖影像組學(xué)特征的列線圖可有效預(yù)測初發(fā)前列腺癌骨轉(zhuǎn)移。