肺癌是人类肿瘤的头号杀手。虽然近年来靶向治疗和免疫治疗取得了突破性进展,但患者仍会出现耐药、免疫治疗无效等情况。在单细胞水平上研究可以精确了解患者体内、以及患者间的异质性,提高细胞差异性研究的分辨率,帮助更好指导治疗。
-
单细胞测序数据(10x Genomics):
入组5例未经治疗的NSCLC患者,其中2例腺癌,2例鳞癌,1例大细胞癌。采集每例患者的3份肿瘤组织样本、以及1份同一肺叶下远端正常肺组织样本,进行单细胞测序。
入组3例未经治疗的NSCLC患者作为验证集。
-
外部数据
TCGA数据库以及其它研究中的1572名NSCLC RNA-Seq数据及预后信息(1027 LUAD 和 545 LUSC)。
-
共获得了52,698个细胞,其中有39,323细胞来源于肿瘤; 13,375个来源于正常组织。
-
52,698个细胞中覆盖了约2000万个转录本。
-
利用2192个高变基因进行主成分分析,确定了8个不同的亚群:内皮细胞、成纤维细胞、B细胞、骨髓细胞、T细胞、肺泡细胞、上皮细胞、癌细胞。
-
进一步对每个亚群再次聚类,划分为52种不同的基质细胞亚群和12种癌细胞亚群。
-
内皮细胞进一步划分为6个亚型:
-
淋巴管内皮细胞 (PDPN 和 PROX1)
-
血管内皮细胞 (FLT1+)
-
2个来源于肿瘤的内皮细胞 (IGFBP3+ 和 SPRY1+)
-
2个来源于非肿瘤的内皮细胞 (MT2A+ 和 EDNRB+)
-
肿瘤来源的内皮细胞转录本表达丰度高于正常组织来源的转录本。
-
通路分析发现肿瘤来源的内皮细胞的代谢通路上调,而免疫活性的通路下调。其中下调的基因包括抗原呈递基因 (CCL2)、趋化基因 (ICAM1)以及免疫归巢基因(HLA-DQA1, HLA-DQB1)。该结果说明肿瘤内皮细胞受到免疫微环境的影响,下调抗原呈递和免疫细胞归巢活性,进而有利于肿瘤的免疫耐受。
-
在肿瘤的淋巴内皮细胞中,利用SCENIC的方法,发现了可能导致肿瘤特异性表型以及削弱免疫活性的基因调控对。
-
成纤维细胞可进一步划分为6个亚型:
-
肿瘤来源的细胞特异性表达了胶原蛋白和细胞外基质的分子(COL10A1,COL4A1)
-
而正常组织来源的成纤维细胞中有更高的弹性蛋白表达、较低的胶原蛋白表达
-
通路分析发现其中一个亚型具有上皮-间质转换(EMT)的特征,伴随着细胞外基质蛋白marker的高表达以及TGF-β基因的表达。
-
利用SCENIC寻找每个成纤维细胞亚型特异的基因调控对,FOXO1和MSC调控的基因在Cluster1中下调,而MEF2C和ELK3调控的基因在Cluster2中上调,说明不同来源的成纤维转录因子调控对的获悉不同,进而导致不同亚型的成纤维细胞有其特异的功能。

-
B 细胞可进一步划分为9个亚型:
-
滤泡型B细胞 (MS4A1, CXCR4, HLA-DRs)
-
浆细胞(immunoglobulin gamma)
-
MALT边缘区B细胞 (immunoglobulins A,M, JCHAIN)
-
大部分的B细胞都来源于肿瘤。
-
在滤泡型B细胞亚群中,负责氧化磷酸化、细胞增殖以及生物产能的基因marker在肿瘤中表达减弱,而肿瘤相关的滤泡型B细胞转录本表达普遍更低,说明在肿瘤中滤泡B细胞出现了耗竭。
-
髓系细胞可进一步划分为12个亚群:
-
粒细胞 (S100A12+)
-
朗格汉斯细胞 (CD207+)
-
单核来源的树突状细胞 (FCGR3A+,CYTIP+)
-
交叉呈递的树突状细胞 (CLEC9A+,XCR1+)
-
CD163+和CD68+的巨噬细胞
-
大多数亚群的细胞都来源于非肿瘤组织。
-
朗格汉斯细胞簇在肿瘤和非肿瘤中数目类似。
-
通路分析在发现来源于肿瘤的巨噬细胞簇中炎症反应减弱、TNF-a诱导的细胞增殖能力不足,活性氧产生的通路下调,呈现出M2-like的巨噬细胞特征。
-
利用SCENIC,在肿瘤中发现了负责M2样巨噬细胞极化的基因调控对。
-
T细胞可进一步划分为9个亚群:
-
Treg细胞 (FOXP3+)
-
NK细胞及NKT细胞 (FGFBP2)
-
CD8+T细胞 (CD8A+)
-
CD4+T细胞 (CD4+)
-
通路分析发现肿瘤来源的T细胞和正常组织来源的T细胞类似,都表现出糖酵解功能增加和氧化磷酸化功能减弱。
-
CD8+T细胞的亚群中,表现出较强的增殖能力、移植排斥的反应,Myc活性等,且免疫检查点抑制剂的表达也较高。
-
肺泡细胞可进一步划分9个亚群:
-
I型肺泡上皮细胞 (AGER, CAV1)
-
II型肺泡上皮细胞 (SFTPC, ABCA3)
-
棒状细胞 (SCGB1A1)
-
基底细胞 (KRT15)
-
其中I型肺泡上皮细胞、II型肺泡上皮细胞,以及呼吸道表皮细胞只肿瘤中出现。
-
缺氧与是肿瘤行为的表现,包括血管生成、转移、代谢重编程以及免疫逃逸。
-
研究人员进一步按照肿瘤的部位将样本分为:肿瘤核心区域、肿瘤边缘区域、以及中间区域。发现肿瘤核心区域有更多的癌细胞。
-
关注不同部位的缺氧情况,发现肿瘤核心区域表现出最强的缺氧现象,其次是中间部位、边缘部位。
-
将不同部位的肿瘤细胞与每种类型细胞数目做线性拟合,发现大多数的基质细胞亚型富集在肿瘤边缘,并且处于常氧的环境中。

-
分析LUAD和LUSC中特异性的基质细胞亚型,发现LUAD来源于AT2的细胞而LUSC来源于基底细胞。
-
观察不同分析下,LUSC和LUAD的基质亚群特征,发现晚期的LUSC间质成分减少,而LUAD在不同分期下未发现明显规律。
-
通过TCGA数据库以及其它外部的RNASeq数据,发现在LUSC中基质marker基因的表达可以反映LUSC的预后状态。
-
I型肺泡上皮细胞、呼吸道上皮细胞以及交叉树突状细胞的marker基因高表达预示着LUSC患者生存较差。而CD8+T细胞marker基因高表达反映了LUAD患者的生存较差。