Shuzhi Ma1,2†
Zhen Guo
Binbin Ji
Yan Wu
Size Chen- 1Department of Oncology, The First Affiliated Hospital of Guangdong Pharmaceutical University, Guangzhou, China
- 2Department of Pathology, Zhujiang Hospital, Southern Medical University, Guangzhou, China
- 3Academician Workstation, Changsha Medical University, Changsha, China
- 4Hunan Key Laboratory of the Research and Development of Novel Pharmaceutical Preparations, Changsha Medical University, Changsha, China
- 5Geneis (Beijing) Co., Ltd., Beijing, China
- 6Guangdong Provincial Engineering Research Center for Esophageal Cancer Precise Therapy, The First Affiliated Hospital of Guangdong Pharmaceutical University, Guangzhou, China
- 7Central Laboratory, The First Affiliated Hospital of Guangdong Pharmaceutical University, Guangzhou, China
Background: Recurrence is still a major obstacle to the successful treatment of gliomas. Understanding the underlying mechanisms of recurrence may help for developing new drugs to combat gliomas recurrence. This study provides a strategy to discover new drugs for recurrent gliomas based on drug perturbation induced gene expression changes.
Methods: The RNA-seq data of 511 low grade gliomas primary tumor samples (LGG-P), 18 low grade gliomas recurrent tumor samples (LGG-R), 155 glioblastoma multiforme primary tumor samples (GBM-P), and 13 glioblastoma multiforme recurrent tumor samples (GBM-R) were downloaded from TCGA database. DESeq2, key driver analysis and weighted gene correlation network analysis (WGCNA) were conducted to identify differentially expressed genes (DEGs), key driver genes and coexpression networks between LGG-P vs LGG-R, GBM-P vs GBM-R pairs. Then, the CREEDS database was used to find potential drugs that could reverse the DEGs and key drivers.
Results: We identified 75 upregulated and 130 downregulated genes between LGG-P and LGG-R samples, which were mainly enriched in human papillomavirus (HPV) infection, PI3K-Akt signaling pathway, Wnt signaling pathway, and ECM-receptor interaction. A total of 262 key driver genes were obtained with frizzled class receptor 8 (FZD8), guanine nucleotide-binding protein subunit gamma-12 (GNG12), and G protein subunit β2 (GNB2) as the top hub genes. By screening the CREEDS database, we got 4 drugs (Paclitaxel, 6-benzyladenine, Erlotinib, Cidofovir) that could downregulate the expression of up-regulated genes and 5 drugs (Fenofibrate, Oxaliplatin, Bilirubin, Nutlins, Valproic acid) that could upregulate the expression of down-regulated genes. These drugs may have a potential in combating recurrence of gliomas.
Conclusion: We proposed a time-saving strategy based on drug perturbation induced gene expression changes to find new drugs that may have a potential to treat recurrent gliomas.
Introduction
Gliomas are the most common type of central nervous system (CNS) tumors, which are composed of various distinct subtype tumors (Galbraith and Snuderl, 2021). Differed from non-CNS neoplasms, the gliomas’ grading system is quite complicate. The latest 2021 WHO Classification of Tumors of the Central Nervous System (WHO CNS5) has integrated certain molecular markers and histological features for more accurate staging of gliomas (Louis et al., 2021). Briefly, the gliomas can be classified as two categories, low-grade gliomas (grade 1–2) with a relatively benign slow-growing feature and favorable prognosis, and high-grade gliomas (grade 3–4) with highly infiltrative ability and malignant form. Glioblastoma multiforme (GBM) accounts for approximately 55% of gliomas and is considered as the most aggressive type of gliomas with a 5-years survival rate less than 5% and a median survival time of 12–15 months (Ostrom et al., 2015a; Ostrom et al., 2015b). Currently, the conventional therapeutic regimen for gliomas is surgical resection followed by radiotherapy and chemotherapy. However, the curative effect is far from satisfaction and recurrence is still the major obstacle to the success of chemoradiotherapy since the majority of GBM would experience recurrence within 6.2 months after diagnosis (Bähr et al., 2009; King and Benhabbour, 2021). Therefore, it is imperative to find new therapeutic target and novel therapeutic strategies for patient with gliomas.
Over the past decades, molecular biomarkers have gained important value in providing diagnostic information and therapeutic target for gliomas. Bevacizumab, an inhibitor of vascular endothelial growth factor (VEGF), was approved to treat recurrent GBM by the Food and Drug Administration (FDA) in March 2009 (Friedman et al., 2009; Kreisl et al., 2009). By targeting VEGF, Bevacizumab inhibits angiogenesis and blocks the nutrient supply, which ultimately impedes the growth and metastasis of GBM. In addition, methylation of the O6 -methylguanine-DNA methyltransferase (MGMT) promoter might serve as a predictive marker for temozolomide (TMZ) treatment response of GBM (Hegi et al., 2008). Poly (ADP ribose) polymerase (PARP) inhibitors can increase tumor sensitivity to TMZ chemotherapy and synergize with radiation therapy (Lesueur et al., 2018). A novel nano-compounds encapsulating wild-type p53 (SGT-53) could enhance the inhibitory effects of TMZ on TMZ-resistant GBM cells (Kim et al., 2015). Moreover, mutations/deregulation in the platelet-derived growth factor receptor alpha (PDGFRα), telomerase reverse transcriptase (TERT), epidermal growth factor receptor (EGFR), c-Myc, phosphatase and tensin homolog (PTEN), serine/threonine-protein kinase (BRAF) are frequently observed in glioma, which have become attractive markers for targeted therapy [Mukasa et al., 2010; Sampson et al., 2010; Killela et al., 2013; Johnson et al., 2014; Liu et al., 2020a; Liu et al., 2020b)].
Though a handful of molecular biomarkers have been discovered, the targeted therapies in clinical trials displayed limited curative effect for gliomas (Wu et al., 2021). This may be attributed to the inter- and intra-heterogeneity in driver mutations and plasticity of gliomas. The recurrent tumor might have a totally distinct gene expression signature in comparison with the primary tumor. Notably, 90% of druggable targets identified at initial diagnosis of gilomas are differentially expressed in a recurrent tumor (Schäfer et al., 2019). Ideally, a patient would select the specific drug according to his/her own molecular genetic feature and change the drugs over the course as the tumor evolves. With the advent of multiomics era, it is becoming possible.
In the era of big data, genome-wide molecular profiling at genome, transcriptome, proteome, and metabolome level have revealed comprehensive landscapes for all major types of gliomas. This has not only enriched our understanding of the molecular mechanism of gliomas pathogenesis and progression, but also broadened our ideas in discovering new therapeutic drugs. As we all know, new drug developing is a time-consuming course with high capital input and low yield, and drug repositioning based on computational tools could largely shorten the process (Liu et al., 2016; Xu et al., 2019; Zhou et al., 2019; Yang et al., 2020a; Liu et al., 2020; Tang et al., 2020; Zhou et al., 2020; Peng et al., 2021). In the present study, we proposed a fast, economical, and comprehensive strategy to find old drugs with new function for combating recurrence of gliomas. Our results on TCGA data suggested that this strategy provides a new direction in discovering drugs and brings hope for people in treating gliomas.
Materials and Methods
Data Collecting and Grouping
The RNA sequencing data of gliomas samples were downloaded from TCGA database. The samples were divided into 4 groups: low grade gliomas primary tumor samples (LGG-P, n = 511), low grade gliomas recurrent tumor samples (LGG-R, n = 18), GBM primary tumor samples (GBM-P, n = 155), and GBM recurrent tumor samples (GBM-R, n = 13).
Screening of DEGs in Gliomas
The R package DESeq2 was applied to identify DEGs in the following data pairs: LGG-P vs LGG-R, GBM-P vs GBM-R. |log2fold change (FC)| ≥ 2, false discovery rate (FDR) < 0.5 and adjusted p value < 0.001 were set as threshold. R package clusterProfiler was used for Gene Ontology (GO) enrichment analysis and calculations. The enriched pathways in the up- or down-regulated gene set were generated using R package ggplot2.
Weighted Gene Correlation Network Analysis and Key Driver Analysis
Network-based methods have been widely used to analyze the associations between various biological entities (Chen et al., 2018; Peng et al., 2018; Peng et al., 2020; Zhang et al., 2021). The R package WGCNA was used to construct a weighted gene co-expression network. The key driver analysis was performed using a software package described by Yang et al.(Yang et al., 2016). The first step was to generate a subnetwork NG, which is located within 2 steps of nodes in a given gene set. Next, the dynamic neighborhood search (DNS) was used to find the gene within 2 steps of each gene in NG. Lastly, by taking the gene set in the first step as the background, the hypergeometric test is carried out to calculate the enrichment value between the gene set in the second step and the input gene set. p value < 0.05 for DEG and p value < 0.01 for subnet were set as threshold in the key driver analysis.
Drug Discovery
Based on the DEGs and key drivers, the CREEDS database was used to find potential drugs. The CREEDS database contains 906 drug perturbation gene expression signatures collected from GEO database (Wang et al., 2016). We screened the drug-gene pairs to identify potential drugs that could reverse the expression of DEGs and key drivers in gliomas. p value < 10–10 was set as the threshold.
Results
A Computational Framework to Identify Biomarkers for Glioma Recurrence and Potential Drugs Targeting Them
We proposed a computational framework biomarker identification and drug discovering for glioma recurrence (Figure 1). Firstly, the differentially expressed genes (DEGs) of primary gliomas samples and recurrent gliomas samples were identified from RNA sequencing data downloaded from The Cancer Genome Atlas (TCGA) database. Secondly, weighted gene correlation network (WGCNA) analysis and key driver analysis were conducted to find co-expression modules and key driver genes. Thirdly, the CREEDS database was applied to find potential drugs that could reverse the DEGs and key drivers. We then applied this framework to the downloaded TCGA data and identified important genes involving in glioma recurrence and potential drugs targeting them.
FIGURE 1. The framework of this study. The RNA-seq data was used to find DEGs, which was followed by GO enrichment, WGCNA, and key driver analysis. Potential drugs that could reverse the DEGs and key drivers were screened through the CREEDS database.
Many DEGs Were Identified Between LGG-P and LGG-R, and Between GBM-P and GBM-R
We conducted a comprehensive analysis of the DEGs between LGG-P and LGG-R, and between GBM-P and GBM-R. The difference between GBM-P and GBM-R samples was not significant as we only obtained 2 upregulated and 29 downregulated genes. A total of 205 DEGs with 75 upregulated and 130 downregulated genes were identified between LGG-P and LGG-R samples. The specific details of each DEGs were shown in Supplementary Table S1. Since the number of DEGs between GBM-P and GBM-R was not large enough, we chose LGG-P and LGG-R pairs for further study. We randomly selected 25 samples to draw the heat map and the top 10 differentially expressed genes were shown in Figure 2.
FIGURE 2. The heat map of the DEGs between LGG-P and LGG-R samples. Dao, mdh1, slc1a3, CNOT4, meox2, pik3c2g, LGR6, ahcyl2, SPARC and LTF were the top 10 differentially upregulated genes; LMNTD2-AS1, MAFA, COL2A1, ONECUT3, MEOX1, BNC1, LMNTD2, GEN1, Wnt11 and trim7 were the top 10 differentially down regulated genes.
GO analysis was utilized to annotate the function of DEGs between LGG-P and LGG-R. The upregulated genes could not be enriched owing to the relatively large p value. The downregulated genes were predominantly enriched in DNA-binding transcription activator activity, extracellular matrix structural constituent, growth factor binding, and platelet-derived growth factor binding for the molecular function (MF) category (Figure 3A). For the cellular component (CC) category, the downregulated genes were correlated with extracellular matrix, collagen-containing extracellular matrix, endoplasmic reticulum lumen, and collagen trimer (Figure 3B). For the biological process (BP) category, the downregulated genes were mainly involved in skeletal system development, extracellular matrix organization, extracellular structure organization, and connective tissue development (Figure 3C). KEGG enrichment analysis was further performed to explore the underling pathological pathways for LGG. As shown in Figure 4, the enrichment pathways include human papillomavirus infection, PI3K-Akt signaling pathway, Wnt signaling pathway, ECM-receptor interaction, and proteoglycans in cancer.
FIGURE 3. GO analysis of the downregulated genes between LGG-P and LGG-R. (A) molecular function category; (B) cellular component category; (C) biological process category. The X-axis is the ratio of differentially expressed genes enriched in the corresponding pathway, and the Y-axis is the name of the pathway.
FIGURE 4. KEGG analysis of the downregulated genes between LGG-P and LGG-R. The DEGs were mainly enriched in protein digestion and absorption, human papillomavirus infection, ECM-receptor interaction, PI3K-Akt signaling pathway, and proteoglycans in cancer pathways.
Coexpression Analysis Revealed Chemical Synaptic Transmission Pathway and T Cell Activation Pathway as Key Modules for Glioma Recurrence
To better understand the function of differentially expressed genes, WGCNA was carried out to identify highly correlated gene clusters. Genes with zero expression were deleted from all samples, and the samples of some separated groups were removed from the hierarchical clustering results. WGCNA finally yielded 150 significant gene modules in LGG-P group and 65 gene modules in LGG-R group. Since there are too many genes to visualize, we randomly selected 400 genes to construct a topological overlapping heat map (Figure 5A) and performed functional enrichment analysis. As shown in Figure 5B and Figure 5C, the highly coexpression genes were mainly enriched in chemical synaptic transmission pathway and T cell activation pathway.
FIGURE 5. Highly correlated gene clusters were identified by WGCNA. (A) Topological overlapping heat map of 400 genes. (B) GO enrichment analysis of module 5. (C) GO enrichment analysis of module 8.
Many Genes Including FZD8 and GNG12 Were Identified as Key Driver Genes for Glioma Recurrence
Key driver gene was considered as the hub gene that connected the up- or downregulated genes. For LGG-P vs LGG-R pair, we obtained 2 key drivers in the upregulated gene set and 260 key drivers in the downregulated gene set. The detailed information of the key drivers can be found in Supplementary Table S2. The most siginificant key drivers and their corresponding subnetwork were shown in Figure 6. We demonstrated that frizzled class receptor 8 (FZD8) and guanine nucleotide-binding protein subunit gamma-12 (GNG12) were two of the hub gene of the downregulated genes. FZD8 is a G protein-coupled receptor protein that plays an important role in β-catenin signaling pathway and regulates cancer invasion and metastasis (Li et al., 2017). GNG12 is a member of the G protein family and participate in a handful of trans-membrane signal transducer pathways (Yuan et al., 2021). G protein subunit β2 (GNB2) is a hub gene of up regulated genes and belongs to the guanine nucleotide-binding proteins family. GNB2 may activate the canonical G protein signaling and involved in cancer initiation and progression (O’hayre et al., 2014). The functions and implications of these hub genes in cancers will be discussed further.
FIGURE 6. A subnetwork of the key drivers that connected the up and down regulated genes. Purple indicates upregulated genes, red indicates key drivers of upregulated genes, green indicates downregulated genes, and blue indicates key drivers of downregulated genes.
A Few Drugs Including Paclitaxel and Fenofibrate Were Identified as Potential Drugs for Preventing Glioma Recurrence
Drugs that have a potential to reverse the expression of DEGs may be valuable for further treatment of gliomas. We obtained 26 drugs that could perturb the expression of up regulated gene sets and 50 drugs that could perturb the expression of down regulated gene sets by CREEDS database. In consideration of drug profile and previous studies, we focused on 4 perturbation drugs of up regulated genes (Paclitaxel, 6-benzyladenine, Erlotinib, Cidofovir) and 5 perturbation drugs of down regulated genes (Fenofibrate, Oxaliplatin, Bilirubin, Nutlins, Valproic acid) (Table 1). These drugs may provide new insight into preventing recurrence of LGG.
TABLE 1. Potential drugs for the treatment of LGG recurrence.
Discussion
Gliomas are highly malignant tumors and recurrence is still the main obstacle to treatment. Though researchers have identified a handful of biomarkers for glioma, the therapeutic effect of targeted drugs is far from satisfaction. Currently, our study provides a fast, economical, and comprehensive method for finding potential drugs to treat gliomas.
Based on RNA sequencing data of LGG-P and LGG-R samples from TCGA, we yielded 75 upregulated and 130 downregulated genes, which were predominantly correlated with human papillomavirus (HPV) infection, PI3K-Akt signaling pathway, Wnt signaling pathway, ECM-receptor interaction. HPV infection is a major cause of cervical cancer, and associated with several epithelial malignancies, including oral cavity, anal, oropharyngeal, penile, vulvar, vaginal, and laryngeal cancers (Lu et al., 2020). Until now, there is no direct evidence indicates that HPV infection is involved in gliomas. We speculate HPV infection may interfere host immune system and participate in LGG recurrence. PI3K-Akt signaling pathway and Wnt signaling pathway are two of the canonical signaling transduction pathways in various cancers. PI3K/Akt pathway controls cell fate by regulating cell growth, apoptosis, angiogenesis, metabolism, autophagy, and chemotherapy resistance of gliomas (Shahcheraghi et al., 2020). Activation of PI3K-Akt pathway is associated with migration and invasion of glioblastoma cells (Huang et al., 2018). Wnt/beta-catenin signaling pathway plays a vital role in ionizing radiation-induced invasion of glioblastoma cells (Dong et al., 2015). It has been reported that high level of beta-catenin was associated with a poor prognosis in glioblastoma patients (Gao et al., 2017). Currently, a number of PI3K inhibitors and wnt inhibitors have entered clinical trials for gliomas treatment, such as BKM120, XL147 and XL765 (Lee et al., 2016; Zhao et al., 2017). ECM-receptor interaction pathway mediates cell migration by regulating neovascularization and diffuse infiltration of tumor cells (Cui et al., 2018). Previous study also indicated ECM-receptor interaction pathway was abnormal in development and survival of glioblastoma (Bo et al., 2015; Yang et al., 2020b).
We also identified several key driver genes that contributed to LGG recurrence, including FZD8, GNG12, GNB2. FZD8 could activate the β-catenin pathway and play a vital role in cancer invasion and metastasis (Chen et al., 2020). Aberrant expression of FZD8 has been reported in gastric cancer, prostate cancer, renal cell carcinoma, lung cancer, pancreatic adenocarcinoma, and overexpression of FZD8 was considered to promote tumor metastasis (Li et al., 2017; Yang et al., 2017; Liu et al., 2019; Chen et al., 2020; Li et al., 2021). In addition, overexpression of FZD8 leaded to chemotherapy resistance in breast cancer patients (Yin et al., 2013). GNG12 acted as an important modulator or transducer in various transmembrane signaling systems. Researchers have demonstrated that GNG12 could regulate cancer cell proliferation, inflammatory response, and immune response via activating the mTORC1 pathway and NF-κB signaling pathway (Larson et al., 2010; Luo et al., 2018; Li et al., 2020). GNB2 was involved in cancer initiation and progression by activating AKT/mTOR pathway, MAPK pathway, and Hippo signaling pathway. Mutations of GNB2 may result in targeted kinase inhibitors resistance to numerous types of cancer (Yoda et al., 2015). The roles of FZD8, GNG12, GNB2 have not been fully illustrated in gliomas and needs further investigation. These key driver genes may help for understanding the pathogenesis of for LGG recurrence and shed new insight for developing new drugs.
However, developing new drugs from a molecular biomarker is a great project and still has a long way to go. In the present study, a total of 9 drugs with potential therapeutic effect against LGG recurrence were selected through a drug-gene perturbation method. Paclitaxel is a natural anticancer drug that has been widely used in the therapy of breast cancer, ovarian cancer, lung cancer, and several head and neck cancers. Paclitaxel binds to tubulin, promotes its assembly with microtubules and inhibits dissociation, which finally prevents mitosis and hinders cell cycle progression (Zhu and Chen, 2019). Erlotinib is used to treat some types of lung cancer and advanced or metastatic pancreatic cancer in clinical. Erlotinib blocks EGFR pathway by inhibiting tyrosine kinase activity and impeding cell proliferation, apoptosis, angiogenesis, invasions, and metastasis. Benzylaminopurine is a first generation cytokinin that stimulates cell division and inhibits respiratory kinase, leading to plant growth and development. Cidofovir is used to treat cytomegalovirus (CMV) infection through inhibiting viral DNA polymerase (Alcamo et al., 2020). Cidofovir also has a strong activity against herpes simplex virus (HSV), varicella zoster virus (VZV), adenovirus (AV), and human papillomavirus (HPV). Fenofibrate is widely used as a lipid-lowering drug through activating PPARα-RXR signal and transcription of lipid metabolism related genes. Numerous evidence has indicated that fenofibrate might exert anticancer effects through regulating cell apoptosis, cell-cycle arrest, invasion, and migration (Lian et al., 2018). Oxaliplatin is a third-generation platinum analog that has been widely used as the first-line drugs for metastatic colorectal cancer. It blocks DNA replication by binding to DNA and forming cross-linked DNA adducts, which consequently leading to cancer cell death (Mauri et al., 2020). Bilirubin is an endogenous metabolite from haem. Recent studies have indicated bilirubin levels may serve as biomarker for several cancers and vascular disease (Horsfall et al., 2020; Seyed Khoei et al., 2020). Nutlins is a small molecule that could displace p53 protein from p53/MDM2 complex, thereby preventing the degradation of p53. It has been revealed that Nutlins could induce p53 dependent cell cycle arrest and apoptosis in a number of tumors (Impicciatore et al., 2010). Valproic acid (VPA), a histone deacetylase (HDAC) inhibitor, is widely used to treat epilepsia, bipolar disorders, migraine, and schizophrenia. In addition, VPA may exert anti-tumor activity by regulating cell proliferation, apoptosis, differentiation, adhesion, invasion, migration, angiogenesis, and inflammation (Michaelis et al., 2007). As shown above, some candidate drugs were already commonly used in clinical, some candidates only showed preclinical antitumor activity. Weather these drugs/agents could prevent LGG recurrence still needs further preclinical and clinical trial validation.
In conclusion, we conducted a comprehensive analysis of LGG-P and LGG-R samples to find DEGs and key driver genes. By using a drug-gene perturbation method, a serious potential drugs/agents were screened to treat LGG recurrence. However, the exact effect of these drugs on glioma recurrence needs further experimental data for verification. Besides, independent dataset with paired primary tumor and recurrent tumor samples is needed to validate the findings in the future. This study may broaden our understanding of the molecular mechanism of LGG recurrence and provide new sights for drug discovery.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
SC contributed to conception and design of the study. BW, MY, and XY organized the database. SM, ZG, and BW performed the statistical analysis. SM and ZG wrote the first draft of the manuscript. MY, XY, BJ, and YW wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was supported by the Natural Science Foundation of China (No. 82104309), Natural Science Foundation of Hunan Province (No.2021JJ40639), Outstanding Youth Project of Hunan Provincial Education Department (No.20B075). Science and Technology Program of Guangzhou, China (Grant numbers: 201904010047); Science and Technology Planning Project of Guangdong Province of China (Grant numbers: 2020A0505100058); Guangdong Educational Committee (Key Project of Regular institutions of higher learning of Guangdong Province (Grant numbers: 2019KZDXM024).
Conflict of Interest
BW, MY, XY, BJ, and YW were employed by the company Geneis (Beijing) Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, orclaim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.832627/full#supplementary-material
Supplementary Table S1 | DEGs discovered in LGG-P and LGG-R samples.
Supplementary Table S2 | Key drivers discovered in LGG-P and LGG-R samples.
References
Alcamo, A. M., Wolf, M. S., Alessi, L. J., Chong, H. J., Green, M., Williams, J. V., et al. (2020). Successful Use of Cidofovir in an Immunocompetent Child with Severe Adenoviral Sepsis. Pediatrics 145. doi:10.1542/peds.2019-1632
Bähr, O., Herrlinger, U., Weller, M., and Steinbach, J. P. (2009). Very Late Relapses in Glioblastoma Long-Term Survivors. J. Neurol. 256, 1756–1758. doi:10.1007/s00415-009-5167-6
Bo, L. J., Wei, B., Li, Z. H., Wang, Z. F., Gao, Z., and Miao, Z. (2015). Bioinformatics Analysis of miRNA Expression Profile between Primary and Recurrent Glioblastoma. Eur. Rev. Med. Pharmacol. Sci. 19, 3579–3586.
Chen, W., Liu, Z., Mai, W., Xiao, Y., You, X., and Qin, L. (2020). FZD8 Indicates a Poor Prognosis and Promotes Gastric Cancer Invasion and Metastasis via B-Catenin Signaling Pathway. Ann. Clin. Lab. Sci. 50, 13–23.
Chen, X., Yin, J., Qu, J., and Huang, L. (2018). MDHGI: Matrix Decomposition and Heterogeneous Graph Inference for miRNA-Disease Association Prediction. Plos Comput. Biol. 14 (8), e1006418. doi:10.1371/journal.pcbi.1006418
Cui, X., Morales, R.-T. T., Qian, W., Wang, H., Gagner, J.-P., Dolgalev, I., et al. (2018). Hacking Macrophage-Associated Immunosuppression for Regulating Glioblastoma Angiogenesis. Biomaterials 161, 164–178. doi:10.1016/j.biomaterials.2018.01.053
Dong, Z., Zhou, L., Han, N., Zhang, M., and Lyu, X. (2015). Wnt/β-catenin Pathway Involvement in Ionizing Radiation-Induced Invasion of U87 Glioblastoma Cells. Strahlenther Onkol 191, 672–680. doi:10.1007/s00066-015-0858-7
Friedman, H. S., Prados, M. D., Wen, P. Y., Mikkelsen, T., Schiff, D., Abrey, L. E., et al. (2009). Bevacizumab Alone and in Combination with Irinotecan in Recurrent Glioblastoma. Jco 27, 4733–4740. doi:10.1200/jco.2008.19.8721
Galbraith, K., and Snuderl, M. (2021). Molecular Pathology of Gliomas. Surg. Pathol. Clin. 14, 379–386. doi:10.1016/j.path.2021.05.003
Gao, L., Chen, B., Li, J., Yang, F., Cen, X., Liao, Z., et al. (2017). Wnt/β-catenin Signaling Pathway Inhibits the Proliferation and Apoptosis of U87 Glioma Cells via Different Mechanisms. PLoS One 12, e0181346. doi:10.1371/journal.pone.0181346
Hegi, M. E., Liu, L., Herman, J. G., Stupp, R., Wick, W., Weller, M., et al. (2008). Correlation of O6-Methylguanine Methyltransferase (MGMT) Promoter Methylation with Clinical Outcomes in Glioblastoma and Clinical Strategies to Modulate MGMT Activity. Jco 26, 4189–4199. doi:10.1200/jco.2007.11.5964
Horsfall, L. J., Burgess, S., Hall, I., and Nazareth, I. (2020). Genetically Raised Serum Bilirubin Levels and Lung Cancer: a Cohort Study and Mendelian Randomisation Using UK Biobank. Thorax 75, 955–964. doi:10.1136/thoraxjnl-2020-214756
Huang, W., Ding, X., Ye, H., Wang, J., Shao, J., and Huang, T. (2018). Hypoxia Enhances the Migration and Invasion of Human Glioblastoma U87 Cells through PI3K/Akt/mTOR/HIF-1α Pathway. Neuroreport 29, 1578–1585. doi:10.1097/wnr.0000000000001156
Impicciatore, G., Sancilio, S., Miscia, S., and Di Pietro, R. (2010). Nutlins and Ionizing Radiation in Cancer Therapy. Cpd 16, 1427–1442. doi:10.2174/138161210791033932
Johnson, B. E., Mazor, T., Hong, C., Barnes, M., Aihara, K., Mclean, C. Y., et al. (2014). Mutational Analysis Reveals the Origin and Therapy-Driven Evolution of Recurrent Glioma. Science 343, 189–193. doi:10.1126/science.1239947
Killela, P. J., Reitman, Z. J., Jiao, Y., Bettegowda, C., Agrawal, N., Diaz, L. A., et al. (2013). TERT Promoter Mutations Occur Frequently in Gliomas and a Subset of Tumors Derived from Cells with Low Rates of Self-Renewal. Proc. Natl. Acad. Sci. 110, 6021–6026. doi:10.1073/pnas.1303607110
Kim, S.-S., Rait, A., Kim, E., Pirollo, K. F., and Chang, E. H. (2015). A Tumor-Targeting P53 Nanodelivery System Limits Chemoresistance to Temozolomide Prolonging Survival in a Mouse Model of Glioblastoma Multiforme. Nanomedicine: Nanotechnology, Biol. Med. 11, 301–311. doi:10.1016/j.nano.2014.09.005
King, J. L., and Benhabbour, S. R. (2021). Glioblastoma Multiforme-A Look at the Past and a Glance at the Future. Pharmaceutics 13. doi:10.3390/pharmaceutics13071053
Kreisl, T. N., Kim, L., Moore, K., Duic, P., Royce, C., Stroud, I., et al. (2009). Phase II Trial of Single-Agent Bevacizumab Followed by Bevacizumab Plus Irinotecan at Tumor Progression in Recurrent Glioblastoma. Jco 27, 740–745. doi:10.1200/jco.2008.16.3055
Larson, K. C., Draper, M. P., Lipko, M., and Dabrowski, M. (2010). Gng12 Is a Novel Negative Regulator of LPS-Induced Inflammation in the Microglial Cell Line BV-2. Inflamm. Res. 59, 15–22. doi:10.1007/s00011-009-0062-2
Lee, Y., Lee, J.-K., Ahn, S. H., Lee, J., and Nam, D.-H. (2016). WNT Signaling in Glioblastoma and Therapeutic Opportunities. Lab. Invest. 96, 137–150. doi:10.1038/labinvest.2015.140
Lesueur, P., Chevalier, F., El-Habr, E. A., Junier, M.-P., Chneiweiss, H., Castera, L., et al. (2018). Radiosensitization Effect of Talazoparib, a Parp Inhibitor, on Glioblastoma Stem Cells Exposed to Low and High Linear Energy Transfer Radiation. Sci. Rep. 8, 3664. doi:10.1038/s41598-018-22022-4
Li, J., Jin, C., Zou, C., Qiao, X., Ma, P., Hu, D., et al. (2020). GNG12 Regulates PD‐L1 Expression by Activating NF‐κB Signaling in Pancreatic Ductal Adenocarcinoma. FEBS Open Bio 10, 278–287. doi:10.1002/2211-5463.12784
Li, Q., Ye, L., Zhang, X., Wang, M., Lin, C., Huang, S., et al. (2017). FZD8, a Target of P53, Promotes Bone Metastasis in Prostate Cancer by Activating Canonical Wnt/β-Catenin Signaling. Cancer Lett. 402, 166–176. doi:10.1016/j.canlet.2017.05.029
Li, Y., Liu, Z., and Zhang, Y. (2021). Expression and Prognostic Impact of FZDs in Pancreatic Adenocarcinoma. BMC Gastroenterol. 21, 79. doi:10.1186/s12876-021-01643-6
Lian, X., Wang, G., Zhou, H., Zheng, Z., Fu, Y., and Cai, L. (2018). Anticancer Properties of Fenofibrate: A Repurposing Use. J. Cancer 9, 1527–1537. doi:10.7150/jca.24488
Liu, C., Wei, D., Xiang, J., Ren, F., Huang, L., Lang, J., et al. (2020). An Improved Anticancer Drug-Response Prediction Based on an Ensemble Method Integrating Matrix Completion and Ridge Regression. Mol. Ther. - Nucleic Acids 21, 676–686. doi:10.1016/j.omtn.2020.07.003
Liu, F., Peng, L., Tian, G., Yang, J., Chen, H., Hu, Q., et al. (2020a). Identifying Small Molecule-miRNA Associations Based on Credible Negative Sample Selection and Random Walk. Front. Bioeng. Biotechnol. 8, 131. doi:10.3389/fbioe.2020.00131
Liu, R., Chen, Y., Shou, T., Hu, J., and Qing, C. (2019). miRNA-99b-5p Targets FZD8 to Inhibit Non-small Cell Lung Cancer Proliferation, Migration and Invasion. Ott 12, 2615–2621. doi:10.2147/ott.s199196
Liu, X., Lang, J., Li, S., Wang, Y., Peng, L., Wang, W., et al. (2020b). Fragment Enrichment of Circulating Tumor DNA with Low-Frequency Mutations. Front. Genet. 11, 147. doi:10.3389/fgene.2020.00147
Liu, X., Yang, J., Zhang, Y., Fang, Y., Wang, F., Wang, J., et al. (2016). A Systematic Study on Drug-Response Associated Genes Using Baseline Gene Expressions of the Cancer Cell Line Encyclopedia. Sci. Rep. 6, 22811. doi:10.1038/srep22811
Louis, D. N., Perry, A., Wesseling, P., Brat, D. J., Cree, I. A., Figarella-Branger, D., et al. (2021). The 2021 WHO Classification of Tumors of the Central Nervous System: a Summary. Neuro Oncol. 23, 1231–1251. doi:10.1093/neuonc/noab106
Lu, Y., Li, P., Luo, G., Liu, D., and Zou, H. (2020). Cancer Attributable to Human Papillomavirus Infection in China: Burden and Trends. Cancer 126, 3719–3732. doi:10.1002/cncr.32986
Luo, C., Zhao, S., Dai, W., Zheng, N., and Wang, J. (2018). Proteomic Analyses Reveal GNG12 Regulates Cell Growth and Casein Synthesis by Activating the Leu-Mediated mTORC1 Signaling Pathway. Biochim. Biophys. Acta (Bba) - Proteins Proteomics 1866, 1092–1101. doi:10.1016/j.bbapap.2018.08.013
Mauri, G., Gori, V., Bonazzina, E., Amatu, A., Tosi, F., Bencardino, K., et al. (2020). Oxaliplatin Retreatment in Metastatic Colorectal Cancer: Systematic Review and Future Research Opportunities. Cancer Treat. Rev. 91, 102112. doi:10.1016/j.ctrv.2020.102112
Michaelis, M., Doerr, H., and Cinatl Jr., J. (2007). Valproic Acid as Anti-cancer Drug. Cpd 13, 3378–3393. doi:10.2174/138161207782360528
Mukasa, A., Wykosky, J., Ligon, K. L., Chin, L., Cavenee, W. K., and Furnari, F. (2010). Mutant EGFR Is Required for Maintenance of Glioma Growth In Vivo, and its Ablation Leads to Escape from Receptor Dependence. Proc. Natl. Acad. Sci. 107, 2616–2621. doi:10.1073/pnas.0914356107
O'hayre, M., Degese, M. S., and Gutkind, J. S. (2014). Novel Insights into G Protein and G Protein-Coupled Receptor Signaling in Cancer. Curr. Opin. Cel Biol 27, 126–135. doi:10.1016/j.ceb.2014.01.005
Ostrom, Q. T., Gittleman, H., Fulop, J., Liu, M., Blanda, R., Kromer, C., et al. (2015a). CBTRUS Statistical Report: Primary Brain and Central Nervous System Tumors Diagnosed in the United States in 2008-2012. Neuro Oncol. 17 Suppl 4 (Suppl. 4), iv1. iv1-iv62. doi:10.1093/neuonc/nov189
Ostrom, Q. T., Gittleman, H., Stetson, L., Virk, S. M., and Barnholtz-Sloan, J. S. (2015b). Epidemiology of Gliomas. Cancer Treat. Res. 163, 1–14. doi:10.1007/978-3-319-12048-5_1
Peng, L.-H., Sun, C.-N., Guan, N.-N., Li, J.-Q., and Chen, X. (20182018). HNMDA: Heterogeneous Network-Based miRNA-Disease Association Prediction. Mol. Genet. Genomics 293 (4), 983–995. doi:10.1007/s00438-018-1438-1
Peng, L. H., Shen, L., Tian, X. F., Liu, F. X., Wang, J. J., Tian, G., et al. (20212010). Prioritizing Antiviral Drugs against SARS-CoV-2 by Integrating Viral Complete Genome Sequences and Drug Chemical Structures. Scientific reportsImmunologic Escape after Prolonged Progression-free Survival with Epidermal Growth Factor Receptor Variant III Peptide Vaccination in Patients with Newly Diagnosed Glioblastoma. J. Clin. Oncol. 1128 (1), 14722–114729.
Peng, L. H., Tian, X. F., Shen, L., Kuang, M., Tian, G., Yang, J. L., et al. (2020). Identifying Effective Antiviral Drugs against SARS-CoV-2 by Drug Repositioning through Virus-Drug Association Prediction[J]. Front. Genet. 11, 1072. doi:10.3389/fgene.2020.577387
Schäfer, N., Gielen, G. H., Rauschenbach, L., Kebir, S., Till, A., Reinartz, R., et al. (2019). Longitudinal Heterogeneity in Glioblastoma: Moving Targets in Recurrent versus Primary Tumors. J. Transl Med. 17, 96. doi:10.1186/s12967-019-1846-y
Seyed Khoei, N., Jenab, M., Murphy, N., Banbury, B. L., Carreras-Torres, R., Viallon, V., et al. (2020). Circulating Bilirubin Levels and Risk of Colorectal Cancer: Serological and Mendelian Randomization Analyses. BMC Med. 18, 229. doi:10.1186/s12916-020-01703-w
Shahcheraghi, S. H., Tchokonte-Nana, V., Lotfi, M., Lotfi, M., Ghorbani, A., and Sadeghnia, H. R. (2020). Wnt/beta-catenin and PI3K/Akt/mTOR Signaling Pathways in Glioblastoma: Two Main Targets for Drug Design: A Review. Cpd 26, 1729–1741. doi:10.2174/1381612826666200131100630
Tang, X., Cai, L., Meng, Y., Xu, J., Lu, C., and Yang, J. (2020). Indicator Regularized Non-negative Matrix Factorization Method-Based Drug Repurposing for COVID-19. Front. Immunol. 11, 603615. doi:10.3389/fimmu.2020.603615
Wang, Z., Monteiro, C. D., Jagodnik, K. M., Fernandez, N. F., Gundersen, G. W., Rouillard, A. D., et al. (2016). Extraction and Analysis of Signatures from the Gene Expression Omnibus by the Crowd. Nat. Commun. 7, 12846. doi:10.1038/ncomms12846
Wu, W., Klockow, J. L., Zhang, M., Lafortune, F., Chang, E., Jin, L., et al. (2021). Glioblastoma Multiforme (GBM): An Overview of Current Therapies and Mechanisms of Resistance. Pharmacol. Res. 171, 105780. doi:10.1016/j.phrs.2021.105780
Xu, X., Long, H., Xi, B., Ji, B., Li, Z., Dang, Y., et al. (2019). Molecular Network-Based Drug Prediction in Thyroid Cancer. Int. J. Mol. Sci. 20. doi:10.3390/ijms20020263
Yang, J. a., Wang, L., Xu, Z., Wu, L., Liu, B., Wang, J., et al. (2020b). Integrated Analysis to Evaluate the Prognostic Value of Signature mRNAs in Glioblastoma Multiforme. Front. Genet. 11, 253. doi:10.3389/fgene.2020.00253
Yang, J., Huang, T., Song, W.-m., Petralia, F., Mobbs, C. V., Zhang, B., et al. (2016). Discover the Network Mechanisms Underlying the Connections between Aging and Age-Related Diseases. Sci. Rep. 6, 32566. doi:10.1038/srep32566
Yang, J., Peng, S., Zhang, B., Houten, S., Schadt, E., Zhu, J., et al. (2020a). Human Geroprotector Discovery by Targeting the Converging Subnetworks of Aging and Age-Related Diseases. Geroscience 42, 353–372. doi:10.1007/s11357-019-00106-x
Yang, Q., Wang, Y., Pan, X., Ye, J., Gan, S., Qu, F., et al. (2017). Frizzled 8 Promotes the Cell Proliferation and Metastasis of Renal Cell Carcinoma. Oncotarget 8, 78989–79002. doi:10.18632/oncotarget.20742
Yin, S., Xu, L., Bonfil, R. D., Banerjee, S., Sarkar, F. H., Sethi, S., et al. (2013). Tumor-initiating Cells and FZD8 Play a Major Role in Drug Resistance in Triple-Negative Breast Cancer. Mol. Cancer Ther. 12, 491–498. doi:10.1158/1535-7163.mct-12-1090
Yoda, A., Adelmant, G., Tamburini, J., Chapuy, B., Shindoh, N., Yoda, Y., et al. (2015). Mutations in G Protein β Subunits Promote Transformation and Kinase Inhibitor Resistance. Nat. Med. 21, 71–75. doi:10.1038/nm.3751
Yuan, J., Yuan, Z., Ye, A., Wu, T., Jia, J., Guo, J., et al. (2021). Low GNG12 Expression Predicts Adverse Outcomes: A Potential Therapeutic Target for Osteosarcoma. Front. Immunol. 12, 758845. doi:10.3389/fimmu.2021.758845
Zhang, L., Yang, P., Feng, H., Zhao, Q., and Liu, H. (20212021). Using Network Distance Analysis to Predict lncRNA-miRNA Interactions. Interdiscip. Sci. Comput. Life Sci. 13 (3), 535–545. doi:10.1007/s12539-021-00458-z
Zhao, H.-f., Wang, J., Shao, W., Wu, C.-p., Chen, Z.-p., To, S.-s. T., et al. (2017). Recent Advances in the Use of PI3K Inhibitors for Glioblastoma Multiforme: Current Preclinical and Clinical Development. Mol. Cancer 16, 100. doi:10.1186/s12943-017-0670-3
Zhou, L., Li, Z., Yang, J., Tian, G., Liu, F., Wen, H., et al. (2019). Revealing Drug-Target Interactions with Computational Models and Algorithms. Molecules 24 (9), 1714. doi:10.3390/molecules24091714
Zhou, L., Wang, J., Liu, G., Lu, Q., Dong, R., Tian, G., et al. (2020). Probing Antiviral Drugs against SARS-CoV-2 through Virus-Drug Association Prediction Based on the KATZ Method. Genomics 112 (6), 4427–4434. doi:10.1016/j.ygeno.2020.07.044
Keywords: low grade gliomas, RNA-seq, differentially expressed genes, WGCNA, key driver genes, drug discovery
Citation: Ma S, Guo Z, Wang B, Yang M, Yuan X, Ji B, Wu Y and Chen S (2022) A Computational Framework to Identify Biomarkers for Glioma Recurrence and Potential Drugs Targeting Them. Front. Genet. 12:832627. doi: 10.3389/fgene.2021.832627
Received: 10 December 2021; Accepted: 29 December 2021;
Published: 17 January 2022.
Edited by:
Lihong Peng, Hunan University of Technology, ChinaCopyright © 2022 Ma, Guo, Wang, Yang, Yuan, Ji, Wu and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Size Chen, Y2hlbnNpemVAZ2RwdS5lZHUuY24=
†These authors have contributed equally to this work