- Open Access
An integrative bioinformatics investigation and experimental validation of critically involved genes in high-grade gliomas
Diagnostic Pathology volume 17, Article number: 73 (2022)
Lack of knowledge around underlying mechanisms of gliomas mandates intense research efforts to improve the disease outcomes. Identification of high-grade gliomas pathogenesis which is known for poor prognosis and low survival is of particular importance. Distinguishing the differentially expressed genes is one of the core approaches to clarify the causative factors.
Microarray datasets of the treatment-naïve gliomas were provided from the Gene Expression Omnibus considering the similar platform and batch effect removal. Interacting recovery of the top differentially expressed genes was performed on the STRING and Cytoscape platforms. Kaplan–Meier analysis was piloted using RNA sequencing data and the survival rate of glioma patients was checked considering selected genes. To validate the bioinformatics results, the gene expression was elucidated by real-time RT-qPCR in a series of low and high-grade fresh tumor samples.
We identified 323 up-regulated and 253 down-regulated genes. The top 20 network analysis indicated that PTX3, TIMP1, CHI3L1, LTF and IGFBP3 comprise a crucial role in gliomas progression. The survival was inversely linked to the levels of all selected genes. Further analysis of RNA sequencing data indicated a significant increase in all five genes in high-grade tumors. Among them, PTX3, TIMP1 and LTF did not show any change in low-grade versus controls. Real-time RT-qPCR confirmed the in-silico results and revealed significantly higher expression of selected genes in high-grade samples compared to low-grade.
Our results highlighted the role of PTX3 and TIMP1 which were previously considered in glioma tumorigenesis as well as LTF as a new potential biomarker.
Gliomas, the most prevalent and aggressive type of brain tumors, are classified by the world health organization (WHO) as grade I-IV according to clinical and histopathological characteristics. Grade I comprise benign and relatively low-risk gliomas. Grade II is known as low-grade gliomas (LGG), displaying well-differentiated tumor cells with a better prognosis. WHO grade III-IV classifications include high-grade gliomas (HGG). They are characterized by undifferentiated cells and consist of anaplastic gliomas (grade III) and glioblastoma multiforme (GBM) (grade IV) .
LGG, which generally affects young adults with an average age of 40, eventually leads to GBM and death with a survival of fewer than ten years. In recent years, a new classification has been proposed by the WHO regarding the status of the isocitrate dehydrogenase (IDH) gene. Tumors bearing mutations in IDH are sub-classified into secondary GBM, indicating that they are originated from a LGG tumor . However, about 90% of GBM are de novo cases with poor prognosis and a median survival of around one year .
Despite the significant advances in understanding gliomas genomic alterations, the molecular basis of HGG remains to be explored. The high molecular heterogeneity of gliomas could explain the complexity of its prognosis prediction .
The clinical aggressiveness within gliomas has been explored by the development of high‐throughput technologies . Glioma-related bioinformatics studies mostly conducted by using the online available microarray datasets such as Gene Expression Omnibus (GEO) , The Cancer Genome Atlas (TCGA)  and ArrayExpress . However, shortcomings in sample selection, including de-novo versus secondary GBM and treatment-naïve versus treatment-experienced have led to heterogeneous outcomes of differentially expressed genes (DEGs). Diverse data processing methods and technological platforms could be mentioned as additional confounding factors and may result in controversies. Integrated bioinformatics and co-expression networks analysis assist network-based gene screening and strengthen the statistical analysis .
In the present study, by the integrative and network-based approach, we analyzed DEGs between adult LGG and HGG. We also determined the involved signaling pathways and potential molecular interactions that are significantly associated with WHO grading and prognosis. DEGs revealing major alterations in their network interactions between LGG and HGG were selected for further validation in patient’s tumor samples.
Materials and methods
Data integration, batch effect removal
GSE4290, GSE15824, GSE19728, GSE43378, and GSE51062 microarray expression datasets were collected from the GEO database (http://www.ncbi.nih.gov/geo) using the GEO query package of R software. All of the selected microarrays had been performed on treatment-naïve adult gliomas by the GPL570 detection platform. The datasets comprised a total of 345 samples, including HGG, LGG, glioma cell lines, and normal brain tissues. The samples of the normal brain, cell lines, unknown subtypes and secondary or recurrent tumors were excluded from the study and ultimately 301 treatment-naïve adult samples including 238 HGG and 63 LGG were recruited for the subsequent analysis. The average expression value of all mapped probe sets was considered for each gene. The preprocessed expression datasets were merged by the R program into one global expression matrices according to the gene symbol. Using the combat method in R, the technical heterogeneity or batch effect was removed across the datasets.
The limma package of R was utilized to determine the DEGs between LGG and HGG in a linear model. To decrease the false positive rate and false discovery rate the cut-off criteria of |log2FC|> 1 and an adjusted P-value < 0.001 were selected. The results were presented in a volcano plot of log-fold changes on the X-axis versus adjusted statistical significance on the Y-axis.
Functional annotation and pathway analyses
Gene Ontology (GO) enrichment analysis and pathway enrichment analysis based on and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were conducted by Database for Annotation, Visualization, and Integrated Discovery (DAVID) online tool (https://david.ncifcrf.gov/). The p-value ˂0.05 was considered statistically significant.
Construction of protein–protein interaction network
The STRING online tool (https://string-db.org)  was constructed to generate the protein–protein interaction (PPI) network of the top 20 up and down-regulated DEGs. To recognize the target genes, the PPI network was visualized by the network analyzer in the Cytoscape software (version 3.8.2, http://www.cytoscape.org/) .
Top DEGs validation by glioma RNA sequencing data
DEGs with significant protein–protein interactions were selected for further assessments through the Gene Expression Profiling and Interactive Analyses (GEPIA) online RNA sequencing-based database (http://gepia.cancer-pku.cn) . This database represents the gene expression variations in cancer samples based on data from TCGA and The Genotype-Tissue Expression (GTEx). Kaplan–Meier and Cox proportional hazard analysis were conducted to check the survival rate of glioma patients by subgrouping into the low and high expression of selected genes.
By the comparison between GBM, LGG and healthy tissue, the expression levels of selected genes were presented as mean (standard deviation) in boxplots. A p-value ˂0.01 was considered statistically significant.
Patient’s samples collection
To validate the in-vitro expression of selected hub genes, we obtained a series of fresh LGG and HGG samples. The research was approved by the Ethics Committees of Tarbiat Modares University and performed under the Helsinki declaration of 1975, as revised in 2013. Following informed consent, 11 LGG and 11 HGG tissue samples were collected from Shariati hospital, Tehran, Iran. A part of the tumor samples was sent for pathological analysis.
Gene expression analysis
Total mRNA of the samples was extracted by the All-In-One kit (Bio Basic, Canada). The reverse transcription‑quantitative polymerase chain reaction (RT-qPCR) was performed on extracted mRNA as described before . For normalizing target genes, we selected eukaryotic translation initiation factor 2B subunit alpha (EIF2B1) as the reference gene . By evaluation on GEPIA, EIF2B1 showed higher expression and lower fluctuation compared to other recommended reference genes such as C-terminal binding protein 1 (CTBP1), mitochondrial ribosomal protein L9 (MRPL19) and TATA box-binding protein (TBP).
The primer sequences of selected and reference genes were as follow: forward PTX3 5'-TGCATCTCCTTGCGATTCTGT-3', reverse PTX3 5'-AGCTTGTCCCATTCCGAGTG-3', forward CHI3L1 5'-ATGATGTGACGCTCTACGGC-3', reverse CHI3L1 5'-ACTCTGGGTGTTGGAGGCTA-3', forward IGFBP3 5'-CAGAATATGGTCCCTGCCGTAG-3', reverse IGFBP3 5'-TTTGGAAGGGCGACACTGC-3', forward TIMP1 5'-GCTTCTGGCATCCTGTTGTTG-3', reverse TIMP1 5'-GGTGGTCTGGTTGACTTCTGG-3', forward LTF 5'-ATCGCCCTGGTGCTGAAAG-3', reverse LTF 5'- GGGTCACTGCTTTGTTGGGA-3', forward IEF2B1 5'-CGGACGTTGCTGGAGTTCTT-3', reverse EIF2B1 5'-AGGCAAGACTGATGAAGCGG-3'.
To perform statistical analysis, the Mann–Whitney test was applied for the comparison of the groups. GraphPad Prism8 was used to generate plots signifying the levels of gene expression. A difference of p-value ˂0.05 was considered statistically significant.
Identification of DEGs in merged datasets
The principal component analysis (PCA) as an unsupervised clustering approach was constructed to classify LGG and HGG samples before and after batch effect corrections. The proper discrepancy between LGG and HGG samples was detected before adjusting the batch effect. Admissible clustering was reached following batch effect removal (Fig. 1). According to the PCA results, five samples were considered an outlier and were removed from the subsequent analysis. By merging the datasets, a total of 22,189 common genes were considered for further analysis. A total of 576 DEGs were identified comprising 323 upregulated and 253 down-regulated genes. The results were illustrated as a volcano plot in Fig. 2, in which the up and down-regulated genes are represented in green and red dots, respectively.
Functional annotation and pathway analyses
To discover the biological functions of the identified DEGs, the up-regulated and down-regulated genes were separately enriched in three categories regarding the molecular function, biological process, and cellular component. The up-regulated DEGs in the molecular function class were mainly enriched in extracellular matrix (ECM) structural constituent, heparin, integrin and PDGF binding. These DEGs were also significantly enriched in the ECM organization, angiogenesis and cell division in the biological process. In the cellular component category, they were enriched in the ECM, extracellular space, and extracellular region. According to the KEGG pathway analysis, the up-regulated DEGs were frequently enriched in ECM-receptor interaction and focal adhesion (Fig. 3).
Construction of PPI network
The top 20 most significant up and down-regulated DEGs are shown in Table 1. The identified top 20 DEGs were mapped into the STRING database to find the interactions between their corresponding proteins. The construction of DEGs networks indicated significant interactions between pentraxin 3 (PTX3), chitinase-3-like protein 1 (CHI3L1), insulin-like growth factor-binding protein 3 (IGFBP3), tissue inhibitor of matrix metalloproteinases 1 (TIMP1) and lactotransferrin (LTF). As illustrated by Cytoscape software in Fig. 4, the 25 genes of the top 20 DEGs interact to each other, amongst them TIMP1, CHI3L1 and PTX3 have the most interface.
In-silico validation of candidate genes
The expression levels of candidate genes were further assessed in 163 HGG, 518 LGG and 207 healthy patients’ samples on the GEPIA web server. Initially, the correlation between gene expression and overall survival of glioma were assessed by Kaplan–Meier analysis. The patient’s survival in high mRNA expression levels of all five selected genes was less than that in the low expression group (Fig. 5).
By dividing glioma patients into GBM and LGG groups, the expression level of candidate genes in each of these groups was compared with normal samples. The data indicated a significant increase in all five genes in GBM compared to LGG and control tissues. Gene expression levels of IGFBP3 and CHI3L1 significantly differed between LGG and normal brain tissue. However, the expression of PTX3, TIMP1 and LTF were not altered in LGG versus controls indicating their specificity in GBM pathogenesis (Fig. 6).
The treatment-naïve glioma cases included eight male and fourteen female patients, with an average age of 52 ± 15 years old. The LGG samples consisted of five cases of diffuse oligodendroglioma (grade II) and six cases of diffuse astrocytoma (grade II), and the HGG samples included five cases of anaplastic astrocytoma (grade III) and six cases of glioblastoma multiform (grade IV).
In vitro validation of gene expression by RT-qPCR
We conducted an RT-qPCR and confirmed the in-silico expression data of PTX3, IGFBP3, TIMP1, CHI3L1 and LTF. Normalized results of the expression, relative to the expression level of EIF2B1 mRNA indicated significantly higher expression of selected genes in HGG compared to LGG (p-value < 0.05) (Fig. 7). The highest fold changes of 2.33 and 2.23 were calculated for TIMP1 and LTF, respectively. CHI3L1 and IGFBP3 exhibited almost a twofold increase, while PTX3 disclosed fold change of 1.67 in HGG compared to LGG.
We could not detect significant expression differences related to patients’ age, sex, tumor size and/or tumor location.
HGG is considered a challenging subtype of brain tumors with a poor prognosis, aggressive features and lack of targeted therapies. Although many efforts have been made to develop effective treatments, tumor heterogeneity and unknown contexts remain limitations of discovered therapies . Identification of the key genes and pathways involved in HGG pathogenesis is of particular importance to develop new therapeutic approaches. In recent decades, high throughput technologies produce huge data that are helpful for a well understanding of the molecular basis of tumorigenesis.
Here, we analyzed multiple expression array datasets to compare gene expression profiles of HGG versus LGG to provide the critically involved genes. Our results elucidated that the PTX3 gene acts as a master regulator in the HGG gene regulatory network by connectivity to four other DEGs. Compared to LGG, HGG showed increased PTX3 mRNA levels in patients’ tumor samples. PTX3 is mainly involved in immune responses and inflammation; however also plays various roles in multiple molecular mechanisms such as cell cycle, angiogenesis, metastasis, and genomic instability . Dysregulation of PTX3 is associated with the early event in oncogenesis and may be acting as an initiating factor. By in vitro inactivation of PTX3, GBM cells migration and invasion were significantly eliminated .
Among the genes that showed a strong link with PTX, the TMP1 was already known to be associated with deteriorating and poor prognosis of glioma . Aaberg-Jessen et al. investigated the prognostic potential of TIMP1 in combination with its cell surface binding protein, CD63, correlated to tumor grade and overall survival of glioma patients. However, no additional prognostic significance of TIMP1 was observed by including the CD63. They also analyzed DEGs comparing the highest and lowest CD63 mRNA level tumors using the TCGA dataset. They reported TIMP1, PTX3 and CHI3L1 among the upregulated DEGs involved in the regulation of cell survival and cellular movement .
In our study, we observed that mRNA expression of CHI3L1, also known as YKL-40, was significantly higher in HGG either by in silico or in vitro analysis. It is expressed as a mesenchymal marker in the most aggressive and challenging subtype of GBM associated with particularly poor outcome . CHI3L1 is a secreted glycoprotein and is suggested as a promising serum marker for GBM diagnosis . However, it has more diagnostic value in recurrent GBM than newly diagnosed patients .
Among up-regulated genes found in our study, LTF transcript was marked as the first ranked (Table 1); however, its promoting role in tumorigenesis is not evident. As an iron-binding protein, LTF has a wide spectrum of protective activities such as anti-microbial and immunity. Studies in nasopharyngeal carcinoma and prostate cancer suggested a tumor suppressor function for LTF [23, 24].
Using GEPIA, we found that the expression of PTX3, TIMP1 and LTF were specificity high in GBM, whereas they were not altered in LGG compared to normal controls. One of the suggested mechanisms that can be shared between their functions is the inflammatory process. Correlation to inflammatory in glioma tumorigenic processes has been demonstrated by different studies. Increased expression of PTX3 under the influence of transcriptional factor CCAAT / enhancer-binding protein delta (CEBPD) has been shown to inhibit the phagocytosis of dying neurons . CEBPD also regulates the stemness of glioma cells by activating platelet-derived growth factor subunit A (PDGFA) expression due to inflammatory stimulation. As mentioned before, inflammation is one of the important processes in the development of GBM and response to treatment [26, 27]. It should be noted that LTF was previously reported as putative interactors of CEBPD .
We also found IGFBP3 among the top 20 DEGs that showed strong involvement with PTX3. IGFBP3 plays a direct role in signaling and cell growth pathways in association with insulin-like growth factor (IGF) 1 and 2. Although studies have suggested a tumor-suppressive role for IGFBP-3, ample evidence in various cancers indicates that its high levels were associated with disease aggressiveness . In a recent study on hepatocellular carcinoma, IGFBP3 was found highly expressed associated with disease poor prognosis . By analyzing TCGA and GTEx data, we observed a significant IGFBP3 increase in HGG compared to normal brain tissues. However, IGFBP3 gene expression was also higher in LGG than in normal controls. It may be considered that IGFBP3 stimulates the tumorigenesis in both LGG and HGG probably via different pathways.
As a conclusion, we identified robust differentially expressed genes associated with HGG and highlighted some important features of HGG tumors such as inflammation and stemness. These findings could enhance our knowledge about the molecular mechanisms underlying HGG tumors and provided reliable biomarkers for novel prognostic and therapeutic targets.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
World health organization
Gene Expression Omnibus
The Cancer Genome Atlas
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Database for Annotation, Visualization, and Integrated Discovery
Search Tool for the Retrieval of Interacting Genes/Proteins
Gene Expression Profiling and Interactive Analyses
Reverse transcription‑quantitative polymerase chain reaction
Eukaryotic translation initiation factor 2B subunit alpha
C-terminal binding protein 1
Mitochondrial ribosomal protein L9
TATA box-binding protein
Principal component analysis
Chitinase-3-like protein 1
Insulin-like growth factor-binding protein 3
Tissue inhibitor of matrix metalloproteinases 1
CCAAT / enhancer-binding protein delta
Platelet-derived growth factor subunit A
Insulin-like growth factor
Transcripts per million
Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, Ohgaki H, Wiestler OD, Kleihues P, Ellison DW. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. 2016;131:803–20.
Claus EB, Walsh KM, Wiencke JK, Molinaro AM, Wiemels JL, Schildkraut JM, Bondy ML, Berger M, Jenkins R, Wrensch M. Survival and low-grade glioma: the emergence of genetic information. Neurosurg Focus. 2015;38:E6.
Ohgaki H, Kleihues P. The definition of primary and secondary glioblastoma. Clin Cancer Res. 2013;19:764–72.
Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, Zheng S, Chakravarty D, Sanborn JZ, Berman SH, et al. The somatic genomic landscape of glioblastoma. Cell. 2013;155:462–77.
Alvarez MJ, Subramaniam PS, Tang LH, Grunn A, Aburi M, Rieckhof G, Komissarova EV, Hagan EA, Bodei L, Clemons PA, et al. A precision oncology approach to the pharmacological targeting of mechanistic dependencies in neuroendocrine tumors. Nat Genet. 2018;50:979–89.
Clough E, Barrett T. The gene expression omnibus database. Methods Mol Biol. 2016;1418:93–110.
Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol (Pozn). 2015;19:A68-77.
Brazma A, Parkinson H, Sarkans U, Shojatalab M, Vilo J, Abeygunawardena N, Holloway E, Kapushesky M, Kemmeren P, Lara GG, et al. ArrayExpress–a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2003;31:68–71.
Wang H, Diaz AK, Shaw TI, Li Y, Niu M, Cho JH, Paugh BS, Zhang Y, Sifford J, Bai B, et al. Deep multiomics profiling of brain tumors identifies signaling networks downstream of cancer driver genes. Nat Commun. 2019;10:3718.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, Jensen LJ. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41:D808-815.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
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-w102.
Ahmadi-Beni R, Vand-Rajabpour F, Ahmadifard M, Daneshpazhooh M, Noormohammadpour P, Rahmati J, Hesari KK, Yaseri M, Tabrizi M. Decreased Sox2 messenger RNA expression in basal cell carcinoma. Indian J Dermatol. 2020;65:178–82.
Bertoni APS, Iser IC, de Campos RP, Wink MR. Normalization in human glioma tissue. Methods Mol Biol. 2020;2065:175–90.
Ahmadi-Beni R, Najafi A, Savar SM, Mohebbi N, Khoshnevisan A. Role of dimethyl fumarate in the treatment of glioblastoma multiforme: a review article. Iran J Neurol. 2019;18:127–33.
Agrawal A, Singh PP, Bottazzi B, Garlanda C, Mantovani A. Pattern recognition by pentraxins. Adv Exp Med Biol. 2009;653:98–116.
Tung JN, Ko CP, Yang SF, Cheng CW, Chen PN, Chang CY, Lin CL, Yang TF, Hsieh YH, Chen KC. Inhibition of pentraxin 3 in glioma cells impairs proliferation and invasion in vitro and in vivo. J Neurooncol. 2016;129:201–9.
Groft LL, Muzik H, Rewcastle NB, Johnston RN, Knäuper V, Lafleur MA, Forsyth PA, Edwards DR. Differential expression and localization of TIMP-1 and TIMP-4 in human gliomas. Br J Cancer. 2001;85:55–63.
Aaberg-Jessen C, Sørensen MD, Matos A, Moreira JM, Brünner N, Knudsen A, Kristensen BW. Co-expression of TIMP-1 and its cell surface binding partner CD63 in glioblastomas. BMC Cancer. 2018;18:270.
Gao YF, Zhu T, Mao CX, Liu ZX, Wang ZB, Mao XY, Li L, Yin JY, Zhou HH, Liu ZQ. PPIC, EMP3 and CHI3L1 are novel prognostic markers for high grade glioma. Int J Mol Sci. 2016;17:1808.
Iwamoto FM, Hormigo A. Unveiling YKL-40, from serum marker to target therapy in glioblastoma. Front Oncol. 2014;4:90.
Holst CB, Christensen IJ, Skjøth-Rasmussen J, Hamerlik P, Poulsen HS, Johansen JS. Systemic immune modulation in gliomas: prognostic value of plasma IL-6, YKL-40, and genetic variation in YKL-40. Front Oncol. 2020;10:478.
Shaheduzzaman S, Vishwanath A, Furusato B, Cullen J, Chen Y, Bañez L, Nau M, Ravindranath L, Kim KH, Mohammed A, et al. Silencing of Lactotransferrin expression by methylation in prostate cancer progression. Cancer Biol Ther. 2007;6:1088–95.
Deng M, Zhang W, Tang H, Ye Q, Liao Q, Zhou Y, Wu M, Xiong W, Zheng Y, Guo X, et al. Lactotransferrin acts as a tumor suppressor in nasopharyngeal carcinoma by repressing AKT through multiple mechanisms. Oncogene. 2013;32:4273–83.
Ko CY, Chang LH, Lee YC, Sterneck E, Cheng CP, Chen SH, Huang AM, Tseng JT, Wang JM. CCAAT/enhancer binding protein delta (CEBPD) elevating PTX3 expression inhibits macrophage-mediated phagocytosis of dying neuron cells. Neurobiol Aging. 2012;33(422):e411-425.
Wang S-M, Lin H-Y, Chen Y-L, Hsu T-I, Chuang J-Y, Kao T-J, Ko C-Y. CCAAT/enhancer-binding protein delta regulates the stemness of glioma stem-like cells through activating PDGFA expression upon inflammatory stimulation. J Neuroinflammation. 2019;16:146.
Ahmadi-Beni R, Khoshnevisan A. An overview of crucial genes involved in stemness of glioblastoma multiforme. Neurochem J. 2017;11:259–65.
Chih DY, Park DJ, Gross M, Idos G, Vuong PT, Hirama T, Chumakov AM, Said J, Koeffler HP. Protein partners of C/EBPepsilon. Exp Hematol. 2004;32:1173–81.
Yu H, Levesque MA, Khosravi MJ, Papanastasiou-Diamandi A, Clark GM, Diamandis EP. Insulin-like growth factor-binding protein-3 and breast cancer survival. Int J Cancer. 1998;79:624–8.
Bai Q, Liu H, Guo H, Lin H, Song X, Jin Y, Liu Y, Guo H, Liang S, Song R, et al. Identification of hub genes associated with development and microenvironment of hepatocellular carcinoma by weighted gene co-expression network analysis and differential gene expression analysis. Front Genet. 2020;11:615308.
Ethics approval and consent to participate
All the patients were filled informed consent forms which were approved by the ethic committee of Tarbiat Modares University (IR.MODARES.REC.1399.135).
Consent for publication
There are no competing interests declared by the authors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ahmadi-Beni, R., Shahbazi, S. & Khoshnevisan, A. An integrative bioinformatics investigation and experimental validation of critically involved genes in high-grade gliomas. Diagn Pathol 17, 73 (2022). https://doi.org/10.1186/s13000-022-01253-0
- Adult glioma
- Integrated bioinformatics
- GEO data
- Differentially expressed genes