<i>CKAP2L</i>, as an Independent Risk Factor, Closely Related to the Prognosis of Glioma

Zhu, Li; Zheng, Yanlei; Hu, Ronghua; Hu, Chenchen

doi:https://doi.org/10.1155/2021/5486131

BioMed Research International

On this page

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

Research Article | Open Access

Volume 2021 | Article ID 5486131 | https://doi.org/10.1155/2021/5486131

CKAP2L, as an Independent Risk Factor, Closely Related to the Prognosis of Glioma

Li Zhu,¹Yanlei Zheng,¹Ronghua Hu,¹and Chenchen Hu¹

Academic Editor: Jafri Malin Abdullah

Received24 Aug 2020

Revised02 Aug 2021

Accepted05 Aug 2021

Published28 Sept 2021

Abstract

Recent studies have found that cytoskeleton-associated protein 2 like (CKAP2L), an important oncogene, is involved in the biological behavior of many malignant tumors, but its function in the malignant course of glioma has not been confirmed. The main purpose of this study was to clarify the relationship between prognostic clinical characteristics of glioma patients and CKAP2L expression using data collected from the GEPIA, HPA, CGGA, TCGA, and GEO databases. CKAP2L expression was significantly increased in glioma. Further, Kaplan-Meier plots revealed that increased expression of CKAP2L was associated with shorter survival time of glioma patients in datasets retrieved from multiple databases. Cox regression analysis indicated that CKAP2L can serve as an independent risk factor but also has relatively reliable diagnostic value for the prognosis of glioma patients. The results of gene set enrichment analysis suggested that CKAP2L may play a regulatory role through the cell cycle, homologous recombination, and N-glycan biosynthesis cell signaling pathways. Several drugs with potential inhibitory effects on CKAP2L were identified in the CMap database that may have therapeutic effects on glioma. Finally, knockdown of CKAP2L inhibited the proliferation and invasion of cells by reducing the expression level of cell cycle-related proteins. This is the first study to demonstrate that high CKAP2L expression leads to poor prognosis in glioma patients, providing a novel target for diagnosis and treatment of glioma.

1. Introduction

Gliomas are primary malignant tumors derived from neural progenitor cells and glial cells that account for the vast majority of intracranial tumors [1, 2]. Because of the high degree of malignancy and invasive growth of glioma cells, gliomas lead to severe physiological dysfunction in patients and eventual death, reducing patient quality of life and increasing the burden on the healthcare system [3]. To improve the prognosis of patients with glioma, many scientists have made unremitting efforts to form the standardized therapeutic principles that have been developed to maximize tumor resection, reduce tumor load, and administer adjuvant radiotherapy and chemotherapy [4]. Despite these efforts, the prognosis of glioma patients has only partially improved, and survival rates remain poor [3]. One of the most important factors affecting the prognosis of glioma patients may be the lack of a sensitive and highly specific biomarker for targeted therapy.

During the process of tumorigenesis and development, a large number of tumor suppressor genes are known to be inactivated, whereas oncogenes are activated. The abnormal expression of these genes plays an extremely important role in regulating the growth of tumor cells, including gliomas [5]. With the development of molecular medicine and a deeper understanding of tumor pathological mechanisms at the molecular level, many of these abnormally expressed genes have been used as biomarkers in the diagnosis and treatment of tumors. In an effort to identify more reliable biomarkers, researchers have established a variety of publicly available databases, including TCGA and CGGA, with sufficient sample sizes and comprehensive clinical information. Mining this data provides guidance for molecular typing and drug target development for glioma, laying the foundation for development of a whole chain of precision medicine for glioma [6, 7]. To date, many effective biomarkers have been identified through database analysis methods for the diagnosis and targeted treatment of glioma, including CD276, HLA-F, and CMTM6 [8–10]. However, the complex pathological mechanisms of glioma may have resulted in many potential targets remaining unexplored, leading to lack of progress in improving the prognosis of glioma patients.

Cytoskeleton-associated protein 2 like (CKAP2L) plays an important regulatory role in the mitochondrial spin of neural progenitor cells. CKAP2L mutation is related to defects in the spindle tissue, including mitotic spindle defects, chromosome hysteresis, and other mitotic instability characteristics, that are involved in the formation and development of cancer [11]. Furthermore, a large number of previous studies have confirmed that CKAP2L is involved in the regulation of a variety of cancers. For example, Xiong et al. determined that abnormal expression of CKAP2L was highly correlated with poor prognosis of lung cancer and may be regulated by the MAPK signaling pathway [12]. Moreover, Wang and He reported that increased expression of CKAP2L in hepatocellular carcinoma (HCC) cells promoted their migration and invasion, which was associated with poor prognosis of HCC [13]. In addition, CKAP2L is reportedly involved in the occurrence and development of other malignant tumors, such as breast and prostate cancers, and is associated with poor prognosis of these cancers [14, 15]. However, it is worth noting that the relationship between CKAP2L and the prognostic features of gliomas has not been elucidated and revealed.

Therefore, this study was aimed at exploring the previously unknown relationship between CKAP2L and the biological characteristics of glioma patients, based on clinical sample analysis of multiple databases, and at confirming the oncogenic role of CKAP2L in glioma. We believe that the results of this study will provide a new target for the diagnosis and treatment of glioma, thus improving the prognosis of glioma patients.

2. Materials and Methods

2.1. Data Collection

Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia.cancer-pku.cn/) is a public data platform containing gene expression levels and patient characteristics for the vast majority of human tumors [16]. Using this database, we collected CKAP2L expression data for 163 glioblastomas, 518 low-grade gliomas, and 207 normal brain tissue samples.

Gene Expression Omnibus (GEO) is a public repository of functional genomic data obtained via high-throughput and next-generation sequencing [17], which has been widely used for bioinformatics analysis of tumor characteristics. We retrieved three glioma microarray datasets from the GEO database. GSE4290, comprised mainly of American samples, contained 23 normal brain samples and 77 glioma samples. GSE50161, comprised mainly of American samples, contained 13 normal brain samples and 34 glioma samples. GSE116520, comprised mainly of Indian samples, contained eight normal brain samples and 34 glioma samples.

The Cancer Genome Atlas (TCGA; http://cancergenome.nih.gov/) is an internationally renowned public data platform containing large-scale genome sequencing analysis of more than 30 human tumors from mainly white and black races [18] and has conducted large-scale genome sequencing analysis, which contains the main white and black races in the database. We obtained mRNA sequences and clinical information for 653 patients with glioma from TCGA database (Table S1).

The Chinese Glioma Genome Atlas (CGGA) is a public data platform focusing on brain tumor research that contains various types of data for more than 2000 tissue samples from Chinese patients. Excluding samples with incomplete clinical information, we obtained an RNA-seq dataset of 748 glioma samples and RNA-microarray data from 268 glioma samples (Table S2 and Table S3) [19].

According to the 2016 World Health Organization Classification of Tumors of the Central Nervous System, we obtained 40 glioma tissue samples and 5 nontumor brain tissue samples from patients at the Hubei Cancer Hospital, including 17 cases of grade III glioma and 23 cases of grade IV glioma. The clinical and molecular characteristics of the 40 glioma patients and their corresponding follow-up characteristics were collected (Table S4). The tissue samples were stored at 80°C for total RNA extraction. The study protocol was performed according to the Declaration of Helsinki and was approved by the Ethics Committee of Hubei Cancer Hospital.

2.2. Gene Set Enrichment Analysis for CKAP2L

Gene set enrichment analysis (GSEA) is a commonly used biological information tool [20] that can simultaneously analyze the expression profiles of multiple genes, enabling accurate identification of a target gene’s function mechanism. In this study, we divided the data obtained from TCGA and CGGA databases into two groups according to expression levels of CKAP2L. GSEA 4.0. jar software was employed to analyze the cell signaling pathways associated with CKAP2L. KEGG pathways were selected after 1000 iterations.

2.3. Cell Culture

Normal human astrocytes and glioma cell lines LN229, T98, U251, and A172 were purchased from Procell Life Science & Technology Co, Ltd. (Hyderabad, India). All cells were incubated in Dulbecco’s modified Eagle’s medium (Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (Hyclone, Logan, UT, USA) at 37°C under 5% CO₂ in a humidified incubator. Expression of CKAP2L was measured in all five cell lines, as described below, and the cell line with the highest CKAP2L expression was selected for further experiments. Accordingly, U251 cells were transfected with lentiviruses containing control shRNA or shRNA against CKAP2L (5-GATCCGCAAACAAAGAGAACTTGCTCGATATTTCAAGAGAATATCGAGCAAGTTCTCTTTGTTTGTTTTTTGGAAG-3), which were purchased from GenePharma (Suzhou, China), and screened using puromycin.

2.4. RT-qPCR Analysis

A Total RNA Kit I (Omega Bio-tek, Norcross, GA, USA) was used to extract the total RNA from cells and brain samples, and a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) was used to confirm the mass of total extracted RNA. Then, cDNA was reverse transcribed from total RNA using NovoScript Plus All-in-one 1st Strand cDNA Synthesis SuperMix (Novoprotein, Shanghai, China). The expression of CKAP2L was analyzed using RT-qPCR. The primer sequences were as follows: forward 5-GGGAAAACTGAAGAGCCAAAACA-3 and reverse 5-AGGTTTGACAGGCAAAACAACA-3. Expression levels were calculated using the 2^−ΔΔCT method.

2.5. MTT Assay

U251 cells were plated at a density of cells/well in 96-well cell culture plates after transfection with shRNA for 24 h. Cells were treated with 20 μL MTT reagent (5 mg/mL) at selected times (0, 24, 48, 72, and 96 h) and then incubated for 4 h. Subsequently, 150 μL DMSO was used to dissolve the formazan crystals in each well, followed by incubation for 15 min. Finally, the absorbance was measured at 490 nm using a microplate reader.

2.6. Wound Healing Assay

U251 cells were plated at a density of cells/well in 6-well cell culture plates after transfection with shRNA for 24 h. After incubation for another 24 h at 37°C under 5% CO₂, the U251 cells formed a monolayer, which was linearly scratched with a sterile pipette tip. The detached cells were washed off, and serum-free cell culture medium was added to the treated wells. Cell migration into the wound was measured after incubation for 24 h, and the distance was calculated using ImageJ software (National Institutes of Health, Bethesda, MD, USA).

2.7. Western Blotting

Total protein was extracted from U251 cells treated with shRNA, the proteins were then fractionated by 12% SDS-PAGE and transferred onto a PVDF membrane, and the membrane was blocked with 5% nonfat milk. Subsequently, the membrane was incubated with primary specific antibodies for cell cycle-related proteins, including CDK4 (1 : 1000, Proteintech, Chicago, IL, USA), CCNB1 (1 : 1000, Proteintech), CCND1 (1 : 1000, Proteintech), CCNE1 (1 : 1000, Proteintech), and GAPDH (1 : 1000, Proteintech), followed by three washes with TBST for 5 min. Then, horseradish peroxidase-conjugated secondary antibodies were applied for 1 h at room temperature, followed by three washes with TBST for 10 min. Finally, enhanced chemiluminescence was used to detect the cell cycle-related proteins, and signals were captured by Image Lab software.

2.8. CMap

The Connectivity Map (CMap; https://portals.broadinstitute.org/cmap/) database is used to compare changes in gene expression profiles of cultured human cells treated with small-molecule drugs to enable prediction of drugs to treat target diseases [21]. We used Pearson correlation coefficients to obtain the top 10 genes positively (upregulated) and negatively (downregulated) correlated with CKAP2L, based on the RNA-seq dataset of 748 glioma samples obtained from the CGGA database. Potential small-molecule drugs with an inhibitory effect on CKAP2L were identified in CMap when and . Finally, the 3D and 2D structures of these drugs and their molecular formulas were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [22].

2.9. Statistical Analysis

All data were statistically analyzed using R language software (version 3.6.1). The comparison of CKAP2L expression in tumor and normal brain tissues was performed using the limma package in R software according to the cut-off standard ( and ). Unpaired -tests were applied to analyze the difference in CKAP2L expression between the two groups, with considered significant. Kaplan Meier curves and Cox regression analysis (univariate and multivariate) were used to determine the impact of CKAP2L on the prognosis and diagnostic value of glioma patients. Pearson correlation coefficients were used to identify genes that were coexpressed with CKAP2L. values < 0.05 were considered statistically significant.

3. Results

3.1. CKAP2L Expression Correlated with Overall Survival in Glioma

GEPIA database analysis revealed abnormal CKAP2L expression in human malignancies. The expression level of CKAP2L in acute myeloid leukemia (LAML) was significantly lower than that in the normal control group, but expression of CKAP2L in most malignant tumors was abnormally increased, including in glioblastoma (GBM), as shown in Figure 1(a). Analysis of three different glioma-related datasets in the GEO database, including dozens of normal brain samples and hundreds of glioma samples, revealed that the expression level of CKAP2L in glioma tissues was significantly higher () than that in normal brain tissue, as shown in Figures 1(b)–1(d).

(a)

(b)

(c)

(d)

To further study the relationship between CKAP2L expression and prognosis of glioma, we analyzed three transcriptome datasets: CGGA microarray, CGGA RNA-seq, and TCGA RNA-seq. As shown in Figures 2(a)–2(c), increased expression of CKAP2L in the three datasets was associated with decreased patient overall survival (). The results were consistent in thousands of tissue samples, reliably indicating that high expression levels of CKAP2L led to poor prognosis in glioma patients.

(a)

(b)

(c)

Receiver operating characteristic (ROC) curve analysis was employed to determine whether CKAP2L has diagnostic value in the prognosis of glioma patients (Figure 3). The area under the ROC curve for one, three, and five years using data obtained from the three transcriptome datasets was >0.7, indicating that CKAP2L has moderate diagnostic value, excluding the CGGA RNA-seq dataset for 1 year (Figures 4(a)–4(c)). These results supported that CKAP2L has a diagnostic value for glioma patients with poor prognosis.

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

3.2. CKAP2L Can Be an Independent Risk Factor of Glioma Patients

To investigate whether CKAP2L can be used as an independent risk factor for poor prognosis of glioma, we analyzed the three transcriptome datasets by univariate and multivariate Cox analyses. The univariate Cox analysis results indicated that increased CKAP2L expression could be a risk factor for poor prognosis in glioma patients in the CGGA RNA-seq (; -1.973, ), CGGA microarray (; -2.140, ), and TCGA RNA-seq (; -1.305, ) datasets. Increased tumor grade in patients with glioma was also a risk factor for poor prognosis in the CGGA RNA-seq (; 95% -3.291, ), CGGA microarray (; -3.100, ), and TCGA RNA-seq (; -5.760, ) datasets. Furthermore, age was also a risk factor for poor prognosis in glioma patients in the CGGA RNA-seq (; -1.960, ), CGGA microarray (; -2.349, ), and TCGA RNA-seq (; -1.083, ) datasets (Figures 3(a), 3(c), and 3(e)).

At the same time, multivariate Cox analysis revealed that increased expression of CKAP2L could be an independent risk factor for poor prognosis in glioma patients in the CGGA RNA-seq (; -1.431, ), CGGA microarray (; -1.671, ), and TCGA RNA-seq (; -1.187, ) datasets. Tumor grade in the CGGA RNA-seq (; -3.398, ), CGGA microarray (; -3.892, ), and TCGA RNA-seq (; -3.721, ) datasets was found to be an independent risk factor for glioma patients (Figures 3(b), 3(d), and 3(f)). Further, the PRS type was an independent risk factor for poor prognosis of glioma patients in the CGGA RNA-seq (; -2.275, ) and CGGA microarray (; -2.086, ) datasets, as shown in Figures 3(b) and 3(d). The above results indicated that CKAP2L can be an independent risk factor associated with poor prognosis of patients with glioma.

3.3. GSEA Predicts Signaling Pathways Associated with CKAP2L and Glioma

GSEA was used to analyze the potential cellular signaling pathways associated with CKAP2L. Analysis of three transcriptome datasets revealed that increased expression of CKAP2L can lead to activation of the p53, cell cycle, homologous recombination, and N-glycan biosynthesis cell signaling pathways (Figure 5 and Table 1). These results indicated that CKAP2L may regulate the occurrence and development of glioma through the above cell signaling pathways.

(a)

(b)

(c)

(d)

3.4. CKAP2L Expression Correlated with Tumor Grade

To further verify the CKAP2L expression profile in glioma, we examined CKAP2L protein levels in glioma data from The Human Protein Atlas, revealing markedly higher CKAP2L protein levels in low-grade and high-grade glioma tissues than in normal brain tissue (Figure 6).

(a)

(b)

(c)

(d)

(e)

(f)

3.5. Validation of Carcinogenic Effect of CKAP2L

We further validated the influence of CKAP2L on patients with glioma by analyzing 40 clinical samples, consisting of both high-and low-grade glioma tissues, and the patient characteristics and outcomes. CKAP2L expression in glioma was significantly increased compared with that in normal brain tissue (Figure 7(a)). Further, high expression of CKAP2L led to shorter survival times for both high-and low-grade glioma patients (Figures 7(c) and 7(d)). Moreover, we examined the expression of CKAP2L in normal human astrocytes and four different glioma cell lines, determining the highest CKAP2L expression in U251 cells (Figure 7(b)). Therefore, we verified that the increased expression level of CKAP2L in glioma in the database resulted in the decrease of overall survival time of glioma patients. Based on the highest expression level of CKAP2L of U251 among the four glioma cell lines, we chose U251 cells to perform subsequent in vitro experiments.

(a)

(b)

(c)

(d)

Using shRNA to knock down the expression of CKAP2L, we then performed MTT and wound healing assays to investigate the effect of CKAP2L on the biological characteristics of glioma cells. As shown in Figures 8(a) and 8(b), downregulation of CKAP2L remarkably inhibited the proliferation and migration of U251 cells. Furthermore, knockdown of CKAP2L significantly reduced the expression of cell cycle-related genes (Figure 8(c)). Therefore, the high expression of CKAP2L may promote the proliferation of glioma cells through the activation of cell cycle signal pathway, which has been further verified.

(a)

(b)

(c)

3.6. Gene Coexpression with CKAP2L and Drugs Targeting CKAP2L

Applying Pearson correlation analysis to the RNA-seq dataset of 748 glioma samples obtained from the CGGA database, we identified 10 genes positively correlated with CKAP2L (BUB1, TTK, KIF14, MELK, ASPM, KIF4A, NCAPG, BUB1B, ESCO2, and SGOL1) and 10 genes negatively correlated with CKAP2L (LDHD, FBXW4, AIFM3, ETNPPL, NOXA1, KNDC1, CYP46A1, LYNX1, MATK, and MYH7B) (Figure 9). Genes positively and negatively correlated with CKAP2L were regarded as upregulated and downregulated, respectively. After screening the genes in the CMap database to find drugs with targeted inhibitory effects on CKAP2L, four potentially therapeutic drugs were identified based on the search criteria, including camptothecin, mercaptopurine, piperidolate, and sanguinarine (Table 2). Finally, we obtained the 2D and 3D structures and molecular formulas of these drugs from PubChem, as shown in Figure 10.

(a)

(b)

(a)

(b)

(c)

(d)

4. Discussion

A large body of literature has confirmed that CKAP2L is involved in regulating many types of malignant tumors, promoting the proliferation of tumor cells, and shortening the survival time of patients. For example, Xu et al. reported that CKAP2L is involved in the biological behavior of prostate cancer and has a significant relationship with its malignant clinical features [15]. Similar conclusions have been reached for other malignant diseases such as hepatocellular carcinoma and lung cancer [12, 13]. However, to date, CKAP2L has not been reportedly associated with glioma.

To determine whether CKAP2L expression differs in glioma, we first investigated the expression of CKAP2L in glioma-related datasets obtained from GEPIA, HPA, and GEO. The results demonstrated that CKAP2L expression was higher in glioma tissues than in normal brain tissues at both the mRNA (Figure 1) and protein (Figure 6) levels. We employed a variety of data fusion analyses, based on thousands of multiethnic tissue samples to improve the scientific rigor and credibility of our results, and further verified the expression of CKAP2L in collected clinical samples and commercial glioma cell lines (Figure 7).

To confirm the impact of abnormally high expression of CKAP2L on the prognosis of glioma patients, we collected patient information from three different TCGA and CGGA transcriptome datasets and 40 clinical glioma patients. Kaplan-Meier curves demonstrated that increased CKAP2L expression was closely associated with poor prognosis of glioma patients. Concurring with our results, studies with liver, prostate, and lung cancers have reported that high expression of CKAP2L is associated with poor prognosis in tumor patients and can be used as a marker for diagnosis and treatment [13–15, 23]. Besides, the clinical molecular characteristics of glioma are closely related to the prognosis of patients. For example, some studies suggest that with the increase of glioma tumor grade, the heterogeneity of tumor cells also increases, and it has a negative correlation with the prognosis of patients [24]. Secondly, due to the infiltrative growth of glioma cells, the recurrence and secondary of glioma are the most important reasons leading to poor prognosis [25] and the incidence rate and mortality of glioma are positively correlated with the age of the patients [24]. Finally, molecular characteristics are closely related to the prognosis of gliomas, especially 1p19q codeletion and IDH mutation which are a good indicator for the prognosis of gliomas [26]. More importantly, in the WHO glioma classification in 2016, these two molecules have been used as a reference indicator of the malignant degree of glioma, and their different states will lead to a greater difference in the prognosis of patients [27]. Therefore, the above reports are consistent with the results of our multivariate analysis, which show that grade, recurrence, secondary, and age were risk factors for the prognosis of patients with glioma, and 1p19q codeletion and IDH mutation had a protective effect on the prognosis of patients. In addition to the above clinical and molecular characteristics of glioma patients, the status of radiotherapy and chemotherapy in our univariate analysis tends to be a risk factor, but the multivariate analysis is actually a protective factor for the prognosis of patients. We speculate that the prognosis of patients with low-grade gliomas is better, but they generally do not always receive chemoradiotherapy. However, patients with glioblastoma generally use chemoradiotherapy, so there is a risk factor in univariate analysis. Therefore, in the results of univariate analysis, the chemoradiotherapy group showed a risk factor. Our multivariate analysis is consistent with previous literature reports, which shows that the prognosis of glioma patients treated with radiotherapy and chemotherapy can improve the prognosis of patients [28]. Therefore, our study findings support that CKAP2L is associated with poor prognosis of glioma, can be used as an independent risk factor, and has diagnostic value.

Although the results of this current study confirmed that CKAP2L is a novel oncogenic gene that significantly reduces the survival time of patients, the molecular mechanism underlying the biological function of CKAP2L remains unknown. GSEA is a scientific, reliable, and widely used method in the biomedical field [19]. Here, GSEA of the above three datasets suggested that CKAP2L may participate in the p53, cell cycle, homologous recombination, and N-glycan biosynthesis signaling pathways. p53 is known to regulate a variety of signaling pathways, including the tumor cell cycle, autophagy, apoptosis, and aging. For example, Zhang et al. demonstrated that p53 can regulate the cell cycle to inhibit pancreatic cancer cell proliferation [29]. Research on the role of N-glycan in tumors has recently attracted increasing attention. N-Glycan can reportedly be used as a diagnostic marker for germline stem cell tumors and is associated with poor prognosis of patients with tumors [30]. Further, studies have reported that blocking N-polysaccharide precursor biosynthesis can inhibit the growth of U87 glioblastoma cells and significantly inhibit glioma size in vivo [31]. GSEA in the current study indicated that CKAP2L may be involved in the abovementioned cell signaling pathways and our subsequent in vitro experiments determined that knockdown of CKAP2L decreased the expression of key cell cycle-related proteins. The MTT and wound healing assays further demonstrated that CKAP2L influenced the proliferative and migratory abilities of glioma cells.

Through gene correlation analysis, we identified genes positively correlated with CKAP2L, including BUB1, TTK, and ASPM, which can promote the development of glioma and treatment resistance. These genes also indirectly support the hypothesis that CKAP2L may promote the pathological process of glioma. For example, BUB1 is a mitotic checkpoint that can promote the proliferation of glioma cells and radioactive resistance and is directly associated with poor prognosis of glioma patients [32]. TTK reportedly plays a role in promoting the growth of gliomas. Indeed, silencing TTK can inhibit the proliferation, invasion, and radiation resistance of glioma cells and other malignant behaviors [33]. CMap analysis predicted four drugs that could negatively regulate CKAP2L. The therapeutic effects of two of these drugs on glioma have been previously reported, including camptothecin, a nanoprodrug with inhibitory effects on U87-MG glioma cell proliferation [34], and sanguinarine, which can reportedly induce ROS-dependent ERK1/2 activation and autophagy cell death in glioma cells [35]. The above findings support that CKAP2L may act synergistically with these identified genes, which may lead to poor prognosis of glioma patients, and the drugs identified via CMap may potentially inhibit CKAP2L and treat glioma.

Although we made every effort to elucidate the relationship between CKAP2L and the prognosis of glioma patients, the use of public databases has certain limitations, such as inconsistent clinical information across different databases. Therefore, our analysis does not include all clinical treatment details and individualized efficacy. However, we must emphasize that public databases also enable access to sufficient numbers of research samples in a short period of time and encompass a variety of ethnicities, which would be difficult to establish within a clinical study. Therefore, the reliability and authenticity of our analysis results were greatly improved through the use of multiple public databases.

5. Conclusions

The study results demonstrate that CKAP2L expression increases with the grade of glioma and is associated with poor prognosis of glioma patients. GSEA revealed that CKAP2L can regulate the development of glioma cells through various signaling pathways, especially that of the cell cycle. Besides, knockdown of CKAP2L inhibited the proliferation and invasion of cells by reducing the expression level of cell cycle-related proteins. Our findings elucidate the possible role of CKAP2L in glioma, identifying a reliable biological biomarker for its diagnosis and treatment.

Data Availability

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Supplementary Materials

Supplementary 1. Table S1: characteristics of patients with glioma based on TCGA.

Supplementary 2. Table S2: characteristics of patients with glioma based on CGGA.

Supplementary 3. Table S3: characteristics of patients with glioma based on CGGA microarray.

Supplementary 4. Table S4: characteristics of patients with glioma based on clinical patients.

References

L. L. Morgan, “The epidemiology of glioma in adults: a "state of the science" review,” Neuro-Oncology, vol. 17, no. 4, pp. 623-624, 2015.
View at: Publisher Site | Google Scholar
T. A. Dolecek, J. M. Propp, N. E. Stroup, and C. Kruchko, “CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2005-2009,” Neuro-Oncology, vol. 14, Supplement 5, pp. v1–49, 2012.
View at: Publisher Site | Google Scholar
L. Xiong, F. Wang, and X. Qi Xie, “Advanced treatment in high-grade gliomas,” Journal of BUON: Official Journal of the Balkan Union of Oncology, vol. 24, no. 2, pp. 424–430, 2019.
View at: Google Scholar
N. A. Bush, S. M. Chang, and M. S. Berger, “Current and future strategies for treatment of glioma,” Neurosurgical Review, vol. 40, no. 1, pp. 1–14, 2017.
View at: Publisher Site | Google Scholar
L. Larsson, “Oncogene- and tumor suppressor gene-mediated suppression of cellular senescence,” Seminars in Cancer Biology, vol. 21, no. 6, pp. 367–376, 2011.
View at: Publisher Site | Google Scholar
K. Tomczak, P. Czerwińska, and M. Wiznerowicz, “The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge,” Contemporary oncology (Poznan, Poland), vol. 19, pp. A68–A77, 2015.
View at: Google Scholar
T. Barrett, S. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, Database issue, pp. D991–D995, 2013.
View at: Publisher Site | Google Scholar
Z. Wang, Z. Wang, C. Zhang et al., “Genetic and clinical characterization of B7-H3 (CD276) expression and epigenetic regulation in diffuse brain glioma,” Cancer Science, vol. 109, no. 9, pp. 2697–2705, 2018.
View at: Publisher Site | Google Scholar
E. Feng, T. Liang, X. Wang et al., “Correlation of alteration of HLA-F expression and clinical characterization in 593 brain glioma samples,” Journal of Neuroinflammation, vol. 16, no. 1, p. 33, 2019.
View at: Publisher Site | Google Scholar
X. Guan, C. Zhang, J. Zhao, G. Sun, Q. Song, and W. Jia, “CMTM6 overexpression is associated with molecular and clinical characteristics of malignancy and predicts poor prognosis in gliomas,” eBioMedicine, vol. 35, pp. 233–243, 2018.
View at: Publisher Site | Google Scholar
T. Yumoto, K. Nakadate, Y. Nakamura et al., “Radmis, a novel mitotic spindle protein that functions in cell division of neural progenitors,” PLoS One, vol. 8, no. 11, article e79895, 2013.
View at: Publisher Site | Google Scholar
G. Xiong, L. Li, X. Chen et al., “Up-regulation of CKAP2L expression promotes lung adenocarcinoma invasion and is associated with poor prognosis,” OncoTargets and Therapy, vol. 12, pp. 1171–1180, 2019.
View at: Publisher Site | Google Scholar
P. Wang and X. He, “Oncogenic and prognostic role of CKAP2L in hepatocellular carcinoma,” International Journal of Clinical and Experimental Pathology, vol. 13, no. 5, pp. 923–933, 2020.
View at: Google Scholar
Y. Fu, Q. Z. Zhou, X. L. Zhang, Z. Z. Wang, and P. Wang, “Identification of hub genes using co-expression network analysis in breast cancer as a tool to predict different stages,” Medical Science Monitor: International Medical Journal of Experimental and Clinical Research, vol. 25, pp. 8873–8890, 2019.
View at: Publisher Site | Google Scholar
N. Xu, S. Chen, T. Lin et al., “Development and validation of hub genes for lymph node metastasis in patients with prostate cancer,” Journal of Cellular and Molecular Medicine, vol. 24, no. 8, pp. 4402–4414, 2020.
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
T. Barrett, D. Troup, S. Wilhite et al., “NCBI GEO: archive for functional genomics data sets--10 years on,” Nucleic Acids Research, vol. 39, no. Database, pp. D1005–D1010, 2011.
View at: Publisher Site | Google Scholar
“The TCGA Legacy,” Cell, vol. 173, no. 2, pp. 281-282, 2018.
View at: Publisher Site | Google Scholar
Z. Liu, F. Shen, H. Wang et al., “Abnormally high expression of HOXA2 as an independent factor for poor prognosis in glioma patients,” Cell Cycle (Georgetown, Tex), vol. 19, no. 13, pp. 1632–1640, 2020.
View at: Publisher Site | Google Scholar
J. Schwarz, J. Payton, R. Rashmi et al., “Pathway-specific analysis of gene expression data identifies the PI3K/Akt pathway as a novel therapeutic target in cervical cancer,” Clinical Cancer Research: An Official Journal of the American Association for Cancer Research, vol. 18, no. 5, pp. 1464–1471, 2012.
View at: Publisher Site | Google Scholar
J. Lamb, E. Crawford, D. Peck et al., “The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease,” Science (New York, NY), vol. 313, no. 5795, pp. 1929–1935, 2006.
View at: Publisher Site | Google Scholar
Z. Liu, X. Wang, G. Yang et al., “Construction of lncRNA-associated ceRNA networks to identify prognostic lncRNA biomarkers for glioblastoma,” Journal of Cellular Biochemistry, vol. 121, no. 7, pp. 3502–3515, 2020.
View at: Publisher Site | Google Scholar
H. Tu, M. Wu, W. Huang, and L. Wang, “Screening of potential biomarkers and their predictive value in early stage non-small cell lung cancer: a bioinformatics analysis,” Translational Lung Cancer Research, vol. 8, no. 6, pp. 797–807, 2019.
View at: Publisher Site | Google Scholar
B. K. Rasmussen, S. Hansen, R. J. Laursen et al., “Epidemiology of glioma: clinical characteristics, symptoms, and predictors of glioma patients grade I-IV in the the Danish Neuro-Oncology Registry,” Journal of Neuro-Oncology, vol. 135, no. 3, pp. 571–579, 2017.
View at: Publisher Site | Google Scholar
H. Ohgaki and P. Kleihues, “The definition of primary and secondary glioblastoma,” Clinical Cancer Research, vol. 19, no. 4, pp. 764–772, 2013.
View at: Publisher Site | Google Scholar
X. Chen, Y. Yan, J. Zhou et al., “Clinical prognostic value of isocitrate dehydrogenase mutation, O-6-methylguanine-DNA methyltransferase promoter methylation, and 1p19q co-deletion in glioma patients,” Annals of Translational Medicine, vol. 7, no. 20, p. 541, 2019.
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
T. A. Lawrie, D. Gillespie, T. Dowswell et al., “Long-term neurocognitive and other side effects of radiotherapy, with or without chemotherapy, for glioma,” Cochrane Database of Systematic Reviews, vol. 8, article CD013047, 2019.
View at: Publisher Site | Google Scholar
H. Zhang, X. Zhang, X. Li et al., “Effect of CCNB1 silencing on cell cycle, senescence, and apoptosis through the p53 signaling pathway in pancreatic cancer,” Journal of Cellular Physiology, vol. 234, no. 1, pp. 619–631, 2018.
View at: Publisher Site | Google Scholar
T. Narita, S. Hatakeyama, T. Yoneyama et al., “Clinical implications of serum N-glycan profiling as a diagnostic and prognostic biomarker in germ-cell tumors,” Cancer Medicine, vol. 6, no. 4, pp. 739–748, 2017.
View at: Publisher Site | Google Scholar
J. Contessa, M. S. Bhojani, H. H. Freeze, B. D. Ross, A. Rehemtulla, and T. S. Lawrence, “Molecular imaging of N-linked glycosylation suggests glycan biosynthesis is a novel target for cancer therapy,” Clinical Cancer Research: An Official Journal of the American Association for Cancer Research, vol. 16, no. 12, pp. 3205–3214, 2010.
View at: Publisher Site | Google Scholar
H. Yu, S. Zhang, A. N. Ibrahim, Z. Deng, and M. Wang, “Serine/threonine kinase BUB1 promotes proliferation and radio-resistance in glioblastoma,” Pathology, Research and Practice, vol. 215, no. 8, p. 152508, 2019.
View at: Publisher Site | Google Scholar
S. Chen, Y. Wang, C. Ni, G. Meng, and X. Sheng, “HLF/miR-132/TTK axis regulates cell proliferation, metastasis and radiosensitivity of glioma cells,” Biomedicine & Pharmacotherapy = Biomedecine & Pharmacotherapie, vol. 83, pp. 898–904, 2016.
View at: Publisher Site | Google Scholar
B. Lee, A. K. Nalla, I. R. Stock, T. C. Shear, K. L. Black, and J. S. Yu, “Oxidative stimuli-responsive nanoprodrug of camptothecin kills glioblastoma cells,” Bioorganic & Medicinal Chemistry Letters, vol. 20, no. 17, pp. 5262–5268, 2010.
View at: Publisher Site | Google Scholar
A. Rahman, F. Thayyullathil, S. Pallichankandy, and S. Galadari, “Hydrogen peroxide/ceramide/Akt signaling axis play a critical role in the antileukemic potential of sanguinarine,” Free Radical Biology & Medicine, vol. 96, pp. 273–289, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Li Zhu 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

290

Downloads

942

Citations