Prognosis-Related Molecular Subtypes and Immune Features Associated with Hepatocellular Carcinoma

與預(yù)后相關(guān)的分子亞型和與肝細胞癌相關(guān)的免疫特征

發(fā)表期刊：Cancers (Basel)

發(fā)表日期：2022 Nov 21

影響因子：6.575

DOI:? 10.3390/cancers14225721

一、研究背景

????????原發(fā)性肝癌是全球癌癥相關(guān)死亡的第四大原因窖式，它表現(xiàn)出顯著的組織學(xué)和生物學(xué)異質(zhì)性抵栈。肝細胞癌（HCC）是最常見的類型，占原發(fā)性肝癌的90%以上扎唾。HCC的主要風(fēng)險因素包括感染丙型肝炎或乙型肝炎病毒盔性、過度飲酒划址、吸煙和飲食。遺傳風(fēng)險因素可能與這些環(huán)境風(fēng)險因素相互作用掖桦，如CTNNB1（編碼β-catenin）本昏、TP53和AXINI基因的突變。

????????由于HCC的異質(zhì)性很強枪汪，身體不同區(qū)域的腫瘤可能代表不同的亞型涌穆。這些亞型在代謝和信號通路上可能有所不同，導(dǎo)致患者生存期的差異雀久。

二宿稀、材料與方法

1、數(shù)據(jù)來源

1）從廣西醫(yī)科大學(xué)附屬腫瘤醫(yī)院收集了六對HCC和配對的癌旁正常組織樣本赖捌，這些樣本是在不同的臨床階段被病理診斷為HCC的患者

2）從TCGA數(shù)據(jù)庫中的肝癌（LIHC）數(shù)據(jù)集中下載了371個原發(fā)腫瘤樣本祝沸、3個復(fù)發(fā)腫瘤樣本和50個正常組織樣本的基因表達譜

3）在GEO數(shù)據(jù)庫，獲得了GSE14520（225個HCC組織和220個非腫瘤組織）越庇、GSE76427（52個相鄰的非腫瘤組織和115個HCC組織）罩锐、GSE25097（268個HCC組織和243個鄰近的非腫瘤組織）、GSE138178（49個HCC組織和配對的鄰近非腫瘤組織）卤唉、GSE84006（38個原發(fā)性HCC組織和配對的鄰近非腫瘤肝組織）涩惑、GSE136319數(shù)據(jù)集（秘魯患者的47個腫瘤和47個非腫瘤肝臟組織的DNA甲基化數(shù)據(jù)）

2、分析流程

流程圖

三桑驱、實驗結(jié)果

01 - HCC中的DEGs及其功能富集度

????????為了確定與HCC預(yù)后相關(guān)的基因竭恬，作者首先利用TCGA跛蛋、GSE76427、GSE25097和GSE14520數(shù)據(jù)集的數(shù)據(jù)進行了差異分析（圖2A,B）痊硕。在四個數(shù)據(jù)集重疊的4036個交叉的DEGs中问芬，有2058個上調(diào)，1272個下調(diào)寿桨，706個DEGs表達方向不同（圖2C）。功能富集分析表明强戴，激活了各種與HCC相關(guān)的途徑亭螟，如P53信號通路、色氨酸代謝以及初級膽汁酸生物合成（圖2D）骑歹。重疊的DEGs可能參與了細胞的氨基酸分解過程预烙、羧酸分解過程和其他代謝過程（圖2E）。

圖2 HCC的DEGs及其在生物過程和KEGG途徑中的富集程度

????????GSEA顯示道媚，DEGs與細胞周期扁掸、錯配修復(fù)和DNA復(fù)制正相關(guān)，但與礦物質(zhì)吸收最域、PPAR信號傳導(dǎo)以及補體和凝血級聯(lián)負相關(guān)谴分。在GSE14520中，預(yù)測HCC患者5年OS的AUC是58%镀脂。Cox回歸和KM分析也顯示牺蹄，3330個重疊的DEG中，有217個與HCC預(yù)后密切相關(guān)薄翅，而使用Metascape的進一步富集分析顯示沙兰，這些預(yù)后基因在小分子分解代謝、小分子生物合成和有絲分裂細胞周期中明顯富集翘魄。

????????為了評估TCGA和GSE14520中預(yù)后相關(guān)基因的診斷價值鼎天，作者計算了它們的AUC值，并從217個AUC>0.80的預(yù)后基因中選出138個（圖3A）暑竟。然后對這138個基因進行單變量Cox分析斋射，得到10個相對重要性大于0.2的生存相關(guān)基因（圖3B）。隨后的LASSO回歸確定了七個特征基因：SORBS2, DHRS1, SLC16A2, RCL1, IGFALS, GNA14和FANCI但荤，預(yù)測5年OS的AUC為0.744（圖3C-E）绩鸣。然后用GMM模型將7個特征基因的表達數(shù)據(jù)整合成3個群組127個組合，并選擇AUC最高的群組來識別預(yù)測HCC復(fù)發(fā)的強勢特征基因纱兑。經(jīng)GMM分類器測定呀闻，127個組合中的1個特征基因的平均準確率為0.9901（圖3F）。通過單變量Cox回歸分析進一步調(diào)查特征基因的獨立預(yù)后價值潜慎，表明FANCI在兩個數(shù)據(jù)集中都與不良的OS顯著相關(guān)（圖3G捡多，H）蓖康，表明它可能作為HCC復(fù)發(fā)的一個新的預(yù)測生物標志物。

圖3 單變量Cox和LASSO回歸

02 - 基于特征基因的預(yù)后風(fēng)險評分作為HCC的預(yù)后工具

????????根據(jù)風(fēng)險評分的中位數(shù)將HCC患者分為高風(fēng)險組和低風(fēng)險組（圖4A,B）垒手，在GSE14520和TCGA數(shù)據(jù)中蒜焊，其預(yù)測1、3科贬、5年OS的AUC值均大于0.65（圖4C,D）泳梆。使用整合了特征基因和臨床病理風(fēng)險因素的nomogram對OS預(yù)測進行了量化（圖4E）。校準圖也顯示榜掌，nomogram與理想模型相比表現(xiàn)良好（圖4F）优妙。

圖4 基于特征基因的風(fēng)險評分可作為HCC的預(yù)后工具

????????為了確定HCC的分子亞型，根據(jù)217個預(yù)后相關(guān)基因?qū)SE14520和TCGA數(shù)據(jù)庫的HCC樣本進行了NMF聚類（圖5A-C）憎账。能夠?qū)CC患者分為兩個分子亞型（圖5C）：預(yù)后差的HCC的C1亞型套硼，和OS良好的C2亞型（圖5E）。還發(fā)現(xiàn)GSE14520的剪影寬度值為0.85（圖5D）胞皱，TCGA的剪影寬度值為0.88邪意，這表明HCC樣本與兩個不同亞型之間有很好的相關(guān)性。

圖5 基于GSE14520數(shù)據(jù)集的HCC亞型確定

03 - HCC亞型對免疫療法和化療藥物的敏感性

????????GSE14520和TCGA中兩種HCC亞型的臨床數(shù)據(jù)分布比較表明反砌，兩種亞型的年齡沒有明顯差異雾鬼，但男性比女性更容易患兩種亞型的疾病。此外宴树，根據(jù)GSE14520（圖6B）和TCGA（圖6C）的數(shù)據(jù)呆贿，C1亞型顯示出比C2亞型更短的生存期（圖6A），以及對CTLA4-R治療更大的反應(yīng)性森渐。相反做入，TCGA中的C2亞型對抗癌藥物ZM.447439和AG.14699的反應(yīng)性明顯大于C1亞型（圖6D）。

圖6 肝細胞癌C1和C2亞型的藥物敏感性和免疫治療反應(yīng)

04 - 干性指數(shù)和FANCI的表達

????????根據(jù)干性指數(shù)值對HCC樣本進行排名同衣，顯示其臨床人口學(xué)特征與mRNAi明顯相關(guān)（圖S4A）竟块。FANCI與干性指數(shù)正相關(guān)（圖S4B），它在GSE138178和GSE84006數(shù)據(jù)集中都是上調(diào)的（圖S4C耐齐，D）浪秘。利用Oncomine數(shù)據(jù)庫進一步分析FANCI在各種癌癥類型中的mRNA表達，發(fā)現(xiàn)FANCI在肝癌中與正常組織相比明顯上調(diào)（圖S4E）埠况。這些結(jié)果被TCGA RNA測序數(shù)據(jù)在TIMER中證實耸携，表明FANCI在腫瘤中的水平明顯高于鄰近的正常組織（圖S4F）。

圖S4 臨床人口學(xué)特征與基于mRNAsi和FANCI在HCC的表達有關(guān)

????????為了探索FANCI在HCC預(yù)后中的作用辕翰，作者進行了時間-ROC分析夺衍，結(jié)果顯示，在GSE14520數(shù)據(jù)集中預(yù)測5年OS的AUC最高（圖7A）喜命，而在TCGA中預(yù)測1沟沙、3河劝、5年OS的AUC都大于0.60（圖7B）。FANCI的AUC在所有四個數(shù)據(jù)集中都高于0.90（圖7C）矛紫。富集分析顯示赎瞎，F(xiàn)ANCI與纖維蛋白溶解的調(diào)節(jié)、環(huán)氧酶P450途徑和蛋白激活級聯(lián)呈正相關(guān)颊咬，而與錯配修復(fù)务甥、中心體分離和轉(zhuǎn)體合成呈負相關(guān)（圖7D）。此外喳篇，F(xiàn)ANCI與細胞周期敞临、DNA復(fù)制和蛋白酶體正相關(guān)，而與FOXO杭隙、IL-17和p53信號通路負相關(guān)（圖7E）。

圖7 FANCI與生物和KEGG途徑的時間-ROC分析和關(guān)聯(lián)性

????????還發(fā)現(xiàn)在人類蛋白質(zhì)圖譜數(shù)據(jù)庫中因妙，F(xiàn)ANCI在HCC組織中的表達高于對照組（圖8A）痰憎。同樣，6名HCC患者的轉(zhuǎn)錄組測序結(jié)果顯示攀涵，與鄰近患者相比铣耘，F(xiàn)ANCI在HCC中高表達，特別是在腫瘤進展的患者2（RFS 66天）和5（RFS 248天）（圖8B-D）以故。同時蜗细，F(xiàn)ANCI在以前的研究中被確定為Wnt信號的下流基因，調(diào)節(jié)HCC的早期復(fù)發(fā)怒详。這些結(jié)果提示炉媒，它可能正向調(diào)節(jié)范可尼貧血癥的途徑（圖8E）。

圖8 FANCI在HCC中的潛在調(diào)節(jié)機制

05 - 免疫細胞浸潤昆烁、體細胞突變和DNA甲基化

????????為了探索免疫細胞浸潤在HCC中的潛在臨床意義吊骤，作者確定了所有四個數(shù)據(jù)集的浸潤水平（圖9A）。與對照組相比静尼，HCC樣本中T輔助2型（Th2）白粉、T輔助和前樹突狀細胞明顯上調(diào)，而T細胞和細胞毒性細胞與C1亞型明顯相關(guān)（圖9B）鼠渺。還分析了HCC組織中24種免疫細胞類型之間的相關(guān)性（圖9C）鸭巴，并構(gòu)建了四個集群，顯示了彼此之間的正負相關(guān)（圖9D）拦盹。此外鹃祖，樹突狀細胞和Th12細胞與七個特征基因相關(guān)，而Th2細胞和FANCI之間有明顯的相關(guān)性（圖9E）普舆。CIBERSORT算法還顯示惯豆，大多數(shù)浸潤的免疫細胞是T細胞（圖9E）池磁。

圖9 免疫細胞之間的浸潤和關(guān)聯(lián)性

????????根據(jù)顯示至少20倍覆蓋率的測序數(shù)據(jù)，在364名HCC患者中鑒定了體細胞單核苷酸變異楷兽。發(fā)現(xiàn)突變使腫瘤抑制基因TP53（占所有患者的31%）地熄、AXINI（8%）和RB1（5%）失活（圖S5A）。此外芯杀，發(fā)現(xiàn)FANCI在HCC患者中有兩個位點發(fā)生突變端考，這些突變可能會影響蛋白功能（圖S5B）。進一步發(fā)現(xiàn)FANCI的mRNA水平與蛋白水平呈負相關(guān)揭厚，這意味著該基因會受到甲基化的影響（圖S5D）却特。對GSE136319中秘魯肝癌患者HCC和非腫瘤組織之間差異甲基化的基因分析表明，甲基化修飾被調(diào)查（圖S5C）筛圆。此外裂明，MEXPRESS工具顯示，TCGA中FANCI啟動子區(qū)域的甲基化水平在HCC樣本中明顯高于正常組織太援。

圖S5 肝細胞癌的基因組圖譜和DNA甲基化變化

四闽晦、結(jié)論

????????作者利用生物信息學(xué)比較了腫瘤和非腫瘤組織之間的基因表達，以確定HCC分子亞型和生物標志物來幫助預(yù)測預(yù)后提岔。此外仙蛉，使用單樣本基因集富集分析（ssGSEA）來確定免疫浸潤的水平并描述腫瘤免疫微環(huán)境的特征。本研究確定了HCC患者的兩種分子亞型碱蒙，并檢測了可能作為HCC潛在生物標志物和治療目標的關(guān)鍵基因荠瘪。