MRI影像組學(xué)預(yù)測(cè)較低級(jí)別膠質(zhì)瘤免疫亞型的臨床研究

陳燕生,楊哲,沈卓,李文菲

低級(jí)別膠質(zhì)瘤(lower grade glioma,LGG)是最常見(jiàn)的原發(fā)性中樞神經(jīng)系統(tǒng)腫瘤之一,盡管神經(jīng)外科專(zhuān)家努力改善LGG患者的生存時(shí)間,但部分LGG患者由于腫瘤異質(zhì)性導(dǎo)致預(yù)后差異顯著[1]。有研究證實(shí)免疫療法顯著改善了LGG的預(yù)后,精確地識(shí)別LGG的腫瘤免疫表型對(duì)患者進(jìn)行免疫治療至關(guān)重要[2]。磁共振成像(magnetic resonance imaging,MRI)由于其無(wú)創(chuàng)性,被廣泛用于LGG的術(shù)前檢查,為臨床診斷提供重要信息。大量研究表明影像組學(xué)廣泛用于LGG的分子分型的診斷、預(yù)后預(yù)測(cè)及基因表達(dá)狀態(tài)[3-7],然而關(guān)于LGG不同影像表型的免疫特征研究并未報(bào)道。最新研究報(bào)道了影像組學(xué)與頭頸部腫瘤免疫亞型存在平行關(guān)系,同時(shí)多中心大樣本研究也證實(shí)影像組學(xué)可以預(yù)測(cè)腫瘤的免疫表型進(jìn)而為評(píng)估免疫檢查點(diǎn)如程序性死亡分子1 (programmed cell death protein-1,PD-1)和程序性死亡配體1(programmed death-ligand 1,PD-L1)治療的療效提供了一種有效的方法[8,9]。本研究旨在通過(guò)MRI影像組學(xué)特征預(yù)測(cè)LGG的預(yù)后進(jìn)而分析不同LGG免疫亞型與影像特征的相關(guān)性。

材料與方法

1．研究隊(duì)列

共納入199例患者來(lái)自于 TCIA(www.cancerimagingarchive.net)中的LGG 患者。納入和排除患者的依據(jù)由兩名具有5年神經(jīng)系統(tǒng)診斷經(jīng)驗(yàn)的高年資放射醫(yī)師決定:①術(shù)前未經(jīng)任何治療,包括但不限于化學(xué)治療、放射治療及保守治療;②同時(shí)具有軸位T2、T2FLAIR和T1增強(qiáng)圖像,層厚5 mm;③臨床信息應(yīng)同時(shí)具有性別、年齡、術(shù)后治療、隨訪等基線數(shù)據(jù)資料和腫瘤分級(jí)、亞型、IDH和1p/19q狀態(tài)等分子病理資料;④同時(shí)具有TCGA(www.cancergenome.nih.gov)數(shù)據(jù)庫(kù)的mRNA轉(zhuǎn)錄數(shù)據(jù)文件及Immport(https://immport.niaid.nih.gov)數(shù)據(jù)庫(kù)的免疫基因信息。 通過(guò)對(duì)TCIA影像數(shù)據(jù)患者與TCGA基因數(shù)據(jù)ID進(jìn)行匹配,最終本研究根據(jù)以上納入和排除標(biāo)準(zhǔn)共納入了133例同時(shí)具有影像和免疫相關(guān)基因數(shù)據(jù)的LGG患者。

2．LGG預(yù)后相關(guān)的MRI影像組學(xué)特征提取

由兩名高年資神經(jīng)組放射醫(yī)師在3D-slicer軟件(4.8.1版本)對(duì)LGG的病灶進(jìn)行識(shí)別并進(jìn)行腫瘤分割(圖1)。首先由低年資醫(yī)師在橫軸面增強(qiáng)T1序列中找到病灶位置,逐層勾畫(huà)感興趣區(qū),每個(gè)病例均在依據(jù)T2/FLAIR等序列排除瘤周水腫,利用General Registration插件進(jìn)行融合對(duì)比,進(jìn)一步提高瘤體與瘤周水腫的劃分。多平面重組(multiplanar reconstruction,MPR)圖像上手動(dòng)描繪膠質(zhì)瘤輪廓(圖1a),為了保證三維圖像(圖1b)的準(zhǔn)確性,將橫軸面MPR圖像和相應(yīng)的輪廓重建為冠狀面(圖1c)和矢狀面(圖1d)圖像輪廓作為參考。然后由另一位具有5年工作經(jīng)驗(yàn)的放射專(zhuān)家檢查膠質(zhì)瘤三維輪廓,隨時(shí)對(duì)ROI進(jìn)行修改和微調(diào)?；赑yradiomics軟件包總共提取851個(gè)放射組學(xué)特征,包括形狀特征(n=14)、一階特征(n=162)、灰度共生矩陣特征(n=216)、灰度游程矩陣特征(n=144)、灰度大小區(qū)域矩陣特征(n=144)、相鄰灰度差矩陣特征(n=45)、灰度依賴(lài)矩陣特征(n=126)?；赑yradiomics軟件包提供的算法提取感興趣區(qū)域內(nèi)影像組學(xué)特征。在進(jìn)行后續(xù)分析前,為了提高特征的穩(wěn)定性,納入符合圖像生物標(biāo)志物標(biāo)準(zhǔn)化倡議中提出的特征[10]。采用Z-score方法對(duì)特征進(jìn)行歸一化處理,兩位醫(yī)師分別對(duì)50名隨機(jī)選擇的患者完成了腫瘤分割和影像組學(xué)特征提取。然后醫(yī)師在1個(gè)月后重復(fù)勾畫(huà)了先前患者。為了減少手動(dòng)分割中的組間差異影響,筆者采用Spearman相關(guān)篩選觀察者內(nèi)和觀察者之間穩(wěn)定性較高的特征(intra correlation coefficient,ICC>0.85)。應(yīng)用10倍交叉驗(yàn)證的最小絕對(duì)收縮和選擇算子(least absolute shrinkage and selection operator,LASSO)比例風(fēng)險(xiǎn)回歸建立影像組學(xué)特征模型。根據(jù)最優(yōu)λ值構(gòu)建影像組學(xué)評(píng)分,通過(guò)以下算法計(jì)算每個(gè)患者的風(fēng)險(xiǎn)評(píng)分:影像組學(xué)評(píng)分(Radscore)=(β1×特征1值)+(β2×特征2值)+…+(βn × 特征n值)。根據(jù)臨界值(中位風(fēng)險(xiǎn)評(píng)分)將所有患者分為高危組和低危組。以上構(gòu)建影像組學(xué)評(píng)分方法,本研究團(tuán)隊(duì)在先前研究中詳細(xì)描述[11-12]。繪制Kaplan Meier生存曲線,比較高低風(fēng)險(xiǎn)組LGG患者生存期。

圖1 膠質(zhì)瘤三維感興趣區(qū)分割示意圖。a)橫軸面圖像勾畫(huà)范圍;b)三維體積示意圖;c、d)冠狀面和矢狀面圖像保證腫瘤勾畫(huà)的準(zhǔn)確性。

3．聚類(lèi)分析對(duì)LGG免疫亞型分組

基于先前發(fā)表的相關(guān)研究[13],使用單樣本基因富集分析(single-sample gene set enrichment analysis,ssGSEA)來(lái)量化133例 LGG患者中免疫細(xì)胞、功能或通路的活性或富集水平。根據(jù)ssGSEA的結(jié)果,利用RStudio中的“hclust”將LGG患者分為高免疫組(沙漠炎癥型)和低免疫組(免疫沙漠型),分別計(jì)算兩組中每個(gè)患者的腫瘤純度評(píng)分、免疫評(píng)分和基質(zhì)評(píng)分。最后分析放射組學(xué)評(píng)分與免疫亞型的相關(guān)性。

4．統(tǒng)計(jì)分析

采用R語(yǔ)言(版本4.02)進(jìn)行統(tǒng)計(jì)學(xué)分析和作圖。采用LASSO回歸篩選最佳影像組學(xué)特征構(gòu)建影像組學(xué)評(píng)分。采用Kruskal-Wallis檢驗(yàn)進(jìn)行不同免疫分型之間影像組學(xué)評(píng)分(Radscore)、ESTIMATE評(píng)分、免疫評(píng)分和基質(zhì)評(píng)分、腫瘤純度的差異性分析。Kaplan-Meier生存分析采用log-rank檢驗(yàn)。以P<0.05為差異具有統(tǒng)計(jì)學(xué)意義。

結(jié) 果

本研究共納入133例LGG患者,男68例,女65例,年齡18～75歲,平均(42.0±14.1)歲。異檸檬酸脫氫酶 (isocitrate dehydrogenase,IDH)野生型28例,突變型105例,其中野生型1p19q共缺失68例。經(jīng)過(guò)LASSO回歸篩選最佳特征,最終11個(gè)組學(xué)特征用來(lái)構(gòu)建影像組學(xué)評(píng)分預(yù)測(cè)LGG患者預(yù)后(表1),不同組學(xué)特征與LGG患者臨床病理相關(guān)性熱圖見(jiàn)圖2。膠質(zhì)瘤級(jí)別越高,影像評(píng)分越高。根據(jù)評(píng)分中位值分為高、低風(fēng)險(xiǎn)組,高風(fēng)險(xiǎn)組患者預(yù)后明顯差于低風(fēng)險(xiǎn)組,風(fēng)險(xiǎn)比(hazard ratios,HR)為3.3(95% CI:1.90～5.91,P<0.001),見(jiàn)圖3a。影像組學(xué)評(píng)分模型在預(yù)測(cè)LGG患者1年、3年、5年生存風(fēng)險(xiǎn)的準(zhǔn)確率分別為79%、78%和80%(圖3b)。

表1 構(gòu)建影像組學(xué)評(píng)分的詳細(xì)特征

圖2 預(yù)后相關(guān)組學(xué)特征與臨床病理變量復(fù)合型熱圖。共有11特征構(gòu)建影像組學(xué)評(píng)分,紅色越深代表評(píng)分越高,風(fēng)險(xiǎn)越高,紫色越深代表評(píng)分越低,風(fēng)險(xiǎn)越低。* 代表不同臨床病理特征影像組學(xué)評(píng)分存在顯著差異(*為P<0.05,**為P<0.01,***為P<0.001)。 圖3 a)影像組學(xué)評(píng)分與LGG患者生存期關(guān)系,高風(fēng)險(xiǎn)患者(患者)中位生存期明顯短于低風(fēng)險(xiǎn)組(紫色);b)影像組學(xué)評(píng)分預(yù)測(cè)LGG患者1、3、5年ROC曲線。

本研究根據(jù)ssGSEA來(lái)量化133例LGG患者中免疫細(xì)胞、功能或通路的活性或富集水平,由于7例患者免疫細(xì)胞信息缺失,共有126例患者進(jìn)行聚類(lèi)分析。利用聚類(lèi)分析將LGG患者分為兩種亞型(免疫炎癥型13例和免疫沙漠型113例),見(jiàn)圖4a。本研究還發(fā)現(xiàn)不同的免疫亞型在總生存期中具有不同預(yù)后意義,免疫沙漠型LGG患者預(yù)后較好(圖4b)。每個(gè)患者腫瘤純度、免疫評(píng)分、基質(zhì)評(píng)分和影像組學(xué)評(píng)分與免疫亞型的相關(guān)性熱圖見(jiàn)圖5。本研究發(fā)現(xiàn)不同免疫亞型的腫瘤純度、免疫評(píng)分、基質(zhì)評(píng)分及影像組學(xué)評(píng)分差異具有統(tǒng)計(jì)學(xué)意義(圖6,P<0.05),發(fā)現(xiàn)免疫炎癥型患者影像組學(xué)評(píng)分較高,腫瘤純度較低,異質(zhì)性較高,預(yù)后較差。

圖4 LGG患者不同免疫分型聚類(lèi)分析。a)聚類(lèi)分析將LGG患者分為兩種不同的免疫表型?！凹t色簇”代表免疫細(xì)胞浸潤(rùn)較多的免疫炎癥型腫瘤,“藍(lán)色簇”代表免疫細(xì)胞浸潤(rùn)較少的免疫沙漠型腫瘤;b)不同免疫分型K-M生存曲線,免疫炎癥型預(yù)后明顯差于免疫沙漠型。 圖5 影像組學(xué)評(píng)分、腫瘤純度、預(yù)估評(píng)分、免疫評(píng)分、基質(zhì)評(píng)分與不同免疫分型相關(guān)性熱圖。紅色代表富集分?jǐn)?shù)越高,綠色代表富集分?jǐn)?shù)越低(*為P<0.05,**為P<0.01,***為P<0.001)。

圖6 不同免疫亞型之間各指標(biāo)差異分析(*為P<0.05,**為P<0.01,***為P<0.001)。a)影像組學(xué)評(píng)分;b)免疫評(píng)分;c)基質(zhì)評(píng)分;d)腫瘤純度。

討 論

本研究首次基于放射基因組學(xué)全面分析影像組學(xué)與LGG預(yù)后及不同免疫分型的相關(guān)性?；?1個(gè)組學(xué)構(gòu)建的組學(xué)特征有效將LGG分為高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組,高風(fēng)險(xiǎn)患者預(yù)后明顯差于低風(fēng)險(xiǎn)組。影像組學(xué)評(píng)分與LGG病理類(lèi)型存在顯著相關(guān),腫瘤級(jí)別越高,組學(xué)評(píng)分越高,IDH突變型患者影像組學(xué)評(píng)分低于IDH野生型。為了探索預(yù)后差異的遺傳背景,分析了相關(guān)的轉(zhuǎn)錄組圖譜。聚類(lèi)分析將LGG患者分為兩種不同免疫亞型,LGG主要表現(xiàn)為免疫沙漠型,預(yù)后明顯好于免疫炎癥型。分析不同免疫分型與影像組學(xué)評(píng)分的相關(guān)性,結(jié)果表明影像組學(xué)評(píng)分與免疫炎癥型存在一定平行關(guān)系,免疫炎癥型患者影像組學(xué)評(píng)分較高,異質(zhì)性明顯高于免疫沙漠型。

膠質(zhì)瘤的時(shí)空異質(zhì)性不能局限于術(shù)后病理獲取,由于少量的活檢組織不能全面反映整個(gè)腫瘤的異質(zhì)性,因此無(wú)創(chuàng)性的MRI圖像分析在LGG膠質(zhì)瘤評(píng)估中扮演著重要角色。先前大量研究證實(shí)MRI組學(xué)特征在預(yù)測(cè)低級(jí)別膠質(zhì)瘤IDH分型及預(yù)后均表現(xiàn)出較高的準(zhǔn)確性[7,14],已經(jīng)成為精準(zhǔn)治療的橋梁;本研究基于MRI增強(qiáng)特征構(gòu)建的影像組學(xué)評(píng)分可以有效將LGG分進(jìn)行風(fēng)險(xiǎn)分層并且具有較高的預(yù)后預(yù)測(cè)能力,與先前研究基本一致[15];本研究同時(shí)發(fā)現(xiàn)不同病理分型、分子分型之間影像組學(xué)評(píng)分存在顯著差異,Li等[16]最新開(kāi)發(fā)基于臨床MRI掃描的放射組學(xué)特征可有效鑒別不同級(jí)別的LGG、不同IDH突變和1p/19q共缺失狀態(tài),準(zhǔn)確率為0.81～0.89,與本研究結(jié)果類(lèi)似。然而不同影像特征預(yù)后差異的生物學(xué)機(jī)制鮮見(jiàn)報(bào)道。

放射基因組學(xué)是整合基因組學(xué)、轉(zhuǎn)錄組學(xué)、表觀基因組學(xué)及影像組學(xué)等多組學(xué)數(shù)據(jù)全面分析腫瘤遺傳背景,已經(jīng)逐漸成為研究關(guān)注的熱點(diǎn)[17,18]。Sun等[17]研究團(tuán)隊(duì)納入435名膠質(zhì)母細(xì)胞瘤患者,研究發(fā)現(xiàn)基于13個(gè)影像組學(xué)特征構(gòu)建的影像組學(xué)評(píng)分與總生存期相關(guān)(HR為3.68) ,聚類(lèi)分析發(fā)現(xiàn)四種不同影像表型與免疫、增殖等不同生物學(xué)行為密切相關(guān)。Lai等[19]利用術(shù)前CT影像特征(腫瘤直徑及轉(zhuǎn)移角)可以有效預(yù)測(cè)胃癌染色體不穩(wěn)定狀態(tài),ROC曲線下面積為0.89。最新研究基于TCGA-TICA數(shù)據(jù)庫(kù)利用73名乳腺腫瘤患者乳腺動(dòng)態(tài)對(duì)比增強(qiáng)磁共振成像中提取的放射學(xué)特征來(lái)預(yù)測(cè)10種細(xì)胞類(lèi)型在乳腺癌病變中的浸潤(rùn),證實(shí)乳腺腫瘤患者影像特征和遺傳微環(huán)境的相關(guān)性[20]。這一系列研究表明放射基因組學(xué)具有預(yù)測(cè)腫瘤微環(huán)境的潛力,更多捕捉腫瘤內(nèi)基因水平的變化,影像特征與分子表型間存在一定的關(guān)聯(lián)性。本研究首次聯(lián)合TCGA-TICA數(shù)據(jù)庫(kù)進(jìn)一步探索133例LGG患者影像特征與免疫基因表達(dá)譜之間的關(guān)聯(lián)性,表明不同影像特征表現(xiàn)為不同的免疫亞型,表明不同的影像特征具有不同的生物學(xué)機(jī)制。本研究了發(fā)現(xiàn)大多數(shù)免疫細(xì)胞在LGG中可以獲得相對(duì)較低的 ssGSEA評(píng)分,表明大部分LGG 可能是一種“免疫沙漠型”腫瘤,少部分LGG患者屬于“免疫炎癥型”。還值得注意的是,免疫沙漠型LGG由于免疫細(xì)胞浸潤(rùn)較少,免疫療法對(duì)大部分LGG患者的療效較差[21,22],相反,免疫炎癥型LGG患者對(duì)免疫療法相對(duì)有效。同時(shí)本研究發(fā)現(xiàn)腫瘤純度低、免疫評(píng)分高、基質(zhì)評(píng)分高、ESTIMATE評(píng)分高為特征的LGG患者預(yù)后較差,上述結(jié)果已經(jīng)在相關(guān)研究證實(shí)[23]。最重要的是,研究結(jié)果表明免疫沙漠型LGG患者影像組學(xué)評(píng)分較低,預(yù)后較好,在一定程度上反映了免疫沙漠型患者異質(zhì)性較低,進(jìn)一步說(shuō)明影像特征對(duì)表征免疫亞型具有一定的潛力。

本研究的局限性:①盡管由于樣本量的限制,本研究的結(jié)論是初步的,但如果在更大的隊(duì)列中得到驗(yàn)證,這些發(fā)現(xiàn)對(duì)膠質(zhì)瘤免疫治療將具有重要的臨床意義;②本研究?jī)H僅提取了T1增強(qiáng)序列的影像組學(xué)特征、瘤周水腫和壞死區(qū)特征對(duì)低級(jí)別膠質(zhì)瘤也具有重要意義,另外感興趣區(qū)的勾畫(huà)存在一定的主觀性,未來(lái)自動(dòng)分割技術(shù)將有助于解決該弊端;③低級(jí)別膠質(zhì)瘤患者免疫治療方案缺乏,不能準(zhǔn)確預(yù)測(cè)免疫治療反應(yīng)。

綜上所述,本研究聯(lián)合LGG的MRI影像組學(xué)與基因轉(zhuǎn)錄組學(xué)特征,構(gòu)建影像組學(xué)評(píng)分預(yù)測(cè)LGG預(yù)后,采用聚類(lèi)分析方法評(píng)估不同免疫分型與影像組學(xué)評(píng)分的關(guān)聯(lián)性進(jìn)行嘗試性分析,結(jié)果表明影像組學(xué)評(píng)分與免疫炎癥型存在一定平行關(guān)系,免疫炎癥型患者影像組學(xué)評(píng)分較高,異質(zhì)性明顯高于免疫沙漠型,預(yù)后較差。放射基因組學(xué)分析通過(guò)影像組學(xué)特征無(wú)創(chuàng)性表征不同免疫亞型信息,可為免疫靶向治療開(kāi)拓新的視角。