Single-cell RNA-seq analysis identifies unique chondrocyte subsets and reveals involvement of ferroptosis in human intervertebral disc degeneration

Summary Objective Nucleus pulposus (NP) plays a central role in disc degeneration pathogenesis, however, as a heterogeneous tissue, cell subsets in NP and their corresponding biological process in intervertebral disc degeneration (IVDD) are unreported. Method Nucleus pulposus were isolated from normal control and IVDD, and then subjected to single-cell RNA sequencing (scRNA-seq). Unsupervised clustering of the cells based on the gene expression profiles using the Seurat package and passed to tSNE for clustering visualization. Rat model of disc degeneration was built to validate the pathways identified by scRNA-Seq. Results Seven chondrocyte subsets were revealed in NP based on differential gene expression, among which 4 subsets (C1–C4) were reported for the first time. Furthermore, GO and KEGG analyses discovered that ferroptosis pathways were enriched. Rat model of disc degeneration was built (n = 6/group, control vs. model) to validate the pathways identified by scRNA-Seq. Iron levels of NP were significantly higher in model group than control group (means 0.712 vs. 0.248, respectively, mg/gpro, p = 0.0026), and the levels of Heme Oxygenase 1 (HO-1) were also elevated in model group (means 14.33 vs. 5.16 IOD, respectively, p = 0.0002). However, the levels of ferritin light chain (FTL) were significantly decreased in model group compared to control group (means 26.17 vs. 9.00 FTL+ cell number, respectively, p = 0.0011). Conclusions Novel chondrocyte subsets in nucleus pulposus were discovered through scRNA-Seq, which provided novel insight to understand the pathological change during the development of IVDD. Ferroptosis participated in disc degeneration pathogenesis and it might serve as a new target for intervening IVDD.