基于多转录组测序的胶质母细胞瘤间质样肿瘤亚群缺氧特征及预后研究

陈涛; 杨帆; 王帅; 罗敏; 颜泽萱; 唐采蝶; 甯芸; 杨思思; 曹若菲; 李政波; 方煊宇; 姚小红

doi:10.16016/j.2097-0927.202504111

您当前的位置：

首页 >

文章列表页 >

基于多转录组测序的胶质母细胞瘤间质样肿瘤亚群缺氧特征及预后研究

基础医学 | 更新时间：2025-08-25

- 基于多转录组测序的胶质母细胞瘤间质样肿瘤亚群缺氧特征及预后研究
- Hypoxia characteristics and prognostic significance of glioblastoma MES-like subpopulation through multi-transcriptomics sequencing
- 陆军军医大学学报 2025年47卷第16期页码：1904-1912
- 作者机构：
  
  陆军军医大学（第三军医大学）第一附属医院病理学研究所，肿瘤免疫病理学教育部重点实验室，重庆
- 作者简介：
  
  姚小红，E-mail：yaoxiaohong@tmmu.edu.cn
- 基金信息：
  
  国家自然科学基金面上项目(82273082);重庆市技术创新与应用发展专项重点项目(cstc2021jscx-cylhX0005)
- DOI：10.16016/j.2097-0927.202504111
  中图分类号： R394-33;R730.23;R730.264
- 收稿日期：2025-04-27，
  
  修回日期：2025-06-26，
  
  纸质出版日期：2025-08-30
- 稿件说明：
移动端阅览
陈涛，杨帆，王帅，等.基于多转录组测序的胶质母细胞瘤间质样肿瘤亚群缺氧特征及预后研究［J］.陆军军医大学学报，2025，47（16）：1904-1912.DOI：10.16016/j.2097-0927.202504111.

CHEN T，YANG F，WANG S，et al.Hypoxia characteristics and prognostic significance of glioblastoma MES-like subpopulation through multi-transcriptomics sequencing［J］.J Army Med Univ，2025，47（16）：1904-1912.DOI：10.16016/j.2097-0927.202504111.
陈涛，杨帆，王帅，等.基于多转录组测序的胶质母细胞瘤间质样肿瘤亚群缺氧特征及预后研究［J］.陆军军医大学学报，2025，47（16）：1904-1912.DOI：10.16016/j.2097-0927.202504111. DOI：

CHEN T，YANG F，WANG S，et al.Hypoxia characteristics and prognostic significance of glioblastoma MES-like subpopulation through multi-transcriptomics sequencing［J］.J Army Med Univ，2025，47（16）：1904-1912.DOI：10.16016/j.2097-0927.202504111. DOI：

摘要

目的

通过多转录组测序数据分析，解析胶质母细胞瘤（glioblastoma， GBM）的肿瘤细胞异质性组成特点和生物学功能，以探究GBM间质样（mesenchymal-like， MES-like）肿瘤亚群的缺氧特点、空间定位特征以及对GBM恶性生物学行为的影响。

方法

分析公共数据库18例GBM单细胞转录组测序（single-cell RNA sequencing， scRNA-seq）数据，以及GBM普通转录组测序（bulk RNA sequencing，bulk RNA-seq）和空间转录组测序（spatial transcriptomics，ST）数据，通过单细胞转录组定义GBM肿瘤亚群，基因本体论（Gene Ontology， GO）和基因集富集分析（gene set enrichment analysis， GSEA）进行功能富集。使用BayesPrism反卷积方法评估普通转录组的细胞亚群占比和定位，并收集2015年至2023年就诊于西南医院神经外科患者的12例符合胶质母细胞瘤病理学诊断标准的石蜡样本（男性10例，女性2例，平均年龄为53.50岁，中位年龄为54.50岁）进行免疫荧光验证。pySCENIC预测肿瘤细胞亚群特异转录因子。

结果

GBM肿瘤细胞具有高度异质性，主要分为4种亚群，即星形胶质细胞样（astrocyte-like，AC-like）、神经前体细胞样（neural progenitor-like，NPC-like）、少突前体细胞样（oligodendrocyte progenitor-like，OPC-like）和MES-like肿瘤亚群。差异基因分析发现MES-like肿瘤亚群高表达血管内皮生长因子A（vascular endothelial growth Factor A，VEGFA）、肾上腺髓质素（adrenomedullin， ADM）、N-myc下游调控基因1（N-myc downstream regulated 1，NDRG1）、胰岛素样生长因子结合蛋白 5（insulin like growth factor binding protein 5，IGFBP5）、A-激酶锚定蛋白 12（A-kinase anchoring protein 12，AKAP12）（

0.001）。pySCENIC转录因子预测发现，MES-like肿瘤亚群的高活性转录因子为AT-rich 交互域 3A（AT-rich interaction domain 3A，ARID3A）、FOS like 2， AP-1 转录因子亚基（FOS Like 2， AP-1 transcription factor subunit，FOSL2）、内皮 PAS 结构域蛋白 1（end

othelial PAS domain protein 1，EPAS1）、CCAAT增强子结合蛋白Delta（CCAAT enhancer binding protein delta，CEBPD）、CCAAT 增强子结合蛋白β（CCAAT enhancer binding protein beta，CEBPB）（

0.05）。GO和GSEA富集发现，MES-like肿瘤亚群富集于缺氧相关通路，如细胞对低氧水平的反应（NES=2.437，

0.001）。BayesPrism反卷积显示，MES-like肿瘤亚群主要存在于假栅栏样坏死区和坏死周围。GBM石蜡样本免疫荧光证实CD44

（CD44 antigen，CD44抗原）MES-like肿瘤亚群主要存在于HIF1α

high

（Hypoxia Inducible Factor 1 Subunit Alpha，缺氧诱导因子1α）的缺氧区域（

0.01）。多因素Cox回归显示MES-like肿瘤亚群与患者不良预后显著相关（

=1.71，95%

：1.38~2.11，

0.001）。

结论

GBM肿瘤细胞具有高度异质性，主要分为AC-like、NPC-like、OPC-like和MES-like肿瘤亚群，其MES-like肿瘤亚群具有缺氧特点，主要定位于假栅栏样坏死区和坏死周围，具有促进GBM恶性生物学行为的作用。

Abstract

Objective‍　‍To analyze the composition characteristics and biological functions of tumor cell subpopulations in glioblastoma (GBM) through multi-transcriptomics sequencing technology

and explore the hypoxia characteristics and spatial localization features of the mesenchymal-like (MES-like) tumor cell subpopulation in GBM and the influence on malignant biological behaviors. Methods‍　‍Multi-transcriptomics sequencing data

including single-cell RNA sequencing (scRNA-seq) data (18 patients)

bulk RNA sequencing (bulk RNA-seq) and spatial transcriptomics (ST) data of GBM

were employed to define cell subpopulations in GBM

and Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) were utilized to analyze their functions. The proportions and locations of cell subpopulations in bulk RNA-seq data were evaluated with BayesPrism deconvolution. Immunofluorescence assay was conducted for verification on 12 paraffin samples of GBM from patients who visited the neurosurgical department of our hospital from 2015 to 2023 and met the pathological diagnostic criteria for GBM (10 males and 2 females

at an average age of 53.50 years and a median age of 54.50 years). pySCENIC was applied to predict sp

ecific transcription factors of tumor cell subpopulations. Results‍　‍Tumor cells in GBM were highly heterogeneous

and could be mainly divided into 4 subpopulations: astrocyte-like (AC-like)

neural progenitor-like (NPC-like)

oligodendrocyte progenitor-like (OPC-like) and MES-like. Differential gene analysis found that the MES-like tumor cells highly expressed vascular endothelial growth factor A (VEGFA)

adrenomedullin (ADM)

N-myc downstream regulated 1 (NDRG1)

insulin like growth factor binding protein 5 (IGFBP5)

and A-kinase anchoring protein 12 (AKAP12) (

0.001). pySCENIC transcription factor prediction found that the high-active transcription factors of the MES-like tumor cells were AT-rich interaction domain 3A (ARID3A)

FOS like 2

AP-1 transcription factor subunit (FOSL2)

endothelial PAS domain protein 1 (EPAS1)

CCAAT enhancer binding protein delta (CEBPD)

and CCAAT enhancer binding protein beta (CEBPB) (

0.05). GO and GSEA enrichment analyses found that the MES-like tumor cells were enriched in hypoxia-related pathways

especially the pathway of cell responses to hypoxia levels (NES=2.437

0.001). BayesPrism deconvolution showed that the MES-like tumor cells mainly existed in PAN(Pseudopalisading cells around necrosis) and perinecrotic zone. Immunofluorescence assay confirmed CD44

(CD44 antigen) MES-like tumor cells were mainly located in hypoxia areas with highly expression of hypoxia inducible factor 1 subunit alpha (HIF1α) (

0.01). Multivariate Cox regression analysis indicated that the MES-like tumor cells were significantly correlated with the adverse prognosis of GBM patients (

=1.71

95%

: 1.38~2.11

0.001). Conclusion‍　‍Tumor cells in GBM are of highly heterogeneity. They could be mainly divided into 4 subpopulations: AC-like

NPC-like

OPC-like and MES-like. MES-like tumor cells

mainly locating in PAN and perinecrotic zone

are charact

erized by hypoxia

which can promote the malignant progression of GBM.

关键词

Keywords

references

TAN A C ， ASHLEY D M ， LÓPEZ G Y ， et al . Management of glioblastoma： State of the art and future directions ［J］. CA Cancer J Clin ， 2020 ， 70 （ 4 ）： 299 - 312 . DOI： 10.3322/caac.21613 http://dx.doi.org/10.3322/caac.21613 .

国家卫生健康委员会医政医管局，中国抗癌协会脑胶质瘤专业委员会，中国医师协会脑胶质瘤专业委员会，等 . 脑胶质瘤诊疗指南（2022版）［J］. 中华神经外科杂志， 2022 ， 38 （ 8 ）： 757 - 777 . DOI： 10.3760/cma.j.cn112050-20220510-00239 http://dx.doi.org/10.3760/cma.j.cn112050-20220510-00239 .

Bureau of Medical Administration of National Health Commission of the People’s Republic of China ， et al . Guidelines for diagnosis and treatment of glioma （2022 edition）［J］. Chin J Neurosurg ， 2022 ， 38 （ 8 ）： 757 - 777 . DOI： 10.3760/cma.j.cn112050-20220510-00239 http://dx.doi.org/10.3760/cma.j.cn112050-20220510-00239 .

PARADA L F ， DIRKS P B ， WECHSLER-REYA R J . Brain tumor stem cells remain in play ［J］. J Clin Oncol ， 2017 ， 35 （ 21 ）： 2428 - 2431 . DOI： 10.1200/JCO.2017.73.9540 http://dx.doi.org/10.1200/JCO.2017.73.9540 .

DU R F ， SANIN A Y ， SHI W J ， et al . Muscarinic receptor drug trihexyphenidyl can alter growth of mesenchymal glioblastoma in vivo ［J］. Front Pharmacol ， 2024 ， 15 ： 1468920 . DOI： 10.3389/fphar.2024.1468920 http://dx.doi.org/10.3389/fphar.2024.1468920 .

ZHAO J X ， CUI X T ， ZHAN Q ， et al . CRISPR-Cas9 library screening combined with an exosome-targeted delivery system addresses tumorigenesis/TMZ resistance in the mesenchymal subtype of glioblastoma ［J］. Theranostics ， 2024 ， 14 （ 7 ）： 2835 - 2855 . DOI： 10.7150/thno.92703 http://dx.doi.org/10.7150/thno.92703 .

Hong M ， Tao S ， Zhang L ， et al . RNA sequencing： new technologies and applications in cancer research ［J］. J Hematol Oncol ， 2020 ， 13 （ 1 ）： 166 . DOI： 10.1186/s13045-020-01005-x http://dx.doi.org/10.1186/s13045-020-01005-x .

BOXER E ， FEIGIN N ， TSCHERNICHOVSKY R ， et al . Emerging clinical applications of single-cell RNA sequencing in oncology ［J］. Nat Rev Clin Oncol ， 2025 ， 22 （ 5 ）： 315 - 326 . DOI： 10.1038/s41571-025-01003-3 http://dx.doi.org/10.1038/s41571-025-01003-3 .

ZORMPAS E ， QUEEN R ， COMBER A ， et al . Mapping the transcriptome： Realizing the full potential of spatial data analysis ［J］. Cell ， 2023 ， 186 （ 26 ）： 5677 - 5689 . DOI： 10.1016/j.cell.2023.11.003 http://dx.doi.org/10.1016/j.cell.2023.11.003 .

WEINSTEIN J N ， COLLISSON E A ， MILLS G B ， et al . The cancer genome atlas pan-cancer analysis project ［J］. Nat Genet ， 2013 ， 45 （ 10 ）： 1113 - 1120 . DOI： 10.1038/ng.2764 http://dx.doi.org/10.1038/ng.2764 .

ZHAO Z ， ZHANG K N ， WANG Q W ， et al . Chinese glioma genome atlas （CGGA）： a comprehensive resource with functional genomic data from Chinese glioma patients ［J］. Genomics Proteomics Bioinformatics ， 2021 ， 19 （ 1 ）： 1 - 12 . DOI： 10.1016/j.gpb.2020.10.005 http://dx.doi.org/10.1016/j.gpb.2020.10.005 .

PUCHALSKI R B ， SHAH N ， MILLER J ， et al . An anatomic transcriptional atlas of human glioblastoma ［J］. Science ， 2018 ， 360 （ 6389 ）： 660 - 663 . DOI： 10.1126/science.aaf2666 http://dx.doi.org/10.1126/science.aaf2666 .

KUECKELHAUS J ， FRERICH S ， KADA-BENOTMANE J ， et al . Inferring histology-associated gene expression gradients in spatial transcriptomic studies ［J］. Nat Commun ， 2024 ， 15 （ 1 ）： 7280 . DOI： 10.1038/s41467-024-50904-x http://dx.doi.org/10.1038/s41467-024-50904-x .

CHU T ， WANG Z ， PE’ER D ， et al . Cell type and gene expression deconvolution with BayesPrism enables Bayesian integrative analysis across bulk and single-cell RNA sequencing in oncology ［J］. Nat Cancer ， 2022 ， 3 （ 4 ）： 505 - 517 . DOI： 10.1038/s43018-022-00356-3 http://dx.doi.org/10.1038/s43018-022-00356-3 .

WANG W Y ， LI T R ， CHENG Y ， et al . Identification of hypoxic macrophages in glioblastoma with therapeutic potential for vasculature normalization ［J］. Cancer Cell ， 2024 ， 42 （ 5 ）： 815 - 832.e12 . DOI： 10.1016/j.ccell.2024.03.013 http://dx.doi.org/10.1016/j.ccell.2024.03.013 .

NEFTEL C ， LAFFY J ， FILBIN M G ， et al . An integrative model of cellular states， plasticity， and genetics for glioblastoma ［J］. Cell ， 2019 ， 178 （ 4 ）： 835 - 849.e21 . DOI： 10.1016/j.cell.2019.06.024 http://dx.doi.org/10.1016/j.cell.2019.06.024 .

CHEN L M ， XIE X S ， WANG T T ， et al . ARL13B promotes angiogenesis and glioma growth by activating VEGFA-VEGFR2 signaling ［J］. Neuro Oncol ， 2023 ， 25 （ 5 ）： 871 - 885 . DOI： 10.1093/neuonc/noac245 http://dx.doi.org/10.1093/neuonc/noac245 .

SUBRAMANIAN A ， TAMAYO P ， MOOTHA V K ， et al . Gene set enrichment analysis： a knowledge-based approach for interpreting genome-wide expression profiles ［J］. Proc Natl Acad Sci USA ， 2005 ， 102 （ 43 ）： 15545 - 15550 . DOI： 10.1073/pnas.0506580102 http://dx.doi.org/10.1073/pnas.0506580102 .

VAN DE SANDE B ， FLERIN C ， DAVIE K ， et al . A scalable SCENIC workflow for single-cell gene regulatory network analysis ［J］. Nat Protoc ， 2020 ， 15 （ 7 ）： 2247 - 2276 . DOI： 10.1038/s41596-020-0336-2 http://dx.doi.org/10.1038/s41596-020-0336-2 .

HUANG H D ， SHAH H ， HAO J ， et al . Long non-coding RNA lung cancer-associated transcript-1 promotes glioblastoma progression by enhancing Hypoxia-inducible factor 1 alpha activity ［J］. Neuro Oncol ， 2024 ， 26 （ 8 ）： 1388 - 1401 . DOI： 10.1093/neuonc/noae036 http://dx.doi.org/10.1093/neuonc/noae036 .

LIU D ， ZHU H T ， CHENG L D ， et al . Hypoxia-induced galectin-8 maintains stemness in glioma stem cells via autophagy regulation ［J］. Neuro Oncol ， 2024 ， 26 （ 5 ）： 872 - 888 . DOI： 10.1093/neuonc/noad264 http://dx.doi.org/10.1093/neuonc/noad264 .

PIETRAS A ， KATZ A M ， EKSTRÖM E J ， et al . Osteopontin-CD44 signaling in the glioma perivascular niche enhances cancer stem cell phenotypes and promotes aggressive tumor growth ［J］. Cell Stem Cell ， 2014 ， 14 （ 3 ）： 357 - 369 . DOI： 10.1016/j.stem.2014.01.005 http://dx.doi.org/10.1016/j.stem.2014.01.005 .

MAO X G ， XUE X Y ， WANG L ， et al . CDH5 is specifically activated in glioblastoma stemlike cells and contributes to vasculogenic mimicry induced by hypoxia ［J］. Neuro Oncol ， 2013 ， 15 （ 7 ）： 865 - 879 . DOI： 10.1093/neuonc/not029 http://dx.doi.org/10.1093/neuonc/not029 .

KANG W ， MO Z X ， LI W S ， et al . Heterogeneity and individualized treatment of microenvironment in glioblastoma （Review）［J］. Oncol Rep ， 2023 ， 50 （ 6 ）： 217 . DOI： 10.3892/or.2023.8654 http://dx.doi.org/10.3892/or.2023.8654 .

LI J B ， SHIH L K ， MARKWELL S M ， et al . A clinically relevant model and method to study necrosis as a driving force in glioma restructuring and progression ［J］. Proc Natl Acad Sci USA ， 2025 ， 122 （ 7 ）： e2416024122 . DOI： 10.1073/pnas.2416024122 http://dx.doi.org/10.1073/pnas.2416024122 .

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

缺氧微环境下HIF-1α与circ-UBE2G1相互调控促进甲状腺癌转移

胶质瘤内肿瘤相关巨噬细胞的缺氧转录表型及细胞超微结构改变

缺氧通过TAZ-H2A.Z轴调控H2A.Z染色质沉积促进肺泡上皮细胞增殖的研究

缺氧特异性诱导肺血管内皮细胞黏附分子3表达以促进单核细胞黏附

SOX6-SOCS3通过STAT3信号通路调控缺氧诱导的人肺动脉平滑肌细胞增殖