hello，大家好瑰钮，我們的分析也慢慢進(jìn)入深水區(qū)了冒滩，很多方法需要聯(lián)合起來使用，如何靈活運用各種方法浪谴，就需要我們的經(jīng)驗和智慧了开睡。今天分享的文獻(xiàn)在Multi-regional characterisation of renal cell carcinoma and microenvironment at single cell resolution，幾乎用到了我們常見的所有的單細(xì)胞的分析思路和方法苟耻，當(dāng)然了篇恒，生物學(xué)意義也進(jìn)行了很好的詮釋，尤其是利用單細(xì)胞數(shù)據(jù)calling Mutation凶杖，大家多多學(xué)習(xí)和積累胁艰。
放一張我覺的最關(guān)鍵的圖
1.png
Abstract
腫瘤行為取決于癌細(xì)胞的致癌特性及其多細(xì)胞相互作用。 通過從 12 名腎腫瘤患者獲得的 270,000 個單細(xì)胞轉(zhuǎn)錄組和 100 個微解剖的完整外顯子組檢查這些依賴性(現(xiàn)在的單細(xì)胞數(shù)據(jù)都是以大數(shù)據(jù)量著稱)官卡。 從腫瘤核心蝗茁、腫瘤-正常界面、正常周圍組織和外周血的多個區(qū)域取樣組織(這個時候空間轉(zhuǎn)錄組更具有說服力)寻咒。 發(fā)現(xiàn) CD8+ T 細(xì)胞克隆型的主要空間位置在很大程度上定義了衰竭狀態(tài)哮翘，而體細(xì)胞瘤內(nèi)異質(zhì)性無法解釋克隆型異質(zhì)性。 從單細(xì)胞 RNA 測序數(shù)據(jù)中calling從頭突變使我們能夠譜系追蹤和推斷細(xì)胞的克隆性毛秘。 發(fā)現(xiàn)了六個區(qū)分腫瘤細(xì)胞功能的meta程序饭寺。 在腫瘤-正常界面富集的上皮-間質(zhì)轉(zhuǎn)化meta程序似乎通過巨噬細(xì)胞表達(dá)的 IL1B 進(jìn)行調(diào)節(jié)，可能形成治療靶點叫挟。
Introduction
Clear cell renal cell carcinoma (ccRCC) 是腎細(xì)胞癌 (RCC) 最常見的亞型艰匙，約占 RCC 病例的 75% 和the majority of deaths from kidney cancer。許多努力已經(jīng)表征了 ccRCC 的基因組landscape抹恳，揭示了重要的驅(qū)動事件员凝，例如 VHL 的雙等位基因失活（最常見的是伴隨著染色體 3p 的丟失和 VHL 中的突變/表觀遺傳沉默），隨后是染色質(zhì)重塑和組蛋白修飾相關(guān)的突變基因 PBRM1奋献、BAP1 和 SETD2健霹。已經(jīng)系統(tǒng)地研究了 ccRCC 演變過程中啟動事件的時間旺上。正如之前的多區(qū)域外顯子組測序研究所揭示的那樣，隨后突變事件的腫瘤內(nèi)異質(zhì)性 (ITH) 似乎是 ccRCC 的一個顯著特征糖埋。相比之下宣吱，ccRCC 在轉(zhuǎn)錄水平上的 ITH 不太清楚，部分原因是包含腫瘤微環(huán)境 (TME) 的多細(xì)胞生態(tài)系統(tǒng)的復(fù)雜性瞳别。特別是征候，ccRCC TME 中惡性和非惡性細(xì)胞的表型異質(zhì)性及其與空間定位的關(guān)系仍然不清楚。
ccRCC 是一種免疫細(xì)胞大量浸潤的癌癥類型祟敛。免疫檢查點阻斷 (ICB) 治療已被證明可有效提高患者的生存率疤坝，突出了探索 ccRCC 免疫微環(huán)境的重要性。以前的研究使用批量測序分析了 ccRCC 的免疫狀況垒棋，這限制了它們剖析不同免疫細(xì)胞群的能力卒煞。使用質(zhì)譜流式細(xì)胞術(shù)的 ccRCC 綜合免疫圖譜揭示了 ccRCC 腫瘤生態(tài)系統(tǒng)中的免疫細(xì)胞多樣性痪宰。單細(xì)胞 RNA 測序 (scRNA-seq) 的最新進(jìn)展及其在癌癥研究中的應(yīng)用徹底改變了我們對腫瘤細(xì)胞表型異質(zhì)性叼架、腫瘤免疫landscape、TME 的復(fù)雜性和可塑性以及 TME 中細(xì)胞間通訊的理解衣撬。具體而言乖订，在 ccRCC 中，最近的一項 scRNA-seq 研究提供了支持其起源于近端腎小管細(xì)胞的證據(jù)具练。其他研究利用 scRNA-seq 研究 ccRCC 的免疫狀況乍构，主要關(guān)注 ICB 治療相關(guān)隊列和不同疾病階段，揭示與治療效果或疾病進(jìn)展相關(guān)的關(guān)鍵特征扛点。
在考慮 TME 的異質(zhì)性時哥遮，感興趣的空間區(qū)域擴展了與突變 ITH 相關(guān)的區(qū)域。更廣泛的目標(biāo)區(qū)域包括循環(huán)血液陵究、腫瘤-正常界面或腫瘤假包膜（代表腫瘤與相鄰正常腎臟之間的邊界）眠饮、相鄰正常腎臟和腎周脂肪組織。假包膜的纖維結(jié)締組織似乎在空間上限制了生長铜邮，侵襲與腫瘤分期和分級相關(guān)仪召。由于 RCC 中的obesity paradox，腎周肥胖引起了人們的興趣松蒜，其中肥胖是診斷腎癌的最強風(fēng)險因素之一扔茅，但也與腫瘤學(xué)結(jié)果的改善有關(guān)。了解 ccRCC 在腫瘤細(xì)胞秸苗、各種免疫/基質(zhì)細(xì)胞方面的空間異質(zhì)性和進(jìn)化召娜，以及它們在更廣泛的 TME 中的相互作用仍然缺乏。為了解決這個問題惊楼，對 12 名患者進(jìn)行了基于多區(qū)域的 scRNA-seq玖瘸，對外周血吐句、正常腎臟、腫瘤核心的四個不同空間區(qū)域和腫瘤-正常界面進(jìn)行了采樣店读，同時對激光捕獲顯微切割進(jìn)行了局部詳盡的外顯子組測序(LCM) 衍生的腫瘤樣本嗦枢。
Results
Multi-region based genomic and single-cell transcriptomic profiling of RCC
對 12 名接受放射學(xué)診斷腎腫瘤手術(shù)切除的患者進(jìn)行了多區(qū)域基因組和單細(xì)胞轉(zhuǎn)錄組分析。 在組織病理學(xué)檢查后屯断，12 名患者中有 10 名的腫瘤被評估為 ccRCC文虏，1 名（PD47172）為嗜酸細(xì)胞瘤，1 名（PD44714）為大的良性厚壁囊腫殖演。 在每位患者中氧秘，從外周血、正常腎臟趴久、腫瘤核心的四個不同空間區(qū)域以及腫瘤-正常界面采集組織丸相。 此外，如果有的話彼棍，還從腎周脂肪灭忠、正常腎上腺、腎上腺轉(zhuǎn)移和腫瘤血栓中取樣座硕。
2.png
在可以從這些樣本中提取足夠數(shù)量的活單細(xì)胞的情況下弛作，使用 10x 平臺進(jìn)行了基于液滴的 5' scRNA-seq 和 T 細(xì)胞受體 (TCR) 富集。 同時华匾，在執(zhí)行全外顯子組測序 (WES) 之前映琳，使用 LCM 從每個含有腫瘤組織的區(qū)域解剖每個患者的微活檢樣本。 LCM 方法使我們能夠通過詢問亞毫米大小的活檢來探索 ITH 的極限蜘拉。 基于 WES 數(shù)據(jù)萨西，確定了在 ccRCC 中報告為復(fù)發(fā)/驅(qū)動事件的基因組改變。 9 名 ccRCC 患者中有 7 名（一名 ccRCC 患者沒有數(shù)據(jù)）具有 VHL 突變旭旭，4 名具有 PBRM1 突變谎脯，3 名攜帶 BAP1 突變。 在所有 9 名患者中都檢測到染色體 3p 的拷貝數(shù)丟失您机。 嗜酸細(xì)胞瘤攜帶整個 1 號染色體的特征性拷貝數(shù)丟失穿肄。
使用 scRNA-seq，在嚴(yán)格的質(zhì)量控制后從大約 270,000 個細(xì)胞中捕獲轉(zhuǎn)錄組际看，根據(jù)典型標(biāo)記基因的表達(dá)可以將其大致分為 12 種主要細(xì)胞類型
3.png
由于合適的單細(xì)胞分離方案咸产，T 細(xì)胞（表達(dá) CD3D、CD3G 和 CD3E）在數(shù)據(jù)中最為豐富仲闽，其次是骨髓細(xì)胞（表達(dá) CD68 和 CD14）和自然殺傷 (NK) 細(xì)胞（表達(dá) NCR1 和 KLRC1 )脑溢。還捕獲了大量其他免疫細(xì)胞群：B 細(xì)胞（表達(dá) CD79A 和 MS4A1）、漿細(xì)胞（表達(dá) IGKC 和 IGHG1）、肥大細(xì)胞（表達(dá) TPSAB1 和 TPSB2）和漿細(xì)胞樣樹突細(xì)胞（pDC屑彻；表達(dá) IRF7 和 LILRA4） .除了免疫區(qū)室验庙，數(shù)據(jù)中還鑒定了四種基質(zhì)細(xì)胞類型：內(nèi)皮細(xì)胞 (EC)（表達(dá) PLVAP 和 TIMP3）、成纖維細(xì)胞（表達(dá) ACTA2 和 TAGLN）社牲、近端腎小管 (PT) 細(xì)胞（表達(dá) SLC22A8 和 GPX3）和非 PT 上皮細(xì)胞（表達(dá) DEFB1 和 UMOD）粪薛。從 scRNAseq 數(shù)據(jù)推斷，RCC 腫瘤細(xì)胞在特異性表達(dá) CA9 并在其基因組中具有廣泛拷貝數(shù)變異 (CNV) 的cluster內(nèi)被鑒定出來搏恤。接下來违寿，我們研究了 12 種主要細(xì)胞類型的起源組織，并觀察了這些細(xì)胞類型的不同組織分布（分析中合并了四個區(qū)域）
%1S6}NI9`AYD_JVVL}%I0D8.png
image.png
- 注：(A) Heatmap illustrating the clinical features (left panel) and the genomic landscape (middle panel) and copy number profiles (right panel) of the tumours sequenced. (B) UMAP, (C) cell type, and (D) sampled region depicting the interpatient variability of scRNA-seq results. (E) UMAP showing marker gene expression for all cells. (F) Copy number inference based on scRNA-seq data.
進(jìn)一步研究中涵蓋的主要細(xì)胞區(qū)室進(jìn)行了亞群聚類分析熟空。 NK 細(xì)胞的亞群產(chǎn)生了 14 個具有差異表達(dá)基因 (DEG) 和異質(zhì)組織來源的cluster藤巢。在這些cluster中，確定了眾所周知的 CD56 (NCAM1) 和 CD16 (FCGR3A) 表達(dá)群體息罗。先天淋巴細(xì)胞 (ILC) cluster的特征在于 IL7R 和 FXYD7 的表達(dá)掂咒。兩個 NK cluster（cluster 2 和 6）顯示干擾素γ（IFNG）的高表達(dá)，cluster 6 也高表達(dá)細(xì)胞因子 CCL4L2迈喉，并可能在正常腎上腺中富集绍刮。確定了一些特征較少的 NK cluster，例如cluster 4弊添，它特異性表達(dá) KRT81 和 KRT86录淡。這種 NK 細(xì)胞亞群之前曾在肝細(xì)胞癌中報道過捌木，但其功能尚不清楚油坝。 B/漿細(xì)胞區(qū)室被分為 13 個cluster，其中確定了眾所周知的主要 B 細(xì)胞群刨裆，如幼稚澈圈、轉(zhuǎn)換記憶和非轉(zhuǎn)換記憶 B 細(xì)胞》校活化的 B 細(xì)胞（AREGhigh 和 RHOBhigh cluster）可能在腫瘤中富集瞬女，而 AREGhigh cluster在腎周脂肪中更富集。與主要細(xì)胞類型的外周血組織分布相比努潘，發(fā)現(xiàn)血漿 IgA诽偷、IgG 和循環(huán)細(xì)胞（分別表達(dá) IGHA1、IGHG1 和 MKI67）在組織中富集疯坤。
image.png
- 注：(A&E) UMAPs, (B&F) dotplots of marker genes, (C&G) tissue distribution, and (D&H) regional enrichment of NK and B cells.
CD8+ T cell lineages reveal evolutionary trajectories and the influence of spatial location on exhaustion
T 細(xì)胞包含數(shù)據(jù)中最豐富的細(xì)胞類型报慕，大致分為兩個主要部分：CD4+ 和 CD8+ T 細(xì)胞（包括 gamma delta T 細(xì)胞）。 CD8+ T 細(xì)胞區(qū)室的亞聚類導(dǎo)致鑒定出 18 個具有不同 DEG 和異質(zhì)組織位置的cluster压怠。
image.png
總的來說眠冈，根據(jù)典型標(biāo)記基因的表達(dá)確定了典型的 CD8+ T 細(xì)胞cluster，它們代表了不同的 T 細(xì)胞功能狀態(tài)菌瘫，包括幼稚蜗顽、效應(yīng)布卡、記憶、前功能障礙和功能障礙雇盖。發(fā)現(xiàn)高表達(dá) LEF1 和 CCR7 等基因的幼稚/中樞記憶 (CM) CD8+ T 細(xì)胞主要在外周血中富集忿等。鑒定了常駐記憶 (RM) T 細(xì)胞，因為它們表達(dá)組織駐留標(biāo)志物（即 ITGAE 和 CD69）崔挖，并且主要在正常腎臟中富集这弧。特別是，cluster 10 還高表達(dá) CXCL13虚汛，它可能在 B 細(xì)胞募集和三級淋巴結(jié)構(gòu)的形成中發(fā)揮潛在作用匾浪。我們發(fā)現(xiàn)cluster 6 高度表達(dá) FGFBP2 和 CX3CR1，并且在外周血中大量富集卷哩，因此該cluster可能代表最近激活的效應(yīng)記憶 T 細(xì)胞（CD8+ T_EMRA）蛋辈。基于包括 LAG3将谊、TIGIT冷溶、PDCD1、HAVCR2 和 CTLA4 在內(nèi)的基因表達(dá)升高尊浓，鑒定了兩個耗盡的 T 細(xì)胞cluster（cluster 7 和 8）逞频。有趣的是，我們發(fā)現(xiàn)cluster 8 具有最高的 LAG3 表達(dá)栋齿，并特異性表達(dá)免疫抑制細(xì)胞因子 IL10苗胀。該cluster可能代表具有極高效應(yīng)子和功能障礙水平的 CD8+ T 細(xì)胞，它們通過產(chǎn)生 IL-10 發(fā)揮調(diào)節(jié)功能瓦堵。發(fā)現(xiàn)粘膜相關(guān)不變性 T (MAIT) 細(xì)胞高表達(dá) TRAV1-2 和 IL7R基协。鑒定出兩個循環(huán)細(xì)胞cluster：一個（表達(dá)MCM5和PCNA）代表細(xì)胞周期G1/S期的細(xì)胞，另一個（表達(dá)TOP2A和MKI67）代表G2/M期的細(xì)胞菇用。除了傳統(tǒng)的 CD8+ T 細(xì)胞cluster澜驮，我們還鑒定了兩個 gamma delta T 細(xì)胞cluster：gdT_Vd1（表達(dá) TRDV1）和 gdT_Vd2（表達(dá) TRDV2）。我們還對 CD4+ T 細(xì)胞群進(jìn)行了亞群聚類分析惋鸥，揭示了各種亞型杂穷，例如 CD4+ na?ve/CM 和 CD4+ 調(diào)節(jié)性 T 細(xì)胞 (Tregs) 及其不同的組織分布.
image.png
- 注：(A) UMAP of the subclustering result, (B) UMAP and (C) barplot showing tissue distribution, and (D) regional enrichment of CD8+ T cell population. (E) UMAP showing the RNA velocity result and (F) exhaustion score across the pseudotime trajectory of CD8+ T cells. (G) UMAP of the sub-clustering result, (H) dotplot showing marker gene expression, and (I) UMAP of tissue distribution of CD4+ T cells. (J) Barplot showing the comparison of TCR clonal expansion between CD8+ and CD4+ T cells, breaking down into three different locations: blood, tumour and other regions
接下來，我們使用 Monocle 3 和 RNA 速度分析對 CD8+ T 細(xì)胞進(jìn)行了偽時間軌跡分析卦绣，不包括 gamma delta T 和循環(huán)cluster耐量。沿著假時間軌跡，我們發(fā)現(xiàn)細(xì)胞毒性相關(guān)基因（即 KLRG1迎卤、GNLY 和 GZMH）是逐漸下調(diào)而功能障礙相關(guān)基因（即 CTLA4拴鸵、HAVCR2 和 LAG3）逐漸上調(diào)。典型的 T 細(xì)胞功能障礙前相關(guān)基因（即 CXCR4、GZMK 和 GZMA) 最初被上調(diào)劲藐，然后沿著偽時間軌跡下降八堡。因此，這個偽時間軌跡概括了 CD8+ T 細(xì)胞從細(xì)胞毒性狀態(tài)經(jīng)過功能障礙前狀態(tài)到功能障礙狀態(tài)的進(jìn)展聘芜，同時耗竭程度逐漸升高兄渺。假時間和耗竭評分之間的正相關(guān)也支持這種進(jìn)展。此外汰现，將前 10 個擴展的 TCR 克隆型投影到軌跡上導(dǎo)致觀察到單個 TCR 譜系通常僅限于類似的表型狀態(tài)挂谍，而不是分布在整個軌跡上。在所有腫瘤中，發(fā)現(xiàn) 90% 的具有 23 個或更多細(xì)胞的克隆型被限制在假時間值范圍內(nèi)（Wilcoxon 檢驗，p < 0.05）晚唇。在多位患者中觀察到每個克隆超過 100 個細(xì)胞的高度擴增的 TCR 克隆，其中顯著高達(dá) 30% 的 CD8+ T 細(xì)胞可以來自單個克隆型妄田。相比之下，與 CD8+ 群體中的 TCR 克隆型相比驮捍，CD4+ 群體中的 TCR 克隆型擴展較少疟呐。許多擴增最多的 CD8+ TCR 克隆具有相當(dāng)比例的循環(huán)細(xì)胞，但在較少消耗的克隆型中觀察到的情況除外东且。這一發(fā)現(xiàn)表明 RCC 中高度耗竭的 T 細(xì)胞的增殖并未完全停止启具，類似于先前在黑色素瘤中的發(fā)現(xiàn)。
檢查了血液中檢測到的 TCR 克隆型是否反映了在其他區(qū)域檢測到的那些珊泳。研究發(fā)現(xiàn)鲁冯，無論克隆大小如何，平均耗竭程度（推斷的假時間）和在外周血中檢測到 CD8+ TCR 克隆的概率都具有很強的反相關(guān)性旨椒，以至于在血液中很少檢測到耗竭的克隆型晓褪。這一發(fā)現(xiàn)出乎意料，表明組織駐留的耗竭 CD8+ T 細(xì)胞克隆似乎不會在外周血中再循環(huán)综慎。為了進(jìn)一步說明 T 細(xì)胞耗竭、克隆擴增及其組織分布之間的關(guān)系勤庐，我們根據(jù) CD8+ T 細(xì)胞是單細(xì)胞還是擴增以及它們的主要組織位置（血液示惊、正常組織或腫瘤）對它們進(jìn)行分類。腫瘤中擴增的 T 細(xì)胞進(jìn)一步細(xì)分為出現(xiàn)在所有腫瘤區(qū)域的那些和不出現(xiàn)的（腫瘤同質(zhì)和異質(zhì)）愉镰。值得注意的是米罚，CD8+ T 細(xì)胞的表型狀態(tài)，就耗竭程度而言丈探，表現(xiàn)出對克隆擴增和組織位置的強烈依賴录择。同時，clones that were private to one tumour region were not significantly more exhausted than those shared between different regions 。
image.png
- 注：(A) Dotplot showing marker gene expression defines principal CD8+ cell types. EM, effector memory; Act, activated; EFF, effector; EX, exhausted. (B) UMAP depicting the pseudotime inference of CD8+ cells隘竭。 (C) Expression of canonical exhaustion markers across cells ordered by pseudotime analysis. All marker genes are statistically significant across pseudotime values (q value = 0). (D) UMAP showing the 10 most expanded clones from patient PD43948. Grey dots represent cells outside the 10 most expanded clones. (E) Boxplot depicting the most expanded clonotypes (>100 cells) across all patients, ordered by their mean pseudotime values, showing the median, interquartile range and outlier pseudotime values (top panel); barplots showing the maximum expansion for the most expanded region (middle panel) and percentage of cycling cells (lower panel). (F) The probability of detecting a given TCR clone in peripheral blood as a function of minimal clone size and mean pseudotime value of the clone. (G) Mean pseudotime values based on the categorisation of clonotypes according to their principal region of enrichment.
Spatial localisation rather than intra-tumoural heterogeneity primarily influences CD8+ clonotypic heterogeneity
使用從 WES 數(shù)據(jù)中calling的體細(xì)胞突變塘秦，構(gòu)建了系統(tǒng)發(fā)育樹來闡明我們研究中腫瘤中的克隆進(jìn)化和 ITH。 總體而言动看，我們發(fā)現(xiàn)所有腫瘤克隆共享一個長主干但有短分支
image.png
- 注：Reconstructed phylogenies from WES of multi-regional LCM biopsies. Each node represents a mutant clone present in one or more of the biopsies.
這表明大多數(shù)體細(xì)胞突變在個體腫瘤中普遍存在尊剔，只有少數(shù)個體突變被檢測到。大多數(shù)檢測到的驅(qū)動突變和關(guān)鍵 CNV（即 VHL 突變和染色體 3p 雜合性丟失 (LOH)）由個體腫瘤內(nèi)的所有腫瘤克隆共享菱皆，因此位于系統(tǒng)發(fā)育樹的樹干上须误。此外，根據(jù)變異等位基因頻率分布仇轻，我們測序的絕大多數(shù) LCM 樣本都出現(xiàn)了克隆性京痢。綜上所述，WES 顯示我們隊列中腫瘤的 ITH 范圍是有限的篷店。先前的研究通過比較在來自腫瘤不同空間定位的樣本中檢測到的體細(xì)胞突變历造，廣泛研究了各種癌癥的腫瘤內(nèi)遺傳異質(zhì)性。然而船庇，體細(xì)胞異質(zhì)性對不同空間定位的局部腫瘤微環(huán)境的影響吭产，尤其是抗腫瘤免疫反應(yīng)，在很大程度上仍未得到表征鸭轮。在這里臣淤，我們系統(tǒng)地比較了個體腫瘤中 CD8+ T 細(xì)胞的體細(xì)胞突變、空間定位和 TCR 克隆型之間的關(guān)系窃爷。
image.png
- 注：Comparison of WES derived phylogenies (left) with geographic location (centre) and CD8+ TCR clonotype expansion (right). Colours reference somatic clones to spatial localisation. Each column in the right panel represents a TCR clonotype, those with significant regional enrichment are highlighted in red. a, c, d, and e represents four different regions of the tumour core; g, tumour-normal interface; n, normal kidney; b, peripheral blood; h, normal adrenal gland.
出乎意料的是邑蒋，發(fā)現(xiàn) T 細(xì)胞克隆型在不同的空間定位中具有高度異質(zhì)性和差異性，即使在僅觀察到體細(xì)胞突變異質(zhì)性可忽略不計的腫瘤區(qū)域也是如此按厘。 在腫瘤細(xì)胞上產(chǎn)生新抗原的體細(xì)胞突變被認(rèn)為是抗原呈遞時 T 細(xì)胞克隆擴增的驅(qū)動因素医吊。 我們的發(fā)現(xiàn)表明 TCR 克隆擴增的異質(zhì)性更多地與 T 細(xì)胞在組織中的不同空間定位相關(guān)，而不是與體細(xì)胞突變的 ITH 相關(guān)逮京。 為了正式檢驗這一點卿堂，我們計算了 T 細(xì)胞克隆型距離與 1) 突變距離之間的相關(guān)性； 2）空間定位距離
image.png
- 注：Scatter plot of the Mantel correlation between tree distances. x-axis represents the correlation coefficient between WES derived clones and TCR clonotype distances. y-axis represents the correlation coefficient between spatial localisation and TCR clonotype distances.懒棉、
通過比較這兩種相關(guān)性草描，我們發(fā)現(xiàn) CD8+ T 細(xì)胞中的 TCR 異質(zhì)性與空間定位而非體細(xì)胞異質(zhì)性的相關(guān)性更強（配對 Wilcoxon 檢驗，p < 0.05）策严。
Precise de novo somatic mutation calling from scRNA-seq data
從轉(zhuǎn)錄組序列檢測單個細(xì)胞內(nèi)的體細(xì)胞突變可能有助于推斷它們的克隆關(guān)系穗慕。 盡管理論上可以從 scRNAseq 數(shù)據(jù)中calling突變，但目前沒有高精度的方法可用妻导。 在這里逛绵，我們開發(fā)了一種算法/pipeline來從 scRNA-seq 數(shù)據(jù)執(zhí)行從頭體細(xì)胞突變calling怀各。 簡而言之，我們在初始過濾步驟之前使用 bcftools 從單細(xì)胞 BAM 文件中calling突變术浪，以去除單個細(xì)胞中存在的突變和不同細(xì)胞譜系之間共享的突變瓢对。 在應(yīng)用二項式過濾器之前，我們檢查了初始步驟中calling的所有位點的參考和變異等位基因計數(shù)添吗。 然后沥曹，我們在對變體進(jìn)行注釋后應(yīng)用了一組最終的過濾指標(biāo)
image.png
- 注：Schematic of the de novo mutation calling algorithm.
為了對我們的突變calling方法進(jìn)行基準(zhǔn)測試，我們首先比較了從腫瘤細(xì)胞的 scRNA-seq 數(shù)據(jù)calling的體細(xì)胞突變與從腫瘤 WES 數(shù)據(jù)calling的那些突變碟联。所有檢測到的突變都被歸類為真陽性（由突變calling檢測到妓美，或未被突變calling者calling但在原始數(shù)據(jù)中有足夠的支持證據(jù)）、假陽性鲤孵、假陰性或不確定（其中突變在沒有足夠的 WES 覆蓋來驗證呼叫）壶栋。總體而言普监，我們的方法取得了良好的性能贵试，精度為 0.64（僅考慮外顯子突變時為 0.70），靈敏度為 0.53凯正。我們還能夠在 CD8+ T 細(xì)胞中對該方法進(jìn)行基準(zhǔn)測試毙玻，結(jié)果表明 84% 的被calling突變僅限于單個 TCR 克隆。這證實了預(yù)期的發(fā)現(xiàn)廊散，即在 CD8+ T 細(xì)胞中calling的大多數(shù)突變僅限于克隆型桑滩，因為在胸腺成熟前 T 細(xì)胞克隆之間可以共享的突變數(shù)量非常有限。
image.png
- 注：Benchmarking results for scRNA-seq derived calls against WES data for each patient donor
使用這些突變calling允睹，我們研究了不同細(xì)胞類型表達(dá)的突變數(shù)量运准，這可能有助于揭示它們的克隆擴增程度。我們計算了具有一個缭受、兩個胁澳、三個或三個以上突變的細(xì)胞的比例。我們需要來自每個細(xì)胞譜系和患者的至少 100 個細(xì)胞來解釋稀有細(xì)胞群中缺乏辨別力的問題米者。正如預(yù)期的那樣韭畸，表達(dá)所謂突變的細(xì)胞數(shù)量最多的譜系是腫瘤細(xì)胞，這主要由譜系的已知克隆結(jié)構(gòu)來解釋塘雳，但也由于與正常細(xì)胞類型相比突變負(fù)擔(dān)增加的可能性陆盘。出于類似的原因，基質(zhì)細(xì)胞通常不會有可辨別數(shù)量的細(xì)胞具有不止一種被稱為突變的細(xì)胞败明。然而，我們觀察到大量表達(dá)突變的骨髓細(xì)胞太防，表明這些細(xì)胞中有相當(dāng)大的比例是無性系相關(guān)的妻顶。隨后是成纖維細(xì)胞和 CD8+ T 細(xì)胞（我們知道這些細(xì)胞是根據(jù) TCR 測序結(jié)果進(jìn)行克隆擴增的）酸员。極少比例的 CD4+ T 細(xì)胞表達(dá)突變，與基于 TCR 分析的低克隆度一致讳嘱。
image.png
- 注：Somatic mutation analysis and the relationship with TCR heterogeneity. (A) Density plots showing the VAF distribution of mutations called from WES data from LCM biopsies of tumour samples. (B) Comparison of WES derived phylogenies (left) with geographic location (centre) and CD8+ TCR clonotype expansion (right) for each patient with matching data. Colours reference somatic clones to spatial localisation. Each column in the right panel represents a TCR clonotype, those with significant regional enrichment are highlighted in red. a, c, d, and e represents four different regions of the tumour core; g, tumour-normal interface; f, perinephric fat; n, normal kidney; b, peripheral blood; h, normal adrenal gland; i, adrenal metastasis; t, thrombus. (C) Bar chart showing benchmarking results of the number of called mutations in CD8+ T cells derived from scRNA-seq data that are restricted to individual TCR clonotypes. (D) Comparison of the proportion of cells with one, two, three, and more than three mutations across the major cell types.
Regional characterisation and evolution of myeloid populations
Overall, we captured transcriptomes from 50,603 myeloid cells in our dataset, which were categorised into 19 clusters in sub-clustering analysis
image.png
我們首先根據(jù)典型標(biāo)志物的表達(dá)和細(xì)胞的組織起源對這些cluster進(jìn)行了廣泛的注釋
image.png
cluster 1幔嗦、2、3 和 4 主要存在于血液中沥潭，CD14 高表達(dá)但缺乏 FCGR3A 表達(dá)邀泉，因此代表循環(huán)經(jīng)典單核細(xì)胞。 cluster 5 代表循環(huán)非經(jīng)典單核細(xì)胞钝鸽，具有 FCGR3A 高表達(dá)但缺乏 CD14 表達(dá)汇恤。 我們確定了三個樹突狀細(xì)胞 (DC) cluster：pDC、1 型和 2 型常規(guī) DC（cDC1 和 cDC2）拔恰，分別以 JCHAIN因谎、CLEC9A 和 CD1C 的特異性表達(dá)為特征。 與其他區(qū)域相比颜懊，cDC1 在三個 DC cluster中顯示出在腫瘤核心中的富集财岔。 我們發(fā)現(xiàn)以 TPSAB1 特異表達(dá)為特征的肥大細(xì)胞可能在腫瘤核心中富集，這與之前的報道一致河爹。 值得注意的是匠璧，我們根據(jù) CD163 和 C1QC 的高表達(dá)確定了 9 個巨噬細(xì)胞cluster（cluster 6-8、11-16）咸这，反映了 RCC 中巨噬細(xì)胞群的顯著異質(zhì)性夷恍。
image.png
- 注：Sub-clustering analysis of myeloid cell compartment. (A) UMAP depicting interpatient variation, (B) tissue distribution, and (C) patient contribution to the different regions sampled. (D) UMAPs showing marker gene expression in myeloid cells. (E) Scatter plot of M1 versus M2 polarisation for all of the macrophages showing individual cells (grey), and mean values for each subtype (coloured). (F) Heatmap plotting the similarity of myeloid clusters derived in our study (G) Violin plots showing expression levels of SPP1 and GPNMB across all myeloid cell subtypes.
為了進(jìn)一步表征我們數(shù)據(jù)集中的異質(zhì)巨噬細(xì)胞群，我們探索了 DEG 和九個巨噬細(xì)胞cluster的組織富集炊苫。 我們發(fā)現(xiàn)與其他正常組織相比裁厅，六個巨噬細(xì)胞cluster（cluster 11-16）優(yōu)先富集在腫瘤核心/界面中，因此被定義為腫瘤相關(guān)巨噬細(xì)胞（TAM）侨艾。 其余三個cluster（cluster 6执虹、7、8）在正常組織中富集唠梨，被認(rèn)為是組織駐留巨噬細(xì)胞 (TR Mac)袋励。 在六個 TAM cluster中，MHC-II TAM（cluster 14）高表達(dá) HLA-DRB5当叭、APOE 和 APOC1茬故，與腫瘤-正常界面相比，腫瘤核心更富集蚁鳖。 相比之下磺芭，其他五個 TAM cluster在腫瘤核心和界面中顯示出相當(dāng)程度的富集。 促炎性 TAM（clster 11）高表達(dá)趨化因子 CXCL9/10 和 NLRP3 炎癥小體組裝激活劑 GBP1/5醉箕，表現(xiàn)出 M1 極化的主要特征钾腺。
image.png
FN1+ TAM（cluster 15）高表達(dá)纖連蛋白 1 (FN1) 和清道夫受體 MARCO徙垫，此前曾報道其為腎癌中的特定巨噬細(xì)胞亞群。我們發(fā)現(xiàn) FN1+ TAM 可能在 ccRCC 中促腫瘤放棒，這反映在髓源性抑制細(xì)胞 (MDSC) 特征和 M2 極化基因的高表達(dá)上姻报。我們在我們的數(shù)據(jù)集中發(fā)現(xiàn)了一個 SPP1+ TAM cluster（cluster 16），該cluster已在各種癌癥類型中有所報道间螟，但在腎癌中不存在吴旋。有趣的是，我們發(fā)現(xiàn) ccRCC 中的 SPP1+ TAM 表達(dá) GPNMB厢破，并顯示出與先前研究確定的 GPNMB+ TAM 高度相似荣瑟。考慮到我們還確定了一個 GPNMB+ TAM cluster（cluster 15）并且可以在多個 TAM cluster中檢測到 GPNMB 的表達(dá)溉奕，這一發(fā)現(xiàn)表明 SPP1+ TAM 可能代表 GPNMB+ TAM 的一個subset褂傀。除了表達(dá) SPP1，我們發(fā)現(xiàn) SPP1+ TAM 還表達(dá) TREM2 并具有高血管生成評分加勤。 TREM2+ 巨噬細(xì)胞與各種生物和病理過程有關(guān)仙辟，例如肥胖和癌癥。
在三個 TR Mac cluster中鳄梅，TR Mac.2 高表達(dá)白細(xì)胞介素 IL1B 和表皮生長因子受體配體 AREG叠国，這可能反映了其在體內(nèi)平衡組織修復(fù)中的可能作用。 TR Mac.3戴尸，顯示SEPP1和MRC1的高表達(dá)粟焊，并且在正常腎上腺中極其富集。 有趣的是孙蒙，TR Mac.3 表現(xiàn)出極高的 M2 表達(dá)和吞噬特征项棠，并顯示出與促腫瘤 TAM cluster（即 FN1+ TAM）相似的通路激活。 在該數(shù)據(jù)集中挎峦，我們無法清楚地區(qū)分胚胎接種的組織巨噬細(xì)胞與單核細(xì)胞衍生的組織巨噬細(xì)胞香追。
image.png
- 注：Heatmap showing the results of pathway enrichments of macrophage subsets using GSVA analysis
接下來，我們探索了在我們的研究中確定的不同 TR Mac 和 TAM 集群的潛在起源坦胶。 使用 RNA 速度分析透典，我們發(fā)現(xiàn)從循環(huán)單核細(xì)胞到組織中的巨噬細(xì)胞有兩個明顯的定向流動：（1）經(jīng)典 mono.3 到 TR Mac.2 和（2）非經(jīng)典單核細(xì)胞到 TR Mac.1。 TR Mac.1 和 TR Mac.2 然后可能在組織中產(chǎn)生其他巨噬細(xì)胞
image.png
- 注：UMAP with superimposed RNA velocity analysis of the monocyte and macrophage subsets with zoomed windows highlighting possible directional flows from monocytes to macrophages.
為了確定巨噬細(xì)胞亞群與循環(huán)單核細(xì)胞的關(guān)系顿苇，我們利用體細(xì)胞突變進(jìn)行譜系追蹤峭咒，其方式類似于通過共享 TCR 克隆型確定 T 細(xì)胞表型狀態(tài)的關(guān)系。 在這里纪岁，我們構(gòu)建了一個鄰接樹來描述不同單核細(xì)胞和巨噬細(xì)胞cluster的關(guān)系凑队，利用從 scRNA-seq 數(shù)據(jù)calling的體細(xì)胞突變
image.png
- 注：Neighbour-joining tree depicting the relationship of different monocyte and macrophage clusters, utilising the somatic mutations called from scRNA-seq data. The numbers of supporting votes in bootstrapping (100 times) are labelled.
我們發(fā)現(xiàn)循環(huán)單核細(xì)胞與組織中的巨噬細(xì)胞分離，與其他經(jīng)典單核細(xì)胞相比幔翰，非經(jīng)典單核細(xì)胞（cluster 5）與組織中的巨噬細(xì)胞顯示出更密切的關(guān)系顽决。 我們的數(shù)據(jù)支持代表循環(huán)單核細(xì)胞和巨噬細(xì)胞之間中間狀態(tài)的非經(jīng)典單核細(xì)胞短条，大多數(shù)巨噬細(xì)胞似乎來自單核細(xì)胞祖細(xì)胞而不是卵黃囊起源导匣。
Regional heterogeneity of endothelial cells, fibroblasts and epithelial cells
我們在數(shù)據(jù)集中觀察到異質(zhì)基質(zhì)細(xì)胞群才菠。 內(nèi)皮細(xì)胞 (EC) 區(qū)室的亞cluster顯示 11 個具有不同 DEG 和組織位置偏好的cluster。
image.png
Pericytes（cluster 11）贡定，以 RGS5 和 TAGLN 的表達(dá)為特征赋访，優(yōu)先富集在腫瘤核心中。發(fā)現(xiàn)一個小cluster（cluster 9）在腎周脂肪中極度富集缓待，并高度表達(dá) TFF3 和 PDPN蚓耽，因此代表淋巴 EC。其余 9 個cluster代表血管 ECs旋炒，其中cluster 10 代表循環(huán) EC步悠，高表達(dá) TOP2A 和 MKI67。我們確定了三個潛在的腫瘤相關(guān) EC cluster：膠原蛋白 EC瘫镇、IGFBP3+ EC 和 ACKR1+ EC鼎兽，因為它們在腫瘤組織中顯示出相當(dāng)多的富集。在這三個cluster中铣除，發(fā)現(xiàn)膠原蛋白 EC（表達(dá) COL4A1 和 COL15A1）在界面中更加豐富谚咬，這可能通過細(xì)胞外基質(zhì) (ECM) 產(chǎn)生與其他細(xì)胞相互作用。 ACKR1+ EC 特異性表達(dá)非典型趨化因子受體 ACKR1尚粘，其支持免疫細(xì)胞的粘附和組織遷移择卦。我們發(fā)現(xiàn) IGFBP3+ EC 也表達(dá)高水平的免疫抑制酶 IDO1，暗示其在 TME 中的免疫調(diào)節(jié)作用郎嫁。我們發(fā)現(xiàn) CRHBP+ EC 和 IGF2+ EC 優(yōu)先富集在正常腎組織中秉继，而 DNASEL3+ EC 顯示在正常腎上腺中富集
我們在sub-clustering分析中確定了九個成纖維細(xì)胞（Fibro）cluster。 與 EC 類似泽铛，我們發(fā)現(xiàn)了一組成纖維細(xì)胞（cluster 8）高表達(dá)的膠原相關(guān)基因（COL1A1 和 COL6A2）并優(yōu)先富集在interface尚辑。
image.png
這表明不同的產(chǎn)生 ECM 的基質(zhì)細(xì)胞傾向于在界面中富集和共定位，發(fā)揮不同的功能厚宰，包括細(xì)胞外環(huán)境重塑和細(xì)胞間相互作用腌巾。 MMP 纖維的特點是基質(zhì)金屬蛋白酶 MMP2、補體因子 CFD 和 lumican LUM 的高表達(dá)铲觉。 發(fā)現(xiàn) MMP 纖維富含腎周脂肪和腎上腺澈蝙。 Cluster 7 高度表達(dá) MYH11、SNCG 和 RERGL撵幽，因此被認(rèn)為是平滑肌細(xì)胞 (SMC)灯荧，如成纖維細(xì)胞。 發(fā)現(xiàn) SMC 樣纖維在正常腎組織中富集
數(shù)據(jù)集中的正常上皮細(xì)胞群表現(xiàn)出預(yù)期的多樣性盐杂，并在sub-clustering分析中分為 13 個聚類逗载。 我們鑒定了兩個近端腎小管 (PT) 細(xì)胞cluster（均表達(dá) NAT8 和 PDZK1IP1）：cluster 1 具有更高的金屬硫蛋白 MT1H 和 MT1G 表達(dá)哆窿，因此可能代表 PT3 細(xì)胞，而cluster 2 可能代表 PT1/2 細(xì)胞厉斟，因為它表現(xiàn)出更高的 PT2 標(biāo)記 SLC22A6挚躯，和 PT1 標(biāo)記 SLC5A12 和 SLC5A2。
image.png
兩個 Henle 環(huán)（LoH）cluster：cluster 3 代表升細(xì)肢（ATL）細(xì)胞擦秽，表達(dá) CLDN3 和 TACSTD2码荔； cluster 8 代表厚升肢 (TAL) 細(xì)胞，其特征在于 SLC12A1 和 UMOD 的表達(dá)感挥。 確定了三個集合管 (CD) 上皮細(xì)胞cluster：A 型插入細(xì)胞（表達(dá) SLC4A1）缩搅、B 型插入細(xì)胞（表達(dá) SLC4A9）和 CD 主細(xì)胞（表達(dá) AQP2 和 FXYD4）。 我們鑒定了特異性表達(dá) SLC12A3 的遠(yuǎn)曲小管 (DCT) 細(xì)胞（7 cluster）触幼、高表達(dá) SLC8A1 和 CALB1 的連接小管細(xì)胞（10 cluster）硼瓣、顯示 PSCA 和 KRT17 特異性表達(dá)的盆腔尿路上皮細(xì)胞（13 cluster）和足細(xì)胞（11 cluster） ) 專門表達(dá) PTGDS 和 PTPRO。
RCC expression meta-programmes show differential abundance at the tumour-normal interface and affects prognosis
為了探索腫瘤細(xì)胞群中的腫瘤內(nèi)表達(dá)異質(zhì)性置谦，我們首先使用非負(fù)矩陣分解（NMF）定義了由每個腫瘤中的共表達(dá)基因組成的腫瘤內(nèi)表達(dá)程序堂鲤。這些表達(dá)程序代表了每個腫瘤中僅由腫瘤細(xì)胞亞群高度表達(dá)的基因模塊，如代表性腫瘤 PD45816 中的 NMF 結(jié)果所示霉祸≈郏總的來說，我們從 10 個 ccRCC 腫瘤中剖析了 45 個腫瘤內(nèi)表達(dá)程序丝蹭。其中一些程序雖然是個體腫瘤中的亞群事件慢宗，但被發(fā)現(xiàn)由不同的腫瘤共享，因此被定義為 ccRCC 中腫瘤細(xì)胞表達(dá)的meta程序奔穿。通過聚類分析確定了六個meta程序镜沽。第一個meta程序 (MP1) 的特征是 FOS 和 JUN 等基因的表達(dá)，因此代表了腫瘤細(xì)胞中與應(yīng)激反應(yīng)相關(guān)的特征贱田。 MP2 由近端腎小管 (PT) 細(xì)胞特異性表達(dá)的基因 (即 NAT8 和 ACSM2B) 組成缅茉。腫瘤細(xì)胞中 PT 特征的存在證實了先前的發(fā)現(xiàn)，即 PT 細(xì)胞是 ccRCC 的細(xì)胞類型男摧。有趣的是蔬墩，我們發(fā)現(xiàn)第三個meta程序 (MP3) 富含 TGFBI 和 MT2A 等與上皮間質(zhì)轉(zhuǎn)化 (EMT) 相關(guān)的基因。這表明 MP3 可以概括 ccRCC 中的 EMT 過程耗拓，這在之前的 RCC scRNA-seq 研究中尚未報道拇颅。 MP4 由 NEAT1 和 HCG18 等非編碼 RNA 基因組成，可能反映了一些壓力或細(xì)胞死亡 (CD) 相關(guān)的細(xì)胞狀態(tài)乔询。 MP5 的特點是表達(dá) MHC-II 相關(guān)基因樟插，如 CD74 和 HLA-DRA。 MP6中發(fā)現(xiàn)TOP2A、MKI67等基因黄锤，說明該meta程序與腫瘤細(xì)胞增殖有關(guān)搪缨。
接下來，我們整合了來自十個腫瘤的腫瘤細(xì)胞鸵熟，通過去除批次效應(yīng)來減輕患者間的異質(zhì)性副编。通過sub-clustering和 DEG 分析，我們驗證了腫瘤細(xì)胞中六個meta程序的存在旅赢。我們計算了使用 NMF 破譯的六個meta程序的基因評分齿桃，并將它們映射到腫瘤細(xì)胞的 UMAP 上。這再次反映了在腫瘤細(xì)胞中亞群的meta程序的表達(dá)煮盼。有趣的是，我們發(fā)現(xiàn) PT 和 EMT 程序的表達(dá)顯示出反轉(zhuǎn)模式带污，這通過在 TCGA 樣本的大量 RNA-seq 數(shù)據(jù)中計算的 PT 和 EMT 分?jǐn)?shù)之間的反相關(guān)性得到進(jìn)一步證實僵控。此外，我們發(fā)現(xiàn)與腫瘤核心相比鱼冀，EMThigh 腫瘤細(xì)胞在腫瘤-正常界面（腫瘤的前沿）更豐富报破，這反映了 EMT 狀態(tài)代表腫瘤細(xì)胞更具侵襲性和遷移性的事實。 PT/EMT 程序的異質(zhì)表達(dá)與腫瘤細(xì)胞的空間位置偏好相結(jié)合千绪，以個體腫瘤為例充易。例如，在腫瘤 PD45815 中荸型，我們發(fā)現(xiàn)當(dāng)通過 EMT 評分對腫瘤細(xì)胞進(jìn)行排序時盹靴，PT 程序是反向表達(dá)的。同時瑞妇，EMThigh腫瘤細(xì)胞傾向于定位在界面稿静，而PThigh細(xì)胞相對更富集在腫瘤核心的R1和R2區(qū)域。最后辕狰，通過對 TCGA 樣本的大量 RNAseq 數(shù)據(jù)進(jìn)行評分改备，我們發(fā)現(xiàn)根據(jù) TCGA 研究顯示出最差預(yù)后的 TCGA 分子亞型 m3，與其他亞型相比蔓倍，顯示出顯著更高的 EMT 評分但較低的 PT 評分悬钳。這一發(fā)現(xiàn)表明我們的meta程序可以成為患者生存的潛在指標(biāo).
image.png
-
注：RCC cell expression programmes, regional enrichment and prognosis. (A) Heatmap showing expression programmes derived in a representative patient using NMF. (B) Heatmap depicting shared expression meta-programmes across all patients. (C) UMAP representing clusters of tumour cell population. (D) Relative expression scores of meta-programmes in each RCC cell cluster (left) and the distributions of cells with different meta-programmes in tumour core and tumour-normal interface. (E) Cells from patient donors PD45815 and PD45816, ranked by decreasing EMT score with corresponding PT score and cell location. (F) Box plots showing the EMT and PT scores of TCGA samples in different molecular subtypes. ***, P < 0.001 (Wilcoxon rank sum test)
image.png
Cellular interactions associated with spatial location reveal biological insights and promising therapeutic targets
為了表征 RCC 微環(huán)境中不同空間位置的細(xì)胞間通訊，我們使用 CellPhoneDB 評估了正常腎臟偶翅、腫瘤-正常界面和腫瘤核心中主要細(xì)胞類型之間的細(xì)胞間相互作用默勾。有趣的是，界面和腫瘤核心中 12 種主要細(xì)胞類型之間發(fā)生的細(xì)胞間相互作用的數(shù)量相對可比倒堕，即使我們排除了那些涉及腫瘤細(xì)胞的相互作用灾测，也比正常腎組織中的數(shù)量高出約兩倍.接下來，我們比較了在正常腎臟、腫瘤-正常界面和腫瘤核心的不同細(xì)胞類型上表達(dá)的特定配體-受體對介導(dǎo)的細(xì)胞-細(xì)胞相互作用媳搪。首先铭段，該分析揭示了腫瘤與正常微環(huán)境之間相互作用的顯著差異。例如秦爆，在 CD8+ T 細(xì)胞和 ECs 之間序愚，我們發(fā)現(xiàn) EC 募集（CCL5-ACKR1）和免疫抑制（LGALS9-HAVCR2）信號在界面和腫瘤核心中更加活躍，反映了腫瘤中血管生成和免疫抑制水平的升高等限。到正常組織爸吮。其次，我們始終觀察到腫瘤邊緣和核心之間的差異望门。例如形娇，由 PVR-TIGIT 介導(dǎo)的潛在免疫抑制相互作用分別在腫瘤核心中更活躍，而細(xì)胞生長/遷移相關(guān)相互作用 IGF-IGF2R 在界面中更活躍筹误。第三桐早，界面處的腫瘤細(xì)胞表達(dá)轉(zhuǎn)錄物，預(yù)測其介導(dǎo)獨特的細(xì)胞間相互作用厨剪，可能促進(jìn)腫瘤發(fā)生哄酝。例如，在界面中祷膳，與腫瘤核心中的相比陶衅，髓系細(xì)胞表達(dá)增強的細(xì)胞遷移相關(guān)信號，可以與腫瘤細(xì)胞相互作用（THBS1-整合素α3β1）直晨。
image.png
- 注： (A) Comparison of the number of cell–cell interactions between cells present in the adjacent normal kidney, tumour-normal interface, and the tumour core. (B) Dot plots representing mean expression levels and significance of important differentially expressed ligand receptor pairs by spatial region.
我們的結(jié)果表明搀军，具有高表達(dá) EMT 特征的腫瘤細(xì)胞（EMThigh 腫瘤細(xì)胞）優(yōu)先定位于腫瘤的前緣。這促使我們探索界面處是否存在任何可能促進(jìn)腫瘤細(xì)胞 EMT 的活躍細(xì)胞間相互作用抡秆。我們使用 NicheNet 分析將來自 TME 中細(xì)胞的配體與腫瘤細(xì)胞中的 EMT 程序聯(lián)系起來奕巍。從該分析中，我們發(fā)現(xiàn)巨噬細(xì)胞表達(dá)的配體可能調(diào)節(jié)在腫瘤細(xì)胞上表達(dá)的大量 EMT 基因儒士。特別是的止，巨噬細(xì)胞來源的 IL1B 對這些 EMT 基因顯示出高度和廣泛的調(diào)節(jié)潛力，可能是通過腫瘤細(xì)胞中表達(dá)的受體 IL1R1着撩。有趣的是诅福，我們發(fā)現(xiàn) IL1B 由 TR Mac.2 特異性表達(dá)，其再次在腫瘤正常界面優(yōu)先富集拖叙∶ト螅總之，我們的研究結(jié)果表明薯鳍，表達(dá) IL1B 的巨噬細(xì)胞 (TR Mac.2) 優(yōu)先存在于腫瘤-正常界面咖气，通過產(chǎn)生 IL1B 上調(diào)腫瘤前沿腫瘤細(xì)胞中的 EMT 程序。這種致癌途徑可能最終促進(jìn)腫瘤細(xì)胞的遷移和侵襲
image.png
- 注： (C) Heatmap depicting the potential regulation of genes expressed by the EMT meta-programme and ligands expressed by macrophages. (D) Schematic illustrating the importance of IL1B signalling between IL1B+ macrophages and RCC cells in promoting EMT.
DISCUSSION
我們使用基于多區(qū)域的基因組和單細(xì)胞轉(zhuǎn)錄組測序來表征 ccRCC 的表型異質(zhì)性和多細(xì)胞生態(tài)系統(tǒng)。 總體而言崩溪，我們的研究描繪了 ccRCC 的 TME 的綜合圖譜以及 ccRCC 中已建立的 ITH浅役，包括腫瘤細(xì)胞和免疫/基質(zhì)細(xì)胞的表型分類，以及它們在 TME 中的細(xì)胞間通訊伶唯，主要與其地理定位相關(guān)觉既。
擴增的 CD8+ TCR 克隆型內(nèi)的細(xì)胞在很大程度上受到耗竭評分的限制谣旁。最近在黑色素瘤中報道了類似的觀察結(jié)果丽柿。克隆型的表型限制可能主要與給定克隆的時間成熟有關(guān)蛛壳，而不是環(huán)境因素粹断，因為個體腫瘤在完全不同的狀態(tài)下都具有克隆型符欠。除了表型限制之外，我們還發(fā)現(xiàn)擴增的 TCR 克隆型在一個或多個宏觀腫瘤活檢中也經(jīng)常受到空間限制姿染。由于我們的研究中觀察到的體細(xì)胞突變的 ITH 有限背亥，因此無法通過暴露于不同的突變相關(guān)新抗原來完全解釋 TCR 克隆擴增的這種空間限制。我們無法在 TME 中定義任何其他可以預(yù)測這種克隆型空間異質(zhì)性的因素悬赏。感知到的 T 細(xì)胞克隆型的隨機定位可能是物理和環(huán)境因素驅(qū)動細(xì)胞從外周循環(huán)到腫瘤駐留的初始遷移的結(jié)果÷玻可以采用縱向取樣策略或方法來確定準(zhǔn)確的 T 細(xì)胞系統(tǒng)發(fā)育闽颇，以詢問 T 細(xì)胞擴增和遷移的精確時間。我們對擴增 TCR 克隆區(qū)域限制的觀察是否在其他癌癥類型中更廣泛地被發(fā)現(xiàn)寄锐，可能需要部署類似采樣和測序策略的額外研究兵多。
外周 TCR 用于非侵入性癌癥檢測和監(jiān)測的效用顯示出前景，尤其是在循環(huán)腫瘤 DNA 片段稀缺的 RCC 中橄仆。 雖然我們發(fā)現(xiàn)在血液和腫瘤區(qū)域都存在許多擴增克隆型剩膘，但我們觀察到衰竭程度與外周血中檢測到 TCR 克隆的概率呈負(fù)相關(guān)，以至于在血液中很少檢測到衰竭克隆型盆顾。 研究結(jié)果表明怠褐，一旦 T 細(xì)胞克隆浸潤到腫瘤中并經(jīng)歷從激活到功能障礙的表型轉(zhuǎn)變，它們就很少再循環(huán)您宪，這可能是由于 CD69 所證明的組織駐留表型奈懒。 因此，對腫瘤反應(yīng)性 TCR 的外周采樣更有可能檢測到耗盡的腫瘤駐留克隆的前因宪巨，而不是那些目前在腫瘤中活躍的克隆.
作者開發(fā)了一種策略 (deSCeRNAMut)磷杏，可以根據(jù)基于液滴的 scRNA-seq 數(shù)據(jù)準(zhǔn)確檢測不同細(xì)胞群中的體細(xì)胞突變。使用許多過濾指標(biāo)（包括不同細(xì)胞類型譜系之間共享胚胎后突變的不可信性）消除了缺乏一致覆蓋捏卓、低讀取深度和容易出錯的測序讀取的主要挑戰(zhàn)极祸。我們檢測到我們數(shù)據(jù)集中不同細(xì)胞譜系的體細(xì)胞突變，并確定了骨髓細(xì)胞中的高度克隆擴增。我們使用從這種方法calling的體細(xì)胞突變來構(gòu)建巨噬細(xì)胞和單核細(xì)胞中的相鄰連接樹遥金，以推斷非經(jīng)典單核細(xì)胞可能是腎癌中循環(huán)單核細(xì)胞和大多數(shù)組織駐留巨噬細(xì)胞之間的中間狀態(tài)浴捆。我們設(shè)想在未來，使用空間成像技術(shù)來可視化一系列細(xì)胞類型中表達(dá)基因的突變將有助于破譯多細(xì)胞 TME 的系統(tǒng)發(fā)育組織汰规。
EMT meta程序是通過每個患者的腫瘤細(xì)胞亞群的表達(dá)來定義的汤功，并且在我們的研究中由多個 ccRCC 腫瘤共享。 更豐富的腫瘤細(xì)胞群和幫助規(guī)避具有挑戰(zhàn)性的批次變化的方法的使用使我們能夠發(fā)現(xiàn)這個以前未報告的特征溜哮。 EMT 程序在 ccRCC 腫瘤細(xì)胞中的表達(dá)與 PT 程序（一種上皮表達(dá)特征）的表達(dá)呈負(fù)相關(guān)滔金。 同時，ccRCC中的EMThigh腫瘤細(xì)胞傾向于定位于腫瘤-正常界面茂嗓，這是腫瘤的前沿和遷移邊緣餐茵。 這些發(fā)現(xiàn)與頭頸癌的 scRNA-seq 研究中報道的結(jié)果相似，反映了 EMT 的定義特征：細(xì)胞中上皮特征的喪失述吸，有利于促進(jìn)其遷移和侵襲能力忿族。
分析細(xì)胞間相互作用揭示了與 ccRCC 的 TME 中不同空間定位相關(guān)的預(yù)測細(xì)胞間通訊的異質(zhì)性。特別是蝌矛，通過使用 NicheNet 連接配體和目標(biāo)基因道批，我們發(fā)現(xiàn) IL1B，由在腫瘤-正常界面 (TR Mac.2) 富集的組織駐留巨噬細(xì)胞亞群特異性表達(dá)入撒，可能促進(jìn)腫瘤細(xì)胞進(jìn)行 EMT隆豹。據(jù)報道，IL1B 的表達(dá)與 RCC 的腫瘤分期呈正相關(guān)茅逮，并且與招募到癌癥基因組圖譜的患者中 RCC 患者的較差預(yù)后相關(guān)璃赡。此外，在 RCC 中抑制 IL1B 已顯示在 RCC 的同源小鼠模型中誘導(dǎo)腫瘤消退献雅。 IL1B 阻斷劑也被證明可以減少動脈粥樣硬化患者肺癌的發(fā)生率碉考，目前正在幾項臨床試驗中研究其使用。在我們的數(shù)據(jù)中挺身，我們表明導(dǎo)致 IL1B 在 RCC 中不利作用的潛在機制通過巨噬細(xì)胞衍生的 IL1B 信號傳導(dǎo)促進(jìn) EMT 起作用侯谁。利用這一途徑可能在治療上是有用的。
Method
scRNA-seq merge and QC
image.png
Cell type sub-clustering and annotation
image.png
image.png
Pseudotime inference, TCR analysis
image.png
De Novo Mutation Calling from scRNA-seq Data 重點
image.png
image.png
image.png
Inferring copy number variations based on scRNA-seq data
image.png
image.png
Cell subtype abundance in different tissues
image.png
Correlation between spatial, somatic RCC evolution and TCR clonotype evolution
image.png
image.png
Gene set enrichment analysis and gene signature scoring in macrophage population
image.png
RNA velocity analysis
image.png
image.png
Lineage tracing using scRNA-seq called somatic mutations
image.png
image.png
Deciphering intra-tumour expression programmes and meta-programmes
image.png
image.png
Integrating and analyzing tumour cells from different patients
image.png
細(xì)胞通訊分析
image.png
生活很好瞒渠，有你更好
