Abstract
Ferroptosis and immune infiltration play an important role in the pathogenesis of intervertebral disc degeneration (IDD). However, there is still a lack of comprehensive analysis on the interaction between ferroptosis-related genes (FRGs) and immune microenvironment in IDD patients. Therefore, this study aims to explore the correlation between FRGs characteristics and immune infiltration in the progression of IDD. The expression profiles (GSE56081 and GSE70362) and FRGs were downloaded from the comprehensive gene expression omnibus (GEO) and FerrDb database, respectively, and the differences were analyzed using R. The intersection of IDD related differential genes (DEGs) and FRGs was taken as differentially expressed FRGs (DE-FRGs) and GO and KEGG enrichment analysis was conducted. Then, we used least absolute shrinkage and selection operator (LASSO) regression algorithm and support vector machine (SVM) algorithm to screen feature genes and draw ROC curve judge the diagnostic value of key DE-FRGs. Then CIBERSORT algorithm is used to evaluate the infiltration of immune cells and analyze the correlation between key DE-FRGs and immune infiltration. Based on the analysis results, we conducted single gene GSEA analysis on key DE-FRGs. RT-PCR and immunohistochemistry further verified the clinical value of the results of biochemical analysis and screening. Seven key DE-FRGs were screened, including the upregulated genes NOX4 and PIR, and the downregulated genes TIMM9, ATF3, ENPP2, FADS2 and TFAP2A. Single gene GSEA analysis further elucidates the role of DE-FRGs in IDD associated with ferroptosis. Correlation analysis showed that seven key DE-FRGs were closely related to immune infiltration in the development of IDD. Finally, RT-PCR and immunohistochemical staining showed that NOX4, ENPP2, FADS2 and TFAP2A were statistically significant differences. In this study, we explored the connection between ferroptosis related characteristics and immune infiltration in IDD, and confirmed that NOX4, ENPP2, FADS2, and TFAP2A may become biomarkers and potential therapeutic targets for IDD.