Liang Ding, Hao Shi, Koon-Kiu Yan, Yogesh Dhungana, Sivaraman Natarajan, J. Easton, J. Chiang, C. Tinkle, Xiaoyan Zhu, Liming Cai, S. Baker, H. Chi, Jiyang Yu
{"title":"Abstract 237: Inferring spatial organization of tumor microenvironment from single-cell RNA sequencing data using graph embedding","authors":"Liang Ding, Hao Shi, Koon-Kiu Yan, Yogesh Dhungana, Sivaraman Natarajan, J. Easton, J. Chiang, C. Tinkle, Xiaoyan Zhu, Liming Cai, S. Baker, H. Chi, Jiyang Yu","doi":"10.1158/1538-7445.AM2021-237","DOIUrl":null,"url":null,"abstract":"Spatial heterogeneity of diverse cellular components in the tumor microenvironment (TME) plays a critical role in the reprogramming of tumor initiation, growth, invasion, metastasis, and response to therapies. Systematic knowledge of TME spatial organization with regards to immune infiltration and tumor resource distribution is of high clinical significance. High throughput single-cell RNA sequencing (scRNA-seq) has become a revolutionary approach for studying cell composition and the development of TME. However, the spatial information of cells is lost as the tissue must be dissociated before the sequencing is performed. While various spatial techniques are emerging, their applicability is still rather limited. To address this challenge computationally, we develop a novel de novo framework to reconstruct TME spatial organization from scRNA-seq data. We hypothesized that cell spatial organization in a microenvironment is mainly determined by cell identity and interactions between different cells. In particular, the spatial organization of structural cells and immune cells follow different mechanisms. Neighboring structural cells, which share similar whole transcriptome profiles, form a scaffold of the TME; immune cells, whose activities are influenced by the structural cells, migrate in the scaffold to interact with structural cells and exert their functions. The algorithm models the scaffold of structural cells using adaptive nearest neighbor graph by taking the cell density estimation into the consideration, where the nearest neighbor graph was further augmented by inserting immune cells into the appropriate locations of the scaffold according to the LR similarities. To reconstruct 3D spatial organization while preserving the cell topology represented by the graph, we employed a graph embedding strategy to minimize the discrepancy between the graph topology and the embedded 3D space. We evaluated the framework on two diffuse intrinsic pontine gliomas (DIPG) samples from a mouse model with coupled scRNA-seq and spatial transcriptome (ST, 10x Visium platform) data. The predicted spatial organization successfully separates the major cell types. The T-cell infiltrated tumor, verified by the T-cell spatial spots of the ST image, is well recapitulated. We deconvoluted the ST data by integrating the scRNA-seq data using SPOTlight. The neighborhood enrichment distributions of predicted spatial organization and the spot deconvoluted ST data show high consistency as measured by Kullback-Leibler divergence. We found heterogeneous neighborhood composition of CD8+ T-cells, indicating diverse clonality and functions with respect to their locations in the TME. Citation Format: Liang Ding, Hao Shi, Koon-Kiu Yan, Yogesh Dhungana, Sivaraman Natarajan, John Easton, Jason Chiang, Christopher L. Tinkle, Xiaoyan Zhu, Liming Cai, Suzanne J. Baker, Hongbo Chi, Jiyang Yu. Inferring spatial organization of tumor microenvironment from single-cell RNA sequencing data using graph embedding [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 237.","PeriodicalId":73617,"journal":{"name":"Journal of bioinformatics and systems biology : Open access","volume":"26 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2021-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of bioinformatics and systems biology : Open access","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1158/1538-7445.AM2021-237","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Spatial heterogeneity of diverse cellular components in the tumor microenvironment (TME) plays a critical role in the reprogramming of tumor initiation, growth, invasion, metastasis, and response to therapies. Systematic knowledge of TME spatial organization with regards to immune infiltration and tumor resource distribution is of high clinical significance. High throughput single-cell RNA sequencing (scRNA-seq) has become a revolutionary approach for studying cell composition and the development of TME. However, the spatial information of cells is lost as the tissue must be dissociated before the sequencing is performed. While various spatial techniques are emerging, their applicability is still rather limited. To address this challenge computationally, we develop a novel de novo framework to reconstruct TME spatial organization from scRNA-seq data. We hypothesized that cell spatial organization in a microenvironment is mainly determined by cell identity and interactions between different cells. In particular, the spatial organization of structural cells and immune cells follow different mechanisms. Neighboring structural cells, which share similar whole transcriptome profiles, form a scaffold of the TME; immune cells, whose activities are influenced by the structural cells, migrate in the scaffold to interact with structural cells and exert their functions. The algorithm models the scaffold of structural cells using adaptive nearest neighbor graph by taking the cell density estimation into the consideration, where the nearest neighbor graph was further augmented by inserting immune cells into the appropriate locations of the scaffold according to the LR similarities. To reconstruct 3D spatial organization while preserving the cell topology represented by the graph, we employed a graph embedding strategy to minimize the discrepancy between the graph topology and the embedded 3D space. We evaluated the framework on two diffuse intrinsic pontine gliomas (DIPG) samples from a mouse model with coupled scRNA-seq and spatial transcriptome (ST, 10x Visium platform) data. The predicted spatial organization successfully separates the major cell types. The T-cell infiltrated tumor, verified by the T-cell spatial spots of the ST image, is well recapitulated. We deconvoluted the ST data by integrating the scRNA-seq data using SPOTlight. The neighborhood enrichment distributions of predicted spatial organization and the spot deconvoluted ST data show high consistency as measured by Kullback-Leibler divergence. We found heterogeneous neighborhood composition of CD8+ T-cells, indicating diverse clonality and functions with respect to their locations in the TME. Citation Format: Liang Ding, Hao Shi, Koon-Kiu Yan, Yogesh Dhungana, Sivaraman Natarajan, John Easton, Jason Chiang, Christopher L. Tinkle, Xiaoyan Zhu, Liming Cai, Suzanne J. Baker, Hongbo Chi, Jiyang Yu. Inferring spatial organization of tumor microenvironment from single-cell RNA sequencing data using graph embedding [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 237.
肿瘤微环境(tumor microenvironment, TME)中不同细胞组分的空间异质性在肿瘤发生、生长、侵袭、转移和治疗反应的重编程中起着至关重要的作用。系统了解TME的空间组织与免疫浸润和肿瘤资源分布的关系具有重要的临床意义。高通量单细胞RNA测序(scRNA-seq)已成为研究细胞组成和TME发展的革命性方法。然而,细胞的空间信息丢失,因为组织必须在测序之前解离。虽然各种空间技术不断涌现,但其适用性仍然相当有限。为了在计算上解决这一挑战,我们开发了一个新的从头框架,从scRNA-seq数据中重建TME空间组织。我们假设微环境中的细胞空间组织主要由细胞身份和不同细胞之间的相互作用决定。特别是结构细胞和免疫细胞的空间组织遵循不同的机制。邻近的结构细胞,具有相似的整个转录组谱,形成TME的支架;免疫细胞的活动受到结构细胞的影响,在支架内迁移,与结构细胞相互作用,发挥其功能。该算法在考虑细胞密度估计的基础上,采用自适应最近邻图对结构细胞的支架进行建模,并根据LR相似性在支架的适当位置插入免疫细胞,进一步增强最近邻图。为了在保留图所表示的细胞拓扑的同时重建三维空间组织,我们采用了图嵌入策略来最小化图拓扑与嵌入的三维空间之间的差异。我们利用scRNA-seq和空间转录组(ST, 10倍Visium平台)数据对来自小鼠模型的两个弥漫性内生性脑桥胶质瘤(DIPG)样本进行了框架评估。预测的空间组织成功地分离了主要的细胞类型。t细胞浸润的肿瘤,通过ST图像的t细胞空间斑点验证,被很好地再现。我们通过使用SPOTlight整合scRNA-seq数据对ST数据进行反卷积。通过Kullback-Leibler散度测量,预测空间组织的邻域富集分布与点反卷积的ST数据具有较高的一致性。我们发现CD8+ t细胞的异质邻域组成,表明它们在TME中的位置具有不同的克隆性和功能。引用格式:丁亮,史浩,严坤桥,Yogesh Dhungana, Sivaraman Natarajan, John Easton, Jason Chiang, Christopher L. Tinkle,朱晓燕,蔡黎明,Suzanne J. Baker,迟洪波,余纪阳。利用图嵌入从单细胞RNA测序数据推断肿瘤微环境的空间组织[摘要]。见:美国癌症研究协会2021年年会论文集;2021年4月10日至15日和5月17日至21日。费城(PA): AACR;癌症杂志,2021;81(13 -增刊):摘要第237期。