Construction of a competitive endogenous RNA network and analysis of potential regulatory axis targets in glioblastoma

Abstract Background Glioblastoma is the most common primary malignant brain tumor. Because of the limited understanding of its pathogenesis, the prognosis of glioblastoma remains poor. This study was conducted to explore potential competing endogenous RNA (ceRNA) network chains and biomarkers in glioblastoma by performing integrated bioinformatics analysis. Methods Transcriptome expression data from The Cancer Genome Atlas database and Gene Expression Omnibus were analyzed to identify differentially expressed genes between glioblastoma and normal tissues. Biological pathways potentially associated with the differentially expressed genes were explored by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, and a protein-protein interaction network was established using the STRING database and Cytoscape. Survival analysis using Gene Expression Profiling Interactive Analysis was based on the Kaplan–Meier curve method. A ceRNA network chain was established using the intersection method to align data from four databases (miRTarBase, miRcode, TargetScan, and lncBace2.0), and expression differences and correlations were verified by quantitative reverse-transcription polymerase chain reaction analysis and by determining the Pearson correlation coefficient. Additionally, an MTS assay and the wound-healing and transwell assays were performed to evaluate the effects of complement C1s (C1S) on the viability and migration and invasion abilities of glioblastoma cells, respectively. Results We detected 2842 differentially expressed (DE) mRNAs, 2577 DE long non-coding RNAs (lncRNAs), and 309 DE microRNAs (miRNAs) that were dysregulated in glioblastoma. The final ceRNA network consisted of six specific lncRNAs, four miRNAs, and four mRNAs. Among them, four DE mRNAs and one DE lncRNA were correlated with overall survival (p