Search

A growing of evidence has showed that patients with osteoarthritis (OA) had a higher coronavirus 2019 (COVID-19) infection rate and a poorer prognosis after infected it. Additionally, scientists have also discovered that COVID-19 infection might cause pathological changes in the musculoskeletal system. However, its mechanism is still not fully elucidated. This study aims to further explore the sharing pathogenesis of patients with both OA and COVID-19 infection and find candidate drugs. Gene expression profiles of OA (GSE51588) and COVID-19 (GSE147507) were obtained from the Gene Expression Omnibus (GEO) database. The common differentially expressed genes (DEGs) for both OA and COVID-19 were identified and several hub genes were extracted from them. Then gene and pathway enrichment analysis of the DEGs were performed; protein–protein interaction (PPI) network, transcription factor (TF)-gene regulatory network, TF-miRNA regulatory network and gene-disease association network were constructed based on the DEGs and hub genes. Finally, we predicted several candidate molecular drugs related to hub genes using DSigDB database. The receiver operating characteristic curve (ROC) was applied to evaluate the accuracy of hub genes in the diagnosis of both OA and COVID-19. In total, 83 overlapping DEGs were identified and selected for subsequent analyses. CXCR4, EGR2, ENO1, FASN, GATA6, HIST1H3H, HIST1H4H, HIST1H4I, HIST1H4K, MTHFD2, PDK1, TUBA4A, TUBB1 and TUBB3 were screened out as hub genes, and some showed preferable values as diagnostic markers for both OA and COVID-19. Several candidate molecular drugs, which are related to the hug genes, were identified. These sharing pathways and hub genes may provide new ideas for further mechanistic studies and guide more individual-based effective treatments for OA patients with COVID-19 infection. [ABSTRACT FROM AUTHOR]

The overall study design is illustrated. [Display omitted] • A four-gene signature (ZBTB16, TNFSF11, SCRG1 and KDELR3) was identified to be involved OA progression based on WGCNA, LASSO, SVM-RFE and RF algorithms combined with bioinformatic analyses. • The immune infiltration analyses revealed that hub genes were most correlated with regulatory T cell, natural killer cell and CD8 + T cell. Although osteoarthritis (OA) is one of the most prevalent joint disorders, effective biomarkers to diagnose OA are still unavailable. This study aimed to acquire some key synovial biomarkers (hub genes) and analyze their correlation with immune infiltration in OA. Gene expression profiles and clinical characteristics of OA and healthy synovial samples were retrieved from the Gene Expression Omnibus (GEO) database. Hub genes for OA were mined based on a combination of weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) algorithms. A diagnostic nomogram model for OA prediction was developed based on the hub genes. Receiver operating characteristic curves (ROC) were performed to confirm the abnormal expression of hub genes in the experimemtal and validation datasets. qRT-PCR using patients' samples were conducted as well. In addition, the infiltration level of 28 immune cells in the expression profile and their relationship with hub genes were analyzed using single-sample GSEA (ssGSEA). 4 hub genes (ZBTB16, TNFSF11, SCRG1 and KDELR3) were obtained by WGCNA, lasso, SVM-RFE, RF algorithms as potential biomarkers for OA. The immune infiltration analyses revealed that hub genes were most correlated with regulatory T cell and natural killer cell. A machine learning model to diagnose OA based on ZBTB16, TNFSF11, SCRG1 and KDELR3 using synovial tissue was constructed, providing theoretical foundation and guideline for diagnostic and treatment targets in OA. [ABSTRACT FROM AUTHOR]