Identification of sulfur metabolism-related gene signature in osteoarthritis and TM9SF2’s sustenance effect on M2 macrophages' phagocytic activity

Abstract

Abstract Background Osteoarthritis (OA) is a chronic and low-grade inflammatory disease associated with metabolism disorder and multiple cell death types in the synovial tissues. Sulfur metabolism has not been studied in OA. Methods First, we calculated the single sample gene set enrichment analysis score of sulfur metabolism-associated annotations (i.e., cysteine metabolism process, regulation of sulfur metabolism process, and disulfidptosis) between healthy and synovial samples from patients with OA. Sulfur metabolism-related differentially expressed genes (DEGs) were analyzed in OA. Least absolute shrinkage and selection operator COX regression were used to identify the sulfur metabolism-associated gene signature for diagnosing OA. Correlation and immune cell deconvolution analyses were used to explore the correlated functions and cell specificity of the signature gene, TM9SF2. TM9SF2’s effect on the phagocytosis of macrophages M2 was analyzed by coculturing macrophages with IgG-coated beads or apoptotic Jurkat cells. Results A diagnostic six gene signature (i.e., MTHFD1, PDK4, TM9SF2, POU4F1, HOXA2, NCKAP1) was identified based on the ten DEGs, validated using GSE12021 and GSE1919 datasets. TM9SF2 was upregulated in the synovial tissues of OA at both mRNA and protein levels. The relationship between TM9SF2 and several functional annotations, such as antigen processing and presentation, lysosome, phagosome, Fcγ-mediated phagocytosis, and tyrosine metabolism, was identified. TM9SF2 and macrophages M2 were significantly correlated. After silencing TM9SF2 in THP-1-derived macrophages M2, a significantly reduced phagocytosis and attenuated activation of PLC-γ1 were observed. Conclusion A sulfur metabolism-associated six-gene signature for OA diagnosis was constructed and upregulation of the phagocytosis-associated gene, TM9SF2, was identified. The findings are expected to deepen our understanding of the molecular mechanism underlying OA development and be used as potential therapeutic targets.