Abstract
The nucleus pulposus (NP) plays a vital role in intervertebral disc degeneration (IVDD). Previous studies have revealed cellular heterogeneity in NP tissue during IVDD progression. However, the cellular and molecular alterations of diverse cell clusters during IVDD remain to be fully elucidated. NP tissues were isolated from patients with different grades of IVDD undergoing discectomy, and then subjected to single-cell RNA sequencing (scRNA-seq). Cell subsets were identified based on unbiased clustering of gene expression profiles. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to determine the molecular features of diverse cell clusters. Monocle analysis was used to illustrate the differentiation trajectories of chondrocytes. Additionally, CellPhoneDB analysis revealed potential interactions between chondrocytes and other cells during IVDD. Based on the expression profiles of 47,610 individual cells, eight putative clusters including chondrocytes, endothelial cells, fibroblasts, macrophages, mural cells, osteoclasts, proliferating stromal cells and T cells were identified. The chondrocyte cluster was classified into three subsets, C1-C3, which were associated with stress-resistance, fibrosis and inflammatory responses, respectively. Pseudo-time trajectories suggested that chondrocytes gradually differentiated into fibroblasts during IVDD. Immune cells including cDC2s, macrophages and monocytes were identified. Further analysis showed that chondrocytes might communicate with immune cells via the MIF, TNFSF9, SPP1 and CCL4L2 signaling pathways. In addition, we found that invading endothelial cells might interact with chondrocytes through the COL4A1, CXCL12, VEGFA and SEMA3E signaling pathways. Our results reveal the cellular complexity and phenotypic characteristics of NP tissues at single-cell resolution, which will contribute to the in-depth investigation of preventative and regenerative strategies for IVDD.