The novel roles of virus infection-associated gene CDKN1A in chemoresistance and immune infiltration of glioblastoma

Kuan Hu; Juanni Li; Geting Wu; Lei Zhou; Xiang Wang; Yuanliang Yan; Zhijie Xu

doi:10.18632/aging.202519

Research Paper Volume 13, Issue 5 pp 6662—6680

The novel roles of virus infection-associated gene CDKN1A in chemoresistance and immune infiltration of glioblastoma

Kuan Hu^{1,
,} , Juanni Li^{2,
,} , Geting Wu^2, , Lei Zhou^3, , Xiang Wang^4, , Yuanliang Yan^4,5, , Zhijie Xu^2,5, ,

¹ Department of Hepatobiliary Surgery, Xiangya Hospital, Central South University, Changsha 410008, Hunan, China
² Department of Pathology, Xiangya Hospital, Central South University, Changsha 410008, Hunan, China
³ Department of Anesthesiology, Third Xiangya Hospital of Central South University, Changsha 410008, Hunan, China
⁴ Department of Pharmacy, National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha 410008, Hunan, China
⁵ National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha 410008, Hunan, China

* Equal contribution

Received: July 16, 2020 Accepted: December 29, 2020 Published: February 17, 2021

https://doi.org/10.18632/aging.202519
How to Cite

Copyright: © 2021 Hu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Chemoresistance is a common limitation for successful treatment of glioblastoma multiforme (GBM). Recently, virus infections have been demonstrated to be associated with tumorigenesis and chemoresistance in tumors. However, the role of infection-related genes in GBM haven’t been clearly demonstrated. Here, we explored the roles and mechanisms of human T-lymphotropic virus type-1 (HTLV-1) infections in tumorigenesis and chemoresistance in GBM. Four candidate genes, CDKN1A, MSX1, MYC and CHEK2, were identified to be the codifferentially expressed genes between three temozolomide (TMZ) chemotherapy datasets and one HTLV-1 infection gene set. Next, comprehensive bioinformatics data from several databases indicated that only CDKN1A was significantly upregulated in both GBM tissues and cells and showed the greatest prognostic value in GBM patients. Clinical data identified the correlations between CDKN1A expression and clinicopathological parameters of GBM patients. Moreover, CDKN1A was found to be involved in AKT-mediated TMZ resistance of glioma cells. In addition, KEGG analysis of CDKN1A-associated coexpression genes showed that CDKN1A was potentially involved in complement and coagulation cascades pathways in GBM. Finally, TISIDB database was used to investigate the role of CDKN1A in tumor-immune system interactions in GBM. These findings enhanced our understanding of the roles of CDKN1A in tumorigenesis and therapy response in GBM.

Introduction

Glioma is the most common primary intracranial tumor and represents approximately 50% of primary brain tumors that cause significant mortality and morbidity [1–3]. Despite very intensive treatments, glioblastoma multiforme (GBM), which is the most aggressive and devastating glioma histology, has a median overall survival of less than 2 years [4]. Surgery, chemotherapy and radiotherapy are the conventional treatments for GBM, and chemoresistance is a common limitation of successful GBM treatment [5, 6]. Recently, studies have indicated that virus infection-associated chemoresistance has been reported in several cancers [7, 8] but their role in GBM treatment has not been clearly demonstrated. Hence, exploring the roles and mechanisms of virus infection-associated chemoresistance in GBM is a critical need and may offer promise for a novel therapeutic biomarker in GBM.

Cyclin dependent kinase inhibitor 1A (CDKN1A), a cell cycle inhibitor, is directly controlled by p53-dependent or independent pathways and is involved in terminal differentiation, stem cell renewal, apoptosis and cell migration. It has become increasingly clear that CDKN1A can function as an oncogene or as a classic tumor suppressor [9]. Serving as an oncogene, cytoplasmic CDKN1A was reported to promote oncogenic transformation of HER-2-positive breast cancer cells [10]. In contrast, CDKN1A inhibition by LincRNAFEZF1-AS1 promoted cell proliferation in gastric cancer [11]. Moreover, CDKN1A was also reported to play a critical role in the immune microenvironment in tumors. Price JG et al. demonstrated that CDKN1A can regulate Langerhans cell survival and promote Treg cell generation upon exposure to ionizing radiation in cutaneous tumors [12]. However, the detailed function and mechanism of CDKN1A in the tumorigenesis and immune infiltration of GBM have not yet been investigated.

Several studies have reported that some key factors related to the immune response were obviously altered in GBM and subsequently resulted in tumor immune evasion [13–15]. In addition to the classic treatments (e.g. surgery, chemotherapy and radiotherapy), immune therapy is increasingly considered a particularly promising treatment modality for tumors that stimulates the immune system and activates specific immune cells to attack tumor cells [16, 17]. Several parameters related to the immune system have been identified to predict prognoses for some glioma histologies [18, 19]. However, comprehensive research on the immune microenvironment of GBM is still rare and needs further exploration.

The purpose of the present study was to investigate the detailed roles and mechanisms of human T-lymphotropic virus type-1 (HTLV-1) infection-related genes for GBM chemotherapy. By analyzing the available data from several public databases, the HTLV-1 infection-related gene, CDKN1A, was found to influence the GBM chemotherapy response. Moreover, higher expression of CDKN1A was identified in GBM tissues and cell lines and was correlated with poor prognoses of GBM patients. Furthermore, by colony formation assay, transient transfection assay and western blot, we found that CDKN1A was highly expressed in TMZ-resistant glioma cells and involved in AKT-mediated TMZ resistance of glioma cells.

Results

Identification of differentially expressed genes between the cases and controls

Recently, the development of resistance to TMZ has become a common limitation for successful GBM treatment [20]. To screen the co-DEGs between the untreated group and TMZ-treated group, we analyzed the gene expression profiles from three publicly available datasets regarding TMZ chemotherapy from the GEO platform, namely, GSE43452, GSE65363 and GSE80729. Using the screening criteria of p-value < 0.05, we identified 1659 genes in GSE43452, 3962 genes in GSE65363 and 1858 genes in GSE80729 (Table 1). Next, using the Venn analysis provided by FunRich, 74 genes were found to be significantly codifferentially expressed in all three datasets (Figure 1A, Supplementary Table 1). The co-DEGs from three datasets were presumed to have an impact on the chemotherapy response of GBM.

Table 1. The main characteristics of 3 selected studies on gene expression profiling by microarray.

GEOa datasets

Platform

Samples size

DEGsb

Co-DEGs

Submission
date

References

Control

Treatment

GSE43452

GPL10558

1659

Jan 11, 2013

[46]

GSE65363

GPL570

3962

Jan 28, 2015

N/A c

GSE80729

GPL10558

1858

Apr 27, 2016

[47]

Note: ^aGEO: Gene Expression Omnibus datasets.

^bDEGs: differentially expressed genes.

^cFurther information can be found at https://www.ncbi.nlm.nih.gov/geo/.

Venn diagrams of gene expression microarray datasets. (A) The 74 co-DEGs in three publicly available datasets including GSE43452, GSE65363 and GSE80729. (B) CDKN1A, MSX1, MYC and CHEK2 were consistently identified between 74 co-DEGs and one HTLV-1 infection-related gene set. Each rectangle represents a dataset. The number in each overlapping region represents the number of differentially expressed genes. The intersection in the middlemost area represents the number of genes that were consistently differentially expressed in all these datasets.

Figure 1. Venn diagrams of gene expression microarray datasets. (A) The 74 co-DEGs in three publicly available datasets including GSE43452, GSE65363 and GSE80729. (B) CDKN1A, MSX1, MYC and CHEK2 were consistently identified between 74 co-DEGs and one HTLV-1 infection-related gene set. Each rectangle represents a dataset. The number in each overlapping region represents the number of differentially expressed genes. The intersection in the middlemost area represents the number of genes that were consistently differentially expressed in all these datasets.

Virus infections have been reported to be involved in the development of GBM [21, 22]. Next, we explored the roles of HTLV-1 infections on the development and chemoresistance of GBM. The HTLV-1 infection-related gene set was extracted from the MalaCards database and four HTLV-1 infection-related genes, namely, CDKN1A, MSX1, MYC and CHEK2, were identified in the aforementioned co-DEGs of the three datasets (Figure 1B). These four genes were hypothesized to have an influence on virus infection-associated chemoresistance in GBM.

CDKN1A shows the greatest prognostic value in GBM

The associations between the expression levels of CDKN1A, MSX1, MYC and CHEK2 and prognosis in GBM patients were analyzed using the GEPIA database. Notably, high CDKN1A expression was marginally associated with poor prognosis in GBM patients (Figure 2A), but the other three genes did not exhibit significant prognostic value (Figure 2B, 2C). Similar findings were identified by the GlioVis database. As shown in Supplementary Figure 1, CDKN1A expression levels significantly impacted prognosis in GBM but MSX1, MYC and CHEK did not have obvious associations with survival in GBM. Based on these results, CDKN1A, as the only gene that showed obvious prognostic significance in GBM, was selected for further study.

Prognostic values of CDKN1A, MSX1, MYC and CHEK2 in GBM. (A–D) Kaplan-Meier analysis of overall survival between the samples with high expression of the four genes and those with low expression in GBM by using the GEPIA database.

Figure 2. Prognostic values of CDKN1A, MSX1, MYC and CHEK2 in GBM. (A–D) Kaplan-Meier analysis of overall survival between the samples with high expression of the four genes and those with low expression in GBM by using the GEPIA database.

We then analyzed the correlation between CDKN1A expression levels and clinicopathological parameters of GBM patients. The clinical data of GBM patients were downloaded from the Wanderer database. As shown in Table 2, CDKN1A expression was significantly associated with age (p = 0.030) and vital status (p = 0.011).

Table 2. Correlation of CDKN1A with clinicopathological parameters in GBM.

Characteristics

Case

CDKN1A

Mean ± SD

P-value

Gender

Female

11.35 ± 1.27

Male

11.20 ± 0.94

0.451^b

Age

≥ 50

11.43 ± 0.92

< 50

10.34 ± 1.51

0.030^b

Race

White

11.26 ± 1.03

Black or African American

11.40 ± 2.60

0.829^b

Ethnicity

Not Hispanic or Latino

11.20 ± 1.15

Hispanic or Latino

11.30 ± 0.00

0.955^b

Vital status

Alive

10.78 ± 1.11

Dead

11.71 ± 0.86

0.011^b

Tumor status

With tumor

11.30 ± 1.18

Tumor free

11.28 ± 0.68

0.858^b

Performance status timing

Pre-adjuvant therapy

11.29 ± 1.20

Post-adjuvant therapy

10.52 ± 2.19

Pre-operative

11.39 ± 0.52

0.828^a

Significant P values are underlined.

^aKruskal–Wallis rank test.

^bMann–Whitney U test.

CDKN1A is upregulated in GBM tissues and cell lines and impacts the treatment outcomes of GBM patients

The expression profiles of CDKN1A were analyzed using two independent bioinformatics databases, GEPIA and GlioVis. First, the analysis results from the GEPIA database indicated that CDKN1A mRNA expression was higher in GBM tissues than that in noncancerous tissues (Figure 3A). Next, data from the GlioVis database revealed that CDKN1A mRNA expression levels were significantly upregulated in GBM tissues (Figure 3B), which was consistent with those obtained using the GEPIA database. Both databases demonstrated that CDKN1A mRNA expression levels were upregulated in GBM tissues. Meanwhile, we evaluated the protein expression levels of CDKN1A in GBM patients by analyzing the immunohistochemical data from the Human Protein Atlas. As shown in Figure 3C, significantly elevated levels of CDKN1A was found in the glioma tissues. Further analysis using the CCLE database to study CDKN1A expression profiles in GBM cell lines and the heatmap revealed elevated CDKN1A expression levels in most GBM cell lines (Figure 3D).

Analysis of the CDKN1A expression levels in GBM tissues and cell lines. (A, B) The mRNA expression of CDKN1A in GBM tissues was detected by using the GEPIA and GlioVis databases. (C) the Human Protein Atlas project showed representative immunohistochemical images of CDKN1A in GBM tissues compared with surrounding normal tissues. (D) The expression levels of CDKN1A in GBM cell lines was detected by using the CCLE database.

Figure 3. Analysis of the CDKN1A expression levels in GBM tissues and cell lines. (A, B) The mRNA expression of CDKN1A in GBM tissues was detected by using the GEPIA and GlioVis databases. (C) the Human Protein Atlas project showed representative immunohistochemical images of CDKN1A in GBM tissues compared with surrounding normal tissues. (D) The expression levels of CDKN1A in GBM cell lines was detected by using the CCLE database.

To further determine the effect of CDKN1A on the treatment outcomes of GBM patients, we checked the expression level of CDKN1A in three microarray datasets related to TMZ chemotherapy. The data from these three datasets, GSE43452, GSE65363 and GSE80729, all showed that treatment with the anticancer agent TMZ clearly upregulated CDKN1A expression in human GBM cells (Figure 4A–4C). In addition, the expression levels of CDKN1A might be negative correlation with TMZ activity in 79 glioma cells from CellMinerCDB [23] (Figure 4D, Supplementary Table 4).

The effect of CDKN1A on the therapeutic response of GBM patients. (A–C) GSE43452, GSE65363 and GSE80729 are three microarray datasets related to TMZ chemotherapy and are employed to identify the impacts of CDKN1A expression levels on GBM therapy. (D) the expression levels of CDKN1A might be negative correlation with TMZ activity in 79 glioma cells from CellMinerCDB.

Figure 4. The effect of CDKN1A on the therapeutic response of GBM patients. (A–C) GSE43452, GSE65363 and GSE80729 are three microarray datasets related to TMZ chemotherapy and are employed to identify the impacts of CDKN1A expression levels on GBM therapy. (D) the expression levels of CDKN1A might be negative correlation with TMZ activity in 79 glioma cells from CellMinerCDB.

The above results indicated that CDKN1A expression levels were upregulated in both GBM tissues and cell lines and may influence the chemotherapy responses of GBM patients.

CDKN1A was involved in TMZ resistance of glioma cells

To verify the promoting effect of CDKN1A on cell chemoresistance, we detect CDKN1A expression in TMZ-resistant glioma cell lines (U87-R and T98G-R cells) and their parental cell lines (U87 and T98G cells). First, using the colony formation assay, we proved that U87-R and T98G-R cells were indeed significantly more resistant to the therapy of TMZ compared with U87 and T98G cells, respectively (Figure 5A–5D). After this, CDKN1A and p-CDKN1A expression was further found to be higher in U87-R and T98G-R cells compared with that in U87 and T98G cells, respectively (Figure 5E–5G). These data thus suggested that CDKN1A may played a critical role in TMZ resistance of GBM cells.

CDKN1A and p-CDKN1A were highly expressed in TMZ-resistant glioma cells. (A, B) Compared with U87 cells, U87-R cells were significantly more resistant to the therapy of TMZ by using the colony formation assay. (C, D) Compared with T98G cells, T98G-R cells were significantly more resistant to the therapy of TMZ by using the colony formation assay. (E–G) Compared with U87 and T98G cells, the protein expression levels of CDKN1A and p-CDKN1A were significantly higher in U87-R and T98G-R cells, respectively. The results were presented as means ± SD (n = 3 for each panel). Statistical significance was concluded at *P

Figure 5. CDKN1A and p-CDKN1A were highly expressed in TMZ-resistant glioma cells. (A, B) Compared with U87 cells, U87-R cells were significantly more resistant to the therapy of TMZ by using the colony formation assay. (C, D) Compared with T98G cells, T98G-R cells were significantly more resistant to the therapy of TMZ by using the colony formation assay. (E–G) Compared with U87 and T98G cells, the protein expression levels of CDKN1A and p-CDKN1A were significantly higher in U87-R and T98G-R cells, respectively. The results were presented as means ± SD (n = 3 for each panel). Statistical significance was concluded at *P < 0.05.

Studies have shown that AKT signaling pathway plays an important role in chemoresistance of GBM [24]. Thus, we then investigated whether the effect of CDKN1A on cell chemoresistance was associated with AKT activity. Western blot analysis showed that the expression of p-AKT, a constitutively active form of AKT, was significantly higher in TMZ-resistant glioma cells compared with that in their parental cells, whereas the expression of total AKT has no change, indicating that AKT played a critical role in chemoresistance of GBM (Figure 6A). We then treated TMZ-resistant glioma cells with the specific inhibitor to AKT, MK2206, and found that MK2206 treatment significantly suppressed phosphorylation of CDKN1A at Ser473 (Figure 6B). However, the down-regulated expression of CDKN1A had no effect on AKT expression in TMZ-resistant glioma cells, suggesting that CDKN1A was a downstream regulatory molecule of AKT (Figure 6C). Next, we treated TMZ-resistant glioma cells with both MK2206 and siCDKN1A, as shown in Figure 6D, 6E, the inhibitory effect on p-CDKN1A expression was remarkably increased by combinational treatment. We performed a colony formation assay and observed that TMZ-resistant glioma cells treated with both MK2206 and siCDKN1A had a significantly lower number of colony formation compared with the control group and either of the individual treatment groups (Figure 6F, 6G). The above results indicated that CDKN1A was involved in AKT-mediated TMZ resistance of glioma cells.

CDKN1A was involved in TMZ resistance of glioma cells. (A) Western blot for phospho-Akt (Ser-473), Akt in U87, U87-R T98G, and T98G-R cells. (B) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R or T98G-R cells treated with DMSO and the specific inhibitor to AKT MK2206. (C) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R or T98G-R cells treated with siCtrl and siCDKN1A. (D) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. (E) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in T98G-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. (F, G) Colony formation assay of U87-R or T98G-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. The results were presented as means ± SD (n = 3 for each panel). Statistical significance was concluded at *P

Figure 6. CDKN1A was involved in TMZ resistance of glioma cells. (A) Western blot for phospho-Akt (Ser-473), Akt in U87, U87-R T98G, and T98G-R cells. (B) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R or T98G-R cells treated with DMSO and the specific inhibitor to AKT MK2206. (C) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R or T98G-R cells treated with siCtrl and siCDKN1A. (D) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in U87-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. (E) Western blot for phospho-Akt (Ser-473), Akt, phospho-CDKN1A, and CDKN1A in T98G-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. (F, G) Colony formation assay of U87-R or T98G-R cells treated with DMSO, MK2206, siCDKN1A and MK2206 plus siCDKN1A. The results were presented as means ± SD (n = 3 for each panel). Statistical significance was concluded at *P < 0.05, **P < 0.01.

Functional enrichment analysis of CDKN1A-associated coexpression genes

To further understand the potential role of CDKN1A in GBM development, we performed functional enrichment annotation analysis of its coexpressed genes. We extracted the datasets (TCGA-GBM) from TCGA to screen DEGs that interact with CDKN1A from the GlioVis database. As shown in Figure 7A and Supplementary Table 3, 99 genes were acquired using screening criteria of |LogFC| ≥1.3 and P≤0.05. Then, a PPI network of the genes that coexpressed with CDKN1A was created with the STRING database and Cytoscape software (Figure 7B). To further understand the biological functions of these coexpressed genes, KEGG pathway analyses were conducted by using the DAVID database. The results showed that these coexpressed genes were functionally enriched in several pathways which were mostly those involving complement and coagulation cascades (Figure 7C).

Figure 7. Functional enrichment analysis of CDKN1A-associated coexpression genes. (A) The coexpression genes of CDKN1A were shown as a heatmap via GlioVis. (B) The PPI network of CDKN1A-associated coexpression genes was created by the STRING and Cytoscape software. (C) The significant KEGG pathways associated with the CDKN1A coexpression genes were identified using the DAVID database.

Regulation of immune molecules by CDKN1A

Guglietta S et al. reported that the complement and coagulation cascade pathways were associated with protumorigenic phenotypes of immune cells and protection of tumor cells from immune attack which ultimately favor the development and metastasis of tumors [25]. Therefore, we investigated the associations between CDKN1A expression and lymphocytes and immunomodulators using the TISIDB database. Figure 8A shows the positive correlation between CDKN1A expression and tumor-infiltrating lymphocytes (TILs) in GBM patients. Additionally, the lymphocytes that exhibited the most significant correlations included central memory CD4 T cells (Tcm_CD4, Spearman: ρ= 0.543), plasmacytoid dendritic cells (pDC, Spearman: ρ= 0.508), activated dendritic cells (Act_DC, Spearman: ρ= 0.508) and type-1 T helper cells (Th1, Spearman: ρ= 0.454) (Figure 8B). Immunomodulators have been classified into three types of molecules which include immunoinhibitors, immunostimulators and major histocompatibility complex (MHC) molecules. Figure 8C shows the correlation between CDKN1A expression and immunostimulators; the immunostimulators displaying the most significant correlations included TNFSF13 (Spearman: ρ= 0.429), CD276 (Spearman: ρ= 0.420), TNFRSF14 (Spearman: ρ= 0.420) and CD86 (Spearman: ρ= 0.351) (Figure 8D). Figure 8E shows the correlations between CDKN1A expression and immunoinhibitors; those immunoinhibitors displaying the most significant correlations included TGFB1 (Spearman: ρ= 0.498), PVRL2 (Spearman: ρ= 0.451), CD274 (Spearman: ρ= 0.420) and IL10RB (Spearman: ρ= 0.417) (Figure 8F). Figure 8G shows the correlations between CDKN1A expression and MHC molecules; the MHC molecules that showed the most significant correlations included HLA-A (Spearman: ρ= 0.507), HLA-B (Spearman: ρ= 0.488), TAP1 (Spearman: ρ= 0.438), and TAPBP (Spearman: ρ= 0.437) (Figure 8H). We speculate that since CDKN1A was significantly correlated with various types of tumor-infiltrating lymphocytes, immunoinhibitors, immunostimulators and MHC molecules in GBM, it might exert a more significant effect on immune fingerprinting in GBM.

Correlation of CDKN1A expression with lymphocytes and immunomodulators in GBM. (A) The correlation between CDKN1A expression and TILs. (B) The top four TILs showing the most significant correlations with CDKN1A expression. (C) The correlation between CDKN1A expression and immunostimulators. (D) The top four immunostimulators showing the most significant correlations with CDKN1A expression. (E) The correlation between CDKN1A expression and immunoinhibitors. (F) The top four immunoinhibitors showing the most significant correlation with CDKN1A expression. (G) The correlation between CDKN1A expression and MHC molecules. (H) The top four MHC molecules showing the most significant correlation with CDKN1A expression.

Figure 8. Correlation of CDKN1A expression with lymphocytes and immunomodulators in GBM. (A) The correlation between CDKN1A expression and TILs. (B) The top four TILs showing the most significant correlations with CDKN1A expression. (C) The correlation between CDKN1A expression and immunostimulators. (D) The top four immunostimulators showing the most significant correlations with CDKN1A expression. (E) The correlation between CDKN1A expression and immunoinhibitors. (F) The top four immunoinhibitors showing the most significant correlation with CDKN1A expression. (G) The correlation between CDKN1A expression and MHC molecules. (H) The top four MHC molecules showing the most significant correlation with CDKN1A expression.

Discussion

Our research group aimed to investigate critical and novel biomarkers involved in the development of virus infection-associated chemoresistance of GBM. Intriguingly, four genes including CDKN1A, MSX1, MYC and CHEK2 were identified by screening co-DEGs in three TMZ chemotherapy-related databases and one HTLV-1 infection-related gene set. However, only CDKN1A expression showed a significant correlation with poor prognosis in GBM patients with remarkable consistency between two independent bioinformatics platforms. CDKN1A, which was also found to be upregulated in GBM tissues and cell lines and impacted treatment outcomes of GBM patients, was selected for further in-depth investigations of its biological processes and signaling pathways as well as its correlations with immune regulation.

CDKN1A is a cyclin-dependent kinase inhibitor [26, 27], and an important versatile cell cycle regulator that is involved in cell migration and autophagy and is often deregulated in various cancers [9]. CDKN1A appears to exhibit dual-role behavior and acts either as an oncogene or tumor suppressor depending on the cell types or cellular localization [28]. A complex of cytoplasmic CDKN1A together with ectopic expressed cytosol-resided p53/CDKN1A complex suppressed cell invasion by targeting Bcl-2 family proteins in non-small cell lung cancer [29]. It has also been shown that the wild-type p53/CDKN1A complex can induce slug protein degradation and further suppress cell invasion [30]. In an opposite manner, Akt-activated CDKN1A was found to accelerate tumor onset and promote lung metastasis in vivo [31]. Cytoplasmic CDKN1A was reported to promote cell migration and invasion abilities in gastric cancer [32]. Similarly, Kreis NN et al. revealed that suppression of CDKN1A could inhibit migration and invasion in various cancer cell lines [33]. In some cases, CDKN1A phosphorylation by activated AKT1 prevents the nuclear translocation of CDKN1A and retains it in the cytoplasm, which is crucial for the pro-survival functions of CDKN1A [34, 35]. Cytoplasmic CDKN1A exhibits its oncogenic function might be dependent on non-traditional cytoplasmic targets of the protein. However, the detailed roles of CDKN1A in human brain tumors, especially GBM, have rarely been studied. In this study, we first demonstrated the prognostic value and potential roles of CDKN1A in GBM biology. We found that CDKN1A was involved in AKT-mediated TMZ resistance of glioma cells.

Recently, increasing evidence has indicated that tumor microenvironments and immune infiltration play a vital role on tumor development and chemoresistance [36–39]. Immunotherapy has emerged as a powerful tool for tumor treatment that targets and attacks tumor cells by stimulating the body's own immune system to recognize tumor cells and activate specific immune cells [40]. Glioblastoma is a very lethal form of human brain cancer with an extremely poor prognosis, which makes the development of novel therapeutic strategies targeting GBM of paramount importance. Immune therapy with the combination of conventional treatment (such as surgery, chemotherapy, and radiation) represents a particularly promising approach [41]. In this study, the correlation between CDKN1A and the immune system was assessed with the TISIDB database, and the results revealed that CDKN1A had the highest correlation with tumor-infiltrating lymphocytes including Tcm_CD4, pDC, Act_DC and Th1. Furthermore, CDKN1A had the most significant correlations with immunostimulators (such as TNFSF13, CD276, TNFRSF14 and CD86), immunoinhibitors (such as TGFB1, PVRL2, CD274 and IL10RB) and MHC molecules (such as HLA-A, HLA-B, TAP1 and TAPBP). DC-based vaccination is a basic form of immunotherapy and is critical for initiating and boosting anti-GBM immunity. Audencel, a DC-based cancer vaccine, was found to significantly upregulate Th1-related immunovariables and had effects on the immune system of GBM patients [42]. CD276 (B7-H3) is an immune checkpoint molecule and plays critical roles in T-cell suppression in GBM, which leads to better understanding of target pathways for immunotherapy in GBM [43]. As the most well-known recent biomarker in the cancer research field, PD-L1, also named CD274, was found to be upregulated in IDH1/2 wild-type GBM [44]. Blockade of the PD-L1/PD-1 axis has been reported to reduce Treg expansion and further improve T cell function and thus prevent immunosuppression in GBM [45]. Together, these findings suggest that CDKN1A, which is associated with these immune molecules, plays a vital role in immune escape in GBM microenvironments and can potentially serve as an immunotherapeutic target for GBM.

Nevertheless, there were several limitations in the present study that need to be mentioned. First, patient ethnicities in the TCGA database were mainly white and black, and more studies including other ethnic groups should be investigated further. Additionally, most of the datasets analyzed here were collected retrospectively, so further prospective studies are needed. Moreover, although attractive findings were identified in this study, more functional and mechanistic experiments and large-scale clinical trials are still required to confirm the clinical application value of CDKN1A.

In summary, this is the first report that CDKN1A is associated with the HTLV-1 infection-related chemoresistance of GBM and shows a significant correlation with poor prognosis. Moreover, CDKN1A is upregulated in GBM and might serve as a promising biomarker in the treatment of GBM patients. Furthermore, CDKN1A was found to be involved in AKT-mediated TMZ resistance of glioma cells. In addition, CDKN1A expression is related to tumor-infiltrating lymphocytes and immunomodulators. Therefore, our findings suggest that CDKN1A likely plays a critical role in immune cell infiltration and is a promising prognostic biomarker in patients with GBM.

Materials and Methods

Data acquisition and reanalysis using different bioinformatics tools

The roles of HTLV-1 infection-associated chemoresistance of GBM were explored using several bioinformatics resources which are summarized in Table 3. Three temozolomide (TMZ) therapeutic transcriptome microarray datasets, GSE43452 [46], GSE65363 and GSE80729 [47], were downloaded from the gene expression omnibus (GEO) database (Table 1) [48]. Codifferentially expressed genes (co-DEGs) were identified in these three datasets using a Venn diagram. Next, we downloaded a comprehensive HTLV-1 infection related gene set from the MalaCards database (Supplementary Table 2) [49]. CDKN1A and other three HTLV-1 infection related genes were identified in the co-DEGs of the three GEO datasets and a corresponding Venn diagram was created by FunRich software [50].

Table 3. The main bioinformatics tools used to analyze the role of CDKN1A in GBM.

Database	Samples	URL	References
MalaCards	Tissues	https://www.malacards.org/	[49]
FunRich	–	http://www.funrich.org	[50]
GEPIA	Tissues	http://gepia.cancer-pku.cn/	[51]
GlioVis	Tissues	http://gliovis.bioinfo.cnio.es/	[52]
Wanderer	Tissues	http://maplab.imppc.org/wanderer/	[53]
THPA	Tissues	http://www.proteinatlas.org/	[54]
CCLE	Cell lines	https://portals.broadinstitute.org/ccle/	[55]
STRING	–	http://string-db.org/	[56]
Cytoscape	–	https://cytoscape.org/	[58]
DIVID	Tissues	https://david.ncifcrf.gov/	[60]
TISIDB	Tissues	http://cis.hku.hk/TISIDB	[61]

Gene expression profiling interactive analysis (GEPIA) and GlioVis were used to explore the relationships between the values of the four HTLV-1 infection-related genes and GBM prognoses. GEPIA, an interactive web server, can analyze RNA sequencing expression from TCGA and the GTEx projects [51]. GlioVis is a web-based tool and contains a large collection of brain tumor entries [52]. We then estimated the relevance among the CDKN1A expression and clinicopathological parameters of GBM patients using clinical data downloaded from Wanderer, which is an interactive viewer containing gene expression profiles from TCGA [53].

The expression levels of CDKN1A in GBM tissues and cell lines were analyzed in the following four databases. The GEPIA and GlioVis databases were also used to analyze the mRNA expression levels of CDKN1A in GBM tissues. The Human Protein Atlas (HPA) was used to analyze the protein expression levels of CDKN1A in GBM tissues [54]. An encyclopedia of cancer cell lines (CCLE) [55] was used to analyze CDKN1A expression in different GBM cell lines. Chemotherapy-related datasets, including GSE43452, GSE65363 and GSE80729, were used to explore the impact of CDKN1A expression on the chemotherapy response of GBM.

Next, CDKN1A-associated coexpression genes in GBM pathology were downloaded from GlioVis. We constructed a protein-protein interaction (PPI) of these coexpression genes from STRING which is a database that integrates known and predicted PPI networks from many organisms [56, 57]. Then, Cytoscape software, a tool that visually integrates networks with phenotypes and gene expression profiles was used to perform detailed visualization [58, 59]. Moreover, we utilized the database for annotation, visualization and integrated discovery (DAVID) to conduct a KEGG pathway analysis of the coexpression genes with CDKN1A in GBM [60].

The TISIDB database is an integrated repository portal for tumor-immune system interactions [61]. In this study, using the TISIDB database, we analyzed the correlation between CDKN1A expression in GBM patients and tumor-infiltrating lymphocytes (TILs) and immunomodulators.

Cells and reagents

The TMZ-resistant glioma cell lines (U87-R and T98G-R) and their parental cell lines (U87 and T98G) were established and cultured as previously described [20]. MK2206 was purchased from Selleck Chemicals and dissolved in dimethylsulfoxide (DMSO) (Sigma, USA). The exposed concentrations of TMZ, MK2206 were 200 and 5mM, respectively.

Colony formation assay

For the colony formation assay, the methods were described previously [62].

Protein extraction and quantification

Protein extraction was carried out as previously described [63]. Protein concentration was determined by DC (detergent compatible) protein assay (Bio-Rad Laboratories, USA) according to manufacturer’s instructions.

Western blot

Protein samples were resolved by SDS-PAGE, transferred to polyvinylidene difluoride membrane, and hybridized with antibodies specific to CDKN1A (2946; Cell Signaling Technology), p-CDKN1A (PA5-37519; Invitrogen, USA), AKT (9272s, Cell Signaling Technology, USA), p-AKT (4060, Cell Signaling Technology, USA), GAPDH (sc-47724; Santa Cruz, USA). The blots were developed by the enhanced chemiluminescence reagent (Thermo Scientific Pierce ECL, USA), visualization of the protein bands was conducted in the ChemiDoc XRS system (Bio-Rad, USA).

Transient transfection

For transient transfections, the two TMZ-resistant glioma cell lines (U87-R and T98G-R) were transfected with siRNAs for CDKN1A (siCDKN1A, 5, - GAUGUCCGUCAGAACCCAUGCGGCA-3’) using Lipofectamine 3000 reagent (Invitrogen, USA) according to the manufacturer's instruction. After the indicated incubation times, the cells were harvested and analyzed.

Statistical analyses

Statistical analyses were performed with SPSS 12.0 software (IBM Analytics, USA). All experiments were performed in at least triplicate with mean ± SD subjected to Student’s t-test. Kaplan-Meier analysis was performed to analyze survival rates for GBM. The differential mRNA expression between cancer and noncancer tissues or between the control group and treatment group were analyzed using Student’s t-test. The associations between CDKN1A expression and clinicopathologic characteristics in GBM patients were assessed using the Kruskal–Wallis rank test or Mann–Whitney U test. Correlations between genes were analyzed using Pearson's correlation coefficient. *p < 0.05, and **p < 0.01 were defined as statistically significant.

Ethical statement

The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. None of the data have been previously published or appeared in copyrighted form elsewhere, and not previously published or unpublished data were cited in this paper. No ethics approval was required for this bioinformatics article, as it did not involve patients or patient data.

Supplementary Materials

Supplementary Figure 1

Supplementary Tables 1, 2 and 3

Supplementary Table 4

Author Contributions

Conception and design: YL Yan and ZJ Xu. Administrative support: GT Wu, L Zhou and X Wang. Provision of study materials or patients: GT Wu, L Zhou and X Wang. Collection and assembly of data: JN Li, K Hu, YL Yan and ZJ Xu. Data analysis and interpretation: JN Li and K Hu. Manuscript writing: All authors. Final approval of manuscript: All authors.

Conflicts of Interest

All authors have completed the ICMJE uniform disclosure form. The authors have no conflicts of interest to declare.

Funding

This study is supported by grants from the Natural Science Foundation of Hunan Province (2018JJ3820, 2019JJ50932, 2020JJ5934), National Natural Science Foundation of China (81703036, 81803035, 81702824), China Postdoctoral Science Foundation (2020M672521, 2017M610510), and Dr. Yuanliang Yan is right now a Postdoctoral Fellow at the Center for Molecular Medicine, Xiangya Hospital, Central South University (248485).

References

1. Zhao M, van Straten D, Broekman ML, Préat V, Schiffelers RM. Nanocarrier-based drug combination therapy for glioblastoma. Theranostics. 2020; 10:1355–72. https://doi.org/10.7150/thno.38147 [PubMed]
2. Bi J, Chowdhry S, Wu S, Zhang W, Masui K, Mischel PS. Altered cellular metabolism in gliomas - an emerging landscape of actionable co-dependency targets. Nat Rev Cancer. 2020; 20:57–70. https://doi.org/10.1038/s41568-019-0226-5 [PubMed]
3. Chen X, Yan Y, Zhou J, Huo L, Qian L, Zeng S, Li Z, Wei J, Xu Z, Gong Z. Clinical prognostic value of isocitrate dehydrogenase mutation, O-6-methylguanine-DNA methyltransferase promoter methylation, and 1p19q co-deletion in glioma patients. Ann Transl Med. 2019; 7:541. https://doi.org/10.21037/atm.2019.09.126 [PubMed]
4. Mazurek M, Litak J, Kamieniak P, Kulesza B, Jonak K, Baj J, Grochowski C. Metformin as potential therapy for high-grade glioma. Cancers (Basel). 2020; 12:210. https://doi.org/10.3390/cancers12010210 [PubMed]
5. Yan Y, Xu Z, Chen X, Wang X, Zeng S, Zhao Z, Qian L, Li Z, Wei J, Huo L, Li X, Gong Z, Sun L. Novel function of lncRNA ADAMTS9-AS2 in promoting temozolomide resistance in glioblastoma via upregulating the FUS/MDM2 ubiquitination axis. Front Cell Dev Biol. 2019; 7:217. https://doi.org/10.3389/fcell.2019.00217 [PubMed]
6. Yan Y, Xu Z, Dai S, Qian L, Sun L, Gong Z. Targeting autophagy to sensitive glioma to temozolomide treatment. J Exp Clin Cancer Res. 2016; 35:23. https://doi.org/10.1186/s13046-016-0303-5 [PubMed]
7. Shin HJ, Kim DN, Lee SK. Association between Epstein-Barr virus infection and chemoresistance to docetaxel in gastric carcinoma. Mol Cells. 2011; 32:173–79. https://doi.org/10.1007/s10059-011-0066-y [PubMed]
8. Tomono T, Kajita M, Yano K, Ogihara T. Adenovirus vector infection of non-small-cell lung cancer cells is a trigger for multi-drug resistance mediated by P-glycoprotein. Biochem Biophys Res Commun. 2016; 476:183–87. https://doi.org/10.1016/j.bbrc.2016.05.070 [PubMed]
9. Kreis NN, Louwen F, Yuan J. The multifaceted p21 (Cip1/Waf1/ CDKN1A) in cell differentiation, migration and cancer therapy. Cancers (Basel). 2019; 11:1220. https://doi.org/10.3390/cancers11091220 [PubMed]
10. Winters ZE, Leek RD, Bradburn MJ, Norbury CJ, Harris AL. Cytoplasmic p21WAF1/CIP1 expression is correlated with HER-2/ neu in breast cancer and is an independent predictor of prognosis. Breast Cancer Res. 2003; 5:R242–49. https://doi.org/10.1186/bcr654 [PubMed]
11. Liu YW, Xia R, Lu K, Xie M, Yang F, Sun M, De W, Wang C, Ji G. LincRNAFEZF1-AS1 represses p21 expression to promote gastric cancer proliferation through LSD1-mediated H3K4me2 demethylation. Mol Cancer. 2017; 16:39. https://doi.org/10.1186/s12943-017-0588-9 [PubMed]
12. Price JG, Idoyaga J, Salmon H, Hogstad B, Bigarella CL, Ghaffari S, Leboeuf M, Merad M. CDKN1A regulates langerhans cell survival and promotes treg cell generation upon exposure to ionizing irradiation. Nat Immunol. 2015; 16:1060–68. https://doi.org/10.1038/ni.3270 [PubMed]
13. Crane CA, Austgen K, Haberthur K, Hofmann C, Moyes KW, Avanesyan L, Fong L, Campbell MJ, Cooper S, Oakes SA, Parsa AT, Lanier LL. Immune evasion mediated by tumor-derived lactate dehydrogenase induction of NKG2D ligands on myeloid cells in glioblastoma patients. Proc Natl Acad Sci USA. 2014; 111:12823–28. https://doi.org/10.1073/pnas.1413933111 [PubMed]
14. Finocchiaro G. TLRgeting evasion of immune pathways in glioblastoma. Cell Stem Cell. 2017; 20:422–24. https://doi.org/10.1016/j.stem.2017.03.018 [PubMed]
15. Gomez GG, Kruse CA. Mechanisms of Malignant glioma immune resistance and sources of immunosuppression. Gene Ther Mol Biol. 2006; 10:133–46. [PubMed]
16. Muppa P, Parrilha Terra SB, Sharma A, Mansfield AS, Aubry MC, Bhinge K, Asiedu MK, de Andrade M, Janaki N, Murphy SJ, Nasir A, Van Keulen V, Vasmatzis G, et al. Immune cell infiltration may be a key determinant of long-term survival in small cell lung cancer. J Thorac Oncol. 2019; 14:1286–95. https://doi.org/10.1016/j.jtho.2019.03.028 [PubMed]
17. Choi BD, Yu X, Castano AP, Darr H, Henderson DB, Bouffard AA, Larson RC, Scarfò I, Bailey SR, Gerhard GM, Frigault MJ, Leick MB, Schmidts A, et al. CRISPR-Cas9 disruption of PD-1 enhances activity of universal EGFRvIII CAR T cells in a preclinical model of human glioblastoma. J Immunother Cancer. 2019; 7:304. https://doi.org/10.1186/s40425-019-0806-7 [PubMed]
18. Berghoff AS, Kiesel B, Widhalm G, Rajky O, Ricken G, Wöhrer A, Dieckmann K, Filipits M, Brandstetter A, Weller M, Kurscheid S, Hegi ME, Zielinski CC, et al. Programmed death ligand 1 expression and tumor-infiltrating lymphocytes in glioblastoma. Neuro Oncol. 2015; 17:1064–75. https://doi.org/10.1093/neuonc/nou307 [PubMed]
19. Li G, Wang Z, Zhang C, Liu X, Cai J, Wang Z, Hu H, Wu F, Bao Z, Liu Y, Zhao L, Liang T, Yang F, et al. Molecular and clinical characterization of TIM-3 in glioma through 1,024 samples. Oncoimmunology. 2017; 6:e1328339. https://doi.org/10.1080/2162402X.2017.1328339 [PubMed]
20. Dai S, Yan Y, Xu Z, Zeng S, Qian L, Huo L, Li X, Sun L, Gong Z. SCD1 Confers Temozolomide Resistance to Human Glioma Cells via the Akt/GSK3β/β-Catenin Signaling Axis. Front Pharmacol. 2018; 8:960. https://doi.org/10.3389/fphar.2017.00960 [PubMed]
21. Burton C, Das A, McDonald D, Vandergrift WA 3rd, Patel SJ, Cachia D, Bartee E. Oncolytic myxoma virus synergizes with standard of care for treatment of glioblastoma multiforme. Oncolytic Virother. 2018; 7:107–16. https://doi.org/10.2147/OV.S179335 [PubMed]
22. Foster H, Ulasov IV, Cobbs CS. Human cytomegalovirus-mediated immunomodulation: effects on glioblastoma progression. Biochim Biophys Acta Rev Cancer. 2017; 1868:273–76. https://doi.org/10.1016/j.bbcan.2017.05.006 [PubMed]
23. Rajapakse VN, Luna A, Yamade M, Loman L, Varma S, Sunshine M, Iorio F, Sousa FG, Elloumi F, Aladjem MI, Thomas A, Sander C, Kohn KW, et al. CellMinerCDB for integrative cross-database genomics and pharmacogenomics analyses of cancer cell lines. iScience. 2018; 10:247–64. https://doi.org/10.1016/j.isci.2018.11.029 [PubMed]
24. Majewska E, Szeliga M. AKT/GSK3β signaling in glioblastoma. Neurochem Res. 2017; 42:918–24. https://doi.org/10.1007/s11064-016-2044-4 [PubMed]
25. Guglietta S, Rescigno M. Hypercoagulation and complement: connected players in tumor development and metastases. Semin Immunol. 2016; 28:578–86. https://doi.org/10.1016/j.smim.2016.10.011 [PubMed]
26. Fulop GA, Kiss T, Tarantini S, Balasubramanian P, Yabluchanskiy A, Farkas E, Bari F, Ungvari Z, Csiszar A. Nrf2 deficiency in aged mice exacerbates cellular senescence promoting cerebrovascular inflammation. Geroscience. 2018; 40:513–21. https://doi.org/10.1007/s11357-018-0047-6 [PubMed]
27. Kiss T, Giles CB, Tarantini S, Yabluchanskiy A, Balasubramanian P, Gautam T, Csipo T, Nyúl-Tóth Á, Lipecz A, Szabo C, Farkas E, Wren JD, Csiszar A, Ungvari Z. Nicotinamide mononucleotide (NMN) supplementation promotes anti-aging miRNA expression profile in the aorta of aged mice, predicting epigenetic rejuvenation and anti-atherogenic effects. Geroscience. 2019; 41:419–39. https://doi.org/10.1007/s11357-019-00095-x [PubMed]
28. Dutto I, Tillhon M, Cazzalini O, Stivala LA, Prosperi E. Biology of the cell cycle inhibitor p21(CDKN1A): molecular mechanisms and relevance in chemical toxicology. Arch Toxicol. 2015; 89:155–78. https://doi.org/10.1007/s00204-014-1430-4 [PubMed]
29. Kim EM, Jung CH, Kim J, Hwang SG, Park JK, Um HD. The p53/p21 complex regulates cancer cell invasion and apoptosis by targeting Bcl-2 family proteins. Cancer Res. 2017; 77:3092–100. https://doi.org/10.1158/0008-5472.CAN-16-2098 [PubMed]
30. Kim J, Bae S, An S, Park JK, Kim EM, Hwang SG, Kim WJ, Um HD. Cooperative actions of p21WAF1 and p53 induce slug protein degradation and suppress cell invasion. EMBO Rep. 2014; 15:1062–68. https://doi.org/10.15252/embr.201438587 [PubMed]
31. Cheng X, Xia W, Yang JY, Hsu JL, Chou CK, Sun HL, Wyszomierski SL, Mills GB, Muller WJ, Yu D, Hung MC. Activation of p21(CIP1/WAF1) in mammary epithelium accelerates mammary tumorigenesis and promotes lung metastasis. Biochem Biophys Res Commun. 2010; 403:103–07. https://doi.org/10.1016/j.bbrc.2010.10.126 [PubMed]
32. Huang Y, Wang W, Chen Y, Huang Y, Zhang J, He S, Tan Y, Qiang F, Li A, Røe OD, Wang S, Zhou Y, Zhou J. The opposite prognostic significance of nuclear and cytoplasmic p21 expression in resectable gastric cancer patients. J Gastroenterol. 2014; 49:1441–52. https://doi.org/10.1007/s00535-013-0900-4 [PubMed]
33. Kreis NN, Friemel A, Ritter A, Roth S, Rolle U, Louwen F, Yuan J. Function of p21 (Cip1/Waf1/ CDKN1A) in migration and invasion of cancer and trophoblastic cells. Cancers (Basel). 2019; 11:989. https://doi.org/10.3390/cancers11070989 [PubMed]
34. Li Y, Dowbenko D, Lasky LA. AKT/PKB phosphorylation of p21Cip/WAF1 enhances protein stability of p21Cip/WAF1 and promotes cell survival. J Biol Chem. 2002; 277:11352–61. https://doi.org/10.1074/jbc.M109062200 [PubMed]
35. Zhou BP, Liao Y, Xia W, Spohn B, Lee MH, Hung MC. Cytoplasmic localization of p21Cip1/WAF1 by Akt-induced phosphorylation in HER-2/neu-overexpressing cells. Nat Cell Biol. 2001; 3:245–52. https://doi.org/10.1038/35060032 [PubMed]
36. Yan Y, Liu W, Liu M, Gong Z, Xu Z. Immune cell infiltration influences long-term survivorship of patients with SCLC. J Thorac Oncol. 2019; 14:e241. https://doi.org/10.1016/j.jtho.2019.06.017 [PubMed]
37. Santín-Márquez R, Alarcón-Aguilar A, López-Diazguerrero NE, Chondrogianni N, Königsberg M. Sulforaphane - role in aging and neurodegeneration. Geroscience. 2019; 41:655–70. https://doi.org/10.1007/s11357-019-00061-7 [PubMed]
38. Royce GH, Brown-Borg HM, Deepa SS. The potential role of necroptosis in inflammaging and aging. Geroscience. 2019; 41:795–811. https://doi.org/10.1007/s11357-019-00131-w [PubMed]
39. Rožman P. The potential of non-myeloablative heterochronous autologous hematopoietic stem cell transplantation for extending a healthy life span. Geroscience. 2018; 40:221–42. https://doi.org/10.1007/s11357-018-0027-x [PubMed]
40. Yan S, Zhao P, Yu T, Gu N. Current applications and future prospects of nanotechnology in cancer immunotherapy. Cancer Biol Med. 2019; 16:486–97. https://doi.org/10.20892/j.issn.2095-3941.2018.0493 [PubMed]
41. Yan Y, Chen X, Wei J, Gong Z, Xu Z. Immunotherapy combinations in patients with small cell lung cancers. J Thorac Oncol. 2019; 14:e244–45. https://doi.org/10.1016/j.jtho.2019.05.021 [PubMed]
42. Erhart F, Buchroithner J, Reitermaier R, Fischhuber K, Klingenbrunner S, Sloma I, Hibsh D, Kozol R, Efroni S, Ricken G, Wöhrer A, Haberler C, Hainfellner J, et al. Immunological analysis of phase II glioblastoma dendritic cell vaccine (audencel) trial: immune system characteristics influence outcome and audencel up-regulates Th1-related immunovariables. Acta Neuropathol Commun. 2018; 6:135. https://doi.org/10.1186/s40478-018-0621-2 [PubMed]
43. Takashima Y, Kawaguchi A, Hayano A, Yamanaka R. CD276 and the gene signature composed of GATA3 and LGALS3 enable prognosis prediction of glioblastoma multiforme. PLoS One. 2019; 14:e0216825. https://doi.org/10.1371/journal.pone.0216825 [PubMed]
44. Heiland DH, Haaker G, Delev D, Mercas B, Masalha W, Heynckes S, Gäbelein A, Pfeifer D, Carro MS, Weyerbrock A, Prinz M, Schnell O. Comprehensive analysis of PD-L1 expression in glioblastoma multiforme. Oncotarget. 2017; 8:42214–25. https://doi.org/10.18632/oncotarget.15031 [PubMed]
45. DiDomenico J, Lamano JB, Oyon D, Li Y, Veliceasa D, Kaur G, Ampie L, Choy W, Lamano JB, Bloch O. The immune checkpoint protein PD-L1 induces and maintains regulatory T cells in glioblastoma. Oncoimmunology. 2018; 7:e1448329. https://doi.org/10.1080/2162402X.2018.1448329 [PubMed]
46. Huang G, Ho B, Conroy J, Liu S, Qiang H, Golubovskaya V. The microarray gene profiling analysis of glioblastoma cancer cells reveals genes affected by FAK inhibitor Y15 and combination of Y15 and temozolomide. Anticancer Agents Med Chem. 2014; 14:9–17. https://doi.org/10.2174/18715206113139990141 [PubMed]
47. Wu L, Bernal GM, Cahill KE, Pytel P, Fitzpatrick CA, Mashek H, Weichselbaum RR, Yamini B. BCL3 expression promotes resistance to alkylating chemotherapy in gliomas. Sci Transl Med. 2018; 10:eaar2238. https://doi.org/10.1126/scitranslmed.aar2238 [PubMed]
48. Barrett T, Edgar R. Reannotation of array probes at NCBI’s GEO database. Nat Methods. 2008; 5:117. https://doi.org/10.1038/nmeth0208-117b [PubMed]
49. Rappaport N, Twik M, Plaschkes I, Nudel R, Iny Stein T, Levitt J, Gershoni M, Morrey CP, Safran M, Lancet D. MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search. Nucleic Acids Res. 2017; 45:D877–87. https://doi.org/10.1093/nar/gkw1012 [PubMed]
50. Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, Mouradov D, Sieber OM, Simpson RJ, Salim A, Bacic A, Hill AF, Stroud DA, et al. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015; 15:2597–601. https://doi.org/10.1002/pmic.201400515 [PubMed]
51. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017; 45:W98–102. https://doi.org/10.1093/nar/gkx247 [PubMed]
52. Bowman RL, Wang Q, Carro A, Verhaak RG, Squatrito M. GlioVis data portal for visualization and analysis of brain tumor expression datasets. Neuro Oncol. 2017; 19:139–41. https://doi.org/10.1093/neuonc/now247 [PubMed]
53. Díez-Villanueva A, Mallona I, Peinado MA. Wanderer, an interactive viewer to explore DNA methylation and gene expression data in human cancer. Epigenetics Chromatin. 2015; 8:22. https://doi.org/10.1186/s13072-015-0014-8 [PubMed]
54. Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson Å, Kampf C, Sjöstedt E, Asplund A, Olsson I, Edlund K, Lundberg E, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015; 347:1260419. https://doi.org/10.1126/science.1260419 [PubMed]
55. Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehár J, Kryukov GV, Sonkin D, Reddy A, Liu M, Murray L, et al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012; 483:603–07. https://doi.org/10.1038/nature11003 [PubMed]
56. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, Jensen LJ, von Mering C. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017; 45:D362–68. https://doi.org/10.1093/nar/gkw937 [PubMed]
57. Wu G, Zhao Z, Yan Y, Zhou Y, Wei J, Chen X, Lin W, Ou C, Li J, Wang X, Xiong K, Zhou J, Xu Z. CPS1 expression and its prognostic significance in lung adenocarcinoma. Ann Transl Med. 2020; 8:341. https://doi.org/10.21037/atm.2020.02.146 [PubMed]
58. Reimand J, Isserlin R, Voisin V, Kucera M, Tannus-Lopes C, Rostamianfar A, Wadi L, Meyer M, Wong J, Xu C, Merico D, Bader GD. Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, cytoscape and EnrichmentMap. Nat Protoc. 2019; 14:482–517. https://doi.org/10.1038/s41596-018-0103-9 [PubMed]
59. Wu G, Yan Y, Zhou Y, Wang X, Wei J, Chen X, Lin W, Ou C, Zhou J, Xu Z. Expression and clinical significance of CPS1 in glioblastoma multiforme. Curr Res Transl Med. 2019; 67:123–28. https://doi.org/10.1016/j.retram.2019.08.003 [PubMed]
60. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009; 4:44–57. https://doi.org/10.1038/nprot.2008.211 [PubMed]
61. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SS, Wu WC, Chu KC, Wong CY, Lau CY, Chen I, Chan NW, Zhang J. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019; 35:4200–02. https://doi.org/10.1093/bioinformatics/btz210 [PubMed]
62. Li J, Li L, Li Z, Gong G, Chen P, Liu H, Wang J, Liu Y, Wu X. The role of miR-205 in the VEGF-mediated promotion of human ovarian cancer cell invasion. Gynecol Oncol. 2015; 137:125–33. https://doi.org/10.1016/j.ygyno.2015.01.531 [PubMed]
63. Li J, Alvero AB, Nuti S, Tedja R, Roberts CM, Pitruzzello M, Li Y, Xiao Q, Zhang S, Gan Y, Wu X, Mor G, Yin G. CBX7 binds the E-box to inhibit TWIST-1 function and inhibit tumorigenicity and metastatic potential. Oncogene. 2020; 39:3965–79. https://doi.org/10.1038/s41388-020-1269-5 [PubMed]

Research Paper Volume 13, Issue 5 pp 6662—6680