Single-cell transcriptomics reveals heterogeneity and intercellular crosstalk in human intervertebral disc degeneration

Fibrocyte enrichment and myofibroblastic adaptation causes nucleus pulposus fibrosis and associates with disc degeneration severity

Fibrotic remodeling of nucleus pulposus (NP) leads to structural and mechanical anomalies of intervertebral discs that prone to degeneration, leading to low back pain incidence and disability. Emergence of fibroblastic cells in disc degeneration has been reported, yet their nature and origin remain elusive. In this study, we performed an integrative analysis of multiple single-cell RNA sequencing datasets to interrogate the cellular heterogeneity and fibroblast-like entities in degenerative human NP specimens. We found that disc degeneration severity is associated with an enrichment of fibrocyte phenotype, characterized by CD45 and collagen I dual positivity, and expression of myofibroblast marker α-smooth muscle actin. Refined clustering and classification distinguished the fibrocyte-like populations as subtypes in the NP cells - and immunocytes-clusters, expressing disc degeneration markers HTRA1 and ANGPTL4 and genes related to response to TGF-β. In injury-induced mouse disc degeneration model, fibrocytes were found recruited into the NP undergoing fibrosis and adopted a myofibroblast phenotype. Depleting the fibrocytes in CD11b-DTR mice in which myeloid-derived lineages were ablated by diphtheria toxin could markedly attenuate fibrous modeling and myofibroblast formation in the NP of the degenerative discs, and prevent disc height loss and histomorphological abnormalities. Marker analysis supports that disc degeneration progression is dependent on a function of CD45+COL1A1+ and αSMA+ cells. Our findings reveal that myeloid-derived fibrocytes play a pivotal role in NP fibrosis and may therefore be a target for modifying disc degeneration and promoting its repair.

Specific plasma biomarker signatures associated with patients undergoing surgery for back pain

BACKGROUND CONTEXT

Intervertebral disc degeneration (IDD) affects numerous people worldwide. The role of inflammation is increasingly recognized but remains incompletely resolved. Peripheral molecules could access neovascularized degenerated discs and contribute to the ongoing pathology.

PURPOSE

To assess a large array of plasma molecules in patients with IDD to identify biomarkers associated with specific spinal pathologies and prognostic biomarkers for the surgery outcome.

DESIGN

Prospective observational study combining clinical data and plasma measures.

PATIENT SAMPLE

Plasma samples were collected just before surgery. Extensive clinical data (age, sex, smoking status, Modic score, glomerular filtration rate, etc.) were extracted from clinical files from 83 patients with IDD undergoing spine surgery.

OUTCOME MEASURES

Recovery 2 months postsurgery as assessed by the treating neurosurgeon.

METHODS

Over 40 biological molecules were measured in patients' plasma using multiplex assays. Statistical analyses were performed to identify associations between biological and clinical characteristics (age, sex, Body Mass Index (BMI), smoking status, herniated disc, radiculopathy, myelopathy, stenosis, MODIC score, etc.) and plasma levels of biological molecules.

RESULTS

Plasma levels of Neurofilament Light chain (NfL) were significantly elevated in patients with myelopathy and spinal stenosis compared to herniated disc. Plasma levels of C- reactive protein (CRP), Neurofilament Light chain (NfL), and Serum Amyloid A (SAA) were negatively associated, while CCL22 levels were positively associated with an efficient recovery 2 months postsurgery.

CONCLUSIONS

Our results show that CRP and CCL22 plasma levels combined with the age of the IDD patient can predict the 2-month postsurgery recovery (Area Under the Curve [AUC]=0.883). Moreover, NfL could become a valuable monitoring tool for patients with spinal cord injuries.

SPP1-ITGα5/β1 Accelerates Calcification of Nucleus Pulposus Cells by Inhibiting Mitophagy via Ubiquitin-Dependent PINK1/PARKIN Pathway Blockade

Low back pain (LBP) caused by nucleus pulposus degeneration and calcification leads to great economic and social burden worldwide. Unexpectedly, no previous studies have demonstrated the association and the underlying mechanism between nucleus pulposus tissue degeneration and calcification formation. Secreted Phosphoprotein 1 (SPP1) exerts crucial functions in bone matrix mineralization and calcium deposition. Here, a novel function of SPP1 is reported, namely that it can aggravate nucleus pulposus cells (NPs) degeneration by negatively regulating extracellular matrix homeostasis. The degenerated NPs have a higher mineralization potential, which is achieved by SPP1. Mechanistically, SPP1 can accelerate the degeneration of nucleus pulposus cells by activating integrin α5β1 (ITGα5/β1), aggravating mitochondrial damage and inhibiting mitophagy. SPP1-ITGα5/β1 axis inhibits mitophagy by PINK1/PARKIN pathway blockade. In conclusion, SPP1 activates ITGα5/β1 to inhibit mitophagy, accelerates NPs degeneration, and induces calcification, thereby leading to intervertebral disc degeneration (IVDD) and calcification, identifying the potentially unknown mechanism and relationship between IVDD and calcification. Important insights are provided into the role of SPP1 in nucleus pulposus calcification in IVDD by inducing nucleus pulposus cell senescence through inhibition of mitophagy and may help develop potential new strategies for IVDD treatment.

Inactivation of Tnf-α/Tnfr signaling attenuates progression of intervertebral disc degeneration in mice

Background

Intervertebral disc degeneration (IVDD) is a major cause of low back pain (LBP), worsened by chronic inflammatory processes associated with aging. Tumor necrosis factor alpha (Tnf-α) and its receptors, Tnf receptor type 1 (Tnfr1) and Tnf receptor type 2 (Tnfr2), are upregulated in IVDD. However, its pathologic mechanisms remain poorly defined.

Methods

To investigate the role of Tnfr in IVDD, we generated global Tnfr1/2 double knockout (KO) mice and age-matched control C57BL/6 male mice, and analyzed intervertebral disc (IVD)-related phenotypes of both genotypes under physiological conditions, aging, and lumbar spine instability (LSI) model through histological and immunofluorescence analyses and μCT imaging. Expression levels of key extracellular matrix (ECM) proteins in aged and LSI mice, especially markers of cell proliferation and apoptosis, were evaluated in aged (21-month-old) mice.

Results

At 4 months, KO and control mice showed no marked differences of IVDD-related parameters. However, at 21 months of age, the loss of Tnfr expression significantly alleviated IVDD-like phenotypes, including a significant increase in height of the nucleus pulposus (NPs) and reductions of endplates (EPs) porosity and histopathological scores, when compared to controls. Tnfr deficiency promoted anabolic metabolism of the ECM proteins and suppressed ECM catabolism. Tnfr loss largely inhibited hypertrophic differentiation, and, in the meantime, suppressed cell apoptosis and cellular senescence in the annulus fibrosis, NP, and EP tissues without affecting cell proliferation. Similar results were observed in the LSI model, where Tnfr deficiency significantly alleviated IVDD and enhanced ECM anabolic metabolism while suppressing catabolism.

Conclusion

The deletion of Tnfr mitigates age-related and LSI-induced IVDD, as evidenced by preserved IVD structure, and improved ECM integrity. These findings suggest a crucial role of Tnf-α/Tnfr signaling in IVDD pathogenesis in mice. Targeting this pathway may be a novel strategy for IVDD prevention and treatment.

Cartilage Endplate-Targeted Engineered Exosome Releasing and Acid Neutralizing Hydrogel Reverses Intervertebral Disc Degeneration

Cartilage endplate cell (CEPC) and nucleus pulposus cell (NPC) inflammation are critical factors that contribute to intervertebral disc degeneration (IVDD). Recent evidence indicated that iron ion influx, reactive oxygen species (ROS), and the cGAS-STING pathway are involved in CEPC inflammatory degeneration. Moreover, cytokines produced by degenerating CEPCs and lactic acid accumulation within the microenvironment significantly contribute to NPC inflammation. Consequently, simultaneous alleviation of CEPC inflammation and correction of the acidic microenvironment are anticipated to reverse IVDD. Herein, CEPC-targeted engineered exosomes loaded with salvianolic acid A are incorporated into a CaCO3/chitosan hydrogel, forming a composite gel, CAP-sEXOs@Gel. Notably, CAP-sEXOs@Gel shows long local retention, realizes the slow release of CAP-sEXOs and specific uptake by CEPCs. After uptake by CEPCs, CAP-sEXOs reduce intracellular iron ion and ROS by inhibiting hypoxia-inducible factor-2α (HIF-2α)/TfR1 expression. Iron ion influx and ROS inhibition contribute to the maintenance of normal mitochondrial function and reduced mtDNA leakage, suppresing the cGAS-STING pathway. Additionally, the CaCO3 component of CAP-sEXOs@Gel neutralizes H+, thereby alleviating NPC inflammation. Collectively, this novel composite hydrogel demonstrates the ability to concurrently inhibit CEPC and NPC inflammation, thereby presenting a promising therapeutic approach for IVDD.

Integrated analysis of single-cell transcriptome and structural biology approach reveals the dynamics changes of NP subtypes and roles of Menaquinone in attenuating intervertebral disc degeneration

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.

Single-cell RNA sequencing reveals the MIF/ACKR3 receptor-ligand interaction between neutrophils and nucleus pulposus cells in intervertebral disc degeneration

OBJECTIVES

To unravel the heterogeneity and function of microenvironmental neutrophils during intervertebral disc degeneration (IDD).

METHODS

Single-cell RNA sequencing (scRNA-seq) was utilized to dissect the cellular landscape of neutrophils in intervertebral disc (IVD) tissues and their crosstalk with nucleus pulposus cells (NPCs). The expression levels of macrophage migration inhibitory factor (MIF) and ACKR3 in IVD tissues were detected. The MIF/ACKR3 axis was identified and its effects on IDD were investigated in vitro and in vivo.

RESULTS

We sequenced here 71520 single cells from 5 control and 9 degenerated IVD samples using scRNA-seq. We identified a unique cluster of neutrophils abundant in degenerated IVD tissues that highly expressed MIF and was functionally enriched in extracellular matrix organization (ECMO). Cell-to-cell communication analyses showed that this ECMO-neutrophil subpopulation was closely interacted with an effector NPCs subtype, which displayed high expression of ACKR3. Further analyses revealed that MIF was positively correlated with ACKR3 and functioned via directly binding to ACKR3 on effector NPCs. MIF inhibition attenuated degenerative changes of NPCs and extracellular matrix, which could be partially reversed by ACKR3 overexpression. Clinically, a significant correlation of high MIF/ACKR3 expression with advanced IDD grade was observed. Furthermore, we also found a positive association between MIF+ ECMO-neutrophil counts and ACKR3+ effector NPCs density as well as higher expression of the MIF/ACKR3 signaling in areas where these two cell types were neighbors.

CONCLUSIONS

These data suggest that ECMO-neutrophil promotes IDD progression by their communication with NPCs via the MIF/ACKR3 axis, which may shed light on therapeutic strategies.

Pip5k1γ promotes anabolism of nucleus pulposus cells and intervertebral disc homeostasis by activating CaMKII-Ampk pathway in aged mice

Degenerative disc disease (DDD) represents a significant global health challenge, yet its underlying molecular mechanisms remain elusive. This study aimed to investigate the role of type 1 phosphatidylinositol 4-phosphate 5-kinase (Pip5k1) in intervertebral disc (IVD) homeostasis and disease. All three Pip5k1 isoforms, namely Pip5k1α, Pip5k1β, and Pip5k1γ, were detectable in mouse and human IVD tissues, with Pip5k1γ displaying a highest expression in nucleus pulposus (NP) cells. The expression of Pip5k1γ was significantly down-regulated in the NP cells of aged mice and patients with severe DDD. To determine whether Pip5k1γ expression is required for disc homeostasis, we generated a Pip5k1γfl/fl; AggrecanCreERT2 mouse model for the conditional knockout of the Pip5k1γ gene in aggrecan-expressing IVD cells. Our findings revealed that the conditional deletion of Pip5k1γ did not affect the disc structure or cellular composition in 5-month-old adult mice. However, in aged (15-month-old) mice, this deletion led to several severe degenerative disc defects, including decreased NP cellularity, spontaneous fibrosis and cleft formation, and a loss of the boundary between NP and annulus fibrosus. At the molecular level, the absence of Pip5k1γ reduced the anabolism of NP cells without markedly affecting their catabolic or anti-catabolic activities. Moreover, the loss of Pip5k1γ significantly dampened the activation of the protective Ampk pathway in NP cells, thereby accelerating NP cell senescence. Notably, Pip5k1γ deficiency blunted the effectiveness of metformin, a potent Ampk activator, in activating the Ampk pathway and mitigating lumbar spine instability (LSI)-induced disc lesions in mice. Overall, our study unveils a novel role for Pip5k1γ in promoting anabolism and maintaining disc homeostasis, suggesting it as a potential therapeutic target for DDD.

Investigating the characteristics of mild intervertebral disc degeneration at various age stages using single-cell genomics

Introduction: Intervertebral disc degeneration often occurs in the elderly population, but in recent years, there has been an increasing incidence of disc degeneration in younger individuals, primarily with mild degeneration. Methods: In order to explore the underlying mechanisms of disc degeneration in both young and aging individuals, we collected four types of nucleus pulposus (NP) single-cell sequencing samples for analysis based on Pfirrmann grading: normal-young (NY) (Grade I), normal-old (NO) (Grade I), mild degenerative-young (MY) (Grade II-III), and mild degenerative-old (MO) (Grade II-III). Results: We found that most NP cells in NO and MY samples exhibited oxidative stress, which may be important pathogenic factors in NO and MY groups. On the other hand, NP cells in MO group exhibited endoplasmic reticulum stress. In terms of inflammation, myeloid cells were mainly present in the degenerative group, with the MY group showing a stronger immune response compared to the MO group. Interestingly, dendritic cells in the myeloid lineage played a critical role in the process of mild degeneration. Discussion: Our study investigated the molecular mechanisms of intervertebral disc degeneration from an age perspective, providing insights for improving treatment strategies for patients with disc degeneration at different age groups.

Wharton's jelly-derived multifunctional hydrogels: New tools to promote intervertebral disc regeneration in vitro and ex vivo

The degeneration of intervertebral disc (IVD) is a disease of the entire joint between two vertebrae in the spine caused by loss of extracellular matrix (ECM) integrity, to date with no cure. The various regenerative approaches proposed so far have led to very limited successes. An emerging opportunity arises from the use of decellularized ECM as a scaffolding material that, directly or in combination with other materials, has greatly facilitated the advancement of tissue engineering. Here we focused on the decellularized matrix obtained from human umbilical cord Wharton's jelly (DWJ) which retains several structural and bioactive molecules very similar to those of the IVD ECM. However, being a viscous gel, DWJ has limited ability to retain ordered structural features when considered as architecture scaffold. To overcome this limitation, we produced DWJ-based multifunctional hydrogels, in the form of 3D millicylinders containing different percentages of alginate, a seaweed-derived polysaccharide, and gelatin, denatured collagen, which may impart mechanical integrity to the biologically active DWJ. The developed protocol, based on a freezing step, leads to the consolidation of the entire polymeric dispersion mixture, followed by an ionic gelation step and a freeze-drying process. Finally, a porous, stable, easily storable, and suitable matrix for ex vivo experiments was obtained. The properties of the millicylinders (Wharton's jelly millicylinders [WJMs]) were then tested in culture of degenerated IVD cells isolated from disc tissues of patients undergoing surgical discectomy. We found that WJMs with the highest percentage of DWJ were effective in supporting cell migration, restoration of the IVD phenotype (increased expression of Collagen type 2, aggrecan, Sox9 and FOXO3a), anti-inflammatory action, and stem cell activity of resident progenitor/notochordal cells (increased number of CD24 positive cells). We are confident that the DWJ-based formulations proposed here can provide adequate stimuli to the cells present in the degenerated IVD to restart the anabolic machinery.

Silencing of secreted phosphoprotein 1 attenuates sciatic nerve injury-induced neuropathic pain: Regulating extracellular signal-regulated kinase and neuroinflammatory signaling pathways

BACKGROUND

Neuropathic pain (NP) is a chronic pathological pain that affects the quality of life and is a huge medical burden for affected patients. In this study, we aimed to explore the effects of secreted phosphoprotein 1 (SPP1) on NP.

METHODS

We established a chronic constriction injury (CCI) rat model, knocked down SPP1 via an intrathecal injection, and/or activated the extracellular signal-regulated kinase (ERK) pathway with insulin-like growth factor 1 (IGF-1) treatment. Pain behaviors, including paw withdrawal threshold (PWT), paw withdrawal latency (PWL), lifting number, and frequency, were assessed. After sacrificing rats, the L4-L5 dorsal root ganglion was collected. Then, SPP1 levels were determined using quantitative polymerase chain reaction (qPCR) and western blot analysis. The levels of interleukin (IL)-1β, tumor necrosis factor (TNF)-α, IL-6, IL-10, epidermal growth factor (EGF), vascular endothelial growth factor (VEGF), and transforming growth factor (TGF)-β were determined using qPCR and enzyme-linked immunosorbent assay. The levels of ERK pathway factors were determined via western blot analysis.

RESULTS

We found that CCI decreased PWT and PWL, increased the lifting number and frequency, and upregulated SPP1 levels. The loss of SPP1 reversed these CCI-induced effects. Additionally, CCI upregulated IL-1β, TNF-α, IL-6, EGF, and VEGF levels, downregulated TGF-β levels, and activated the ERK pathway, while silencing of SPP1 abrogated these CCI-induced effects. Moreover, IGF-1 treatment reversed the effects of SPP1 loss.

CONCLUSIONS

The data indicate that silencing SPP1 attenuates NP via inactivation of the ERK pathway, suggesting that SPP1 may be a promising target for NP treatment.

Transcriptional profiling of human cartilage endplate cells identifies novel genes and cell clusters underlying degenerated and non-degenerated phenotypes

BACKGROUND

Low back pain is a leading cause of disability worldwide and is frequently attributed to intervertebral disc (IVD) degeneration. Though the contributions of the adjacent cartilage endplates (CEP) to IVD degeneration are well documented, the phenotype and functions of the resident CEP cells are critically understudied. To better characterize CEP cell phenotype and possible mechanisms of CEP degeneration, bulk and single-cell RNA sequencing of non-degenerated and degenerated CEP cells were performed.

METHODS

Human lumbar CEP cells from degenerated (Thompson grade ≥ 4) and non-degenerated (Thompson grade ≤ 2) discs were expanded for bulk (N=4 non-degenerated, N=4 degenerated) and single-cell (N=1 non-degenerated, N=1 degenerated) RNA sequencing. Genes identified from bulk RNA sequencing were categorized by function and their expression in non-degenerated and degenerated CEP cells were compared. A PubMed literature review was also performed to determine which genes were previously identified and studied in the CEP, IVD, and other cartilaginous tissues. For single-cell RNA sequencing, different cell clusters were resolved using unsupervised clustering and functional annotation. Differential gene expression analysis and Gene Ontology, respectively, were used to compare gene expression and functional enrichment between cell clusters, as well as between non-degenerated and degenerated CEP samples.

RESULTS

Bulk RNA sequencing revealed 38 genes were significantly upregulated and 15 genes were significantly downregulated in degenerated CEP cells relative to non-degenerated cells (|fold change| ≥ 1.5). Of these, only 2 genes were previously studied in CEP cells, and 31 were previously studied in the IVD and other cartilaginous tissues. Single-cell RNA sequencing revealed 11 unique cell clusters, including multiple chondrocyte and progenitor subpopulations with distinct gene expression and functional profiles. Analysis of genes in the bulk RNA sequencing dataset showed that progenitor cell clusters from both samples were enriched in "non-degenerated" genes but not "degenerated" genes. For both bulk- and single-cell analyses, gene expression and pathway enrichment analyses highlighted several pathways that may regulate CEP degeneration, including transcriptional regulation, translational regulation, intracellular transport, and mitochondrial dysfunction.

CONCLUSIONS

This thorough analysis using RNA sequencing methods highlighted numerous differences between non-degenerated and degenerated CEP cells, the phenotypic heterogeneity of CEP cells, and several pathways of interest that may be relevant in CEP degeneration.

Association of the FGF4L2 retrogene with fibrocartilaginous embolic myelopathy in dogs

BACKGROUND

Fibrocartilaginous embolic myelopathy (FCE) is a well-documented condition in dogs although rarely reported in chondrodystrophic breeds. Genetic associations have not been defined.

OBJECTIVES

Define the association of the chondrodystrophy-associated FGF4L2 retrogene with histopathologically confirmed cases of FCE.

ANIMALS

Ninety-eight dogs with a histopathologic diagnosis of FCE.

METHODS

Retrospective multicenter study. Dogs were genotyped for the FGF4L2 and FGF4L1 retrogenes using DNA extracted from formalin-fixed, paraffin-embedded tissue. Associations between breed, FCE and retrogene status were investigated with reference to a hospital population and known breed and general population allele frequencies.

RESULTS

FGF4L2 genotype was defined in 89 FCE cases. Fibrocartilaginous embolic myelopathy was present in 22 dogs from FGF4L2-segregating breeds with allele frequencies of ≥5%; however, all dogs were wild type. Two Labrador retrievers with FCE carried FGF4L2 alleles. Frequency of the FGF4L2 allele was significantly (P < .001) and negatively associated with FCE relative to predicted hospital-population dogs. FCE was overrepresented in Boxer, Great Dane, Yorkshire Terrier, Bernese Mountain Dog, Miniature Schnauzer, Rottweiler, and Shetland Sheepdog breeds.

CONCLUSIONS AND CLINICAL IMPORTANCE

Study data based on genotypically and histopathologically defined cases support the historical observation that FCE is uncommon in chondrodystrophic dog breeds. FGF4 plays an important role in angiogenesis and vascular integrity; anatomical studies comparing chondrodystrophic and non-chondrodystrophic dogs might provide insight into the pathogenesis of FCE.

Gene Expression Profiles Perturbed by Injury to the Mouse Intervertebral Disc

OBJECTIVES

Back pain subsequent to intervertebral disc injury is a common clinical problem. Previous work examining early molecular changes post injury mainly used a candidate marker approach. In this study, gene expression in the injured and intact mouse tail intervertebral discs was determined with a nonbiased whole transcriptome approach.

DESIGN

Mouse tail intervertebral disc injury was induced by a needle puncture. Whole murine transcriptome was determined by RNASeq. Transcriptomes of injured intervertebral discs were compared with those of intact controls by bioinformatic methods.

RESULTS

Among the 18,078 murine genes examined, 592 genes were differentially expressed ( P.adj < 0.01). Novel genes upregulated in injured compared with intact intervertebral discs included Chl1, Lum , etc. Ontology study of upregulated genes revealed that leukocyte migration was the most enriched biological process, and network analysis showed that Tnfa had the most protein-protein interactions. Novel downregulated genes in the injured intervertebral discs included 4833412C05Rik , Myoc , etc . The most enriched downregulated pathways were related to cytoskeletal organization.

CONCLUSIONS

Novel genes highly regulated after disc injury were identified with an unbiased approach; they may serve as biomarkers of injury and response to treatments in future experiments. Enriched biological pathways and molecules with high numbers of connections may be targets for treatments after injury.