Identification of Immune-Related Risk Genes in Osteoarthritis Based on Bioinformatics Analysis and Machine Learning

In this research, we aimed to perform a comprehensive bioinformatic analysis of immune cell infiltration in osteoarthritic cartilage and synovium and identify potential risk genes. Datasets were downloaded from the Gene Expression Omnibus database. We integrated the datasets, removed the batch effects and analyzed immune cell infiltration along with differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA) was used to identify the positively correlated gene modules. LASSO (least absolute shrinkage and selection operator)-cox regression analysis was performed to screen the characteristic genes. The intersection of the DEGs, characteristic genes and module genes was identified as the risk genes. The WGCNA analysis demonstrates that the blue module was highly correlated and statistically significant as well as enriched in immune-related signaling pathways and biological functions in the KEGG and GO enrichment. LASSO-cox regression analysis screened 11 characteristic genes from the hub genes of the blue module. After the DEG, characteristic gene and immune-related gene datasets were intersected, three genes, PTGS1, HLA-DMB and GPR137B, were identified as the risk genes in this research. In this research, we identified three risk genes related to the immune system in osteoarthritis and provide a feasible approach to drug development in the future.

[1]  Haiyuan Lv,et al.  Identification of the osteoarthritis signature gene PDK1 by machine learning and its regulatory mechanisms on chondrocyte autophagy and apoptosis , 2023, Frontiers in Immunology.

[2]  Anyu Zeng,et al.  Identification of SLC3A2 as a Potential Therapeutic Target of Osteoarthritis Involved in Ferroptosis by Integrating Bioinformatics, Clinical Factors and Experiments , 2022, Cells.

[3]  Yihao Liang,et al.  Identification of Biomarkers Associated with Diagnosis of Osteoarthritis Patients Based on Bioinformatics and Machine Learning , 2022, Journal of immunology research.

[4]  Yiqin Zhou,et al.  Identification and development of a novel 5-gene diagnostic model based on immune infiltration analysis of osteoarthritis , 2021, Journal of translational medicine.

[5]  Letian Shan,et al.  Identification and validation of hub genes of synovial tissue for patients with osteoarthritis and rheumatoid arthritis , 2021, Hereditas.

[6]  H. Lorenz,et al.  Proinflammatory T cell polarization is already present in patients with early knee osteoarthritis , 2021, Arthritis Research & Therapy.

[7]  Wei Wang,et al.  Integration of gene interaction information into a reweighted Lasso-Cox model for accurate survival prediction , 2020, Bioinform..

[8]  Dongqiang Zeng,et al.  IOBR: Multi-Omics Immuno-Oncology Biological Research to Decode Tumor Microenvironment and Signatures , 2020, bioRxiv.

[9]  Chao Guo,et al.  The landscape of gene co-expression modules correlating with prognostic genetic abnormalities in AML , 2020, Journal of translational medicine.

[10]  S. Gregory,et al.  Synovial cell cross-talk with cartilage plays a major role in the pathogenesis of osteoarthritis , 2020, Scientific Reports.

[11]  A. Barbero,et al.  Anti-Inflammatory and Chondroprotective Effects of Vanillic Acid and Epimedin C in Human Osteoarthritic Chondrocytes , 2020, Biomolecules.

[12]  D. Bueno,et al.  Macrophage: A Potential Target on Cartilage Regeneration , 2020, Frontiers in Immunology.

[13]  X. Bai,et al.  Macrophages regulate the progression of osteoarthritis. , 2020, Osteoarthritis and cartilage.

[14]  Richard J. Miller,et al.  The Innate Immune Response As A Mediator Of Osteoarthritis Pain. , 2019, Osteoarthritis and cartilage.

[15]  J. Hadfield,et al.  RNA sequencing: the teenage years , 2019, Nature Reviews Genetics.

[16]  W. Robinson,et al.  IgE-mediated mast cell activation promotes inflammation and cartilage destruction in osteoarthritis , 2019, eLife.

[17]  Cuijiang Wang,et al.  Naproxen attenuates osteoarthritis progression through inhibiting the expression of prostaglandinl‐endoperoxide synthase 1 , 2018, Journal of cellular physiology.

[18]  K. Harrar,et al.  Combining GLCM with LBP features for knee osteoarthritis prediction: Data from the Osteoarthritis initiative , 2018, EAI Endorsed Trans. Scalable Inf. Syst..

[19]  Philip G Conaghan,et al.  Synovitis in osteoarthritis: current understanding with therapeutic implications , 2017, Arthritis Research & Therapy.

[20]  D Loeuille,et al.  Subchondral tibial bone texture analysis predicts knee osteoarthritis progression: data from the Osteoarthritis Initiative: Tibial bone texture & knee OA progression. , 2017, Osteoarthritis and cartilage.

[21]  G. Bernardini,et al.  CXCR3/CXCL10 Axis Regulates Neutrophil–NK Cell Cross-Talk Determining the Severity of Experimental Osteoarthritis , 2017, The Journal of Immunology.

[22]  Pornpimol Charoentong,et al.  Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade , 2016, bioRxiv.

[23]  Qian Wang,et al.  Low-grade inflammation as a key mediator of the pathogenesis of osteoarthritis , 2016, Nature Reviews Rheumatology.

[24]  Richard F. Loeser,et al.  Ageing and the pathogenesis of osteoarthritis , 2016, Nature Reviews Rheumatology.

[25]  C. Manferdini,et al.  From osteoarthritic synovium to synovial-derived cells characterization: synovial macrophages are key effector cells , 2016, Arthritis Research & Therapy.

[26]  Daniel J. Gaffney,et al.  A survey of best practices for RNA-seq data analysis , 2016, Genome Biology.

[27]  Ash A. Alizadeh,et al.  Robust enumeration of cell subsets from tissue expression profiles , 2015, Nature Methods.

[28]  V. Kraus,et al.  The Role of Innate Immunity in Osteoarthritis: When Our First Line of Defense Goes On the Offensive , 2015, The Journal of Rheumatology.

[29]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[30]  P. Richette,et al.  Biologic agents in osteoarthritis: hopes and disappointments , 2013, Nature Reviews Rheumatology.

[31]  P. Low,et al.  Direct in vivo evidence of activated macrophages in human osteoarthritis. , 2013, Osteoarthritis and cartilage.

[32]  Hugues Bersini,et al.  Unlocking the potential of publicly available microarray data using inSilicoDb and inSilicoMerging R/Bioconductor packages , 2012, BMC Bioinformatics.

[33]  J. Schoones,et al.  Synovial inflammation, immune cells and their cytokines in osteoarthritis: a review. , 2012, Osteoarthritis and cartilage.

[34]  Annette Lee,et al.  Genome-Wide Association Study of Rheumatoid Arthritis in Koreans: Population-Specific Loci as Well as Overlap With European Susceptibility Loci , 2011 .

[35]  Ali Guermazi,et al.  The role of the meniscus in knee osteoarthritis: a cause or consequence? , 2009, Radiologic clinics of North America.

[36]  S. Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[37]  S. Abramson,et al.  Current concepts in the pathogenesis of osteoarthritis. , 2008, Osteoarthritis and cartilage.

[38]  R. Lister,et al.  Highly Integrated Single-Base Resolution Maps of the Epigenome in Arabidopsis , 2008, Cell.

[39]  G. Castellani,et al.  Inflammaging and anti-inflammaging: A systemic perspective on aging and longevity emerged from studies in humans , 2007, Mechanisms of Ageing and Development.

[40]  B. Combe,et al.  HLA-DMA*0103 and HLA-DMB*0104 alleles as novel prognostic factors in rheumatoid arthritis , 2004, Annals of the rheumatic diseases.

[41]  L. Sundblad,et al.  Protein patterns in synovial fluid and serum in rheumatoid arthritis and osteoarthritis. , 1959, Arthritis and rheumatism.

[42]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

[43]  C. Brinckerhoff,et al.  Matrix metalloproteinases: role in arthritis. , 2006, Frontiers in bioscience : a journal and virtual library.

[44]  R. Villar.  Knee replacement. , 1990, The Practitioner.