預(yù)后模型到底做了什么發(fā)到十分呢?沒錯,他就是結(jié)合了臨床!
關(guān)鍵秘方只需要:一點點簡單的臨床知識+更多一點點的生信分析!想到就是賺到,快來學(xué)習(xí),打開思路吧!
今天分享的文章“基于MVI相關(guān)基因的肝細胞腫瘤預(yù)后模型的鑒定和驗證 Identification and Validation of a Prognostic Model Based on Three MVI-Related Genes in Hepatocellular Carcinoma”,這篇文章于2022年發(fā)表(Int J Biol Sci , IF 10.75)。
微血管浸潤( MVI)是指顯微鏡下鑒定的tumor侵犯血管,與tumor血道轉(zhuǎn)移密切相關(guān),是HCC肝內(nèi)或遠處轉(zhuǎn)移,術(shù)后復(fù)發(fā)和不良預(yù)后的獨立危險因素。目前仍沒MVI相關(guān)分子標(biāo)志物指導(dǎo)肝細胞癌診斷和治療。作者探索了新的MVI(微血管侵犯)相關(guān)生物標(biāo)志物,并建立了預(yù)測HCC預(yù)后的風(fēng)險評分模型,可以為HCC患者個性化治療提供基礎(chǔ)。
一、背景
肝細胞癌是世界上最常見的消化系統(tǒng)惡性腫瘤之一,發(fā)病隱匿,易復(fù)發(fā)和轉(zhuǎn)移,由于缺乏有效的預(yù)測和治療策略,預(yù)后往往較差,即使經(jīng)過治愈性治療,許多HCC患者仍會在5年內(nèi)出現(xiàn)腫瘤復(fù)發(fā)。因此,迫切需要開發(fā)新的預(yù)后評估方法來預(yù)測患者的臨床預(yù)后。
微血管浸潤(MVI)是指顯微鏡下可見血管腔內(nèi)內(nèi)皮細胞排列的腫瘤細胞團。腫瘤細胞可存在于肝門靜脈或肝靜脈,是肝內(nèi)或遠處轉(zhuǎn)移的潛在因素,并可促使術(shù)后復(fù)發(fā),MVI是一個獨立的組織病理學(xué)預(yù)后因素,與HCC患者在所有階段的生存相關(guān),還是術(shù)后早期復(fù)發(fā)和生存的關(guān)鍵因素和重要指標(biāo)。近年來,雖然MVI在肝癌治療選擇和預(yù)后評估中的臨床價值越來越受到重視,但由于MVI的診斷主要基于術(shù)后病理檢查,目前仍沒有推薦的MVI相關(guān)分子標(biāo)志物來預(yù)測HCC的預(yù)后,限制了MVI在指導(dǎo)肝癌診斷和治療中的應(yīng)用。
在本研究中,作者探索了新的MVI相關(guān)生物標(biāo)志物,并建立了預(yù)測HCC預(yù)后的風(fēng)險評分模型,旨在為HCC患者提供合適的治療方案。術(shù)前準(zhǔn)確預(yù)測MVI,有利于醫(yī)生在做出更合理的、個性化的手術(shù)范圍選擇術(shù)后隨訪策略制定,實現(xiàn)基于腫瘤生物學(xué)行為的個體化治療,有望改善患者的預(yù)后。
二、方法
1.數(shù)據(jù)采集和預(yù)處理
收集了兩組HCC RNA-seq數(shù)據(jù)。HCC的全基因組mRNA表達和臨床病理信息作為訓(xùn)練隊列從(TCGA LIHC下載)。共收集了95例MVI癌組織,210例無MVI癌組織和58例與臨床數(shù)據(jù)相匹配的正常組織。于國際癌癥基因組協(xié)會(ICGC LIRI-JP)下載驗證隊列。
2.TCGA隊列中差異表達MVI相關(guān)基因的鑒定
有或沒有MVI的癌組織分別用于與正常組織的差異分析。將DEG(具有MVI的癌組織與正常組織)中的差異基因作為MVI相關(guān)的差異表達基因。
3.Go和KEGG分析
將153個MVI相關(guān)基因應(yīng)用于GO和KEGG途徑富集分析。
4.HCC的預(yù)后signature開發(fā)和評估
基于訓(xùn)練集構(gòu)建預(yù)后signature,然后在驗證集中驗證其預(yù)測性能。首先進行單變量Cox比例風(fēng)險回歸分析確定了預(yù)后相關(guān)基因,以評估訓(xùn)練集中MVI相關(guān)基因與OS之間的相關(guān)性。利用最小絕對收縮和選擇算子(LASSO)懲罰Cox比例風(fēng)險回歸來減少模型的基因并限制解決過擬合問題的復(fù)雜性。逐步Cox回歸分析用于確定用于構(gòu)建風(fēng)險模型以預(yù)測HCC患者的預(yù)后。將中位風(fēng)險評分設(shè)置為臨界值,將患者分為高風(fēng)險組和低風(fēng)險組。Kaplan–Meier生存曲線用于低危組和高危組之間的生存比較。此外,建立了時間依賴性受試者工作曲線(包括1年,3年和5年生存率)來反映signature的敏感性和特異性。
5.signature的臨床特征
應(yīng)用Wilcoxon檢驗分析風(fēng)險評分與臨床特征(年齡,性別,階段,等級和分化等級)的關(guān)聯(lián),并且使用相同的方法將signature應(yīng)用于具有不同臨床特征的患者亞組,以區(qū)分低風(fēng)險和高風(fēng)險病例。
6.識別構(gòu)造signature的獨立預(yù)后作用
使用Cox回歸模型方法進行了多變量分析,以驗證其獨立于后價值。對具有不同臨床特征(包括年齡,性別和分期)的患者的生存率進行亞組分析,以進一步評估特征的預(yù)后表現(xiàn)。
7.預(yù)后基因signature中基因的外部驗證
基于驗證組,構(gòu)建Kaplan-Meier曲線來驗證基因signature的預(yù)測價值,多變量Cox分析評估基因signature在驗證組中的獨立預(yù)后作用。
三、結(jié)果
1.MVI陽性和MVI陰性HCC組織中的差異表達基因
比較MVI和非MVI組中上調(diào)和下調(diào)的DEG(與非tumor組織相比),得到差異表達 MVI相關(guān)基因(MVIRGs),用于后續(xù)建模。
MVIRGs的功能分析:主要參與血管內(nèi)皮的生長和微管的聚集,并且主要參與代謝相關(guān)途徑、tumor轉(zhuǎn)錄障礙和DNA復(fù)制途徑。(圖1)
2.預(yù)后模型的構(gòu)建和識別
單因素Cox篩選出與與預(yù)后相關(guān)的MVIRG,基于這些MVIRG進行LOSSO和Cox回歸,構(gòu)建基于MVI的預(yù)后模型,最后選擇三個關(guān)鍵的MVIRG(DBF4、ARG2和SLC16A3)構(gòu)建模型(圖2)。
3.三種MVIRGs的表達與HCC患者臨床參數(shù)之間的相關(guān)性。
為了進一步驗證該模型,作者分析了三種MVIRGs的表達與HCC患者臨床參數(shù)之間的相關(guān)性。KM曲線顯示,DBF4和SLC16A3基因的高表達與HCC患者的低OS顯著相關(guān),DBF4高表達、SLC16A3高表達與HCC的tumor分級存在顯著相關(guān)性(圖3)。
4.分析模型計算的風(fēng)險評分與臨床數(shù)據(jù)之間的相關(guān)性
高風(fēng)險評分與腫瘤分級、分期、T分期顯著相關(guān);多元Cox顯示風(fēng)險評分和分期是OS的獨立危險因素(圖4)。
5.使用附加數(shù)據(jù)隊列(ICGC)確認預(yù)后模型
將ICGC隊列中的HCC患者分為低或高風(fēng)險組,KM分析顯示兩組之間存在顯著差異;多因素Cox分析還顯示,風(fēng)險評分和分期是HCC患者預(yù)后的獨立危險因素,驗證了預(yù)測模型具有良好的預(yù)測價值(圖5)。
6.臨床組織樣本的驗證
HCC組織和成對的相鄰組織相比,顯示與3個MVIRGs(DBF4、ARG2和SLC16A3)相對應(yīng)的特征蛋白在腫瘤組織中染色更強烈,陽性區(qū)域更多(免疫組化);HCC樣本中三種MVIRGs(DBF4、ARG2和SLC16A3)的表達高于相鄰組織(qRT-PCR);隨訪數(shù)據(jù)進行分析,MVIRGs高表達者RFS和OS更低(圖6)。
7.列線圖模型的構(gòu)建與識別
將常規(guī)臨床參數(shù)引入并構(gòu)建列線圖模型,c-index為0.697,校準(zhǔn)曲線驗證了該模型具有良好的預(yù)測值(圖7)。
四、討論
MVI是HCC患者生存的重要預(yù)測因子,基于此,作者將MVI相關(guān)基因與常規(guī)臨床參數(shù)相結(jié)合,構(gòu)建了一個預(yù)后模型。富集分析表明,這些MVIRGs主要富集血管內(nèi)皮的生長和微管的積累,主要通過代謝途徑,腫瘤轉(zhuǎn)錄紊亂途徑和DNA復(fù)制途徑影響HCC的預(yù)后。這些結(jié)果與以前的研究結(jié)果一致。例如,MVI意味著在內(nèi)皮細胞排列的血管腔中可見腫瘤細胞簇,這表明它與血管內(nèi)皮的生長和微管的積累有關(guān)。此外,MVI與上皮-間質(zhì)轉(zhuǎn)化(EMT)密切相關(guān),EMT涉及多種機制,包括代謝變化,轉(zhuǎn)錄調(diào)控和表觀遺傳異常,表明MVI也與這些途徑有關(guān)。
用于建模的三個關(guān)鍵的MVIRGs(DBF4,ARG2和SLC16A3)與腫瘤密切相關(guān),其中DBF4編碼的蛋白質(zhì)具有模塊化結(jié)構(gòu),在DNA復(fù)制、復(fù)制檢查點激活、減數(shù)分裂、有絲分裂退出、跨病變合成和組蛋白穩(wěn)態(tài)中起著關(guān)鍵作用。過去的一些研究曾指出,較高的ASK/Dbf4表達黑色素瘤與較低的無復(fù)發(fā)生存率相關(guān),抑制Cdc7/Dbf4激酶活性會影響癌細胞中MCM2上的特異性磷酸化位點。ARG2(精氨酸酶2)是一種蛋白質(zhì)編碼基因,其活性的增加通常與更晚期的疾病和不良的臨床預(yù)后有關(guān),此外與腫瘤免疫反應(yīng)密切相關(guān)。最近的研究發(fā)現(xiàn),ARG2的過度表達是多種癌癥類型的預(yù)后不良因素,包括胰腺癌、甲狀腺腫瘤、胃癌、神經(jīng)母細胞瘤、頭頸部鱗狀細胞癌和急性髓系白血?。ˋML)。其他的一些研究發(fā)現(xiàn),ARG2的缺失可導(dǎo)致氨的積累和肥胖相關(guān)胰腺癌生長的抑制,線粒體ARG2是CD8+T細胞功能和抗腫瘤功效的細胞自主調(diào)節(jié)因子。有趣的是,一些學(xué)者發(fā)現(xiàn)ARG2還具有抑癌作用。它可以通過生物合成輔因子磷酸吡哆醛的消耗和多胺毒性的增加來抑制腎癌的進展。
SLC16A3(溶質(zhì)載體家族16成員3)是蛋白質(zhì)編碼基因。由SLC16A3基因編碼的單羧酸轉(zhuǎn)運蛋白4(MCT4)可以催化單羧酸的質(zhì)子連接轉(zhuǎn)運,單羧酸參與體內(nèi)許多代謝過程,并產(chǎn)生抗凋亡作用。MCT4在糖酵解組織中顯著存在,如缺氧癌細胞,在某些癌細胞中過表達,并在癌細胞生長和增殖中起關(guān)鍵作用。已有的研究表明,SLC16A3是肺腺癌患者預(yù)后不良和轉(zhuǎn)移的獨立指標(biāo),是胰腺癌代謝過程的關(guān)鍵調(diào)節(jié)因子,SLC16A3對腫瘤診斷和預(yù)后具有良好的預(yù)測能力,并且參與甲狀腺癌和透明細胞腎細胞癌中SLC16A3啟動子DNA甲基化和/或MCT4蛋白水平的表觀遺傳調(diào)控機制。
綜上,這三個關(guān)鍵的MVIRGs(DBF4,ARG2和SLC16A3)與腫瘤及其預(yù)后密切相關(guān),這也在一定程度上證明了選擇這三個MVIRGs建立預(yù)后模型的正確性。