Bioinformatic analysis of long non-coding RNA-associated competing endogenous RNA network in adrenocortical carcinoma
Ontology highlight
ABSTRACT: Background
Adrenocortical carcinoma (ACC) is a malignant tumor with poor prognosis and unclear pathogenesis. This study aimed to explore the role of long non-coding RNAs (lncRNAs) in ACC. Methods
We obtained the lncRNA expression profiles of 10 ACC samples and 6 normal control samples from the GEO database and identified differentially expressed RNAs using the limma package in R. Results
We obtained a total of 391 differentially expressed lncRNAs (DElncRNAs) and 1,313 differentially expressed mRNAs (DEmRNAs) between ACC samples and normal control samples. Using Cytoscape v3.7.0, we then constructed a lncRNA-miRNA-mRNA (competing endogenous RNA, or ceRNA) network consisting of 87 lncRNAs, 31 miRNAs, and 78 mRNAs. Applying GO and KEGG enrichment analysis for 78 mRNAs in the ceRNA network, we identified 9 GO terms and 21 significantly enriched pathways. A PPI network was constructed using STRING online tools and Cytoscape v3.7.0, identifying 10 key genes. Finally, through Kaplan-Meier survival analysis, we identified five lncRNAs (LINC00887, MEIS1-AS2, MIR29B2CHG, MIR503HG, and SREBF2-AS1) associated with prognosis in patients with ACC. Conclusions
In summary, we constructed a ceRNA network and propose a new method for lncRNA research in ACC. Our results provide new clues for further exploration of lncRNAs in the pathogenesis of ACC.
SUBMITTER: Zhou Y
PROVIDER: S-EPMC8798464 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Objective: To construct ceRNA network and identify pivotal competing endogenous RNAs (ceRNAs) in adrenocortical carcinoma (ACC) using ceRNA network analysis. Methods: The RNA sequencing expression data of 77 ACCs in TCGA were obtained from GEPIA. Cancer specific ceRNAs, cancer specific microRNAs (miRNAs), and cancer specific messenger RNAs (mRNAs) were identified. The interaction of cancer specific miRNAs with cancer specific ceRNAs and cancer specific mRNAs were predicted. CeRNA network was constructed and visualized by Cytoscape 3.7.0 software. The genes in ceRNA network regulated GO terms and regulated pathways were performed by function analysis. Survival analysis of pivotal ceRNAs was performed for the pivotal lncRNAs. Result: Twenty-eight cancer specific ceRNAs, 149 cancer specific miRNAs, and 104 mRNAs were identified. CeRNA network was constructed including 10 ceRNAs, 35 miRNAs, and 34 mRNAs. The genes in ceRNA network regulated GO terms and were classified into three groups: cellular component (CC), molecular function (MF), and biological process (BP). The genes in ceRNA network regulated the following pathways: leukocyte transendothelial migration, and proteoglycans in cancer. Survival analysis showed that CTB-63M22.1 and RP1-241P17.4 were significantly associated with ACC patient disease free survival and overall survival. Conclusion: This study has constructed ceRNA networks in ACC. The study provides a set of pivotal ceRNAs for future investigation into the molecular mechanisms.
| S-EPMC6529643 | biostudies-literature
Project description:BackgroundThis study aimed to develop multi-RNA-based models using a competing endogenous RNA (ceRNA) regulatory network to provide survival risk prediction in head and neck squamous cell carcinoma (HNSCC).MethodsAll long non-coding RNA (lncRNA), microRNA (miRNA), and mRNA expression data and clinicopathological features related to HNSCC were derived from The Cancer Genome Atlas. Differentially expressed RNAs were calculated using R. Prognostic factors were identified using univariate Cox regression analysis. Functional analysis was performed using GO, KEGG pathways, and PPI network. Based on the results, we derived a risk signature and compared high- and low-risk subgroups using LASSO regression analysis. Survival analysis and the relationship between risk signature and clinicopathological features were performed using log-rank tests and Cox regression analysis. A ceRNA regulatory network was constructed, and prognostic lncRNAs and miRNA expression levels were validated in vitro and in vivo.ResultsA list of 207 lncRNAs, 18 miRNAs and 362 mRNAs related to overall survival was established. Five lncRNAs (HOTTIP, LINC00460, RMST, SFTA1P, and TM4SF19-AS1), one miRNA (hsa-miR-206), and one mRNA (STC2) were used to construct the ceRNA network. Three prognostic models contained 13 lncRNAs, eight miRNAs, and 17 mRNAs, which correlated with the patient status, disease-free survival (DFS), stage, grade, T stage, N stage, TP53 mutation status, angiolymphatic invasion, HPV status, and extracapsular spread. KEGG pathway analysis revealed significant enrichment of "Transcriptional misregulation in cancer" and "Neuroactive ligand-receptor interaction." In addition, HOTTIP, LINC00460, miR-206 and STC2 were validated in GTEx data, GEO microarrays and six HNSCC cell lines.ConclusionsOur findings clarify the interaction of ceRNA regulatory networks and crucial clinicopathological features. These results show that prognostic biomarkers can be identified by constructing multi-RNA-based prognostic models, which can be used for survival risk prediction in patients with HNSCC.
| S-EPMC7500352 | biostudies-literature
Project description:Schizophrenia (SCZ) is a serious psychiatric condition with a 1% lifetime risk. SCZ is one of the top ten global causes of disabilities. Despite numerous attempts to understand the function of genetic factors in SCZ development, genetic components in SCZ pathophysiology remain unknown. The competing endogenous RNA (ceRNA) network has been demonstrated to be involved in the development of many kinds of diseases. The ceRNA hypothesis states that cross-talks between coding and non-coding RNAs, including long non-coding RNAs (lncRNAs), via miRNA complementary sequences known as miRNA response elements, creates a large regulatory network across the transcriptome. In the present study, we developed a lncRNA-related ceRNA network to elucidate molecular regulatory mechanisms involved in SCZ. Microarray datasets associated with brain regions (GSE53987) and lymphoblasts (LBs) derived from peripheral blood (sample set B from GSE73129) of SCZ patients and control subjects containing information about both mRNAs and lncRNAs were downloaded from the Gene Expression Omnibus database. The GSE53987 comprised 48 brain samples taken from SCZ patients (15 HPC: hippocampus, 15 BA46: Brodmann area 46, 18 STR: striatum) and 55 brain samples taken from control subjects (18 HPC, 19 BA46, 18 STR). The sample set B of GSE73129 comprised 30 LB samples (15 patients with SCZ and 15 controls). Differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) were identified using the limma package of the R software. Using DIANA-LncBase, Human MicroRNA Disease Database (HMDD), and miRTarBase, the lncRNA- associated ceRNA network was generated. Pathway enrichment of DEmRNAs was performed using the Enrichr tool. We developed a protein-protein interaction network of DEmRNAs and identified the top five hub genes by the use of STRING and Cytoscape, respectively. Eventually, the hub genes, DElncRNAs, and predictive miRNAs were chosen to reconstruct the subceRNA networks. Our bioinformatics analysis showed that twelve key DEmRNAs, including BDNF, VEGFA, FGF2, FOS, CD44, SOX2, NRAS, SPARC, ZFP36, FGG, ELAVL1, and STARD13, participate in the ceRNA network in SCZ. We also identified DLX6-AS1, NEAT1, MINCR, LINC01094, DLGAP1-AS1, BABAM2-AS1, PAX8-AS1, ZFHX4-AS1, XIST, and MALAT1 as key DElncRNAs regulating the genes mentioned above. Furthermore, expression of 15 DEmRNAs (e.g., ADM and HLA-DRB1) and one DElncRNA (XIST) were changed in both the brain and LB, suggesting that they could be regarded as candidates for future biomarker studies. The study indicated that ceRNAs could be research candidates for investigating SCZ molecular pathways.
| S-EPMC8709859 | biostudies-literature
Project description:Long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNAs), interacting with microRNAs (miRNAs) and playing an important role in tumor progression. However, the role of lncRNA-mediated ceRNAs in glioma remains largely unknown. The present study aimed to identify novel lncRNAs and their associated function in glioma. RNA sequencing and corresponding clinical data from patients with glioma were obtained from The Cancer Genome Atlas. A total of 598 glioma tissues and 5 normal brain tissues were analyzed in the present study. The differentially expressed (DE) lncRNAs, mRNAs and miRNAs were identified using R packages and were used to construct a ceRNA network. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to investigate the biological functions of the DEmRNAs. Kaplan-Meier curve analysis was performed to investigate the association between DElncRNA expression and patient outcome. A total of 752 DElncRNAs, 2,079 DEmRNAs and 113 DEmiRNAs were identified between glioma and normal tissues. A lncRNA-miRNA-mRNA ceRNA network consisting of 61 lncRNAs, 12 miRNAs and 92 mRNAs was constructed. Survival analysis indicated that 36 DElncRNAs, 72 DEmRNAs and 3 DEmiRNAs were associated with overall survival in patients with glioma. The present study identified novel lncRNAs associated with survival prognosis and may facilitate further investigation of lncRNA-mediated ceRNA regulatory mechanisms in glioma.
| S-EPMC7436175 | biostudies-literature
Project description:Gastric cancer (GC) is a common type of cancer, and identification of novel diagnostic biomarkers associated with this disease is important. The present study aimed to identify novel diagnostic biomarkers associated with the prognosis of GC, using an integrated bioinformatics approach. Differentially expressed long non-coding RNAs (lncRNAs) associated with GC were identified using Gene Expression Omnibus datasets (GSE58828, GSE72305 and GSE99416) and The Cancer Genome Atlas database. A competing endogenous RNA network that incorporated five lncRNAs [long intergenic non-protein coding RNA 501 (LINC00501), LINC00365, SOX21 antisense divergent transcript 1 (SOX21-AS1), GK intronic transcript 1 (GK-IT1) and DLEU7 antisense RNA 1 (DLEU7-AS1)], 29 microRNAs and 114 mRNAs was constructed. Gene Ontology and protein-protein interaction network analyses revealed that these lncRNAs may be involved in 'biological regulation', 'metabolic process', 'cell communication', 'developmental process', 'cell proliferation', 'reproduction' and the 'cell cycle'. The results of receiver operating characteristic curve analysis demonstrated that LINC00501 (AUC=0.819), LINC00365 (AUC=0.580), SOX21-AS1 (AUC=0.736), GK-IT1 (AUC=0.823) and DLEU7-AS1 (AUC=0.932) had the potential to become valuable diagnostic biomarkers for GC. Associations with clinicopathological characteristics demonstrated that LINC00501 expression was significantly associated with sex (P=0.015) and tumor grade (P=0.022). Furthermore, LINC00365 expression was significantly associated with lymph node metastasis (P=0.025). Gene set enrichment analysis revealed that LINC00501, LINC00365 and SOX21-AS1 were enriched in signaling pathways associated with GC. Reverse transcription-quantitative PCR analysis demonstrated that LINC00501 expression (P=0.043) was significantly upregulated in GC tissues, whereas the expression levels of LINC00365 (P=0.033) and SOX21-AS1 (P=0.037) were significantly downregulated in GC tissues. Taken together, the results of the present study suggest that LINC00501, LINC00365, SOX21-AS1, GK-IT1 and DLEU7-AS1 may be used as novel diagnostic biomarkers for GC, and may be functionally associated with GC development and progression.
| S-EPMC8063256 | biostudies-literature
Project description:BackgroundOsteosarcoma is the primary bone malignant neoplasm that often develops metastasis. Increasing evidences have shown that non-coding RNAs (ncRNAs) relate to the progression of osteosarcoma. However, the ncRNAs' roles in osteosarcoma metastasis are still unknown.MethodsDifferentially expressed (DE) RNAs were identified from Gene Expression Omnibus (GEO) database. Protein-protein interaction (PPI) of DE messenger RNAs (DEmRNAs) was built through STRING database. The target mRNAs and long ncRNAs (lncRNAs) of microRNAs (miRNA) were predicted through miRDB, Targetscan and Genecode databases, which then cross-checked with previously obtained DERNAs to construct competing endogenous RNA (ceRNA) network. All networks were visualized via Cytoscape and the hub RNAs were screened out through Cytoscape plug-in Cytohubba. The gene functional and pathway analyses were performed through DAVID and MirPath databases. The survival analyses of hub RNAs were obtained through Kaplan-Meier (KM) survival curves.ResultsFive hundred sixty-four DEmRNAs, 16 DElncRNAs and 22 DEmiRNAs were screened out. GO functional and KEGG pathway analyses showed that DERNAs were significantly associated with tumor metastasis. The ceRNA network including 6 lncRNAs, 55 mRNAs and 20 miRNAs were constructed and the top 10 hub RNAs were obtained. Above all, PI3K/AKT signaling pathway was identified as the most important osteosarcoma metastasis-associated pathway and its hub ceRNA module was constructed. The survival analyses showed that the RNAs in hub ceRNA module closely related to osteosarcoma patients' prognosis.ConclusionsThe current study provided a new perspective on osteosarcoma metastasis. More importantly, the RNAs in hub ceRNA module might act as the novel therapeutic targets and prognostic factors for osteosarcoma patients.
| S-EPMC7887822 | biostudies-literature
Project description:Emerging evidence indicates that long noncoding RNAs (lncRNAs) play important roles in the regulation of cell differentiation by acting as competing endogenous RNA (ceRNA). However, the regulatory mechanisms of lncRNA and the lncRNA-associated ceRNA network involved in adipogenic differentiation of chicken preadipocytes remain elusive. Here, we first constructed the chicken preadipocyte in vitro induction model. Then, we identified differentially expressed lncRNAs (DELs), miRNAs (DEMis), and mRNAs (DEMs) between differentiated and undifferentiated preadipocytes. Furthermore, we constructed the lncRNA associated ceRNA network by gene expression correlation analysis and target prediction of DELs, DEMis, and DEMs. Finally, we determined twelve candidate lncRNA-miRNA-mRNA interactions from the lncRNA associated ceRNA network. Eight out of the twelve interactions were validated by RT-qPCR, indicating their potential role in the regulation of chicken preadipocytes differentiation. Among the eight interactions, TCONS_00026544-gga-miR-128-1-5p-RASD1, TCONS_00055280-gga-miR-135a-5p-JAM3, TCONS_00055280-gga-miR-135a-5p-GPR133, TCONS_00055280-gga-miR-135a-5p-CLDN1, and TCONS_00055280-gga-miR-135a-5p-TMEM123 may promote adipogenic differentiation of chicken preadipocytes while TCONS_00057272-gga-miR-146a-3p-FOXO6, TCONS_00057242-gga-miR-6615-3p-FOXO6, and TCONS_00057242-gga-miR-6615-3p-ENSGALT00000043224 have the opposite effects. Our results not only provide novel insights into ceRNA roles of lncRNAs in chicken preadipocytes differentiation and but also contribute to a better understanding of chicken fat deposition.
| S-EPMC6826404 | biostudies-literature
Project description:Background:Long non-coding RNAs (lncRNAs) play crucial roles in the initiation and progression of renal cell carcinoma (RCC) by competing in binding to miRNAs, and related competitive endogenous RNA (ceRNA) networks have been constructed in several cancers. However, the coexpression network has been poorly explored in RCC. Methods:We collected RCC RNA expression profile data and relevant clinical features from The Cancer Genome Atlas (TCGA). A cluster analysis was explored to show different lncRNA expression patterns. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and gene set enrichment analysis (GSEA) were performed to analyze the functions of the intersecting mRNAs. Targetscan and miRanda bioinformatics algorithms were used to predict potential relationships among RNAs. Univariate Cox proportional hazards regression was conducted to determine the RNA expression levels and survival times. Results:Bioinformatics analysis revealed that the expression profiles of hundreds of aberrantly expressed lncRNAs, miRNAs, and mRNAs were significantly changed between different stages of tumors and non-tumor groups. By combining the data predicted by databases with intersection RNAs, a ceRNA network consisting of 106 lncRNAs, 26 miRNAs and 69 mRNAs was established. Additionally, a protein interaction network revealed the main hub nodes (VEGFA, NTRK2, DLG2, E2F2, MYB and RUNX1). Furthermore, 63 lncRNAs, four miRNAs and 31 mRNAs were significantly associated with overall survival. Conclusion:Our results identified cancer-specific lncRNAs and constructed a ceRNA network for RCC. A survival analysis related to the RNAs revealed candidate biomarkers for further study in RCC.
| S-EPMC6054097 | biostudies-literature
Project description:Accumulating evidence has emerged revealing that noncoding RNAs (ncRNAs) play essential roles in the occurrence and development of hepatocellular carcinoma (HCC). However, the complicated regulatory interactions among various ncRNAs in the development of HCC are not entirely understood. The newly discovered mechanism of competing endogenous RNAs (ceRNAs) uncovered regulatory interactions among different varieties of RNAs. In recent years, a growing number of studies have suggested that ncRNAs, including long ncRNAs, circular RNAs and pseudogenes, play major roles in the biological functions of the ceRNA network in HCC. These ncRNAs can share microRNA response elements to affect microRNA affinity with target RNAs, thus regulating gene expression at the transcriptional level and both physiological and pathological processes. The ncRNAs that function as ceRNAs are involved in diverse biological processes in HCC cells, such as tumor cell proliferation, epithelial-mesenchymal transition, invasion, metastasis and chemoresistance. Based on these findings, ncRNAs that act as ceRNAs may be promising candidates for clinical diagnosis and treatments. In this review, we discuss the mechanisms and research methods of ceRNA networks. We also reviewed the recent advances in studying the roles of ncRNAs as ceRNAs in HCC and highlight possible directions and possibilities of ceRNAs as diagnostic biomarkers or therapeutic targets. Finally, the limitations, gaps in knowledge and opportunities for future research are also discussed.
| S-EPMC7552539 | biostudies-literature
Project description:Gastric cancer (GC) is one of the most frequently occurring life-threatening malignancies worldwide. Due to its high mortality rate, the discovery of putative biomarkers that may be sensitive and specific to GC is of seminal importance. Long non-coding RNAs (lncRNAs) are non-translatable RNAs whose transcript length exceeds 200 base pairs. The dysregulation of lncRNA expression plays a key role in tumorigenesis and development. In the present study, the expression profiles of lncRNAs, microRNAs and mRNAs of 361 GC tissues (and 32 normal gastric tissues) were downloaded from The Cancer Genome Atlas database. Furthermore, differentially expressed RNAs were analyzed by the DEseq package. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses confirmed some significant dysregulated signaling pathways and target RNAs. As a result, an lncRNA-associated competing endogenous RNA (ceRNA) network was constructed. Kaplan-Meier analysis of the differentially expressed RNAs associated with GC pathogenesis confirmed that the lncRNAs PVT1, HAND2-AS1 and ZNF667-AS1 were potentially associated with the prognosis of GC (P<0.05). The present study suggests the mechanism of ceRNA networks in GC, and further demonstrates that aberrant lncRNA expression may be used as an effective diagnostic tool (or target) for the prognosis of GC.
| S-EPMC7039062 | biostudies-literature