儲海瑞，何劍，鄭屹峰，丁建平

（1.湖州市中心醫(yī)院放射科，浙江 湖州 313000；2.杭州師范大學(xué)附屬醫(yī)院放射科，浙江 杭州 310006）

腦膜瘤在顱內(nèi)原發(fā)腫瘤中發(fā)病率較高，僅次于膠質(zhì)瘤[1-2]。根據(jù)2021 年WHO 中樞神經(jīng)系統(tǒng)腫瘤的分類，腦膜瘤分為WHO Ⅰ、Ⅱ、Ⅲ級[3]。腦膜瘤的治療手段以手術(shù)切除為主，歐洲神經(jīng)腫瘤學(xué)會腦膜瘤診治指南中特別建議[4]：WHO Ⅰ級腦膜瘤如不伴有明顯臨床癥狀可定期隨訪，如伴有認(rèn)知障礙等癥狀則需要手術(shù)切除，術(shù)后每年進(jìn)行一次MR 檢查，隨訪5 年；WHO Ⅱ級或WHO Ⅲ級腦膜瘤則需要強(qiáng)制手術(shù)，并可根據(jù)病情輔助放療，或嘗試細(xì)胞毒藥物化療、激素治療、靶向治療等[5]，治療后每間隔半年一次MR 檢查隨訪5年，之后改為間隔1 年。由此可見，腦膜瘤的分級和臨床診療決策密切相關(guān)。因此，我們將利用多維度MR 特征共同構(gòu)建用于腦膜瘤分級的列線圖預(yù)測模型，希望通過量化評分提高對腦膜瘤的術(shù)前分級診斷能力。

1 資料與方法

1.1 患者資料

回顧性分析我院自2015 年1 月—2021 年12月符合以下標(biāo)準(zhǔn)的患者資料。納入標(biāo)準(zhǔn)：①M(fèi)R 檢查發(fā)現(xiàn)顱內(nèi)占位性病變，經(jīng)手術(shù)病理診斷為腦膜瘤并有明確病理分級的患者；②患者在行MR 檢查之前，顱內(nèi)病變沒有接受過任何治療（包括手術(shù)、穿刺活檢、放療、化療或激素治療等）。排除標(biāo)準(zhǔn)：①因各種原因造成的病例資料不全或失訪；②由于運(yùn)動或磁敏感偽影使圖像質(zhì)量欠佳不能達(dá)到研究要求。

共納入腦膜瘤患者127 例。訓(xùn)練集（n=89 例）包含低級別腦膜瘤71 例、高級別腦膜瘤18例，其中男性22例，女性67例，平均年齡（58±12）歲。驗(yàn)證集（n=38 例）包含低級別腦膜瘤30 例、高級別腦膜瘤8例，其中男性12例，女性26例，平均年齡（57±11）歲。性別、發(fā)病年齡在兩組間均無顯著性差異。

1.2 MRI 圖像采集與分析

采用 GE Discovery MR750 3.0T 磁共振掃描儀，8 通道頭部相控陣線圈，T1WI 序列參數(shù)為：TR 500 ms，TE 10 ms，NEX 2，矩陣320×256，層厚5 mm，F(xiàn)OV 24 cm。T2WI 參數(shù)：TR 5 200 ms，TE 95 ms，NEX 2，矩陣320×256，層厚5 mm，F(xiàn)OV 24 cm。使用釓噴酸葡胺對比劑進(jìn)行T1WI 序列增強(qiáng)掃描。

由2 位具有8 年和20 年工作經(jīng)驗(yàn)的放射科診斷醫(yī)師采用盲法評估每位患者的MR 特征，記錄結(jié)果不一致時經(jīng)討論達(dá)成共識。參考既往研究[6-7]結(jié)果中對鑒別不同級別腦膜瘤有價值的MR 特征，本次納入評估的常規(guī)MR 特征包括腫瘤最大徑、形態(tài)、囊變壞死、瘤周水腫（腫瘤周圍是否存在水腫帶）、瘤腦界面、鄰近骨質(zhì)改變。

1.3 影像組學(xué)分析

將原始影像資料導(dǎo)入ITK-SNAP（version 3.6，http://www.ITK SNAP.org）軟件，由上述兩位醫(yī)師在不知病理結(jié)果的情況下，對增強(qiáng)T1WI 序列顯示腫瘤的所有層面逐層勾畫感興趣區(qū)（ROI），ROI 盡量貼近腫瘤邊緣（圖1）。之后將原始圖像及分割的腫瘤ROI 文件同步導(dǎo)入Analysis Kit 軟件，提取圖像ROI內(nèi)的定量紋理參數(shù)。采用最小絕對收縮與選擇算子算法（Least absolute shrinkage and selection operator，LASSO）對特征進(jìn)行篩選，對篩選出的特征采用邏輯回歸的方法構(gòu)建影像組學(xué)標(biāo)簽并計算每位患者的影像組學(xué)分?jǐn)?shù)。

圖1 在增強(qiáng)T1WI 序列上勾畫腫瘤ROI。Figure 1.The tumor ROI was delineated on the enhanced T1WI.

1.4 列線圖模型的構(gòu)建和評估

在訓(xùn)練集納入具有獨(dú)立預(yù)測效能的MR 影像特征和影像組學(xué)分?jǐn)?shù)構(gòu)建列線圖預(yù)測模型（Model 1），利用具有預(yù)測效能的MR 影像特征構(gòu)建常規(guī)MR 特征預(yù)測模型（Model 2），利用影像組學(xué)分?jǐn)?shù)構(gòu)建影像組學(xué)預(yù)測模型（Model 3）。采用受試者工作特征（ROC）曲線評價3 個模型的預(yù)測效能，通過曲線下面積（Area under the curve，AUC）進(jìn)行比較，并在驗(yàn)證集進(jìn)行驗(yàn)證。

1.5 統(tǒng)計學(xué)分析

采用SPSS 23.0 和R 軟件（Version：3.6.1）進(jìn)行統(tǒng)計學(xué)分析。計量資料符合正態(tài)分布的采用獨(dú)立樣本t 檢驗(yàn)，非正態(tài)分布的采用Mann-Whitney U 檢驗(yàn)。分類資料以頻數(shù)表示并采用χ2檢驗(yàn)或Fisher 確切概率法。使用R 軟件的“rms”軟件包構(gòu)建列線圖模型，使用“pROC”軟件包繪制ROC 曲線并比較不同模型的AUC 值。

2 結(jié)果

2.1 常規(guī)MR 影像特征

低級別和高級別腦膜瘤的腫瘤最大徑、形態(tài)、囊變壞死、瘤周水腫在訓(xùn)練集和驗(yàn)證集均有統(tǒng)計學(xué)差異（P＜0.05），見表1。訓(xùn)練集的多因素回歸分析顯示腫瘤形態(tài)、囊變壞死、瘤周水腫是具有獨(dú)立預(yù)測效能的常規(guī)MR 因子，見表2。

表1 訓(xùn)練集和驗(yàn)證集的常規(guī)MRI 特征比較

表2 訓(xùn)練集多因素邏輯回歸分析

2.2 影像組學(xué)特征

采用Analysis Kit（AK）軟件通過量化分析，共提取396 個定量影像特征參數(shù)，包括直方圖特征、形態(tài)學(xué)特征、灰度共生矩陣特征、游程矩陣特征、灰度區(qū)域大小矩陣特征5 大類。使用LASSO 降維算法選擇篩選得到8 個有較大預(yù)測價值的特征參數(shù)，其中包括直方圖特征3 個、形態(tài)學(xué)特征1 個、灰度共生矩陣特征2 個、游程矩陣特征1 個、灰度區(qū)域大小矩陣特征1個，見表3。

表3 影像組學(xué)特征及系數(shù)值

2.3 列線圖模型

我們將腦膜瘤形態(tài)、囊變壞死、瘤周水腫和影像組學(xué)評分作為構(gòu)建列線圖模型的指標(biāo)，并成功繪制出用于預(yù)測腦膜瘤分級的列線圖（圖2）。列線圖預(yù)測模型（Model 1）在訓(xùn)練集的AUC 為0.965（95%CI：0.927～1.000），靈敏度為83.3%，特異度為98.6%；在驗(yàn)證集的AUC 為0.975（95%CI：0.931～1.000），靈敏度為100%，特異度為86.7%。常規(guī)MR特征預(yù)測模型（Model 2）在訓(xùn)練集的AUC 為0.817（95%CI：0.717～0.916），靈敏度為88.9%，特異度為56.3%；在驗(yàn)證集的AUC 為0.850（95%CI：0.703～0.997），靈敏度為62.5%，特異度為93.3%。影像組學(xué)預(yù)測模型（Model 3）在訓(xùn)練集的AUC 為0.909（95%CI：0.833～0.985），靈敏度為97.2%，特異度為72.2%；在驗(yàn)證集的AUC 為0.896（95%CI：0.785～1.000），靈敏度為100%，特異度為66.7%。結(jié)果顯示，列線圖模型的AUC 值明顯提高（圖3）。

圖2 預(yù)測腦膜瘤分級的列線圖。Figure 2.The nomogram for predicting meningiomas grading.

圖3 不同模型的ROC 曲線,圖3a 為訓(xùn)練集，圖3b 為驗(yàn)證集。Figure 3.ROC curves of different models,Figure 3a is of the training set and Figure 3b is of the validation set.

3 討論

3.1 構(gòu)建列線圖的各個參數(shù)與腫瘤之間的相關(guān)性

本研究通過統(tǒng)計學(xué)分析發(fā)現(xiàn)腦膜瘤形態(tài)、囊變壞死、瘤周水腫在不同級別的腦膜瘤之間有顯著差異性，與既往研究結(jié)果一致[6-7]。腫瘤細(xì)胞增生活躍易導(dǎo)致腫瘤向各個方向生長速率有差異，進(jìn)而致使腫瘤形態(tài)各異。同時，腫瘤生長過快易造成腫瘤內(nèi)部缺血缺氧導(dǎo)致壞死，進(jìn)而出現(xiàn)液化囊變，說明腫瘤形態(tài)和液化壞死與腫瘤的增殖快慢關(guān)聯(lián)[8-9]。有研究表明瘤周水腫具有預(yù)測腦膜瘤侵襲性的能力[10]。腦膜瘤雖起源于腦外，但WHO Ⅱ級腦膜瘤在生長過程中可浸潤周圍組織，使得腫瘤周圍血-腦屏障遭到破壞，并分泌較多的血管內(nèi)皮生長因子（Vascular endothelial growth factor，VEGF），VEGF 進(jìn)入瘤周誘導(dǎo)液體成分積聚于周圍膠質(zhì)細(xì)胞內(nèi)，從而出現(xiàn)瘤周水腫[11]。同時腫瘤對周圍血管壓迫可導(dǎo)致局部血供障礙，腫瘤細(xì)胞在局部組織環(huán)境缺氧的刺激下分泌更多VEGF，加劇瘤周水腫。

隨著人工智能的飛速進(jìn)展，影像組學(xué)的方法越來越廣泛的應(yīng)用于腫瘤的分級、病理分型、鑒別診斷、療效評估、預(yù)后預(yù)測等方面。虞芯儀等人通過對101 例經(jīng)病理證實(shí)的腦膜瘤進(jìn)行直方圖和灰度共生矩陣紋理分析，發(fā)現(xiàn)多個紋理參數(shù)在鑒別低級別與高級別腦膜瘤方面具有統(tǒng)計學(xué)意義，且偏度值、熵值具有較好的診斷效能[12]。李曉欣等[13]之前的一項研究表明，基于腫瘤全域測量的MR 信號強(qiáng)度直方圖參數(shù)有助于腦膜瘤分級，均勻性是最佳影像學(xué)分級診斷指標(biāo)。本次研究中，LASSO 最終選取的最佳紋理特征中有3 個直方圖特征，其中包括均勻性、頻率、最小強(qiáng)度，均勻性在所有參數(shù)中系數(shù)值最大。高級別腦膜瘤的實(shí)質(zhì)成分較復(fù)雜、不均勻，有較大的異質(zhì)性。因此，均勻性對不同級別的腦膜瘤表現(xiàn)出較高的鑒別效能?；叶裙采仃囂卣髦邪? 個熵類特征，熵反映了腫瘤內(nèi)部的空間分布情況，與腫瘤異質(zhì)性有一定的關(guān)系[14]。除此之外，我們還提取了更多、更全面的高階特征參數(shù)。游程矩陣主要反映紋理的粗糙程度和方向性，區(qū)域大小矩陣有效地描述了紋理的均勻性。這些高階特征能更好的反映不同級別腫瘤內(nèi)部的空間異質(zhì)性，為腦膜瘤分級診斷提供更豐富的信息[15]。

3.2 列線圖在預(yù)測腦膜瘤分級中的應(yīng)用價值

列線圖是一種復(fù)雜計算公式的圖像表現(xiàn)，它通過整合能夠決定結(jié)局變量的危險因素，將這些因素賦值并圖像化，最后通過計算分?jǐn)?shù)預(yù)測事件發(fā)生的可能性[16]。隨著精準(zhǔn)醫(yī)療的日益飛速發(fā)展，列線圖作為一種個體化預(yù)測方法，能針對多種疾病定制相應(yīng)的預(yù)測模型，對不同病情的病人進(jìn)行個性化評估并將結(jié)果可視化，已被越來越廣泛地運(yùn)用于疾病的鑒別診斷及分級預(yù)測，日益成為現(xiàn)代醫(yī)療決策的重要組成部分[17]。

本研究通過分析發(fā)現(xiàn)腫瘤形態(tài)不規(guī)則、囊變壞死、瘤周水腫被認(rèn)為是鑒別高低級別腦膜瘤的關(guān)鍵MR 特征，對腦膜瘤的分級具有一定診斷價值。此外，影像組學(xué)特征更好地反映不同級別腦膜瘤內(nèi)部的空間異質(zhì)性，對腦膜瘤分級診斷具有一定的應(yīng)用前景。我們通過將常規(guī)MR 特征和影像組學(xué)評分進(jìn)行整合，并以具有刻度的線段按各自比例平鋪陳列于列線圖之上，通過計算每個參數(shù)對應(yīng)的分?jǐn)?shù)值總和，并找到總分標(biāo)尺下方相對應(yīng)的預(yù)測風(fēng)險值，便可快速獲得針對腦膜瘤分級的預(yù)測概率。

3.3 本項研究的局限性

本研究存在以下幾點(diǎn)不足：①本研究的樣本量偏小，由于發(fā)病率、臨床診療策略等因素的限制，WHO Ⅱ級和Ⅲ級腦膜瘤的病例數(shù)量十分有限；②本研究為單中心回顧性研究，還缺乏外部驗(yàn)證及多中心大樣本數(shù)據(jù)的多重驗(yàn)證，來保證列線圖模型實(shí)際應(yīng)用于臨床的可靠性。因此，尋求多中心的合作，建立用于腦膜瘤研究的大樣本數(shù)據(jù)庫，將是未來的努力方向。

綜上所述，基于MR 圖像的影像組學(xué)列線圖模型具有良好的預(yù)測效能，對腦膜瘤的術(shù)前分級預(yù)測具有一定的指導(dǎo)價值。