Crosstalk in competing endogenous RNA network reveals the complex molecular mechanism underlying lung cancer.
Ontology highlight
ABSTRACT:
SUBMITTER: Jin X
PROVIDER: S-EPMC5710922 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:microRNAs (miRNAs) regulate gene expression at post-transcriptional level by repressing target RNA molecules. Competition to bind miRNAs tends in turn to correlate their targets, establishing effective RNA-RNA interactions that can influence expression levels, buffer fluctuations and promote signal propagation. Such a potential has been characterized mathematically for small motifs both at steady state and away from stationarity. Experimental evidence, on the other hand, suggests that competing endogenous RNA (ceRNA) crosstalk is rather weak. Extended miRNA-RNA networks could however favour the integration of many crosstalk interactions, leading to significant large-scale effects in spite of the weakness of individual links. To clarify the extent to which crosstalk is sustained by the miRNA interactome, we have studied its emergent systemic features in silico in large-scale miRNA-RNA network reconstructions. We show that, although generically weak, system-level crosstalk patterns (i) are enhanced by transcriptional heterogeneities, (ii) can achieve high-intensity even for RNAs that are not co-regulated, (iii) are robust to variability in transcription rates, and (iv) are significantly non-local, i.e. correlate weakly with miRNA-RNA interaction parameters. Furthermore, RNA levels are generically more stable when crosstalk is strongest. As some of these features appear to be encoded in the network's topology, crosstalk may functionally be favoured by natural selection. These results suggest that, besides their repressive role, miRNAs mediate a weak but resilient and context-independent network of cross-regulatory interactions that interconnect the transcriptome, stabilize expression levels and support system-level responses.
| S-EPMC6853376 | biostudies-literature
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:BackgroundAcute lung injury (ALI) is a fatal syndrome frequently induced by lipopolysaccharide (LPS) released from the bacterial cell wall. LPS could also trigger autophagy of lung bronchial epithelial cell to relieve the inflammation, while the overwhelming LPS would impair the balance of autophagy consequently inducing serious lung injury.MethodsWe observed the autophagy variation of 16HBE, human bronchial epithelial cell, under exposure to different concentrations of LPS through western blot, immunofluorescence staining, and electron microscopy. Eight strands of 16HBE were divided into two groups upon 1000 ng/ml LPS stimulation or not, which were sent to be sequenced at whole transcriptome. Subsequently, we analyzed the sequencing data in functional enrichment, pathway analysis, and candidate gene selection and constructed a hsa-miR-663b-related competing endogenous RNA (ceRNA) network.ResultsWe set a series of concentrations of LPS to stimulate 16HBE and observed the variation of autophagy in related protein expression and autophagosome count. We found that the effective concentration of LPS was 1000 ng/ml at 12 hours of exposure and sequenced the 1000 ng/ml LPS-stimulated 16HBE. As a result, a total of 750 differentially expressed genes (DEGs), 449 differentially expressed lncRNAs (DElncRNAs), 76 differentially expressed circRNAs (DEcircRNAs), and 127 differentially expressed miRNAs (DEmiRNAs) were identified. We constructed the protein-protein interaction (PPI) network to visualize the interaction between DEGs and located 36 genes to comprehend the core discrepancy between LPS-stimulated 16HBE and the negative control group. In combined analysis of differentially expressed RNAs (DERNAs), we analyzed all the targeted relationships of ceRNA in DERNAs and figured hsa-miR-663b as a central mediator in the ceRNA network to play when LPS induced the variation of autophagy in 16HBE.ConclusionOur research indicated that the hsa-miR-663b-related ceRNA network may contribute to the key regulatory mechanism in LPS-induced changes of autophagy and ALI.
| S-EPMC8553468 | biostudies-literature
Project description:Increasing evidence has shown that noncoding RNAs play significant roles in the initiation, progression, and metastasis of tumours via participating in competing endogenous RNA (ceRNA) networks. However, the survival-associated ceRNA in lung adenocarcinoma (LUAD) remains poorly understood. In this study, we aimed to investigate the regulatory mechanisms underlying ceRNA in LUAD to identify novel prognostic factors. mRNA, lncRNA, and miRNA sequencing data obtained from the GDC data portal were utilized to identify differentially expressed (DE) RNAs. Survival-related RNAs were recognized using univariate Kaplan-Meier survival analysis. We performed functional enrichment analysis of survival-related mRNAs using the clusterProfiler package of R and STRING. lncRNA-miRNA and miRNA-mRNA interactions were predicted based on miRcode, Starbase, and miRanda. Subsequently, the survival-associated ceRNA network was constructed for LUAD. Multivariate Cox regression analysis was used to identify prognostic factors. Finally, we acquired 15 DE miRNAs, 49 DE lncRNAs, and 843 DE mRNAs associated with significant overall survival. Functional enrichment analysis indicated that survival-related DE mRNAs were enriched in cell cycle. The survival-associated lncRNA-miRNA-mRNA ceRNA network was constructed using five miRNAs, 49 mRNAs, and 21 lncRNAs. Furthermore, seven hub RNAs (LINC01936, miR-20a-5p, miR-31-5p, TNS1, TGFBR2, SMAD7, and NEDD4L) were identified based on the ceRNA network. LINC01936 and miR-31-5p were found to be significant using the multifactorial Cox regression model. In conclusion, we successfully constructed a survival-related lncRNA-miRNA-mRNA ceRNA regulatory network in LUAD and identified seven hub RNAs, which provide novel insights into the regulatory molecular mechanisms associated with survival of LUAD, and identified two independent prognostic predictors for LUAD.
| S-EPMC7895565 | biostudies-literature
Project description:BackgroundEsophageal cancer (ESCA) constitutes one of the most common cancers worldwide. The identification of potential biomarkers is important to improving the diagnostic accuracy and treatment efficiency for patients with ESCA. In this study, we aimed to identify biomarkers related to ESCA progression through a comprehensive analysis of long non-coding RNAs (lncRNAs), microRNA (miRNAs), and mRNA expression profiles in ESCA.MethodsDifferentially expressed lncRNAs, miRNAs, and mRNAs (DElncRNAs, DEmiRNAs, and DEmRNAs, respectively) in ESCA samples compared with normal controls were obtained. A competing endogenous RNA (ceRNA) network consisting of interacting DElncRNAs, DEmiRNAs, and DEmRNAs was constructed using a combination of the miRCode and TargetScan databases. Relationships between RNAs in the ceRNA network and overall survival in patients with EC were explored through another independent ESCA dataset from The Cancer Genome Atlas.ResultsA total of 1,014 DElncRNAs, 3,677 DEmRNAs, and 35 DEmiRNAs were identified in ESCA samples compared with normal samples. Functional enrichment analysis indicated that the DEmRNAs were involved in cell activity, inflammatory response, and oxygen metabolism-related biological processes. A ceRNA network containing 5 DEmiRNAs, 582 DEmRNAs and 764 DElncRNAs was obtained. In the survival analysis, 39 genes were found to be significantly associated with overall survival in patients with EC, including GOLGA7, NFYB, TOP1, and TMTC3.ConclusionsOur study constructed a ceRNA network for ESCA for the first time, which will be helpful for the disease's diagnosis and treatment.
| S-EPMC8506737 | biostudies-literature
Project description:Evidence is increasingly indicating that circular RNAs (circRNAs) are closely involved in tumorigenesis and cancer progression. However, the function and application of circRNAs in lung adenocarcinoma (LUAD) are still unknown. In this study, we constructed a circRNA-associated competitive endogenous RNA (ceRNA) network to investigate the regulatory mechanism of LUAD procession and further constructed a prognostic signature to predict overall survival for LUAD patients. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were selected to construct the ceRNA network. Based on the TargetScan prediction tool and Pearson correlation coefficient, we constructed a circRNA-associated ceRNA network including 11 DEcircRNAs, 8 DEmiRNAs, and 49 DEmRNAs. GO and KEGG enrichment indicated that the ceRNA network might be involved in the regulation of GTPase activity and endothelial cell differentiation. After removing the discrete points, a PPI network containing 12 DEmRNAs was constructed. Univariate Cox regression analysis showed that three DEmRNAs were significantly associated with overall survival. Therefore, we constructed a three-gene prognostic signature for LUAD patients using the LASSO method in the TCGA-LUAD training cohort. By applying the signature, patients could be categorized into the high-risk or low-risk subgroups with significant survival differences (HR: 1.62, 95% CI: 1.12-2.35, log-rank p = 0.009). The prognostic performance was confirmed in an independent GEO cohort (GSE42127, HR: 2.59, 95% CI: 1.32-5.10, log-rank p = 0.004). Multivariate Cox regression analysis proved that the three-gene signature was an independent prognostic factor. Combining the three-gene signature with clinical characters, a nomogram was constructed. The primary and external verification C-indexes were 0.717 and 0.716, respectively. The calibration curves for the probability of 3- and 5-year OS showed significant agreement between nomogram predictions and actual observations. Our findings provided a deeper understanding of the circRNA-associated ceRNA regulatory mechanism in LUAD pathogenesis and further constructed a useful prognostic signature to guide personalized treatment of LUAD patients.
| S-EPMC8421160 | biostudies-literature
Project description:BackgroundThe establishment of immunotherapy has led to a new era in oncotherapy. But the signature of immune-related genes (IRGs) in LUAD remains to be elucidated. Here we use integrated analysis to identify IRGs roles in immune signature and detect their relationship with competing endogenous RNA (ceRNA) networks in LUAD progression.MethodsBy analyzing the RNA-seq data from different platforms, we recognized the differentially expressed genes (DEGs) of each platform and screened out the top 20 hub IRGs related to immune responses. Then, we applied the CIBERSORT algorithm to explore the landscape of tumor-infiltrating immune cells (TILs) in LUAD and their connection with hub genes. Next, we predicted and validated the upstream miRNAs and lncRNAs according to their expression and prognostic roles. Finally, we constructed and validated an immune-related ceRNA network by co-expression analysis.ResultsA total of 71 IRGs were identified among 248 DEGs, which play key roles in immune responses. CIBERSORT analysis showed that six hub genes were closely related to TILs, such as SPP1 and naive B cells (R = -0.17), TEK and resting mast cells (R = 0.37). Stepwise prediction and validation from mRNA to lncRNA, including 6 hub genes, 5 miRNAs, and 9 lncRNAs, were applied to construct a ceRNA network. Ultimately, we confirmed the TMPO-AS1/miR-126-5p/SPP1 and CARD8-AS1/miR-21-5p/TEK as immune-related ceRNA networks in LUAD progression.ConclusionWe elucidated two immune-related ceRNA networks in LUAD progression, which can be considered as immunotherapy targets for this disease.
| S-EPMC8209499 | biostudies-literature
Project description:BackgroundThe events in early HIV infection (EHI) are important determinants of disease severity and progression rate to AIDS, but the mechanisms of pathogenesis in EHI have not been fully understood. Circular RNAs (circRNAs) have been verified as "microRNA sponges" that regulate gene expression through competing endogenous RNA (ceRNA) networks, but circRNA expression profiles and their contribution to EHI pathogenesis are still unclear.MethodsTwo different libraries were constructed with RNA from human peripheral blood mononuclear cells from 3 HARRT-naive EHI patients and 3 healthy controls (HCs). The complete transcriptomes were sequenced with RNA sequencing (RNA-Seq) and miRNA sequencing (miRNA-Seq). The differentially expressed (DE) RNAs were validated with RT-qPCR. The circRNA profile and circRNA-associated-ceRNA network in EHI were analyzed with the integrated data of RNA-Seq and miRNA-Seq. Gene ontology (GO) analysis was used to annotate the circRNAs involved in the circRNA-associated-ceRNA networks.ResultsA total of 1365 circRNAs, 30 miRNAs, and 2049 mRNAs were differentially expressed between HARRT-naive EHI patients and HCs. A ceRNA network was constructed with 516 DE circRNAs and 903 DE mRNAs that shared miR response elements with 21 DE miRNAs. GO analysis demonstrated the multiple roles of the circRNAs enriched in EHI with circRNA-associated-ceRNA networks, such as immune response, inflammatory response and defense responses to virus, 67 circRNAs were revealed to be potentially involved in HIV-1 replication through regulating the expression of CCNK, CDKN1A and IL-15.ConclusionsThis study, for the first time, revealed a large circRNA profile and complex pathogenesis roles of circRNAs in EHI. A group of enriched circRNAs and associated circRNA-associated-ceRNA networks might contribute to HIV replication regulation and provide novel potential targets for both the pathogenesis of EHI and antiviral therapy.
| S-EPMC6264784 | biostudies-literature
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.
| S-EPMC6059698 | biostudies-literature
Project description:A new form of circuitry for gene regulation has been identified in which RNAs can crosstalk by competing for shared microRNAs (miRNAs). Such competing endogenous RNAs (ceRNAs) form a network via shared miRNA response elements (MREs) to antagonize miRNA function. We previously reported natural antisense RNA (AS) as an important modulator of interferon-α1 (IFN-α1) mRNA levels by promoting IFN-α1 mRNA stability. We show that IFN-α1 AS forms a ceRNA network with specific IFN-α AS (IFN-α7/-α8/-α10/-α14) and mRNA (IFN-α8/-α10/-α14/-α17) subtypes from the IFN-α gene (IFNA) family to antagonize miRNA-1270 (miR-1270), thereby modulating IFN-α1 mRNA levels. Bioinformatic analysis demonstrated that IFN-α1 AS harbors multiple miR-1270 MREs (MRE-1270s), whose presence was substantiated by miR-1270 overexpression and transfection of antimiR-1270. The antimiR-1270, complementary to the miR-1270 seed region, revealed that IFN-α1 AS likely shares the MRE-1270 with IFN-α1 mRNA and specific IFN-α AS and mRNA subtypes. Subsequent bioinformatic analysis for MRE-1270s showed that IFN-α1 AS and other RNA subtypes shared the 6-mer MRE-1270 site. Further MRE-mapping demonstrated that the total number of MRE-1270s in IFN-α1 AS accounted for approximately 30 % of the miR-1270 population. AntimiR-1270 transfection also caused specific de-repression of five cellular mRNAs, including that of CAPRIN1. These results suggest that IFN-α1 AS, together with specific IFN-α AS and mRNA subtypes, as well as the five cellular mRNAs, participate as competing molecules in the ceRNA network against miR-1270. This coordinated regulatory architecture suggests a vital function for the innate immune system in maintaining precise physiological type I IFN levels via post-transcriptional regulatory mechanisms.
| S-EPMC4477080 | biostudies-literature