Molecular Features of Glioma Determined and Validated Using Combined TCGA and GTEx Data Analyses

Abstract

Gliomas are among the most common intracranial tumors which originated from neuroepithelial cells. Increasing evidence has revealed that long noncoding RNA (lncRNA)-microRNA (miRNA)-mRNA module regulation and tumor-infiltrating immune cells play important regulatory roles in the occurrence and progression of gliomas. However, the precise underlying molecular mechanisms remain largely unknown. Data on gliomas in The Cancer Genome Atlas lack normal control samples; to overcome this limitation, we combined 665 The Cancer Genome Atlas glioma RNA sequence datasets with 188 Genotype-Tissue Expression normal brain RNA sequences to construct an expression matrix profile after normalization. We systematically analyzed the expression of mRNAs, lncRNAs, and miRNAs between gliomas and normal brain tissues. Kaplan–Meier survival analyses were conducted to screen differentially expressed mRNAs, lncRNAs, and miRNAs. A prognostic miRNA-related competitive endogenous RNA network was constructed, and the core subnetworks were filtered using 6 miRNAs, 3 lncRNAs, and 11 mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to investigate the biological functions of significantly dysregulated mRNAs. Co-expression network analysis was performed to analyze and screen the core genes. Furthermore, single-sample Gene Set Enrichment Analysis and immune checkpoint gene expression analysis were performed, as co-expression analysis indicated immune gene dysregulation in glioma. Finally, the expression of representative dysregulated genes was validated in U87 cells at the transcriptional level, establishing a foundation for further research. We identified 7017 mRNAs, 437 lncRNAs, and 9 miRNAs that were differentially expressed in gliomas. Kaplan–Meier survival analysis revealed 5684 mRNAs, 61 lncRNAs, and 7 miRNAs with potential as prognostic signatures in patients with glioma. The hub subnetwork of the competing endogenous RNA network between PART1-hsa-mir-25-SLC12A5/TACC2/BSN/TLN2/ZDHHC8 was screened out. Gene co-expression network, single-sample Gene Set Enrichment Analysis, and immune checkpoint expression analysis demonstrated that tumor-infiltrating immune cells are closely related to gliomas. We identified novel potential biomarkers to predict survival and therapeutic targets for patients with gliomas based on a large-scale sample. Importantly, we filtered pivotal genes that provide valuable information for further exploration of the molecular mechanisms underlying glioma tumorigenesis and progression.