Dataset Information

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma.

ABSTRACT: RNAs may act as competing endogenous RNAs (ceRNAs), a critical mechanism in determining gene expression regulations in many cancers. However, the roles of ceRNAs in thyroid carcinoma remains elusive. In this study, we have developed a novel pipeline called Molecular Network-based Identification of ceRNA (MNIceRNA) to identify ceRNAs in thyroid carcinoma. MNIceRNA first constructs micro RNA (miRNA)-messenger RNA (mRNA)long non-coding RNA (lncRNA) networks from miRcode database and weighted correlation network analysis (WGCNA), based on which to identify key drivers of differentially expressed RNAs between normal and tumor samples. It then infers ceRNAs of the identified key drivers using the long non-coding competing endogenous database (lnCeDB). We applied the pipeline into The Cancer Genome Atlas (TCGA) thyroid carcinoma data. As a result, 598 lncRNAs, 1025 mRNAs, and 90 microRNA (miRNAs) were inferred to be differentially expressed between normal and thyroid cancer samples. We then obtained eight key driver miRNAs, among which hsa-mir-221 and hsa-mir-222 were key driver RNAs identified by both miRNA-mRNA-lncRNA and WGCNA network. In addition, hsa-mir-375 was inferred to be significant for patients' survival with 34 associated ceRNAs, among which RUNX2, DUSP6 and SEMA3D are known oncogenes regulating cellular proliferation and differentiation in thyroid cancer. These ceRNAs are critical in revealing the secrets behind thyroid cancer progression and may serve as future therapeutic biomarkers.

SUBMITTER: Lu M

PROVIDER: S-EPMC5793195 | biostudies-literature | 2018 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma.

Lu Minjia M Xu Xingyu X Xi Baohang B Dai Qi Q Li Chenli C Su Li L Zhou Xiaonan X Tang Min M Yao Yuhua Y Yang Jialiang J

Genes 20180119 1

RNAs may act as competing endogenous RNAs (ceRNAs), a critical mechanism in determining gene expression regulations in many cancers. However, the roles of ceRNAs in thyroid carcinoma remains elusive. In this study, we have developed a novel pipeline called Molecular Network-based Identification of ceRNA (MNIceRNA) to identify ceRNAs in thyroid carcinoma. MNIceRNA first constructs micro RNA (miRNA)-messenger RNA (mRNA)long non-coding RNA (lncRNA) networks from miRcode database and weighted correl ...[more]

PMID: 29351231

Similar Datasets

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:BackgroundHepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide. Recently, competing endogenous RNAs (ceRNA) have revealed a significant role in the progression of HCC. Herein, we aimed to construct a ceRNA network to identify potential biomarkers and illustrate its correlation with immune infiltration in HCC.MethodsRNA sequencing data and clinical traits of HCC patients were downloaded from TCGA. The limma R package was used to identify differentially expressed (DE) RNAs. The predicted prognostic model was established using univariate and multivariate Cox regression. A K-M curve, TISIDB and GEPIA website were utilized for survival analysis. Functional annotation was determined using Enrichr and Reactome. Protein-to-protein network analysis was implemented using SRTNG and Cytoscape. Hub gene expression was validated by quantitative polymerase chain reaction, Oncomine and the Hunan Protein Atlas database. Immune infiltration was analyzed by TIMMER, and Drugbank was exploited to identify bioactive compounds.ResultsThe predicted model that was established revealed significant efficacy with 3- and 5-years of the area under ROC at 0.804 and 0.744, respectively. Eleven DEmiRNAs were screened out by a K-M survival analysis. Then, we constructed a ceRNA network, including 56 DElncRNAs, 6 DEmiRNAs, and 28 DEmRNAs. The 28 DEmRNAs were enriched in cancer-related pathways, for example, the TNF signaling pathway. Moreover, six hub genes, CEP55, DEPDC1, KIF23, CLSPN, MYBL2, and RACGAP1, were all overexpressed in HCC tissues and independently correlated with survival rate. Furthermore, expression of hub genes was related to immune cell infiltration in HCC, including B cells, CD8+ T cells, CD4+ T cells, monocytes, macrophages, neutrophils, and dendritic cells.ConclusionThe findings from this study demonstrate that CEP55, DEPDC1, KIF23, CLSPN, MYBL2, and RACGAP1 are closely associated with prognosis and immune infiltration, representing potential therapeutic targets or prognostic biomarkers in HCC.

Project description:Phospholipase C epsilon 1 (PLCE1) and the competing endogenous RNA (ceRNA) network are crucial for tumorigenesis and the progression of esophageal squamous cell carcinoma (ESCC). However, whether PLCE1 can regulate the ceRNA network in ESCC has not been clarified. In the present study, we aimed to identify the PLCE1‑regulated ceRNA network and further elucidate the regulatory mechanisms by which ESCC is promoted. Microarray analysis was used to identify differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) from three pairs of samples of PLCE‑silenced Eca109 and control Eca109 cells. Next, the ceRNA regulatory network was established and visualized in Cytoscape, and functional enrichment analysis was performed to analyze DEGs from ceRNAs. Protein‑protein interaction (PPI) networks among the DEGs were established by the STRING database to screen hub genes. Kaplan‑Meier survival analysis was used to validate hub genes. Finally, PLCE1‑related hub gene/lncRNA/miRNA axes were also constructed based on the ceRNA network. A total of 105 DELs and 346 DEGs were found to be dysregulated in the microarray data (|log2FC| >1.5, adjusted P<0.05). We constructed a PLCE1‑regulated ceRNA network that incorporated 12 lncRNAs, 43 miRNAs, and 169 mRNAs. Functional enrichment analysis indicated that the DEGs might be associated with ESCC onset and development. A PPI network was established, and 9 hub genes [WD and tetratricopeptide repeats 1 (WDTC1), heat shock protein family A (Hsp70) member 5 (HSPA5), N‑ethylmaleimide sensitive factor, vesicle fusing ATPase (NSF), fibroblast growth factor 2 (FGF2), cyclin dependent kinase inhibitor 1A (CDKN1A or P21), bone morphogenetic protein 2 (BMP2), complement C3 (C3), GM2 ganglioside activator (GM2A) and discs large MAGUK scaffold protein 4 (DLG4)] were determined from the network. Kaplan‑Meier survival analysis validated four hub genes (BMP2, CDKN1A, GM2A, and DLG4) that were treated as prognostic factors. Ultimately, hub gene/lncRNA/miRNA subnetworks were obtained based on the 4 hub genes, 13 DEmiRNAs, and 10 DELs. In conclusion, the PLCE1‑regulated ceRNA contributes to the onset and progression of ESCC and the underlying molecular mechanisms may provide insights into personalized prognosis and new therapies for ESCC patients.

Project description:Liver cancer is the most common cancer and is the epitome of a recalcitrant cancer. Increasing evidence shown that long noncoding RNAs (lncRNA) were associated with cancer-related death and could function as a competing endogenous RNA (ceRNA). To explore regulatory roles and potential prognostic biomarkers of lncRNA for liver cancer, RNA-sequencing expression data were downloaded from the TCGA database and GEO database. A total of 357 patients were randomly divided into a discovery group and a validation group, of which 313 patients can obtain clinical data. In discovery phrase, 58 lncRNAs, 16 miRNAs, and 34 mRNAs were screened to construct the ceRNA network based on 252 patients employed from discovery group. Univariate and multivariate Cox hazard regression analysis model revealed that five lncRNAs (AATK-AS1, C10orf91, LINC00162, LINC00200, and LINC00501) from 58 lncRNAs were formulated to predict the overall survival (OS). We used the value of gene expression and regression coefficients to construct a risk score based on the five lncRNAs. Next, we validated our model in the GSE116174 dataset (n = 64) and the validation group (n = 94) from TCGA database. Subgroup analysis suggest that the five lncRNAs played critical parts in early stage in cancer from both discovery and validation groups. The five lncRNAs were also found to be associated with immune cells infiltration including CD4+ memory activated, NK cells activated and mast cells activated, then the results were also validated according to the validation group. Furthermore, KEGG pathway enrichment analysis showed that these nine coexpressed modules using the method of WGCNA, and many of these pathways are associated with the development and progression of disease. At last, the transcription factor binding sites (TFBS) of the five lncRNAs were predicted, which help us to understand the potential mechanism that the TFBS adjusted the ceRNA network. In summary, the ceRNA regulatory network may contribute to a better understanding of liver cancer mechanism and provide potential prognostic biomarkers and therapeutic targets.

Project description:The purpose of this work was to extract key players such as mRNAs and long non-coding RNA (lncRNAs) in the etiopathogenesis of osteosarcoma (OS). The sequencing analyses (mRNAs and lncRNAs) of OS were conducted followed by differentially expressed mRNAs and lncRNAs (DEmRNAs and DElncRNAs) identification between U-2OS cells with has-miR-590-5p overexpression and negative control cells. Following this, the co-expression and functional enrichment analyses of DEmRNAs and DElncRNAs were carried out. Also, the miRNAs-DElncRNAs-DEmRNAs regulatory network was constructed with DElncRNAs-miRNAs and DElncRNAs-DEmRNAs pairs after the target gene analysis of miRNA. In addition, the ceRNA-has-miR-590-5p was further extracted based on the has-miR-590-5p-DElncRNAs and DElncRNAs-DEmRNAs interactions. Finally, the results of the bioinformatics analysis was verified by reverse-transcription polymerase chain reaction (RT-PCR). Totally, 980 DEmRNAs (539 up-regulated DEmRNAs and 441 down-regulated DEmRNAs) and 682 DElncRNAs (352 up-regulated DElncRNAs and 330 down-regulated DElncRNAs) were extracted between cells with hsa-miR-590-5p overexpression and normal cells. The functional analyses suggested that up-regulated genes were significantly enriched in several GO terms such as signal transduction and cytokine-cytokine receptor interaction pathway while down-regulated genes (SCUBE3, HIST1H4E and EDIL3) were associated with calcium ion binding, cell surface function and nucleosome assembly. Additionally, the miRNAs-DEmRNAs-DEmRNAs network represented 220 pairs among 41 miRNAs, 38 DElncRNAs and 61 DEmRNAs. Furthermore, the ceRNA-hsa-miR-590-5p network consisted of 70 interaction pairs including hsa-miR-590-5p-SCUBE3-CTB-113D17.1, hsa-miR-590-5p-EDIL3-CTB-113D17.1 and hsa-miR-590-5p-HIST1H4E-CTB-113D17.1) among hsa-miR-590-5p, 30 DEmRNAs and 4 down-regulated DElncRNAs. Meanwhile, the RT-PCR results incidated that compared with the blank (KB) and negative control (NC) group, the mRNA expression of SCUBE3, HIST1H4E, and EDIL3 were significantly descreased in mimics group (P value <0.05). The lncRNA CTB-113D17.1 might implicate with OS development probably via serving as a hsa-miR-590-5p sponge to regulate gene targets (SCUBE3, EDIL3 and HIST1H4E), which will facilitate the deep understandings of OS progression.

Dataset Information

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma.

Publications

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets