Dataset Information

Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers.

ABSTRACT: 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.

SUBMITTER: Chiu YC

PROVIDER: S-EPMC4416191 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers.

Chiu Yu-Chiao YC Hsiao Tzu-Hung TH Chen Yidong Y Chuang Eric Y EY

BMC genomics 20150421

<h4>Background</h4>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 a ...[more]

PMID: 25917195

Similar Datasets

Project description:BackgroundOsteoarthritis (OA) is one of the most common diseases in elderly people; however, the correlation between molecular alterations and the occurrence and progression of OA are still not well understood. We conducted this study to investigate the molecular changes in OA via the competing endogenous ribonucleic acid (ceRNA) network.MethodsWe downloaded the messenger RNA (mRNA) data set, GSE48556, the microRNA (miRNA) data set, GSE105027, and the long non-coding (lncRNA) data set, GSE126963 from the Gene Expression Omnibus (GEO) database, and examined the differentially expressed genes (DEGs) in these data sets. Further, we constructed a ceRNA network of the differentially expressed miRNAs, mRNAs, and lncRNAs. To determine the biological functions of the ceRNA network, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses. Finally, we conducted an immune cell infiltration analysisusing single-sample gene set enrichment analysis to examine the abundance of immune cells in healthy and OA patients, and compared the infiltration of 28 immune cells between the healthy and OA samples. We also analyzed the relationship between the abundance of immune cells and mRNA expression levels in the ceRNA network.ResultsUltimately hsa-mir-425-3p, dual specificity phosphatase 1, and 24 lncRNAs were identified in the ceRNA network. The functional enrichment analyses showed that these lncRNAs, miRNAs, and mRNAs are involved in various significant biological process, such as the regulation of leukocyte migration, Mitogen-Activated Protein (MAP) kinase tyrosine/serine/threonine phosphatase activity, the interleukin-17 signaling pathway, the tumor necrosis factor signaling pathway, and osteoclast differentiation, and can also have a strong effect on immune cell infiltration.ConclusionsThe dual-specificity phosphatase 1-specific ceRNA network can be used as a diagnostic tool to assess the progression of OA patients.

Project description: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.

Project description:The present study aimed to screen all types of RNAs involved in the development of papillary thyroid carcinoma (PTC). RNA‑sequencing data of PTC and normal samples were used for screening differentially expressed (DE) microRNAs (DE‑miRNAs), long non‑coding RNAs (DE‑lncRNAs) and genes (DEGs). Subsequently, lncRNA‑miRNA, miRNA‑gene (that is, miRNA‑mRNA) and gene‑gene interaction pairs were extracted and used to construct regulatory networks. Feature genes in the miRNA‑mRNA network were identified by topological analysis and recursive feature elimination analysis. A support vector machine (SVM) classifier was built using 15 feature genes, and its classification effect was validated using two microarray data sets that were downloaded from the Gene Expression Omnibus (GEO) database. In addition, Gene Ontology function and Kyoto Encyclopedia Genes and Genomes pathway enrichment analyses were conducted for genes identified in the ceRNA network. A total of 506 samples, including 447 tumor samples and 59 normal samples, were obtained from The Cancer Genome Atlas (TCGA); 16 DE‑lncRNAs, 917 DEGs and 30 DE‑miRNAs were screened. The miRNA‑mRNA regulatory network comprised 353 nodes and 577 interactions. From these data, 15 feature genes with high predictive precision (>95%) were extracted from the network and were used to form an SVM classifier with an accuracy of 96.05% (486/506) for PTC samples downloaded from TCGA, and accuracies of 96.81 and 98.46% for GEO downloaded data sets. The ceRNA regulatory network comprised 596 lines (or interactions) and 365 nodes. Genes in the ceRNA network were significantly enriched in 'neuron development', 'differentiation', 'neuroactive ligand‑receptor interaction', 'metabolism of xenobiotics by cytochrome P450', 'drug metabolism' and 'cytokine‑cytokine receptor interaction' pathways. Hox transcript antisense RNA, miRNA‑206 and kallikrein‑related peptidase 10 were nodes in the ceRNA regulatory network of the selected feature gene, and they may serve import roles in the development of PTC.

Project description:Background: Glioblastoma multiforme (GBM), the most common and aggressive human malignant brain tumor, is notorious for its limited treatment options and poor prognosis. MicroRNAs (miRNAs) are found to be involved in tumorigenesis of GBM. However, a comprehensive miRNA-mRNA regulatory network has still not been established. Methods: A miRNA microarray dataset (GSE90603) was obtained from GEO database. Then, we employed GEO2R tool to perform differential expression analysis. Potential transcription factors and target genes of screened differentially expressed miRNAs (DE-miRNAs) were predicted. The GBM mRNA dataset were downloaded from TCGA database for identifying differentially expressed genes (DEGs). Next, GO annotation and KEGG pathway enrichment analysis was conducted. PPI network was then established, and hub genes were identified via Cytoscape software. The expression and prognostic roles of hub genes was further evaluated. Results: Total 33 DE-miRNAs, consisting of 10 upregulated DE-miRNAs and 23 downregulated DE-miRNAs, were screened. SP1 was predicted to potentially regulate most of screened DE-miRNAs. Three thousand and twenty seven and 3,879 predicted target genes were obtained for upregulated and downregulated DE-miRNAs, respectively. Subsequently, 1,715 upregulated DEGs and 1,259 downregulated DEGs were identified. Then, 149 and 295 potential downregulated and upregulated genes commonly appeared in target genes of DE-miRNAs and DEGs were selected for GO annotation and KEGG pathway enrichment analysis. The downregulated genes were significantly enriched in cGMP-PKG signaling pathway and calcium signaling pathway whereas the upregulated genes were enriched in pathways in cancer and PI3K-Akt signaling pathway. Construction and analysis of PPI network showed that STXBP1 and TP53 were recognized as hub genes with the highest connectivity degrees. Expression analytic result of the top 20 hub genes in GBM using GEPIA database was generally identical with previous differential expression analysis for TCGA data. EGFR, PPP3CB, and MYO5A expression was significantly associated with patients' OS. Conclusions: In this study, we established a potential GBM-related miRNA-mRNA regulatory network, which explores a comprehensive understanding of the molecular mechanisms and provides key clues in seeking novel therapeutic targets for GBM. In the future, more experiments need to be performed to validate our current findings.

Project description:BackgroundGlioblastoma (GBM) has a high incidence rate, invasive growth, and easy recurrence, and the current therapeutic effect is less than satisfying. Pyroptosis plays an important role in morbidity and progress of GBM. Meanwhile, the tumor microenvironment (TME) is involved in the progress and treatment tolerance of GBM. In the present study, we analyzed prognosis model, immunocyte infiltration characterization, and competing endogenous RNA (ceRNA) network of GBM on the basis of pyroptosis-related genes (PRGs).MethodsThe transcriptome and clinical data of 155 patients with GBM and 120 normal subjects were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx). Lasso (Least absolute shrinkage and selection operator) Cox expression analysis was used in predicting prognostic markers, and its predictive ability was tested using a nomogram. A prognostic risk score formula was constructed, and CIBERSORT, ssGSEA algorithm, Tumor IMmune Estimation Resource (TIMER), and TISIDB database were used in evaluating the immunocyte infiltration characterization and tumor immune response of differential risk samples. A ceRNA network was constructed with Starbase, mirtarbase, and lncbase, and the mechanism of this regulatory axis was explored using Gene Set Enrichment Analysis (GSEA).ResultsFive PRGs (CASP3, NLRP2, TP63, GZMB, and CASP9) were identified as the independent prognostic biomarkers of GBM. Prognostic risk score formula analysis showed that the low-risk group had obvious survival advantage compared with the high-risk group, and significant differences in immunocyte infiltration and immune related function score were found. In addition, a ceRNA network of messenger RNA (CASP3, TP63)-microRNA (hsa-miR-519c-5p)-long noncoding RNA (GABPB1-AS1) was established. GSEA analysis showed that the regulatory axis played a considerable role in the extracellular matrix (ECM) and immune inflammatory response.ConclusionsPyroptosis and TME-related independent prognostic markers were screened in this study, and a prognosis risk score formula was established for the first time according to the prognosis PRGs. TME immunocyte infiltration characterization and immune response were assessed using ssGSEA, CIBERSORT algorithm, TIMER, and TISIDB database. Besides a ceRNA network was built up. This study not only laid foundations for further exploring pyroptosis and TME in improving prognosis of GBM, but also provided a new idea for more effective guidance on clinical immunotherapy to patients and developing new immunotherapeutic drugs.

Project description:BackgroundPancreatic cancer is an insidious and heterogeneous malignancy with poor prognosis that is often locally unresectable. Therefore, determining the underlying mechanisms and effective prognostic indicators of pancreatic cancer may help optimize clinical management. This study was conducted to develop a prognostic model for pancreatic cancer based on a competing endogenous RNA (ceRNA) network.MethodsWe obtained transcriptomic data and corresponding clinicopathological information of pancreatic cancer samples from The Cancer Genome Atlas (TCGA) database (training set). Based on the ceRNA interaction network, we screened candidate genes to build prediction models. Univariate Cox regression analysis was performed to screen for genes associated with prognosis, and least absolute shrinkage and selection operator (LASSO) regression analysis was conducted to construct a predictive model. A receiver operating characteristic (ROC) curve was drawn, and the C-index was calculated to evaluate the accuracy of the prediction model. Furthermore, we downloaded transcriptomic data and related clinical information of pancreatic cancer samples from the Gene Expression Omnibus database (validation set) to evaluate the robustness of our prediction model.ResultsEight genes (ANLN, FHDC1, LY6D, SMAD6, ACKR4, RAB27B, AUNIP, and GPRIN3) were used to construct the prediction model, which was confirmed as an independent predictor for evaluating the prognosis of patients with pancreatic cancer through univariate and multivariate Cox regression analysis. By plotting the decision curve, we found that the risk score model is an independent predictor has the greatest impact on survival compared to pathological stage and targeted molecular therapy.ConclusionsAn eight-gene prediction model was constructed for effectively and independently predicting the prognosis of patients with pancreatic cancer. These eight genes identified show potential as diagnostic and therapeutic targets.

Dataset Information

Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers.

Publications

Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets