Dataset Information

Identification of key genes and small molecule drugs in osteoarthritis by integrated bioinformatics analysis.

ABSTRACT:

Background

Osteoarthritis (OA) is a common joint degenerative disease that can affect multiple joints. Genetic events may play an important regulatory role in the early stages of the disease, but the specific mechanisms have not yet been fully elucidated. The main purpose of this study was to screen for disease-causing hub genes and effective small molecule drugs to reveal the pathogenesis of OA and to develop novel drugs for treatment.

Methods

Two gene expression profile datasets, GSE55235 and GSE55457, were integrated and further analyzed. The consistently differentially expressed genes (DEGs) were identified, and functional annotation and pathway analysis of these genes were performed with GO and KEGG. A protein-protein interaction network (PPI) of the DEGs was generated using STRING, and potential small molecule drug screening was performed on the connectivity map (CMap).

Results

A total of 158 consistently differentially expressed genes were identified from the two profile datasets. The functions of these DEGs are mainly related to the TNF signaling pathway, osteoclast differentiation, MAPK signaling pathway and so on. The PPI network contains 127 nodes and 1802 edges, and the ten hub genes were interleukin 6 (IL6), vascular endothelial growth factor A (VEGFA)and so on. 7 small molecule drugs were identified as potential interactors with these hubs.

Conclusions

This study explains the disorder of expression in the pathological process of OA at transcriptome, which will help to understand the pathogenesis of OA.

SUBMITTER: Liu Z

PROVIDER: S-EPMC10009689 | biostudies-literature | 2023 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification of key genes and small molecule drugs in osteoarthritis by integrated bioinformatics analysis.

Liu Zhendong Z Wang Hongbo H Cheng Xingbo X Zhang Jiangfen J Gao Yanzheng Y

Biochemistry and biophysics reports 20230305

<h4>Background</h4>Osteoarthritis (OA) is a common joint degenerative disease that can affect multiple joints. Genetic events may play an important regulatory role in the early stages of the disease, but the specific mechanisms have not yet been fully elucidated. The main purpose of this study was to screen for disease-causing hub genes and effective small molecule drugs to reveal the pathogenesis of OA and to develop novel drugs for treatment.<h4>Methods</h4>Two gene expression profile datasets ...[more]

PMID: 36923006

Similar Datasets

Project description:Knee osteoarthritis (KOA) is a common chronic joint disease globally. Synovial inflammation plays a pivotal role in its pathogenesis, preceding cartilage damage. Identifying biomarkers in osteoarthritic synovial tissues holds promise for early diagnosis and targeted interventions. Gene expression profiles were obtained from the Gene Expression Omnibus database. Subsequent analyses included differential expression gene (DEG) analysis and weighted gene co-expression network analysis (WGCNA) on the combined datasets. We performed functional enrichment analysis on the overlapping genes between DEGs and module genes and constructed a protein-protein interaction network. Using Cytoscape software, we identified hub genes related to the disease and conducted gene set enrichment analysis on these hub genes. The CIBERSORT algorithm was employed to evaluate the correlation between hub genes and the abundance of immune cells within tissues. Finally, Mendelian randomization analysis was utilized to assess the potential of these hub genes as biomarkers. We identified 46 differentially expressed genes (DEGs), comprising 20 upregulated and 26 downregulated genes. Using WGCNA, we constructed a gene co-expression network and selected the most relevant modules, resulting in 24 intersecting genes with the DEGs. KEGG enrichment analysis of the intersecting genes identified the IL-17 signaling pathway, associated with inflammation, as the most significant pathway. Cytoscape software was utilized to rank the candidate genes, with JUN, ATF3, FOSB, NR4A2, and IL6 emerging as the top five based on the Degree algorithm. A nomogram model incorporating these five genes, supported by ROC curve analysis, validated their diagnostic efficacy. Immune infiltration and correlation analysis revealed that macrophages were significantly associated with JUN (p < 0.01), FOSB (p < 0.01), and NR4A2 (p < 0.05). Additionally, T follicular helper cells showed significant associations with ATF3 (p < 0.05), FOSB (p < 0.05), and JUN (p < 0.05). Mendelian randomization analysis provided strong evidence linking JUN (IVW: OR = 0.910, p = 0.005) and IL6 (IVW: OR = 1.024, p = 0.026) with KOA. Through the utilization of various bioinformatics analysis methods, we have pinpointed key hub genes relevant to knee osteoarthritis. These findings hold promise for advancing pre-symptomatic diagnostic strategies and enhancing our understanding of the biological underpinnings behind knee osteoarthritis susceptibility genes.

Project description:Objective: The currently established diagnostic and prognostic tools for diabetic kidney disease (DKD) have limitations, which demands the necessity to find new genes and pathways associated with diagnosis and treatment. Our study aims to reveal the gene expression alteration and discover critical genes involved in the development of DKD, thus providing novel diagnostic molecular markers and therapeutic targets. Materials and methods: The differences of infiltrating immune cells within kidney were compared between healthy living donors and DKD patients. Besides, differentially expressed genes (DEGs) within kidney from healthy living donor, early stage DKD and advanced stage DKD samples were detected. Furthermore, the weighted co-expressed network (WGCNA) and protein-protein interaction (PPI) network were constructed, followed by recognition of core hub genes and module analysis. Receiver operating characteristic (ROC) curve analysis was implemented to determine the diagnostic value of hub genes, correlation analysis was employed to explore the association between hub genes and infiltrating immune cells, and certain hub genes was validated by quantitative real-time PCR and immunohistochemistry staining in cultured tubule cells and diabetic mice kidney. Finally, the candidate small molecules as potential drugs to treat DKD were anticipated through utilizing virtual screening and molecular docking investigation. Results: Our study revealed significantly higher proportion of infiltrating immune cells within kidney from DKD patients via probing the immune landscape by single-cell transcriptomics. Besides, 126 commonly shared DEGs identified among three group samples were enriched in immune biological process. In addition, the ROC curve analysis demonstrated the strong diagnostic accuracy of recognized hub genes (NFKB1, DYRK2, ATAD2, YAP1, and CHD3) from PPI network. Correlation analysis further confirmed the positive association between these hub genes with infiltrating natural killer cells. More importantly, the mRNA transcripts and protein abundance of YAP1 were significantly higher in high glucose-treated renal tubule cells and diabetic mice kidney, and the small molecules exhibiting the best binding affinities with YAP1 were predicted and acquired. Conclusion: Our findings for the first time indicate that NFKB1, DYRK2, ATAD2, YAP1, and CHD3 might be potential novel biomarkers and therapeutic targets for DKD, providing insights into the molecular mechanisms underlying the pathogenesis of DKD.

Project description:The present study aimed to identify key genes involved in osteoarthritis (OA). Based on a bioinformatics analysis of five gene expression profiling datasets (GSE55457, GSE55235, GSE82107, GSE12021 and GSE1919), differentially expressed genes (DEGs) in OA were identified. Subsequently, a protein-protein interaction (PPI) network was constructed and its topological structure was analyzed. In addition, key genes in OA were identified following a principal component analysis (PCA) based on the DEGs in the PPI network. Finally, the functions and pathways enriched by these key genes were also analyzed. The PPI network consisted of 241 nodes and 576 interactives, including a total of 171 upregulated DEGs [e.g., aspartylglucosaminidase (AGA), CD58 and CD86] and a total of 70 downregulated DEGs (e.g., acetyl-CoA carboxylase β and dihydropyrimidine dehydrogenase). The PPI network complied with an attribute of scale-free small-world network. After PCA, 47 key genes were identified, including β-1,4-galactosyltransferase-1 (B4GALT1), AGA, CD58, CD86, ezrin, and eukaryotic translation initiation factor 4 γ 1 (EIF4G1). Subsequently, the 47 key genes were identified to be enriched in 13 Gene Ontology (GO) terms and 2 Kyoto Encyclopedia of Genes and Genomes pathways, with the GO terms involving B4GALT1 including positive regulation of developmental processes, protein amino acid terminal glycosylation and protein amino acid terminal N-glycosylation. In addition, B4GALT1 and EIF4G1 were confirmed to be downregulated in OA samples compared with healthy controls, but only EIF4G1 was determined to be significantly downregulated in OA samples, as determined via a meta-analysis of the 5 abovementioned datasets. In conclusion, B4GALT1 and EIF4G1 were indicated to have significant roles in OA, and B4GALT1 may be involved in positive regulation of developmental processes, protein amino acid terminal glycosylation and protein amino acid terminal N-glycosylation. The present study may enhance the current understanding of the molecular mechanisms of OA and provide novel therapeutic targets.

Project description:BackgroundAtherosclerosis (AS) is a chronic inflammatory disease that might induce severe cardiovascular events, such as myocardial infarction and cerebral infarction. These risk factors in the pathogenesis of AS remain uncertain and further research is needed. This study aims to explore the potential molecular mechanisms of AS by bioinformatics analyses.MethodsGSE100927 gene expression profiles, including 69 AS samples and 35 healthy controls, were downloaded from Gene Expression Omnibus database and indenfied for key genes and pathways in AS.ResultsA total of 443 differentially expressed genes (DEGs) between control and AS were identified, including 323 down-regulated genes and 120 up-regulated genes. The Gene ontology terms enriched by the up-regulated DEGs were associated with the regulation of leukocyte activation, endocytic vesicle, and cytokine binding, while the down-regulated DEGs were associated with negative regulation of cell growth, extracellular matrix, and G protein-coupled receptor binding. KEGG pathway analysis showed that the up-regulated DEGs were enriched in Osteoclast differentiation and Phagosome, while the down-regulated DEGs were enriched in vascular smooth muscle contraction and cGMP-PKG signaling pathway. Using the modular analysis of Cytoscape, we identified 3 modules mainly involved in Leishmaniasis and Osteoclast differentiation. The GSEA analysis showed the up-regulated gene sets were enriched in the ribosome, ascorbated metabolism, and propanoate metabolism. The LASSO Cox regression analysis showed the top 3 genes were TNF, CX3CR1, and COL1R1. Finally, we found these immune cells were conferred significantly higher infiltrating density in the AS group.ConclusionsOur data showed the pathway of Osteoclast differentiation and Leishmaniasis was involved in the AS process and we developed a three-gene model base on the prognosis of AS. These findings clarified the gene regulatory network of AS and may provide a novel target for AS therapy.

Project description:Colorectal cancer (CRC) is a life-threatening disease with a poor prognosis. Therefore, it is crucial to identify molecular prognostic biomarkers for CRC. The present study aimed to identify potential key genes that could be used to predict the prognosis of patients with CRC. Three CRC microarray datasets (GSE20916, GSE73360 and GSE44861) were downloaded from the Gene Expression Omnibus (GEO) database, and one dataset was obtained from The Cancer Genome Atlas (TCGA) database. The three GEO datasets were analyzed to detect differentially expressed genes (DEGs) using the BRB-ArrayTools software. Functional and pathway enrichment analyses of these DEGs were performed using the Database for Annotation, Visualization and Integrated Discovery tool. A protein-protein interaction (PPI) network of DEGs was constructed, hub genes were extracted, and modules of the PPI network were analyzed. To investigate the prognostic values of the hub genes in CRC, data from the CRC datasets of TCGA were used to perform the survival analyses based on the sample splitting method and Cox regression model. Correlation among the hub genes was evaluated using Spearman's correlation analysis. In the three GEO datasets, a total of 105 common DEGs were identified, including 51 down- and 54 up-regulated genes in CRC compared with normal colorectal tissues. A PPI network consisting of 100 DEGs and 551 edges was constructed, and 44 nodes were identified as hub genes. Among these 44 genes, the four hub genes TIMP metallopeptidase inhibitor 1 (TIMP1), solute carrier family 4 member 4 (SLC4A4), aldo-keto reductase family 1 member B10 (AKR1B10) and ATP binding cassette subfamily E member 1 (ABCE1) were associated with overall survival (OS) in patients with CRC. Three significant modules were extracted from the PPI network. The hub gene TIMP1 was present in Module 1, ABCE1 was involved in Module 2 and SLC4A4 was identified in Module 3. Univariate analysis revealed that TIMP1, SLC4A4, AKR1B10 and ABCE1 were associated with the OS of patients with CRC. Multivariate analysis demonstrated that SLC4A4 may be an independent prognostic factor associated with OS. Furthermore, the results from correlation analysis revealed that there was no correlation between TIMP1, SLC4A4 and ABCE1, whereas AKR1B10 was positively correlated with SLC4A4. In conclusion, the four key genes TIMP1, SLC4A4, AKR1B10 and ABCE1 associated with the OS of patients with CRC were identified by integrated bioinformatics analysis. These key genes may be used as prognostic biomarkers to predict the survival of patients with CRC, and may therefore represent novel therapeutic targets for CRC.

Dataset Information

Identification of key genes and small molecule drugs in osteoarthritis by integrated bioinformatics analysis.

Background

Methods

Results

Conclusions

Publications

Identification of key genes and small molecule drugs in osteoarthritis by integrated bioinformatics analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets