掃碼聯(lián)系客服
各位小伙伴，大家好呀！ 今天給大家?guī)淼氖且黄趩渭?xì)胞RNA測(cè)序數(shù)據(jù)分析人類胃癌器官特異性轉(zhuǎn)移轉(zhuǎn)錄異質(zhì)性的文章；文章分析思路整體來說都是偏常規(guī)的分析，但是亮點(diǎn)在于作者將腫瘤微環(huán)境中的每一種免疫細(xì)胞都進(jìn)行了細(xì)致入微的刻畫，非常適合大家在分析數(shù)據(jù)時(shí)開拓自己的思路，快來和小編一起看看吧~
研究背景
胃癌(GC)是全球第二大癌癥死亡原因，每年有103萬例，而亞洲胃癌發(fā)病率最高，占全球病例的70%。大多數(shù)GC患者被診斷為晚期和轉(zhuǎn)移期，再進(jìn)行藥物治療時(shí)往往患者已經(jīng)產(chǎn)生了耐藥性，因此，手術(shù)治療已經(jīng)不再是治療胃癌的唯一方式。
腫瘤細(xì)胞可通過血液、淋巴系統(tǒng)和腹腔內(nèi)擴(kuò)散侵入遠(yuǎn)處部位以進(jìn)行轉(zhuǎn)移，在轉(zhuǎn)移過程中，腫瘤細(xì)胞通過激活特定的基因和通路，向特定的器官進(jìn)化適應(yīng)度并表現(xiàn)出趨向性，不同器官的轉(zhuǎn)移是胃癌診斷和治療面臨的最困難和最緊迫的問題。因此，精確繪制器官特異性轉(zhuǎn)移特征是至關(guān)重要的，特別是對(duì)于不同轉(zhuǎn)移的特異性治療策略的開發(fā)以及臨床診斷的潛在生物標(biāo)志物的識(shí)別。然而，胃癌向不同器官轉(zhuǎn)移的特點(diǎn)在很大程度上尚不清楚。
因此，作者收集了包括轉(zhuǎn)移和肺轉(zhuǎn)移共6名胃癌患者進(jìn)行單細(xì)胞RNA測(cè)序，刻畫了胃癌患者腫瘤間和腫瘤內(nèi)部微環(huán)境異質(zhì)性的問題，更好地分析了胃癌不同器官轉(zhuǎn)移的模式。
二、主要結(jié)果
Result1: 刻畫原發(fā)性和轉(zhuǎn)移性GC的單細(xì)胞轉(zhuǎn)錄圖譜
作者對(duì)6名患者的10份新鮮人體組織樣本進(jìn)行了scRNA序列分析，包括三份原發(fā)腫瘤樣本（PT；即PT1、PT2和PT3）、一份相鄰的非腫瘤樣本（NT；即NT1）和六份轉(zhuǎn)移樣本（M）。在轉(zhuǎn)移樣本中，作者在胃鏡檢查、活檢或手術(shù)切除中獲得了兩個(gè)肝轉(zhuǎn)移樣本（Li：即Li1和Li2）、兩個(gè)淋巴結(jié)轉(zhuǎn)移樣本（LN：即LN1和LN2）、一個(gè)腹膜轉(zhuǎn)移樣本（P：P1）和一個(gè)卵巢轉(zhuǎn)移樣本（O：O1）。其中，PT1和Li1來自患者1；PT2和NT1來自患者2；O1來自患者3；PT3、Li2和LN1來自患者4；LN2來自患者5；P1來自患者6。PT1、PT2、PT3、Li1、Li2、LN1和P1是腸道組織學(xué)GC樣品，而LN2和O1是組織混合型GC樣本。
在經(jīng)過質(zhì)量篩選后，作者共檢測(cè)到42968個(gè)細(xì)胞以用于下游分析（圖1A、B）。在進(jìn)行降維和聚類以及運(yùn)用Marker基因進(jìn)行細(xì)胞注釋后共確定了七種細(xì)胞類型，分別為上皮細(xì)胞、基質(zhì)細(xì)胞、T細(xì)胞、B細(xì)胞、自然殺傷細(xì)胞和髓系細(xì)胞等（圖1C）并展示了注釋每種細(xì)胞類型所運(yùn)用到的Marker（圖1D、E）。作者對(duì)不同的細(xì)胞類型進(jìn)行數(shù)目上的統(tǒng)計(jì)，發(fā)現(xiàn)由于不同細(xì)胞類型的分離效率不同，免疫細(xì)胞在所有樣本中所占比例最大（90.4%），尤其是淋巴結(jié)和肝臟。
最后作者對(duì)于原發(fā)性腫瘤和轉(zhuǎn)移性腫瘤進(jìn)行了分析，發(fā)現(xiàn)每種細(xì)胞的比例差異很大（圖1F），揭示了轉(zhuǎn)移性腫瘤與原發(fā)性腫瘤相比較腫瘤間的異質(zhì)性。
圖1：GC原發(fā)腫瘤、轉(zhuǎn)移灶和鄰近非腫瘤樣本的單細(xì)胞轉(zhuǎn)錄組譜
Result2:惡性上皮細(xì)胞的四個(gè)亞群特征與腫瘤轉(zhuǎn)移和生存時(shí)間的關(guān)系
作者借鑒之前發(fā)表在GUT文章上面的一個(gè)工作，通過計(jì)算每個(gè)惡性細(xì)胞中的腫瘤患者差異基因的表達(dá)豐度（score值）以此定義每個(gè)細(xì)胞的惡性程度。運(yùn)用該方法作者共識(shí)別到1615個(gè)惡性細(xì)胞和128個(gè)非惡性上皮細(xì)胞。其中大部分的非惡性上皮細(xì)胞來自NT（圖2A）。此外，作者發(fā)現(xiàn)與非惡性上皮細(xì)胞相比，惡性上皮細(xì)胞中幾種腫瘤特異性基因（CLDN4、CLDN7、TFF3）的表達(dá)顯著上調(diào)，而非惡性上皮細(xì)胞強(qiáng)烈表達(dá)與胃粘液和消化酶分泌相關(guān)的幾個(gè)基因（MUC5AC、GKN1、PGC、LIPF），驗(yàn)證了注釋細(xì)胞惡性程度的準(zhǔn)確性。接著作者對(duì)PT和M樣本中惡性細(xì)胞進(jìn)行聚類分析，共顯示了四個(gè)亞聚類G0-G3（圖2B），同時(shí)對(duì)不同樣本以及不同轉(zhuǎn)移瘤中檢測(cè)到惡性細(xì)胞的數(shù)目進(jìn)行統(tǒng)計(jì)以展示惡性細(xì)胞的轉(zhuǎn)錄異質(zhì)性（圖2C）。接下來，作者進(jìn)一步通過功能分析等定義每個(gè)惡性細(xì)胞亞群的潛在功能。
具體來說：作者發(fā)現(xiàn)血管生成的關(guān)鍵因子（VEGFA）在G0細(xì)胞中顯著上調(diào)（圖2D），運(yùn)用IPA(Upon Ingenuity pathway analysis)分析，發(fā)現(xiàn)G0顯著富集于血管內(nèi)皮細(xì)胞的血管生成和粘附等信號(hào)通路（圖2E）。惡性細(xì)胞亞群G1特異性表達(dá)上皮-間充質(zhì)轉(zhuǎn)化（EMT）相關(guān)基因SRGN、VIM和LAPTM5A，說明惡性細(xì)胞亞群G1具有上皮間質(zhì)轉(zhuǎn)化的功能（圖2E）。惡性細(xì)胞亞群G2主要存在于LN3和PT患者中（圖2C）。作者在進(jìn)行IPA分析后發(fā)現(xiàn)G2細(xì)胞中的p53信號(hào)被激活，這將誘導(dǎo)細(xì)胞凋亡并抑制腫瘤生長(zhǎng)。因此，G2細(xì)胞可能在GC患者中的呈現(xiàn)休眠狀態(tài)。惡性細(xì)胞亞群G3細(xì)胞主要見于PT和M樣本中。在G3細(xì)胞中，腫瘤惡性相關(guān)基因RASD1、PIK3R1、JUN和FOS顯著上調(diào)（圖2D），SPINK1癌癥途徑和HGF信號(hào)通路被激活，表明G3促進(jìn)GC的惡性進(jìn)展。生長(zhǎng)因子途徑對(duì)EMT的調(diào)節(jié)也顯著增強(qiáng)（圖2D），且E M T相關(guān)基因（ZEB2、VIM和ID2）顯著上調(diào)，表明G3細(xì)胞中誘導(dǎo)了EMT進(jìn)展。同時(shí)作者發(fā)現(xiàn)與惡性上皮細(xì)胞的其他亞群相比，G3中CD44（圖2F）、PROM1（CD133）、LGR5等癌干細(xì)胞（CSCs）標(biāo)記物的表達(dá)也顯著上調(diào)，說明腫瘤干細(xì)胞在癌細(xì)胞增殖、遷移和轉(zhuǎn)移中發(fā)揮重要作用。
之后，作者運(yùn)用Monocle進(jìn)行細(xì)胞擬時(shí)序分析，可以看出惡性上皮細(xì)胞的動(dòng)態(tài)特征和異質(zhì)性（圖2G）。具體而言，在最初階段作者觀察到具有EMT誘導(dǎo)的CSC表型的G3細(xì)胞，隨后是G1細(xì)胞。休眠樣細(xì)胞G2和促血管生成G0細(xì)胞位于不同的軌跡分支，意味著惡性細(xì)胞進(jìn)入了不同的分化狀態(tài)。
最后作者進(jìn)一步使用來自TCGA胃腺癌數(shù)據(jù)集，進(jìn)行惡性細(xì)胞亞群特征基因的生存分析：我們發(fā)現(xiàn)高表達(dá)G0、G1、G3特征基因的患者預(yù)后較差(圖2H-J)。
圖2：胃癌中惡性上皮細(xì)胞的四個(gè)亞群及其與轉(zhuǎn)移相關(guān)的特征
Result3: 細(xì)胞T和B細(xì)胞在GC進(jìn)展過程中介導(dǎo)各種免疫反應(yīng)
作者通過細(xì)胞Marker基因注釋后共獲得24 448個(gè)T細(xì)胞，進(jìn)一步降維聚類獲得11個(gè)亞群。其中包括5個(gè)CD8+亞群、5個(gè)CD4+亞群和1個(gè)未知亞群(圖3A、3B)。
CD4+亞群由naive CD4+ T細(xì)胞，調(diào)節(jié)性T細(xì)胞(Tregs)，耗竭CD4+ T細(xì)胞，效應(yīng)記憶CD4+ T細(xì)胞(CD4+ TEM)和GADD45B+ Th1 CD4+ T細(xì)胞組成。CD8+亞群由naive CD8+ T細(xì)胞，細(xì)胞毒性CD8+ T細(xì)胞(GNLY, GZMB)，耗竭CD8+ T細(xì)胞，效應(yīng)記憶CD8+ T細(xì)胞和mucosal-associated invariant T細(xì)胞 (MAIT)組成（圖3A）。
如圖3B所示，腫瘤樣本(PT和M)中CD4+ T細(xì)胞和treg細(xì)胞的耗散程度較高，而M中CD4+ TEM、CD8+ TEM和細(xì)胞毒性CD8+ T細(xì)胞的比例與NT相比如預(yù)期的略有降低，表明腫瘤進(jìn)展和轉(zhuǎn)移過程中產(chǎn)生了免疫抑制情況。此外作者通過比較不同種類的樣本中的MAIT細(xì)胞比例，發(fā)現(xiàn)在轉(zhuǎn)移樣本（M）中，MAIT細(xì)胞不表達(dá)GZMB（編碼一種可以殺死感染細(xì)胞的蛋白質(zhì)）但高表達(dá)KLRG1（T細(xì)胞功能失調(diào)相關(guān)的基因）(圖3C)，與生物學(xué)先驗(yàn)知識(shí)背景相呼應(yīng)。免疫熒光染色顯示，在PT和M中均可觀察到MAIT細(xì)胞，在M中數(shù)量較多(圖3D)。
為了識(shí)別T細(xì)胞亞群的特征并分析他們之間的關(guān)聯(lián)，作者進(jìn)一步進(jìn)行了細(xì)胞通訊的分析，檢測(cè)細(xì)胞毒性/耗竭性CD8+ T細(xì)胞和NK細(xì)胞中KLRC1/KLRC2的表達(dá)水平：作者發(fā)現(xiàn)在腫瘤內(nèi)CD8+ T細(xì)胞和NK細(xì)胞中KLRC1表達(dá)均上調(diào)，同時(shí)細(xì)胞通訊分析還預(yù)測(cè)了每個(gè)患者不同細(xì)胞類型中HLA-E-KLRC1/KLRC2受配體的表達(dá)豐度 (圖3E)。作者進(jìn)一步將不同患者的CD8+ T細(xì)胞中的PDCD1表達(dá)與檢查點(diǎn)阻斷治療進(jìn)行關(guān)聯(lián)分析，以此推斷TIGIT和LAG3為患者1，2，5，6的共抑制受體(圖3F)。接著，作者對(duì)7708個(gè)B細(xì)胞進(jìn)行細(xì)胞亞群的注釋，共注釋為5個(gè)B細(xì)胞亞群（圖3G），不同樣本中B細(xì)胞類型的比例差異很大（圖3H），表明胃原發(fā)腫瘤和不同轉(zhuǎn)移灶中腫瘤微環(huán)境中的免疫異質(zhì)性問題。最后在本研究中，作者發(fā)現(xiàn)了一種細(xì)胞類型并命名為T細(xì)胞樣B細(xì)胞（T cell like-B cell），在M中表達(dá)最高，T cell like-B cell相關(guān)DEGs的通路分析顯示，白細(xì)胞/淋巴細(xì)胞的結(jié)合和NK細(xì)胞/淋巴細(xì)胞的細(xì)胞毒性顯著增強(qiáng)，說明T cell like-B cell與其他細(xì)胞亞型的相互作用(圖3I)。
圖3：GC進(jìn)展過程中由T和B細(xì)胞介導(dǎo)的各種免疫反應(yīng)。
Result4：內(nèi)皮細(xì)胞促進(jìn)血管生成并產(chǎn)生免疫抵抗
在本研究中，399個(gè)細(xì)胞被鑒定為內(nèi)皮細(xì)胞，其典型標(biāo)記基因表達(dá)一致。對(duì)其降維聚類后共分為4個(gè)亞型，分別為E0-E3(圖4A)。其中E0細(xì)胞在PT和M中特異性表達(dá)IGFBP5、STC1和IGFBP3，而在NT中不表達(dá)IGFBP3(圖4B)。促血管生成因子VEGFA和TGF-β是IGFBP5的上游調(diào)控因子。IPA顯示，與胃癌中腫瘤侵襲密切相關(guān)的mTOR和IGF-1信號(hào)在E0細(xì)胞中均被正調(diào)控，因此作者定義E0細(xì)胞亞群可能增加腫瘤細(xì)胞的侵襲和遷移(圖4C和D)。同時(shí)，IPA還顯示E1細(xì)胞與T細(xì)胞衰竭信號(hào)通路的調(diào)節(jié)密切相關(guān)(圖4D)。E2細(xì)胞在NT細(xì)胞中占比最高(88%)，其活性低于其他三個(gè)亞群，在本文章定義為正常內(nèi)皮細(xì)胞(圖4C和D)。E3細(xì)胞中白色脂肪組織褐變通路被激活，其中VEGF受體編碼基因NRP1和成纖維細(xì)胞生長(zhǎng)因子受體編碼基因FGFR1表達(dá)顯著上調(diào)(圖4D)，因此定義E3細(xì)胞亞群為參與血管生成的細(xì)胞亞群。最后作者又運(yùn)用TCGA胃癌的隊(duì)列數(shù)據(jù)進(jìn)行驗(yàn)證，進(jìn)一步說明高表達(dá)E1和E3細(xì)胞亞群特征的患者的預(yù)后較差(圖4F)。
圖4：內(nèi)皮細(xì)胞促進(jìn)GC血管生成并產(chǎn)生免疫抵抗
Result5：炎性癌成纖維細(xì)胞(iCAFs)與GC生長(zhǎng)和侵襲有關(guān)
作者將889個(gè)成纖維細(xì)胞無監(jiān)督聚類為兩個(gè)不同的亞型(圖5A-5C)。其中F0（高表達(dá)PDGFRA和CXCL12）是主要的成纖維細(xì)胞類型，在PT和M占比較高(圖5D)：高度表達(dá)RGS5和ACTA2，類似于膀胱尿路上皮癌中描述的iCAFs。F1細(xì)胞 (圖5D)被鑒定為肌癌相關(guān)成纖維細(xì)胞(mCAFs)。與F0類似，F(xiàn)1細(xì)胞在不同樣品中的比例差異較大(圖5B)。作者通過GO進(jìn)行功能富集分析，發(fā)現(xiàn)iCAF和mCAF均富集于細(xì)胞外基質(zhì)和細(xì)胞-底物粘附富集這一功能；但mCAFs特異性富集于肌肉相關(guān)過程；而iCAF則富集于細(xì)胞外基質(zhì)分解、白細(xì)胞遷移調(diào)節(jié)、以及細(xì)胞生長(zhǎng)和血管生成的調(diào)節(jié)(圖5E)。MMP2和MMP14被認(rèn)為促進(jìn)細(xì)胞外基質(zhì)降解和腫瘤細(xì)胞侵襲的基因，這兩個(gè)基因主要在iCAFs中表達(dá)，影響top富集的生物功能(圖5F)，并能促進(jìn)細(xì)胞外基質(zhì)降解和腫瘤細(xì)胞侵襲。CXCL12也是iCAFs中的特異性標(biāo)記物(圖5D)，其在CAFs中的高表達(dá)可促進(jìn)GC細(xì)胞侵襲和EMT。為了確認(rèn)iCAFs是否能調(diào)節(jié)細(xì)胞生長(zhǎng)和血管生成，作者進(jìn)一步分析了iCAFs和mCAFs中各種生長(zhǎng)因子編碼基因(VEGF、IGF、FGF家族)的表達(dá) (圖5G)，提示iCAFs可以促進(jìn)GC腫瘤生長(zhǎng)。最后通過細(xì)胞通訊分析(圖5I)，配體-受體對(duì)CXCL12-CXCR4參與了iCAFs與T、B或髓系細(xì)胞之間的相互作用，經(jīng)過文獻(xiàn)證實(shí)這是GC中免疫抵抗的一種機(jī)制。最后，作者將iCAF與生存進(jìn)行關(guān)聯(lián)分析，發(fā)現(xiàn)與低表達(dá)患者相比，高表達(dá)iCAF特征基因的患者患者的總體生存率更差（p<0.05），進(jìn)一步驗(yàn)證了iCAF與細(xì)胞生長(zhǎng)和血管生成相關(guān)這一結(jié)論。(圖5H)
圖5：CAF可在GC的TME中識(shí)別，并與腫瘤侵襲相關(guān)。
Result 6: GC中髓細(xì)胞譜系的多樣性
作者通過無監(jiān)督聚類和細(xì)胞marker基因的注釋，共確定了1801個(gè)單核吞噬細(xì)胞系統(tǒng)細(xì)胞（326個(gè)樹突狀細(xì)胞（DC）、894個(gè)巨噬細(xì)胞和581個(gè)單細(xì)胞）、139個(gè)肥大細(xì)胞和3579個(gè)中性粒細(xì)胞（圖6A）。作者進(jìn)一步對(duì)每種細(xì)胞亞群又進(jìn)行了細(xì)分，其中樹突狀細(xì)胞DC被重新分類為經(jīng)典DC1（cDC1）、經(jīng)典DC2（cDC2）、漿細(xì)胞樣DC（pDCs）和活化DC，經(jīng)過分類別的數(shù)量統(tǒng)計(jì)，cDC1和cDC2在PT和M中富集，而NT中活化DC的比例較大（圖6B）。有趣的是，在NT中很少發(fā)現(xiàn)pDC，但在LN和一些PT中出現(xiàn)，分析原因?yàn)閜DC高度表達(dá)顆粒酶基因GZMB和白細(xì)胞免疫球蛋白樣受體家族基因LILRA4和LILRB4，低表達(dá)CD86、CD80、CD83和LAMP3等基因，因此顯示出免疫抑制表型（圖6C）。
作者對(duì)Ma0-Ma3四個(gè)巨噬細(xì)胞亞群計(jì)算之前較為常見的M0和M1打分，進(jìn)一步對(duì)巨噬細(xì)胞進(jìn)行了詳細(xì)的劃分(圖6D和E)。Ma0聚集在PT中，高表達(dá)的誘導(dǎo)血管生成和癌細(xì)胞遷移的重要巨噬細(xì)胞基因VEGFA和SPARC(圖6F)，說明Ma0可能是一種與血管生成相關(guān)的巨噬細(xì)胞類型。Ma1細(xì)胞表達(dá)的DEGs包括HLA-DPB1和HLA-DPA1 (MHC-II)(圖6F)，它們與促炎表型相關(guān)。Ma2中FABP5高表達(dá)顯示其與免疫抑制密切相關(guān)。
單核細(xì)胞分為兩種已知類型，CD14+CD16（FCGR3A）? 經(jīng)典單核細(xì)胞亞型Mo1和Mo3和CD14+CD16+細(xì)胞亞型Mo0和Mo2(圖6G-I)。PT含有大量Mo0、Mo1和Mo3單核細(xì)胞。相反，M特別是O和Li在Mo2細(xì)胞中富集（圖6H）。IPA顯示Mo2細(xì)胞具有增強(qiáng)腫瘤細(xì)胞吞噬和細(xì)胞死亡的功能，因此，Mo2細(xì)胞可能對(duì)腫瘤細(xì)胞產(chǎn)生細(xì)胞毒性作用。
中性粒細(xì)胞亞簇顯示六個(gè)簇（圖6J和K）。其中LYZ編碼的溶菌酶是Paneth細(xì)胞分泌的一種抗菌肽，與自噬和維持免疫穩(wěn)態(tài)有關(guān)。作者同時(shí)進(jìn)行了IPA分析，富集結(jié)果顯示N0細(xì)胞吞噬功能增強(qiáng)（圖6L），并進(jìn)行實(shí)驗(yàn)驗(yàn)證，通過RNAscope和標(biāo)記基因LYZ和NAMPT進(jìn)一步證實(shí)了GC中存在Paneth樣細(xì)胞（圖6M）。同時(shí)作者對(duì)髓源性抑制細(xì)胞（MDSC）進(jìn)行了進(jìn)一步亞型的劃分，進(jìn)一步刻畫了胃癌患者腫瘤異質(zhì)性的問題。
圖6：GC進(jìn)展期髓樣細(xì)胞的富集情況
Result 7:淋巴結(jié)源性耗竭性CD8+T細(xì)胞的20個(gè)基因特征被證實(shí)可以預(yù)測(cè)淋巴結(jié)轉(zhuǎn)移
腫瘤浸潤(rùn)性T細(xì)胞和B淋巴細(xì)胞在癌癥發(fā)展和轉(zhuǎn)移中起著關(guān)鍵作用；了解這些免疫細(xì)胞在不同器官GC轉(zhuǎn)移中的浸潤(rùn)狀態(tài)，可以為器官特異性轉(zhuǎn)移制定特定的診斷和治療策略。作者采用無監(jiān)督方法對(duì)T細(xì)胞和B淋巴細(xì)胞的每個(gè)亞群進(jìn)行重新聚類，探討GC中淋巴細(xì)胞的器官特異性轉(zhuǎn)移特征。如圖7A TSNE降維結(jié)果所示，在大多數(shù)T細(xì)胞亞群和一些B細(xì)胞亞群中，不同的器官轉(zhuǎn)移在轉(zhuǎn)錄上是異質(zhì)的（每一種疾病樣本類型中包含不同的亞型）。
早期GC中，淋巴結(jié)轉(zhuǎn)移的發(fā)生率極高，在單細(xì)胞分辨率下識(shí)別淋巴結(jié)轉(zhuǎn)移的特征具有重要意義。因此，作者通過差異分析，識(shí)別了20個(gè)（top 20）GC中生成淋巴結(jié)衍生（lymph node-derived）T細(xì)胞和B細(xì)胞亞群的基因表達(dá)特征并為每個(gè)淋巴結(jié)衍生的T細(xì)胞和B細(xì)胞亞群添加了20個(gè)基因標(biāo)簽。如圖7B所示，在耗竭性CD8+T細(xì)胞中，淋巴結(jié)衍生亞群中前20個(gè)上調(diào)的DEG在LN中的表達(dá)高于其他轉(zhuǎn)移瘤，同時(shí)與其他原發(fā)腫瘤相比，PT3（淋巴結(jié)轉(zhuǎn)移患者的原發(fā)腫瘤樣本）中的表達(dá)也更高，表明20個(gè)基因特征可以反映從原發(fā)部位向淋巴結(jié)轉(zhuǎn)移的模式；作者進(jìn)一步將20個(gè)signature與患者的預(yù)后進(jìn)行關(guān)聯(lián)，說明了這20個(gè)signature對(duì)于患者生存時(shí)間的影響（圖7D）。
圖7：不同器官轉(zhuǎn)移的T和B淋巴細(xì)胞亞群異質(zhì)性探究
自此這篇文章的主題內(nèi)容就介紹完畢啦，總體來說，作者對(duì)胃癌特異性器官轉(zhuǎn)移的6名患者進(jìn)行了單細(xì)胞RNA-seq測(cè)序，通過對(duì)腫瘤微環(huán)境中每種細(xì)胞類型的刻畫，識(shí)別了胃癌患者中，原發(fā)腫瘤和器官轉(zhuǎn)移的異質(zhì)性問題，結(jié)合TCGA組學(xué)數(shù)據(jù)進(jìn)一步進(jìn)行了驗(yàn)證；最終識(shí)別到了可以作為胃癌特異性性器官轉(zhuǎn)移的腫瘤浸潤(rùn)性T細(xì)胞和B淋巴細(xì)胞的signatures。整體的文章思路清晰，尤其是在GC腫瘤微環(huán)境中免疫細(xì)胞的刻畫上細(xì)致入微，可見，數(shù)據(jù)量并不是發(fā)好文章的唯一法寶，與生物學(xué)意義緊密結(jié)合，也同樣可以盡顯工作的意義。
