影響因子：27.7

研究概述：NK細(xì)胞是固有免疫淋巴細(xì)胞喧半，能直接殺死病毒感染的細(xì)胞奴迅，并調(diào)控適應(yīng)性免疫反應(yīng)的早期階段，在免疫反應(yīng)中發(fā)揮重要作用挺据。NK細(xì)胞通常分為兩類：CD56bright和CD56dim取具，但這種分類過于簡化脖隶。通過單細(xì)胞技術(shù)的分析，健康個(gè)體中的NK細(xì)胞具有高達(dá)10萬個(gè)獨(dú)特亞群的廣泛表型多樣性暇检。這種多樣性源于隨機(jī)表達(dá)的产阱、胚系編碼的激活和抑制受體的組合，這些受體與HLA I類結(jié)合块仆，調(diào)節(jié)NK細(xì)胞功能构蹬。另一層次的NK細(xì)胞多樣性是從初始CD56bright NK細(xì)胞到終末分化的適應(yīng)性CD62L?NKG2C+CD57+KIR+CD56dim NK細(xì)胞的連續(xù)分化，與巨細(xì)胞病毒CMV感染有關(guān)悔据。在本篇研究中庄敛，作者建立了NK細(xì)胞單細(xì)胞轉(zhuǎn)錄參考圖譜，解析了血液和正常組織中NK細(xì)胞分化期間的基因表達(dá)趨勢蜜暑。通過參考圖譜分析各種條件下的細(xì)胞差異基因铐姚，識(shí)別出新的NK細(xì)胞狀態(tài)。作者使用來自44,640個(gè)PB-NK細(xì)胞和27,732個(gè)組織駐留TrNKNK細(xì)胞的數(shù)據(jù)肛捍，研究了來自427名患者的7種實(shí)體瘤中的腫瘤浸潤TiNK細(xì)胞。發(fā)現(xiàn)TrNK和TiNK細(xì)胞具有明顯的組織駐留特征之众。在泛癌圖譜中拙毫，作者確認(rèn)了六種NK細(xì)胞功能狀態(tài)，其中CD56bright和細(xì)胞毒性效應(yīng)CD56dim狀態(tài)在多種腫瘤類型中普遍富集棺禾∽禾悖患者生存數(shù)據(jù)分析顯示，效應(yīng)CD56dim與CD56bright狀態(tài)的高比例與骨肉瘤和黑色素瘤患者的良好生存相關(guān)膘婶。這揭示了TME如何導(dǎo)致NK細(xì)胞功能障礙缺前，并可能為設(shè)計(jì)新的NK細(xì)胞治療策略提供靈感。具體結(jié)果如下：

研究結(jié)果：

NK細(xì)胞亞群注釋

為了建立一個(gè)腫瘤浸潤NK細(xì)胞的泛癌圖譜悬襟，作者首先在轉(zhuǎn)錄水平上定義了NK細(xì)胞的分化衅码。作者對(duì)七名健康供體的總NK細(xì)胞群體進(jìn)行了單細(xì)胞RNA測序，并整合了五個(gè)公開的單細(xì)胞數(shù)據(jù)集成一個(gè)均勻的人群（圖a）脊岳。上文提及CD56bright和CD56dim NK細(xì)胞代表兩個(gè)不同的分化階段逝段，因此作者基于CD56 bright/ dim signature來對(duì)單細(xì)胞數(shù)據(jù)進(jìn)行評(píng)分，以識(shí)別CD56bright和CD56dim NK細(xì)胞（圖b）割捅。CD56dim NK細(xì)胞異質(zhì)性較大的群體通常根據(jù)抑制性和激活性受體的選擇性表達(dá)來表征功能上不同的亞群奶躯。因此，為了識(shí)別與這些功能階段相關(guān)的預(yù)測基因特征亿驾，作者從兩名供體中分選了相等數(shù)量的CD56bright NK細(xì)胞和四個(gè)CD56dim NK細(xì)胞亞群（NKG2A+KIR-CD57-嘹黔，NKG2A-self-KIR+CD57-，NKG2A-nonself-KIR+CD57-莫瞬，NKG2A-self-KIR+CD57+或NKG2A-/+self-KIR+CD57+NKG2C+）（圖c）儡蔓。這五個(gè)轉(zhuǎn)錄上不同的NK細(xì)胞亞群被重新命名以反映它們的成熟階段：‘CD56bright’醉锄、‘早期CD56dim’、‘中期CD56dim’浙值、‘晚期CD56dim’和‘適應(yīng)性CD56dim’（圖c）恳不。作者接下來訓(xùn)練了一個(gè)半監(jiān)督模型--scANVI，利用識(shí)別的NK細(xì)胞亞群基因特征來預(yù)測和推斷整理的整體NK細(xì)胞亞群注釋开呐。結(jié)果發(fā)現(xiàn)烟勋，轉(zhuǎn)錄上不同的亞群（CD56bright和適應(yīng)性NK細(xì)胞這兩個(gè)亞群）注釋準(zhǔn)確性高，而轉(zhuǎn)錄重疊較高的亞群注釋準(zhǔn)確性略低（圖d）筐付。使用該模型卵惦，作者在亞群水平注釋包含12個(gè)患者的23,253個(gè)單細(xì)胞轉(zhuǎn)錄組的總NK細(xì)胞數(shù)據(jù)集。scANVI模型捕獲了亞群的轉(zhuǎn)錄特征瓦戚，并用于識(shí)別差異表達(dá)的基因（圖e-f）沮尿。為了驗(yàn)證這一注釋模型，作者又對(duì)總NK細(xì)胞數(shù)據(jù)集（12個(gè)患者）進(jìn)行了無偏聚類分析（Leiden）较解，識(shí)別了與前面注釋的5個(gè)NK細(xì)胞亞群緊密匹配的5個(gè)簇（圖g）畜疾。部分中期CD56dim注釋的NK細(xì)胞與晚期CD56dim注釋的NK細(xì)胞在聚類4中聚集，可能對(duì)應(yīng)于群體中更成熟的細(xì)胞（圖h）印衔。有/無適應(yīng)性NK細(xì)胞擴(kuò)展的患者對(duì)比中（NK細(xì)胞可以在遇到某些病毒感染如CMV后經(jīng)歷適應(yīng)性擴(kuò)展）啡捶，晚期CD56dim和適應(yīng)性NK細(xì)胞亞群的頻率也會(huì)相應(yīng)發(fā)生變化（圖i）。

特定的基因調(diào)控網(wǎng)絡(luò)

為了在轉(zhuǎn)錄水平上解析驅(qū)動(dòng)NK細(xì)胞分化的調(diào)控基因路徑奸焙，作者使用Palantir和RNA velocity進(jìn)行了擬時(shí)間分析瞎暑。作者基于最低CD56dim評(píng)分定義起始細(xì)胞，識(shí)別了兩個(gè)終末期的細(xì)胞（橙色）与帆，預(yù)測屬于晚期CD56dim和適應(yīng)性CD56dim細(xì)胞亞群（圖a）了赌。作者為了驗(yàn)證這一軌跡，利用scVelo27中的動(dòng)態(tài)模型計(jì)算剪接與未剪接RNA速度玄糟，在沒有預(yù)定義類似上面起始細(xì)胞的情況下推斷擬時(shí)間勿她，結(jié)果確認(rèn)了從CD56bright NK細(xì)胞亞群開始并終止于適應(yīng)性CD56dim細(xì)胞亞群的軌跡（圖b）。既然識(shí)別了兩種不同的終末期的細(xì)胞：晚期CD56dim和適應(yīng)性CD56dim細(xì)胞亞群茶凳。作者對(duì)這兩細(xì)胞進(jìn)一步分析其來源嫂拴。在兩種來源類型中，早期CD56dim NK細(xì)胞形成了CD56bright和其余CD56dim亞群之間的連接環(huán)節(jié)（圖c,d）贮喧。雖然適應(yīng)性分化中NK細(xì)胞向中期CD56dim細(xì)胞進(jìn)展筒狠；常規(guī)分化中則轉(zhuǎn)向中期/晚期CD56dim群體（圖c,d）。在建立了NK細(xì)胞分化的時(shí)間軸后箱沦，作者利用廣義加性模型（GAMs）計(jì)算每個(gè)基因隨時(shí)間變化的基因表達(dá)趨勢辩恼，這些基因聚集成五個(gè)不同的趨勢（圖e）。在常規(guī)和適應(yīng)性分化中表達(dá)差異的基因（如KLRC2、CD52和IL32）在常規(guī)晚期CD56dim分化中聚集于趨勢1灶伊，在適應(yīng)性分化中聚集于趨勢4（圖e）疆前。基于雙分化模型聘萨，作者構(gòu)建了按五個(gè)基因趨勢分層的基因調(diào)控網(wǎng)絡(luò)（GRNs）竹椒，識(shí)別了跨偽時(shí)間的主導(dǎo)轉(zhuǎn)錄因子（TFs）及其已知的下游靶基因（圖f）。趨勢1由從CD56bright到CD56dim細(xì)胞分化過程中下調(diào)的基因主導(dǎo)米辐，包括先前報(bào)道的TFs（MYC胸完，LEF1，RUNX2）翘贮，涉及Notch信號(hào)通路的RBPJ赊窥，視黃酸受體（RXRA）和調(diào)節(jié)ID2表達(dá)的TFs（HOXA9，HOXA10）狸页。

通過遷移學(xué)習(xí)生成泛癌圖譜

在外周血（PB）中轉(zhuǎn)錄定義NK細(xì)胞分化后锨能，作者繼續(xù)使用scVI訓(xùn)練第二個(gè)模型M2，結(jié)合136名患者的6個(gè)健康組織（前列腺芍耘、肺址遇、胰腺、皮膚齿穗、乳腺傲隶、大腦）scRNA-seq數(shù)據(jù)，生成健康人群的PB-NK + TrNK參考圖譜（圖a）窃页。使用scANVI對(duì)組織特定數(shù)據(jù)集進(jìn)行整合和注釋，并在泛組織水平（圖b）和各個(gè)組織內(nèi)識(shí)別感興趣的免疫亞群复濒。接著脖卖，作者對(duì)來自427名腫瘤患者患者在泛癌水平（圖c,d）和各個(gè)腫瘤類型內(nèi)（包括7個(gè)實(shí)體腫瘤：前列腺癌、肺癌巧颈、黑色素瘤畦木、胰腺癌、乳腺癌砸泛、膠質(zhì)母細(xì)胞瘤和骨肉瘤）進(jìn)行注釋和整合十籍。為了評(píng)估注釋NK細(xì)胞在組織和腫瘤數(shù)據(jù)集中的組織駐留狀態(tài)，作者利用文獻(xiàn)中組織駐留特征以及作者自己的圖譜衍生出了組織駐留（atlas-TR）NK細(xì)胞的特征（圖e）唇礁。結(jié)果發(fā)現(xiàn)勾栗，相比CD56dim NK細(xì)胞，CD56bright NK細(xì)胞在正常組織和腫瘤中總體上對(duì)組織駐留征評(píng)分較高（圖e）盏筐。最后围俘，作者使用遷移學(xué)習(xí)（scArches44），將CD56bright和CD56dim注釋的TiNK細(xì)胞映射到參考圖譜（PB-NK，TrNK）界牡，生成最終模型（M3）--泛癌NK圖譜（圖f）簿寂。PB、組織和腫瘤中的CD56bright和CD56dim亞群聚集在一起（圖g,h）宿亡。

不同組織和腫瘤中NK細(xì)胞亞群頻率的變化

對(duì)健康組織和腫瘤注釋的免疫亞群進(jìn)行泛癌比較（圖b,d）顯示常遂，健康組織的血漿細(xì)胞和初始B細(xì)胞比例增加，而腫瘤組織的CD56dim NK細(xì)胞挽荠、單核細(xì)胞克胳、樹突狀細(xì)胞、NK T細(xì)胞坤按、效應(yīng)記憶/效應(yīng)T輔助細(xì)胞毯欣、效應(yīng)記憶/再表達(dá)CD45RA的細(xì)胞毒性T細(xì)胞和駐留記憶細(xì)胞毒性T細(xì)胞的比例減少（圖a）。在乳腺癌中臭脓，總免疫細(xì)胞中CD56bright NK細(xì)胞比例增加酗钞。而在黑色素瘤中，CD56dim NK細(xì)胞比例增加来累，但在非小細(xì)胞肺癌和乳腺癌中比例減少（圖a-c）砚作。作者進(jìn)一步使用第一部分構(gòu)建好的亞群訓(xùn)練模型（M1）對(duì)NK細(xì)胞進(jìn)行亞群水平注釋（圖d,e）。大多數(shù)腫瘤類型（圖d,e）中健康組織與腫瘤之間存在CD56bright比例的偏斜嘹锁，通過流式細(xì)胞術(shù)在非小細(xì)胞肺癌隊(duì)列中得到同樣的驗(yàn)證（圖f）葫录。在腫瘤浸潤NK細(xì)胞（TiNKs）中，中期CD56dim群體普遍減少（圖d,e）领猾。對(duì)非小細(xì)胞肺癌隊(duì)列CD56dim群體的蛋白質(zhì)水平注釋還發(fā)現(xiàn)米同，與健康血液對(duì)照相比，早期和中期CD56dim亞群減少摔竿，晚期CD56dim亞群略有增加（圖g）面粮。

NK細(xì)胞六種功能上不同的細(xì)胞狀態(tài)

作者實(shí)施了一種無偏方法（Milo47）通過識(shí)別6,932個(gè)獨(dú)立的細(xì)胞鄰域來確定泛癌NK細(xì)胞圖譜中的細(xì)胞狀態(tài)。這6,932個(gè)鄰域被分為6個(gè)不同的鄰域組继低，并測試了TiNK細(xì)胞和參考NK細(xì)胞（Ref-NK細(xì)胞）之間的差異豐富度（圖a）熬苍。鄰域組1和2由在TiNK細(xì)胞中顯著豐富的鄰域組成，而組6包括在Ref-NK細(xì)胞中豐富的鄰域（圖b）袁翁。接下來柴底，作者使用注釋模型（M1）可視化每個(gè)組內(nèi)NK細(xì)胞亞群的分布。組1和2富集了CD56bright細(xì)胞粱胜，而組3-6主要由CD56dim NK細(xì)胞亞群組成（圖c）柄驻。先前識(shí)別的PB-NK細(xì)胞分化的主導(dǎo)TF調(diào)控子確認(rèn)組1和2為兩個(gè)CD56bright狀態(tài)，而組3-6為四個(gè)CD56dim NK細(xì)胞狀態(tài)（圖d）年柠。接著凿歼，作者依據(jù)細(xì)胞狀態(tài)特異性GRNs褪迟、DEGs、基因集富集分析和特征評(píng)分注釋了以下的NK細(xì)胞狀態(tài)信息：應(yīng)激CD56bright（組1）答憔、典型CD56bright（組2）味赃、效應(yīng)CD56dim（組3）、適應(yīng)性CD56dim（組4）虐拓、活化CD56dim（組5）和典型CD56dim（組6）（圖e–n）心俗。比較應(yīng)激與典型CD56bright狀態(tài)（組1對(duì)組2）發(fā)現(xiàn)：ATF3調(diào)控子（細(xì)胞應(yīng)激反應(yīng)）、MAFF調(diào)控子（缺氧誘導(dǎo)）和眾多熱休克蛋白表達(dá)增加（圖e,g）蓉驹。應(yīng)激CD56bright細(xì)胞狀態(tài)在免疫抑制通路（轉(zhuǎn)化生長因子TGFβ信號(hào)城榛、缺氧、活性氧ROS）中得分較高态兴，并表現(xiàn)出糖酵解狠持、膽固醇穩(wěn)態(tài)、脂肪酸代謝激活增加（圖g,j–l）瞻润。此外喘垂，低NK細(xì)胞細(xì)胞毒性評(píng)分表明應(yīng)激CD56bright這一細(xì)胞狀態(tài)功能降低，該狀態(tài)在所有七種腫瘤類型中普遍表達(dá)較高（圖m,o）绍撞。與TME中CD56bright細(xì)胞浸潤增加一致正勒，典型CD56bright細(xì)胞狀態(tài)在與健康組織相比的五種腫瘤類型中也表現(xiàn)出豐富度增加，兩個(gè)CD56bright組均表現(xiàn)出免疫調(diào)節(jié)分子（如XCL1傻铣、XCL2和IFNG）表達(dá)增加（圖n–o）章贞。在所有腫瘤類型中，CD56dim狀態(tài)中的效應(yīng)狀態(tài)最常見非洲，富集了頂端連接鸭限、肌動(dòng)蛋白和細(xì)胞骨架相關(guān)基因以及效應(yīng)分子（圖h）。該狀態(tài)在表型上富集了中期和晚期CD56dim NK細(xì)胞亞群两踏，在NK細(xì)胞毒性和氧化磷酸化中得分較高里覆，重要的是，在免疫抑制中得分較低（圖i,k–m）缆瓣。活化CD56dim狀態(tài)的特點(diǎn)是缺氧增加虹统、營養(yǎng)運(yùn)輸?shù)鞍咨险{(diào)和mTORC1–Myc軸上調(diào)（圖i,k）弓坞。最后，典型CD56dim狀態(tài)表現(xiàn)出低應(yīng)激評(píng)分和高細(xì)胞毒性評(píng)分车荔，并與IFN渡冻、腫瘤壞死因子TNF和JAK/STAT信號(hào)相關(guān)（圖i–j,m）。

TME中特定狀態(tài)的信號(hào)傳導(dǎo)與功能的聯(lián)系

為了闡明TME對(duì)識(shí)別出的六種功能狀態(tài)的NK細(xì)胞的影響忧便，作者使用CellChat推斷細(xì)胞間通訊族吻。組1和組2 NK細(xì)胞狀態(tài)在四條主要通訊路徑中在所有腫瘤類型中均表現(xiàn)出較強(qiáng)的輸入信號(hào)（圖a）。在非小細(xì)胞肺癌中，高表達(dá)的CD44超歌、CXCR4和CD74在組1和組2 NK細(xì)胞上砍艾，大多來自成纖維細(xì)胞、內(nèi)皮細(xì)胞巍举、腫瘤細(xì)胞和巨噬細(xì)胞的信號(hào)（COLLAGEN脆荷、MIF、LAMININ）匯聚在這些細(xì)胞上懊悯，促進(jìn)了輸入信號(hào)的增強(qiáng)（圖b,c）蜓谋。此外，組1和組2 NK細(xì)胞主要通過主要組織相容性復(fù)合物I（MHC-I）（HLA-E/KLRC1）通路接收抑制性輸入炭分，這歸因于這些細(xì)胞狀態(tài)中KLRC1的高表達(dá)（圖a,d）桃焕。因此，組1和組2細(xì)胞狀態(tài)更易接受通過上調(diào)的CD44捧毛、CXCR4观堂、CD74和KLRC1所介導(dǎo)的TME誘導(dǎo)的免疫抑制信號(hào)。為了理解NK細(xì)胞如何通過免疫調(diào)節(jié)作用影響TME岖妄，作者將分析重點(diǎn)放在主要由NK細(xì)胞進(jìn)行的輸出信號(hào)傳導(dǎo)上型将。作者確定了三條信號(hào)路徑（CC趨化因子配體（CCL）、蛋白酶激活抑制劑（PARs）荐虐、IFN-II）七兜，通過這些路徑，NK細(xì)胞主要與樹突狀細(xì)胞福扬、巨噬細(xì)胞腕铸、成纖維細(xì)胞和內(nèi)皮細(xì)胞進(jìn)行通訊（圖e,f）。效應(yīng)NK細(xì)胞狀態(tài)（組3）在轉(zhuǎn)錄水平上高表達(dá)的顆粒酶A可誘導(dǎo)TME中表達(dá)F2R的細(xì)胞（如成纖維細(xì)胞）的凋亡（圖g）铛碑。與健康血液對(duì)照相比狠裹，非小細(xì)胞肺癌患者的中央腫瘤樣本中的CD56dim NK細(xì)胞顆粒酶A的表達(dá)頻率和強(qiáng)度均減少，提示NK細(xì)胞在腫瘤內(nèi)釋放顆粒酶A（圖h,i）汽烦。通過HLA-E軸的抑制信號(hào)顯著抑制CD56bright和CD56dim NK細(xì)胞的脫顆粒和顆粒酶B釋放涛菠，通過與預(yù)先用IFN-γ刺激以上調(diào)HLA-E表達(dá)的A549靶細(xì)胞共培養(yǎng)證明這一點(diǎn)（圖j,k）。使用抗NKG2A抗體阻斷NKG2A–HLA-E軸撇吞，顯著恢復(fù)了功能俗冻，包括脫顆粒和顆粒酶B釋放（圖j,k）。

細(xì)胞狀態(tài)比例預(yù)測患者結(jié)局

作者將黑色素瘤SKCM牍颈、肺癌NSCLC和膠質(zhì)瘤GBM的空間RNA測序組織切片進(jìn)行解卷積迄薄，以識(shí)別這些數(shù)據(jù)集中的細(xì)胞類型（圖a）。SKCM煮岁、NSCLC和GBM的主要免疫亞群組成分析在腫瘤類型之間存在很大差異讥蔽，但在測序技術(shù)（scRNA-seq與空間-seq）之間高度一致（圖b）涣易。因?yàn)镹K細(xì)胞比例在SKCM最高（圖b），作者進(jìn)一步將注釋的NK細(xì)胞細(xì)分為CD56bright和CD56dim亞群（圖c）及細(xì)分細(xì)胞狀態(tài)（圖d）冶伞。效應(yīng)（組3）和典型（組6）CD56dim狀態(tài)在與NK細(xì)胞毒性相關(guān)的基因中得分較高（圖e）新症。類似地，應(yīng)激反應(yīng)相關(guān)基因及免疫抑制相關(guān)基因（ROS碰缔，缺氧）在應(yīng)激CD56bright（組1）狀態(tài)中得分最高（圖f,g）账劲。鑒于在血液、組織和實(shí)體腫瘤中的scRNA-seq和空間-seq數(shù)據(jù)集中識(shí)別出的六種功能狀態(tài)的NK細(xì)胞金抡，作者繼續(xù)測試這些細(xì)胞狀態(tài)的臨床相關(guān)性瀑焦，使用包含生存數(shù)據(jù)的TCGA數(shù)據(jù)。結(jié)果發(fā)現(xiàn)較高的效應(yīng)CD56dim：應(yīng)激CD56bright NK狀態(tài)特征比率預(yù)測了SARC和SKCM患者的更好生存率（圖h）梗肝。

研究總結(jié)：

本文研究了腫瘤浸潤自然殺傷（NK）細(xì)胞的泛癌特征及其在腫瘤微環(huán)境中的功能榛瓮。研究通過生成健康人血液和組織來源的NK細(xì)胞單細(xì)胞轉(zhuǎn)錄圖譜，分析了腫瘤微環(huán)境誘導(dǎo)的擾動(dòng)巫击。通過將39個(gè)數(shù)據(jù)集中的腫瘤浸潤NK細(xì)胞轉(zhuǎn)錄組整合到參考圖譜中禀晓，發(fā)現(xiàn)了六種功能上不同的NK細(xì)胞狀態(tài)，包括易受腫瘤微環(huán)境免疫抑制的功能障礙CD56bright狀態(tài)和抵抗腫瘤微環(huán)境的細(xì)胞毒性CD56dim狀態(tài)坝锰。文章揭示了NK細(xì)胞在不同分化階段的特定基因調(diào)控網(wǎng)絡(luò)粹懒，并通過遷移學(xué)習(xí)生成了一個(gè)泛癌NK細(xì)胞圖譜，顯示了組織駐留和腫瘤浸潤NK細(xì)胞的特征顷级。研究表明凫乖，不同腫瘤類型中NK細(xì)胞亞群頻率的變化與患者預(yù)后密切相關(guān)，效應(yīng)CD56dim狀態(tài)的高比例與黑色素瘤和骨肉瘤患者的生存改善有關(guān)弓颈。這些發(fā)現(xiàn)為NK細(xì)胞在抗腫瘤免疫中的作用及其作為治療靶點(diǎn)提供了新的見解帽芽。