Next Article in Journal
Influence of Mutations and N-Glycosylation Sites in the Receptor-Binding Domain (RBD) and the Membrane Protein of SARS-CoV-2 Variants of Concern on Antibody Binding in ELISA
Previous Article in Journal
Ovotestis Isolation and Cryopreservation of Nesiohelix samarangae (Oriental Snail) as a Snail Model for Conserving Other Endangered Snail Species
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades

by
Marco A. Orda
1,2,
Peter Matthew Paul T. Fowler
1,3 and
Lemmuel L. Tayo
1,3,*
1
School of Chemical, Biological, and Materials Engineering and Sciences, Mapúa University, Manila City 1002, Philippines
2
School of Graduate Studies, Mapúa University, Manila City 1002, Philippines
3
Department of Biology, School of Health Sciences, Mapúa University, Makati City 1203, Philippines
*
Author to whom correspondence should be addressed.
Biology 2024, 13(4), 206; https://doi.org/10.3390/biology13040206
Submission received: 12 February 2024 / Revised: 7 March 2024 / Accepted: 21 March 2024 / Published: 23 March 2024
(This article belongs to the Section Biochemistry and Molecular Biology)

Abstract

:

Simple Summary

Currently, gliomas still stand as one of the global health concerns due to its resistance to treatment strategies, resulting in poor survival rate and high mortality despite sustained efforts. Even though chemotherapeutic and targeted drugs for glioma exist, its heterogeneity per grade and cell origin poses as a significant contributing factor to its negative response to treatment. With this, our study utilized a systems biology approach to investigate correlations between various grades of glioma using DNA microarray profiles of samples that contain pilocytic astrocytoma, oligodendroglioma, anaplastic astrocytoma, and glioblastoma multiforme gene expression data. The highly preserved gene clusters (modules) and their corresponding hub genes were identified, which provided further insights into the signaling pathways and cellular processes involved. Drug repurposing was also performed based on the upregulated and downregulated hub genes. This approach determined potential drug candidates that are known to have mediating effects on the signaling pathways, thereby highlighting the possibility of the gene network as a potential therapeutic avenue.

Abstract

Gliomas have displayed significant challenges in oncology due to their high degree of invasiveness, recurrence, and resistance to treatment strategies. In this work, the key hub genes mainly associated with different grades of glioma, which were represented by pilocytic astrocytoma (PA), oligodendroglioma (OG), anaplastic astrocytoma (AA), and glioblastoma multiforme (GBM), were identified through weighted gene co-expression network analysis (WGCNA) of microarray datasets retrieved from the Gene Expression Omnibus (GEO) database. Through this, four highly correlated modules were observed to be present across the PA (GSE50161), OG (GSE4290), AA (GSE43378), and GBM (GSE36245) datasets. The functional annotation and pathway enrichment analysis done through the Database for Annotation, Visualization, and Integrated Discovery (DAVID) showed that the modules and hub genes identified were mainly involved in signal transduction, transcription regulation, and protein binding, which collectively deregulate several signaling pathways, mainly PI3K/Akt and metabolic pathways. The involvement of several hub genes primarily linked to other signaling pathways, including the cAMP, MAPK/ERK, Wnt/β-catenin, and calcium signaling pathways, indicates potential interconnectivity and influence on the PI3K/Akt pathway and, subsequently, glioma severity. The Drug Repurposing Encyclopedia (DRE) was used to screen for potential drugs based on the up- and downregulated hub genes, wherein the synthetic progestin hormones norgestimate and ethisterone were the top drug candidates. This shows the potential neuroprotective effect of progesterone against glioma due to its influence on EGFR expression and other signaling pathways. Aside from these, several experimental and approved drug candidates were also identified, which include an adrenergic receptor antagonist, a PPAR-γ receptor agonist, a CDK inhibitor, a sodium channel blocker, a bradykinin receptor antagonist, and a dopamine receptor agonist, which further highlights the gene network as a potential therapeutic avenue for glioma.

Graphical Abstract

1. Introduction

Glioma is a type of brain tumor that originates from glial cells. A common type of glioma, astrocytoma, emerges from astrocytes, cells that primarily provide structural support to neurons and regulate the blood–brain barrier. Astrocytomas are often considered to be more aggressive than other types of glioma due to their high cell proliferation, invasiveness, and recurrence [1]. Astrocytomas can be further classified based on WHO grades of increasing malignancy: Grade I: pilocytic astrocytoma (PA); Grade II: diffuse astrocytoma (DA); Grade III: anaplastic astrocytoma (AA); and Grade IV: glioblastoma multiforme (GBM) [2]. On the other hand, oligodendrogliomas are less-common gliomas that originate from oligodendrocytes, which produce and maintain the myelin sheath around the neuronal axons for nerve signal transmission. Like astrocytomas, they are classified based on WHO grades, wherein oligodendroglioma (OG) is Grade II and anaplastic oligodendroglioma (AO) is Grade III [3]. They are noted for their specific genetic alterations that affect their prognosis and treatment, specifically 1p/19q codeletion, as these alterations are usually present in both grades; in contrast, different grades of astrocytomas contain key sets of genetic alterations [2,4]. Despite their differences, astrocytomas and oligodendrogliomas may share similar genetic alterations and clinical symptoms [5].
Along with surgery and radiation therapy, chemotherapeutic and targeted drugs, particularly temozolomide and bevacizumab, are used as treatment for glioma [6]. Several studies have also indicated the use of inhibitors for the epidermal growth factor receptor (EGFR) due to the strong involvement of the downstream PI3K/Akt pathway in glioma cell proliferation [7,8]. However, resistance against temozolomide and the toxicity of bevacizumab have been documented [9,10], while EGFR inhibitors tend to struggle in penetrating the blood–brain barrier [11]. Hence, the poor prognosis of the five-year survival rate of patients with GBM has not improved significantly despite sustained efforts [12,13]. As a result, systems biology approaches, which excel in identifying gene network hubs that play central roles in regulating overall cellular processes, are currently being considered. One of these approaches is the weighted gene co-expression network analysis (WGCNA), which is used to characterize the patterns of gene association found in the gene expression data of several diseases. This method enables identifying potential therapeutic targets linked to particular biological processes or illnesses by categorizing genes into modules according to their expression profiles, offering a new point of view for the comprehensive understanding of the molecular foundations of related diseases. Furthermore, the identification of hub genes in the modules can also be done to provide insight into the main pathways and processes that are mutated. Lastly, these hub genes can then be subjected to drug repurposing techniques, which has emerged as a novel approach in the field of drug discovery due to its potential to save time and costs, as it enables the screening of approved and experimental drugs that can interfere with the mutated genes.
In this study, gene expression datasets for four different grades of glioma, represented by PA (GSE50161), OG (GSE4290), AA (GSE43378), and GBM (GSE36245), were acquired from the Gene Expression Omnibus (GEO) and were then used to perform a meta-analysis approach through WGCNA to determine highly conserved modules among all the datasets. Four such modules were identified and then subjected to functional annotation and pathway enrichment analysis, where it was found that they were highly involved in signal transduction, transcription regulation, protein binding, and metabolic changes. In each module, the protein–protein association network and key hub genes were determined, which further allowed for screening of potential already-existing drugs. With this bioinformatics approach, this study aims to offer further insights for understanding the pivotal points at which different glioma grades align genetically to pave the way for the enhancement of existing treatment strategies or the development of novel ones.

2. Materials and Methods

2.1. Dataset Acquisition and Preparation

The microarray gene expression datasets were sourced from the National Center for Biotechnology Information: Gene Expression Omnibus (NCBI GEO) [https://www.ncbi.nlm.nih.gov/geo/ accessed on 7 December 2023] online database for WGCNA analysis. The datasets gathered contained expression data taken from primary tumor samples of patients diagnosed with pilocytic astrocytoma (PA), oligodendroglioma (OG), anaplastic astrocytoma (AA), and glioblastoma multiforme (GBM). Furthermore, only the datasets carried out through GPL570-HG-U133 Plus 2 Affymetrix Human Genome U133 Plus 2.0 Array were selected in order to prevent inconsistencies arising from probe design, normalization issues, and batch effects [14,15,16]. This initial screening resulted in the exclusion of datasets for diffuse astrocytoma (DA) and anaplastic oligodendroglioma (AO) from further investigation. A total of 113 samples were obtained. Table 1 summarizes the details and distribution of samples for the four tumor types.
All screened raw data underwent background correction, quantile normalization, and log-2 transformation through the robust multi-array average (RMA) method in the affy package in Bioconductor using R [http://www.bioconductor.org accessed on 19 December 2023]. A boxplot of expression values and sample clustering dendrograms were generated to visually inspect the resulting data for outliers. After pre-processing, control probes were removed to eliminate non-biological variation and focus on the samples containing primary tumor expression. Then, the expression data were filtered to select the genes whose mean and variance throughout all the samples in each dataset were higher than the 20% percentile cut-off. Probe IDs were converted into gene symbols through the AnnotationDbi function using the hgu133plus2.db database for biological interpretation and analysis. Lastly, only the probes that were present in all datasets were used, and samples with no values after log-2 transformation were filtered out using the WGCNA R package’s goodSamplesGenesMs function.

2.2. Weighted Gene Co-Expression Network Analysis (WGCNA)

2.2.1. Scale-Free Network Approximation

Using the pickSoftThreshold function of the WGCNA R package, a plot of the scale-free topology fit versus power index (1–20) was generated. The soft-thresholding power (β) was then selected as the lowest power for which the scale-free topology criterion is satisfied. In a scale-free topology, the distribution of gene degrees (number of gene connections) follows the power-law distribution. This is often evaluated based on how well the co-expression network conforms to the linear relationship between log of connectivity and log of connectivity probability—by plotting the scale-free topology fit versus soft-thresholding power and selecting the power at which the scale-free topology fit reaches a plateau or a reasonably high value that indicates a good fit [17]. Hence, an approximate straight-line relationship was plotted using the values for soft connectivity (k) to evaluate the chosen soft-thresholding power.

2.2.2. TOM-Based Network Construction and Module Identification

Using Pearson’s correlation with the network type “signed,” the adjacency matrices were generated for network construction. The Topological Overlap Measure (TOM) dissimilarities were then calculated using the soft-thresholding power to raise the adjacency matrices. This allowed the emphasis of the strong correlations and the downweighting of weak correlations in performing fast hierarchal clustering of genes using the flashClust function. To produce gene dendrograms using the highly correlated genes, the hclust function was used to perform hierarchal clustering based on the distance matrix of the expression profiles with the “average” method. Gene modules were identified based on hierarchal clustering through the cutreeHybrid function, in which the deep split parameter was tested from 0 to 3 for branch splitting sensitivity [18]. Stable modules that remain consistent across different values of a deep split parameter indicate that the clustering algorithm has accurately identified significant clusters within the datasets, while varying results would suggest additional validation of the parameters used, such as the soft-thresholding power.

2.2.3. Module Preservation Analysis

Module preservation analysis measures the extent to which the connectivity patterns of the reference network are preserved in other networks to evaluate the biological relevance of the gene modules across different diseases or disease states. The modulepreservation function from the WGCNA R package was used to analyze gene co-expression network preservation of PA, OG, AA, and GBM modules qualitatively and quantitatively with the network type “signed,” the number of permutations set to 100, and a minimum module size of 30 genes. After then, the eigengene-based connectivity (kME), a measure of the correlation between the expression profile of each gene and its gene co-expression network, was calculated to identify the genes within each module using the moduleEigengenes function. Module membership was associated to the correlation between the gene expression profile and its module eigengene.

2.3. Functional Annotation and Pathway Enrichment

The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used to perform functional annotation clustering on the modules of interest [19]. Within DAVID Gene Ontology (GO), the databases for biological processes (BP), cellular components (CC), and molecular functions (MF) were utilized. This database contains information on the basic features and activities of the provided genes based on the proteins they encode. The classification stringency was set to medium, and only significant GO terms (p adj. < 0.05) that have enrichment scores greater than 1.3 were chosen for analysis. For pathway enrichment analysis, searches through the Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed. KEGG terms that scored significantly and clustered with the selected GO terms were then used to provide details on the overall functions of the modules of interest [20].

2.4. Identification of Protein-Protein Interaction (PPI) and Hub Genes

In generating PPI networks, genes within each of the highly preserved modules were considered for potential associated protein–protein interactions through the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database [21]. PPI networks with a minimum interaction score of 0.7 (high confidence) were built for each module of interest. The generated networks were then imported into Cytoscape to find hub genes based on degree centrality using the Cytohubba feature, which measures the number of interactions that the gene has within the network that correlates with its essentiality [22,23]. In this case, the top 10 genes with the highest interactions within the modules were considered hub genes.

2.5. Signature-Based Drug Repurposing

The top 10 hub genes identified based on their degree centrality were first classified as either being “upregulated” or “downregulated” using GEO2R [https://www.ncbi.nlm.nih.gov/geo/geo2r/ accessed on 23 December 2023] prior to searching for potential drugs. The Drug Repurposing Encyclopedia (DRE) webserver was used to search for available drugs that can be repurposed based on their transcription profiles from the Molecular Signatures Database (MSigDB) and Connectivity Map (CMap) [24,25]. The groups of upregulated and downregulated genes were submitted for drug repurposing tests based on molecular signature screening. Experimental drugs without an indicated mechanism of action were excluded from the analysis, and only candidates with false discovery rates (FDR) of less than 0.05 were considered.

3. Results

3.1. Weighted Gene Co-Expression Network Analysis (WGCNA)

3.1.1. Data Pre-Processing and Approximation of Scale-Free Networks

After the data preparation and filtering stages (Figure A1), a total of 25,466 genes remained. All samples for each dataset were included in the weighted gene co-expression network analysis, as no significant outliers were identified based on the sample clustering dendrograms of the datasets in Figure A2. Figure 1 shows the variation of the scale-free topology fit index versus the soft threshold (β) values ranging from 1 to 20. At β > 10, it can be observed that the index across all the datasets has stabilized and has steadily plateaued onwards. This indicates that the scale-free topology fit for β > 10 is not any more significantly affected by the increase in power and that a robust scale-free structure has been achieved. With this, it was decided that the scale-free topologies for all datasets would be approximated as β = 10. The approximation of the optimal soft-thresholding power for the scale-free topology fit is crucial, as it guarantees robustness, efficiency, and resilience to errors of the gene network. Lower soft-thresholding power leads to denser gene networks with more connections, supporting the selection of the least power at which the network indices have stabilized. This ensures that the resulting modules will contain biologically relevant relationships between genes.
Figure 2 shows the corresponding log-log plot for the GBM dataset, constructed by plotting the connectivity frequencies corresponding to each value of connectivity at a soft-thresholding power of 10. This dataset exhibited the most accurate representation of a scale-free network by achieving the highest R2 value of 0.9 among the datasets, which suggests that gene co-expression networks derived from this dataset are more likely to have related expression patterns, features, or information for confident interpretation of the network regarding its biological significance. Hence, the GBM dataset was used as the reference dataset for WGCNA.

3.1.2. TOM-Based Network Construction and Module Identification

Several approaches to conducting meta-analysis through WGCNA exist, one of which involves designating a single dataset as a reference and then projecting the module eigengenes of other datasets onto this reference. The choice of the reference dataset plays a crucial role in determining the network’s robustness, considering factors such as the dataset’s sample size, scale-free network approximation, and the resolution of clustering in the gene dendrograms based on the TOM gene dendrograms [22]. Following these criteria, the GBM dataset was chosen as the reference for the meta-analysis, and the identified modules are shown in Figure 3. Aside from this, it aligns with the fact that GBM is the end stage of glioma, which makes it suitable to be referenced to glioma datasets of lower grades, as it may indicate modules that are involved in glioma progression from the initial and latter stages. It can be observed from the sensitivity plots that the same sets of modules remained despite the increasing sensitivity, which may be due to strong correlation of the genes, having fewer significant gene clusters, or both (Figure S1). All in all, WGCNA identified a total of 21 gene co-expression modules classified arbitrarily by different colors: green-yellow (270), gold (100), red (492), midnight blue (140), yellow (1215), purple (327), tan (265), cyan (180), light green (80), light cyan (119), black (409), magenta (351), pink (323), gray60 (108), light yellow (32), green (1135), brown (1775), salmon (256), turquoise (2500), gray (2500), and blue (2500). These modules serve as clusters of related genes that exhibit coordinated expression patterns across the glioma grades, which can provide further insights into the underlying molecular foundations and regulatory networks in glioma.

3.2. Module Preservation Analysis

The preservation of the modules identified within the gene co-expression network of GBM was assessed in the datasets for PA, OG, and AA using the Z-summary score, which measures the connectivity of the modules that indicate the gain of preservation; a higher Z score indicated better preservation. As shown in Figure 4, there are few significant variations between the preserved modules across the datasets. In addition, most modules retain their relative position in the graphs but with more modules surpassing the threshold per dataset, as demonstrated by the increase in the Z-score values of the modules from the PA to the AA datasets. Biologically, this incidence is consistent with the severity of glioma, since the module eigengenes projected on the GBM dataset reflected increasing preservation from PA to AA, demonstrating the similarity of genes towards GBM. It may also be speculated that this relates to glioma progression. This results in the AA dataset having the highest number of significantly preserved modules, although only the modules with preservation scores higher than Z = 10 across all datasets—represented by the green, yellow, brown, and gray modules—were identified as modules of interest since the analysis specifically focused only on shared networks across the different glioma grades.
In addition, the analysis of module membership included the application of kME (eigengene-based connectivity) to measure the connectivity of each gene within a specific module. This was determined by calculating the correlation between the gene’s expression profile and the module eigengene, which functions as a representative expression profile for that module [26]. After that, the genes in each module were ranked according to their kME values in each dataset (Figure S2). The genes with the highest maximum rank across all datasets were chosen for the subsequent functional annotation and pathway enrichment analysis. Such genes consistently showed strong correlation within the particular modules across different glioma grades, indicating their critical roles in the biological processes that these modules represent.

3.3. Functional Annotation and Pathway Enrichment

The DAVID webserver was used to cluster functional annotations based on the top genes in each module. For each module, the top-enriched Gene Ontology (GO) terms are shown in Figure 5a–c. Figure 5d shows KEGG pathways that were prioritized only if they were enriched in the same cluster as the top GO keywords. Significant enrichment scores were found in several clusters, which were noteworthy since they included links to the progression of glioma. Particularly, the green and brown modules were involved in changes in signaling pathways directly affecting metabolism, while the yellow and gray modules mainly affect transcription regulation. All modules were involved in protein binding. These findings may indicate that the green and brown modules are linked with the same or closely related metabolic pathways, as they both operate in the plasma membrane. On the other hand, pathway enrichment results have shown the involvement of the yellow and gray modules in the PI3K/Akt pathway, although possibly in different processes since the gray module could be operating in a specific cytoplasmic organelle as compared to the yellow module, which was indicated to be in the free cytoplasm.

3.4. Identification of Protein–Protein Interaction Networks and Hub Genes

PPI networks were built for each module of interest using the STRING database to further investigate any possible connections between the different proteins corresponding to the genes within the discovered modules. To guarantee the accuracy of the anticipated interactions, a high confidence level of 0.7 was established [21]. In order to determine probable hub genes, each PPI network was imported into Cytoscape, and its network and gene (node) scores were calculated. In Cytoscape, CytoHubba was used to rank the genes within each network through the topological algorithm degree centrality, and the ten genes with the highest rankings for each network were considered as hub genes (Figure 6). The color of the genes represents the ranking of the hub genes, with the strongest red being the highest ranking. These identified hub genes require further analysis to determine their possible roles in the main pathways and processes associated with glioma due to their high interaction scores within the PPI networks of their corresponding modules, which indicates their crucial role in the gene network.

3.5. Signature-Based Drug Repurposing

In order to search for drugs that can reverse the expression of the top 10 hub genes from each module, transcriptional and molecular signatures were analyzed through the DRE webserver. Table 2 lists the top therapeutic candidates obtained together with the relevant mechanisms of action. Consideration was given to the top five potential drug candidates with the lowest false discovery rate (FDR) and Tau values, as compounds with low FDR tend to produce fewer false positive drug discovery results, and compounds with more negative Tau values induce more effective gene expression changes, countering the provided gene signature. For the upregulated hub genes, norgestimate, phentolamine, GW0742, olomoucine, and ambroxol were the top-ranked therapeutic candidates, whereas ethisterone, noscapine, nomegestrol, carmoxirole, and oxcarbazepine were the most promising drug candidates for the downregulated hub genes.

4. Discussion

4.1. Gene Co-Expression Modules across the Datasets

The WHO categorized the different grades of glioma based on histological features, molecular characteristics, and clinical behaviors. The vast heterogeneity of gliomas per grade and cell origin is one of the significant factors contributing to its difficulty responding well to treatment [27]. Through WGCNA, this study focused on identifying points at which the distinct glioma grades align genetically, since several studies have demonstrated the effectiveness of this approach in determining highly correlated genes in related diseases from a systems biology point of view [28,29,30]. Furthermore, it is able to determine pathways involved in disease networks due to its ability to investigate the interplay among the highly preserved genes. These reasons supported the use of WGCNA to elaborate the disease network present in different glioma grades, represented by PA, OG, AA, and GBM datasets. Since the module eigengenes of other datasets were projected onto the GBM dataset, it can be speculated that the increasing trend in the module preservation analysis from the PA dataset to the AA dataset (Figure 4) was related to the fact that GBM is at the highest severity of the disease. This can also be observed in Figure A3, where the ranked connectivity of GBM across the datasets increases from the PA dataset to the AA dataset. With the recent discovery of the possible lineage conversion of oligodendroglioma to astrocytoma [31], the network’s relation to glioma progression was further investigated. Figure 7 shows the top KEGG pathways per module (Table A2) that may potentially serve as pivotal points of glioma progression, as it demonstrates the increasing preservation of modules per glioma grade with reference to GBM. It is important to note that the modules that did not reach the threshold of high preservation (z > 10) may still induce effects in other signaling pathways depending on the PPI networks. Particularly, the deregulation of the PI3K/Akt pathway is known to cause changes in MAPK/ERK, calcium signaling, and cAMP signaling pathways; Ras and Rap1 proteins are known regulators of the MAPK/ERK pathway [32,33,34,35,36]. On the other hand, metabolite biomarkers in neurodegenerative diseases have been discovered in recent studies, showing their connection to metabolic pathways [37,38]. Moving on, high preservation of the Wnt/β-catenin signaling pathway in AA suggests its tendency to migrate [39], and the deregulation of the mRNA surveillance pathway poses challenges in neurodegeneration [40]. Lastly, phosphatidic acid is known to regulate the PI3K/Akt signaling pathway through mTOR [41]. Hence, it can be hypothesized that there are existing synergistic processes, cross-talks, or even signal amplification events involved in these pathways, which may be relevant to glioma progression. This study may open new avenues for future explorations in glioma research for studies involving computational, wet laboratory, and clinical experimentations.
However, only the co-expression modules green, yellow, brown, and gray were highly preserved in all the datasets and subjected to further analysis. The functional annotation and pathway enrichment analysis was able to deduce that the green and brown modules are primarily associated with signal transduction deregulation in metabolic pathways, while the yellow and gray modules are mainly involved with cell proliferation of glioma, specifically the PI3K/Akt pathway (Table A3, Table A4, Table A5 and Table A6).

4.2. Module Hub Genes and Their Protein Functions

4.2.1. Involvement of PI3K/Akt Pathway and Other Signaling Pathways

The hub genes of each highly preserved module have been identified and are listed along with their protein functions in Table A1. Notably, the EGFR gene was present across the green, brown, and gray modules while retaining its high interaction score, implying the strong involvement of the PI3K/Akt pathway, as implied by the pathway enrichment results. It has been seen as an important signaling pathway in the progression of low-grade to high-grade glioma since it is commonly upregulated through several types of genetic or post-translational mutations [42]. With this, many have made use of EGFR inhibitors, as most of the EGFR mutations in cancer cells cause overactivation of the PI3K-Akt pathway. However, the lack of blood–brain barrier permeability posed as its major obstacle [43,44]. Aside from EGFR, AKT1 was also included in the upregulated hub genes, indicating the pathway’s relevance to the present study. Overexpression of the AKT1 gene has been one of the common characteristics of glioma cells, as it promotes viability and malignancy [45]. On the other hand, the PRKACA and PRKACB genes that may indirectly inhibit the activation of the PI3K/Akt pathway through the cAMP signaling pathway were downregulated. This is consistent with the need to prevent apoptosis and ensure the survival of the gliomas.
The inclusion of hub genes primarily involved in other signaling pathways that affect the PI3K/Akt pathway further supports the existing pathway interconnectivity in glioma. It was found that KRAS and GRB2 in the MAPK/ERK signaling pathway were upregulated across the datasets. Recent studies have indicated the strong cooperativity between the PI3K/AKT and MAPK/ERK signaling pathways in glioma progression due to them being regulators of cell proliferation [32,34]. Aside from this, the CTNNB1 and FN1 genes in the Wnt/β-catenin signaling pathway were also upregulated. The relationship between the Wnt/β-catenin and PI3K/Akt signaling pathways has been known to induce glioma malignancy [39]. Lastly, it was observed that the calcium signaling pathway is significantly involved, represented by the CACNA1C, CACNG2, CALML3, CALM1, and CAMK2A hub genes. Its dysregulation affects the PI3K/Akt pathway by increasing calcium influx, which further activates biological processes such as growth, gene expression, and neurotransmitter release, and it is claimed that the calcium signaling pathway is another relevant factor contributing to the motility and invasion of gliomas [33,46].

4.2.2. Deregulation of Cellular Processes in Glioma

Several upregulated hub genes relating to DNA repair have been observed in the gray module, such as EXO1, H3C12, and ATM. In cancer cells, the upregulation of genes involved in DNA repair may indicate increased DNA damage that may possibly be due to genomic instability, increased replication stress, and resistance to DNA-damaging therapies [47]. Specifically, high expression of ATM gene has been shown to be linked with radio-resistance of glioma [48]. On the other hand, upregulated hub genes in the yellow module partake mainly on cell cycle regulation. The TP53 gene, which encodes for the p53 protein, which is a main regulator of both the cell cycle and DNA repair and replication, was found to be upregulated across the datasets. It functions by inducing the expression of DNA repair genes while inhibiting cell cycle progression in response to DNA damage to ensure genomic integrity, although cancer cells normally observe upregulated expression of TP53 due to mutations that enable them to become chemo-resistant [49]. In general, the upregulation of genes involved in the cell cycle may promote cell cycle progression, genomic instability, and enhanced proliferation, which could have an impact on the behavior of tumors and their response to therapeutics. Across the modules, hub genes that are involved in immune response were also observed. Among these, genes that code for clusters of differentiation proteins such as CD4, CD2, and CD8A were found to be downregulated and observed in the green module, while CD8A was also present in the brown and gray modules. The downregulation of these genes may indicate the impairment of immune surveillance and effector functions of T cells and natural killer cells that overall modulate immune responses against anomalies, including tumors. Figure 8 shows a graphical summary of the hub genes and the potential cellular processes affected due to their involvement in glioma. Furthermore, the involvement of specific hub genes in several cellular processes supports the existing pathway interconnectivity.

4.3. Metabolic Reprogramming of Glioma Cells

Collectively, the presented results indicate that the gene network constructed is primarily involved with metabolic reprogramming of glioma, as the modules were found to be concerned with the PI3K/Akt and metabolic pathways. Metabolic reprogramming is a major hallmark in glioma, as it is when tumor cells alter their metabolism into aerobic glycolysis to generate the necessary energy to invade surrounding cells and adapt to different microenvironments [50,51]. Aside from metabolizing glucose for energy production, cancer cells are able to utilize lipids as raw materials for energetic currency [52]. In particular, the green and brown modules were highly involved with metabolic pathways and operate in cytosol, which aligns with malignant glioma cells also being dependent on fatty acid oxidation [53,54]. Furthermore, these metabolic changes are known to be regulated by phosphoinositide 3-kinase/protein kinase (PI3Ks), the main ligand for the PI3K/Akt pathway, as it can directly regulate the TCA cycle, wherein increased PI3K activation also increases AKT activity, which promotes citrate conversion to acetyl-CoA for fatty acid synthesis [51].

4.4. Signature-Based Drug Repurposing

The involvement of several signaling pathways affecting the PI3K/Akt pathway creates avenues into which therapeutics may be introduced. Through the DRE webserver, potential drug candidates were screened, as shown in Table 3, along with their corresponding pathways and affected processes based on their mechanism of action. Norgestimate, ethisterone, and nomegestrol are synthetic progestin hormones that mimic the effects of progesterone to activate physiological responses through binding with progesterone receptors [55,56,57]. The influence of sex hormones on glioma progression was observed in women, who are otherwise less likely to develop glioma than men but who become more prone to developing glioma after menopause due to decreased estrogen production [58,59,60]. Progesterone, on the other hand, can indirectly modulate estrogen production by interfering with estrogen receptor activity or releasing gonadotropins through feedback mechanisms [61]. Despite producing pro-tumorigenic effects at low doses, high doses of progesterone have demonstrated therapeutic effects against GBM by downregulating the activity of progesterone receptor B (PR-B), a protein that promotes tumor cell growth, or through alteration of detoxification mechanisms, stress, immune response, and glucose metabolism [62,63,64,65,66,67]. Furthermore, a significant decrease in the expression levels of EGFR, Akt, phospho-Akt, mTOR, and phospho-mTOR has been observed following high doses of progesterone, indicating its potential attenuating effect against GBM through the PI3K/Akt and metabolic pathways [68,69,70]. Aside from this, progesterone has a known regulatory effect in the significant signaling pathways previously mentioned, in cAMP, MAPK/ERK, Wnt/β-catenin, and calcium signaling pathway in other types of cancer [71,72,73,74], which further demonstrates that the use of synthetic progestin hormones against glioma progression can be a potential area of interest.
Phentolamine, an adrenergic receptor antagonist drug, also scored significantly, along with other approved and experimental drugs, as shown in Table 2. Previous studies have shown the regulating effect of adrenergic receptor antagonist drugs against cancer initiation and progression due to their influence on the cAMP/PKA, MAPK/ERK, and PI3K/Akt pathways [75,76]. Another drug candidate is GW0742, as its potential modulating effect towards fatty acid oxidation has become an area of interest for cancer treatment research [77,78]. A CDK inhibitor, olomoucine, was also included as a potential drug candidate due to its regulation of the cell cycle. CDK inhibitors promote cell cycle arrest and apoptosis and decrease cell proliferation in glioma [79]. Ambroxol and oxcarbazepine were found to be possible drug candidates due to their disrupting effect on the depolarization process and ion homeostasis, in which calcium signaling is also affected [80]. On the other hand, noscapine, a bradykinin receptor antagonist, also interferes with the calcium signaling pathway by reducing calcium production and, subsequently, glioma cell migration [81]. Lastly, a dopamine receptor agonist, carmoxirole, was included due to its potential to modulate cAMP production, with dopamine receptor D1 inhibiting GBM tumorigenicity and dopamine receptor D2 modulating GBM survival and death [82,83].

5. Conclusions

The present study was able to successfully determine the four highly conserved modules in different glioma grades that were represented by PA (GSE50161), OG (GSE4290), AA (GSE43378), and GBM (GSE36245). Through WGCNA, the top 10 hub genes per module were identified and used to screen for potential drug candidates. Furthermore, the functional annotation and pathway enrichment analysis paved the way for deducing the complex molecular basis for glioma, focusing on the PI3K/Akt pathway, signal transduction deregulation, and metabolic reprogramming. Based on the hub genes, several signaling pathways associated with the PI3K/Akt pathway, such as the cAMP signaling pathway, MAPK/ERK, Wnt/β-catenin, and calcium signaling, have also been deregulated, which often results in metabolic reprogramming. The involvement of several signaling pathways in deregulating the PI3K/Akt pathway and vice versa further indicates the possible interconnection of signaling pathways in glioma severity. The synthetic progestin hormones norgestimate and ethisterone were the top drug candidates identified in the transcriptional signature-based drug repurposing approach. Aside from these, several experimental and approved drug candidates have also been identified, including an adrenergic receptor antagonist, a PPAR-γ receptor agonist, a CDK inhibitor, a sodium channel blocker, and a bradykinin receptor antagonist, which further demonstrates that the gene network is a potential therapeutic avenue for glioma. Our study posts limitations since all were conducted in silico; however, the theoretical data presented can be used as a platform to explore avenues in designing drug cocktails or drug candidates for glioma. Furthermore, this study provides insights on the interconnectivities of molecular signaling, which can be further evaluated using wet laboratory experimentations.

Supplementary Materials

The following supporting information can be downloaded at www.mdpi.com/article/10.3390/biology13040206/s1: Figure S1: TOM-based dissimilarity clustering of (a) PA, (b) OG, (c) AA, and (d) GBM datasets; Figure S2: Summary of in-module connectivities of genes in highly preserved modules across all datasets referenced in GBM and plotted against (a) PA, (b) OG, and (c) AA datasets.

Author Contributions

Conceptualization, M.A.O. and L.L.T.; methodology, M.A.O. and L.L.T.; validation, P.M.P.T.F. and L.L.T.; formal analysis, M.A.O.; investigation, M.A.O. and P.M.P.T.F.; resources, L.L.T.; data curation, M.A.O.; writing—original draft preparation, M.A.O.; writing—review and editing, P.M.P.T.F. and L.L.T.; visualization, M.A.O.; supervision, L.L.T.; funding acquisition, L.L.T.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The gene microarray datasets used for the study are openly available in the NCBI Gene Expression Omnibus (GEO) database under the accession IDs GSE50161, GSE4290, GSE43378, and GSE36245 at https://www.ncbi.nlm.nih.gov/geo/ accessed on 7 December 2023.

Acknowledgments

The authors would like to express their gratitude to Ryan Christian Mailem and Mathew Timothy Manuel for their guidance in the making of the study.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Figure A1. Boxplot of microarray datasets (a) PA, (b) OG, and (c) AA and (d) GBM datasets after data preparation and filtering.
Figure A1. Boxplot of microarray datasets (a) PA, (b) OG, and (c) AA and (d) GBM datasets after data preparation and filtering.
Biology 13 00206 g0a1
Figure A2. Sample clustering dendrogram of (a) PA, (b) OG, (c) AA, and (d) GBM datasets. The ends of each branch represent the individual samples in the datasets.
Figure A2. Sample clustering dendrogram of (a) PA, (b) OG, (c) AA, and (d) GBM datasets. The ends of each branch represent the individual samples in the datasets.
Biology 13 00206 g0a2
Figure A3. Ranked expression plots and ranked connectivity plots for dataset comparison: (a,b) PA vs. OG; (c,d) PA vs. AA; (e,f) OG vs. AA; (g,h) GBM vs. PA; (i,j) GBM vs. AA; and (k,l) GBM vs. OG.
Figure A3. Ranked expression plots and ranked connectivity plots for dataset comparison: (a,b) PA vs. OG; (c,d) PA vs. AA; (e,f) OG vs. AA; (g,h) GBM vs. PA; (i,j) GBM vs. AA; and (k,l) GBM vs. OG.
Biology 13 00206 g0a3
Table A1. Summary of module hub genes and their protein functions obtained from the GeneCards database (www.genecards.org (accessed 27 December 2023)) [84].
Table A1. Summary of module hub genes and their protein functions obtained from the GeneCards database (www.genecards.org (accessed 27 December 2023)) [84].
GeneProteinFunctionModule/s
EGFRepidermal growth factor receptoractivates MAPK and PI3K/AKT pathways; promotes cell proliferation and survivalGreen, brown, and grey
CD8ACD8 subunit alphabinds to MHC class I; facilitates T cell-mediated cytotoxicity against infected or transformed cellsGreen, brown, and grey
CALML3calmodulin like 3regulates calmodulin-dependent signaling; potential roles in intracellular calcium homeostasisGreen and brown
PRKACAprotein kinase cAMP-activated catalytic subunit alphaphosphorylates target proteins; involved in cAMP signaling; regulates various cellular processesGreen and grey
CACNA1Ccalcium voltage-gated channel subunit alpha1 Cmediates calcium ion influx; regulates neuronal excitability and cardiac functionBrown
CAMK2Acalcium/calmodulin dependent protein kinase II alpharegulates synaptic plasticity and memory formationBrown
PRKACBprotein kinase cAMP-activated catalytic subunit betainvolved in cAMP signaling; regulates various cellular processesBrown
CACNG2calcium voltage-gated channel auxiliary subunit gamma 2modulates channel function; regulates calcium ion influx and neuronal excitabilityBrown
KRASKRAS proto-oncogene, GTPaseactivates MAPK pathway; regulates cell proliferation, survival, and differentiationBrown
TNFtumor necrosis factorregulates inflammation, apoptosis, and immune cell functionsBrown
GRIN2Bglutamate ionotropic receptor NMDA type subunit 2Bmediates synaptic transmission; regulates synaptic plasticity, learning, and memoryBrown
ITGAMintegrin subunit alpha Mregulates immune cell adhesion, migration, and phagocytosisGreen
HSP90AA1heat shock protein 90 alpha family class A member 1assists in protein folding and stabilizationGreen
GRB2growth factor receptor bound protein 2links receptor tyrosine kinases to Ras/MAPK pathway activationGreen
TLR4Toll-like receptor 4detects lipopolysaccharides; activates innate immune responsesGreen
CD4CD4 moleculebinds to MHC class II; facilitates antigen recognition and T cell activationGreen
CD2CD2 moleculeregulates immune cell interactions; facilitates T cell adhesion and co-stimulationGreen
CALM1calmodulin 1modulates calcium-dependent signaling pathways; regulates cellular responses to calcium ionsGrey
EXO1exonuclease 1involved in DNA mismatch repair and recombination; maintains genomic stabilityGrey
KITKIT proto-oncogene, receptor tyrosine kinaseactivates PI3K/AKT and MAPK pathways; regulates cell proliferation, survival, and differentiationGrey
CREB1cAMP responsive element binding protein 1binds to cAMP response elements; regulates gene expression in response to cAMP signalingGrey
GNASG-protein alpha subunitactivates adenylyl cyclase; regulates cellular signaling pathways via GPCRsGrey
ATMATM serine/threonine kinaseregulates DNA damage response; activates cell cycle checkpoints and DNA repair mechanismsGrey
H3C12H3 clustered histone 12packaging and organizing of DNA into nucleosomeGrey
FN1fibronectin 1mediates cell adhesion and migration; regulates tissue remodeling and repairYellow
AKT1AKT serine/threonine kinase 1activates mTOR and other downstream pathways; regulates cell survival, growth, and metabolismYellow
BUB1BBUB1 mitotic checkpoint serine/threonine kinase Bregulates chromosome alignment and segregation; ensures genomic stability during cell divisionYellow
CCNA2cyclin A2forms complexes with CDKs; regulates G1/S and G2/M transitions of the cell cycleYellow
CTNNB1catenin beta 1component of adherens junctions; regulates cell-cell adhesion and Wnt signaling pathwayYellow
CDC20cell division cycle 20facilitates ubiquitination and degradation of cell cycle regulators; controls mitotic progressionYellow
TP53tumor protein p53activates DNA repair or apoptosis in response to DNA damage; regulates cell cycle checkpointsYellow
CCNB1cyclin B1forms complexes with CDK1; regulates the G2/M transition of the cell cycleYellow
TOP2ADNA topoisomerase II alpharegulates DNA topology during replication and transcription; targeted in cancer therapyYellow
CDK1cyclin dependent kinase 1forms complexes with cyclins; regulates cell cycle transitions and mitotic entryYellow
Table A2. Top KEGG pathway analysis results in identified modules obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20].
Table A2. Top KEGG pathway analysis results in identified modules obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20].
ModuleTermCountAdj. p-Value
Greenhsa01100 metabolic pathways1941.4 × 10−3
Brownhsa01100 metabolic pathways1721.3 × 10−2
Yellowhsa04151 PI3K/Akt signaling pathway 596.3 × 10−4
Grayhsa04151 PI3K/Akt signaling pathway562.9 × 10−3
Purplehsa04020 calcium signaling pathway652.8 × 10−13
Tanhsa04015 Rap1 signaling pathway292.0 × 10−2
Cyanhsa04020 calcium signaling pathway513.6 × 10−7
Blackhsa04020 calcium signaling pathway553.9 × 10−9
Magentahsa04014 Ras signaling pathway349.2 × 10−4
Pinkhsa04024 cAMP signaling pathway458.1 × 10−6
Salmonhsa04010 MAPK signaling pathway526.1 × 10−6
Turquoisehsa04010 MAPK signaling pathway381.3 × 10−2
Bluehsa04020 calcium signaling pathway681.9 × 10−10
Gray60hsa04010 MAPK signaling pathway531.2 × 10−5
Midnight Bluehsa04151 PI3K/Akt signaling pathway514.2 × 10−3
Redhsa04115 p53 signaling pathway211.4 × 10−5
Light Cyanhsa05022 pathways of neurodegeneration-multiple diseases592.3 × 10−2
Green-yellowhsa03015 mRNA surveillance pathway318.8 × 10−9
Light Yellowhsa04310 Wnt signaling pathway259.8 × 10−3
Light Greenhsa05200 pathways in cancer751.2 × 10−3
Goldhsa04072 phospholipase D signaling pathway265.1 × 10−3
Table A3. Top GO and pathway analysis results in brown module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Table A3. Top GO and pathway analysis results in brown module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Category TermCount Adj. p-Value
BPGO:0007165 signal transduction150 7.8 × 10−5
GO:0035556 intracellular signal transduction694.4 × 10−6
GO:0006468 protein phosphorylation643.1 × 10−5
GO:0006508 proteolysis622.5 × 10−4
GO:0007268 chemical synaptic transmission569.4 × 10−11
CCGO:0005886 plasma membrane6673.8 × 10−29
GO:0016020 integral component of membrane5884.4 × 10−12
GO:0005737 cytoplasm5599.7 × 10−5
GO:0005829 cytosol5432.0 × 10−4
GO:0016020 membrane4241.9 × 10−12
MFGO:0005515 protein binding11925.4 × 10−4
GO:0005524 ATP binding1801.1 × 10−4
GO:0042802 identical protein binding1746.4 × 10−2
GO:0005509 calcium ion binding1055.3 × 10−6
GO:0004712 protein serine/threonine/tyrosine kinase activity681.4 × 10−5
KEGGhsa01100 metabolic pathways1721.3 × 10−2
hsa05200 pathways in cancer689.3 × 10−3
hsa05022 pathways of neurodegeneration-multiple diseases602.0 × 10−2
hsa04020 calcium signaling pathway492.1 × 10−6
hsa04010 MAPK signaling pathway484.3 × 10−4
Table A4. Top GO and pathway analysis results in green module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Table A4. Top GO and pathway analysis results in green module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Category TermCount Adj. p-Value
BPGO:0007165 signal transduction1814.7 × 10−10
GO:0030154 cell differentiation891.1 × 10−3
GO:0007155 cell adhesion873.6 × 10−7
GO:0045087 innate immune response832.6 × 10−4
GO:0006915 apoptotic process754.2 × 10−3
CCGO:0005886 plasma membrane7191.0 × 10−37
GO:0005737 cytoplasm6304.9 × 10−11
GO:0016020 integral component of membrane5953.3 × 10−10
GO:0005829 cytosol5694.2 × 10−5
GO:0016020 membrane4312.6 × 10−11
MFGO:0005515 protein binding12563.3 × 10−7
GO:0042802 identical protein binding2144.2 × 10−6
GO:0005524 ATP binding1724.2 × 10−3
GO:0005509 calcium ion binding1181.0 × 10−8
GO:0042803 protein homodimerization activity932.3 × 10−3
KEGGhsa01100 metabolic pathways1941.4 × 10−3
hsa05200 pathways in cancer821.6 × 10−4
hsa04010 MAPK signaling561.5 × 10−5
hsa04020 calcium signaling537.3 × 10−7
hsa04151 PI3K/Akt signaling pathway511.6 × 10−2
Table A5. Top GO and pathway analysis results in gray module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Table A5. Top GO and pathway analysis results in gray module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Category TermCount Adj. p-Value
BPGO:0045944 positive regulation of transcription by RNA polymerase II160 3.9 × 10−6
GO:0007165 signal transduction1582.0 × 10−4
GO:0000122 negative regulation of transcription by RNA polymerase II1169.2 × 10−3
GO:0045893 positive regulation of DNA-templated transcription971.6 × 10−4
GO:0007155 cell adhesion855.9 × 10−6
CCGO:0005737 cytoplasm6287.0 × 10−9
GO:0005886 plasma membrane6195.8 × 10−12
GO:0005829 cytosol5961.3 × 10−6
GO:0005634 nucleus5882.7 × 10−2
GO:0016020 integral component of membrane5216.2 × 10−2
MFGO:0005515 protein binding13078.6 × 10−9
GO:0046872 metal ion binding2882.3 × 10−2
GO:0005524 ATP binding2171.1 × 10−9
GO:0005509 calcium ion binding1105.5 × 10−6
GO:0004712 protein serine/threonine/tyrosine kinase activity781.3 × 10−7
KEGGhsa04151 PI3K/Akt signaling pathway562.9 × 10−3
hsa04010 MAPK signaling pathway497.3 × 10−4
hsa04020 calcium signaling pathway482.0 × 10−5
hsa04015 Rap1 signaling pathway384.4 × 10−4
hsa04024 cAMP signaling pathway365.5 × 10−3
Table A6. Top GO and pathway analysis results in yellow module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Table A6. Top GO and pathway analysis results in yellow module obtained from the Gene Ontology Consortium (http://www.geneontology.org (accessed 27 December 2023)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg (accessed 27 December 2023)) [20,85].
Category TermCount Adj. p-Value
BPGO:0045944 positive regulation of transcription by RNA polymerase II155 2.3 × 10−6
GO:0000122 negative regulation of transcription by RNA polymerase II1204.8 × 10−4
GO:0045893 positive regulation of DNA-templated transcription1039.5 × 10−7
GO:0006355 regulation of DNA-templated transcription1036.9 × 10−2
GO:0006915 apoptotic process994.0 × 10−9
CCGO:0005829 cytosol6813.1 × 10−29
GO:0005634 nucleus6797.3 × 10−19
GO:0005737 cytoplasm6231.5 × 10−13
GO:0005654 nucleoplasm5548.4 × 10−36
GO:0016020 membrane4751.1 × 10−24
MFGO:0005515 protein binding13963.1 × 10−46
GO:0046872 metal ion binding2687.0 × 10−2
GO:0005524 ATP binding2321.3 × 10−15
GO:0042802 identical protein binding2274.2 × 10−9
GO:0003723 RNA binding2235.1 × 10−15
KEGGhsa04151 PI3K/Akt signaling pathway 596.3 × 10−4
hsa05200 pathways in cancer532.7 × 10−3
hsa04141 protein processing in endoplasmic reticulum483.1 × 10−10
hsa05205 proteoglycans in cancer447.7 × 10−6
hsa04010 MAPK signaling pathway442.5 × 10−2

References

  1. de Almeida Sassi, F.; Lunardi Brunetto, A.; Schwartsmann, G.; Roesler, R.; Abujamra, A.L. Glioma Revisited: From Neurogenesis and Cancer Stem Cells to the Epigenetic Regulation of the Niche. J. Oncol. 2012, 2012, 537861. [Google Scholar] [CrossRef]
  2. Torp, S.H.; Solheim, O.; Skjulsvik, A.J. The WHO 2021 Classification of Central Nervous System Tumours: A Practical Update on What Neurosurgeons Need to Know—A Minireview. Acta Neurochir. 2022, 164, 2453–2464. [Google Scholar] [CrossRef]
  3. van den Bent, M.J.; Brandes, A.A.; Taphoorn, M.J.B.; Kros, J.M.; Kouwenhoven, M.C.M.; Delattre, J.-Y.; Bernsen, H.J.J.A.; Frenay, M.; Tijssen, C.C.; Grisold, W.; et al. Adjuvant Procarbazine, Lomustine, and Vincristine Chemotherapy in Newly Diagnosed Anaplastic Oligodendroglioma: Long-Term Follow-Up of EORTC Brain Tumor Group Study 26951. J. Clin. Oncol. 2012, 31, 344–350. [Google Scholar] [CrossRef]
  4. Johnson, D.R.; Diehn, F.E.; Giannini, C.; Jenkins, R.B.; Jenkins, S.M.; Parney, I.F.; Kaufmann, T.J. Genetically Defined Oligodendroglioma Is Characterized by Indistinct Tumor Borders at MRI. Am. J. Neuroradiol. 2017, 38, 678. [Google Scholar] [CrossRef]
  5. Cohen, A.L.; Holmen, S.L.; Colman, H. IDH1 and IDH2 Mutations in Gliomas. Curr. Neurol. Neurosci. Rep. 2013, 13, 345. [Google Scholar] [CrossRef]
  6. Carter, T.C.; Medina-Flores, R.; Lawler, B.E. Glioblastoma Treatment with Temozolomide and Bevacizumab and Overall Survival in a Rural Tertiary Healthcare Practice. BioMed Res. Int. 2018, 2018, 6204676. [Google Scholar] [CrossRef]
  7. Pan, P.C.; Magge, R.S. Mechanisms of EGFR Resistance in Glioblastoma. Int. J. Mol. Sci. 2020, 21, 8471. [Google Scholar] [CrossRef]
  8. Liu, X.; Chen, X.; Shi, L.; Shan, Q.; Cao, Q.; Yue, C.; Li, H.; Li, S.; Wang, J.; Gao, S.; et al. The Third-Generation EGFR Inhibitor AZD9291 Overcomes Primary Resistance by Continuously Blocking ERK Signaling in Glioblastoma. J. Exp. Clin. Cancer Res. 2019, 38, 219. [Google Scholar] [CrossRef]
  9. Nasir, S.; Nazir, S.; Hanif, R.; Javed, A. Glioblastoma Multiforme: Probing Solutions to Systemic Toxicity towards High-Dose Chemotherapy and Inflammatory Influence in Resistance against Temozolomide. Pharmaceutics 2023, 15, 687. [Google Scholar] [CrossRef]
  10. Xie, X.; Bao, S.; Zhao, H.; Li, L.; Fu, X. Efficacy and Safety of Bevacizumab for Treating Glioblastoma: A Systematic Review and Meta-Analysis of Phase II and III Randomized Controlled Trials. Cancer Investig. 2023, 41, 305–317. [Google Scholar] [CrossRef]
  11. Reardon, D.A.; Nabors, L.B.; Mason, W.P.; Perry, J.R.; Shapiro, W.; Kavan, P.; Mathieu, D.; Phuphanich, S.; Cseh, A.; Fu, Y.; et al. Phase I/Randomized Phase II Study of Afatinib, an Irreversible ErbB Family Blocker, with or without Protracted Temozolomide in Adults with Recurrent Glioblastoma. Neuro Oncol. 2015, 17, 430–439. [Google Scholar] [CrossRef]
  12. Ostrom, Q.T.; Cioffi, G.; Waite, K.; Kruchko, C.; Barnholtz-Sloan, J.S. CBTRUS Statistical Report: Primary Brain and Other Central Nervous System Tumors Diagnosed in the United States in 2014–2018. Neuro Oncol. 2021, 23, iii1–iii105. [Google Scholar] [CrossRef]
  13. Low, J.T.; Ostrom, Q.T.; Cioffi, G.; Neff, C.; Waite, K.A.; Kruchko, C.; Barnholtz-Sloan, J.S. Primary Brain and Other Central Nervous System Tumors in the United States (2014–2018): A Summary of the CBTRUS Statistical Report for Clinicians. Neuro-Oncol. Pract. 2022, 9, 165–182. [Google Scholar] [CrossRef]
  14. Liu, H.; Bebu, I.; Li, X. Microarray Probes and Probe Sets. Front. Biosci. (Elite Ed.) 2010, 2, 325–338. [Google Scholar] [CrossRef]
  15. Junet, V.; Farrés, J.; Mas, J.M.; Daura, X. CuBlock: A Cross-Platform Normalization Method for gene-Expression Microarrays. Bioinformatics 2021, 37, 2365–2373. [Google Scholar] [CrossRef]
  16. Yu, Y.; Zhang, N.; Mai, Y.; Ren, L.; Chen, Q.; Cao, Z.; Chen, Q.; Liu, Y.; Hou, W.; Yang, J.; et al. Correcting Batch Effects in Large-Scale Multiomics Studies Using a Reference-Material-Based Ratio Method. Genome Biol. 2023, 24, 201. [Google Scholar] [CrossRef]
  17. Zhu, Y.; Yang, X.; Zu, Y. Integrated Analysis of WGCNA and Machine Learning Identified Biomarkers in Dilated Cardiomyopathy with Heart. Front. Cell. Dev. Biol. 2022, 10, 1089915. [Google Scholar] [CrossRef]
  18. Langfelder, P.; Zhang, B.; Horvath, S. Defining Clusters from a Hierarchical Cluster Tree: The Dynamic Tree Cut Package for R. Bioinformatics 2008, 24, 719–720. [Google Scholar] [CrossRef]
  19. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A Web Server for Functional Enrichment Analysis and Functional Annotation of Gene Lists (2021 Update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef]
  20. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  21. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein–Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  22. Chin, C.-H.; Chen, S.-H.; Wu, H.-H.; Ho, C.-W.; Ko, M.-T.; Lin, C.-Y. CytoHubba: Identifying Hub Objects and Sub-Networks from Complex Interactome. BMC Syst. Biol. 2014, 8, S11. [Google Scholar] [CrossRef]
  23. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  24. Li, X.; Pan, L.; Sanchez-Burgos, L.; Hühn, D.; Fernandez-Capetillo, O. The Drug Repurposing Encyclopedia (DRE): A Web Server for Systematic Drug Repurposing across 20 Organisms. bioRxiv 2023. [Google Scholar] [CrossRef]
  25. Lamb, J.; Crawford, E.D.; Peck, D.; Modell, J.W.; Blat, I.C.; Wrobel, M.J.; Lerner, J.; Brunet, J.-P.; Subramanian, A.; Ross, K.N.; et al. The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease. Science 2006, 313, 1929–1935. [Google Scholar] [CrossRef] [PubMed]
  26. Miller, J.A.; Horvath, S.; Geschwind, D.H. Divergence of Human and Mouse Brain Transcriptome Highlights Alzheimer Disease Pathways. Proc. Natl. Acad. Sci. USA 2010, 107, 12698–12703. [Google Scholar] [CrossRef] [PubMed]
  27. Hutóczki, G.; Virga, J.; Birkó, Z.; Klekner, A. Novel Concepts of Glioblastoma Therapy Concerning Its Heterogeneity. Int. J. Mol. Sci. 2021, 22, 10005. [Google Scholar] [CrossRef] [PubMed]
  28. Mailem, R.C.; Tayo, L.L. Identification of Hub Genes and Key Pathways in TNF-α and IFN-γ Induced Cytokine Storms via Bioinformatics. In Proceedings of the 2022 10th International Conference on Bioinformatics and Computational Biology, ICBCB 2022, Hangzhou, China, 13–15 May 2022; pp. 6–12. [Google Scholar]
  29. Mailem, R.C.; Tayo, L.L. Drug Repurposing Using Gene Co-Expression and Module Preservation Analysis in Acute Respiratory Distress Syndrome (ARDS), Systemic Inflammatory Response Syndrome (SIRS), Sepsis, and COVID-19. Biology 2022, 11, 1827. [Google Scholar] [CrossRef] [PubMed]
  30. Manuel, M.T.A.; Tayo, L.L. Navigating the Gene Co-Expression Network and Drug Repurposing Opportunities for Brain Disorders Associated with Neurocognitive Impairment. Brain Sci. 2023, 13, 1564. [Google Scholar] [CrossRef] [PubMed]
  31. Kim, J.-H.; Jang, W.-Y.; Jung, T.-Y.; Jung, S.; Kim, K.-K.; Kim, H.-S.; Kim, E.-H.; Lee, M.-C.; Moon, K.-S.; Lee, K.-H. Recurrent Glioma with Lineage Conversion from oligodendroglioma to Astrocytoma in Two Cases. Front. Oncol. 2019, 9, 828. [Google Scholar] [CrossRef]
  32. Vitucci, M.; Karpinich, N.O.; Bash, R.E.; Werneke, A.M.; Schmid, R.S.; White, K.K.; McNeill, R.S.; Huff, B.; Wang, S.; Van Dyke, T.; et al. Cooperativity between MAPK and PI3K Signaling Activation Is Required for Glioblastoma Pathogenesis. Neuro Oncol. 2013, 15, 1317–1329. [Google Scholar] [CrossRef]
  33. Simms, B.A.; Zamponi, G.W. Neuronal Voltage-Gated Calcium Channels: Structure, Function, and Dysfunction. Neuron 2014, 82, 24–45. [Google Scholar] [CrossRef]
  34. Zhao, H.; Wang, J.; Shao, W.; Wu, C.; Chen, Z.; To, S.T.; Li, W. Recent Advances in the Use of PI3K Inhibitors for Glioblastoma Multiforme: Current Preclinical and Clinical Development. Mol. Cancer 2017, 16, 100. [Google Scholar] [CrossRef]
  35. Gündüz, D.; Troidl, C.; Tanislav, C.; Rohrbach, S.; Hamm, C.; Aslam, M. Role of PI3K/Akt and MEK/ERK Signalling in cAMP/Epac-Mediated Endothelial Barrier Stabilisation. Front. Physiol. 2019, 10, 1387. [Google Scholar] [CrossRef]
  36. York, R.D.; Yao, H.; Dillon, T.; Ellig, C.L.; Eckert, S.P.; McCleskey, E.W.; Stork, P.J. Rap1 Mediates Sustained MAP Kinase Activation Induced by Nerve Factor. Nature 1998, 392, 622–626. [Google Scholar] [CrossRef]
  37. Kori, M.; Aydn, B.; Unal, S.; Arga, K.Y.; Kazan, D. Metabolic Biomarkers and Neurodegeneration: A Pathway Enrichment of Alzheimer’s Disease, Parkinson’s Disease, and amyotrophic Lateral Sclerosis. OMICS 2016, 20, 645–661. [Google Scholar] [CrossRef]
  38. Sahebjam, S.; McNamara, M.G.; Mason, W.P. Emerging Biomarkers in Anaplastic Oligodendroglioma: Implications for Clinical Investigation and Patient Management. CNS Oncol. 2013, 2, 351–358. [Google Scholar] [CrossRef]
  39. Barzegar Behrooz, A.; Talaie, Z.; Jusheghani, F.; Łos, M.J.; Klonisch, T.; Ghavami, S. Wnt and PI3K/Akt/MTOR Survival Pathways as Therapeutic Targets in Glioblastoma. Int. J. Mol. Sci. 2022, 23, 1353. [Google Scholar] [CrossRef]
  40. Weskamp, K.; Barmada, S.J. RNA Degradation in Neurodegenerative Disease. Adv. Neurobiol. 2018, 20, 103–142. [Google Scholar]
  41. Bruntz, R.C.; Lindsley, C.W.; Brown, H.A. Phospholipase D Signaling Pathways and Phosphatidic Acid as therapeutic Targets in Cancer. Pharmacol. Rev. 2014, 66, 1033–1079. [Google Scholar] [CrossRef]
  42. Robinson, G.; Parker, M.; Kranenburg, T.A.; Lu, C.; Chen, X.; Ding, L.; Phoenix, T.N.; Hedlund, E.; Wei, L.; Zhu, X.; et al. Novel Mutations Target Distinct Subgroups of Medulloblastoma. Nature 2012, 488, 43–48. [Google Scholar] [CrossRef]
  43. Westphal, M.; Maire, C.L.; Lamszus, K. EGFR as a Target for Glioblastoma Treatment: An Unfulfilled Promise. CNS Drugs 2017, 31, 723–735. [Google Scholar] [CrossRef]
  44. Kizilbash, S.H. Why Has Targeting EGFR Aberrations in Glioblastoma Therapy Had Limited Success? Expert. Rev. Anticancer. Ther. 2022, 22, 1261–1263. [Google Scholar] [CrossRef]
  45. Schlegel, J.; Piontek, G.; Budde, B.; Neff, F.; Kraus, A. The Akt/Protein Kinase B-Dependent Anti-Apoptotic Pathway and the Mitogen-Activated Protein Kinase Cascade Are Alternatively Activated in Human Glioblastoma Multiforme. Cancer Lett. 2000, 158, 103–108. [Google Scholar] [CrossRef]
  46. So, J.-S.; Kim, H.; Han, K.-S. Mechanisms of Invasion in Glioblastoma: Extracellular Matrix, Ca2+ Signaling, and Glutamate. Front. Cell. Neurosci. 2021, 15, 663092. [Google Scholar] [CrossRef]
  47. Li, L.; Guan, Y.; Chen, X.; Yang, J.; Cheng, Y. DNA Repair Pathways in Cancer Therapy and Resistance. Front. Pharmacol. 2021, 11, 629266. [Google Scholar] [CrossRef]
  48. Li, Y.; Li, L.; Li, B.; Wu, Z.; Wu, Y.; Wang, Y.; Jin, F.; Li, D.; Ma, H.; Wang, D. Silencing of Ataxia-Telangiectasia Mutated by SiRNA Enhances the in Vitro and in Vivo Radiosensitivity of Glioma. Oncol. Rep. 2016, 35, 3303–3312. [Google Scholar] [CrossRef]
  49. Ozaki, T.; Nakagawara, A. Role of P53 in Cell Death and Human Cancers. Cancers 2011, 3, 994–1013. [Google Scholar] [CrossRef]
  50. Garcia, J.H.; Jain, S.; Aghi, M.K. Metabolic Drivers of Invasion in Glioblastoma. Front. Cell. Dev. Biol. 2021, 9, 683276. [Google Scholar] [CrossRef]
  51. Agnihotri, S.; Zadeh, G. Metabolic Reprogramming in Glioblastoma: The Influence of Cancer Metabolism on Epigenetics and Unanswered Questions. Neuro Oncol. 2016, 18, 160–172. [Google Scholar] [CrossRef]
  52. Strickland, M.; Stoll, E.A. Metabolic Reprogramming in Glioma. Front. Cell. Dev. Biol. 2017, 5, 43. [Google Scholar] [CrossRef]
  53. Lin, H.; Patel, S.; Affleck, V.S.; Wilson, I.; Turnbull, D.M.; Joshi, A.R.; Maxwell, R.; Stoll, E.A. Fatty Acid Oxidation Is Required for the Respiration and Proliferation of Malignant Glioma Cells. Neuro Oncol. 2017, 19, 43–54. [Google Scholar] [CrossRef]
  54. Ward, P.S.; Thompson, C.B. Metabolic Reprogramming: A Cancer Hallmark Even Warburg Did Not Anticipate. Cancer Cell 2012, 21, 297–308. [Google Scholar] [CrossRef]
  55. Carrera, A.R.M.; Eleazar, E.G.; Caparanga, A.R.; Tayo, L.L. Theoretical Studies on the Quantitative Structure–Toxicity Relationship of Polychlorinated Biphenyl Congeners Reveal High Affinity Binding to Multiple Human Nuclear Receptors. Toxics 2024, 12, 49. [Google Scholar] [CrossRef] [PubMed]
  56. Hahn, D.W.; Allen, G.O.; McGuire, J.L. The Pharmacological Profile of Norgestimate, a New Orally Active Progestin. Contraception 1977, 16, 541–553. [Google Scholar] [CrossRef]
  57. McRobb, L.; Handelsman, D.J.; Kazlauskas, R.; Wilkinson, S.; McLeod, M.D.; Heather, A.K. Structure–Activity Relationships of Synthetic Progestins in a Yeast-Based in Vitro Androgen Bioassay. J. Steroid Biochem. Mol. Biol. 2008, 110, 39–47. [Google Scholar] [CrossRef]
  58. Huang, K.; Whelan, E.A.; Ruder, A.M.; Ward, E.M.; Deddens, J.A.; Davis-King, K.E.; Carreón, T.; Waters, M.A.; Butler, M.A.; Calvert, G.M.; et al. Reproductive Factors and Risk of Glioma in Women. Cancer Epidemiol. Biomark. Prev. 2004, 13, 1583–1588. [Google Scholar] [CrossRef]
  59. McKinley, B.P.; Michalek, A.M.; Fenstermaker, R.A.; Plunkett, R.J. The Impact of Age and Gender on the Incidence of Glial Tumors in New York State from 1976–1995. J. Neurosurg. 2000, 93, 932–939. [Google Scholar] [CrossRef]
  60. Le Rhun, E.; Weller, M. Sex-Specific Aspects of Epidemiology, Molecular Genetics and outcome: Primary Brain Tumours. ESMO Open 2020, 5, e001034. [Google Scholar] [CrossRef] [PubMed]
  61. Skinner, D.C.; Evans, N.P.; Delaleu, B.; Goodman, R.L.; Bouchard, P.; Caraty, A. The Negative Feedback Actions of Progesterone On-Releasing Hormone Secretion Are Transduced by the classical Progesterone Receptor. Proc. Natl. Acad. Sci. USA 1998, 95, 10978–10983. [Google Scholar] [CrossRef] [PubMed]
  62. Altinoz, M.A.; Ozpinar, A.; Elmaci, I. Reproductive Epidemiology of Glial Tumors May Reveal Novel: High-Dose Progestins or Progesterone Antagonists as endocrino-Immune Modifiers against Glioma. Neurosurg. Rev. 2019, 42, 351–369. [Google Scholar] [CrossRef] [PubMed]
  63. Cheng, Y.; Yeung, W.L.; Zhang, P.D.; Li, N.; Kiang, M.Y.; Leung, G.K.K. Progesterone Is More Effective Than Dexamethasone in Prolonging Overall Survival and Preserving Neurologic Function in Experimental Animals with Orthotopic Glioblastoma Allografts. World Neurosurg. 2019, 125, e497–e507. [Google Scholar] [CrossRef]
  64. Altinoz, M.A.; Ucal, Y.; Yilmaz, M.C.; Kiris, İ.; Ozisik, O.; Sezerman, U.; Ozpinar, A.; Elmaci, İ. Progesterone at High Doses Reduces the Growth of U87 and A172 Glioblastoma Cells: Proteomic Changes Regarding and Immunity. Cancer Med. 2020, 9, 5767–5780. [Google Scholar] [CrossRef]
  65. Altinoz, M.A.; Yilmaz, M.; Ucal, Y.; Ozpinar, A.; Baskan, O.; Elmaci, I. P04.54 High Dose Progesterone Induces Growth Inhibition in Human U87 and A172 Glioblastoma Cell with Concomitant Changes in Mitochondrial and Cytoskeleton Proteins. Neuro Oncol. 2018, 20, iii291–iii292. [Google Scholar] [CrossRef]
  66. Germán-Castelán, L.; Manjarrez-Marmolejo, J.; González-Arenas, A.; González-Morán, M.G.; Camacho-Arroyo, I. Progesterone Induces the Growth and Infiltration of Human Astrocytoma Cells Implanted in the Cerebral Cortex of the Rat. BioMed Res. Int. 2014, 2014, 393174. [Google Scholar] [CrossRef] [PubMed]
  67. Cabrera-Muñoz, E.; Hernández-Hernández, O.T.; Camacho-Arroyo, I. Role of Progesterone in Human Astrocytomas Growth. Curr. Top. Med. Chem. 2011, 11, 1663–1667. [Google Scholar] [PubMed]
  68. Hernández-Hernández, O.T.; Camacho-Arroyo, I. Regulation of Gene Expression by Progesterone in Cancer Cells: Effects on Cyclin D1, EGFR and VEGF. Mini Rev. Med. Chem. 2013, 13, 635–642. [Google Scholar] [CrossRef]
  69. Atif, F.; Yousuf, S.; Stein, D.G. Anti-Tumor Effects of Progesterone in Human Glioblastoma Multiforme: Role of PI3K/Akt/MTOR Signaling. J. Steroid Biochem. Mol. Biol. 2015, 146, 62–73. [Google Scholar] [CrossRef]
  70. Atif, F.; Yousuf, S.; Espinosa-Garcia, C.; Sergeeva, E.; Stein, D.G. Progesterone Treatment Attenuates Glycolytic Metabolism and Induces Senescence in Glioblastoma. Sci. Rep. 2019, 9, 988. [Google Scholar] [CrossRef]
  71. Migliaccio, A.; Piccolo, D.; Castoria, G.; Di Domenico, M.; Bilancio, A.; Lombardi, M.; Gong, W.; Beato, M.; Auricchio, F. Activation of the Src/P21ras/Erk Pathway by Progesterone Receptor via Cross–Talk with Estrogen Receptor. EMBO J. 1998, 17, 2008–2018. [Google Scholar] [CrossRef]
  72. Lu, C.-C.; Tsai, S.-C. The Cyclic AMP-Dependent Protein Kinase A Pathway Is Involved in Progesterone Effects on Calcitonin Secretion from TT Cells. Life Sci. 2007, 81, 1411–1420. [Google Scholar] [CrossRef]
  73. Wang, Y.; Hanifi-Moghaddam, P.; Hanekamp, E.E.; Kloosterboer, H.J.; Franken, P.; Veldscholte, J.; van Doorn, H.C.; Ewing, P.C.; Kim, J.J.; Grootegoed, J.A.; et al. Progesterone Inhibition of Wnt/β-Catenin Signaling in Normal Endometrium and Endometrial Cancer. Clin. Cancer Res. 2009, 15, 5784–5793. [Google Scholar] [CrossRef]
  74. Luoma, J.I.; Stern, C.M.; Mermelstein, P.G. Progesterone Inhibition of Neuronal Calcium Signaling Underlies Aspects of Progesterone-Mediated Neuroprotection. J. Steroid Biochem. Mol. Biol. 2012, 131, 30–36. [Google Scholar] [CrossRef]
  75. Schuller, H.M. Beta-Adrenergic Signaling, a Novel Target for Cancer Therapy? Oncotarget 2010, 1, 466–469. [Google Scholar] [CrossRef]
  76. Braadland, P.R.; Ramberg, H.; Grytli, H.H.; Taskén, K.A. Î2-Adrenergic Receptor Signaling in Prostate Cancer. Front. Oncol. 2014, 4, 375. [Google Scholar] [CrossRef]
  77. Chi, T.; Wang, M.; Wang, X.; Yang, K.; Xie, F.; Liao, Z.; Wei, P. PPAR-γ Modulators as Current and Potential Cancer. Front. Oncol. 2021, 11, 737776. [Google Scholar] [CrossRef]
  78. Melbourne, J.K.; Thompson, K.R.; Peng, H.; Nixon, K. Chapter Eight—Its Complicated: The Relationship between Alcohol and Microglia in the Search for Novel Pharmacotherapeutic Targets for Alcohol Use Disorders. In Molecular Basis of Neuropsychiatric Disorders: From Bench to Bedside; Rahman, S., Ed.; Progress in Molecular Biology and Translational Science; Academic Press: Cambridge, MA, USA, 2019; Volume 167, pp. 179–221. [Google Scholar]
  79. Pucci, B.; Kasten, M.; Giordano, A. Cell Cycle and Apoptosis. Neoplasia 2000, 2, 291–299. [Google Scholar] [CrossRef]
  80. Griffin, M.; Khan, R.; Basu, S.; Smith, S. Ion Channels as Therapeutic Targets in High Grade Gliomas. Cancers 2020, 12, 3068. [Google Scholar] [CrossRef]
  81. Seifert, S.; Sontheimer, H. Bradykinin Enhances Invasion of Malignant Glioma into the Brain Parenchyma by Inducing Cells to Undergo Amoeboid Migration. J. Physiol. 2014, 592, 5109–5127. [Google Scholar] [CrossRef]
  82. Yang, K.; Wei, M.; Yang, Z.; Fu, Z.; Xu, R.; Cheng, C.; Chen, X.; Chen, S.; Dammer, E.; Le, W. Activation of Dopamine Receptor D1 Inhibits Glioblastoma Tumorigenicity by Regulating Autophagic Activity. Cell. Oncol. 2020, 43, 1175–1190. [Google Scholar] [CrossRef]
  83. Jeon, H.-M.; Oh, Y.T.; Shin, Y.J.; Chang, N.; Kim, D.; Woo, D.; Yeup, Y.; Joo, K.M.; Jo, H.; Yang, H.; et al. Dopamine Receptor D2 Regulates Glioblastoma Survival and death through MET and Death Receptor 4/5. Neoplasia 2023, 39, 100894. [Google Scholar] [CrossRef] [PubMed]
  84. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinform. 2016, 54, 1–30. [Google Scholar] [CrossRef] [PubMed]
  85. Balakrishnan, R.; Harris, M.A.; Huntley, R.; Van Auken, K.; Cherry, J.M. A Guide to Best Practices for Gene Ontology (GO) Manual Annotation. Database 2013, 2013, bat054. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Network indices to approximate scale-free topology through average number of connections per gene in the networks. Numbers represent the scale free topology model fit per dataset at given soft-thresholding power. The green box indicates the point at which the index across the datasets has stabilized (β = 10).
Figure 1. Network indices to approximate scale-free topology through average number of connections per gene in the networks. Numbers represent the scale free topology model fit per dataset at given soft-thresholding power. The green box indicates the point at which the index across the datasets has stabilized (β = 10).
Biology 13 00206 g001
Figure 2. Approximate linear relationship for the GBM dataset at β = 10 presented by log-log plot.
Figure 2. Approximate linear relationship for the GBM dataset at β = 10 presented by log-log plot.
Biology 13 00206 g002
Figure 3. Dendrogram of gene clustering on TOM-based dissimilarity and module split sensitivities of GBM dataset. The assorted colors represent identified gene co-expression modules corresponding to the portion of the dendrogram directly above.
Figure 3. Dendrogram of gene clustering on TOM-based dissimilarity and module split sensitivities of GBM dataset. The assorted colors represent identified gene co-expression modules corresponding to the portion of the dendrogram directly above.
Biology 13 00206 g003
Figure 4. Module preservation analysis performed on the gene co-expression modules from the GBM network in (a) pilocytic astrocytoma, (b) oligodendroglioma, and (c) anaplastic astrocytoma datasets. The modules in the yellow partition pertain to highly preserved modules (Z > 10); green, yellow, brown, and gray exhibited high preservation in all datasets.
Figure 4. Module preservation analysis performed on the gene co-expression modules from the GBM network in (a) pilocytic astrocytoma, (b) oligodendroglioma, and (c) anaplastic astrocytoma datasets. The modules in the yellow partition pertain to highly preserved modules (Z > 10); green, yellow, brown, and gray exhibited high preservation in all datasets.
Biology 13 00206 g004
Figure 5. Top-enriched terms for the green, yellow, brown, and gray modules in terms of (a) biological processes, (b) cellular component, (c) molecular functions, and (d) KEGG pathways.
Figure 5. Top-enriched terms for the green, yellow, brown, and gray modules in terms of (a) biological processes, (b) cellular component, (c) molecular functions, and (d) KEGG pathways.
Biology 13 00206 g005
Figure 6. The identified top 10 hub gene networks based on the PPI networks of the (a) green, (b) yellow, (c) brown, and (d) gray modules. Nodes represent the hub genes, and edges represent interactions between connected nodes. Colors represent ranking with each network, from red (high) to yellow (low).
Figure 6. The identified top 10 hub gene networks based on the PPI networks of the (a) green, (b) yellow, (c) brown, and (d) gray modules. Nodes represent the hub genes, and edges represent interactions between connected nodes. Colors represent ranking with each network, from red (high) to yellow (low).
Biology 13 00206 g006
Figure 7. KEGG pathways of modules with increasing preservation per glioma grade (I, II, III, and IV) based on module preservation and pathway enrichment analysis results. Signaling pathways within the range of glioma grade are highly preserved (z > 10).
Figure 7. KEGG pathways of modules with increasing preservation per glioma grade (I, II, III, and IV) based on module preservation and pathway enrichment analysis results. Signaling pathways within the range of glioma grade are highly preserved (z > 10).
Biology 13 00206 g007
Figure 8. Cellular processes interconnectivity based on hub gene functions.
Figure 8. Cellular processes interconnectivity based on hub gene functions.
Biology 13 00206 g008
Table 1. Summary of GEO Datasets.
Table 1. Summary of GEO Datasets.
Accession No.GSE50161GSE4290GSE43378GSE36245
ConditionPilocytic Astrocytoma OligodendrogliomaAnaplastic AstrocytomaGlioblastoma Multiforme
TypeExpression Profiling by Array
PlatformGPL570-HG-U133 Plus 2 Affymetrix Human Genome U133 Plus 2.0 Array
SourcePrimary Brain Tumor Samples
No. of Samples14361944
Table 2. Top five drug candidates for the upregulated and downregulated hub genes.
Table 2. Top five drug candidates for the upregulated and downregulated hub genes.
ExpressionGenesDrugMechanismTauFDR
UpregulatedKRAS, CCNB1, BUB1B, KIT, TP53, EGFR, ATM, EXO1, GNAS, CDC20, TOP2A, HSP90AA1, FN1, H3C12, GRIN2B, GRB2, CCNA2, CDK1, CALM1, CALML3, CAMK2A, CREB1, TNF, AKT1, CTNNB1, and ITGAMNorgestimate
Phentolamine
GW0742
Olomoucine
Ambroxol
Progesterone receptor agonist
Adrenergic receptor antagonist
PPAR receptor agonist
CDK inhibitor
Sodium channel blocker
−99.8
−99.7
−99.5
−99.4
−99.3
7.75 × 10−3
6.53 × 10−3
5.05 × 10−4
7.01 × 10−3
8.35 × 10−3
DownregulatedCD4, CACNA1C, CACNG2, CD2, PRKACA, PRKACB, TLR4, and CD8AEthisterone
Noscapine
Nomegestrol
Carmoxirole
Oxcarbazepine
Progestogen hormone
Bradykinin receptor antagonist
Progestogen hormone
Dopamine receptor agonist
Sodium channel blocker
−99.7
−99.6
−99.5
−99.5
−99.5
2.36 × 10−3
6.22 × 10−5
7.89 × 10−4
7.40 × 10−3
9.94 × 10−3
Table 3. Potential drug candidates and their corresponding pathways or processes impeded.
Table 3. Potential drug candidates and their corresponding pathways or processes impeded.
DrugStatusPathway/ProcessReference
NorgestimateApprovedPI3K/Akt Pathway[68,69,70]
PhentolamineApprovedcAMP Signaling, MAPK/ERK, and PI3K/Akt Pathway[75,76]
GW0742ExperimentalLipid Metabolic Pathways[77,78]
OlomoucineApprovedCell Cycle [79]
AmbroxolApprovedCalcium Signaling Pathway[80]
EthisteroneApprovedPI3K/Akt Pathway[68,69,70]
NoscapineApprovedCalcium Signaling Pathway[81]
NomegestrolApprovedPI3K/Akt Pathway[68,69,70]
CarmoxiroleExperimentalcAMP Signaling Pathway[82,83]
OxcarbazepineApprovedCalcium Signaling Pathway[80]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Orda, M.A.; Fowler, P.M.P.T.; Tayo, L.L. Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades. Biology 2024, 13, 206. https://doi.org/10.3390/biology13040206

AMA Style

Orda MA, Fowler PMPT, Tayo LL. Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades. Biology. 2024; 13(4):206. https://doi.org/10.3390/biology13040206

Chicago/Turabian Style

Orda, Marco A., Peter Matthew Paul T. Fowler, and Lemmuel L. Tayo. 2024. "Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades" Biology 13, no. 4: 206. https://doi.org/10.3390/biology13040206

APA Style

Orda, M. A., Fowler, P. M. P. T., & Tayo, L. L. (2024). Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades. Biology, 13(4), 206. https://doi.org/10.3390/biology13040206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop