Integrated analysis of single-cell and bulk RNA sequencing data identifies the characteristics of ferroptosis in lumbar disc herniation

Abstract

Lumbar disc herniation (LDH) is the most common condition associated with low back pain, and it adversely impacts individuals' health. Ferroptosis has recently emerged as a novel factor in the pathogenesis of LDH; however, the specific impacts of ferroptosis on LDH have not been fully elucidated. Ferroptosis-related differentially expressed genes (FRDEGs) were identified from the transcriptomic datasets of LDH. Gene set enrichment analysis (GSEA) was conducted to identify biological mechanism and related pathways. LASSO and SVM-RFE algorithms were applied to detect signature genes. Function of the signature gene was confirmed by RT-qPCR. The CIBERSORT algorithm was used to compare immune infiltration between LDH and normal samples. Correlation analysis between MYB and immune cells was analyzed using the Pearson method. Additionally, we used scRNA-seq to dissect cell clusters and cellular interactions. AUCell scoring was used to analyze the ferroptosis scores of different cell types. We found that MYB, a highly expressed ferroptosis-related gene, was associated with LDH By leveraging bioinformatics analysis. In immune infiltration analysis, the abundance of monocytes and macrophages varied significantly between the LDH group and disc spondylolisthesis (DS) group. MYB was correlated with most immune cells. GSEA revealed MYB was significantly enriched in immune-related pathways. Furthermore, scRNA-seq analysis revealed the presence of eight distinct cell types. AUCell analysis showed that macrophages had a high ferroptosis score. Cell trajectory analysis revealed that chondrocyte 1 was at the beginning of the trajectory, while calcification inhibiting chondrocytes and fibrochondrocytes accumulated along the middle and tail end of the trajectory, respectively. Cell-cell communication analysis identified chondrocyte 1 had an extensive communication network with other clusters and interacted with nucleus pulposus through collagen signaling pathway. Our analysis demonstrated that MYB may be a potential therapeutic target for LDH. This study provides a resource for the orthopedics community that will facilitate additional discoveries directedly toward understanding the pathogenesis process of LDH.

Keywords: Cell–cell communication; Ferroptosis; Lumbar disc herniation; Single-cell RNA sequencing; Transcriptome profiling analysis.