Identification of osteoarthritis-characteristic genes and immunological micro-environment features through bioinformatics and machine learning-based approaches

BMC Med Genomics. 2023 Oct 7;16(1):236. doi: 10.1186/s12920-023-01672-y.

Abstract

Background: Osteoarthritis (OA) is a multifaceted chronic joint disease characterized by complex mechanisms. It has a detrimental impact on the quality of life for individuals in the middle-aged and elderly population while also imposing a significant socioeconomic burden. At present, there remains a lack of comprehensive understanding regarding the pathophysiology of OA. The objective of this study was to examine the genes, functional pathways, and immune infiltration characteristics associated with the development and advancement of OA.

Methods: The Gene Expression Omnibus (GEO) database was utilized to acquire gene expression profiles. The R software was employed to conduct the screening of differentially expressed genes (DEGs) and perform enrichment analysis on these genes. The OA-characteristic genes were identified using the Weighted Gene Co-expression Network Analysis (WGCNA) and the Lasso algorithm. In addition, the infiltration levels of immune cells in cartilage were assessed using single-sample gene set enrichment analysis (ssGSEA). Subsequently, a correlation analysis was conducted to examine the relationship between immune cells and the OA-characteristic genes.

Results: A total of 80 DEGs were identified. As determined by functional enrichment, these DEGs were associated with chondrocyte metabolism, apoptosis, and inflammation. Three OA-characteristic genes were identified using WGCNA and the lasso algorithm, and their expression levels were then validated using the verification set. Finally, the analysis of immune cell infiltration revealed that T cells and B cells were primarily associated with OA. In addition, Tspan2, HtrA1 demonstrated a correlation with some of the infiltrating immune cells.

Conclusions: The findings of an extensive bioinformatics analysis revealed that OA is correlated with a variety of distinct genes, functional pathways, and processes involving immune cell infiltration. The present study has successfully identified characteristic genes and functional pathways that hold potential as biomarkers for guiding drug treatment and facilitating molecular-level research on OA.

Keywords: Bioinformatics; Biomarkers; Cartilage; Herc5; HtrA1; Osteoarthritis; Tspan2.

MeSH terms

  • Aged
  • Computational Biology
  • High-Temperature Requirement A Serine Peptidase 1
  • Humans
  • Inflammation
  • Machine Learning
  • Middle Aged
  • Osteoarthritis* / genetics
  • Quality of Life*

Substances

  • High-Temperature Requirement A Serine Peptidase 1