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.