Identification of Osteoarthritis Inflamm-Aging Biomarkers by Integrating Bioinformatic Analysis and Machine Learning Strategies and the Clinical Validation

Objective To identify inflamm-aging related biomarkers in osteoarthritis (OA).

Methods Microarray gene profiles of young and aging OA patients were obtained from the Gene Expression Omnibus (GEO) database and aging-related genes (ARGs) were obtained from the Human Aging Genome Resource (HAGR) database. The differentially expressed genes of young OA and older OA patients were screened and then intersected with ARGs to obtain the aging-related genes of OA. Enrichment analysis was performed to reveal the potential mechanisms of aging-related markers in OA. Three machine learning methods were used to identify core senescence markers of OA and the receiver operating characteristic (ROC) curve was used to assess their diagnostic performance. Peripheral blood mononuclear cells were collected from clinical OA patients to verify the expression of senescence-associated secretory phenotype (SASP) factors and senescence markers.

Results A total of 45 senescence-related markers were obtained, which were mainly involved in the regulation of cellular senescence, the cell cycle, inflammatory response, etc. Through the screening with the three machine learning methods, 5 core senescence biomarkers, including FOXO3, MCL1, SIRT3, STAG1, and S100A13, were obtained. A total of 20 cases of normal controls and 40 cases of OA patients, including 20 cases in the young patient group and 20 in the elderly patient group, were enrolled. Compared with those of the young patient group, C-reactive protein (CRP), interleukin (IL)-6, and IL-1β levels increased and IL-4 levels decreased in the elderly OA patient group (P<0.01); FOXO3, MCL1, and SIRT3 mRNA expression decreased and STAG1 and S100A13 mRNA expression increased (P<0.01). Pearson correlation analysis demonstrated that the selected markers were associated with some indicators, including erythrocyte sedimentation rate (ESR), IL-1β, IL-4, CRP, and IL-6. The area under the ROC curve of the 5 core aging genes was always greater than 0.8 and the C-index of the calibration curve in the nomogram prediction model was 0.755, which suggested the good calibration ability of the model.

Conclusion FOXO3, MCL1, SIRT3, STAG1, and S100A13 may serve as novel diagnostic biomolecular markers and potential therapeutic targets for OA inflamm-aging.