Spatially multicellular variability of intervertebral disc degeneration by comparative single-cell analysis

Previous studies have revealed cellular heterogeneity in intervertebral discs (IVDs). However, the cellular and molecular alteration patterns of cell populations during degenerative progression remain to be fully elucidated. To illustrate the cellular and molecular alteration of cell populations in intervertebral disc degeneration (IDD), we perform single cell RNA sequencing on cells from four anatomic sites of healthy and degenerative goat IVDs. EGLN3(+) StressCs, TGFBR3(+) HomCs and GPRC5A(+) RegCs exhibit the characteristics associated with resistance to stress, maintaining homeostasis and repairing, respectively. The frequencies and signatures of these cell clusters fluctuate with IDD. Notably, the chondrogenic differentiation programme of PROCR+ progenitor cells is altered by IDD, while notochord cells turn to stemness exhaustion. In addition, we characterise CAV1(+) endothelial cells that communicate with chondrocytes through multiple signalling pathways in degenerative IVDs. Our comprehensive analysis identifies the variability of key cell clusters and critical regulatory networks responding to IDD, which will facilitate in-depth investigation of therapeutic strategies for IDD.