Identification of Candidate Genetic Markers and a Novel 4-genes Diagnostic Model in Osteoarthritis through Integrating Multiple Microarray Data

Identification of key hub genes in knee osteoarthritis through integrated bioinformatics analysis

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.

Identification and Verification of Endoplasmic Reticulum Stress-Related Genes as Novel Signatures for Osteoarthritis Diagnosis and Therapy: A Bioinformatics Analysis-Oriented Pilot Study

BACKGROUND AND PURPOSE

Endoplasmic reticulum stress (ERS) has been reported to be closely associated with the development of osteoarthritis (OA), but the underlying mechanisms are not fully delineated. The present study was designed to investigate the involvement of ERS-related genes in regulating OA progression.

METHODS

The expression profiles of OA patients and normal people were downloaded from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) in datasets GSE55457 and GSE55235 were screened and identified by R software with the construction of the protein-protein interaction (PPI) networks. Through the STRING and Venn diagram analysis, hub ERS-related genes were obtained. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed. Biomarkers with high diagnostic values of osteoarthritis (OA) were studied. The hematoxylin and eosin (H&E) staining and micro-CT were applied to evaluate the establishment of the OA model. The expression levels of biomarkers were validated with the use of reverse transcription‑quantitative polymerase chain reaction (RT-qPCR) and western blot. Finally, we evaluated the correlations of hub ERS-related genes with the immune infiltration cells via the CIBERSORT algorithm.

RESULTS

A total of 60 downregulated and 52 upregulated DEGs were identified, and the following GO and KEGG pathway analyses verified that those DEGs were mainly enriched in biological process (BP), cellular component (CC), molecular function (MF), and inflammation-associated signal pathways. Interestingly, among all the DEGs, six ER stress-associated genes, including activating transcription factor 3 (ATF3), DEAD-Box Helicase 3 X-Linked (DDX3X), AP-1 transcription factor subunit (JUN), eukaryotic initiation factor 4 (EIF4A1), KDEL endoplasmic reticulum protein retention receptor 3 (KDELR3), and vascular endothelial growth factor A (VEGFA), were found to be closely associated with OA progression, and the following RT-qPCR and Western Blot analysis confirmed that DDX3X, JUN, and VEGFA were upregulated, whereas KDELR3, EIF4A1, and ATF3 were downregulated in OA rats tissues compared to the normal tissues, which were in accordance with our bioinformatics findings. Furthermore, our receiver operating characteristic (ROC) curve analysis verified that the above six ER stress-associated genes could be used as ideal biomarkers for OA diagnosis and those genes also potentially regulated immune responses by influencing the biological functions of mast cells and macrophages.

CONCLUSION

Collectively, the present study firstly identified six ER stress-associated genes (ATF3, DDX3X, JUN, EIF4A1, KDELR3, and VEGFA) that may play critical role in regulating the progression of OA.

ATF3 as a potential diagnostic marker of early-stage osteoarthritis and its correlation with immune infiltration through bioinformatics analysis

AIMS

This study aimed, through bioinformatics analysis, to identify the potential diagnostic markers of osteoarthritis, and analyze the role of immune infiltration in synovial tissue.

METHODS

The gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified by R software. Functional enrichment analyses were performed and protein-protein interaction networks (PPI) were constructed. Then the hub genes were screened. Biomarkers with high value for the diagnosis of early osteoarthritis (OA) were validated by GEO datasets. Finally, the CIBERSORT algorithm was used to evaluate the immune infiltration between early-stage OA and end-stage OA, and the correlation between the diagnostic marker and infiltrating immune cells was analyzed.

RESULTS

A total of 88 DEGs were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DEGs were significantly enriched in leucocyte migration and interleukin (IL)-17 signalling pathways. Disease ontology (DO) indicated that DEGs were mostly enriched in rheumatoid arthritis. Six hub genes including FosB proto-oncogene, AP-1 transcription factor subunit (FOSB); C-X-C motif chemokine ligand 2 (CXCL2); CXCL8; IL-6; Jun proto-oncogene, AP-1 transcription factor subunit (JUN); and Activating transcription factor 3 (ATF3) were identified and verified by GEO datasets. ATF3 (area under the curve = 0.975) turned out to be a potential biomarker for the diagnosis of early OA. Several infiltrating immune cells varied significantly between early-stage OA and end-stage OA, such as resting NK cells (p = 0.016), resting dendritic cells (p = 0.043), and plasma cells (p = 0.043). Additionally, ATF3 was significantly correlated with resting NK cells (p = 0.034), resting dendritic cells (p = 0.026), and regulatory T cells (Tregs, p = 0.018).

CONCLUSION

ATF3 may be a potential diagnostic marker for early diagnosis and treatment of OA, and immune cell infiltration provides new perspectives for understanding the mechanism during OA progression.Cite this article: Bone Joint Res 2022;11(9):679-689.