Expression, Prognosis, and Immune Infiltrates Analyses of E2Fs in Human Brain and CNS Cancer

Liao, Pan; Han, Shuzhen; Qu, Honglan

doi:https://doi.org/10.1155/2020/6281635

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Abbreviations Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 6281635 | https://doi.org/10.1155/2020/6281635

Expression, Prognosis, and Immune Infiltrates Analyses of E2Fs in Human Brain and CNS Cancer

Pan Liao,¹Shuzhen Han,²and Honglan Qu³

Academic Editor: Bilal Alatas

Received23 Apr 2020

Accepted06 Aug 2020

Published15 Dec 2020

Abstract

Objective. We investigated the expression patterns, potential functions, unique prognostic value, and potential therapeutic targets of E2Fs in brain and CNS cancer and tumor-infiltrating immune cell microenvironments. Methods. We analyzed E2F mRNA expression levels in diverse cancer types via Oncomine and GEPIA databases, respectively. Moreover, we evaluated the prognostic values using GEPIA database and TCGAportal database and the correlation of E2F expression with immune infiltration and the correlation between immune cell infiltration and GBM and LGG prognosis via TIMER database. Then, cBioPortal, GeneMANIA, and DAVID databases were used for mutation analysis, PPI network analysis of coexpressed gene, and functional enrichment analysis. Results. E2F1-8 expression increased in most cancers, including brain and CNS cancer. Higher expression in E2F1, 2, 4, 6, 7, and 8 indicated poor OS of LGG. Higher E2F3–6 and E2F1–8 expressions correlated with poor prognosis and increased immune infiltration levels in CD8+ T cells, macrophages, neutrophils, and DCs in GBM and CD8+ T cells, B cells, CD4+ T cells, neutrophils, macrophages, and DCs in LGG, respectively. Conclusion. E2F1–8 and E2F2–8 could be hopeful prognostic biomarkers of GBM and LGG, respectively. E2F3–6 and E2F1–8 could be likely therapeutic targets in patients with immune cell infiltration of GBM and LGG, respectively.

1. Introduction

E2Fs, as a set of genes encoding transcription factor family (TFs) in higher eukaryotes, play an essential role in cell cycle regulation and DNA synthesis [1]. In some human malignancies, E2F activator expression is deactivated and is deregulated. A report demonstrating that increased expression in E2F3a played a significant role in the development of glioma. This study suggests that E2Fs could promote cancer development [2].

Worldwide, brain and CNS cancer are commonly diagnosed as cancers and the causes of tumor death [3]. Almost 80% of malignant brain cancers are gliomas [4]. In 2016, based on data from groundbreaking research [5], the WHO improved the correct classification of glioma subtypes by integrating tumor morphology and molecular genetic information [6]. The clinical therapy of these cancers includes surgery, radiotherapy, and chemotherapy, which is far from sufficient in combating cancer development. Tumor-infiltrating lymphocytes are an independent predictor of sentinel lymph node status, and survival in cancers [7] and immunotherapy has been considered a promising direction because of its ability to penetrate the blood-brain barrier to treat these tumors since the first discovery of lymphatic in CNS [8]. However, current immunotherapies, for example, anti-CTLA4, demonstrated poor clinical efficacy in brain tumors [9]. CAR T cells also still face substantial obstacles in treating brain tumors [10]. Besides, after the diagnosis of malignant brain and CNS tumors, the overall 5-year relative survival rate was 32.1%, of which the survival rate was only 4.9% after the diagnosis of GBM (glioblastoma multiforme) [4], and metastasis often leads to a poor prognosis. Due to the heterogeneity of tumors, currently, there are some limitations in the biomarkers that can predict prognosis, so it is necessary to search for new biomarkers in this field as prognostic indicators, to improve prognosis. Also, there is an urgent need to identify new therapeutic targets in the brain and CNS cancer so as to effectively improve the accuracy of treatment.

E2F plays a complex and unique role in human brain and CNS tumors. E2F1/DP1 complex was reported to inhibit the mismatch repair proteins MSH2, MSH6, and EXO1, as well as the homologous recombinant protein RAD51, leading to cell apoptosis [11]. Overexpression of E2F1 blocked the proliferation repression caused by miR-342-3p or miR-377 in glioma cells [12], and E2F1 overexpression conferred resistance to cisplatin in glioma cells [13]. Studies have shown that downregulation of E2F2 inhibits glioma cell growth [14]. Overexpression of E2F2 promoted the progression of NSCLC cells [15]. Increased expression in E2F3a played a significant function in glioma development [2]. Change in human GBM cell cycle is related to the decrease of E2F4 level [16]. The silencing of E2F5 inhibited the proliferative ability of GBM cells [17]. Moreover, the combination of E2F6 and MATN1-AS1, which negatively targeted RELA, inhibited the MAPK signaling pathway, thus inhibiting the proliferation and invasion of GBM cells [18]. E2F7 could induce autophagy of glioma [19]. Targeting E2F8 can also regulate the development in glioma [20]. However, the underlying mechanisms of E2F activation or inhibition and the distinct functions of E2Fs and tumor immunology in brain and CNS tumor have not been fully elucidated.

As far as we know, bioinformatics analysis and immune infiltrate analysis have not been used to study the function of E2Fs from tumors of the brain and CNS. As an important part of biological and biomedical research, RNA and DNA research has undergone revolutionary changes due to microarray technology [21]. In addition, immunotherapy, a revolution in cancer treatment, has become a promising strategy [22]. Based on thousands of gene expression or copy number variations published online, we analyzed the expression patterns, potential functions, and prognostic value of E2Fs in brain and CNS cancer. Furthermore, we analyzed the interrelationship between E2Fs and immune cell in different cancer microenvironments to identify underlying therapeutic targets of brain and CNS tumor.

2. Materials and Methods

2.1. GEPIA Dataset Analysis

GEPIA was used to analyze the RNA sequencing expression data from 9,736 tumors and 8,587 normal samples from TCGA and the Genotype-Tissue Expression (GTEx) projects. GEPIA has many functional modules, including differential expression, survival, and correlation analysis [23].

2.2. Oncomine Dataset Analysis

The transcription levels of E2Fs of various tumors and coexpression genes of E2Fs in brain and CNS cancer were determined through analysis in Oncomine cancer microarray datasets.

2.3. The Cancer Genome Atlas Data [24] and cBioPortal Datasets Analyses

cBioPortal currently contains 225 cancer studies to explore multidimensional tumor genomics data sets [25]. We used cBioPortal to analyze E2F alterations in the glioblastoma multiforme (TCGA, Firehose Legacy) dataset, including data from 604 cases with pathology reports, and the brain lower-grade glioma (TCGA, PanCancer Atlas) dataset, including data from 514 cases with pathology reports. The search parameters included mutations, putative copy number alterations (CNAs) from genomic identification of significant targets in cancer (GISTIC), mRNA expression scores (RNA-seq v.2 RSEM), and protein expression scores (reverse phase protein array (RPPA)). Additionally, Kaplan–Meier survival curves were drawn in cBioPortal to evaluate the influence of gene alterations of E2Fs on the overall survival of GBM and LGG patients.

2.4. TCGAportal Database Analysis

TCGAportal was used to further verify the prognostic value of mRNA level of E2F factors in GBM and LGG patients.

2.5. TIMER Database Analysis

TIMER database includes 10,897 samples across 32 cancer types from The Cancer Genome Atlas (TCGA) to estimate the abundance of immune infiltrates [26]. We analyzed the correlation of the E2F expression with the abundance of immune infiltrates, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, via gene modules. In addition, correlations between immune cell infiltration and GBM and LGG prognosis were explored via correlation modules. Lastly, we assessed how the E2F expression correlated with the expression of particular immune infiltrating cell subset markers.

2.6. GeneMANIA Database Analysis

GeneMANIA is a flexible web to construct protein-protein interaction networks, generating hypotheses on gene role, exploring gene list, while prioritizing genes [27]. We visualized the gene network and predicted the functions of E2F coexpression genes via GeneMANIA.

3. Results

3.1. Expression Levels of E2Fs of Brain and CNS Tumor

Oncomine differential expression analysis revealed that the E2F5 expression was remarkably upregulated in brain and CNS tumor in five datasets (Figure 1).

In the French Brain Statistics [28], E2F5 was overexpressed in anaplastic oligodendroglioma versus normal tissues with a fold change of 2.425 (Table 1). In the Sun Brain Statistics (Table 1) [29], E2F5 was also overexpressed in anaplastic astrocytoma with a fold change of 2.545, in oligodendroglioma with a fold change of 2.580, and in glioblastoma with a fold change of 2.877. In the Murat Brain Statistics (Table 1) [30], E2F5 was also overexpressed in glioblastoma with a fold change of 2.016. The Murat Brain Statistics (Table 1) [30] showed that E2F7 was also increased in glioblastoma () compared to normal samples. In addition, in the Sun Brain Statistics (Table 1) [29], E2F7 was overexpressed in glioblastoma (), and in the French Brain Statistics (Table 1) [28], the increased E2F7 expression was found in anaplastic oligodendroglioma compared to normal sample (). In the Sun Brain Statistics (Table 1) [29], E2F8 was also overexpressed in anaplastic astrocytoma (), glioblastoma (), and oligodendroglioma () compared to normal samples.

GEPIA differential expression analysis indicated that the E2F1-8 expressions were higher in GBM tissue compared with normal tissue and that E2F2-8 was higher in LGG (Figure 2).

(a)

(b)

3.2. The Correlations between Different E2Fs in Brain and CNS Cancer

The cBioPortal online tool calculated Pearson correlations. The results showed that the following E2Fs showed significantly positive correlations: E2F1 with E2F2, E2F3, E2F7, and E2F8; E2F2 with E2F1, E2F3, E2F7, and E2F8; E2F3 with E2F1 and E2F2; E2F7 with E2F1, E2F2, and E2F8; and E2F8 with E2F1, E2F2, and E2F7 (Figure 3).

3.3. Gene Alteration of E2Fs in Brain and CNS Cancer Tissue from cBioPortal

cBioPortal analyses revealed that E2Fs were altered in 197 (33·8%) of all the 591 sequenced cases with GBM, which included 3 mutation, 8 cases of amplification, 3 cases of deep deletion, 29 cases of mRNA downregulation, 97 cases of mRNA upregulation, and 5 cases of multiple alterations. Statistical results of the gene alteration frequency indicated that mRNA upregulation occupied the overwhelming majority of alteration types (Figure 4(a)). E2Fs were altered in 324 (43%) of all the 514 sequenced cases with LGG, which included 1 mutation, 6 cases of amplification, 3 cases of deep deletion, 10 cases of mRNA downregulation, 103 cases of mRNA upregulation, and 4 cases of multiple alterations. Among the LGG cases with gene alteration of E2Fs, mRNA upregulation occupied the overwhelming majority of alteration types (Figure 4(b)). Kaplan–Meier overall survival curves showed that gene alterations of E2F1, E2F2, E2F4, and E2F6-8 were remarkably related to overall survival (OS) () among LGG (Figure 4(d)). No significant correlation was identified between gene alteration of E2Fs and OS in GBM patients (Figure 4(c)).

(a)

(b)

(c)

(d)

3.4. The Prognostic Value of E2Fs in GBM and LGG Patients

GEPIA and TCGAportal survival analysis both revealed that decreased E2F1, E2F2, E2F4, and E2F6-8 mRNA levels were remarkably related to overall survival (OS) () among LGG. No significant correlation was identified between E2F mRNA levels and OS in GBM patients (Figure 5).

(a)

(b)

3.5. Prognostic Relevance of Immune Cell Infiltration in GBM and LGG

TIMER analysis showed that the survival time of LGG patients with low B cells, macrophage, CD4+ T cells, neutrophil, and DCs infiltrations was significantly higher than that of LGG patients with high infiltration. However, there were no meaningful findings about GBM (Figure 6).

3.6. Correlations of E2F Expression with Immune Infiltration Levels in GBM and LGG

TIMER analysis revealed that the E2F1 expression is remarkably positively related to tumor purity (, ) in GBM (Figure 7(a)). And the E2F1 expression demonstrated a very weak interrelation with CD8+ T cell infiltration levels in LGG (, ). Also, the E2F2 expression was significantly positively related to cancer purity (, ) and remarkably negatively related to macrophages (, ), neutrophils (, ), and dendritic cells (, ) but no significant correlation with CD4+ T cells in GBM (Figure 7(b)). E2F2 expression level was remarkably positively related to infiltrating levels of B cells (, ), CD4+ T cells (, ), macrophages (, ), neutrophils (, ), and dendritic cells (, ) but no significant correlation with tumor purity (, ) in LGG (Figure 7(b)). The E2F3 expression was remarkably positively related to tumor purity (, ) and neutrophils (, ) in GBM (Figure 7(c)). The E2F3 expression was remarkably positively related to tumor purity (, ) and infiltrating levels of B cells (, ), CD8+ T cells (, ), and DCs (, ) but no significant correlation with macrophages (, ) or neutrophils (, ) in LGG (Figure 7(c)). The E2F4 expression demonstrated a very weak interrelation with dendritic cell (, ) infiltration levels in GBM (, ). The E2F4 expression was remarkably positively related to infiltrating levels of B cells (, ), CD4+ T cells (, ), macrophages (, ), neutrophils (, ), and dendritic cells (, ) in LGG (Figure 7(d)). The E2F5 expression was remarkably positively related to tumor purity (, ) and neutrophils (, ) in GBM (Figure 7(e)). The E2F5 expression was remarkably positively related to tumor purity (, ), B cells (, ), CD8+ T cells (, ), and DCs (, ) but no significant correlation with macrophages (, ) in LGG (Figure 7(e)). The E2F6 expression was remarkably positively related to tumor purity (, ), CD8+ T cells (, ), and macrophages (, ) but no significant correlation with neutrophils (, ) in GBM (Figure 7(f)). The E2F6 expression was remarkably positively related to tumor purity (, ), B cells (, ), macrophages (, ), CD8+ T cells (, ), neutrophils (, ), and DCs (, ) but no significant correlation with CD4+ T cells in LGG (Figure 7(f)). The E2F7 expression was remarkably positively related to tumor purity (, ) in GBM (Figure 7(g)). The E2F7 expression was remarkably positively related to tumor purity (, ), B cells (, ), CD8+ T cells (, ), macrophages (, ), and DCs (, ) in LGG (Figure 7(g)). The E2F8 expression was remarkably positively related to tumor purity (, ) in GBM. The E2F8 expression was remarkably positively related to infiltrating levels of B cells (, ), CD8+ T cells (, ), CD4+ T cells (, ), macrophages (, ), neutrophils (, ), and dendritic cells (, ) in LGG (Figure 7(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 7

Correlation of the E2F expression with immune infiltration level in GBM and LGG: (a) correlation of the E2F1 expression with immune infiltration level in GBM and LGG; (b) correlation of the E2F2 expression with immune infiltration level in GBM and LGG; (c) correlation of the E2F3 expression with immune infiltration level in GBM and LGG; (d) correlation of the E2F4 expression with immune infiltration level in GBM and LGG; (e) correlation of the E2F5 expression with immune infiltration level in GBM and LGG; (f) correlation of the E2F6 expression with immune infiltration level in GBM and LGG; (g) correlation of the E2F7 expression with immune infiltration level in GBM and LGG; (h) correlation of the E2F8 expression with immune infiltration level in GBM and LGG.

3.7. Assessment of the Correlation between E2Fs and Immune Marker Expression

In particular, TIMER analysis revealed that E2F1 was significantly correlated with monocyte markers (CD86), TAM markers (CCL2, IL10), and M2 macrophage markers (VSIG4, MS4A4A) in GBM (; Table 2). E2F2 was significantly correlated with monocyte markers (CD86), TAM markers (CCL2, IL10), and M2 macrophage markers (VSIG4) in GBM (; Table 2). E2F5 was significantly correlated with monocyte markers (CD115) and M1 macrophage markers (IRF5) in GBM (; Table 2). E2F6 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CD68, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (VSIG4, MS4A4A) in GBM (; Table 2). E2F7 was significantly correlated with TAM markers (IL10) in GBM (; Table 2). E2F8 was significantly correlated with M2 macrophage markers (VSIG4) in GBM (; Table 2). Meanwhile, E2F1 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CCL2), M1 macrophage markers (INOS), and M2 macrophage markers (VSIG4) in LGG (; Table 2). E2F2 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CD68, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, VSIG4, MS4A4A) in LGG (; Table 2). E2F4 was significantly correlated with monocyte markers (CD86), TAM markers (CD68), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, VSIG4, MS4A4A) in LGG (; Table 2). E2F8 was significantly correlated with monocyte markers (CD68), TAM markers (CCL2, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, MS4A4A) in LGG (; Table 2).

3.8. Functional Enrichment Analysis of Coexpressed Genes Correlated with E2Fs in Brain and CNS Cancer

Oncomine analyses revealed that E2F1 was positively related to CDK2, RAB3B, FANCI, NUP88, CKS1B, TMEM194A, UBE2S, GINS1, SMC1A, TOPBP1, and USP1. E2F2 was positively corrected with ASAP3, TCEA3, ID3, RPL11, TCEB3, C1orf128, LYPLA2, GALE, HMGCL, SFRS13A, and PNRC2. E2F3 was positively corrected with ARHGAP11A, MYB, BLM, WEE1, CDC7, EXOSC8, HIST1H4C, HIST1H4I, HIST1H4B, IGF2BP, and RRP1B. E2F4 was positively corrected with ELMO3, LRRC29, EXOC3L, KIAA0895L, NOL3, HSF4, FBXL8, TRADD, B3GNT9, FHOD1, and TMEM208. E2F5 was positively corrected with CCDC85B, NDN, SMAD2, TFPI2, DNALI1, CDH11, FGF2, SLC39A6, MBD2, SMAD2, and TGFA. E2F6 was positively corrected with POLA1, STOML2, RPA2, RPA1, JMJD6, CCT7, THOC4, SF3B3, CSTF2, CIAO1, and EIF2B1. E2F7 was positively corrected with CDKN3, CEP55, KIF4A, SHCBP1, CCNB1, CENPK, ECT2, DEPDC1B, NEK2, NCAPG, and KIF14. E2F8 was positively corrected with KIFC1, CENPM, MYBL2, ERCC6L, AURKB, TK1, MCM10, CDCA3, NCAPH, NUSAP1, and MKI67 (Figure 8).

Next, we constructed a network of E2Fs and E2F-correlated coexpressed genes via GeneMANIA. The results showed that the biological functions of the enriched gene set, including DNA replication, chromosome segregation, mitosis, regulation of cell division, cell cycle phase, and regulation of DNA metabolic process, were intimately associated with E2F alterations (Figure 9).

The GO (gene ontology) and KEGG pathways in the DAVID database were analyzed to predict the functions of E2Fs and E2F-correlated coexpressed genes. Based on biological processes, cell components, and molecular functions, GO enrichment analyses predicted the functions of target host genes. We found that cell cycle phase, mitotic cell cycle, DNA replication, cell cycle checkpoint, and mRNA transport were remarkably regulated by E2F alterations in brain and CNS cancer (Figure 10(a)). Chromosome, replication, transcription factor complex, and protein serine/threonine kinase activity were also significantly controlled by these E2F alterations (Figures 10(b) and 10(c)). The corresponding genes are known to be associated with cell cycle.

(a)

(b)

(c)

(d)

KEGG analysis can define the pathways related to the functions of E2F alterations and coexpressed genes correlated with E2Fs. Nine pathways related to the functions of E2F alterations in brain and CNS cancer were found via KEGG analysis (Figure 10(d)). Among these pathways, cfa04110: cell cycle, cfa04350: TGF-beta signaling pathway, ptr05200: pathways in cancer, and cfa04115: p53 signaling pathway were involved in the tumorigenesis and pathogenesis about brain and CNS cancer (Figures 11(a) and 11(b)).

(a)

(b)

4. Discussion

The function of E2F activators in tumorigenesis and their prognostic relevance has been partially demonstrated in several cancers [31–33]; however, further bioinformatics analysis and analysis of immune infiltrates in brain and CNS cancers have not been performed. This is the first study that explored the mRNA expressions of different E2F factors in brain and CNS cancer and investigated their prognostic relevance and correlation with various tumor-infiltrating immune cells. Our findings may help improve treatment of patients with brain and CNS cancers.

E2F1 has been studied extensively in tumors of the brain and CNS [11–13, 34–38]. E2F1/DP1 complex was reported to inhibit the mismatch repair proteins MSH2, MSH6, and EXO1, as well as the homologous recombinant protein RAD51, leading to cell apoptosis [11]. Overexpression of E2F1 blocked the proliferation repression caused by miR-342-3p or miR-377 in glioma cells [12], and E2F1 overexpression conferred resistance to cisplatin in glioma cells [13]. MicroRNA-138 inhibited the development of GBM by targeting E2F1 [34]. PPARα inhibited growth of glioma cells through the E2F1/miR-19a feedback loop [35]. miRNA-320 and miRNA-329 inhibited glioma proliferation by targeting E2F1 [36, 37]. The ECT2/PSMD14/PTTG1 axis promoted glioma proliferation by stabilizing E2F1 [38]. In this study, analysis of Oncomine dataset and GEPIA dataset revealed that elevated expression of E2F1 in human GBM. Analysis of GEPIA and TCGAportal datasets revealed the prognostic value of E2F1 in patients with GBM and LGG. High E2F1 expression and gene alteration of E2F1 in LGG was associated with poor OS; however, the E2F1 expression and gene alteration of E2F1 in GBM showed no significant correlation with OS. In TIMER datasets, the E2F1 expression in GBM showed a positive association with tumor purity while the E2F1 expression in LGG demonstrated a very weak interrelation with CD8+ T cells infiltration level. Specifically, E2F1 was significantly correlated with monocyte markers (CD86), TAM markers (CCL2, IL10), and M2 macrophage markers (VSIG4, MS4A4A) in GBM. E2F1 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CCL2), M1 macrophage markers (INOS), and M2 macrophage markers (VSIG4) in LGG. This suggests that E2F1 plays a role in regulating monocyte polarization in GBM and LGG.

Downregulation of E2F2 was shown to inhibit glioma cell growth in vitro and in vivo by inducing cell cycle arrest in G0/G1 [14]. Overexpression of E2F2 was found to promote invasive growth of NSCLC cells [15]. miR-125b regulated glioblastoma progress via E2F2 [39]. In our study, we observed overexpression of E2F2 in GBM and LGG tissues. However, while high E2F2 expression and gene alteration of E2F2 in LGG was related to low overall survival (consistent with its role as an oncogene), the E2F2 expression and gene alteration of E2F1 showed no significant correlation with OS in the context of GBM. Furthermore, the E2F2 expression was positively related to cancer purity and negatively related to macrophages, neutrophils, and DCs of GBM; however, the E2F2 expression in GBM showed no significant correlation with CD4+ T cells. The E2F2 expression in LGG showed a positive correlation with infiltrating levels of B cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells but not with tumor purity. Specifically, E2F2 was significantly correlated with monocyte markers (CD86), TAM markers (CCL2, IL10), and M2 macrophage markers (VSIG4) in GBM. E2F2 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CD68, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, VSIG4, MS4A4A) in LGG. This suggests that E2F2 plays a role in regulating monocyte polarization in GBM and LGG.

E2F3 overexpression is a cancer suppressor event in many types of tumors, including brain and CNS cancers [32, 40, 41]. miR128-1 inhibited GBM growth by targeting E2F3 [40]. The primary factor of miR-195-mediated cell arrest was E2F3 [41]. Increased expression in E2F3a played a significant role in the development of glioma [2]. Interestingly, SNHG5 promoted the development of glioma by targeting E2F3 [42]. In our study, the E2F3 expression was upregulated in GBM and LGG. The E2F3 expression and gene alteration of E2F3 in LGG and GBM showed no correlation with survival outcomes. However, the E2F3 expression was positively related to tumor purity and neutrophil infiltration remarkably in GBM. About LGG, the E2F3 expression showed a positive correlation with tumor purity and infiltrating levels of B cells, CD8+ T cells, and dendritic cells but not with macrophages and neutrophils.

TFs play a vital role in inhibiting proliferation-related genes. As a member of the TFs E2F family, E2F4 is rich in nonproliferating and differentiated cells [43]. Decreased levels of E2F4 in human GBM cells induced G2/M cell cycle arrest with a concomitant 2-fold decrease with S phase cell [16]. In our study, the expression of E2F4 in GBM and LGG was higher than that in normal tissue. High E2F4 expression in LGG was related to poor overall survival, while the E2F4 expression and gene alteration of E2F4 of GBM showed no correlation with OS. In addition, the E2F4 expression in GBM demonstrated a very weak interrelation with dendritic cell infiltration. E2F4 expression level in LGG was positively related to CD4+ T cells, neutrophil, B cells, macrophage, and DCs. Specifically, E2F4 was significantly correlated with monocyte markers (CD86), TAM markers (CD68), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, VSIG4, MS4A4A) in LGG. This suggests that E2F4 plays a role in regulating monocyte polarization in LGG.

E2F5 was highly expressed in all kinds of tumors, such as prostate tumor [44] and glioblastoma [45]. In addition, silencing of E2F5 inhibited the proliferation of GBM cells and induced cell cycle arrest. More importantly, reintroduction of E2F5 into GBM cells reversed the tumor-suppressive function of miR-1179 [17]. miRNA-129-3p suppressed growth of glioblastoma via targeting E2F5 [45]. Moreover, Let-7c inhibited glioma development by targeting E2F5 [46]. In our study, the expression of E2F5 in GBM and LGG was higher than that in normal tissue. The E2F5 expression and gene alteration of E2F5 in LGG and GBM showed no correlation with OS. However, the E2F5 expression of GBM was positively related to cancer purity and neutrophil infiltration. The E2F5 expression of LGG showed a positive correlation with tumor purity and infiltrating levels of B cells, CD8+ T cells, and dendritic cells but not with infiltrating level of macrophages. Specifically, E2F5 was significantly correlated with monocyte markers (CD115) and M1 macrophage markers (IRF5) in GBM. This suggests that E2F5 plays a role in regulating monocyte polarization in GBM.

The combination of E2F6 and MATN1-AS1 (which negatively targets RELA) was shown to inhibit the MAPK signaling pathway, thus inhibiting the development of GBM [18]. In our study, the expression of E2F6 in GBM and LGG was higher than that in normal tissue. Increased E2F6 expression and gene alteration of E2F6 in LGG was related to poor OS, while the E2F6 expression of GBM showed no correlation with OS. In addition, the E2F6 expression of GBM was positively related to cancer purity, CD8+ T cells, and macrophages but not with neutrophils. The E2F6 expression in LGG showed a significant positive correlation with tumor purity and infiltrating levels of B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells but not with CD4+ T cells. Specifically, E2F6 was significantly correlated with monocyte markers (CD86, CD115), TAM markers (CD68, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (VSIG4, MS4A4A) in GBM. This suggests that E2F6 plays a role in regulating monocyte polarization in GBM.

E2F7 and E2F8 function as transcriptional repressors [47]. E2F7 could induce autophagy of glioma [19]. HOXD-AS1 regulated glioma progress via E2F8 [20]. In our study, the expressions of E2F7 and E2F8 in GBM and LGG were higher than those in normal tissues. High E2F7 and E2F8 expressions and gene alteration of E2F7 and E2F8 in LGG were associated with poor OS. Moreover, the E2F7 and E2F8 expressions were positively related to tumor purity in GBM remarkably. E2F7 expression level in LGG showed a positive correlation with tumor purity and infiltrating levels of B cells, CD8+ T cells, macrophages, and dendritic cells. Similarly, E2F8 expression level in LGG showed a positive correlation with infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. Specifically, E2F7 was significantly correlated with TAM markers (IL10) in GBM. E2F8 was significantly correlated with M2 macrophage markers (VSIG4) in GBM. E2F8 was significantly correlated with monocyte markers (CD68), TAM markers (CCL2, IL10), M1 macrophage markers (IRF5), and M2 macrophage markers (CD163, MS4A4A) in LGG. This suggests that E2F8 plays a role in regulating monocyte polarization.

5. Conclusions

In this study, we systematically analyzed the expressions of E2Fs in brain and CNS cancer tissues, assessed their prognostic value, and investigated their correlation with tumor-infiltrating immune cells. Our findings characterize the heterogeneity and complexity of the molecular biology and the immune infiltrates in brain and CNS cancers. Bioinformatics analysis indicated that the increased expressions of E2F1, 2, and 5–8 in GBM tissues and those of E2F5 and 6 in LGG tissues may play a vital role in oncogenesis. High E2F1–8 expressions may serve as molecular markers to identify high-risk subgroups of patients with GBM. Similarly, high E2F2–8 expressions may help identify high-risk patients in LGG. Our findings suggest that E2F1–8 and E2F2–8 are potential prognostic biomarkers for GBM and LGG, respectively. However, increased expressions of E2F3–6 in GBM correlated with poor prognosis and increased infiltration of CD8+ T cells, macrophages, neutrophils, and DCs. Increased expressions of E2F1-8 in LGG correlated with poor prognosis and increased infiltration of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and DCs. Therefore, E2F3–6 and E2F1–8 are potential therapeutic targets in patients with immune cell infiltration of GBM and LGG, respectively.

Abbreviations

E2Fs:	E2F transcription factors
CNS:	Central nervous system
GBM:	Glioblastoma multiforme
LGG:	Low grade glioma
PPI:	Protein-protein interaction
OS:	Overall survival
TAM:	Tumor-associated macrophages
GO:	Gene ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

No conflict of interest.

Acknowledgments

We would like to thank Hui Sun from the Second Clinical Medical School of Inner Mongolia, University for the Nationalities for helpful suggestions.

References

S. Gaubatz, G. J. Lindeman, S. Ishida et al., “E2F4 and E2F5 play an essential role in pocket protein-mediated G1 control,” Molecular Cell, vol. 6, no. 3, pp. 729–735, 2000.
View at: Publisher Site | Google Scholar
Z.-G. Shen, X.-Z. Liu, C.-X. Chen, and J. M. Lu, “Knockdown of E2F3 inhibits proliferation, migration, and invasion and increases apoptosis in glioma cells,” Oncology Research, vol. 25, no. 9, pp. 1555–1566, 2017.
View at: Publisher Site | Google Scholar
F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
J. S. Barnholtz-Sloan, Q. T. Ostrom, and D. Cote, “Epidemiology of brain tumors,” Neurologic Clinics, vol. 36, no. 3, pp. 395–419, 2018.
View at: Publisher Site | Google Scholar
The Cancer Genome Atlas Research Network, “Comprehensive genomic characterization defines human glioblastoma genes and core pathways,” Nature, vol. 455, no. 7216, pp. 1061–1068, 2008.
View at: Publisher Site | Google Scholar
D. N. Louis, A. Perry, G. Reifenberger et al., “The 2016 World Health Organization classification of tumors of the central nervous system: a summary,” Acta Neuropathologica, vol. 131, no. 6, pp. 803–820, 2016.
View at: Publisher Site | Google Scholar
F. Azimi, R. A. Scolyer, P. Rumcheva et al., “Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma,” Journal of Clinical Oncology : Official Journal of the American Society of Clinical Oncology, vol. 30, no. 21, pp. 2678–2683, 2012.
View at: Publisher Site | Google Scholar
M. Weller, T. Cloughesy, J. R. Perry, and W. Wick, “Standards of care for treatment of recurrent glioblastoma--are we there yet?” Neuro-Oncology, vol. 15, no. 1, pp. 4–27, 2013.
View at: Publisher Site | Google Scholar
S. C. Kurz, L. P. Cabrera, D. Hastie et al., “PD-1 inhibition has only limited clinical benefit in patients with recurrent high-grade glioma,” Neurology, vol. 91, no. 14, pp. e1355–e1359, 2018.
View at: Publisher Site | Google Scholar
K. B. Long, R. M. Young, A. C. Boesteanu et al., “CAR T cell therapy of non-hematopoietic malignancies: detours on the road to clinical success,” Frontiers in Immunology, vol. 9, p. 2740, 2018.
View at: Publisher Site | Google Scholar
D. Aasland, L. Götzinger, L. Hauck et al., “Temozolomide induces senescence and repression of DNA repair pathways in glioblastoma cells via activation of ATR-CHK1, p21, and NF-κB,” Cancer Research, vol. 79, no. 1, pp. 99–113, 2019.
View at: Publisher Site | Google Scholar
Y. Huang and C. Chi, “Glioma cell proliferation is inhibited by miR-342-3p, miR-377 / E2F1 signaling pathway,” Neoplasma, vol. 66, no. 4, pp. 524–531, 2019.
View at: Publisher Site | Google Scholar
F. F. Li, C. Xing, L. L. Wu, and F. Xue, “MiR-205 enhances cisplatin sensitivity of glioma cells by targeting E2F1,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 2, pp. 299–306, 2018.
View at: Publisher Site | Google Scholar
Y. Gao, D. Han, L. Sun et al., “PPARαRegulates the Proliferation of Human Glioma Cells through miR-214 and E2F2,” BioMed Research International, vol. 2018, Article ID 3842753, 10 pages, 2018.
View at: Publisher Site | Google Scholar
X. Zhou and H. Tao, “Overexpression of microRNA-936 suppresses non-small cell lung cancer cell proliferation and invasion via targeting E2F2,” Experimental and Therapeutic Medicine, vol. 16, no. 3, pp. 2696–2702, 2018.
View at: Publisher Site | Google Scholar
D. Yin, H. Zhou, T. Kumagai et al., “Proteasome inhibitor PS-341 causes cell growth arrest and apoptosis in human glioblastoma multiforme (GBM),” Oncogene, vol. 24, no. 3, pp. 344–354, 2005.
View at: Publisher Site | Google Scholar
X. Xu, N. Cai, T. Zhi et al., “MicroRNA-1179 inhibits glioblastoma cell proliferation and cell cycle progression via directly targeting E2F transcription factor 5,” American Journal of Cancer Research, vol. 7, no. 8, pp. 1680–1692, 2017.
View at: Google Scholar
N. Han, L. Yang, X. Zhang et al., “LncRNA MATN1-AS1 prevents glioblastoma cell from proliferation and invasion via RELA regulation and MAPK signaling pathway,” Annals of translational medicine, vol. 7, no. 23, pp. 784–784, 2019.
View at: Publisher Site | Google Scholar
X. Chen, Y. Zhang, Y. Shi et al., “MiR-129 triggers autophagic flux by regulating a novel Notch-1/ E2F7/Beclin-1 axis to impair the viability of human malignant glioma cells,” Oncotarget, vol. 7, no. 8, pp. 9222–9235, 2016.
View at: Publisher Site | Google Scholar
Y. Chen, F. Zhao, D. Cui et al., “HOXD-AS1/miR-130a sponge regulates glioma development by targeting E2F8,” International Journal of Cancer, vol. 142, no. 11, pp. 2313–2322, 2018.
View at: Publisher Site | Google Scholar
S. C. Sealfon and T. T. Chu, “RNA and DNA Microarrays,” in Methods in Molecular Biology, pp. 3–34, Humana Press, Totowa, NJ, 2011.
View at: Google Scholar
Z. K. Ballas, “The 2018 Nobel Prize in Physiology or Medicine: an exemplar of bench to bedside in immunology,” The Journal of Allergy and Clinical Immunology, vol. 142, no. 6, pp. 1752-1753, 2018.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017.
View at: Publisher Site | Google Scholar
The Cancer Genome Atlas Network, “Comprehensive molecular portraits of human breast tumours,” Nature, vol. 490, no. 7418, pp. 61–70, 2012.
View at: Publisher Site | Google Scholar
J. Gao, B. A. Aksoy, U. Dogrusoz et al., “Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal,” Science Signaling, vol. 6, no. 269, p. pl1, 2013.
View at: Publisher Site | Google Scholar
T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
View at: Publisher Site | Google Scholar
D. Warde-Farley, S. L. Donaldson, O. Comes et al., “The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function,” Nucleic Acids Research, vol. 38, suppl_2, pp. W214–W220, 2010.
View at: Publisher Site | Google Scholar
P. J. French, S. M. A. Swagemakers, J. H. A. Nagel et al., “Gene expression profiles associated with treatment response in oligodendrogliomas,” Cancer Research, vol. 65, no. 24, pp. 11335–11344, 2005.
View at: Publisher Site | Google Scholar
L. Sun, A.-M. Hui, Q. Su et al., “Neuronal and glioma-derived stem cell factor induces angiogenesis within the brain,” Cancer Cell, vol. 9, no. 4, pp. 287–300, 2006.
View at: Publisher Site | Google Scholar
A. Murat, E. Migliavacca, T. Gorlia et al., “Stem cell-related “self-renewal” signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma,” Journal of Clinical Oncology : Official Journal of the American Society of Clinical Oncology, vol. 26, no. 18, pp. 3015–3024, 2008.
View at: Publisher Site | Google Scholar
C.-l. Huang, D. Liu, J. Nakano et al., “E2F1 overexpression correlates with thymidylate synthase and survivin gene expressions and tumor proliferation in non small-cell lung cancer,” Clinical cancer research : an official journal of the American Association for Cancer Research, vol. 13, no. 23, pp. 6938–6946, 2007.
View at: Publisher Site | Google Scholar
L. Chen, J. H. Yu, Z. H. Lu, and W. Zhang, “E2F2 induction in related to cell proliferation and poor prognosis in non-small cell lung carcinoma,” International Journal of Clinical and Experimental Pathology, vol. 8, no. 9, pp. 10545–10554, 2015.
View at: Google Scholar
S.-A. Park, J. Platt, J. W. Lee, F. López-Giráldez, R. S. Herbst, and J. S. Koo, “E2F8 as a Novel Therapeutic Target for Lung Cancer,” JNCI: Journal of the National Cancer Institute, vol. 107, no. 9, 2015.
View at: Publisher Site | Google Scholar
S. Qiu, D. Huang, D. Yin et al., “Suppression of tumorigenicity by microRNA-138 through inhibition of EZH2-CDK4/6-pRb-E2F1 signal loop in glioblastoma multiforme,” Biochimica et Biophysica Acta, vol. 1832, no. 10, pp. 1697–1707, 2013.
View at: Publisher Site | Google Scholar
Y. Shi, T. Tao, N. Liu et al., “PPARα, a predictor of patient survival in glioma, inhibits cell growth through the E2F1/miR-19a feedback loop,” Oncotarget, vol. 7, no. 51, pp. 84623–84633, 2016.
View at: Publisher Site | Google Scholar
J. I.-Y. O. N. G. SUN, W. E. I.-Z. H. O. N. G. XIAO, F. E. I. WANG et al., “MicroRNA-320 inhibits cell proliferation in glioma by targeting E2F1,” Molecular Medicine Reports, vol. 12, no. 2, pp. 2355–2359, 2015.
View at: Publisher Site | Google Scholar
B. Xiao, L. Tan, B. He, Z. Liu, and R. Xu, “MiRNA-329 targeting E2F1 inhibits cell proliferation in glioma cells,” Journal of Translational Medicine, vol. 11, no. 1, p. 172, 2013.
View at: Publisher Site | Google Scholar
T. Zhi, K. Jiang, X. Xu et al., “ECT2/PSMD14/PTTG1 axis promotes the proliferation of glioma through stabilizing E2F1,” Neuro-Oncology, vol. 21, no. 4, pp. 462–473, 2019.
View at: Publisher Site | Google Scholar
N. Wu, L. Xiao, X. Zhao et al., “miR-125b regulates the proliferation of glioblastoma stem cells by targeting E2F2,” FEBS Letters, vol. 586, no. 21, pp. 3831–3839, 2012.
View at: Publisher Site | Google Scholar
Z.-N. Shan, R. Tian, M. Zhang et al., “miR128-1 inhibits the growth of glioblastoma multiforme and glioma stem-like cells via targeting BMI1 and E2F3,” Oncotarget, vol. 7, no. 48, pp. 78813–78826, 2016.
View at: Publisher Site | Google Scholar
Q.-Q. Zhang, H. Xu, M.-B. Huang et al., “MicroRNA-195 plays a tumor-suppressor role in human glioblastoma cells by targeting signaling pathways involved in cellular proliferation and invasion,” Neuro-Oncology, vol. 14, no. 3, pp. 278–287, 2012.
View at: Publisher Site | Google Scholar
X. Li, L. Liu, Y. Luo et al., “Long non-coding RNA SNHG5 promotes glioma progression via miR-205/E2F3 axis,” Bioscience Reports, vol. 39, no. 7, 2019.
View at: Publisher Site | Google Scholar
S. Schwemmle and G. P. Pfeifer, “Genomic structure and mutation screening of theE2F4 gene in human tumors,” International Journal of Cancer, vol. 86, no. 5, pp. 672–677, 2000.
View at: Publisher Site | Google Scholar
S.-L. Li, Y. Sui, J. Sun, T. Q. Jiang, and G. Dong, “Identification of tumor suppressive role of microRNA-132 and its target gene in tumorigenesis of prostate cancer,” International Journal of Molecular Medicine, vol. 41, no. 4, pp. 2429–2433, 2018.
View at: Publisher Site | Google Scholar
D.-Z. Fang, Y.-P. Wang, J. Liu et al., “MicroRNA-129-3p suppresses tumor growth by targeting E2F5 in glioblastoma,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 4, pp. 1044–1050, 2018.
View at: Publisher Site | Google Scholar
M. Huang and X. Gong, “Let-7c inhibits the proliferation, invasion, and migration of glioma cells via targeting E2F5,” Oncology Research, vol. 26, no. 7, pp. 1103–1111, 2018.
View at: Publisher Site | Google Scholar
T. Lammens, J. Li, G. Leone, and L. de Veylder, “Atypical E2Fs: new players in the E2F transcription factor family,” Trends in Cell Biology, vol. 19, no. 3, pp. 111–118, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Pan Liao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

577

Downloads

1283

Citations