Abstract
BACKGROUND: Osteoarthritis (OA) is a worldwide chronic disease of the articulating joints. An increasing body of data demonstrates the immune system's involvement in osteoarthritis. The molecular mechanisms of OA are still unclear. This study aimed to search for OA immunerelated hub genes and determine appropriate diagnostic markers to help the detection and treatment of the disease. METHODS: Gene expression data were downloaded from the GEO database. Firstly, we analyzed and identified the differentially expressed genes (DEGs) using R packages. Meanwhile, ssGSEA was used to determine the activation degree of immune-related genes (IRGs), and WGCNA analysis was applied to search for co-expressed gene modules associated with immune cells. Then, critical networks and hub genes were found in the PPI network. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analyzed the biological functions of genes. The ability of the hub genes to differentiate OA from controls was assessed by the area under the ROC curve. A miRNA and transcription factor (TF) regulatory network was constructed according to their relationship with hub genes. Finally, the validation of hub genes was carried out by qPCR. RESULTS: In total, 353 DEGs were identified in OA patients compared with controls, including 222 upregulated and 131 downregulated genes. WGCNA successfully identified 34 main functional modules involved in the pathogenesis of OA. The most crucial functional module involved in OA included 89 genes. 19 immune-related genes were obtained by overlapping DEGs with the darkgrey module. The String database was constructed using the protein-protein interaction (PPI) network of 19 target genes, and 7 hub genes were identified by MCODE. ROC curve showed that 7 hub genes were potential biomarkers of OA. The expression levels of hub genes were validated by qPCR, and the results were consistent with those from bioinformatic analyses. CONCLUSION: Immune-related hub genes, including TYROBP, ITGAM, ITGB2, C1QC, MARCO, C1QB, and TLR8, may play critical roles in OA development. ITGAM had the highest correction on immune cells.