Comprehensive analysis reveals a prognostic and therapeutic biomarker CD3D in the breast carcinoma microenvironment

Breast carcinoma (BRCA) is the most common carcinoma among women worldwide. Despite the great progress achieved in early detection and treatment, morbidity and mortality rates remain high. In the present study, we make a systematic analysis of BRCA using TCGA database by applying CIBERSORT and ESTIMATE computational methods, uncovered CD3D as a prognostic biomarker by intersection analysis of univariate COX and protein–protein interaction (PPI). It revealed that high CD3D expression was strongly associated with poor survival of BRCA, based on The Cancer Genome Atlas (TCGA) database and online websites. Gene Set Enrichment Analysis (GSEA) revealed that the high CD3D expression group was mainly enriched for the immune-related pathways and the low CD3D expression group was mainly enriched for metabolic-related activities. Based on CIBERSORT analysis, the difference test and correlation test suggested that CD3D had a strong correlation with T cells, particularly CD8 + T cells, which indicated that CD3D up-regulation may increase T cell immune infiltration in the TME and induce antitumor immunity by activating T lymphocytes. Furthermore, the correlation analysis showed that CD3D expression had a strongly positive correlation with immune checkpoints, which indicating that the underlying mechanism involves CD3D mediated regulation of T cell functions in BRCA, and single cell RNA-seq analysis revealed that CD3D correlate with CD8 + T cells and it is itself highly expressed in CD8 + T cells. In summary, we identified a prognostic biomarker CD3D in BRCA, which was associated with lymphocyte infiltration, immune checkpoints and could be developed for innovative therapeutics of BRCA.     

Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas

The tumor microenvironment is highly correlated with tumor occurrence, progress, and prognosis. We aimed to investigate the immune-related gene (IRG) expression and immune infiltration pattern in the tumor microenvironment of lower-grade glioma (LGG). We employed the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate immune and stromal scores and identify prognostic IRG based on The Cancer Genome Atlas data set. The potential molecular functions of these genes were explored with the help of functional enrichment analysis and the protein–protein interaction network. Remarkably, three cohorts that were downloaded from the Chinese Glioma Genome Atlas database were analyzed to further verify the prognostic values of these genes. Moreover, the Tumor IMmune Estimation Resource (TIMER) algorithm was used to estimate the abundance of infiltrating immune cells and explore the immune infiltration pattern in LGG. And unsupervised cluster analysis determined three clusters of the immune infiltration pattern and indicated that CD8⁺ T cells and macrophages were significantly associated with LGG outcomes. Altogether, our study identified a list of prognostic IRGs and provided a perspective to explore the immune infiltration pattern in LGG.     

Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of cutaneous melanoma

Cutaneous malignant melanoma (hereafter called melanoma) is one of the most aggressive cancers with increasing incidence and mortality rates worldwide. In this study, we performed a systematic investigation of the tumor microenvironmental and genetic factors associated with melanoma to identify prognostic biomarkers for melanoma. We calculated the immune and stromal scores of melanoma patients from the Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm and found that they were closely associated with patients’ prognosis. Then the differentially expressed genes were obtained based on the immune and stromal scores, and prognostic immune-related genes further identified. Functional analysis and the protein–protein interaction network further revealed that these genes enriched in many immune-related biological processes. In addition, the abundance of six infiltrating immune cells was analyzed using prognostic immune-related genes by TIMER algorithm. The unsupervised clustering analysis using immune-cell proportions revealed eight clusters with distinct survival patterns, suggesting that dendritic cells were most abundant in the microenvironment and CD8⁺ T cells and neutrophils were significantly related to patients’ prognosis. Finally, we validated these genes in three independent cohorts from the Gene Expression Omnibus database. In conclusion, this study comprehensively analyzed the tumor microenvironment and identified prognostic immune-related biomarkers for melanoma.     

Activated breast stromal fibroblasts exhibit myoepithelial and mammary stem cells features

BackgroundActive breast cancer-associated fibroblasts (CAFs) promote tumor growth and spread, and like tumor cells they are also heterogeneous with various molecular sub-types and different pro-tumorigenic capacities.MethodsWe have used immunoblotting as well as quantitative RT-PCR to assess the expression of various epithelial/mesenchymal as well as stemness markers in breast stromal fibroblasts. Immunofluorescence was utilized to assess the level of different myoepithelial and luminal markers at the cellular level. Flow cytometry allowed to determine the proportion of CD44- and ALDH1-positive breast fibroblasts, while sphere formation assay was used to test the ability of these cells to form mammospheres.ResultsWe have shown here that IL-6-dependent activation of breast and skin fibroblasts promotes mesenchymal-to-epithelial transition and stemness in a STAT3- and p16-dependent manner. Interestingly, most primary CAFs isolated from breast cancer patients exhibited such transition and expressed lower levels of the mesenchymal markers N-cadherin and vimentin as compared to their adjacent normal fibroblasts (TCFs) isolated from the same patients. We have also shown that some CAFs and IL-6-activated fibroblasts express high levels of the myoepithelial markers cytokeratin 14 and CD10. Interestingly, 12 CAFs isolated from breast tumors showed higher proportions of CD24^low/CD44^high and ALDH^high cells, compared to their corresponding TCF cells. These CD44^high cells have higher abilities to form mammospheres and to enhance cell proliferation of breast cancer cells in a paracrine manner relative to their corresponding CD44^low cells.ConclusionTogether, the present findings show novel characteristics of active breast stromal fibroblasts, which exhibit additional myoepithelial/progenitor features.     

Reprogramming of cancer-associated fibroblasts combined with immune checkpoint inhibitors: A potential therapeutic strategy for cancers

Activated fibroblasts, namely cancer-associated fibroblasts (CAFs), are highly heterogeneous in phenotypes, functions, and origins. CAFs originated from varieties of cell types, including local resident fibroblasts, epithelial cells, mesenchymal stromal cells, or others. These cells participate in tumor angiogenesis, mechanics, drug access, and immune suppression, with the latter being particularly important. It was difficult to distinguish CAFs by subsets due to their complex origins until the use of scRNA-seq. Reprogramming CAFs with TGFβ-RI inhibitor, a CXCR4 blocker, or other methods increases T cells activation and infiltration, together with a decrease in CAFs recruitment, thus improving the prognosis. As depletion of CAFs can't bring clinical benefit, the combination of reprogramming CAFs and immune checkpoint inhibitors (ICIs) come into consideration. It has shown better outcomes compared with monotherapy respectively in basic/preclinical researches, and needs more data on clinical trials. Combination therapy may be a promising and expecting method for treatment of cancer.     

Biological Characteristics and Genetic Heterogeneity between Carcinoma-Associated Fibroblasts and Their Paired Normal Fibroblasts in Human Breast Cancer

Background

The extensional signals in cross-talk between stromal cells and tumor cells generated from extracellular matrix molecules, soluble factor, and cell-cell adhesion complexes cooperate at the extra- and intracellular level in the tumor microenvironment. CAFs are the primary type of stromal cells in the tumor microenvironment and play a pivotal role in tumorigenesis and development. Hitherto, there is hardly any systematic analysis of the intrinsic relationship between CAFs function and its abnormal signaling pathway. The extreme complexity of CAFs’ features and their role in tumor development are needed to be further investigated.

Methodology/Principal Findings

We primary cultured CAFs and NFs from early stages of breast cancer tissue and identified them using their biomarker by immunohistochemistry for Fibronectin, α-SMA and FAP. Microarray was applied to analyze gene expression profiles of human breast CAFs and the paired NFs. The Up-regulated genes classified by Gene Ontology, signal pathways enriched by DAVID pathway analysis. Abnormal signaling pathways in breast cancer CAFs are involved in cell cycle, cell adhesion, signal transduction and protein transport being reported in CAFs derived from other tumors. Significantly, the altered ATM signaling pathway, a set of cell cycle regulated signaling, and immune associated signaling are identified to be changed in CAFs.

Conclusions/Significance

CAFs have the vigorous ability of proliferation and potential of invasion and migration comparing with NFs. CAFs could promote breast cancer cell invasion under co-culture conditions through up-regulated CCL18 and CXCL12. Consistently with its biologic behavior, the gene expression profiling analyzed by microarray shows that some of key signaling pathways, such as cell cycle, cell adhesion, and secreting factors play an important role in CAFs. The altered ATM signaling pathway is abnormally active in the early stage of breast cancer. The set of immune associated signaling may be involved in tumor cell immune evasion.     

Immune and Inflammatory Cell Composition of Human Lung Cancer Stroma

Recent studies indicate that the abnormal microenvironment of tumors may play a critical role in carcinogenesis, including lung cancer. We comprehensively assessed the number of stromal cells, especially immune/inflammatory cells, in lung cancer and evaluated their infiltration in cancers of different stages, types and metastatic characteristics potential. Immunohistochemical analysis of lung cancer tissue arrays containing normal and lung cancer sections was performed. This analysis was combined with cyto-/histomorphological assessment and quantification of cells to classify/subclassify tumors accurately and to perform a high throughput analysis of stromal cell composition in different types of lung cancer. In human lung cancer sections we observed a significant elevation/infiltration of total-T lymphocytes (CD3⁺), cytotoxic-T cells (CD8⁺), T-helper cells (CD4⁺), B cells (CD20⁺), macrophages (CD68⁺), mast cells (CD117⁺), mononuclear cells (CD11c⁺), plasma cells, activated-T cells (MUM1⁺), B cells, myeloid cells (PD1⁺) and neutrophilic granulocytes (myeloperoxidase⁺) compared with healthy donor specimens. We observed all of these immune cell markers in different types of lung cancers including squamous cell carcinoma, adenocarcinoma, adenosquamous cell carcinoma, small cell carcinoma, papillary adenocarcinoma, metastatic adenocarcinoma, and bronchioloalveolar carcinoma. The numbers of all tumor-associated immune cells (except MUM1⁺ cells) in stage III cancer specimens was significantly greater than those in stage I samples. We observed substantial stage-dependent immune cell infiltration in human lung tumors suggesting that the tumor microenvironment plays a critical role during lung carcinogenesis. Strategies for therapeutic interference with lung cancer microenvironment should consider the complexity of its immune cell composition.     

CA9 transcriptional expression determines prognosis and tumour grade in tongue squamous cell carcinoma patients

CA9 is a member of the carbonic anhydrases’ family, that is often expressed in cancer cells under hypoxic condition. However, the role of CA9 in the molecular mechanisms of tongue squamous cell carcinoma (TSCC) pathogenesis remains unclear. CA9 expression was analysed using the TCGA database, and its influence on survival was performed using Kaplan-Meier, LASSO and COX regression analyses. The correlation between CA9 and immune infiltration was investigated by CIBERSORT and ESTIMATE. Moreover, the relationship between CA9 expression and downstream molecular regulation pathways was analysed by GSEA, GO and WGCNA. CA9 expression correlated with clinical prognosis and tumour grade in TSCC. Moreover, CA9 expression potentially contributes to the regulation of cancer cell differentiation and mediates tumour-associated genes and signalling pathways, including apoptosis, hypoxia, G2M checkpoint, PI3K/AKR/mTOR signalling and TGF-beta signalling pathways. However, the follicular helper T cells, regulatory T cells, immune and stromal scores showed no significance between high and low CA9 expression groups. These findings suggested that CA9 plays a critical role of TSCC prognosis and tumour grade. CA9 expression significantly correlated with the regulation of cell differentiation, various oncogenes and cancer-associated pathways.     

Immune Checkpoint Blockade to Improve Tumor Infiltrating Lymphocytes for Adoptive Cell Therapy

Tumor-infiltrating lymphocytes (TIL) has been associated with improved survival in cancer patients. Within the tumor microenvironment, regulatory cells and expression of co-inhibitory immune checkpoint molecules can lead to the inactivation of TIL. Hence, there is a need to develop strategies that disrupt these negative regulators to achieve robust anti-tumor immune responses. We evaluated the blockade of immune checkpoints and their effect on T cell infiltration and function. We examined the ability of TIL to induce tumor-specific immune responses in vitro and in vivo. TIL isolated from tumor bearing mice were tumor-specific and expressed co-inhibitory immune checkpoint molecules. Administration of monoclonal antibodies against immune checkpoints led to a significant delay in tumor growth. However, anti-PD-L1 antibody treated mice had a significant increase in T cell infiltration and IFN-γ production compared to other groups. Adoptive transfer of in vitro expanded TIL from tumors of anti-PD-L1 antibody treated mice led to a significant delay in tumor growth. Blockade of co-inhibitory immune checkpoints could be an effective strategy to improve TIL infiltration and function.     

A dual immune signature of CD8+ T cells and MMP9 improves the survival of patients with hepatocellular carcinoma

The 5-year survival of hepatocellular carcinoma (HCC) is difficult due to the high recurrence rate and metastasis. Tumor infiltrating immune cells (TICs) and immune-related genes (IRGs) bring hope to improve survival and treatment of HCC patients. However, there are problems in predicting immune signatures and identifying novel therapeutic targets. In the study, the CIBERSORT algorithm was used to evaluate 22 immune cell infiltration patterns in gene expression omnibus (GEO) and the cancer genome atlas (TCGA) data. Eight immune cells were found to have significant infiltration differences between the tumor and normal groups. The CD8+ T cells immune signature was constructed by least absolute shrinkage and selection operator (LASSO) algorithm. The high infiltration level of CD8+ T cells could significantly improve survival of patients. The weighted gene co-expression network analysis (WGCNA) algorithm identified MMP9 was closely related to the overall survival of HCC patients. K-M survival and tROC analysis confirmed that MMP9 had an excellent prognostic prediction. Cox regression showed that a dual immune signature of CD8+ T cells and MMP9 was independent survival factor in HCC. Therefore, a dual prognostic immune signature could improve the survival of patient and may provide a new strategy for the immunotherapy of HCC.     

IL2RA+VSIG4+ tumor-associated macrophage is a key subpopulation of the immunosuppressive microenvironment in anaplastic thyroid cancer

Extensive infiltration of tumor-associated macrophages was correlated poor prognosis in anaplastic thyroid cancer (ATC). However, the heterogeneity and characteristics of the ATC-associated macrophages (ATAMs) in ATC remain far from clear. We combined single-cell RNA-sequencing analysis and gene expression microarray datasets to assess the molecular signature of ATAMs. Compared with normal thyroid-associated macrophages (NTAMs), 778 differentially expressed genes (DEGs) significantly changed in ATAMs compared with NTAMs. These DEGs were correlated with oxidative phosphorylation (M2 phenotype) and phagocytosis (M1 phenotype). Moreover, ATAMs highly expressed pro-tumor genes associated with angiogenesis, fibrosis, metalloprotease activity, and metastasis. Notably, we identified one ATC-specific subset, IL2RA⁺ VSIG4⁺ ATAMs, co-expressed M1 and M2 markers. The infiltration of IL2RA⁺ VSIG4⁺ ATAMs showed strong correlation with BRAF and RAS signaling, and its high infiltration was associated with favorable prognosis in thyroid-cancer patients. IL2RA⁺ VSIG4⁺ ATAMs were associated with increased tumor-infiltrating lymphocytes (B cells, CD8⁺ T cells, Tregs). IL2RA⁺ VSIG4⁺ ATAMs interacted with CD8⁺ T cells and Tregs through immune checkpoints (such as LGALS9_HAVCR2), cytokines (such as CXCL10_CXCR3), and receptors (such as CSF1R_CSF1), thereby forming an immunosuppressive microenvironment. Multiplex immunohistochemistry staining and coculture experiment confirmed that ATC cancer cells were able to induce the polarization of IL2RA⁺ VSIG4⁺ ATAMs. Besides, we identified several novel ATC-specific immune checkpoint genes including the immunosuppressive molecule VSIG4, LAIR1, and LILRB2. Expression of VSIG4 was also significantly correlated with tumor-infiltrating lymphocytes (B cells, CD8⁺ T cells, Tregs). In conclusion, our study revealed an ATC-specific ATAM subset with bifunctional phenotype, which provided a comprehensive insight to delineate the molecular characteristics of ATC-associated macrophages.     

PRF1 is a prognostic marker and correlated with immune infiltration in head and neck squamous cell carcinoma

PurposeHead and neck squamous cell carcinoma (HNSCC) is a highly invasive malignancy with poor survival. Perforin (PRF1) plays essential roles in host immunity. Our research intended to identify the correlations of PRF1 with clinical prognosis and tumor immune infiltration in HNSCC.MethodsWe explored PRF1 expression and its associations with the clinical features of HNSCC via the Tumor Immune Estimation Resource (TIMER), Oncomine and The Cancer Genome Atlas (TCGA) databases. The prognostic value of PRF1 for HNSCC was further explored by Kaplan–Meier plotter and TIMER. Finally, the relation between PRF1 and immune infiltration in HNSCC was estimated via CIBERSORT and TIMER.ResultsPRF1 expression was remarkably elevated in HNSCC and associated with clinical stage and HPV infection. High PRF1 expression predicted favorable outcomes in HNSCC, especially in HPV+ HNSCC. Moreover, higher infiltration of CD8+ T cells and CD4+ T cells were found in the PRF1^high group of HNSCC. PRF1 expression in HNSCC was strongly correlated with infiltrating CD8+ T cells and dendritic cells (DCs), with higher relevance in HPV+ HNSCC.ConclusionOur findings suggested that PRF1 could be a novel prognostic biomarker in HNSCC and that its expression was related to immune cell infiltration, which was impacted by HPV status.Key words: PRF1, prognosis, head and neck squamous cell carcinoma, tumor immune infiltration, HPV     

Overexpression of COL11A1 by Cancer-Associated Fibroblasts: Clinical Relevance of a Stromal Marker in Pancreatic Cancer

Background

The collagen11A1 (COL11A1) gene is overexpressed in pancreatic cancer. The expression of COL11A1 protein could be involved in desmoplastic events in pancreatic cancer, but an antibody that specifically stains the COL11A1 protein is not currently available.

Methods and findings

A total of 54 pancreatic ductal adenocarcinomas (PDAC), 23 chronic pancreatitis (CP) samples, and cultured peritumoral stromal cells of PDAC (passages 3-6) were studied. Normal human pancreas tissue samples were obtained through a cadaveric organ donation program.1) Validation of COL11A1 gene overexpression by q-RT-PCR. Findings: the expression of COL11A1 gene is significantly increased in PDAC samples vs. normal and CP samples.2) Analysis of COL11A1 by immunohistochemistry using highly specific anti-proCOL11A1 antibodies. Findings: anti-proCOL11A1 stains stromal cells/cancer-associated fibroblasts (CAFs) of PDAC but it does not stain chronic benign condition (chronic pancreatitis) stromal cells, epithelial cells, or normal fibroblasts.3) Evaluation of the discrimination ability of the antibody. Findings: anti-proCOL11A1 immunostaining accurately discriminates between PDAC and CP (AUC 0.936, 95% CI 0.851, 0.981).4) Phenotypic characterization of proCOL11A1+ stromal cells co-staining with mesenchymal, epithelial and stellate cell markers on pancreatic tissue samples and cultured peritumoral pancreatic cancer stromal cells. Findings: ProCOL11A1+ cells present co-staining with mesenchymal, stellate and epithelial markers (EMT phenotype) in different proportions.

Conclusions/Significance

Detection of proCOL11A1 through immunostaining with this newly-developed antibody allows for a highly accurate distinction between PDAC and CP. Unlike other available antibodies commonly used to detect CAFs, anti-proCOL11A1 is negative in stromal cells of the normal pancreas and almost absent in benign inflammation. These results strongly suggest that proCOL11A1 is a specific marker for CAFs, and thus, anti-proCOL11A1 is a powerful new tool for cancer research and clinical diagnostics.     

Gene expression and DNA methylation analyses suggest that immune process-related ADCY6 is a prognostic factor of luminal-like breast cancer

Breast cancer is a malignant tumor that seriously threatens women's health, and luminal-like cancer subtypes account for the majority of the cases. The purpose of this study was to investigate the relationships among DNA methylation, gene expression profile, and the tumor-immune microenvironment of luminal-like breast cancer, and to identify the potential key genes that regulate immune cell infiltration in luminal-like breast cancer. The ESTIMATE algorithm was applied to calculate immune scores and stromal scores of patients with breast cancer. Kaplan-Meier curves were generated for survival analysis. The clusterProfile package was used for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The protein-protein interaction (PPI) network was constructed using the STRING database and Cytoscape software. Correlations between ADCY6 expression and immune cell infiltration-related pathways were analyzed by gene set variation analysis. R software was used for the statistical analysis and figure generation. Disease-free survival was higher in the immune score-high group than it was in the immune score-low group, while the stromal score had no correlation with prognosis. There were 515 genes that differed in both gene expression and DNA methylation levels, and these genes were mainly enriched in immune process-related pathways. ADCY6 was enriched in module A of the PPI network. Patients with downregulation and hypermethylation of ADCY6 associated with a better prognosis. ADCY6 expression was negatively correlated with the activation of immune process-related signaling pathways, immune checkpoint receptors, and ligands, except for CLEC4G. DNA methylation was found to be involved in the regulation of the key cellular pathways of luminal-like breast cancer immune cell infiltration. Additionally, ADCY6 was identified as a prognostic factor involved in the DNA methylation-regulated immune processes in luminal-like breast cancer.     

Epigenetic silencing of microRNA-149 in cancer-associated fibroblasts mediates prostaglandin E2/interleukin-6 signaling in the tumor microenvironment

Tumor initiation and growth depend on its microenvironment in which cancer-associated fibroblasts (CAFs) in tumor stroma play an important role. Prostaglandin E2 (PGE2) and interleukin (IL)-6 signal pathways are involved in the crosstalk between tumor and stromal cells. However, how PGE2-mediated signaling modulates this crosstalk remains unclear. Here, we show that microRNA (miR)-149 links PGE2 and IL-6 signaling in mediating the crosstalk between tumor cells and CAFs in gastric cancer (GC). miR-149 inhibited fibroblast activation by targeting IL-6 and miR-149 expression was substantially suppressed in the CAFs of GC. miR-149 negatively regulated CAFs and their effect on GC development both in vitro and in vivo. CAFs enhanced epithelial-to-mesenchymal transition (EMT) and the stem-like properties of GC cells in a miR-149-IL-6-dependent manner. In addition to IL-6, PGE2 receptor 2 (PTGER2/EP2) was revealed as another potential target of miR-149 in fibroblasts. Furthermore, H. pylori infection, a leading cause of human GC, was able to induce cyclooxygenase-2 (COX-2)/PGE2 signaling and to enhance PGE2 production, resulting in the hypermethylation of miR-149 in CAFs and increased IL-6 secretion. Our findings indicate that miR-149 mediates the crosstalk between tumor cells and CAFs in GC and highlight the potential of interfering miRNAs in stromal cells to improve cancer therapy.     

A polysaccharide from Dictyophora indusiata inhibits the immunosuppressive function of cancer‐associated fibroblasts

Reversing the function of cancer‐associated fibroblasts (CAFs) may improve the efficacy of cancer therapy. Here, we isolated a novel polysaccharide from Dictyophora indusiata (ZSP4) and examined its effects on the function of prostate CAFs. The supernatant of prostate CAFs can stimulate the proliferation of immune cells and inhibit the growth of CD4+/CD8+ T cells. However, after ZSP4 stimulation, the functions of prostate CAFs were inhibited. The mechanism experiment shows that ZSP4 can stimulate prostate CAFs by down‐regulating the expression of α‐smooth muscle actin. Polysaccharides extracted from Dictyophora indusiata stimulate the proliferation of immune cells and reverse the immune‐suppressive functions of prostate CAFs, shedding new light on the development of novel anticancer strategies. The endocrine therapy used to treat prostate cancer aims to eliminate androgenic activity from prostatic tissue; these therapies are painful and of poor therapeutic effect. In this study, we found that polysaccharides extracted from Dictyophora indusiata may affect the micro‐environment of tumours and inhibit the growth of the tumours. Our results suggest that polysaccharides may modulate negative immune regulation and enhance antitumour immunity, which is important for clinical therapy.     

Prognostic value and immune cell infiltration of hypoxic phenotype‐related gene signatures in glioblastoma microenvironment

Glioblastoma (GBM) is a malignant intracranial tumour with the highest proportion and lethality. It is characterized by invasiveness and heterogeneity. However, the currently available therapies are not curative. As an essential environmental cue that maintains glioma stem cells, hypoxia is considered the cause of tumour resistance to chemotherapy and radiation. Growing evidence shows that immunotherapy focusing on the tumour microenvironment is an effective treatment for GBM; however, the current clinicopathological features cannot predict the response to immunotherapy and provide accurate guidance for immunotherapy. Based on the ESTIMATE algorithm, GBM cases of The Cancer Genome Atlas (TCGA) data set were classified into high‐ and low‐immune/stromal score groups, and a four‐gene tumour environment‐related model was constructed. This model exhibited good efficiency at forecasting short‐ and long‐term prognosis and could also act as an independent prognostic biomarker. Additionally, this model and four of its genes (CLECL5A, SERPING1, CHI3L1 and C1R) were found to be associated with immune cell infiltration, and further study demonstrated that these four genes might drive the hypoxic phenotype of perinecrotic GBM, which affects hypoxia‐induced glioma stemness. Therefore, these might be important candidates for immunotherapy of GBM and deserve further exploration.