Identification of hub genes and construction of an mRNA-miRNA-lncRNA network of gastric carcinoma using integrated bioinformatics analysis.
Ontology highlight
ABSTRACT: Background
Gastric carcinoma (GC) is one of the most common cancer globally. Despite its worldwide decline in incidence and mortality over the past decades, gastric cancer still has a poor prognosis. However, the key regulators driving this process and their exact mechanisms have not been thoroughly studied. This study aimed to identify hub genes to improve the prognostic prediction of GC and construct a messenger RNA-microRNA-long non-coding RNA(mRNA-miRNA-lncRNA) regulatory network.Methods
The GSE66229 dataset, from the Gene Expression Omnibus (GEO) database, and The Cancer Genome Atlas (TCGA) database were used for the bioinformatic analysis. Differential gene expression analysis methods and Weighted Gene Co-expression Network Analysis (WGCNA) were used to identify a common set of differentially co-expressed genes in GC. The genes were validated using samples from TCGA database and further validation using the online tools GEPIA database and Kaplan-Meier(KM) plotter database. Gene set enrichment analysis(GSEA) was used to identify hub genes related to signaling pathways in GC. The RNAInter database and Cytoscape software were used to construct an mRNA-miRNA-lncRNA network.Results
A total of 12 genes were identified as the common set of differentially co-expressed genes in GC. After verification of these genes, 3 hub genes, namely CTHRC1, FNDC1, and INHBA, were found to be upregulated in tumor and associated with poor GC patient survival. In addition, an mRNA-miRNA-lncRNA regulatory network was established, which included 12 lncRNAs, 5 miRNAs, and the 3 hub genes.Conclusions
In summary, the identification of these hub genes and the establishment of the mRNA-miRNA-lncRNA regulatory network provide new insights into the underlying mechanisms of gastric carcinogenesis. In addition, the identified hub genes, CTHRC1, FNDC1, and INHBA, may serve as novel prognostic biomarkers and therapeutic targets.
SUBMITTER: Wei G
PROVIDER: S-EPMC8718005 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Liver hepatocellular carcinoma (LIHC) is one of the major causes of cancer-related death worldwide with increasing incidences, however there are very few studies about the underlying mechanisms and pathways in the development of LIHC. We obtained LIHC samples from The Cancer Genome Atlas (TCGA) to screen differentially expressed mRNAs, lncRNAs, miRNAs and driver mutations. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Gene ontology enrichment analyses and protein-protein interaction (PPI) network were performed. Moreover, we constructed a competing endogenous lncRNAs-miRNAs-mRNAs network. Finally, cox proportional hazards regression analysis was used to identify important prognostic differentially expressed genes. Total of 1284 mRNAs, 123 lncRNAs, 47 miRNAs were identified within different tissues of LIHC patients. GO analysis indicated that upregulated and downregulated differentially expressed mRNAs (DEmRNAs) were mainly associated with cell division, DNA replication, mitotic sister chromatid segregation and complement activation respectively. Meanwhile, KEGG terms revealed that upregulated and downregulated DEmRNAs were primarily involved in DNA replication, Metabolic pathways, cell cycle and Metabolic pathways, chemical carcinogenesis, retinol metabolism pathway respectively. Among the DERNAs, 542 lncRNAs-miRNAs-mRNAs pairs were predicted to construct a ceRNA regulatory network including 35 DElncRNAs, 26 DEmiRNAs and 112 DEmRNAs. In the Kaplan-Meier analysis, total of 43 mRNAs, 14 lncRNAs and 3 miRNAs were screened out to be significantly correlated with overall survival of LIHC. The mutation signatures were analyzed and its correlation with immune infiltrates were evaluated using the TIMER in LIHC. Among the mutation genes, TTN mutation is often associated with poor immune infiltration and a worse prognosis in LIHC. This work conducted a novel lncRNAs-miRNAs-mRNAs network and mutation signatures for finding potential molecular mechanisms underlying the development of LIHC. The biomarkers also can be used for predicting prognosis of LIHC.
| S-EPMC8051809 | biostudies-literature
Project description:The aim of the present study was to investigate the long non-coding RNA (lncRNA)-microRNA (miRNA)-mRNA regulatory network in gastric cancer (GC) using bioinformatics analysis. Two mRNA gene expression profiles, GSE79973 and GSE54129, and two miRNA expression profiles, GSE93415 and GSE78091, were downloaded from the Gene Expression Omnibus database. The differentially expressed mRNAs (DEMs) and the differentially expressed miRNAs (DEMis) were merged separately. Gene ontology and pathway enrichment analysis were conducted using the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network was then constructed and the 10 top hub genes in the network were analyzed using the Search Tool for the Retrieval of Interacting Genes. The lncRNA-miRNA-mRNA networks were visualized using Cytoscape software. As a result, 158 shared DEMs (40 upregulated and 118 downregulated) were identified from two mRNA datasets. A total of 30 upregulated miRNAs and 1 downregulated miRNA functioned as DEMis. The PPI network consisted of 129 nodes and 572 interactions. The 10 top hub genes were selected by degree using Cytohubba, including Jun proto-oncogene, mitogen-activated protein kinase (MAPK)3, transforming growth factor-β1, Fos proto-oncogene, AP-1 transcription factor subunit, interleukin (IL)-8, MAPK1, RELA proto-oncogene nuclear factor-κB subunit, interferon regulatory factor 7, ubiquitin like modifier and vascular endothelial growth factor A. In the lncRNA-miRNA-mRNA network, a total of 1,215 regulatory associations were constructed using Cytoscape. In conclusion, the present study provides a novel perspective of the molecular mechanisms underlying GC by identifying the lncRNA-miRNA-mRNA regulatory network via bioinformatics analysis.
| S-EPMC6865131 | biostudies-literature
Project description:BackgroundNumerous studies have highlighted that long non-coding RNAs (lncRNAs) can bind to microRNA (miRNA) sites as competing endogenous RNAs (ceRNAs), thereby affecting and regulating the expression of mRNAs and target genes. These lncRNA-associated ceRNAs have been theorized to play a significant role in cancer initiation and progression. However, the roles and functions of the lncRNA-miRNA-mRNA ceRNA network in squamous cell carcinoma of the tongue (SCCT) are still unclear.MethodsThe miRNA, mRNA and lncRNA expression profiles from 138 patients with SCCT were downloaded from The Cancer Genome Atlas database. We identified the differential expression of miRNAs, mRNAs, and lncRNAs using the limma package of R software. We used the clusterProfiler package for GO and KEGG pathway annotations. The survival package was used to estimate survival analysis according to the Kaplan-Meier curve. Finally, the GDCRNATools package was used to construct the lncRNA-miRNA-mRNA ceRNA network.ResultsIn total, 1943 SCCT-specific mRNAs, 107 lncRNAs and 100 miRNAs were explored. Ten mRNAs (CSRP2, CKS2, ADGRG6, MB21D1, GMNN, RIPOR3, RAD51, PCLAF, ORC1, NAGS), 9 lncRNAs (LINC02560, HOXC13?-?AS, FOXD2?-?AS1, AC105277.1, AC099850.3, STARD4?-?AS1, SLC16A1?-?AS1, MIR503HG, MIR100HG) and 8 miRNAs (miR?-?654, miR?-?503, miR?-?450a, miR?-?379, miR?-?369, miR?-?190a, miR?-?101, and let-7c) were found to be significantly associated with overall survival (log-rank p?<?0.05). Based on the analysis of the lncRNA-miRNA-mRNA ceRNA network, one differentially expressed (DE) lncRNA, five DEmiRNAs, and three DEmRNAs were demonstrated to be related to the pathogenesis of SCCT.ConclusionsIn this study, we described the gene regulation by the lncRNA-miRNA-mRNA ceRNA network in the progression of SCCT. We propose a new lncRNA-associated ceRNA that could help in the diagnosis and treatment of SCCT.
| S-EPMC6686570 | biostudies-literature
Project description:BACKGROUND As all we know, gastric cancer (GC) is a highly aggressive disease. Recently, circular RNA (circRNA) was found to play a vital role in regulation of GC. Some circRNAs could regulate messenger RNA (mRNA) expression by functioning as a microRNA (miRNA) sponge. Nevertheless, the circRNA-miRNA-mRNA regulatory network involved GC rarely has been explored and researched. MATERIAL AND METHODS All the differentially expressed circRNAs, miRNAs, and mRNA were derived from Gene Expression Omnibus (GEO) microarray data (GSE78092, GSE89143, GSE93415, and GSE54129). GC level 3 miRNA-sequencing data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Furthermore, a circRNA-miRNA-mRNA regulatory network was constructed by Cytoscape (version 3.6.1). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway revealed the functions and signaling pathways associated with these target genes. Hub genes of protein-protein interaction (PPI) network were identified by STRING database and cytoHubba. RESULTS The regulatory network consists of 3 circRNAs, 22 miRNAs, and 128 mRNAs. Only 3 miRNAs of the network were consistent with the expression of TCGA and were associated with some clinical features. The results of the functional analysis of 128 mRNAs showed that GO analysis and KEGG pathways of inclusion criteria were 49 and 24, respectively. PPI network and Cytoscape showed that the top 10 hub genes were MYC, CTGF, TGFBR2, TGFBR1, SERPINE1, KRAS, ZEB1, THBS1, CDK6, and TNS1; 4 of which were verified by GEPIA based on TCGA. Highly expressed SERPINE1 had a poor OS (over survival) and DFS (disease-free survival), and TGFBR1 expression increased along with the increase of clinical stages. CONCLUSIONS This study looked at a circRNA-miRNA-mRNA regulatory network associated with GC and explored the potential functions of mRNA in the network, then identified a new molecular marker for prediction, prognosis, and therapeutic targets for clinical patients.
| S-EPMC6880644 | biostudies-literature
Project description:Background:Chronic obstructive pulmonary disease (COPD) has become a major cause of morbidity and mortality worldwide. Increasing evidence indicates that aberrantly expressed microRNAs (miRNAs) are involved in the pathogenesis of COPD. However, an integrative exploration of miRNA-mRNA regulatory network in COPD plasma remains lacking. Methods:The microarray datasets GSE24709, GSE61741, and GSE31568 were downloaded from the GEO database and analyzed using GEO2R tool to identify differentially expressed miRNAs (DEMs) between COPD and normal plasma. The consistently changing miRNAs in the three datasets were screened out as candidate DEMs. Potential upstream transcription factors and downstream target genes of candidate DEMs were predicted by FunRich and miRNet, respectively. Next, GO annotation and KEGG pathway enrichment analysis for target genes were performed using DAVID. Then, PPI and DEM-hub gene network were constructed using the STRING database and Cytoscape software. Finally, GSE56768 was used to evaluate the hub gene expressions. Results:A total of nine (six upregulated and three downregulated) DEMs were screened out in the above three datasets. SP1 was predicted to potentially regulate most of the downregulated DEMs, while YY1 and E2F1 could regulate both upregulated and downregulated DEMs. 1139 target genes were then predicted, including 596 upregulated DEM target genes and 543 downregulated DEM target genes. Target genes of DEMs were mainly enriched in PI3K/Akt signaling pathway, mTOR signaling pathway, and autophagy. Through the DEM-hub gene network construction, most of the hub genes were found to be potentially modulated by miR-497-5p, miR-130b-5p, and miR-126-5p. Among the top 12 hub genes, MYC and FOXO1 expressions were consistent with that in the GSE56768 dataset. Conclusion:In the study, potential miRNA-mRNA regulatory network was firstly constructed in COPD plasma, which may provide a new insight into the pathogenesis and treatment of COPD.
| S-EPMC7490070 | biostudies-literature
Project description:BackgroundMicroRNA (miRNA) has been confirmed to be involved in the occurrence, development, and prevention of diabetic nephropathy (DN), but its mechanism of action is still unclear.ObjectiveWith the help of the GEO database, bioinformatics methods are used to explore the miRNA-mRNA regulatory relationship pairs related to diabetic nephropathy and explain their potential mechanisms of action.MethodsThe DN-related miRNA microarray dataset (GSE51674) and mRNA expression dataset (GSE30122) are downloaded through the GEO database, online analysis tool GEO2R is used for data differential expression analysis, TargetScan, miRTarBase, and miRDB databases are used to predict potential downstream target genes regulated by differentially expressed miRNAs, and intersection with differential genes is used to obtain candidate target genes. According to the regulatory relationship between miRNA and mRNA, the miRNA-mRNA relationship pair is clarified, and the miRNA-mRNA regulatory network is constructed using Cytoscape. DAVID is used to perform GO function enrichment analysis and KEGG pathway analysis of candidate target genes. By GeneMANIA prediction of miRNA target genes and coexpressed genes, the protein interaction network is constructed. Results and Conclusions. A total of 67 differentially expressed miRNAs were screened in the experiment, of which 42 were upregulated and 25 were downregulated; a total of 448 differentially expressed mRNAs were screened, of which 93 were upregulated and 355 were downregulated. Using TargetScan, miRTarBase, and miRDB databases to predict downstream targets of differentially expressed miRNAs, 2283 downstream target genes coexisting in 3 databases were predicted to intersect with differentially expressed mRNAs to obtain 96 candidate target genes. Finally, 44 miRNA-mRNA relationship pairs consisting of 12 differentially expressed miRNAs and 27 differentially expressed mRNAs were screened out; further analysis showed that miRNA regulatory network genes may participate in the occurrence and development of diabetic nephropathy through PI3K/Akt, ECM-receptor interaction pathway, and RAS signaling pathway.
| S-EPMC8616647 | biostudies-literature
Project description:The aim of this study was to establish a novel competing endogenous RNA (ceRNA) network able to predict prognosis in patients with triple-negative breast cancer (TNBC). Differential gene expression analysis was performed using the GEO2R tool. Enrichr and STRING were used to conduct protein-protein interaction and pathway enrichment analyses, respectively. Upstream lncRNAs and miRNAs were identified using miRNet and mirTarBase, respectively. Prognostic values, expression, and correlational relationships of mRNAs, lncRNAs, and miRNAs were examined using GEPIA, starBase, and Kaplan-Meier plotter. It total, 860 upregulated and 622 downregulated differentially expressed mRNAs were identified in TNBC. Ten overexpressed and two underexpressed hub genes were screened. Next, 10 key miRNAs upstream of these key hub genes were predicted, of which six upregulated miRNAs were significantly associated with poor prognosis and four downregulated miRNAs were associated with good prognosis in TNBC. NEAT1 and MAL2 were selected as key lncRNAs. An mRNA-miRNA-lncRNA network in TNBC was constructed. Thus, we successfully established a novel mRNA-miRNA-lncRNA regulatory network, each component of which is prognostic for TNBC.
| S-EPMC7835059 | biostudies-literature
Project description:BackgroundMany studies on long chain non-coding RNAs (lncRNAs) are published in recent years. But the roles of lncRNAs in aortic dissection (AD) are still unclear and should be further examined. The present work focused on determining the molecular mechanisms underlying lncRNAs regulation in aortic dissection on the basis of the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network.MethodsThis study collected the lncRNAs (GSE52093), mRNAs (GSE52093) and miRNAs (GSE92427) expression data within human tissue samples with aortic dissection group and normal group based on Gene Expression Omnibus (GEO) database.ResultsThis study identified three differentially expressed lncRNAs (DELs), 19 differentially expressed miRNAs (DEmiRs) and 1046 differentially expressed mRNAs (DEGs) identified regarding aortic dissection. Furthermore, we constructed a lncRNA-miRNA-mRNA network through three lncRNAs (including two with up-regulation and one with down-regulation), five miRNAs (five with up-regulation), as well as 211 mRNAs (including 103 with up-regulation and 108 with down-regulation). Simultaneously, we conducted functional enrichment and pathway analyses on genes within the as-constructed ceRNA network. According to our PPI/ceRNA network and functional enrichment analysis results, four critical genes were found (E2F2, IGF1R, BDNF and PPP2R1B). In addition, E2F2 level was possibly modulated via lncRNA FAM87A-hsa-miR-31-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p. The expression of IGF1R may be regulated by lncRNA FAM87A-hsa-miR-16-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p.ConclusionIn conclusion, the ceRNA interaction axis we identified is a potentially critical target for treating AD. Our results shed more lights on the possible pathogenic mechanism in AD using a lncRNA-associated ceRNA network.
| S-EPMC8495997 | biostudies-literature
Project description:Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA-mRNA and lncRNA-mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein-protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.
| S-EPMC9332185 | biostudies-literature
Project description:Gastric cancer (GC), as an epidemic cancer worldwide, has more than 1 million new cases and an estimated 769,000 deaths worldwide in 2020, ranking fifth and fourth in global morbidity and mortality. In mammals, both miRNAs and transcription factors (TFs) play a partial role in gene expression regulation. The mRNA expression profile and miRNA expression profile of GEO database were screened by GEO2R for differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs). Then, DAVID annotated the functions of DEGs to understand the functions played in biological processes. The prediction of potential target genes of miRNA and key TFs of mRNA was performed by mipathDB V2.0 and CHEA3, respectively, and the gene list comparison was performed to look for overlapping genes coregulated by key TFs and DEMs. Finally, the obtained miRNAs, TF, and overlapping genes were used to construct the miRNA-mRNA-TF regulatory network, which was verified by RT-qPCR. 76 upregulated DEGs, 199 downregulated DEGs, and 3 upregulated miRNAs (miR-199a-3p/miR-199b-3p, miR-125b-5p, and miR-199a-5p) were identified from the expression profiles of mRNA (GSE26899, GSE29998, GSE51575, and GSE13911) and miRNA (GSE93415), respectively. Through database prediction and gene list comparison, it was found that among the 199 downregulated DEGs, 61, 71, and 69 genes were the potential targets of miR-199a-3p/miR-199b-3p, miR-125b-5p, and miR-199a-5p, respectively. 199 downregulated DEGs were used as the gene list for the prediction of key TFs, and the results showed that RFX6 ranked the highest. The potential target overlap genes of miR-199a-3p/miR-199b-3p, miR-125b-5p, and miR-199a-5p were 4 genes (SH3GL2, ATP4B, CTSE, and SORBS2), 7 genes (SLC7A8, RNASE4, ESRRG, PGC, MUC6, Fam3B, and FMO5), and 6 genes (CHGA, PDK4, TMPRSS2, CLIC6, GPX3, and PSCA), respectively. Finally, we constructed a miRNA-mRNA-TF regulatory network based on the above 17 mRNAs, 3 miRNAs, and 1 TF and verified by RT-qPCR and western blot results that the expression of RFX6 was downregulated in GC tissues. These identified miRNAs, mRNAs, and TF have a certain reference value for further exploration of the regulatory mechanism of GC.
| S-EPMC8616677 | biostudies-literature