大家好呀!今天給大家介紹一篇2022年1月發表在Cancer Cell(IF:31.743)上的文章。作者對77例接受抗PD-1和抗CTLA-4治療的晚期皮膚黑色素瘤患者的全基因組數據,轉錄組數據和甲基化數據進行分析。結果表明TMB,新抗原負荷,IFNγ相關基因表達程序性死亡配體表達,PSMB8低甲基化和T細胞與免疫治療反應有關。結合TMB和IFNγ相關基因表達水平構建的模型可以準確預測免疫治療反應。高TMB和高IFNγ特征對免疫治療反應最好。
Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance
檢查點抑制劑治療黑色素瘤的多組學分析:鑒定預測反應和抵抗的標志物和生物紊亂標志物
流程圖:

結果:
1.接受抗PD-1免疫治療的黑色素瘤患者的臨床特征
作者收集抗PD-1免疫治療的77例黑色素瘤患者的樣本,包括單獨PD1和PD1聯合IPI治療。免疫治療中位時間為5.2個月,腫瘤組織主要來自皮下和區域淋巴結。77例患者中有66例患者發生于非肢端皮膚部位,11例發生隱匿性原發性黑色素瘤(表1)。年齡中位數為58,51例為男性。

2.高TMB,低結構變異負荷和高新抗原負荷與反應較好有關
作者首先研究了每個樣本的SNV,SV和CNA與免疫治療反應的關系(圖1A)。高TMB與免疫治療反應顯著相關(圖1B)。高SV和高CNA與免疫反應較差有關(圖1C和1D)。當根據免疫治療類型進行分類是,PD1治療反應與TMB顯著相關而IPI+PD1治療反應與SV顯著相關。高TMB和低SV與免疫治療的PFS改善顯著相關。此外,TMB與新抗原負荷具有良好相關性,新抗原負荷較高時反映較好(圖1E和1F)。

3.基因突變與免疫治療反應和抗性無關
接下來,作者研究免疫治療反應是否與特定的基因突變有關。作者發現TCGA分子亞型與免疫治療反應無關(圖1A)。作者對SNV和Indel進行分析然而并沒有發現具有統計學意義的基因。此外,作者還對已有報道與反應有關的基因中尋找特定突變然而與反應較好或較差沒有相關性。
4.干擾素和T細胞受體信號通路富集于免疫治療反應組但沒有特征與反應不了有關
對53例樣本進行RNAseq,對反應良好和反應不良進行差異分析鑒定到205個DEGs。在反應良好組中大量免疫基因上調表達包括CD274,STAT1,IFNγ,GZMA,PRF1,GZMB和CD8A。富集分析表明,反應良好組中只有IFNγ反應通路顯著富集(圖2A)。
對43個樣本進行Illumina EPIC甲基化測序,差異分析僅鑒定到5個差異甲基化CpG位點。由于這5個CpG位點均不在基因啟動子區域,作者接下來重點關注205個DEGs的甲基化水平。有5647個CpG位點注釋到192個差異基因中,173個基因中有1671個CpG位點與表達顯著相關,有41個基因中有5個或更多的CpG位點與表達負相關。其中,PSMB8的CpG位點與表達負相關,且在反應良好和反應不良之間具有顯著差異。
接下來,作者分析了與免疫治療反應有關的基因的表達水平。有一部分基因在反應良好和反應不良中差異表達(圖2B)。還有一些特征在反應良好和反應不良中具有顯著差異包括6-基因IFNγ特征(圖2C),18-基因IFNγ特征(圖2D),效應T細胞特征(圖2E),結合IFNγ/效應T細胞特征(圖2F),趨化因子特征(圖2G)和TIDE特征(圖2H)。

5.CD8 T細胞,M1巨噬細胞和細胞溶解活性對應答者免疫微環境的影響
作者使用CIBERSORT評價免疫細胞浸潤水平,鑒定反應良好和反應不良組的浸潤水平具有顯著差異的細胞類型(圖3A)。其中反應良好組中M1巨噬細胞和CD8 T細胞浸潤水平較高(圖3B和3C)。使用methylCibersort的甲基化水平評估免疫細胞浸潤水平,CD8 T細胞在反應良好組中浸潤水平較高(圖3D)。基于PRF1和GZMA的表達水平評估細胞溶解活性,細胞溶解國興在CD8 T細胞激活時上調,在反應良好組同樣上調(圖3E)。PRF1和GZMA的低甲基化與反應良好有關(圖3F)。免疫組化分析表明CD8 T細胞與反應良好有關(圖3G),瘤內CD8 細胞表達增加(圖3H)。此外,PD-L1表達水平在反應良好組中較高(圖3J)。

6.TMB和IFNγ-6打分可以預測免疫治療反應
由于TMB,IFNγ-6和CD8 T細胞與免疫治療反應有關,因此作者認為高TMB與免疫應答高有關,作者接下來研究TMB與IFNγ-6和CD8 T細胞的相關性。作者對53例樣本構建預測免疫治療反應的回歸模型。首先,作者基于TMB,SV,腫瘤純度和特征基因表達水平等構建單因素模型,其中TMB的作用最強。接下來,作者使用TMB和其他因子依次構建多因素模型。TMB結合IFNγ或TMB結合T細胞特征的模型具有相似的預測性能。TMB,IFNγ或TMB結合IFNγ的模型的ROC曲線如圖4A和5A所示。因此,TMB高和IFNγ-6高的患者對免疫治療反應較好(圖4B)且OS較長(圖4C)。接下來,作者使用一個包括60例黑色素瘤患者的獨立數據集來驗證TMB和IFNγ-6的性能,ROC曲線如圖4D和5B所示。

7.分析腫瘤中的生物學紊亂
對TMB/IFNγ-6模型不能準確分類的12例樣本進行臨床和分子檢測可以深入了解潛在的耐藥性機制和預測耐藥性反應(圖5A和5C)。12例樣本中有7例接受了PD1治療與可以準確預測的患者沒有相同的分子特征(圖5C)。第8例樣本具有潛在的耐藥生物學機制(JAK3突變)。8個樣本中有4個的TMB高且IFNγ-6低,4個的TMB低且IFNγ-6高。

結論:
作者對77例黑色素瘤患者的轉錄組數據,全基因組數據和DNA甲基化數據進行分析,鑒定與免疫治療反應顯著相關的生物標志物。結果表明,TMB,IFNγ-6和CD8 T細胞與免疫治療反應顯著相關。此外,結合TMB和IFNγ-6構建模型可以準確預測患者的免疫治療反應和總生存期。本篇文章的亮點在于第一,作者的樣本是自己收集的樣本而不是和常規生信分析一樣的TCGA數據集,第二作者的分析數據為多組學數據。如果你手里也有自測數據的話,不妨學習一下這篇文章的研究思路,說不定下一篇高分文章就是你的~~
參考文獻:
Newell F , Silva I , Johansson P A , et al. Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance. 2021.