Dataset Information

Systematic analysis of the long noncoding RNA (lncRNA)-miRNA-mRNA competing endogenous RNA network to identify prognostic biomarkers and the potential regulatory axis in glioblastoma multiforme

ABSTRACT:

Background

Glioblastoma (GBM) is an intracranial brain tumor characterized by a high final lethality rate and recurrence rate, and limited available therapies. With the development of high-throughput sequencing technology, the genomic and transcriptomic features of GBM have been fully characterized. Therefore, our study aimed to identify its underlying genetic mechanisms, thus facilitating the development of novel therapies for GBM.

Methods

Based on the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases, differential expression of RNAs in GBM and control group was analyzed. After constructing the long noncoding RNA (lncRNA)-miRNA-mRNA regulatory network of GBM, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGGs) were performed to analyze related key nodes and the lncRNAs interacting with them. Further univariate Cox regression was conducted to explore independent factors, and then multivariate Cox regression was performed to construct risk prediction models.

Results

We first constructed the lncRNA-miRNA-mRNA regulatory network of GBM and two effective prediction models that included 2 mRNAs [transcription factor 12 (TCF12) and discoidin, CUB and LCCL domain containing 2 (DCBLD2)] and 5 lncRNAs (C10orf25, LINC00343, HOXA transcript antisense RNA, myeloid-specific 1 (HOTAIRM1), FGF12 antisense RNA 2 (FGF12-AS2) and H19). Additionally, we identified several key molecules [TCF12, integrin β3 (ITGB3), high mobility group AT-hook 2 (HMGA2), C10orf25 and LINC00336] closely associated with GBM prognosis. C10orf25/miR-218/DCBLD2 may be an important regulatory pathway in GBM.

Conclusions

Key molecules (TCF12, ITGB3, HMGA2, C10orf25, LINC00336 and H19) that are independent prognostic factors may be possible biomarkers to further optimize GBM prognosis. Two effective prognostic risk models that include 2 mRNAs (TCF12 and DCBLD2) and 5 lncRNAs (C10orf25, LINC00343, HOTAIRM1, FGF12-AS2 and H19) were constructed. C10orf25/miR-218/DCBLD2 may be an important regulatory pathway associated with the pathogenesis of GBM. Our findings contribute to further understanding the pathogenesis of GBM and finding possible candidate genes for prognostic and therapeutic usage with GBM.

SUBMITTER: Li X

PROVIDER: S-EPMC8798109 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundPulpitis is a common inflammatory disease that affects dental pulp. It is important to understand the molecular signals of inflammation and repair associated with this process. Increasing evidence has revealed that long noncoding RNAs (lncRNAs), via competitively sponging microRNAs (miRNAs), can act as competing endogenous RNAs (ceRNAs) to regulate inflammation and reparative responses. The aim of this study was to elucidate the potential roles of lncRNA, miRNA and messenger RNA (mRNA) ceRNA networks in pulpitis tissues compared to normal control tissues.MethodsThe oligo and limma packages were used to identify differentially expressed lncRNAs and mRNAs (DElncRNAs and DEmRNAs, respectively) based on expression profiles in two datasets, GSE92681 and GSE77459, from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were further analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein-protein interaction (PPI) networks and modules were established to screen hub genes using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the Molecular Complex Detection (MCODE) plugin for Cytoscape, respectively. Furthermore, an lncRNA-miRNA-mRNA-hub genes regulatory network was constructed to investigate mechanisms related to the progression and prognosis of pulpitis. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was applied to verify critical lncRNAs that may significantly affect the pathogenesis in inflamed and normal human dental pulp.ResultsA total of 644 upregulated and 264 downregulated differentially expressed genes (DEGs) in pulpitis samples were identified from the GSE77459 dataset, while 8 up- and 19 downregulated probes associated with lncRNA were identified from the GSE92681 dataset. Protein-protein interaction (PPI) based on STRING analysis revealed a network of DEGs containing 4,929 edges and 623 nodes. Upon combined analysis of the constructed PPI network and the MCODE results, 10 hub genes, including IL6, IL8, PTPRC, IL1B, TLR2, ITGAM, CCL2, PIK3CG, ICAM1, and PIK3CD, were detected in the network. Next, a ceRNA regulatory relationship consisting of one lncRNA (PVT1), one miRNA (hsa-miR-455-5p) and two mRNAs (SOCS3 and PLXNC1) was established. Then, we constructed the network in which the regulatory relationship between ceRNA and hub genes was summarized. Finally, our qRT-PCR results confirmed significantly higher levels of PVT1 transcript in inflamed pulp than in normal pulp tissues (p = 0.03).ConclusionOur study identified a novel lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of pulpitis.

Project description:In addition to direct targeting and repressing mRNAs, recent studies reported that microRNAs (miRNAs) can bridge up an alternative layer of post-transcriptional gene regulatory networks. The competing endogenous RNA (ceRNA) regulation depicts the scenario where pairs of genes (ceRNAs) sharing, fully or partially, common binding miRNAs (miRNA program) can establish coexpression through competition for a limited pool of the miRNA program. While the dynamics of ceRNA regulation among cellular conditions have been verified based on in silico and in vitro experiments, comprehensive investigation into the strength of ceRNA regulation in human datasets remains largely unexplored. Furthermore, pan-cancer analysis of ceRNA regulation, to our knowledge, has not been systematically investigated.In the present study we explored optimal conditions for ceRNA regulation, investigated functions governed by ceRNA regulation, and evaluated pan-cancer effects. We started by investigating how essential factors, such as the size of miRNA programs, the number of miRNA program binding sites, and expression levels of miRNA programs and ceRNAs affect the ceRNA regulation capacity in tumors derived from glioblastoma multiforme patients captured by The Cancer Genome Atlas (TCGA). We demonstrated that increased numbers of common targeting miRNAs as well as the abundance of binding sites enhance ceRNA regulation and strengthen coexpression of ceRNA pairs. Also, our investigation revealed that the strength of ceRNA regulation is dependent on expression levels of both miRNA programs and ceRNAs. Through functional annotation analysis, our results indicated that ceRNA regulation is highly associated with essential cellular functions and diseases including cancer. Furthermore, the highly intertwined ceRNA regulatory relationship enables constitutive and effective intra-function regulation of genes in diverse types of cancer.Using gene and microRNA expression datasets from TCGA, we successfully quantified the optimal conditions for ceRNA regulation, which hinge on four essential parameters of ceRNAs. Our analysis suggests optimized ceRNA regulation is related to disease pathways and essential cellular functions. Furthermore, although the strength of ceRNA regulation is dynamic among cancers, its governing functions are stably maintained. The findings of this report contribute to better understanding of ceRNA dynamics and its crucial roles in cancers.