從現(xiàn)在開(kāi)始解讀2018年4月的TCGA終結(jié)系列文章，因?yàn)樯婕暗降拇a知識(shí)不多，我就不親自出手了，會(huì)邀請(qǐng)我的幾個(gè)學(xué)生來(lái)一一解讀，當(dāng)然，也歡迎廣大粉絲投稿，運(yùn)氣的話有機(jī)會(huì)得到我的關(guān)注和指導(dǎo)。

免疫微環(huán)境在腫瘤發(fā)生、侵襲、轉(zhuǎn)移及治療耐受過(guò)程中占據(jù)重要位置，預(yù)示著我們可以從免疫微環(huán)境參與腫瘤治療耐受的這一視角出發(fā)，為尋找新的抗腫瘤途徑開(kāi)辟新的思路，下面我們就一起來(lái)嘗試了解這篇癌癥免疫圖譜文章吧。

2018年4月Immunity雜志上發(fā)表了文章The Immune Landscape of Cancer ，由34個(gè)單位共同合作完成。文章對(duì)TCGA中33種癌癥，超過(guò)10,000個(gè)腫瘤樣本進(jìn)行了免疫原性分析，將所有腫瘤分成6種免疫亞型，即：

wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-β dominant

進(jìn)一步分析：

• 不同亞型間巨噬細(xì)胞或淋巴細(xì)胞特征、Th1:Th2細(xì)胞比例、腫瘤異質(zhì)性程度、非整倍性、新抗原負(fù)荷程度、細(xì)胞增殖、免疫調(diào)節(jié)基因的表達(dá)、預(yù)后等指標(biāo)的差異。

• 與免疫相關(guān)的驅(qū)動(dòng)突變

• 參與腫瘤免疫的細(xì)胞內(nèi)和細(xì)胞間調(diào)控網(wǎng)絡(luò)（調(diào)控網(wǎng)絡(luò)包含轉(zhuǎn)錄、microRNA、拷貝數(shù)、表觀遺傳信息）

分析流程

step1. 下載TCGA平臺(tái)33種癌癥類型的6種分子平臺(tái)數(shù)據(jù)（mRNA、microRNA、外顯子測(cè)序數(shù)據(jù)；甲基化、拷貝數(shù)、蛋白質(zhì)芯片數(shù)據(jù)）

step2. 為全部非血液類癌癥進(jìn)行免疫表達(dá)特征（160種）評(píng)分，將TCGA全部腫瘤聚類成6種免疫亞型

step3. 根據(jù)DNA甲基化、mRNA和H&E圖像分析確定白細(xì)胞類型及比例（用于評(píng)估免疫浸潤(rùn)細(xì)胞的組成）

step4. 生存分析評(píng)估免疫亞型與預(yù)后的關(guān)系

step5. 分別以TCGA腫瘤類型、TCGA定義的分子亞型、6種免疫亞型分組，比較新抗原、病毒RNA表達(dá)、T細(xì)胞受體(TCR)和B細(xì)胞受體(BCR)、免疫調(diào)節(jié)因子(IM)表達(dá)與調(diào)控等特征（評(píng)估影響免疫原性和免疫浸潤(rùn)的因素間的關(guān)系）

免疫原性（immunogenicity)：能引起機(jī)體產(chǎn)生免疫應(yīng)答的物質(zhì)

腫瘤浸潤(rùn)免疫細(xì)胞包括：T cells, B cells, natural killer cells, macrophages, neutrophils, dendritic cells, mast cells, eosinophils, basophils等

step6. 鑒定與免疫浸潤(rùn)相關(guān)的體細(xì)胞變異（包括pathway, 拷貝數(shù)變異，driver mutation）

step7. 評(píng)估性別與祖源是否對(duì)特定腫瘤類型的免疫反應(yīng)有影響

step8. 鑒定控制腫瘤免疫應(yīng)答的細(xì)胞內(nèi)調(diào)控網(wǎng)絡(luò)、參與形成腫瘤免疫微環(huán)境(TME)的細(xì)胞外通訊網(wǎng)絡(luò)

癌癥的免疫亞型

計(jì)算160種免疫特征間的相關(guān)系數(shù)，聚類分析，得到5個(gè)核心模塊免疫表達(dá)特征（巨噬細(xì)胞 , 免疫浸潤(rùn)淋巴細(xì)胞, TGF-β response , IFN-γ response, wound healing ) → 據(jù)此將TCGA中的非血液腫瘤聚成6種免疫亞型(C1-C6) → 6種免疫亞型中這5種免疫特征的分布情況。（圖1A）

6種免疫亞型中腫瘤樣本表現(xiàn)出的其它關(guān)鍵免疫特征（圖1B,C）

每種腫瘤類型中免疫亞型的組成情況，圖1D。

Figure 1. Immune Subtypes in Cancer

腫瘤免疫浸潤(rùn)的組成

不同免疫亞型的主要免疫浸潤(rùn)細(xì)胞分布，圖2A

白細(xì)胞比例(leukocyte fraction,LF)在不同免疫亞型間（圖1C）和不同腫瘤類型間（圖2B）存在差異。

腫瘤基質(zhì)比例與白細(xì)胞比例的相關(guān)性在腫瘤類型中不同，圖2C顯示2種代表性的腫瘤。（該相關(guān)性在免疫亞型的差別見(jiàn)補(bǔ)充圖2B）

基質(zhì)細(xì)胞：器官中的結(jié)締組織細(xì)胞，為器官中實(shí)質(zhì)細(xì)胞提供支持和營(yíng)養(yǎng)。

已知基質(zhì)細(xì)胞和腫瘤細(xì)胞間的相互作用在腫瘤生長(zhǎng)和進(jìn)展中起主要作用?；|(zhì)細(xì)胞可以提供腫瘤細(xì)胞可以生長(zhǎng)的細(xì)胞外基質(zhì)，調(diào)節(jié)細(xì)胞生長(zhǎng)因子，還可能通過(guò)一氧化氮的產(chǎn)生限制T細(xì)胞增殖，從而阻礙免疫能力。

通過(guò)TCGA的H&E染色圖片分析估計(jì)出具有浸潤(rùn)淋巴細(xì)胞的腫瘤區(qū)域的空間比例，按免疫亞型查看，發(fā)現(xiàn)C2亞型的淋巴細(xì)胞浸潤(rùn)區(qū)比例最高，圖2D。

Figure 2. Composition of the Tumor Immune Infiltrate

腫瘤免疫反應(yīng)與預(yù)后

6種免疫亞型OS生存分析（圖3A）：C3預(yù)后最佳；C1,C2雖然免疫成分多，但預(yù)后并不好；C4,C6免疫特征混雜，預(yù)后最差。（PFI結(jié)果見(jiàn)補(bǔ)充圖3A）

計(jì)算5種免疫表達(dá)特征分?jǐn)?shù)與OS的一致性指數(shù)，紅色代表高風(fēng)險(xiǎn)。按免疫亞型和腫瘤類型分組（圖3B），淋巴細(xì)胞特征與C1,C2的預(yù)后改善相關(guān)；5種特征中任何一個(gè)增加都會(huì)導(dǎo)致C3預(yù)后更差。
計(jì)算輔助T細(xì)胞（Th1,Th2,Th17）與OS的一致性指數(shù)，按免疫亞型分組（圖3C），Th17細(xì)胞的增加與多數(shù)亞型的預(yù)后改善相關(guān)。

輔助T細(xì)胞（T helper cells, Th）,是一種T細(xì)胞（白細(xì)胞的一種），它的表面有抗原受體，可以識(shí)別抗原呈遞細(xì)胞的MHC-II 類分子呈獻(xiàn)的抗原片段。 輔助T細(xì)胞主要可區(qū)分為T(mén)h1 Th2 Th17 及 Thαβ等四種。

使用CoxPH建立免疫特征分?jǐn)?shù)高低的生存模型，分?jǐn)?shù)高低與生存顯著相關(guān)，圖3D是驗(yàn)證集的KM-polt，圖3E顯示模型的預(yù)測(cè)準(zhǔn)確度。圖3F說(shuō)明添加腫瘤stage和tpye信息后模型預(yù)測(cè)更準(zhǔn)確。

補(bǔ)充圖3E,F討論了其它免疫特征與OS的關(guān)系，某些腫瘤類型中免疫亞型與OS的關(guān)系

Figure 3. Immune Response and Prognostics

免疫反應(yīng)與體細(xì)胞變異

6中免疫亞型的免疫浸潤(rùn)（白細(xì)胞比例LF）與DNA損傷（包括CNV負(fù)荷、非整倍性、雜合性丟失LOH、同源重組缺陷HRD、腫瘤內(nèi)異質(zhì)性ITH）的相關(guān)系數(shù)熱圖，圖4A。發(fā)現(xiàn)LF與C6和C2的相關(guān)性最強(qiáng)，從整體來(lái)看與非整倍性、LOH、HRD及突變負(fù)荷成正相關(guān)。

圖B顯示平均LF值在1-22號(hào)染色體上的觀察值與預(yù)估值間的差異，左右圖分別為拷貝數(shù)擴(kuò)增和缺失，黑色部分代表差異顯著。揭示了LF與拷貝數(shù)變異的相關(guān)性，如：chr1p擴(kuò)增，LF高于預(yù)估值，缺失則LF低于預(yù)估值，據(jù)此進(jìn)一步將該位置上的重要基因的CNV與免疫浸潤(rùn)相聯(lián)系。

基因的平均LF預(yù)估值：每種腫瘤的平均IF的平均值由每種疾病類型中存在的“擴(kuò)增/缺失”樣品的數(shù)量加權(quán)。

補(bǔ)充圖4探討了腫瘤內(nèi)異質(zhì)性ITH與預(yù)后、LF、免疫亞型間的關(guān)系。

將299個(gè)癌癥驅(qū)動(dòng)基因的突變與免疫亞型相關(guān)聯(lián)，發(fā)現(xiàn)33個(gè)顯著相關(guān)的基因，圖C。

圖D為L(zhǎng)F-驅(qū)動(dòng)基因變異火山圖，x軸為基因變異與LF的多因素相關(guān)性（考慮了腫瘤類型和突變數(shù)量），y軸為顯著性，橙色即為與LF顯著相關(guān)的驅(qū)動(dòng)基因變異。其中23個(gè)驅(qū)動(dòng)基因變異與LF增加相關(guān)，12個(gè)與LF降低相關(guān)。

補(bǔ)充圖4C分析了體細(xì)胞變異（突變+CNVs）對(duì)8種癌信號(hào)通路的影響。

Figure 4. Immune Response and Genome State

人口統(tǒng)計(jì)學(xué)因素及遺傳變異與免疫反應(yīng)

考察了不同腫瘤類型中免疫細(xì)胞與性別和祖源的相關(guān)性，圖4E。發(fā)現(xiàn)在某些腫瘤中，女性的PD-L1表達(dá)高于男性，非洲血統(tǒng)患者PD-L1表達(dá)較低。

免疫原性調(diào)查

文章中的免疫原性主要討論的是新抗原負(fù)荷，即SNV和Indel突變產(chǎn)生的預(yù)測(cè)能與pMHC結(jié)合并誘導(dǎo)免疫反應(yīng)的肽。

MHC proteins(pMHC)：指一些列細(xì)胞表面蛋白，主要作用是結(jié)合抗原并將其呈遞至細(xì)胞表面，使其被T-cell識(shí)別，即幫助免疫系統(tǒng)識(shí)別外來(lái)物質(zhì)。

突變與pMHC：圖5A為與突變數(shù)量相關(guān)的pMHCs數(shù)量分布，標(biāo)注了源于超過(guò)40個(gè)突變的pMHCs。99.8的pMHC由單突變產(chǎn)生，0.2%是由不同突變產(chǎn)生的相同pMHC。

常見(jiàn)pMHC：圖5B為y軸表達(dá)共有pMHC的腫瘤數(shù)，標(biāo)注了產(chǎn)生最常見(jiàn)pMHC的已知癌基因(PIK3CA, KRAS, BRAF等)

補(bǔ)充圖5A,B討論了pMHC數(shù)與預(yù)后的關(guān)系：在多數(shù)腫瘤中pMHC數(shù)與生存無(wú)關(guān)，但在某些免疫亞型中，pMHC數(shù)與PFI預(yù)后相關(guān)。

補(bǔ)充圖5C討論了病毒與腫瘤中免疫特征的關(guān)系，病毒含量會(huì)影響免疫細(xì)胞的組成。

Figure 5. The Tumor-Immune Interface

癌癥中的適應(yīng)性免疫受體庫(kù)

適應(yīng)性免疫Adaptive immune：又稱獲得性免疫，對(duì)特定病原體具有高度特異性。獲得性免疫系統(tǒng)的細(xì)胞是T和B淋巴細(xì)胞。B細(xì)胞主要參與體液免疫，而T細(xì)胞則參與細(xì)胞免疫。B細(xì)胞和T細(xì)胞都攜帶能夠識(shí)別特定靶標(biāo)的受體分子。

T細(xì)胞負(fù)責(zé)識(shí)別“非自身”靶標(biāo)，如病原體，但需要一種被稱為主要組織相容性復(fù)合體（MHC）呈遞抗原之后才能實(shí)現(xiàn)。B細(xì)胞抗原特異性受體則是位于B細(xì)胞表面的抗體，負(fù)責(zé)識(shí)別病原體，不需要經(jīng)過(guò)抗原呈遞。

免疫球蛋白immunoglobulin，又稱抗體，由適應(yīng)性免疫系統(tǒng)的B細(xì)胞分泌。

根據(jù)RNA-seq數(shù)據(jù)評(píng)估TCR(T cell receptor)α和β、免疫球蛋白重鏈和輕鏈庫(kù)。

發(fā)現(xiàn)TCR多樣性值因免疫亞型而不同(圖5C)，C6和C2的多樣性最高。（補(bǔ)充圖5D展示了不同腫瘤類型的TCR多樣性）

評(píng)估CDR3的α和β鏈，以確定具有相同TCR的患者頻率。圖5C為CDR3α-CDR3β pairs，顏色代表p值，圓的大小代表sample數(shù)目，至少在2種腫瘤中出現(xiàn)的α-β pairs有2812個(gè)。

檢驗(yàn)SNV特異性pMHC-CDR3 pairs，圖5D展示了206個(gè)pMHC-CDR3α pairs和196個(gè)pMHC-CDR3β pairs。（說(shuō)明這些患者中抗原受體多，抗原少，多數(shù)T細(xì)胞抗原反應(yīng)是由公共抗原介導(dǎo)）

抗原受體通常由位于重鏈和輕鏈上的兩個(gè)可變結(jié)構(gòu)域組成，每個(gè)可變結(jié)構(gòu)域上有三個(gè)非連續(xù)排列的CDR（CDR1，CDR2和CDR3）

Complementarity-determining regions, CDR: 互補(bǔ)決定區(qū)，是免疫球蛋白（抗體）和T細(xì)胞受體中可變鏈的一部分，分別由B細(xì)胞和T細(xì)胞產(chǎn)生，與其特異性抗原結(jié)合。作為分子中變化最大的部分，CDR對(duì)淋巴細(xì)胞產(chǎn)生的抗原特異性的多樣性至關(guān)重要。

在某些腫瘤類型中，TCR多樣性與PFI改善相關(guān)，見(jiàn)補(bǔ)充圖5F。

免疫球蛋白重鏈和輕鏈的多樣性模式與TCR類似。

免疫調(diào)節(jié)劑的調(diào)控

免疫調(diào)節(jié)劑(Immunomodulators, IM)基因的情況可以反映癌細(xì)胞對(duì)腫瘤微環(huán)境(TME)的調(diào)節(jié)。

IM基因的表達(dá)在免疫亞型中不同（圖6A），圖6B為亞型間差異最大的基因EDNRB和CXCL10的表達(dá)值。

許多IM的甲基化程度與表達(dá)值負(fù)相關(guān)，即表觀遺傳沉默（圖6A），如圖6C的C3亞型中的CD40基因。

294個(gè)miRNA可能是IM基因的調(diào)節(jié)因子，補(bǔ)充圖6C。

拷貝數(shù)改變影響多個(gè)IM，并在免疫亞型中存在差異。C1,C2亞型中IM基因頻繁擴(kuò)增和缺失，C3和C5中IM基因變異較少（圖6A）。KIR2DL3等基因在C5中明顯缺失，CD40等基因常在C1中發(fā)生擴(kuò)增（圖6D）。

圖6A熱圖中的擴(kuò)增/缺失頻率：obs是指特定亞型中IM基因擴(kuò)增/缺失頻率，exp是指全部樣本中擴(kuò)增/缺失頻率的差異，熱圖反映的是obs與exp的差異

Figure 6. Regulation of Immunomodulators

免疫反應(yīng)的網(wǎng)絡(luò)調(diào)控

免疫應(yīng)答由腫瘤細(xì)胞、免疫細(xì)胞、基質(zhì)細(xì)胞中的細(xì)胞內(nèi)分子網(wǎng)絡(luò)狀態(tài)和細(xì)胞外通訊網(wǎng)絡(luò)（包括細(xì)胞間直接互作、通過(guò)可溶性蛋白或細(xì)胞因子通訊）構(gòu)成。

根據(jù)FANTOM5中記錄的信息獲得配體 - 受體，細(xì)胞 - 受體和細(xì)胞 - 配體對(duì)的網(wǎng)絡(luò)，圖A,B,C分別是提取的IFM-γ, TGF-β, T細(xì)胞的細(xì)胞外通訊子網(wǎng)絡(luò)網(wǎng)絡(luò)。展示了免疫亞型間的一些關(guān)聯(lián)，例如C2,C3中，CD4 T細(xì)胞、CD8 T細(xì)胞和NK細(xì)胞都與IFM-γ及CCL5的表達(dá)相關(guān)；C2,C3,C6中，多種細(xì)胞都與TGF-β的表達(dá)相關(guān)……

構(gòu)建泛癌癥轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)，關(guān)聯(lián)基因組事件 - 轉(zhuǎn)錄調(diào)控因子 - 下游靶基因 - 免疫浸潤(rùn)-預(yù)后 。使用Master Regulator(MR)和SYGNAL兩種互補(bǔ)的方法，產(chǎn)生了2個(gè)轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)：

• 圖D是Master Regulator(MR)-Pan-Immune網(wǎng)絡(luò)，26個(gè)MRs(橙色填充)代表基因表達(dá)與LF相關(guān)的hub，連接15個(gè)上游driver events(橙色環(huán))。

• SYGNAL-PamImmune網(wǎng)絡(luò)，包含了171個(gè)富集且與LF相關(guān)的的IMs，主要展示下游調(diào)控網(wǎng)絡(luò)（突變 - 轉(zhuǎn)錄因子 - miRNA間的聯(lián)系）。（圖7E展示部分SYGNAL-PamImmune網(wǎng)絡(luò)）

MRs網(wǎng)絡(luò)圖的構(gòu)建過(guò)程：根據(jù)下游靶基因的表達(dá)變化，推斷影響蛋白質(zhì)活性的轉(zhuǎn)錄因子，為MRs → 根據(jù)蛋白互作尋找與MRs相關(guān)的發(fā)生變異的蛋白cMRs(與LF的相關(guān)系數(shù)高的優(yōu)先) → 尋找與cMRs相關(guān)的mutation和CNVs事件。

SYGNAL網(wǎng)絡(luò)構(gòu)建過(guò)程：用高表達(dá)的基因作為input，得到具有轉(zhuǎn)錄因子和miRNA共調(diào)節(jié)的基因(bicluster) → 過(guò)濾篩選得到171個(gè)與免疫浸潤(rùn)和免疫調(diào)節(jié)均顯著相關(guān)的biclusters→ 添加biclusters的轉(zhuǎn)錄因子和miRNA調(diào)控信息 → 整合這些調(diào)控關(guān)系，連接成網(wǎng)絡(luò)

（網(wǎng)絡(luò)圖構(gòu)建步驟比較復(fù)雜，大致應(yīng)該是這個(gè)意思）

綜合2個(gè)調(diào)控網(wǎng)絡(luò)，MR及SYGNAL中有7個(gè)共享的TFs，調(diào)節(jié)27個(gè)IM基因，推斷AKAP9, HRAS, KRAS和PREX2基因突變影響IMs。

圖7E展示免疫亞型-IM基因表達(dá)調(diào)控-腫瘤類型間的關(guān)系，可以鑒定不同免疫亞型內(nèi)的IM基因調(diào)控因子。例如，C1和C2中的IM基因共同受BCL5B, ETV7, IRF1, IRF2, IRF4, PRDM1, 和SPIB調(diào)控；C3中IM受KLF15 and miR-141-3p調(diào)控……

某些腫瘤中，bicluster的表達(dá)增加與預(yù)后相關(guān)，在KIRC, LGG, LUSC和READ中預(yù)后差，在SKCM中預(yù)后好。

Figure 7. Predicted Networks Modulating the Immune Response to Tumors

結(jié)果總結(jié)

文章利用TCGA數(shù)據(jù)構(gòu)建了癌癥的免疫基因組學(xué)圖譜：

• 根據(jù)5種免疫表達(dá)特征將TCGA的腫瘤分為6種免疫亞型

• 通過(guò)多種方法評(píng)估腫瘤樣本中的免疫組成，包括：根據(jù)基因表達(dá)和甲基化數(shù)據(jù)估計(jì)免疫細(xì)胞組成、預(yù)測(cè)來(lái)自突變和HLA-typing的新抗原-MHC pairs、根據(jù)RNA測(cè)序數(shù)據(jù)評(píng)估BCR和TCR庫(kù)。

• 比較癌癥和免疫亞型的免疫組成

• 鑒定與腫瘤微環(huán)境相關(guān)的體細(xì)胞變異

• 構(gòu)建影響腫瘤微環(huán)境的分子調(diào)控網(wǎng)絡(luò)及細(xì)胞間通訊網(wǎng)絡(luò)

• 聯(lián)系免疫特征與生存，OS和PFI在腫瘤間或一些腫瘤內(nèi)的免疫亞型中存在差異

數(shù)據(jù)和結(jié)果的補(bǔ)充表格下載地址：https://gdc./about-data/publications/panimmune （劃重點(diǎn)?。。。?/p>

文章很清晰，做了pan-cancer的免疫亞型分類，然后把能用到的所有基因組指標(biāo)和免疫指標(biāo)都進(jìn)行了分組比較。文章的數(shù)據(jù)量很大，用到的方法很多，方法介紹也寫(xiě)得很詳細(xì)，需要時(shí)可以再詳細(xì)了解；補(bǔ)充表格的數(shù)據(jù)有很多信息可以直接拿來(lái)用，避免重復(fù)分析一次。