劉欣，楊海濤，王琪琪，周素穎，趙潔

類風(fēng)濕性關(guān)節(jié)炎(rheumatoid arthritis，RA)和痛風(fēng)性關(guān)節(jié)炎(gouty arthritis，GA)是臨床上常見的滑膜性關(guān)節(jié)炎，二者臨床上均以慢性起病，表現(xiàn)為四肢關(guān)節(jié)腫脹、疼痛及功能障礙，晚期導(dǎo)致關(guān)節(jié)不可逆性破壞、功能喪失而致殘。其診斷均需結(jié)合臨床表現(xiàn)、實(shí)驗(yàn)室和影像學(xué)檢查綜合診斷[1-2]。但實(shí)驗(yàn)室檢查不典型者并不少見，血清陰性RA 占所有RA 患者20%左右[3]；臨床上不伴有高尿酸血癥GA 患者逐漸增多[4]。二者常規(guī)MRI都表現(xiàn)為滑膜增生、關(guān)節(jié)積液、骨髓水腫、骨質(zhì)侵蝕以及關(guān)節(jié)周圍軟組織腫脹[5]，并沒有足夠的特異度。因此，無論在臨床上還是影像學(xué)上特別是疾病早期對(duì)二者的鑒別仍有一定難度。

紋理分析能充分挖掘圖像中肉眼無法觀察到的隱藏信息，能反映組織的潛在病理及生理學(xué)特征，在鑒別診斷方面潛力巨大[6-8]。增強(qiáng)MRI 能區(qū)分增生滑膜與關(guān)節(jié)積液，一直是滑膜性關(guān)節(jié)炎的最佳影像學(xué)檢查[5，9]。但RA和GA治療過程中常需行多次MRI增強(qiáng)檢查進(jìn)行療效評(píng)估，而釓劑的腦沉積、腎功能損害等毒副作用已日益受到關(guān)注，使其應(yīng)用受到限制；文獻(xiàn)報(bào)道[10-11]T2WI序列增生滑膜信號(hào)強(qiáng)度低于關(guān)節(jié)積液，增生滑膜與關(guān)節(jié)積液在T2WI-FS 序列上具有良好對(duì)比度。因此，筆者探討能否利用MRI 平掃的T2WI 脂肪抑制單序列圖像紋理分析來協(xié)助鑒別診斷RA與GA。

1 材料與方法

1.1 一般資料

本研究經(jīng)過我院醫(yī)學(xué)倫理委員會(huì)批準(zhǔn)(批準(zhǔn)文號(hào)：2018-092)，免除受試者知情同意?；仡櫺苑治鑫以?016 年1 月至2020 年9 月臨床綜合確診為RA、GA患者的臨床和MRI 資料。納入標(biāo)準(zhǔn)：RA 納入根據(jù)2010 年ACR/EULAR 標(biāo)準(zhǔn)記分表總分≥6；GA 納入根據(jù)2015年ACR/EULAR分類標(biāo)準(zhǔn)記分表總分≥8。排除標(biāo)準(zhǔn)如下：(1)臨床資料不完善；(2)滑膜不易勾畫ROI(滑膜輕度增生或與關(guān)節(jié)積液、軟骨難以區(qū)分)；(3)臨床診斷不明確；(4)圖像質(zhì)量差不能滿足診斷。最終68 例RA (81 個(gè)關(guān)節(jié))、47 例GA (61 個(gè)關(guān)節(jié))納入本次研究，關(guān)節(jié)部位分布如表1所示。

表1 類風(fēng)濕關(guān)節(jié)炎和痛風(fēng)性關(guān)節(jié)炎關(guān)節(jié)部位分布[個(gè)(%)]Tab.1 Distribution of locations of rheumatoid arthritis and gouty arthritis[n(%)]

1.2 儀器與方法

采用GE Signa Hispeed plus 1.5 T (65 個(gè)關(guān)節(jié))、GE Signa HDx 3.0 T (77個(gè)關(guān)節(jié))磁共振儀掃描。MRI掃描根據(jù)檢查部位應(yīng)用相應(yīng)的關(guān)節(jié)表面線圈。各關(guān)節(jié)的掃描序列及參數(shù)總結(jié)如下：T1WI快速自旋回波序列(TR 460～600 ms，TE 9～20 ms)；T2WI 快速自旋回波序列(TR 2800～4200 ms，TE 40～60 ms)；T2WI脂肪抑制快速自旋回波序列(TR 2600～3600 ms，TE 40～80 ms)；FOV 80 mm×80 mm～180 mm×180 mm，矩陣128×128～256×256，層厚3～4 mm，層間距0.5～1.0 mm；分別行冠狀位、矢狀位及軸位掃描。共51 個(gè)關(guān)節(jié)(RA 共45 個(gè)、GA 共6 個(gè))經(jīng)肘靜脈注射對(duì)比劑釓噴酸葡胺(Gd-DTPA)行動(dòng)態(tài)增強(qiáng)掃描，總量0.2 mL/kg 體質(zhì)量，注射速率3 mL/s，與注射同時(shí)進(jìn)行，行增強(qiáng)脂肪抑制T1WI冠矢及軸位掃描。

1.3 圖像分割

所有MR 圖像以DICOM 格式加載到計(jì)算機(jī)使用開源軟件ITK-SNAP (版本3.4)勾畫增生滑膜ROI。在T2WI-FS 序列上增生的滑膜組織界定為關(guān)節(jié)囊內(nèi)信號(hào)強(qiáng)度低于關(guān)節(jié)積液、高于肌肉組織信號(hào)強(qiáng)度的片絮狀影，放大圖像后可見其大致邊界。為確保ROI正確勾畫，首先由2名分別有3年(醫(yī)師1)和5年(醫(yī)師2)骨肌診斷經(jīng)驗(yàn)放射科醫(yī)師對(duì)照MR 增強(qiáng)圖像在T2WI-FS 序列上手動(dòng)勾畫3 個(gè)ROI (滑膜增生明顯處)，所有勾畫的滑膜ROI由1名高級(jí)職稱骨肌診斷醫(yī)師復(fù)核基本無誤后再進(jìn)行其他無增強(qiáng)圖像的T2WI-FS滑膜ROI勾畫，結(jié)果仍由高級(jí)醫(yī)師確認(rèn)符合。勾畫ROI時(shí)盡量放大圖像，避免重疊鄰近骨質(zhì)和關(guān)節(jié)積液。手動(dòng)分割示意圖見圖1、2。利用Python 軟件提取75 個(gè)紋理特征參數(shù)，包括19 個(gè)一階統(tǒng)計(jì)量(first order statistics)特征、24個(gè)灰度共生矩陣(GLCM)特征、16 個(gè)灰度長度矩陣(GLRLM)特征以及16個(gè)灰度形狀矩陣(GLSZM)特征。

圖1 女，39 歲，類風(fēng)濕關(guān)節(jié)炎。A、B：膝關(guān)節(jié)橫軸位T2WI-FS 及其對(duì)應(yīng)手動(dòng)勾畫ROI(紅色)圖像；C、D：膝關(guān)節(jié)橫軸位T1WI/FS 增強(qiáng)及其對(duì)應(yīng)手動(dòng)勾畫的ROI(深紅色)圖像 圖2 男，66 歲，痛風(fēng)性關(guān)節(jié)炎。A、B：腕關(guān)節(jié)冠狀位T2WI-FS 及其對(duì)應(yīng)手動(dòng)勾畫ROI(紅色)圖像；C、D：腕關(guān)節(jié)冠狀位T1WI-FS 增強(qiáng)及其對(duì)應(yīng)手動(dòng)勾畫ROI(深紅色)圖像Fig. 1 Female, 39 years old, rheumatoid arthritis.A, B were axial T2WI-FS of knee image and corresponding manually delineated ROI (red) image. C, D were axial enhanced T1WI-FS image and corresponding manually delineated ROI (dark red) image. Fig. 2 Male, 66 years old, gouty arthritis. A, B were coronal T2WI-FS image and corresponding manually delineated ROI(red)image.C,D were coronal enhanced T1WI-FS image and corresponding manually delineated ROI(dark red)image.

1.4 統(tǒng)計(jì)學(xué)分析

應(yīng)用SPSS 19.0統(tǒng)計(jì)軟件，患者性別定性資料進(jìn)行卡方檢驗(yàn)；年齡定量資料(參數(shù)符合正態(tài)分布及方差齊性檢驗(yàn))采用均數(shù)±標(biāo)準(zhǔn)差描述，進(jìn)行t檢驗(yàn)。采用組內(nèi)相關(guān)系數(shù)(ICC)評(píng)估2 名醫(yī)師提取增生滑膜ROI紋理特征的一致性。對(duì)2名醫(yī)師測得一致性較好(ICC＞0.75)的紋理特征進(jìn)行后續(xù)的統(tǒng)計(jì)分析，符合正態(tài)分布且方差齊性檢驗(yàn)的2組滑膜性關(guān)節(jié)炎紋理特征之間的差異行t檢驗(yàn)；反之，則采用Mann-Whitney U檢驗(yàn)進(jìn)行分析(所有參數(shù)均用2名醫(yī)師測得紋理特征的平均值進(jìn)行分析)。對(duì)各特征集具有統(tǒng)計(jì)學(xué)差異的紋理特征運(yùn)用因子分析方法，選取載荷指標(biāo)較大(＞0.75)的紋理特征，應(yīng)用Logistic 回歸分析行多特征聯(lián)合回歸分析，并進(jìn)行多特征集分析。分析相關(guān)紋理特征ROC 曲線的AUC，并比較其診斷效能。P＜0.05認(rèn)為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 一般資料分析

本組RA 患者年齡22～82 歲；GA 患者年齡21～75歲。如表2所示，兩組間性別比例差異有統(tǒng)計(jì)學(xué)意義，年齡分布差異無統(tǒng)計(jì)學(xué)意義。

表2 兩組患者的一般資料Tab.2 General data of the two groups

2.2 紋理特征提取一致性檢驗(yàn)結(jié)果分析

2名醫(yī)師提取紋理特征的ICC值范圍為0.620～0.915，篩選出一致性較好(ICC＞0.75)的55 個(gè)紋理特征進(jìn)行后續(xù)統(tǒng)計(jì)分析，包括15個(gè)一階統(tǒng)計(jì)量、18個(gè)灰度共生矩陣、12 個(gè)灰度長度矩陣和10 個(gè)灰度形狀矩陣特征。20個(gè)一致性較差(ICC≤0.75)的紋理特征中有8 個(gè)具有統(tǒng)計(jì)學(xué)差異(P＜0.05)，包括2 個(gè)一階統(tǒng)計(jì)量、3 個(gè)灰度共生矩陣、1 個(gè)灰度長度矩陣以及2 個(gè)灰度形狀矩陣特征。

2.3 紋理特征分析

一階統(tǒng)計(jì)量、灰度共生矩陣、灰度長度矩陣、灰度形狀矩陣紋理特征分析分別有10、9、6、5個(gè)紋理特征在RA 和GA 組中具有統(tǒng)計(jì)學(xué)差異(P＜0.05)。所有關(guān)節(jié)特征集聯(lián)合分析如表3所示，一階統(tǒng)計(jì)量特征紋理分析鑒別兩組關(guān)節(jié)滑膜病變的AUC 值最大(AUC=0.898)，僅略低于多特征集聯(lián)合診斷效能(AUC=0.902)。將兩組膝關(guān)節(jié)病例單獨(dú)提取的紋理特征集聯(lián)合分析如表4 所示，多特征集聯(lián)合特征紋理分析AUC 為0.900。各特征集AUC 值最大的紋理特征鑒別效能結(jié)果如表5 所示：一階統(tǒng)計(jì)量特征中的Minimum在所有單個(gè)紋理特征中AUC 值(0.836)最大，敏感度、特異度分別為0.778、0.770。

表3 所有關(guān)節(jié)特征集聯(lián)合ROC曲線分析鑒別效能Tab.3 Differential performance of texture feature sets combined ROC curve analysis of all joints

表4 膝關(guān)節(jié)紋理特征集聯(lián)合分析鑒別效能Tab.4 Differential performance of texture feature sets combined ROC curve analysis of the knees

表5 各特征集AUC值最大的紋理特征鑒別效能Tab.5 Differential performance of texture feature with maximum AUC value of each texture feature set

3 討論

高尿酸血癥是GA較特異的臨床表現(xiàn)，但文獻(xiàn)[4]報(bào)道GA 患者痛風(fēng)急性發(fā)作時(shí)血尿酸水平可為正常水平，而高尿酸血癥并不一定導(dǎo)致痛風(fēng)發(fā)作，尿酸水平正常不能排除痛風(fēng)，這類血尿酸水平不高患者臨床診斷較困難，常需要與RA鑒別。RA常多發(fā)、對(duì)稱侵犯手足小關(guān)節(jié)，但也有RA 以單關(guān)節(jié)起病并持續(xù)數(shù)年的相關(guān)文獻(xiàn)報(bào)道[12]。以往GA 相關(guān)研究多集中于雙能CT (DECT)對(duì)痛風(fēng)石檢出的相關(guān)應(yīng)用，但MRI 可以顯示GA 主要病理改變，并能監(jiān)測RA、GA 的治療反應(yīng)，在臨床應(yīng)用廣泛[13-14]。雖然文獻(xiàn)報(bào)道穿鑿樣骨質(zhì)破壞形式在GA 中較具特異度，但敏感度較低；大多數(shù)GA 滑膜炎T2WI 信號(hào)以等至低信號(hào)為主[4，15]，RA 病變早期滑膜組織由于炎細(xì)胞浸潤及小血管增生在T2WI 常呈高信號(hào)，然而文獻(xiàn)報(bào)道[16]RA病變晚期隨著滑膜纖維性血管翳形成，T2WI 也可呈低信號(hào)，另外二者都可伴有關(guān)節(jié)積液、骨髓水腫、關(guān)節(jié)周圍軟組織腫脹等表現(xiàn)，常規(guī)MRI 并不足以鑒別RA 和GA[17]。紋理分析通過高通量提取影像特征來更客觀和定量地描述病變特征，提供一種相對(duì)無創(chuàng)的鑒別診斷評(píng)價(jià)方式。目前尚未見紋理分析在滑膜性關(guān)節(jié)炎的相關(guān)研究報(bào)道[18]。釓劑由于存在腦沉積、腎功能損害等毒副作用，其應(yīng)用受到限制，而RA和GA患者需行多次MRI增強(qiáng)檢查進(jìn)行療效評(píng)估。近年文獻(xiàn)[19-20]報(bào)道了新興的影像學(xué)方法，如不應(yīng)用對(duì)比劑的脂肪抑制液體衰減反轉(zhuǎn)恢復(fù)(FLAIR-FS)序列和體素內(nèi)不相干運(yùn)動(dòng)擴(kuò)散加權(quán)成像(intravoxel incoherent motion diffusion-weighted imaging，IVIM-DWI)對(duì)于評(píng)估關(guān)節(jié)滑膜炎具有一定價(jià)值。另外文獻(xiàn)報(bào)道[10-11]增生滑膜與關(guān)節(jié)積液在T2WI-FS序列上對(duì)比度良好，本研究試圖利用MR T2WI抑脂單序列紋理分析來鑒別RA和GA。

本組RA 和GA 患者性別分布存在差異性。與以往研究[1，4]一致，本組RA 患者以女性多見，而GA 以男性多見。所有關(guān)節(jié)特征集聯(lián)合ROC 曲線分析結(jié)果顯示一階統(tǒng)計(jì)量、灰度共生矩陣、灰度長度矩陣、灰度形狀矩陣特征紋理分析的AUC 值分別為0.898、0.826、0.831、0.830，一階統(tǒng)計(jì)量紋理特征鑒別診斷效能較好，對(duì)應(yīng)的AUC 值、敏感度及特異度分別為0.898、0.852和0.802。加上膝關(guān)節(jié)單獨(dú)紋理分析鑒別效能也較高，表明基于MRI 常規(guī)T2WI 抑脂序列紋理分析對(duì)于RA 和GA 的有一定鑒別診斷價(jià)值。區(qū)分RA、GA 的主要紋理特征是測量病變信號(hào)的異質(zhì)性，如在一階統(tǒng)計(jì)量特征中的Entropy (熵)和灰度長度矩陣、灰度形狀矩陣中的GLN (灰度非一致性)紋理特征反映了ROI 的異質(zhì)性，其值越大提示灰度均勻性越差?；叶裙采仃囍械穆?lián)合能量(joint energy)度量圖像灰度分布均勻性，值越大表明圖像相鄰灰度分布越均勻。在本研究中，灰度長度矩陣和灰度形狀矩陣單個(gè)紋理特征ROC 曲線AUC 值最大對(duì)應(yīng)的紋理特征均為GLN，與之相應(yīng)GA 的GLN、Entropy 紋理特征值較RA更大，提示GA增生滑膜組織較RA的異質(zhì)性相對(duì)大。RA 和GA 滑膜炎病理組織學(xué)不同，在鏡下RA增生滑膜內(nèi)見較多炎細(xì)胞彌漫性浸潤，伴有淋巴濾泡樣小結(jié)節(jié)形成及漿細(xì)胞浸潤，并可見毛細(xì)血管增生[21]；GA 引起滑膜炎是對(duì)持續(xù)性尿酸鹽晶體沉積的反應(yīng)，其增生滑膜組織主要為結(jié)晶狀物體沉積[16]。綜合比較各特征集紋理分析，發(fā)現(xiàn)從一階統(tǒng)計(jì)量提取的紋理特征診斷效能最高(AUC=0.898)，僅略低于所有特征集聯(lián)合診斷效能(AUC=0.902)。

本研究具有以下不足：第一，不是全部病例均有增強(qiáng)圖像對(duì)照，部分圖像基于T2WI 手動(dòng)勾畫滑膜ROI，難以做到滑膜組織的準(zhǔn)確分割，文獻(xiàn)[22]報(bào)道應(yīng)用深度學(xué)習(xí)技術(shù)可實(shí)現(xiàn)腫瘤組織的精確檢測和自動(dòng)分割，后續(xù)研究可以采用類似半自動(dòng)或全自動(dòng)組織分割方法來更精確分割增生滑膜。第二，本研究為回顧性分析研究，不同場強(qiáng)儀器、不同部位MR圖像掃描參數(shù)無法做到精準(zhǔn)控制，對(duì)圖像本身紋理特征一致性存在偏差，難以做到完全一致。第三，本研究僅僅基于MR圖像紋理分析，并未結(jié)合臨床特征如性別、年齡和實(shí)驗(yàn)室檢查等，有待后續(xù)整理收集相關(guān)資料后運(yùn)用更深入的影像組學(xué)方法進(jìn)行綜合分析。

綜上所述，本研究顯示RA 與GA 的MRI 常規(guī)T2WI抑脂單序列紋理特征存在差異，紋理分析對(duì)兩者具有潛在鑒別價(jià)值；一階統(tǒng)計(jì)量提取的紋理特征在所有特征集中的鑒別診斷效能最高；一階統(tǒng)計(jì)量中的Minimum 為診斷效能最高的單個(gè)紋理特征。這或許能給慢性滑膜性關(guān)節(jié)炎長期隨訪療效評(píng)價(jià)中減少對(duì)比劑使用和副反應(yīng)發(fā)生提供一些借鑒。

作者利益沖突聲明：全體作者均聲明無利益沖突。