Dataset Information

Identification of hub genes and construction of an mRNA-miRNA-lncRNA network of gastric carcinoma using integrated bioinformatics analysis.

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

ACCESS DATA

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.

Project description:Long non-coding RNAs (lncRNAs) are ncRNA transcripts >200 nucleotides that are important genetic regulators. LncRNAs can directly regulate mRNA through a lncRNA-mRNA regulatory mode and can also regulate mRNA through competitive binding to micro (mi)RNA, which is generally known as the competitive endogenous RNA (ceRNA) network. The present study evaluated the functional roles and regulatory networks of lncRNAs in chronic glomerulonephritis (CGN). The proliferative ability of mouse glomerular mesangial cells (GMCs) induced by different concentrations of lipopolysaccharide (LPS) was assessed using the Cell Counting Kit-8 assay, and RNA sequencing (RNA-seq) was performed to identify differentially expressed lncRNAs in LPS-induced GMCs. Based on the sequencing results, six lncRNAs were selected for validation using reverse transcription-quantitative PCR (RT-qPCR). Furthermore, the lncRNA-mRNA regulatory network and the lncRNA-miRNA-mRNA ceRNA network were constructed to assess the role and mechanism of CGN-related lncRNAs. To elucidate the biological functions of lncRNAs, Gene Ontology (GO) biological process term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on all mRNAs involved in the lncRNA-mRNA regulatory network and in the ceRNA network. A total of 1,532 differentially expressed lncRNAs, including 594 upregulated lncRNAs and 938 downregulated lncRNAs, were identified using RNA-seq. The results of RT-qPCR validation were consistent with RNA-seq results. An lncRNA-mRNA regulatory network, including 236 lncRNAs and 556 mRNAs, and a ceRNA network, including 6 lncRNAs, 18 miRNAs and 419 mRNAs, were successfully constructed. The GO biological process term enrichment and KEGG pathway enrichment analyses demonstrated that those lncRNAs were often related to inflammatory response and substance metabolism. The present study identified key CGN-related lncRNAs in LPS-induced GMCs, and further demonstrated a global view of the lncRNA-mRNA regulatory network and ceRNA network involved in CGN. These results offered novel insights into the roles of lncRNAs in the pathogenesis of CGN and identified potential diagnostic biomarkers.

Project description:Asthma is a chronic respiratory disease caused by environment-host interactions. Bronchial epithelial cells (BECs) are the first line of defense against environmental toxins. However, the mechanisms underlying the role of BECs in severe asthma (SA) are not yet fully understood. Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in the regulation of gene expression in the pathogenesis of SA. In this study, bioinformatics was used for the first time to reveal the lncRNA-miRNA-mRNA regulatory network of BECs in SA. Five mRNA datasets of bronchial brushing samples from patients with SA and healthy controls (HC) were downloaded from the Gene Expression Omnibus (GEO) database. A combination of the Venn diagram and robust rank aggregation (RRA) method was used to identify core differentially expressed genes (DEGs). Protein-protein interaction (PPI) analysis of core DEGs was performed to screen hub genes. The miRDB, miRWalk, and ENCORI databases were used to predict the miRNA-mRNA relationships, and the ENCORI and starBase v2.0 databases were used to predict the upstream lncRNAs of the miRNA-mRNA relationships. Four core DEGs were identified: carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5), interleukin-1 receptor type 2 (IL1R2), trefoil factor 3 (TFF3), and vascular endothelial growth factor A (VEGFA). These 4 core DEGs indicated that SA was not significantly associated with sex. Enrichment analysis showed that the MAPK, Rap1, Ras, PI3K-Akt and Calcium signaling pathways may serve as the principal pathways of BECs in SA. A lncRNA-miRNA-mRNA regulatory network of the severe asthmatic bronchial epithelium was constructed. The top 10 competing endogenous RNAs (ceRNAs) were FGD5 antisense RNA 1 (FGD5-AS1), metastasis associated lung adenocarcinoma transcript 1 (MALAT1), X inactive specific transcript (XIST), HLA complex group 18 (HCG18), small nucleolar RNA host gene 16 (SNHG16), has-miR-20b-5p, has-miR-106a-5p, hsa-miR-106b-5p, has-miR-519d-3p and Fms related receptor tyrosine kinase 1 (FLT1). Our study revealed a potential mechanism for the lncRNA-miRNA-mRNA regulatory network in BECs in SA.

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.

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.

Project description:Background: Acute myocardial infarction (AMI) is one of the main fatal diseases of cardiovascular diseases. Circular RNA (circRNA) is a non-coding RNA (ncRNA), which plays a role in cardiovascular disease as a competitive endogenous RNA (ceRNA). However, their role in AMI has not been fully clarified. This study aims to explore the mechanism of circRNA-related ceRNA network in AMI, and to identify the corresponding immune infiltration characteristics. Materials and Methods: The circRNA (GSE160717), miRNA (GSE24548), and mRNA (GSE60993) microarray datasets of AMI were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed circRNAs (DEcircRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were identified by the "limma" package. After integrating the circRNA, miRNA and mRNA interaction, we constructed a circRNA-miRNA-mRNA network. The "clusterProfiler" package and String database were used for functional enrichment analysis and protein-protein interaction (PPI) analysis, respectively. After that, we constructed a circRNA-miRNA-hub gene network and validated the circRNAs and mRNAs using an independent dataset (GSE61144) as well as qRT-PCR. Finally, we used CIBERSORTx database to analyze the immune infiltration characteristics of AMI and the correlation between hub genes and immune cells. Results: Using the "limma" package of the R, 83 DEcircRNAs, 54 DEmiRNAs, and 754 DEmRNAs were identified in the microarray datasets of AMI. Among 83 DEcircRNAs, there are 55 exonic DEcircRNAs. Then, a circRNA-miRNA-mRNA network consists of 21 DEcircRNAs, 11 DEmiRNAs, and 106 DEmRNAs were predicted by the database. After that, 10 hub genes from the PPI network were identified. Then, a new circRNA-miRNA-hub gene network consists of 14 DEcircRNAs, 7 DEmiRNAs, and 9 DEmRNAs was constructed. After that, three key circRNAs (hsa_circ_0009018, hsa_circ_0030569 and hsa_circ_0031017) and three hub genes (BCL6, PTGS2 and PTEN) were identified from the network by qRT-PCR. Finally, immune infiltration analysis showed that hub genes were significantly positively correlated with up-regulated immune cells (neutrophils, macrophages and plasma cells) in AMI. Conclusion: Our study constructed a circRNA-related ceRNA networks in AMI, consists of hsa_circ_0031017/hsa-miR-142-5p/PTEN axis, hsa_circ_0030569/hsa-miR-545/PTGS2 axis and hsa_circ_0009018/hsa-miR-139-3p/BCL6 axis. These three hub genes were significantly positively correlated with up-regulated immune cells (neutrophils, macrophages and plasma cells) in AMI. It helps improve understanding of AMI mechanism and provides future potential therapeutic targets.

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.