Competing endogenous RNA network for esophageal cancer progression
Ontology highlight
ABSTRACT: Background
Esophageal 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. Methods
Differentially 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. Results
A 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. Conclusions
Our study constructed a ceRNA network for ESCA for the first time, which will be helpful for the disease’s diagnosis and treatment.
SUBMITTER: Chen S
PROVIDER: S-EPMC8506737 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BACKGROUND:It is well acknowledged that cancer-related pathways play pivotal roles in the progression of pancreatic cancer (PC). Employing Integrated analysis, we aim to identify the pathway-related ceRNA network associated with PC progression. METHODS:We divided eight GEO datasets into three groups according to their platform, and combined TCGA and GTEx databases as a group. Additionally, we screened out the differentially expressed genes (DEGs) and performed functional enrichment analysis in each group, and recognized the top hub genes in the most enriched pathway. Furthermore, the upstream of miRNAs and lncRNAs were predicted and validated according to their expression and prognostic roles. Finally, the co-expression analysis was applied to identify a pathway-related ceRNA network in the progression of PC. RESULTS:A total of 51 significant pathways that common enriched in all groups were spotted. Enrichment analysis indicated that pathway in cancer was greatly linked with tumor formation and progression. Next, the top 20 hug genes in this pathway were recognized, and stepwise prediction and validation from mRNA to lncRNA, including 11 hub genes, 4 key miRNAs, and 2 key lncRNAs, were applied to identify a meaningful ceRNA network according to ceRNA rules. Ultimately, we identified the PVT1/miR-20b/CCND1 axis as a promising pathway-related ceRNA axis in the progression of PC. CONCLUSION:Overall, we elucidate the pathway-related ceRNA regulatory network of PVT1/miR-20b/CCND1 in the progression of PC, which can be considered as therapeutic targets and encouraging prognostic biomarkers for PC.
| S-EPMC7532576 | biostudies-literature
Project description:As the most common type of cancer in female patients, the morbidity and mortality rates of breast cancer (BC) are high, and its incidence is gradually increasing worldwide. However, the underlying molecular and genetic mechanisms involved in the etiopathogenesis of BC remain unclear. Circular RNAs (circRNAs) are a novel type of non-coding RNAs that have been verified to serve a crucial role in tumorigenesis. However, the majority of functions and mechanisms of circRNAs remain unknown. The present study identified 47 differentially expressed circRNAs in a dataset from Gene Expression Omnibus. Using the cancer-specific circRNA database, the potential microRNA (miRNA) response elements, RNA-binding proteins and open reading frames of the candidate circRNAs were predicted. Combing the predictions of miRNAs and target mRNAs, a competing endogenous RNA network was constructed, which may serve as the theoretical basis for further research. Furthermore, the analyses conducted using Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways indicated that candidate circRNAs may serve a role in transcriptional regulation. Moreover, 20 BC tissue specimens and their paired adjacent normal tissue specimens were used to evaluate the expression levels of the screened circRNAs. Thus, the analyses of the raw microarray data conducted in the present study offer perspectives on the exploration of mechanisms associated with BC tumorigenesis with regard to the circRNA-miRNA-mRNA network.
| S-EPMC6960389 | biostudies-literature
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.
| S-EPMC7859920 | biostudies-literature
Project description:Non-small-cell lung cancer (NSCLC) is the most prevalent histological cancer subtype worldwide. As the majority of patients present with invasive, metastatic disease, it is vital to understand the basis for lung cancer progression. Hmga2 is highly expressed in metastatic lung adenocarcinoma, in which it contributes to cancer progression and metastasis. Here we show that Hmga2 promotes lung cancer progression in mouse and human cells by operating as a competing endogenous RNA (ceRNA) for the let-7 microRNA (miRNA) family. Hmga2 can promote the transformation of lung cancer cells independent of protein-coding function but dependent upon the presence of let-7 sites; this occurs without changes in the levels of let-7 isoforms, suggesting that Hmga2 affects let-7 activity by altering miRNA targeting. These effects are also observed in vivo, where Hmga2 ceRNA activity drives lung cancer growth, invasion and dissemination. Integrated analysis of miRNA target prediction algorithms and metastatic lung cancer gene expression data reveals the TGF-β co-receptor Tgfbr3 (ref. 12) as a putative target of Hmga2 ceRNA function. Tgfbr3 expression is regulated by the Hmga2 ceRNA through differential recruitment to Argonaute 2 (Ago2), and TGF-β signalling driven by Tgfbr3 is important for Hmga2 to promote lung cancer progression. Finally, analysis of NSCLC-patient gene-expression data reveals that HMGA2 and TGFBR3 are coordinately regulated in NSCLC-patient material, a vital corollary to ceRNA function. Taken together, these results suggest that Hmga2 promotes lung carcinogenesis both as a protein-coding gene and as a non-coding RNA; such dual-function regulation of gene-expression networks reflects a novel means by which oncogenes promote disease progression.
| S-EPMC3886898 | biostudies-literature
Project description:Cutaneous melanoma (CM) is the most malignant tumor of skin cancers because of its rapid development and high mortality rate. Long noncoding RNAs (lncRNAs), which play essential roles in the tumorigenesis and metastasis of CM and interplay with microRNAs (miRNAs) and mRNAs, are hopefully considered to be efficient biomarkers to detect deterioration during the progression of CM to improve the prognosis. Bioinformatics analysis was fully applied to predict the vital lncRNAs and the associated miRNAs and mRNAs, which eventually constructed the competing endogenous RNA (ceRNA) network to explain the RNA expression patterns in the progression of CM. Further statistical analysis emphasized the importance of these key genes, which were statistically significantly related to one or few clinical features from the ceRNA network. The results showed the lncRNAs MGC12926 and LINC00937 were verified to be strongly connected with the prognosis of CM patients.
| S-EPMC5911588 | biostudies-literature
Project description:Background: Genetic characteristics remain underutilized for establishing prognostic models for colorectal cancer (CRC). We explored the underlying regulatory mechanisms of circular RNAs (circRNAs) that act as competing endogenous RNAs (ceRNAs) and constructed a gene-based nomogram to predict overall survival (OS) in patients with CRC. Methods: We obtained circRNA expression profiling data from the Gene Expression Omnibus (GEO) database. MicroRNA (miRNA) and mRNA expression profiles, with associated clinical data, were obtained from The Cancer Genome Atlas (TCGA). A ceRNA network was established using Cytoscape. Interactions between differential genes were analyzed, and hub genes were identified using the cytoHubba application. The R package "clusterProfiler" was used to evaluate the Gene Ontology (GO) annotations of the differentially expressed mRNAs and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Database-extracted patients were randomized into a training and validation cohorts. A prognostic model was developed using the training set based on multivariate Cox analyses and was then assessed in the validation set. The accuracy of the model was evaluated using discrimination and calibration plots. Results: Thirteen circRNAs, 62 miRNAs, and 301 mRNAs were used to construct the ceRNA network; 10 hub genes were identified via the PPI network. Next, a circRNA- miRNA hub of gene-regulatory modules was established based on four differentially expressed circRNAs, eight differentially expressed miRNAs, and nine differentially expressed mRNAs (DEmRNAs). GO and KEGG pathway analyses indicated the possible association of DEmRNAs with CRC onset and progression. Multivariate analyses revealed that age, tumor stage, and CXCR5 expression were independent risk factors for OS. A CXCR5-based model was developed to predict the OS of patients with CRC in our training set. Our nomogram showed relatively good accuracy, with C-indices of 0.757 and 0.702 in the training and validation sets, respectively. The areas under the curve of the nomograms predicting 3- and 5-years OS were 0.749 and 0.805 in the training set and 0.706 and 0.779 in the validation set, respectively. Conclusions: Our data suggested that the hsa_circ_00001666/has-mir-1229/CXCR5 axis plays an important role in the pathogenesis of CRC, thereby identifying a potential therapeutic target. The proposed CXCR5-based nomogram may also assist surgeons in devising personalized treatments for patients with this disease.
| S-EPMC6857072 | biostudies-literature
Project description:Purpose:This study was aimed to develop a lncRNA-associated competing endogenous RNA (ceRNA) network to provide further understanding of the ceRNA regulatory mechanism and pathogenesis in colorectal cancer (CRC). Patients and methods:Expression profiles of mRNAs, lncRNAs, and miRNAs, and clinical information for CRC patients were obtained from The Cancer Genome Atlas. The differentially expressed mRNAs, lncRNAs, and miRNAs (referred to as "DEmRNAs", "DElncRNAs", and "DEmiRNAs", respectively) were screened out between 539 CRC samples and 11 normal samples. The interactions between DElncRNAs and DEmiRNAs were predicted by miRcode. The DEmRNAs targeted by the DEmiRNAs were retrieved according to TargetScan, miRTar-Base, and miRDB. The lncRNA-miRNA-mRNA ceRNA network was constructed based on the DEmiRNA-DElncRNA and DEmiRNA-DEmRNA interactions. Functional enrichment analysis revealed the biological processes and pathways of DEmRNAs involved in the development of CRC. Key lncRNAs were further analyzed for their associations with overall survival and clinical features of CRC patients. Results:A total of 1,767 DEmRNAs, 608 DElncRNAs, and 283 DEmiRNAs were identified as CRC-specific RNAs. Three hundred eighty-two DEmiRNA-DElncRNA interactions and 68 DEmiRNA-DEmRNA interactions were recognized according to the relevant databases. The lncRNA-miRNA-mRNA ceRNA network was constructed using 25 DEmiRNAs, 52 DEmRNAs, and 64 DElncRNAs. Two DElncRNAs, five DEmiRNAs, and six DEmRNAs were demonstrated to be related to the prognosis of CRC patients. Four DElncRNAs were found to be associated with clinical features. Twenty-eight Gene Ontology terms and 10 Kyoto Encyclopedia of Genes and Genomes pathways were found to be significantly enriched by the DEmRNAs in the ceRNA network. Conclusion:Our results showed cancer-specific mRNA, lncRNA, and miRNA expression patterns and enabled us to construct an lncRNA-associated ceRNA network that provided new insights into the molecular mechanisms of CRC. Key RNA transcripts related to the overall survival and clinical features were also found with promising potential as biomarkers for diagnosis, survival prediction, and classification of CRC.
| S-EPMC5937496 | 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: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