An integrated bioinformatic analysis of bulk and single-cell sequencing clarifies immune microenvironment and metabolic profiles of lung adenocarcinoma to predict immunotherapy efficacy

Frontiers in Cell and Developmental Biology

Cell and Developmental Biology

Mengling Li¹  Baosen Zhou¹  Chang Zheng² 

[1] Department of Clinical Epidemiology and Center of Evidence-Based Medicine, The First Hospital of China Medical University, Shenyang, China;null;

关键词: lung adenocarcinoma;    immunotherapy;    metabolism;    epithelia-mesenchymal transition;    cancer-associated fibroblasts;    single-cell multi-omics;   

DOI  :  10.3389/fcell.2023.1163314

 received in 2023-02-10, accepted in 2023-03-24,  发布年份 2023

来源: Frontiers

【 摘 要 】

Targeting the tumor microenvironment is increasingly recognized as an effective treatment of advanced lung adenocarcinoma (LUAD). However, few studies have addressed the efficacy of immunotherapy for LUAD. Here, a novel method for predicting immunotherapy efficacy has been proposed, which combines single-cell and bulk sequencing to characterize the immune microenvironment and metabolic profile of LUAD. TCGA bulk dataset was used to cluster two immune subtypes: C1 with “cold” tumor characteristics and C2 with “hot” tumor characteristics, with different prognosis. The Scissor algorithm, which is based on these two immune subtypes, identified GSE131907 single cell dataset into two groups of epithelial cells, labeled as Scissor_C1 and Scissor_C2. The enrichment revealed that Scissor_C1 was characterized by hypoxia, and a hypoxic microenvironment is a potential inducing factor for tumor invasion, metastasis, and immune therapy non-response. Furthermore, single cell analysis was performed to investigate the molecular mechanism of hypoxic microenvironment-induced invasion, metastasis, and immune therapy non-response in LUAD. Notably, Scissor_C1 cells significantly interacted with T cells and cancer-associated fibroblasts (CAF), and exhibited epithelial–mesenchymal transition and immunosuppressive features. CellChat analysis revealed that a hypoxic microenvironment in Scissor_C1elevated TGFβ signaling and induced ANGPTL4 and SEMA3C secretion. Interaction with endothelial cells with ANGPTL4, which increases vascular permeability and achieves distant metastasis across the vascular endothelium. Additionally, interaction of tumor-associated macrophages (TAM) and Scissor_C1 via the EREG/EFGR pathway induces tyrosine kinase inhibitor drug-resistance in patients with LAUD. Thereafter, a subgroup of CAF cells that exhibited same features as those of Scissor_C1 that exert immunosuppressive functions in the tumor microenvironment were identified. Moreover, the key genes (EPHB2 and COL1A1) in the Scissor_C1 gene network were explored and their expressions were verified using immunohistochemistry. Finally, the metabolism dysfunction in cells crosstalk was determined, which is characterized by glutamine secretion by TAM and uptake by Scissor_C1 via SLC38A2 transporter, which may induce glutamine addiction in LUAD cells. Overall, single-cell sequencing clarifies how the tumor microenvironment affects immunotherapy efficacy via molecular mechanisms and biological processes, whereas bulk sequencing explains immunotherapy efficacy based on clinical information.

【 授权许可】

Unknown   
Copyright © 2023 Li, Zhou and Zheng.