劉思野，文 露，侯 靜，聶少麟，周菊梅，曹 芳，盧 強，覃玉卉，于小平

1.中南大學湘雅醫(yī)學院附屬腫瘤醫(yī)院，湖南省腫瘤醫(yī)院放射診斷科，湖南 長沙 410006；

2.中南大學湘雅醫(yī)學院附屬腫瘤醫(yī)院，湖南省腫瘤醫(yī)院結(jié)直腸外科，湖南 長沙 410006；

3.中南大學湘雅醫(yī)學院附屬腫瘤醫(yī)院，湖南省腫瘤醫(yī)院放療科，湖南 長沙 410006；

4.中南大學湘雅醫(yī)學院附屬腫瘤醫(yī)院，湖南省腫瘤醫(yī)院病理科，湖南 長沙 410006

局部進展期直腸癌的標準治療方案是新輔助放化療聯(lián)合全直腸系膜切除術［1-2］，而病理學檢查是評價新輔助放化療效果的金標準。早期評價直腸癌新輔助放化療的病理學反應，可幫助患者實現(xiàn)個性化治療。目前，有不少研究［3-5］采用動態(tài)對比增強磁共振成像（dynamic contrast-enhanced magnetic resonance imaging，DCE-MRI）預測直腸癌患者對新輔助放化療的病理學反應狀態(tài)，但結(jié)論不一，甚至相互矛盾。原因可能在于DCE-MRI參數(shù)值不能反映腫瘤的異質(zhì)性［6］。紋理分析，作為一種新興的醫(yī)學圖像分析方法，可以量化腫瘤異質(zhì)性［7-9］。目前，鮮有研究報道基線DCE-MRI參數(shù)圖像的紋理特征是否有助于預測直腸癌原發(fā)灶新輔助放化療后達到病理學完全緩解（pathological complete response，pCR）。因此，本研究旨在探討基線DCE-MRI紋理特征在評估直腸癌原發(fā)灶放化療病理學反應方面的潛能。

1 資料和方法

1.1 一般資料

本研究經(jīng)中南大學湘雅醫(yī)學院附屬腫瘤醫(yī)院醫(yī)學倫理委員會批準。共收集38例局部進展期（cT3-4N0-2M0）直腸癌患者，包括27例男性患者和11例女性患者，年齡22～73歲，平均年齡51歲。所有患者均在新輔助放化療前經(jīng)結(jié)腸鏡下病理學活檢證實為直腸非黏液型腺癌。

1.2 設備及掃描方法

所有患者于新輔助放化療1～3 d前，依次行盆腔常規(guī)MRI平掃及DCE-MRI掃描。患者取仰臥位，使用美國GE公司的MR360 1.5 T超導型磁共振儀及相控體線圈。常規(guī)平掃檢查：① 軸位T2WI快速自旋回波（fast spin echo，F(xiàn)SE）圖像，重復時間（repetition time，TR）為4 435 ms，回波時間（echo time，TE）為102 ms，層厚5 mm，層距1 mm，視野（ field of view，F(xiàn)OV）380 mm，采集矩陣320×224，激發(fā)次數(shù)（number of excitation，NEX）為2；② 斜軸位高分辨率T2加權(quán)FSE圖像，TR為4 500 ms，TE為102 ms，層厚3 mm，層距0.5 mm，F(xiàn)OV為256 mm，采集矩陣200×200，NEX為4；③ 矢狀位高分辨率T2加權(quán)FSE圖像，TR 4 500 ms，TE 102 ms，層厚3 mm，層距0.5 mm，F(xiàn)OV為256 mm，采集矩陣200×200，NEX為4。DCE-MRI檢查：① 肝臟快速三維容積動態(tài)增強（liver acquisition with volume acceleration，LAVA）-T1WI平掃：翻轉(zhuǎn)角為2o和15o，TR為3 ms，TE為1.3 ms，F(xiàn)OV為380 mm，層厚3 mm，層距0.5 mm；② 多時相LAVA-T1WI動態(tài)增強掃描：翻轉(zhuǎn)角15o，TR為3 ms，TE為1.3 ms，F(xiàn)OV為380 mm，層厚3 mm，層距0.5 mm。每個時相6 s，共掃描56個時相。在多時相動態(tài)增強掃描開始后18 s使用高壓注射器經(jīng)肘部靜脈注射造影劑釓雙胺（劑量：0.1 mmol/kg，流速3.5 mL/s），隨后以20 mL 0.9%的NaCl溶液沖管。

1.3 圖像采集及紋理分析

將原始磁共振數(shù)據(jù)導入至GE Advantage Windows Workstation 4.6（AW 4.6）工作站，使用Cinetool軟件進行后處理。由1名放射診斷科醫(yī)師在LAVA-T1WI動態(tài)增強掃描強化最明顯的一期圖像上參照軸位T2WI圖像逐層沿直腸癌原發(fā)灶外緣手動勾畫感興趣區(qū)域（region of interest，ROI），并自動生成4種DCE-MRI參數(shù)圖，分別為容量轉(zhuǎn)運常數(shù)（volume transfer constant，Ktrans）、速率常數(shù)（rate constant，Kep）、血漿容量分數(shù)（plasma volume fraction，fPV）、血管外細胞外間隙容量分數(shù)（extravascular extracellular volume fraction，Ve）。病灶所有層面的ROI勾畫完成后，軟件自動生成感興趣體積（volume of interest，VOI），并提供相應的DCE-MRI參數(shù)值。

DCE-MRI參數(shù)圖的紋理分析采用MaZda軟件。MaZda可從每個VOI中獲?。孩?9種一級紋理特征，包括均值、方差、偏度、峰度，以及第1、10、50、90和99百分位數(shù)（Pere.1%、Pere.10%、Pere.50%、Pere.90%、Pere.99%）；② 11種二級紋理特征（即共生矩陣特征），包括角二階矩（angular second moment，AngScMom）、對比度（Contrast）、相關度（Correlat）、差分熵（difference entropy，DifEntrp）、差分方差（difference variance，DifVarnc）、熵（Entropy）、逆差矩（inverse difference moment，InvDfMom）、和平均數(shù)（sum average，SumAverg）、和熵（sum entropy，SumEntrp）、平方和（sum of squares，SumOfSqs）、和方差（sum variance，SumVarnc）。

1.4 病理學反應評估

患者行MRI檢查后，接受新輔助放化療，放療總劑量為50 Gy，分25次完成，同期口服卡培他濱1 650 mg/（m2·d）。放化療結(jié)束后8周行全直腸系膜切除術。腫瘤原發(fā)灶的病理學反應狀態(tài)由1名病理科醫(yī)師根據(jù)Dworak腫瘤消退分級（tumor regression grade，TRG）標準進行評估［10］：① TRG 4級，鏡下無存活腫瘤細胞；② TRG 3級，纖維組織中有散在的存活腫瘤細胞；③ TRG 2級，存活腫瘤細胞數(shù)量較多，但少于纖維組織；④ TRG 1級，存活腫瘤細胞多于纖維組織；⑤ TRG 0級，腫瘤細胞對治療無反應。pCR組包括病理學反應為TRG 4級的患者，非pCR組包括病理學反應屬TRG 0～3級的患者。

1.5 統(tǒng)計學處理

采用SPSS 22.0及Medcalc 15.0進行統(tǒng)計學分析。計量資料以±s表示，計數(shù)資料以頻數(shù)表示。P＜0.05為差異有統(tǒng)計學意義。pCR組與非pCR組治療前DCE-MRI參數(shù)值及紋理特征值上的差異（單變量分析）用Mann-Whitney U檢驗比較。pCR的獨立預測因子由多變量logistic回歸分析獲取，共構(gòu)建3個logistic回歸模型，分別為模型1（基于一級紋理特征）、模型2（基于二級紋理特征）及模型3（基于一級及二級紋理特征）。單個紋理特征和3個多變量回歸模型對pCR的預測能力用受試者工作特征（receiver operating characteristic，ROC）曲線評價。采用Fisher確切概率法比較pCR與非pCR兩組之間在T分期、N分期和病理學分化程度上的差異，比較一級與二級紋理特征組之間在紋理特征構(gòu)成比上的差異。紋理特征構(gòu)成比定義為對pCR有預測能力的紋理特征數(shù)占該組紋理特征總數(shù)的百分比。

2 結(jié) 果

腫瘤T分期：cT3期患者有10例，cT4期28例。N分期：cN0、cN1及cN2期分別有11、20和7例。病理分化程度：高、中和低分化患者分別為6、25和7例。T分期、N分期及病理學分化程度在pCR組（8例）與非pCR組（30例）之間差異均無統(tǒng)計學意義（P值分別為0.653、1.000和0.848）。

Ktrans、Kep、fPV及Ve參數(shù)值在pCR組與非pCR組之間差異無統(tǒng)計學意義（表1）。在36個一級紋理特征中，有2個（5.56%）對pCR有預測價值，即均數(shù)Ve（MeanVe）及偏度Ve（SkewnessVe）；在44個二級紋理特征中，有16個（36.36%）對pCR有預測價值（表2）。二級紋理特征構(gòu)成比明顯高于一級紋理特征的構(gòu)成比（P=0.001）。

與非pCR組相比，pCR組的SkewnessVe、相關度Kep（CorrelatKep）、差分方差fPV（difference variance of fPV，DifVarncfPV）、差分方差Kep（diあerence variance of Kep，DifVarncKep）、差分熵fPV（diあerence entropy of fPV，DifEntrpfPV）、差分熵Ve（diあerence entropy of Ve，DifEntrpVe）、熵fPV（EntropyfPV）、熵Ve（EntropyVe）、和平均數(shù)fPV（sum average of fPV，SumAvergfPV）、和熵fPV（sum entropy of fPV，SumEntrpfPV）及和熵Ve（sum entropy of Ve，SumEntrpVe）值更高，而MeanVe、角二階矩fPV（angular second moment of fPV，AngScMomfPV）、角二階矩Ve（angular second moment of Ve，AngScMomVe）、對比度ktrans（Contrastktrans）、逆差矩fPV（inverse diあerence moment of fPV，InvDfMomfPV）、和平均數(shù)ktrans（sum average of Ktrans，SumAvergKtrans）及和平均數(shù)Ve（sum average of Ve，SumAvergVe）值更低（P均＜0.05，表2）。圖1和圖2分別為1例pCR及1例非pCR患者的MRI圖像。

就紋理特征對pCR的預測能力而言，單因素分析時，單個紋理特征的ROC曲線的曲線下面積（area under curve，AUC）為0.571～0.817，靈敏度為33.33%～100.00%，特異度25.00%～100.00%（表2）。在多變量分析中，模型1中pCR的獨立預測因子為偏度ve（P=0.014），此模型的AUC為0.729（95% CI：0.514～0.944），靈敏度90.00%，特異度62.50%；模型2中pCR的獨立預測因子是SumAvergfPV（P=0.003），此模型的AUC值為0.975（95%CI：0.930～1.000），靈敏度93.33%，特異度100.00%；模型1與模型2的AUC值之間差異有統(tǒng)計學意義（P=0.028）；模型3中pCR的獨立預測因子、AUC值、靈敏度及特異度與模型2完全一致。

表1 pCR組與非pCR組之間在DCE-MRI參數(shù)值上的差異

表2 DCE-MRI參數(shù)圖紋理特征的單變量分析結(jié)果

圖1 新輔助放化療后局部進展期直腸癌（cT3N0M0期）患者1影像圖及病理學觀察

圖2 新輔助放化療后局部進展期直腸癌（cT3N0M0期）患者2影像圖及病理學觀察

3 討 論

3.1 一級紋理特征與DCE-MRI參數(shù)值的比較

本研究顯示，在pCR與非pCR兩組之間在直腸癌原發(fā)灶的4個治療前DCE-MRI參數(shù)值上差異無統(tǒng)計學意義，提示基線DCEMRI參數(shù)值對直腸癌新輔助放化療病理學反應狀態(tài)的預測價值有限，與Gollub等［4］和Kim等［11］的研究結(jié)果相似。本研究的一級紋理特征中，MeanVe及SkewnessVe在pCR與非pCR組之間差異有統(tǒng)計學意義，提示一級紋理特征對直腸癌pCR有預測價值。相關研究［12］發(fā)現(xiàn)，SkewnessVe可以幫助鑒別直腸癌放化療后病理學反應良好組（TRG 3～4級）和反應不良組（TRG 0～2級）。以上研究表明，一級紋理特征在直腸癌療效預測方面可能優(yōu)于DCE-MRI參數(shù)值，原因可能在于DCE-MRI參數(shù)值是基于平均值分析，不能反映瘤內(nèi)異質(zhì)性［6］，而腫瘤的異質(zhì)性會影響腫瘤對治療的反應［13］。

3.2 二級紋理特征與一級紋理特征的比較

本研究發(fā)現(xiàn)，在預測pCR方面，二級紋理特征構(gòu)成比明顯高于一級紋理特征構(gòu)成比，提示二級紋理特征可以為預測直腸癌pCR提供更多的指標。本研究也發(fā)現(xiàn)，與基于一級紋理特征構(gòu)建的多變量模型相比，基于二級紋理特征構(gòu)建的多變量模型對pCR的預測效能更高（AUC分別為0.975、0.729，P=0.028）。本研究還發(fā)現(xiàn)，基于一級特征與二級特征共同構(gòu)建的多變量模型，其預測直腸癌原發(fā)灶pCR的獨立預測因子為二級紋理特征，并沒有一級紋理特征，且該模型的AUC值、靈敏度及特異度與單獨基于二級紋理特征構(gòu)建的多變量模型完全一致。綜合本研究的上述發(fā)現(xiàn)，在預測直腸癌原發(fā)灶新輔助放化療的病理學反應方面，二級紋理特征可能要優(yōu)于一級紋理特征。其原因可能在于一級紋理特征僅反映腫瘤MRI圖像的像素及體素灰度頻率分布，而二級紋理特征可反映像素及體素之間的空間與位置關系，故二級紋理特征較一級紋理特征更能充分地反映腫瘤的異質(zhì)性［14］。

3.3 紋理特征的含義及其與異質(zhì)性的關系

本研究中，均值、偏度、AngScMom、Correlat、Contrast、Entropy、DifVarnc、DifEntrp、InvDfMom、SumAverg及SumEntrp值在pCR組與非pCR組中差異有統(tǒng)計學意義。均值反映直方圖的中心趨勢，并不能直接反映組織異質(zhì)性，而偏度表示統(tǒng)計數(shù)據(jù)分布偏斜方向和程度的度量，可直接反映組織異質(zhì)性［15］。AngScMom和InvDfMom是組織的同質(zhì)性的度量［16-17］。Contrast用于測量圖像的局部變化［17］，反映圖像的清晰度和紋理溝紋的深淺程度，Entropy反映圖像紋理的復雜度和混亂度［17］，DifVarnc指體素對間灰度差的變化，Contrast、Entropy及DifVarnc可直接反映組織的異質(zhì)性［17-18］。SumEntrp和DifEntrp分別是相鄰體素灰度總和、相鄰體素灰度差的隨機性度量［17］，二者與異質(zhì)性的關系并不明確。Correlat指圖像灰度線性相關度，SumAverg反映圖像的整體亮度［16-17］，它們與組織異質(zhì)性無直接相關［18］。在本研究中，偏度、Contrast、Entropy及DifVarnc在pCR組與非pCR組中差異有統(tǒng)計學意義，提示pCR組與非pCR組在瘤內(nèi)異質(zhì)性方面有一定差別。

本研究存在的不足：① 樣本量小，這會導致一定的統(tǒng)計學偏倚；② 沒有探討直腸癌原發(fā)灶治療后DCE-MRI的紋理特征與其療效之間的關系。

總之，本研究發(fā)現(xiàn)基線DCE-MRI參數(shù)圖的紋理特征在預測直腸癌原發(fā)灶新輔助放化療的病理學反應方面或許有一定價值，且二級紋理特征的預測價值可能優(yōu)于一級紋理特征。