Integrated Analysis of ceRNA Regulatory Network Associated With Tumor Stage in Cervical Cancer.
Ontology highlight
ABSTRACT: Objective
To analyze the abnormally expressed genes involved in cervical cancer occurrence and development.Materials and methods
Integrated bioinformatics methods were used to analyze differentially expressed (DE) RNAs, including mRNAs, microRNAs (miRNAs), and long non-coding RNAs (lncRNAs), in stage I, II, III, and IV cervical cancer patients from the TCGA database to fully reveal the dynamic changes caused by cervical cancer.Results
First, DE RNAs in cervical cancer tissues from stage I, II, III, and IV patients and normal cervical tissues were identified and divided into different profiles. Several DE RNA profiles were down-regulated or up-regulated in stage I, III, and IV patients. GO and KEGG analysis of DE mRNA profile 1, 2, 4, 5, 6 and 22 which were significantly down-regulated or up-regulated showed that DE mRNAs are involved in cell division, DNA replication, cell adhesion, the positive and negative regulation of RNA polymerase ll promoter transcription. Besides, DE RNA profiles with significant differences in patient stages were analyzed to perform a competing endogenous RNA (ceRNA) regulatory network of lncRNA, miRNA, and mRNA. The protein-protein interaction (PPI) network of DE mRNAs in the ceRNA regulatory network was also constructed. The network had nine central genes (up-regulated genes: CDKN2A, GSK3B, BIRC5, CYCS, MAD2L1; down-regulated genes: PTEN, FOXO3, CCND2, TGFBR2). Survival analysis found that 5 lncRNAs, 9 mRNAs, and 4 miRNAs can be used as prognostic indicators of cervical cancer. Finally, combined with cluster analysis results, we further screened 2 DE RNAs (AMZ2P1 and HDAC5) using clinical samples, suggesting that AMZ2P1, and HDAC5 may act as diagnostic biomarkers for the development of cervical cancer.Conclusion
This research provides new effective targets and reliable biological markers for the diagnosis and prognosis of cervical cancer.
SUBMITTER: Ma X
PROVIDER: S-EPMC8021857 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:This study aimed to identify potential therapeutic targets in osteosarcoma (OS) through the network analysis of competing endogenous RNAs (ceRNAs). The differentially expressed miRNAs (DEMIs) and mRNAs (DEMs) were identified between OS cell lines and human mesenchymal stem cells (hMSCs) from the data deposited under GSE70415 using limma package. Functional analysis of DEMs was performed using DAVID and clusterProfiler to identify significantly enriched Gene Ontology biological processes and KEGG pathways, respectively. The DEMI-DEM interaction network was constructed using Cytoscape. LncRNA-miRNA interactions were predicted using starBase database. The ceRNA regulatory network was constructed by integrating mRNAs, miRNAs, and lncRNAs, and functional enrichment analysis was performed for the genes involved. The analysis revealed a total of 326 DEMs and 54 DEMIs between OS cells and hMSCs. We identified several novel therapeutic targets involved in the progression and metastasis of OS, such as CBX7, RAD9A, SNHG7 and miR-34a-5p. The miRNA, miR-543 (target gene: CBX7) was found to be associated with the pathway Mucin type O-glycan biosynthesis. Using the ceRNA network, we established the following regulatory interactions: NEAT1/miR-543/CBX7, SNHG7/miR-34a-5p/RAD9A, and XIST/miR-34a-5p/RAD9A. CBX7, RAD9A, lncRNA SNHG7, miR-543, and miR-34a-5p may be explored as novel therapeutic targets for treatment of OS.
| S-EPMC7252519 | biostudies-literature
Project description:BackgroundOur previous studies have confirmed that miR-154-5p can regulate pRb expression, and thus, play a tumor suppressor role in HPV16 E7-induced cervical cancer. However, its upstream molecules have not been elucidated in the progression of cervical cancer. This study aimed to explore the role of the miR-154-5p upstream molecule, hsa_circ_0000276 in cervical cancer development and its possible mechanisms of action.MethodsWe detected differences in whole transcriptome expression profiles of cervical squamous carcinoma and tissues adjacent to cervical cancer tissues from patients using microarray technology to predict circular RNAs (circRNAs) with binding sites to miR-154-5p. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to detect the expression of hsa_circ_0000276 (which had the strongest binding capacity to miR-154 and was selected as the target molecule) in cervical cancer tissues, followed by in vitro functional assays. Downstream microRNAs (miRNAs) and mRNAs of hsa_circ_0000276 were identified using transcriptome microarray data and databases, while the protein-protein interaction networks were obtained using STRING. A competing endogenous RNA (ceRNA) network centered on hsa_circ_0000276 was constructed using Cytoscape and GO and KEGG databases. Abnormal expression and prognosis of critical downstream molecules were analyzed using gene databases and molecular experiments. qRT-PCR and western blot analysis was performed to verify the expression of candidate genes.ResultsWe identified 4,001 differentially expressed circRNAs between HPV16-positive cervical squamous carcinoma and benign cervical tissues and 760 circRNAs targeting miR-154-5p, including hsa_circ_0000276. hsa_circ_0000276 and miR-154-5p directly bound, and hsa_circ_0000276 was upregulated, in cervical precancerous lesions and cervical cancer tissues and cells. Silencing hsa_circ_0000276 inhibited G1/S transition and cell proliferation and promoted apoptosis in SiHa and CaSki cells. Bioinformatics analysis showed that the hsa_circ_0000276 ceRNA network included 17 miRNAs and seven mRNAs, and downstream molecules of hsa_circ_0000276 were upregulated in cervical cancer tissues. These downstream molecules were associated with a poor prognosis and affected cervical cancer-associated immune infiltration. Of these, expression of CD47, LDHA, PDIA3, and SLC16A1 was downregulated in sh_hsa_circ_0000276 cells.ConclusionsOur findings show that hsa_circ_0000276 exerts cancer-promoting effects in cervical cancer and is an underlying biomarker for cervical squamous cell carcinoma.
| S-EPMC9999601 | biostudies-literature
Project description:Esophageal cancer (ESCA) is one of the most commonly diagnosed cancers in the world. Tumor immune microenvironment is closely related to tumor prognosis. The present study aimed at analyzing the competing endogenous RNA (ceRNA) network and tumor-infiltrating immune cells in ESCA. The expression profiles of mRNAs, lncRNAs, and miRNAs were downloaded from the Cancer Genome Atlas database. A ceRNA network was established based on the differentially expressed RNAs by Cytoscape. CIBERSORT was applied to estimate the proportion of immune cells in ESCA. Prognosis-associated genes and immune cells were applied to establish prognostic models basing on Lasso and multivariate Cox analyses. The survival curves were constructed with Kaplan-Meier method. The predictive efficacy of the prognostic models was evaluated by the receiver operating characteristic (ROC) curves. The differentially expressed mRNAs, lncRNAs, and miRNAs were identified. We constructed the ceRNA network including 23 lncRNAs, 19 miRNAs, and 147 mRNAs. Five key molecules (HMGB3, HOXC8, HSPA1B, KLHL15, and RUNX3) were identified from the ceRNA network and five significant immune cells (plasma cells, T cells follicular helper, monocytes, dendritic cells activated, and neutrophils) were selected via CIBERSORT. The ROC curves based on key genes and significant immune cells all showed good sensitivity (AUC of 3-year survival: 0.739, AUC of 5-year survival: 0.899, AUC of 3-year survival: 0.824, AUC of 5-year survival: 0.876). There was certain correlation between five immune cells and five key molecules. The present study provides an effective bioinformatics basis for exploring the potential biomarkers of ESCA and predicting its prognosis.
| S-EPMC8164107 | biostudies-literature
Project description:Introduction: Stroke, of which ischemic stroke (IS) is the major type, is the second leading cause of disability and death worldwide. Circular RNAs (circRNAs) are reported to play important role in the physiology and pathology of IS. CircRNAs often act as competing endogenous RNA (ceRNA) to regulate gene expression by acting as miRNA sponges. However, whole transcriptome-wide screenings of circRNA-mediated ceRNA networks associated with IS are still lacking. In the present study, we constructed a circRNA-miRNA-mRNA ceRNA network by whole transcriptome-wide analysis. Methods: CircRNAs, miRNAs and mRNAs expression profiles were downloaded from the Gene Expression Omnibus (GEO) datasets. We identified differentially expressed (DE) circRNAs, miRNAs, and mRNAs in IS patients. StarBase and CircBank databases were used to predict the miRNA targets of DEcircRNAs, and mirDIP database was used to predict the mRNA targets of DEmiRNAs. CircRNA-miRNA pairs and miRNA-mRNA pairs were established. Then, we identified hub genes via protein-protein interaction analysis and constructed a core ceRNA sub-network. Results: In total, 276 DEcircRNAs, 43 DEmiRNAs, and 1926 DEmRNAs were explored. The ceRNA network included 69 circRNAs, 24 miRNAs, and 92 mRNAs. The core ceRNA subnetwork included hsa_circ_0011474, hsa_circ_0023110, CDKN1A, FHL2, RPS2, CDK19, KAT6A, CBX1, BRD4, and ZFHX3. Discussion: In conclusion, we established a novel hsa_circ_0011474 - hsa-miR-20a-5p/hsa-miR-17-5p - CDKN1A ceRNA regulatory axis associated with IS. Our findings provide new insights into the pathogenesis of IS and offer promising diagnostic and predictive biomarkers.
| S-EPMC10267749 | biostudies-literature
Project description:BackgroundBreast cancer is a malignancy and lethal tumor in women. Metastasis of breast cancer is one of the causes of poor prognosis. Increasing evidences have suggested that the competing endogenous RNAs (ceRNAs) were associated with the metastasis of breast cancer. Nonetheless, potential roles of ceRNAs in regulating the metastasis of breast cancer remain unclear.MethodsThe RNA expression (3 levels) and follow-up data of breast cancer and noncancerous tissue samples were downloaded from the Cancer Genome Atlas (TCGA). Differentially expressed and metastasis associated RNAs were identified for functional analysis and constructing the metastasis associated ceRNA network by comprehensively bioinformatic analysis. The Kaplan-Meier (K-M) survival curve was utilized to screen the prognostic RNAs in metastasis associated ceRNA network. Moreover, we further identified the metastasis associated biomarkers with operating characteristic (ROC) curve. Ultimately, the data of Cancer Cell Line Encyclopedia (CCLE, https://portals.broadinstitute.org/ccle) website were selected to obtained the reliable metastasis associated biomarkers.Results1005 mRNAs, 22 miRNAs and 164 lncRNAs were screened as differentially expressed and metastasis associated RNAs. The results of GO function and KEGG pathway enrichment analysis showed that these RNAs are mainly associated with the metabolic processes and stress responses. Next, a metastasis associated ceRNA (including 104 mRNAs, 19 miRNAs, and 16 lncRNAs) network was established, and 12 RNAs were found to be related to the overall survival (OS) of patients. In addition, 3 RNAs (hsa-miR-105-5p, BCAR1, and PANX2) were identified to serve as reliable metastasis associated biomarkers. Eventually, the results of mechanism analysis suggested that BCAR1 might promote the metastasis of breast cancer by facilitating Rap 1 signaling pathway.ConclusionIn the present research, we identified 3 RNAs (hsa-miR-105-5p, BCAR1 and PANX2) might associated with prognosis and metastasis of breast cancer, which might be provide a new perspective for metastasis of breast cancer and contributed to the treatment of breast cancer.
| S-EPMC8158160 | biostudies-literature
Project description:Circular RNAs (circRNAs) have displayed dysregulated expression in several types of cancer. However, the functions of the majority of circRNAs in pancreatic ductal adenocarcinoma (PDAC) remain unknown. The present study aimed to investigate the expression, functions and molecular mechanisms of circRNAs in PDAC. The circRNAs, mRNAs and the microRNA (miRNAs) expression profiles were obtained from three Gene Expression Omnibus microarray datasets, and a circRNA-miRNA-mRNA and circRNA-miRNA-hubgene network was established. The interactions between proteins were analyzed using the Search Tool for the Retrieval of Interacting Genes/Proteins database, and hubgenes were identified using the MCODE plugin. A total of eight differentially expressed circRNAs (DEcircRNAs), 44 differentially expressed miRNAs (DEmiRNAs), and 2,052 differentially expressed mRNAs (DEmRNAs) were identified. The present study successfully constructed a circRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network based on four circRNAs, six miRNAs and 111 mRNAs in PDAC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathways analyses indicated that DEmRNAs may participate in the pathogenesis and progression of PDAC. The protein-protein interaction network and module analysis identified six hubgenes (THBS1, FN1, TIMP3, TGFB2, ITGA1 and ITGA3). Furthermore, the circRNA-miRNA-hubgene regulatory modules were constructed based on the three DEcircRNAs, one DEmiRNAs and five DEmRNAs. In conclusion, the results of the present study improve the current understanding of the pathogenesis of PDAC.
| S-EPMC7039142 | biostudies-literature
Project description:Circular RNA (circRNA) is closely related to tumorigenesis and cancer progression. Yet, the roles of cancer-specific circRNAs in the circRNA-related ceRNA network of breast cancer (BRCA) remain unclear. The aim of this study was to construct a ceRNA network associated with circRNA and to explore new therapeutic and prognostic targets and biomarkers for breast cancer. We downloaded the circRNA expression profile of BRCA from Gene Expression Omnibus (GEO) microarray datasets and downloaded the miRNA and mRNA expression profiles of BRCA from The Cancer Genome Atlas (TCGA) database. Differentially expressed mRNAs (DEmRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed circRNAs (DEcircRNAs) were identified, and a competitive endogenous RNA (ceRNA) regulatory network was constructed based on circRNA-miRNA pairs and miRNA-mRNA pairs. Gene ontology and pathway enrichment analyses were performed on mRNAs regulated by circRNAs in ceRNA networks. Survival analysis and correlation analysis of all mRNAs and miRNAs in the ceRNA network were performed. A total of 72 DEcircRNAs, 158 DEmiRNAs, and 2762 DE mRNAs were identified. The constructed ceRNA network contains 60 circRNA-miRNA pairs and 140 miRNA-mRNA pairs, including 40 circRNAs, 30 miRNAs, and 100 mRNAs. Functional enrichment indicated that DEmRNAs regulated by DEcircRNAs in ceRNA networks were significantly enriched in the PI3K-Akt signaling pathway, microRNAs in cancer, and proteoglycans in cancer. Survival analysis and correlation analysis of all mRNAs and miRNAs in the ceRNA network showed that 13 mRNAs and 6 miRNAs were significantly associated with overall survival, and 48 miRNA-mRNA interaction pairs had a significant negative correlation. A PPI network was established, and 21 hub genes were determined from the network. This study provides an effective bioinformatics basis for further understanding of the molecular mechanisms and predictions of breast cancer. A better understanding of the circRNA-related ceRNA network in BRCA will help identify potential biomarkers for diagnosis and prognosis.
| S-EPMC8712159 | biostudies-literature
Project description:Increasing evidence has experimentally proved the competitive endogenous RNA (ceRNA) hypothesis that long non-coding RNA (lncRNA) can affect the expression of RNA targets by competitively combining microRNA (miRNA) via miRNA response elements. However, an extensive ceRNA network of thyroid carcinoma in a large cohort has not been evaluated. We analyzed the RNAseq and miRNAseq data of 348 cases of primary papillary thyroid cancer (PTC) patients with clinical information downloaded from The Cancer Genome Atlas (TCGA) project to search for potential biomarkers or therapeutic targets. A computational approach was applied to build an lncRNA-miRNA-mRNA regulatory network of PTC. In total, 780 lncRNAs were detected as collectively dysregulated lncRNAs in all 3 PTC variants compared with normal tissues (fold change >2 and false discovery rate <0.05). The interactions among 45 lncRNAs, 13 miRNAs and 86 mRNAs constituted a ceRNA network of PTC. Nine out of the 45 aberrantly expressed lncRNAs were related to the clinical features of PTC patients. However, the expression levels of 3 lncRNAs (LINC00284, RBMS3-AS1 and ZFX-AS1) were identified to be tightly correlated with the patients overall survival (log-rank, P<0.05). The present study identified a list of specific lncRNAs associated with PTC progression and prognosis. This complex ceRNA interaction network in PTC may provide guidance for better understanding the molecular mechanisms underlying PTC.
| S-EPMC5802034 | biostudies-literature
Project description:Long noncoding RNAs (lncRNAs) have important biological functions as competing endogenous RNAs (ceRNAs) in tumors, yet the functions and regulatory mechanisms of lncRNA-related ceRNAs in gastric cancer have not been fully elucidated. In this study, we constructed a lncRNA-miRNA-mRNA ceRNA network and identified potential lncRNA biomarkers in gastric cancer. Basing on the RNA profiles downloaded from The Cancer Genome Atlas (TCGA) platform, the gastric cancer-specific differentially expressed lncRNAs, miRNAs, and mRNAs were screened for constructing a ceRNA network using bioinformatic tools. The enrichment analysis of the biological processes in Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathways was performed on the ceRNA-related DEmRNAs. According to the modularization of protein-protein interaction (PPI) network, we extracted a ceRNA subnetwork and analyzed the correlation between the expression of the lncRNAs involved and specific clinical features of patients. Next, the expression of highly up-regulated in liver cancer (HULC) and RP11-314B1.2 showed significant changes in several pathological processes involved in gastric cancer, and nine lncRNAs were found to be correlated with the overall survival of patients with gastric cancer. Through the univariate and multivariate Cox regression analyses, two lncRNAs (LINC00106 and RP11-999E24.3) were identified and utilized to establish a risk score model for assessing the prognosis of patients. The analysis results were also partially verified using quantitative real-time PCR. The findings from this study indicate that HULC, RP11-314B1.2, LINC00106, and RP11-999E24.3 could be considered as potential therapeutic targets or prognostic biomarkers in gastric cancer, and provide a new perspective for cancer pathogenesis research.
| S-EPMC7050084 | biostudies-literature
Project description:BackgroundGastric cancer (GC) is a malignant tumour that originates in the gastric mucosal epithelium and is associated with high mortality rates worldwide. Long noncoding RNAs (lncRNAs) have been identified to play an important role in the development of various tumours, including GC. Yet, lncRNA biomarkers in a competing endogenous RNA network (ceRNA network) that are used to predict survival prognosis remain lacking. The aim of this study was to construct a ceRNA network and identify the lncRNA signature as prognostic factors for survival prediction.MethodsThe lncRNAs with overall survival significance were used to construct the ceRNA network. Function enrichment, protein-protein interaction, and cluster analysis were performed for dysregulated mRNAs. Multivariate Cox proportional hazards regression was performed to screen the potential prognostic lncRNAs. RT-qPCR was used to measure the relative expression levels of lncRNAs in cell lines. CCK8 assay was used to assess the proliferation of GC cells transfected with sh-lncRNAs.ResultsDifferentially expressed genes were identified including 585 lncRNAs, 144 miRNAs, and 2794 mRNAs. The ceRNA network was constructed using 35 DElncRNAs associated with overall survival of GC patients. Functional analysis revealed that these dysregulated mRNAs were enriched in cancer-related pathways, including TGF-beta, Rap 1, calcium, and the cGMP-PKG signalling pathway. A multivariate Cox regression analysis and cumulative risk score suggested that two of those lncRNAs (LINC01644 and LINC01697) had significant prognostic value. Furthermore, the results indicate that LINC01644 and LINC01697 were upregulated in GC cells. Knockdown of LINC01644 or LINC01697 suppressed the proliferation of GC cells.ConclusionsThe authors identified 2-lncRNA signature in ceRNA regulatory network as prognostic biomarkers for the prediction of GC patient survival and revealed that silencing LINC01644 or LINC01697 inhibited the proliferation of GC cells.
| S-EPMC8479203 | biostudies-literature