Bioinformatics and Integrative Experimental Method to Identifying and Validating Co-Expressed Ferroptosis-Related...

Osteoarthritis (OA) is the most common joint disease worldwide and is the primary cause of disability and chronic pain in older adults.Ferroptosis is a type of programmed cell death characterized by aberrant iron metabolism and reactive oxygen species accumulation; however, its role in OA is not known.

To identify ferroptosis markers co-expressed in articular cartilage and synovium samples from patients with OA, in silico analysis was performed.Signature genes were analyzed and the results were evaluated using a ROC curve prediction model.The biological function, correlation between Signature genes, immune cell infiltration, and ceRNA network analyses were performed. Signature genes and ferroptosis phenotypes were verified through in vivo animal experiments and clinical samples. The expression levels of non-coding RNAs in samples from patients with OA were determined using qRT-PCR. ceRNA network analysis results were confirmed using dual-luciferase assays.

JUN, ATF3, and CDKN1A were identified as OA- and ferroptosis-associated signature genes. GSEA analysis demonstrated an enrichment of these genes in immune and inflammatory responses, and amino acid metabolism. The CIBERSORT algorithm showed a negative correlation between T cells and these signature genes in the cartilage, and a positive correlation in the synovium. Moreover, RP5-894D12.5 and FAM95B1 regulated the expression of JUN, ATF3, and CDKN1A by competitively binding to miR-1972, miR-665, and miR-181a-2-3p. In vivo, GPX4 was downregulated in both OA cartilage and synovium; however, GPX4 and GSH were downregulated, while ferrous ions were upregulated in patient OA cartilage and synovium samples, indicating that ferroptosis was involved in the pathogenesis of OA. Furthermore, JUN, ATF3, and CDKN1A expression was downregulated in both mouse and human OA synovial and cartilage tissues. qRT-PCR demonstrated that miR-1972, RP5-894D12.5, and FAM95B1 were differentially expressed in OA tissues. Targeted interactions between miR-1972 and JUN, and a ceRNA regulatory mechanism between RP5-894D12.5, miR-1972, and JUN were confirmed by dual-luciferase assays.

This study identified JUN, ATF3, and CDKN1A as possible diagnostic biomarkers and therapeutic targets for joint synovitis and OA. Furthermore, our finding indicated that RP5-894D12.5/miR-1972/JUN was a potential ceRNA regulatory axis in OA, providing an insight into the connection between ferroptosis and OA.