Identification of susceptibility modules and characteristic genes to osteoarthritis by WGCNA

The susceptibility modules and characteristic genes of patients with osteoarthritis (OA) were determined by weighted gene co-expression network analysis (WGCNA), and the role of immune cells in OA related microenvironment was analyzed. GSE98918 and GSE117999 data sets are from GEO database. R language was used to conduct difference analysis for the new data set after merging. The formation of gene co-expression network, screening of susceptibility modules and screening of core genes are all through WGCNA. GO and KEGG enrichment analyses were used for Hub genes. The characteristic genes of the disease were obtained by Lasso regression screening. SSGSEA was used to estimate immune cell abundance in sample and a series of correlation analyses were performed. WGCNA was used to form 6 gene co-expression modules. The yellow-green module is identified as the susceptible module of OA. 202 genes were identified as core genes. Finally, RHOT2, FNBP4 and NARF were identified as the characteristic genes of OA. The results showed that the characteristic genes of OA were positively correlated with plasmacytoid dendritic cells, NKT cells and immature dendritic cells, but negatively correlated with active B cells. MDSC were the most abundant immune cells in cartilage. This study identified the Hippo signaling pathway, mTOR signaling pathway, and three characteristic genes (RHOT2, FNBP4, NARF) as being associated with osteoarthritis (OA). These three genes are downregulated in the cartilage of OA patients and may serve as biomarkers for early diagnosis and targeted therapy. Proper regulation of immune cells may aid in the treatment of OA. Future research should focus on developing tools to detect these genes and exploring their therapeutic applications.