RESUMO
Primary glioblastoma(pGBM) is the most malignant tumor of the central nervous system. Radiotherapy, chemotherapy and surgical treatment have little effect on the survival of pGBM patients. The prognosis is often poorly once the tumor recurs. It is urgent to develop new therapies for patients. In recent years, studies have been clarified that miRNA have a powerful regulating effect on the genes. However, the main group of miRNAs in regulating long-term survival specific related genes of pGBM is still unclear. Given that the survival period of most glioma patients is relatively short, studying long-term survival patients with pGBM is of great value for this disease. Our study aim to identify key miRNAs with long-term survival related genes present in pGBM and uncover their potential mechanisms. The gene expression profiles of GSE53733, GSE15824, GSE30563, GSE50161 were obtained from the Gene Expression Omnibus database. Firstly, samples were divided into 3 groups according to its survival time and each group compare to the normal control group. Then we obtained differential expression genes (DEGs) with a long-term survival specific (LTSDEGs) and a short-term survival specific DEGs (STSDEGs). Next, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted with LTSDEGs and STSDEGs together. Moreover, we used the UALCAN database to verify LTSDEGs and STSDEGs, and obtained long-term verified survival specific DEGs(LTVSDEGs) and short-term verified survival specific DEGs(STVSDEGs). Finally, we established the predicted key miRNAs-LTVSDEGs interaction network. The protein expressions of the top 4 LTVSDEGs were verified in the HPA database with immunohistochemical staining. In total, we found 260 genes changed in LTSDEGs and 822 genes changed in STSDEGs. GO and KEGG results shown that the major changes are focused on tumor metabolism. 9 LTVSDEGs and 18 STVSDEGs were verified in UALCAN database. As for protein expression verification in top 4 LTVSDEGs, ZNF630, BLVRB and RPA3 were verified, while TPBG was not detected. We obtained 59 key miRNA from the predicted key miRNAs-LTVSDEGs interaction network. 25 key miRNAs were verified using GSE90603. Finally, we constructed the key miRNAs-LTVSDEGs network using a Sankey diagram, including 25 miRNAs and 7 LTVSDEGs. In conclusion, our study shows that there is a close relationship between metabolic changes and survival in pGBM. Besides, we established a key miRNAs-LTVSDEGs network for pGBM, which could be the key path in prolonging the life of pGBM patients.
