Purine and Pyrimidine-Linked Enzymes and Genes are Strongly Responsible for the Development of Tumors, Particularly Glioblastoma Multiforme

We were aiming to delineate, by the utility of the biological data results, in our investigations the link between the purine and pyrimidine metabolism and development of the glioblastoma. We analyzed the sets of the genes, belonging to the purine and pyrimidine metabolism by the utility of GSEA software as well as MSIgnDB application of the GSEA. The GEO database, GEOR2 tools were serving for the visualization of the genes expression profiles of the disease. The Cancer Proteome Atlas as well as the tools of the data sets were also used to collect and analyze the results. We concluded and came to the following consequential results. 1) Neurogenesis and Glioblastoma are sharing some common genes. 2) Purine and pyrimidine metabolism-linked enzymes and genes are responsible for the upregulation of DNA and mRNA synthesis in the settings of the tumor development. 3) EGFR expression responsible genes, mRNA as well as protein is upregulated during the development of the glioblastoma. 4) GMPS genes are more strongly upregulated in the settings of the glioblastoma than ADSL. 5) PRPS1 is strongly synthetized in neurospheres in contrast to the mature tissue during glioblastoma development.


Introduction
We were aiming to delineate in our investigations, by the utility of the biological data results, the link between the purine and pyrimidine metabolism and development of glioblastoma.
In our previous experiments, we have proved experimentally, by the inhibition of the Xanthine Oxidoreductase (EC 1. 17.3.2), the one of the last and the key regulative enzymes in the chain of the purines catabolism-responsible enzymes, by its classical compound-Allopurinol, it is possible to decrease the speed of purines decay [1] [2]. In addition, we have proved, by impacting the XOR, it is possible to accelerate the cells proliferation or have an influence on the cells' deaths [2] [3] [4].
During our recent experiments we injected primer human glioblastoma cells (island-like spheres groups of the glioblastoma cells culture) into the rats' temporal cortex. At day 21 st AdenosineDeaminase2 (ADA2; EC 3.5.4.4), PRPS1 as well as the XOR activities were measured. Interestingly, activities of the mentioned enzymes were elevated [6].
The purpose of our current study was the globalized analyses of the involvement of the purine and pyrimidine metabolism responsible compounds and enzymes in the processes of the glioblastoma development.

Methods
We analyzed the sets of the genes, belonging to the purine and pyrimidine metabolism by the utility of GSEA software (on-line version, https://www.gsea-msigdb.org/gsea/index.jsp), Molecular Signature Databases (MSIgnDB) application of the GSEA (https://www.gsea-msigdb.org/gsea/msigdb). The GEO database, GEOR2 tools were serving for the visualization of the genes expression profiles of the disease. The Cancer Proteome Atlas as well as the tools of the data sets also were used to collect and analyze the results, including glioblastoma elevated and specific gene profile, mRNA and proteins (https://tcpaportal.org/tcpa/).
The graphs and tables, built by the programs, were replaced with the manually drawn to increase the resolution.

Results
1) First we extracted the genes from data sets, responsible for the purines' and The existence of the numerous cancer-linked genes of the purine and pyrimidine catabolism was providing the fact of strong correlation between metabolic mentioned processes and cancer development. In addition, we drew out the KEGG pathways and sets of the genes of KEGG pathways, which were involving that particular genes as the cancer modules (Table 1). It is necessary to mention, we withdrew from GSEA 98 genes for pyrimidine metabolism, 159-from purine.
After summarizing the purine and pyrimidine genes and collection them in the MSIgnDB window we applied cancer related "filters". Analyzing the results after application of the filters, we concluded, 91 genes from 257 were involved in the processes of cancer development.
We concluded, the predominant parts of those pathways are related to the RNA, DNA syntheses as well as the DNA repair in the settings of the cancer development.
2) Then we isolated the common genes (RRM1, POLR2I, POLE3, POLA2) from the purine metabolism vs gene neighborhoods network and pyrimidine metabolism vs gene neighborhoods and found the KEGG pathways, where these genes were involved, in the settings of the cancer development (Table 2). It's logical and the results are evidencing that in the cancer-linked environment these genes are predominantly responsible for DNA replication.
3) Then we found the specific genes up which were regulated only during the Glioblastoma from the Human protein atlas.
These genes were transferred into the GSEA software for the further analyses. We found, these genes were involved into the data sets of the other sets of the genes, responsible mostly for several biological processes such as formation of the synapses and neurogenesis, which are evidencing about the common mechanism between the regenerative neurogenesis and glioblastoma development. Some sets of glioblastoma development genes are responsible for the establishment of the cell-cell signaling and formation of the synapsis, regulation of the nervous system development, ions transmembrane transporting and, finally, neurogenesis.
Interestingly, in our previous publications, we proved, the activation of the PRPS1 in the settings of experimental stroke modeling might increase the number of the Ki-67 positive cells in the region of the hippocampus [7]. 4) Further, we analyzed the mRNA existence profile from the Cancer Proteome Atlas in the settings of glioblastoma. We found 6 mRNA. MRNA and genes from that particular subpart of the data sets which are Annexin 1, Caveolin 1, CKIT, EGFR, EGFR_pY1068, P16INK4A. These are mostly the structural and EGF related proteins.
It is very well know, EgfR (Egf Receptor) activation initiates formation of the tumors, different types of it, moreover, has a crucial role in progression of the metastasis [8]. EgfR are responsible for the triggering of the following pathways: Ras/Raf signaling cascade, which promotes cell survival and cell proliferation; phosphatidylinositol 3-kinase/Akt signaling cascade, responsible for cell growth, apoptosis resistance, cell invasion and migration; signal transducers and activators of transcription (STAT) pathway, inducing tumor progression, oncogenesis and angiogenesis; phospholipase Cγ signaling regulation of ion channels, which induces cell migration, calcium-mediated signaling; Nck/PAK signaling cascade cell, responsible for the survival and cell migration; Cbl mediated endocytosis, inducing respectfully endocytotic processes [9]- [14].
Interestingly, EGF and EGFR activates the MAP kinase pathways, which in its own turn, by the process of phosphorylation activates the CAD (carbamoylphosphatesynthetase-aspartate carbamoyltransferase-dihydroorotase). It is the multifunctional protein, triggering the mammalian pyrimidine biosynthesis [15] [16]. Zollner et al. proved, by influencing on the pyrimidine metabolism, it is possible to have an impact on the purine metabolic pathways [17]. 5) We worked also with the GEO data set and analyzed the profile of expression of the PRPS1 enzymes responsible for the syntheses of purines and pyrimidines.
We have chosen analysis of primary neurosphere cultures of glioblastoma (GBM) or oligodendroglial (ODG) tumors from 6 patients study (GDS4473).
Patient-derived glioma-propagating cells (GPC) contain karyotypic and gene expression profiles that are found in the primary tumor. From the bottom of the page of the particular study (GDS4473) we chose GEOR tools. By "define" window the comparison groups were divided into the two folders of genes "oligodendroglialglioma" vs "neurospheres" in Homo sapiens. The "Value Distribution" window allowed to compare the quantity related homogeneity of the groups to be compared. "Top 250" allowed the comparison of the genes. Any genes distribution in the groups might be compared from the list of the entire sets of the genes.
Two sets of the results were analyzed regarding PRPS1 gene expression in glioblastoma multiforme and oligodendroglialglioma vsneurospheres in Homo sapiens. The GEO data results are evidencing about the up regulation of PRPS1 in neurospheres fraction in comparison with the already formed and maturated glioblastoma tissue (Figure 1).
By the utility of GEO2R we came to the conclusion, another two key regulative enzymes of purine and pyrimidine metabolism were up regulated in the glioblastoma multiforme tissues in comparison with the normal brain fractions.
Interestingly, from those two enzymes GMPS was more prominently unregulated in comparison with the ADSL (Figure 2).

Discussion
Glioblastoma (GBM) is the most common primary malignant brain tumor, comprising 16% of all primary brain and central nervous system neoplasms [18]. The average age-adjusted incidence rate is 3.2 per 100,000 population Originally, GBMs were thought to be derived solely from glial cells; however, evidence suggests  that they may arise from multiple cell types with neural stem cell-like properties.
These cells are at multiple stages of differentiation from stem cell to neuron, to glia, with phenotypic variations determined, in large part, by molecular alterations in signaling pathways rather than by differences in the cell type of origin [19]. They might differentiate into the normal or atypical cells.

Conclusions
We concluded and came to the following consequential results.
1) Neurogenesis and Glioblastoma are sharing some common genes. 5) PRPS1 is strongly synthetized in neurospheres than in the mature tissue during glioblastoma development.