Intervertebral disc degeneration (IDD) is a progressive and chronic disease, the mechanisms have been studied extensively as a whole, while the cellular heterogeneity of cells in nucleus pulposus (NP) tissues remained controversial for a long time. This study conducted integrated analysis through single-cell sequencing analysis, weighted gene co-expression network analysis (WGCNA), and differential expression analysis, to systematically decipher the longitudinal alterations of distinct NP subtypes, and also analyzed the most essential genes in the development of IDD. Then, this study further conducted structural biology method to discover the potential lead compounds through a suite of advanced approaches like high-throughput screening (HTVS), pharmaceutical characteristics assessment, CDOCKER module as well as molecular dynamics simulation, etc., aiming to ameliorate the progression of IDD. Totally 5 NP subpopulations were identified with distinct biological functions based on their unique gene expression patterns. The predominant dynamics changes mainly involved RegNPs and EffNPs, the RegNPs were mainly aggregated in normal NP tissues and drastically decreased in degenerative NP, while EffNPs, as pathogenic subtype, exhibited opposite phenomenon. Importantly, this study further reported the essential roles of Menaquinone in alleviating degenerative NP cells for the first time, which could provide solid evidence for the application of nutritional therapy in the treatment of IDD. This study combined scRNA-seq, bulk-RNA seq and HTVS techniques to systematically decipher the longitudinal changes of NP subtypes during IDD. EffNPs were considered to be 'chief culprit' in IDD progression, while the novel natural drug Menaquinone could reverse this phenomenon.Communicated by Ramaswamy H. Sarma.
Keywords: Intervertebral disc degeneration; high-throughput virtual screening analysis; longitudinal approach; single-cell transcriptome landscape; trajectory lineage analysis.