Characterization of tumor-associated reactive astrocytes in gliomas by single-cell and bulk tumor sequencing

Chuan-Bao Zhang; Zhi-Liang Wang; Han-Jie Liu; Zheng Wang; Wang Jia

doi:10.3389/fneur.2023.1193844

Characterization of tumor-associated reactive astrocytes in gliomas by single-cell and bulk tumor sequencing

Front Neurol. 2023 Jun 21:14:1193844. doi: 10.3389/fneur.2023.1193844. eCollection 2023.

Authors

Chuan-Bao Zhang^{1

2}, Zhi-Liang Wang^{1

2}, Han-Jie Liu^{1

2}, Zheng Wang^{1

2}, Wang Jia^{1

2}

Affiliations

¹ Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China.
² The Chinese Glioma Genome Atlas (CGGA) Project, Beijing, China.

Abstract

Objective: Astrocytes constitute approximately 30% of cells in gliomas and play important roles in synapse construction and survival. Recently, JAK/STAT pathway activation associated with a new type of astrocyte was reported. However, the implications of these tumor-associated reactive astrocytes (TARAs) in glioma are not known.

Methods: We comprehensively assessed TARAs in gliomas, both in single cells and at the bulk tumor level, by analyzing five independent datasets. First, we analyzed two single-cell RNA sequencing datasets of 35,563 cells from 23 patients to estimate the infiltration level of TARAs in gliomas. Second, we collected clinical information and genomic and transcriptomic data of 1,379 diffuse astrocytoma and glioblastoma samples from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas datasets to evaluate the genomic, transcriptomic and clinical characteristics of TARA infiltration. Third, we downloaded expression profiles of recurrent glioblastoma samples from patients receiving PD-1 inhibitors to analyze the predictive value of TARAs for immune checkpoint inhibition.

Results: Single-cell RNA sequencing data showed TARAs were abundant in the glioma micro-environment (15.7% in the CGGA dataset and 9.1% in the Gene Expression Omnibus GSE141383 dataset, respectively). Bulk tumor sequencing data showed that the extent of TARA infiltration was highly associated with major clinical and molecular features of astrocytic gliomas. Patients with more TARA infiltration were more likely to have MUC16, FLG, and PICK3A mutations, chromosome 9p21.3, 10q23.3, and 13q14.2 deletions and 7p11.2 amplification. Gene Ontology analysis revealed that the high level of astrocyte infiltration was characterized by immune and oncogenic pathways, such as the inflammatory response, positive regulation of the JAK-STAT cascade, positive regulation of NIK/NF-kappa B signaling and the tumor necrosis factor biosynthetic process. Patients with greater TARA infiltration showed inferior prognosis. Meanwhile, the extent of reactive astrocyte infiltration exhibited a predictive value for recurrent glioblastoma patients undergoing anti-PD-1 immune therapy.

Conclusion: TARA infiltration might promote glioma tumor progression and can be used as a diagnostic, predictive and prognostic marker in gliomas. Prevention of TARA infiltration might be a new therapeutic strategy for glioma.

Keywords: glioma; immunotherapy; molecular and genomic features; reactive astrocytes; sequencing.