Dataset Information

Identification of Novel Kinase–Transcription Factor–mRNA–miRNA Regulatory Network in Nasopharyngeal Carcinoma by Bioinformatics Analysis

ABSTRACT:

Purpose

Nasopharyngeal carcinoma (NPC) is one of the most common malignant tumors of the head and neck. This study aimed to investigate the crucial genes and regulatory networks involved in the carcinogenesis of NPC using a bioinformatics approach.

Methods

Five mRNA and two miRNA expression datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and miRNAs (DEMs) between NPC and normal samples were analyzed using R software. The WebGestalt tool was used for functional enrichment analysis, and protein–protein interaction (PPI) network analysis of DEGs was performed using STRING database. Transcription factors (TFs) were predicted using TRRUST and Transcriptional Regulatory Element Database (TRED). Kinases were identified using X2Kgui. The miRNAs of DEGs were predicted using miRWalk database. A kinase–TF–mRNA–miRNA integrated network was constructed, and hub nodes were selected. The hub genes were validated using NPC datasets from the GEO and Oncomine databases. Finally, candidate small-molecule agents were predicted using CMap.

Results

A total of 122 DEGs and 44 DEMs were identified. DEGs were associated with the immune response, leukocyte activation, endoplasmic reticulum stress in GO analysis, and the NF-κB signaling pathway in KEGG analysis. Four significant modules were identified using PPI network analysis. Subsequently, 26 TFs, 73 kinases, and 2499 miRNAs were predicted. The predicted miRNAs were cross-referenced with DEMs, and seven overlapping miRNAs were selected. In the kinase–TF–mRNA–miRNA integrated network, eight genes (PTGS2, FN1, MMP1, PLAU, MMP3, CD19, BMP2, and PIGR) were identified as hub genes. Hub genes were validated with consistent results, indicating the reliability of our findings. Finally, six candidate small-molecule agents (phenoxybenzamine, luteolin, thioguanosine, reserpine, blebbistatin, and camptothecin) were predicted.

Conclusion

We identified DEGs and an NPC regulatory network involving kinases, TFs, mRNAs, and miRNAs, which might provide promising insight into the pathogenesis, treatment, and prognosis of NPC.

SUBMITTER: Gao L

PROVIDER: S-EPMC8566004 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundMedullary thyroid carcinoma (MTC), a thyroid C cell-derived malignancy, is poorly differentiated and more aggressive than papillary, follicular and oncocytic types of thyroid cancer. The current therapeutic options are limited, with a third of population suffering resistance. The differential gene expression pattern among thyroid cancer subtypes remains unclear. This study intended to explore the exclusive gene profile of MTC and construct a comprehensive regulatory network via integrated analysis, to uncover the potential key biomarkers.MethodsMultiple datasets of thyroid and other neuroendocrine tumors were obtained from GEO and TCGA databases. Differentially expressed genes (DEGs) specific in MTC were identified to construct a transcription factor (TF)-mRNA-miRNA network. The impact of the TF-mRNA-miRNA network on tumor immune characteristics and patient survival was further explored by single-sample GSEA (ssGSEA) and ESTIMATE algorithms, as well as univariate combined with multivariate analyses. RT-qPCR, cell viability and apoptosis assays were performed for in vitro validation.ResultsWe identified 81 genes upregulated and 22 downregulated in MTC but not in other types of thyroid tumor compared to the normal thyroid tissue. According to the L1000CDS2 database, potential targeting drugs were found to reverse the expressions of DEGs, with panobinostat (S1030) validated effective for tumor repression in MTC by in vitro experiments. The 103 DEGs exclusively seen in MTC were involved in signal release, muscle contraction, pathways of neurodegeneration diseases, neurotransmitter activity and related amino acid metabolism, and cAMP pathway. Based on the identified 15 hub genes, a TF-mRNA-miRNA linear network, as well as REST-cored coherent feed-forward loop networks, namely REST-KIF5C-miR-223 and REST-CDK5R2-miR-130a were constructed via online prediction and validation by public datasets and our cohort. Hub-gene, TF and miRNA scores in the TF-mRNA-miRNA network were related to immune score, immune cell infiltration and immunotherapeutic molecules in MTC as well as in neuroendocrine tumor of lung and neuroblastoma. Additionally, a high hub-gene score or a low miRNA score indicated good prognoses of neuroendocrine tumors.ConclusionThe present study uncovers underlying molecular mechanisms and potential immunotherapy-related targets for the pathogenesis and drug discovery of MTC.

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: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: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: Osteoarthritis (OA) is a degenerative joint disease that seriously affects the quality of people. Unfortunately, the pathogenesis of OA has not been fully known. Therefore, this study aimed to construct a ceRNA regulatory network related to OA to explore the pathogenesis of OA. Methods: Differentially expressed circRNAs (DEcircRNAs), microRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were obtained from the Gene Expression Omnibus microarray data (GSE175959, GSE105027, and GSE169077). The miRNA response elements and target mRNAs were identified using bioinformatics approaches. Additionally, a circRNA-miRNA-mRNA network was established using Cytoscape version 3.8.0. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of mRNAs in the network were conducted to explore the possible mechanisms underlying OA development. Protein-protein interaction (PPI) analysis was performed to determine the hub genes. Based on the hub genes, a sub network was constructed using Cytoscape 3.8.0 version. Finally, connectivity map (CMap) and drug-gene interaction database (DGIdb) analyses were performed to identify the potential therapeutic targets for OA. Results: Altogether, five DEcircRNAs, 89 DEmiRNAs, and 345 DEmRNAs were identified. Moreover, a circRNA-miRNA-mRNA network was established using three circRNAs, seven miRNAs, and 37 mRNAs. GO and KEGG analyses demonstrated that the mRNAs in the network could be related to the occurrence and development of OA. PPI analysis was performed and six key genes, namely serpin family H member 1 [SERPINH1], collagen type VIII alpha 2 chain [COL8A2], collagen type XV alpha 1 chain [COL15A1], collagen type VI alpha 3 chain [COL6A3], collagen type V alpha 1 chain [COL5A1], and collagen type XI alpha 1 chain [COL11A1], were identified. Furthermore, a circRNA-miRNA-hub gene subnetwork was established in accordance with two circRNAs (hsa_circ_0075320 and hsa_circ_0051428), two miRNAs (hsa-miR-6124 and hsa-miR-1207-5p), and six hub genes (COL11A1, SERPINH1, COL6A3, COL5A1, COL8A2, and COL15A1). Finally, three chemicals (noscapine, diazepam, and TG100-115) based on CMap analysis and two drugs (collagenase Clostridium histolyticum and ocriplasmin) based on DGIdb were discovered as potential treatment options for OA. Conclusion: This study presents novel perspectives on the pathogenesis and treatment of OA based on circRNA-related competitive endogenous RNA regulatory networks.

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.