Abstract
BACKGROUND: Our study aimed to identify potential specific biomarkers for osteoarthritis (OA) and assess their relationship with immune infiltration. METHODS: We utilized data from GSE117999, GSE51588, and GSE57218 as training sets, while GSE114007 served as a validation set, all obtained from the GEO database. First, weighted gene co-expression network analysis (WGCNA) and functional enrichment analysis were performed to identify hub modules and potential functions of genes. We subsequently screened for potential OA biomarkers within the differentially expressed genes (DEGs) of the hub module using machine learning methods. The diagnostic accuracy of the candidate genes was validated. Additionally, single gene analysis and ssGSEA was performed. Then, we explored the relationship between biomarkers and immune cells. Lastly, we employed RT-PCR to validate our results. RESULTS: WGCNA results suggested that the blue module was the most associated with OA and was functionally associated with extracellular matrix (ECM)-related terms. Our analysis identified ALB, HTRA1, DPT, MXRA5, CILP, MPO, and PLAT as potential biomarkers. Notably, HTRA1, DPT, and MXRA5 consistently exhibited increased expression in OA across both training and validation cohorts, demonstrating robust diagnostic potential. The ssGSEA results revealed that abnormal infiltration of DCs, NK cells, Tfh, Th2, and Treg cells might contribute to OA progression. HTRA1, DPT, and MXRA5 showed significant correlation with immune cell infiltration. The RT-PCR results also confirmed these findings. CONCLUSIONS: HTRA1, DPT, and MXRA5 are promising biomarkers for OA. Their overexpression strongly correlates with OA progression and immune cell infiltration.