Pengcheng Liu
Xiang Ren
Beiting Zhang1†
Song Guo
Qiang Fu- 1Department of Orthopedics Surgery, Shanghai General Hospital, Shanghai Jiao Tong University, Shanghai, China
- 2Interdisciplinary Research Center on Biology and Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai, China
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.
Introduction
Low back pain (LBP) is a common clinical condition that has become increasingly prevalent. Statistical reports indicate that the global incidence rate of LBP ranges between 13.1% and 28.5% (Maher et al., 2017). Among the various identified factors, intervertebral disc degeneration (IVD) stands out as a primary contributor to LBP, emerging as a leading cause of disability worldwide and exerting significant strain on healthcare systems and economies (Chou, 2021; Knezevic et al., 2021; Yang et al., 2024). To date, most treatments of IVD degeneration are limited to invasive surgical interventions, such as disc replacement and spinal fusion, or pain management, which often fail to restore the integrity and normal physiological function of degenerated intervertebral discs (Knezevic et al., 2017; Meisel et al., 2019).
The intervertebral disc, a sealed structure situated between the vertebral bodies in the human spine, comprises cartilage plates, fibrous rings, and a nucleus pulposus (NP) (Gantenbein et al., 2023), which help sustain the natural flexibility and normal biomechanical functioning of the spine. Researchers have established that dysfunction in any component of the IVD leads to intervertebral disc degeneration (IVDD), with NP tissue degeneration recognized as the primary cause (Lama et al., 2022). IVDD is a multifactorial process with aging being a significant risk factor (Ngo et al., 2017; Iwakiri et al., 2023; Novais et al., 2024). However, disc degeneration can also manifest independently of the aging process, making it challenging to distinguish between physiological disc aging and degeneration (Soufi et al., 2023).
In recent years, single-cell RNA sequencing (scRNA-seq) has gained prominence as a powerful method for evaluating single-cell gene expression, identifying heterogeneous cell subpopulations within tissues, unveiling rare subpopulations, delineating cellular heterogeneity, discovering new markers, and predicting developmental trajectories. It offers insights into physiological and pathological processes and has been increasingly applied in IVD research. Several studies have identified novel phenotypic biomarkers of NP cells, significantly enhancing our understanding of NP cell biology and development (Fernandes et al., 2020; Ling et al., 2021; Panebianco et al., 2021; Zhang et al., 2021; Cherif et al., 2022; Han et al., 2022; RNA, 2022; Guo et al., 2023a; Chen et al., 2024). However, the heterogeneity of IVDD among young and elderly individuals remains to be fully elucidated.
To address this gap, we collected four types of nucleus pulposus (NP) single-cell sequencing data 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). We conducted in-depth gene ontology term enrichment and pathway analyses and predicted developmental trajectories for each type. Furthermore, we validated our findings through histology and immunofluorescence. The objective of this study was to explore the underlying mechanisms of disc degeneration in both young and aging individuals. These investigations aim to clarify the mechanisms of IVDD from an age perspective and offer insights for improving treatment strategies for patients with disc degeneration across different age groups.
Methods
Patients and tissue samples
The Institutional Review Board of the Shanghai General Hospital approved this study. Degenerative NP tissues were collected from patients undergoing surgery due to IVDD (Supplementary Figure S6). Healthy tissues were collected from patients who underwent surgery for idiopathic scoliosis or vertebral fracture. All samples were trimmed in order to preserve the nucleus pulposus tissue while removing the fibrous annulus. Before sample collection, the degree of IVD degeneration was evaluated according to the Pfirrmann grading system (Pfirrmann et al., 2001) (Supplementary Table S1).
Histology and immunofluorescence
IVD samples from the sample were fixed for 48 h using 4% paraformaldehyde (PFA), dehydrated and embedded in paraffin. Five-micron thick sections were cut from the paraffin blocks. For histological assessment, the paraffin sections were de-paraffinized in graded xylene, rehydrated in graded alcohol solutions, and washed and stained with hematoxylin and eosin (H&E) staining (Servicebio, Wuhan, China) per the manufacturer’s protocols.
For immunofluorescence staining, prepared paraffin sections were incubated in antigen retrieval buffer (Roche Diagnostics) at 37°C for 30 min. After cooling to RT, the sections were washed with PBS three times for 5 min each, followed by permeabilized for 10 min with 0.1% Triton X-100. Then the sections were washed and blocked for 1 h with 1% PBS/BSA and incubated with primary antibodies (anti-SOD2, anti-TOMM20, anti-CD68, anti-Collegen2, anti-HIF1α) in PBS/BSA at 4°C overnight. Afterwards secondary antibodies conjugated with Alexa Fluor 488 and/or 594 were used to incubate cells for 1 h at room temperature. The nucleus was labelled with 10 μg/mL DAPI solution. Pictures were captured using section scanner (Pannoramic MIDI). Images were analyzed using CaseViewer software (C.V2.4).
Single-cell RNA sequencing
Single cell sequencing data was downloaded from NCBI GEO with accession number GSE233666 (Guo et al., 2023b), GSE205535 (RNA, 2022), GSE189916 (Jiang et al., 2022) and CNGBdb with accession number CNP0002664 (Han et al., 2022). In short, the downloaded data was used to construct a Seurat object using the CreateSeuratObject function in Seurat package (v4.3.0) (Stuart et al., 2019), and the DoubletFinder package (v2.0.3) (McGinnis et al., 2019) was utilized to remove doublets. Cells were filtered by gene counts between 200 and 7,500 and UMI counts over 1,000. Number of genes detected per UMI are higher than 0.8 to remove the outlier cells which have less complex RNA species like red blood cells. Cells with over 15% mitochondrial content were removed. For dimension reduction and clustering, we utilized functions from Seurat package. To ensure uniform gene expression measurement, normalization and scaling were performed using NormalizeData and ScaleData methods. We conducted PCA on the top 2,000 variable genes, which were identified using FindVariableFeatures. The batch effects were eliminated in all datasets using the default parameters of the Harmony package (v0.1.1) (Korsunsky et al., 2019). The cells were subsequently classified into clusters using FindNeighbors and FindClusters with the top 20 principle components and a resolution parameter of 0.1. The DimPlot function was employed to generate a 2D Uniform Manifold Approximation and Projection (UMAP) plot using the “UMAP” reduction method. Cell-type identity of each cluster was determined through manual annotation using the expression of canonical markers. DotPlot, FeaturePlot, and Vlnplot functions were utilized to generate dot plots, feature plots, and violin plots, respectively, to visualize the expression of markers used for cell type identification. Differentially expressed genes (DEGs) within each celltype were identified using Seurat’s FindMarkers function, employing a Wilcoxon likelihood ratio test with default parameters. Genes expressed in over 10% of cells within a cluster and exhibiting an average log-fold change greater than 0.25 were considered as DEGs. The GO enrichment analysis with DEGs was conducted using the clusterProfiler package. The strength of signaling pathways was calculated using the AUCell package (v1.20.2) (Deng et al., 2021).
Pseudotime analysis
To determine the cell lineage of subclusters, Slingshot (Street et al., 2018) was used to find lineage by fitting minimum spanning trees based on clusters through the getLineages function. After constructing the smooth lineages, we inferred the pseudo time of each cell by fitting the principal curves through the getCurves function.
Cell communication analysis
The cell-cell communication network was identified by the CellChat R package as previously described (Jin et al., 2021). The cellular communications between different groups in IVD were analyzed. Briefly, we construct cellchat object base on the expression matrix and meta information from Seurat by using createCellChat function. Then, merge cellchat object of different grade sample by mergeCellChat function. After these, we use compareInteractions to compare different interactions between different grades and net Visual_diffInteraction function to visualize results.
CellChat (v1.1.3) was applied for ligand–receptor analysis. The normalized counts and cell-type annotations for each cell were imputed into CellChat to determine the potential ligand–receptor pairs. Interaction pairs with p-value > 0.05 were filtered out from further analysis. The interaction of myeloid cells and NP cells was analyzed. Selected specific pairs were plotted in CellChat with default parameters.
Results
Single-cell sequencing analysis reveals heterogeneity of intervertebral disc cells during degeneration and aging processes
To investigate the transcriptomic changes in the intervertebral disc during normal and mild degeneration in both young and elderly individuals, capturing potential initiating factors of intervertebral disc degeneration (IDD), we collected a total of 10 nucleus pulposus (NP) single-cell sequencing datasets (Supplementary Table S2) based on the Pfirrmann grading. The datasets included samples from normal-young (Grade I, NY), normal-old (Grade I, NO), mild degenerative-young (Grade II-III, MY), and mild degenerative-old (Grade II-III, MO). After quality control processes such as filtering out low-quality cells, doublets, and batch effects, a total of 50,077 high-quality NP cells were obtained for downstream analysis (Figure 1A).
Figure 1. Single-cell RNA sequencing analysis of human IVD cells. (A) Schematic workflow of sc-RNA sequencing analysis and experimental validation. (B) Uni-form Mani-fold Approximation and Projection (UMAP) showing the single-cell transcriptomic landscape of human IVD cells. Six clusters were visualized. (C) Dot plot of marker genes for distinct cell types. Color scale indicates the mean of normalized expression of marker genes in each cell type, and dot size is proportional to the percentage of cells within each cell cluster expressing the marker genes. (D) Distribution of cell clusters (defined in B) in different groups. NY, normal-young (Grade I), NO, normal-old (Grade I), MO, mild degenerative-young (Grade II-III), MY, mild degenerative-old (Grade II-III). (E) Same UMAP visualization as the one in B but split by different groups.
Based on the classical markers, we have initially divided the overall cells into six celltypes (Figures 1B, C), which are NP cell (COL2A1, MIA, ACAN), nucleus pulposus progenitor cell (NPPC) (COL1A1, PDGFRA, PRRX1), division-NPPC (MKI67), myeloid cell (CD68, CSF1R, CD163), endothelial cell (PECAM1, VWF, CD34), and T cell (CD3D, CD3E, CD3G). We found that NP cell were the main components in all groups, with no significant changes between groups. The proportion of NPPC decreased with aging in both normal and degenerated samples, while immune cells were mainly present in degenerated samples. Moreover, MY group had higher levels of immune cell infiltration than the MO group, suggesting that samples from younger patients exhibited more immune cell infiltration in mild degenerative disc disease (Figures 1D, E).
Oxidative stress is an important pathogenic factor in the NP cell of both NO and MY groups
To further understand the heterogeneity within total NP cells (including NPPC, division-NPPC, and NP cell in Figure 1B), we conducted subpopulation analysis. Using Uniform Manifold Approximation and Projection (UMAP) dimensionality reduction analysis, we divided the total NP cells into 10 clusters (Supplementary Figure S1). GO enrichment analysis suggested significant enrichment of the oxidative phosphorylation pathway in cells from cluster 0 and cluster 5. Since the normal intervertebral disc is avascular tissue, the internal environment is characterized by hypoxia. Studies have demonstrated that HIF1A plays a crucial role in enabling NP cells to adapt to hypoxic conditions and in inhibiting the production of reactive oxygen species (ROS). Conversely, normoxic conditions have been shown to result in a downregulation of HIF1A expression in NP cells, which in turn leads to an increase in ROS levels. Furthermore, the reduction of ROS in NP cells under hypoxic conditions is dependent on HIF1A, emphasizing the critical function of HIF1A in the typical survival of NP cells under hypoxic conditions (Silagi et al., 2021; Madhu et al., 2023; Yang et al., 2023; Song et al., 2024). Therefore, combining the expression pattern of HIF1A with the results of GO enrichment analysis, we divided the total NP cells into 8 celltypes, namely, HIF1Ahigh_NP1-2, HIF1Alow_NP1-3, and NPPC1-3 (Figures 2A, B; Supplementary Figure S2). We found that in the NY group, HIF1Ahigh_NP2 cells with high expression of HIF1A were predominant, which is consistent with previous research indicating that HIF1A is highly expressed in NP cells of normal intervertebral discs to adapt to hypoxia. In the NO and MY groups, the main composition consisted of HIF1Alow_NP1 cells with low expression of HIF1A, suggesting that the hypoxic environment within the intervertebral discs may be disrupted in these groups, leading to the suppression of HIF1A expression in the NP cells. It is worth mentioning that in the MO group, HIF1Ahigh_NP1 accounted for the majority. Although these NP cells had high expression of HIF1A, they still exhibited differences from HIF1Ahigh_NP2 cells in normal NP on the UMAP, indicating potential differential expression at other molecular levels (Figure 2C).
Figure 2. Identification of NP cell subpopulations and their distribution in different groups. (A) Expression of HIF1A in human NP cells (including NP cell, NPPC, division-NPPC) projected onto UMAP visualization. NPPC, nucleus pulposus progenitor cell. (B) UMAP plot showing the single-cell transcriptomic landscape of human NP cells. Eight subpopulations were visualized. (C) Distribution of eight NP subpopulations (defined in B) in different groups. (D) Heatmap revealing the scaled expression of DEGs for each NP subpopulation. DEGs, differentially expressed genes. (E) GO (gene ontology) functions of eight NP subpopulations defined in B. (F) Immunofluorescence staining showing the expression of HIF1A, TOMM20, SOD2 in human NP tissues. Scale bar 100 µm. (G) Quantification of HIF1A, TOMM20, SOD2 expression.
To further explore the functional differences among subgroups of cells, we calculated the DEGs of each celltype (Figure 2D) and performed GO enrichment analysis. Of note, HIF1Ahigh_NP2 showed high expression of COL2A1 and NR4A1. COL2A1 encodes type II collagen fibers, which increase intervertebral disc stability (Trefilova et al., 2021), while NR4A1 has been reported to participate in the regulatory process of inhibiting intervertebral disc degeneration (Feng et al., 2017). HIF1Alow_NP1 cells exhibited high expression of MMP3 and SOD2. MMP3 can degrade the extracellular matrix and is believed to promote intervertebral disc degeneration (Goupille et al., 1998; Roberts et al., 2000; Vo et al., 2013), while SOD2 is a pivotal antioxidant enzyme that eliminates intracellular superoxide and can convert oxygen radicals into hydrogen peroxide, thereby alleviating intracellular oxidative stress (Tamagawa et al., 2024). This suggests the presence of oxidative stress pressure within these cells and their potential contribution to intervertebral disc degeneration. HIF1Ahigh_NP1 cells showed high expression of APOD, DNAJB1, and HSPB1. APOD encodes apolipoprotein, suggesting a possible alteration in lipid metabolism pathways in these cells, while DNAJB1 and HSPB1 are known stress-response characteristic genes. Previous studies have shown that an increase in ROS levels can induce elevated expression levels of these two genes (He et al., 2024). This implies that although HIF1Ahigh_NP1 cells have low HIF1A expression, there may still be some stress factors present within these cells. GO enrichment analysis (Figure 2E) indicates that HIF1Ahigh_NP2 cells are mainly associated with ossification, connective tissue development, and cartilage development. In HIF1Alow_NP1 cells, pathways related to mitochondrial translation, aerobic respiration, and oxidative phosphorylation were significantly enriched. Previous research has shown that due to the low oxygen environment within the intervertebral disc, normal nucleus pulposus cells primarily rely on glycolysis as their main metabolic mode, and high-intensity oxidative phosphorylation is an abnormal manifestation (Agrawal et al., 2007; Francisco et al., 2023). Considering the low expression level of HIF1A in these cells, it suggests the presence of higher oxygen levels in the microenvironment that HIF1Alow_NP1 cells encounter, enabling them to engage in oxidative phosphorylation. Previous studies have indicated that high-intensity oxidative phosphorylation can increase ROS levels, leading to oxidative stress (Dimozi et al., 2015; Cheng et al., 2022; Excessive reactive oxygen species, 2022). These cells also exhibited high expression of SOD2, which may serve as a protective mechanism against high intracellular ROS levels, further reflecting the high ROS levels (Tamagawa et al., 2024).
To further explore the oxidative stress pressure of various subtypes of cells, we used the AUCell to analyze the strength of relevant pathways in different cell populations. The results showed that the hypoxia pathway was weaker in HIF1Alow_NP1, while the oxidative phosphorylation, ROS, DNA repair, and ferroptosis pathways were stronger. Interestingly, the opposite trend was observed in the HIF1Alow_NP1 cells for the above-mentioned pathways. It is worth noting that there was no significant difference in the strength of the glycolysis pathway among different cell populations (Supplementary Figure S3A). Since HIF1Alow_NP1 cells are mainly distributed in MY, with fewer in MO, indicating an upregulation of oxidative stress levels in the NP of young degenerative patients, which is a defining feature of this type of disc degeneration. We performed hematoxylin-eosin staining (H&E stain, Supplementary Figure S3B) and immunofluorescence staining on NP tissue samples from MY and MO patients. The results showed lower HIF1A level and higher TOMM20 level in the MY group, with a patchy distribution suggestive of intense oxidative phosphorylation. SOD2 level increased in the MY group (Figures 2F, G), while COL2A1 level decreased (Supplementary Figures S3C, D). In summary, the above results indicate that NP cells from patients in the NO and MY groups exhibit relatively higher levels of oxidative phosphorylation intensity, resulting in a higher degree of oxidative stress within their cells. Although NP cells from patients in the MO group also express stress response characteristic genes, their pathogenic factors are unrelated to oxidative phosphorylation and instead associated with the unfolded proteins and other stressors, indicating endoplasmic reticulum stress.
Pseudotime trajectories analysis reveals distinct differentiation trajectories of nucleus pulposus cells
To investigate how different terminally differentiated NP cells differentiate from NPPC, we conducted pseudotime trajectories analysis using Slingshot. Since NPPC1 possesses strong proliferative ability (Figure 2E; Supplementary Figure S4 cell cycle), it is recognized as the most stem-like NP cell and set as the starting point of differentiation. We found three differentiation trajectories in NP cell differentiation (Figures 3A–D), from NPPC1 to NPPC2, further to NPPC3 and HIF1Alow_NP3, and subsequently branching into three distinct outcomes: HIF1Alow_NP1, HIF1Ahigh_NP1, and HIF1Ahigh_NP2. These three types of cells represent the dominant cell components of the NO and MY, MO, and NY groups, respectively. pseudotime trajectories analysis suggests that they also represent terminal differentiation states within their respective groups. Therefore, it can be inferred that NPPC2, NPPC3, and HIF1Alow_NP3 serve as common transitional states during NPPC1 differentiation into terminally differentiated NP cells. In trajectory 1, HIF1Alow_NP2 represents a unique transitional state in this fate process. Thus, through pseudotime trajectories analysis, we have identified the differentiation trajectories of representative cells in each group, and their fate determination may depend on the environmental factors encountered by HIF1Alow_NP3.
Figure 3. Trajectories of NP subpopulations. (A)Slingshot-based pseudotime trajectories inferred from UMAP embedding of NP subpopulations. (B) showed the trajectory, NPPC1, NPPC2, NPPC3, HIF1Alow_NP3, HIF1Alow_NP1. (C) showed the trajectory, NPPC1, NPPC2, NPPC3, HIF1Alow_NP3, HIF1Ahigh_NP1. (D) showed the trajectory, NPPC1, NPPC2, NPPC3, HIF1Alow_NP3, HIF1Ahigh_NP2.
Dendritic cells play an important role in the IVD
Due to the avascular nature of intervertebral discs, they are in an immune-privileged state. To explore the pathological and physiological mechanisms of immune cells in different groups, we performed subpopulation analysis of myeloid cells. Based on the classic marker and GO enrichment analysis results, we divided myeloid cells into six celltypes (Figures 4A, B; Supplementary Figure S5), namely, monocytes 1-2 (CD14, FCGR3A, FCGR3B), dendritic cells 1-2 (ITGAX, CD1C, CD83), and macrophages 1-2 (CD68, CSF1R, CD163). We found that myeloid cells were mainly distributed in degenerated groups, and the MY group had more myeloid cells. Among the myeloid cells, dendritic cells were the predominant component (Figure 4C). In monocytes, monocyte 1 exhibited strong proliferative activity (MKI67+), while the GO enrichment analysis results for monocyte 2 (Figure 4D) suggested significant enrichment of antigen processing and presentation function, which is an important function of dendritic cells. Although both dendritic cell 1 and 2 showed high expression of classic markers, only dendritic cell 2 had significant enrichment of antigen processing and presentation function, revealing the heterogeneity within dendritic cells. Additionally, we also found significant enrichment of the response to lipopolysaccharide pathway in dendritic cells 1 and 2, suggesting that they exhibit a pro-inflammatory phenotype. The AUCell results indicated a significant increase in the intensity of the inflammation pathway in dendritic cells 1 and 2, while the intensity was relatively low in other myeloid cell populations (Figure 4E).
Figure 4. Identification of myeloid cell subpopulations in human IVD and their distribution in different groups. (A) UMAP plot showing the single-cell transcriptomic landscape of myeloid cells. Six clusters were visualized. (B) Dot plot of marker genes for myeloid cell subpopulations. Color scale indicates the mean of normalized expression of marker genes in each cell type, and dot size is proportional to the percentage of cells within each cell cluster expressing the marker genes. (C) Same UMAP visualization as the one in B but split by different groups. (D) GO (gene ontology) functions of six cluster defined in (A). (E) UMAP visualization AUC scores of the inflammation pathway in myeloid cells.
Through Slingshot pseudotime analysis, we discovered two differentiation trajectories from monocytes to dendritic cells and monocytes to macrophages (Figure 5A). Monocyte 1 exhibited strong proliferative capacity, so we defined it as the starting point of differentiation. Slingshot pseudotime analysis showed two differentiation trajectories. Monocyte1 differentiated into macrophages 1-2 or monocyte1 differentiated into dendritic cells 1-2 via monocyte 2, which served as a transitional state (Figures 5B, C). This suggests that once monocytes begin to differentiate, their proliferative capacity diminishes, and they start exhibiting immune-related functions, particularly enhanced antigen presentation. This indicates a tendency for monocyte 2 to differentiate into dendritic cells, as supported by the higher proportion of dendritic cells in this group, confirming the accuracy of this differentiation trajectory. Dendritic cell 1 represents the transitional state from monocyte 2 to dendritic cell 2. Interestingly, during this transitional state, antigen presentation functionality decreases but reappears in the terminally differentiated dendritic cell 2 (Figure 4D). The role of dendritic cells is to present antigens through MHC-II molecules to facilitate antigen recognition by CD4+ T cells and initiate downstream immune responses. In the immune-privileged state of the intervertebral disc, immune cell infiltration may lead to the recognition of self-antigens present in the NP as foreign substances. Consequently, there is a tendency for a significant differentiation into dendritic cells to meet the demand for antigen presentation. We found more CD68 positive cell in MY group with immunofluorescence staining, which is consist with our single-cell sequencing analysis results (Figures 5D, E). Therefore, in the NP of MY group, dendritic cells are involved in antigen presentation of internal antigens and exhibit pro-inflammatory responses.
Figure 5. Trajectories of myeloid cell subpopulations. (A) Slingshot-based pseudotime trajectories inferred from UMAP embedding of myeloid cell subpopulations. (B) showed the trajectory, monocyte1, macrophage1, macrophage2. (C) showed the trajectory, monocyte1, dendritic cell1, dendritic cell2. (D) Immunofluorescence staining showing the expression of CD68 in human NP tissues. Scale bar 100 µm. (E) Quantification of CD68 expression.
Cellular communication analysis reveals the interactions between various cell types in different groups
To further understand the differences in cellular communication between the groups, we used CellChat to analyze cellular communication among all cell types in the four groups (Figure 6A). The results showed that within the HIF1Ahigh_NP2 subgroup of the NY group, there is a strong mutual communication between NPPC3 and HIF1Ahigh_NP2. In the NO group, self-communication exists only within the HIF1Alow_NP1 subgroup. In the MO group, communication mainly occurs between HIF1Ahigh_NP1, HIF1Alow_NP3, and NPPC2. In the MY group, there is communication of varying intensities between HIF1Alow_NP1, NPPC2, and dendritic cells 1-2.
Figure 6. Overview of the crosstalk networks among the clusters of human IVD cells. (A) The intercellular crosstalk networks among NP cells and myeloid cells in different groups. (B–E) Dot plot showing the communication probability of the indicated ligand-receptor pairs between NP cell subpopulations and myeloid cells subpopulations in different groups.
Bubble plots display communication within the NY group (Figure 6B). HIF1Ahigh_NP2 interacts with NPPC3 through FGF2/FGFR2. FGF2/FGFR2 plays a crucial role in various important biological processes such as cell survival, proliferation, differentiation, and migration. We speculate that this may be a mechanism maintaining normal homeostasis in the NY group. Additionally, GAS6/AXL communication also exists within HIF1Ahigh_NP2, which is a classic efferocytosis mechanism widely found in the process of macrophage engulfing apoptotic cells (Bellan et al., 2019; Vago et al., 2021). Previous studies have indicated that NP cells possess certain phagocytic ability (Nerlich et al., 2002; Jones et al., 2008; Lin et al., 2019). Hence, we hypothesize that the senescent NP cells within the NY group, lacking immune cells, are likely engulfed by the healthy NP cells through the process of phagocytosis, thereby maintaining internal homeostasis. It is noteworthy that POSTN/(ITGAV + ITGB5) and TNFSF11/TNFRSF11B represent the most significant modes of communication in the NY group, primarily occurring between NPPC3 and HIF1Ahigh_NP2 cells. In the NO group (Figure 6C), ANGPTL2/(ITGA5 + ITGB1), FGF2/FGFR1, and NAMPT/(ITGA5 + ITGB1) communication is predominantly observed within HIF1Alow_NP1 cells. In the MO group (Figure 6D), the main mode of communication for HIF1Alow_NP1-3 is through FGF2/FGFR1 with HIF1Ahigh_NP1 cells. In the MY group (Figure 6E), the most significant modes of communication are FGF2/FGFR1 and MIF/ACKR3 within HIF1Alow_NP1. Ligand MDK and receptors LRP1, NCL, SDC2, and SDC4 mediate communication between NPPC2 and HIF1Alow_NP1. Furthermore, dendritic cells 1-2 secrete SPP1 protein that acts on multiple receptors of HIF1Alow_NP1 cells. HIF1Alow_NP1 cells can also affect dendritic cells 1-2 through MIF acting on CD74, CXCR4, and CD44 receptors (Figure 3). This suggests a strong immunomodulatory effect between dendritic cells and HIF1Alow_NP1 within the MY group. Additionally, HIF1Alow_NP1 can act on endothelial cell VEGFR through VEGFA, potentially leading to vascular formation within the intervertebral disc and disrupting its hypoxic environment. In summary, we have discovered distinct characteristics of cellular communication among different groups. These communication modes suggest potential molecular mechanisms for maintaining homeostasis in NP cells and the changes that occur in communication between cell types when homeostasis is disrupted under different circumstances.
Discussion
Intervertebral disc degeneration is often thought to occur mainly in the elderly population. However, there is a gradual increase in the proportion of young patients (Liu et al., 2023). In recent years, single-cell sequencing has been increasingly used in the field of intervertebral disc degeneration (Wang et al., 2023; Jiang et al., 2022; Guo et al., 2023b; Gao et al., 2022; Gan et al., 2021). However, most studies have focused on exploring the severity of degeneration or disc development, without investigating the differences between intervertebral disc degeneration in different age groups. In order to determine whether there are differences in degenerative mechanisms between young and elderly patients, we collected single-cell sequencing data from intervertebral discs of 10 cases, including Normal Young (NY), Normal Old (NO), Mild Young (MY), and Mild Old (MO) types based on the Pfirrmann classification. We conducted a single-cell resolution analysis to elucidate the differences between young and elderly degenerative patients.
During the aging process, the proportion of nucleus pulposus progenitor cells (NPPCs) decreases, indicating a depletion of stem cells. We subdivided NP cells based on the expression levels of hypoxia-inducible factor 1 alpha (HIF1A) and identified a subset in young patients with low HIF1A expression and aerobic respiration, which we termed HIF1Alow_NP1. It is worth noting that this cell type is predominant in the intervertebral discs of normal elderly individuals. Previous studies have shown that NP cells primarily metabolize through glycolysis to adapt to the hypoxic environment in the disc center. It has been indicated that under hypoxic conditions, HIF1A expression in NP cells increases, leading to the inhibition of reactive oxygen species (ROS) generation. Conversely, normoxic conditions induce downregulation of HIF1A expression in NP cells, resulting in increased ROS levels. Furthermore, under hypoxic conditions, the decrease in ROS in NP cells is reliant on HIF1A, emphasising the crucial function of HIF1A in the typical survival of NP cells under hypoxic conditions (Peng et al., 2024).
HIF1Alow_NP1 cells are the primary component in both young degenerative patients and normal elderly disc nucleus pulposus cells. This suggests that the anaerobic environment within their intervertebral discs may be disrupted. Previous studies have shown that oxidative phosphorylation can induce ROS generation (Zhang et al., 2022). Elevated levels of oxidative phosphorylation within HIF1Alow_NP1 cells may lead to increased intracellular ROS generation, inducing oxidative stress damage. High expression of SOD2 and MMP3 was observed in this cell type. SOD2 is an important antioxidant enzyme that clears intracellular superoxide (Bolduc et al., 2019), and its upregulation indirectly indicates elevated intracellular oxidative stress levels. MMP3 can degrade extracellular matrix components, including various types of collagen, and is believed to promote disc degeneration (Lerner et al., 2018). Immunofluorescence staining confirmed a decrease in type II collagen in young patients, which may be associated with the high expression of MMP3 in HIF1Alow_NP1 cells. This reduction in extracellular matrix weakens the disc’s resistance to mechanical pressure, exacerbating cell death within the disc and promoting degeneration (Hu et al., 2020).
Intercellular communication indicates that HIF1Alow_NP1 cells can influence endothelial cells through VEGFA, promoting angiogenesis within the disc and disrupting the hypoxic environment. Therefore, the HIF1Alow_NP1 cell population is a significant characteristic of intervertebral disc degeneration in young patients, and these cells may contribute to disc degeneration. In contrast, elderly patients exhibit a predominance of HIF1Ahigh_NP1 cells within the intervertebral disc. Although these cells share high HIF1A expression with normal NP cells, functional enrichment analysis suggests that they respond to unfolded proteins and other stressors, indicating endoplasmic reticulum stress (Di Conza et al., 2023). Previous studies have also shown that endoplasmic reticulum stress is a significant pathological feature of disc degeneration, leading to apoptosis of NP cells (Liao et al., 2019).
To summarise, intervertebral disc degeneration in young patients is characterised by oxidative stress, whereas in elderly patients, endoplasmic reticulum stress is a significant feature.
The study found that immune cells were only present in degenerative samples, while they were absent in the intervertebral discs of normal and elderly populations. This is consistent with previous reports of the disc being an immune-privileged organ (Sun et al., 2020; Ye et al., 2022). The sequencing data and immunofluorescence staining results indicated a higher proportion of immune cells in young patients compared to elderly patients, which suggests a potential role for immune cells in young patients. Previous studies have indicated that autoimmunity may play an important role in IVDD. In these circumstances, the autoimmune disorder commences following the exposure of antigenic disc components to the circulating blood, resulting in abnormalities in cell-mediated and humoral autoimmune responses and the production of chronic inflammation with low back pain and sciatica (Di Martino et al., 2013). Monocytes within the myeloid cell subset are more likely to differentiate into dendritic cells than macrophages, although both can represent antigen-presenting cells in the process of autoimmune response. Dendritic cells present antigen to naïve CD4+ T cells during the initiation of an immune response (Geiss et al., 2016), whereas macrophages display antigen to previously activated CD4+ T cells at a later stage of the manifested immune response. Additionally, dendritic cells display a stronger pro-inflammatory phenotype (Lan et al., 2022). Studies have shown that a culture medium that favours M1 macrophages can promote intervertebral disc degeneration in rats. This can induce NP cells to express MMP13 and inflammatory factors, suggesting that pro-inflammatory immune cells may have adverse effects on NP cells (Li et al., 2022). Therefore, the pro-inflammatory dendritic cells in this study may further promote disc degeneration. As previously stated, the infiltration of immune cells in the NP area is closely related to the formation of nerves and blood vessels, which may result in an interaction between the NP and the blood circulation. Thus, a recent study has demonstrated that the circulating levels of cartilage oligomeric matrix protein (COMP) and a disintegrin and metalloproteinase with thrombospondin motifs 7 (ADAMTS7) can be examined as promising indicators to detect IVDD (Ding et al., 2024). These biomarkers may help us better dignose IVDD at an early stage.
Future research should aim to elucidate the following points. Why do young and elderly patients present different pathological characteristics in intervertebral disc degeneration? The pathogenic mechanisms in young patients are unclear. Can inhibiting oxidative stress or endoplasmic reticulum stress improve disc degeneration, and what are the key mechanisms and targets? Finally, this project leaves unanswered questions about why monocytes differentiate into dendritic cells and the specific roles that dendritic cells play in intervertebral disc degeneration. These questions require further exploration in future studies.
Data availability statement
The data presented in the study are deposited in the GEO database, accession number (GSE233666, GSE205535, GSE189916), and in the CNGBdb database accession number CNP0002664 available via the following link: https://db.cngb.org/search/project/CNP0002664/.
Ethics statement
The studies involving humans were approved by The Institutional Review Board of the Shanghai General Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from primarily isolated as part of your previous study for which ethical approval was obtained. Written informed consent for participation was not required from the participants or the participant’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
PL: Methodology, Formal Analysis, Data curation, Writing–review and editing, Writing–original draft, Conceptualization. XR: Software, Investigation, Writing–review and editing, Writing–original draft, Formal Analysis, Data curation, Conceptualization. BZ: Methodology, Writing–review and editing, Writing–original draft, Software, Data curation. SG: Validation, Resources, Writing–review and editing, Methodology. QF: Writing–original draft, Supervision, Funding acquisition, Conceptualization, Writing–review and editing, Validation.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported in part by grants from the National Natural Science Foundation of China (NSFC) (No. 82202694 and 82272539).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2024.1409287/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | UMAP visualization showing the human IVD cells were divided into 10 clusters.
SUPPLEMENTARY FIGURE S2 | Violin polt showing HIF1A expression in different NP cell subpopulations. Using wilcoxon test, p < 0.0001 (****).
SUPPLEMENTARY FIGURE S3 | (A) Strength of the glycolysis pathway among different cell populations. (B) H&E staining on NP tissues. (C) Immunofluorescence staining showing the expression of COL2A1 in human NP tissues. Scale bar 100 µm. (D) Quantification of COL2A1 expression.
SUPPLEMENTARY FIGURE S4 | (A) Bar graph representing the proportion of each cell cycle stage in 8 NP cell subpopulations. (B) Bar graph representing the amount of cells in 8 NP cell subpopulations.
SUPPLEMENTARY FIGURE S5 | (A) UMAP visualization showing the immune cells were divided into 6 clusters.
SUPPLEMENTARY FIGURE S6 | (A) The sagittal and horizontal views of the computed tomography image serve to illustrate the area that was chosen for the immunofluorescence procedure. (B) A schematic diagram of the disc represent the area that was chosen for the immunofluorescence.
References
Agrawal, A., Guttapalli, A., Narayan, S., Albert, T. J., Shapiro, I. M., and Risbud, M. V. (2007). Normoxic stabilization of HIF-1alpha drives glycolytic metabolism and regulates aggrecan gene expression in nucleus pulposus cells of the rat intervertebral disk. Am. J. Physiol. Cell. Physiol. 293 (2), C621–C631. doi:10.1152/ajpcell.00538.2006
Bellan, M., Cittone, M. G., Tonello, S., Rigamonti, C., Castello, L. M., Gavelli, F., et al. (2019). Gas6/TAM system: a key modulator of the interplay between inflammation and fibrosis. Int. J. Mol. Sci. 20 (20), 5070. doi:10.3390/ijms20205070
Bolduc, J. A., Collins, J. A., and Loeser, R. F. (2019). Reactive oxygen species, aging and articular cartilage homeostasis. Free Radic. Biol. Med. 132, 73–82. doi:10.1016/j.freeradbiomed.2018.08.038
Chen, F., Lei, L., Chen, S., Zhao, Z., Huang, Y., Jiang, G., et al. (2024). Serglycin secreted by late-stage nucleus pulposus cells is a biomarker of intervertebral disc degeneration. Nat. Commun. 15 (1), 47. doi:10.1038/s41467-023-44313-9
Cheng, F., Yang, H., Cheng, Y., Liu, Y., Hai, Y., and Zhang, Y. (2022). The role of oxidative stress in intervertebral disc cellular senescence. Front. Endocrinol. (Lausanne) 13, 1038171. doi:10.3389/fendo.2022.1038171
Cherif, H., Mannarino, M., Pacis, A. S., Ragoussis, J., Rabau, O., Ouellet, J. A., et al. (2022). Single-cell RNA-seq analysis of cells from degenerating and non-degenerating intervertebral discs from the Same individual reveals new biomarkers for intervertebral disc degeneration. Int. J. Mol. Sci. 23 (7), 3993. doi:10.3390/ijms23073993
Deng, Y., Zheng, Y., Li, D., Hong, Q., Zhang, M., Li, Q., et al. (2021). Expression characteristics of interferon-stimulated genes and possible regulatory mechanisms in lupus patients using transcriptomics analyses. EBioMedicine 70, 103477. doi:10.1016/j.ebiom.2021.103477
Di Conza, G., Ho, P. C., Cubillos-Ruiz, J. R., and Huang, S. C. C. (2023). Control of immune cell function by the unfolded protein response. Nat. Rev. Immunol. 23 (9), 546–562. doi:10.1038/s41577-023-00838-0
Di Martino, A., Merlini, L., and Faldini, C. (2013). Autoimmunity in intervertebral disc herniation: from bench to bedside. Expert Opin. Ther. Targets 17 (12), 1461–1470. doi:10.1517/14728222.2013.834330
Dimozi, A., Mavrogonatou, E., Sklirou, A., and Kletsas, D. (2015). Oxidative stress inhibits the proliferation, induces premature senescence and promotes a catabolic phenotype in human nucleus pulposus intervertebral disc cells. Eur. Cell. Mater 30, 89–102. doi:10.22203/ecm.v030a07
Ding, J. Y., Yan, X., Zhang, R. J., Zhang, H. Q., Kang, L., Jia, C. Y., et al. (2024). Diagnostic value of serum COMP and ADAMTS7 for intervertebral disc degeneration. Eur. J. Med. Res. 29 (1), 196. doi:10.1186/s40001-024-01784-w
Excessive reactive oxygen species (2022). Excessive reactive oxygen species are therapeutic targets for intervertebral disc degeneration - PubMed. Available at: https://pubmed.ncbi.nlm.nih.gov/26542776/.
Feng, G., Zhang, Z., Dang, M., Zhang, X., Doleyres, Y., Song, Y., et al. (2017). Injectable nanofibrous spongy microspheres for NR4A1 plasmid DNA transfection to reverse fibrotic degeneration and support disc regeneration. Biomaterials 131, 86–97. doi:10.1016/j.biomaterials.2017.03.029
Fernandes, L. M., Khan, N. M., Trochez, C. M., Duan, M., Diaz-Hernandez, M. E., Presciutti, S. M., et al. (2020). Single-cell RNA-seq identifies unique transcriptional landscapes of human nucleus pulposus and annulus fibrosus cells. Sci. Rep. 10 (1), 15263. doi:10.1038/s41598-020-72261-7
Francisco, V., Ait Eldjoudi, D., González-Rodríguez, M., Ruiz-Fernández, C., Cordero-Barreal, A., Marques, P., et al. (2023). Metabolomic signature and molecular profile of normal and degenerated human intervertebral disc cells. Spine J. 23 (10), 1549–1562. doi:10.1016/j.spinee.2023.06.005
Gan, Y., He, J., Zhu, J., Xu, Z., Wang, Z., Yan, J., et al. (2021). Spatially defined single-cell transcriptional profiling characterizes diverse chondrocyte subtypes and nucleus pulposus progenitors in human intervertebral discs. Bone Res. 9 (1), 37. doi:10.1038/s41413-021-00163-z
Gantenbein, B., Sun, Z., Liu, Z., and Samartzis, D. (2023). Editorial: immunological imbalance: what is its role in intervertebral disc degeneration? Front. Cell. Dev. Biol. 11, 1196377. doi:10.3389/fcell.2023.1196377
Gao, B., Jiang, B., Xing, W., Xie, Z., Luo, Z., and Zou, W. (2022). Discovery and application of postnatal nucleus pulposus progenitors essential for intervertebral disc homeostasis and degeneration. Adv. Sci. (Weinh) 9 (13), e2104888. doi:10.1002/advs.202104888
Geiss, A., Sobottke, R., Delank, K. S., and Eysel, P. (2016). Plasmacytoid dendritic cells and memory T-cells infiltrate true sequestrations stronger than subligamentous sequestrations: evidence from flow cytometric analysis of disc infiltrates. Eur. Spine J. 25 (5), 1417–1427. doi:10.1007/s00586-015-4325-z
Goupille, P., Jayson, M. I., Valat, J. P., and Freemont, A. J. (1998). Matrix metalloproteinases: the clue to intervertebral disc degeneration? Spine (Phila Pa 1976) 23 (14), 1612–1626. doi:10.1097/00007632-199807150-00021
Guo, S., Yan, M., Li, X., Zhang, S., Liu, Z., Li, K., et al. (2023a). Single-cell RNA-seq analysis reveals that immune cells induce human nucleus pulposus ossification and degeneration. Front. Immunol. 14, 1224627. doi:10.3389/fimmu.2023.1224627
Guo, S., Yan, M., Li, X., Zhang, S., Liu, Z., Li, K., et al. (2023b). Single-cell RNA-seq analysis reveals that immune cells induce human nucleus pulposus ossification and degeneration. Front. Immunol. 14, 1224627. doi:10.3389/fimmu.2023.1224627
Han, S., Zhang, Y., Zhang, X., Zhang, H., Meng, S., Kong, M., et al. (2022). Single-cell RNA sequencing of the nucleus pulposus reveals chondrocyte differentiation and regulation in intervertebral disc degeneration. Front. Cell. Dev. Biol. 10, 824771. doi:10.3389/fcell.2022.824771
He, J., Hou, X., Wu, J., Wang, K., Qi, X., Wei, Z., et al. (2024). Hspb1 protects against severe acute pancreatitis by attenuating apoptosis and ferroptosis via interacting with Anxa2 to restore the antioxidative activity of Prdx1. Int. J. Biol. Sci. 20 (5), 1707–1728. doi:10.7150/ijbs.84494
Hu, B., Zhang, S., Liu, W., Wang, P., Chen, S., Lv, X., et al. (2020). Inhibiting heat shock protein 90 protects nucleus pulposus-derived stem/progenitor cells from compression-induced necroptosis and apoptosis. Front. Cell. Dev. Biol. 8, 685. doi:10.3389/fcell.2020.00685
Iwakiri, K., Sasaki, T., Sotoyama, M., Du, T., Miki, K., and Oyama, F. (2023). Effect of relative weight limit set as a body weight percentage on work-related low back pain among workers. PLoS One 18 (4), e0284465. doi:10.1371/journal.pone.0284465
Jiang, W., Glaeser, J. D., Salehi, K., Kaneda, G., Mathkar, P., Wagner, A., et al. (2022). Single-cell atlas unveils cellular heterogeneity and novel markers in human neonatal and adult intervertebral discs. iScience 25 (7), 104504. doi:10.1016/j.isci.2022.104504
Jin, S., Guerrero-Juarez, C. F., Zhang, L., Chang, I., Ramos, R., Kuan, C. H., et al. (2021). Inference and analysis of cell-cell communication using CellChat. Nat. Commun. 12 (1), 1088. doi:10.1038/s41467-021-21246-9
Jones, P., Gardner, L., Menage, J., Williams, G. T., and Roberts, S. (2008). Intervertebral disc cells as competent phagocytes in vitro: implications for cell death in disc degeneration. Arthritis Res. Ther. 10 (4), R86. doi:10.1186/ar2466
Knezevic, N. N., Candido, K. D., Vlaeyen, J. W. S., Van Zundert, J., and Cohen, S. P. (2021). Low back pain. Lancet. 398 (10294), 78–92. doi:10.1016/S0140-6736(21)00733-9
Knezevic, N. N., Mandalia, S., Raasch, J., Knezevic, I., and Candido, K. D. (2017). Treatment of chronic low back pain - new approaches on the horizon. J. Pain Res. 10, 1111–1123. doi:10.2147/JPR.S132769
Korsunsky, I., Millard, N., Fan, J., Slowikowski, K., Zhang, F., Wei, K., et al. (2019). Fast, sensitive and accurate integration of single-cell data with Harmony. Nat. Methods 16 (12), 1289–1296. doi:10.1038/s41592-019-0619-0
Lama, P., Tewari, J., Adams, M. A., and Le Maitre, C. (2022). Degenerative physiochemical events in the pathological intervertebral disc. Histol. Histopathol. 37 (1), 11–20. doi:10.14670/HH-18-395
Lan, T., Hu, Z., Guo, W., Yan, B., and Zhang, Y. (2022). Development of a novel inflammatory-associated gene signature and immune infiltration patterns in intervertebral disc degeneration. Oxid. Med. Cell. Longev. 2022, 2481071. doi:10.1155/2022/2481071
Lerner, A., Neidhöfer, S., Reuter, S., and Matthias, T. (2018). MMP3 is a reliable marker for disease activity, radiological monitoring, disease outcome predictability, and therapeutic response in rheumatoid arthritis. Best. Pract. Res. Clin. Rheumatol. 32 (4), 550–562. doi:10.1016/j.berh.2019.01.006
Li, X. C., Luo, S. J., Fan, W., Zhou, T. L., Tan, D. Q., Tan, R. X., et al. (2022). Macrophage polarization regulates intervertebral disc degeneration by modulating cell proliferation, inflammation mediator secretion, and extracellular matrix metabolism. Front. Immunol. 13, 922173. doi:10.3389/fimmu.2022.922173
Liao, Z., Luo, R., Li, G., Song, Y., Zhan, S., Zhao, K., et al. (2019). Exosomes from mesenchymal stem cells modulate endoplasmic reticulum stress to protect against nucleus pulposus cell death and ameliorate intervertebral disc degeneration in vivo. Theranostics 9 (14), 4084–4100. doi:10.7150/thno.33638
Lin, Y., Cong, H., Liu, K., Jiao, Y., Yuan, Y., Tang, G., et al. (2019). Microbicidal phagocytosis of nucleus pulposus cells against Staphylococcus aureus via the TLR2/MAPKs signaling pathway. Front. Immunol. 10, 1132. doi:10.3389/fimmu.2019.01132
Ling, Z., Liu, Y., Wang, Z., Zhang, Z., Chen, B., Yang, J., et al. (2021). Single-cell RNA-seq analysis reveals macrophage involved in the progression of human intervertebral disc degeneration. Front. Cell. Dev. Biol. 9, 833420. doi:10.3389/fcell.2021.833420
Liu, X., Jin, L., Jiang, C., Jiang, X., Chen, Z., and Cao, Y. (2023). Characteristics of cervical intervertebral disc signal intensity: an analysis of T2-weighted magnetic resonance imaging in 5843 asymptomatic Chinese subjects. Eur. Spine J. 32 (7), 2415–2424. doi:10.1007/s00586-023-07742-0
Madhu, V., Hernandez-Meadows, M., Boneski, P. K., Qiu, Y., Guntur, A. R., Kurland, I. J., et al. (2023). The mitophagy receptor BNIP3 is critical for the regulation of metabolic homeostasis and mitochondrial function in the nucleus pulposus cells of the intervertebral disc. Autophagy 19 (6), 1821–1843. doi:10.1080/15548627.2022.2162245
Maher, C., Underwood, M., and Buchbinder, R. (2017). Non-specific low back pain. Lancet. 389 (10070), 736–747. doi:10.1016/S0140-6736(16)30970-9
McGinnis, C. S., Murrow, L. M., and Gartner, Z. J. (2019). DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell. Syst. 8 (4), 329–337. doi:10.1016/j.cels.2019.03.003
Meisel, H. J., Agarwal, N., Hsieh, P. C., Skelly, A., Park, J. B., Brodke, D., et al. (2019). Cell therapy for treatment of intervertebral disc degeneration: a systematic review. Glob. Spine J. 9 (1 Suppl. l), 39S-52S–52S. doi:10.1177/2192568219829024
Nerlich, A. G., Weiler, C., Zipperer, J., Narozny, M., and Boos, N. (2002). Immunolocalization of phagocytic cells in normal and degenerated intervertebral discs. Spine (Phila Pa 1976) 27 (22), 2484–2490. doi:10.1097/00007632-200211150-00012
Ngo, K., Pohl, P., Wang, D., Leme, A. S., Lee, J., Di, P., et al. (2017). ADAMTS5 deficiency protects mice from chronic tobacco smoking-induced intervertebral disc degeneration. Spine (Phila Pa 1976) 42 (20), 1521–1528. doi:10.1097/BRS.0000000000002258
Novais, E. J., Narayanan, R., Canseco, J. A., van de Wetering, K., Kepler, C. K., Hilibrand, A. S., et al. (2024). A new perspective on intervertebral disc calcification-from bench to bedside. Bone Res. 12 (1), 3. doi:10.1038/s41413-023-00307-3
Panebianco, C. J., Dave, A., Charytonowicz, D., Sebra, R., and Iatridis, J. C. (2021). Single-cell RNA-sequencing atlas of bovine caudal intervertebral discs: discovery of heterogeneous cell populations with distinct roles in homeostasis. FASEB J. 35 (11), e21919. doi:10.1096/fj.202101149R
Peng, Z., Wang, Z., Guo, S., Tan, B., Cao, R., Cheng, S., et al. (2024). An analysis of the transcriptional landscape in hypoxia-treated primary nucleus pulposus cells. Genes. Dis. 11 (2), 558–560. doi:10.1016/j.gendis.2023.04.009
Pfirrmann, C. W., Metzdorf, A., Zanetti, M., Hodler, J., and Boos, N. (2001). Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976) 26 (17), 1873–1878. doi:10.1097/00007632-200109010-00011
RNA (2022). Single-cell RNA sequencing reveals the difference in human normal and degenerative nucleus pulposus tissue profiles and cellular interactions - PubMed. Available at: https://pubmed.ncbi.nlm.nih.gov/35874809/.
Roberts, S., Caterson, B., Menage, J., Evans, E. H., Jaffray, D. C., and Eisenstein, S. M. (2000). Matrix metalloproteinases and aggrecanase: their role in disorders of the human intervertebral disc. Spine (Phila Pa 1976) 25 (23), 3005–3013. doi:10.1097/00007632-200012010-00007
Silagi, E. S., Schipani, E., Shapiro, I. M., and Risbud, M. V. (2021). The role of HIF proteins in maintaining the metabolic health of the intervertebral disc. Nat. Rev. Rheumatol. 17 (7), 426–439. doi:10.1038/s41584-021-00621-2
Song, C., Hu, P., Peng, R., Li, F., Fang, Z., and Xu, Y. (2024). Bioenergetic dysfunction in the pathogenesis of intervertebral disc degeneration. Pharmacol. Res. 202, 107119. doi:10.1016/j.phrs.2024.107119
Soufi, K. H., Castillo, J. A., Rogdriguez, F. Y., DeMesa, C. J., and Ebinu, J. O. (2023). Potential role for stem cell regenerative therapy as a treatment for degenerative disc disease and low back pain: a systematic review. Int. J. Mol. Sci. 24 (10), 8893. doi:10.3390/ijms24108893
Street, K., Risso, D., Fletcher, R. B., Das, D., Ngai, J., Yosef, N., et al. (2018). Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics. BMC Genomics 19 (1), 477. doi:10.1186/s12864-018-4772-0
Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W. M., et al. (2019). Comprehensive integration of single-cell data. Cell. 177 (7), 1888–1902. doi:10.1016/j.cell.2019.05.031
Sun, Z., Liu, B., and Luo, Z. J. (2020). The immune privilege of the intervertebral disc: implications for intervertebral disc degeneration treatment. Int. J. Med. Sci. 17 (5), 685–692. doi:10.7150/ijms.42238
Tamagawa, S., Sakai, D., Nojiri, H., Nakamura, Y., Warita, T., Matsushita, E., et al. (2024). SOD2 orchestrates redox homeostasis in intervertebral discs: a novel insight into oxidative stress-mediated degeneration and therapeutic potential. Redox Biol. 71, 103091. doi:10.1016/j.redox.2024.103091
Trefilova, V. V., Shnayder, N. A., Petrova, M. M., Kaskaeva, D. S., Tutynina, O. V., Petrov, K. V., et al. (2021). The role of polymorphisms in collagen-encoding genes in intervertebral disc degeneration. Biomolecules 11 (9), 1279. doi:10.3390/biom11091279
Vago, J. P., Amaral, F. A., and van de Loo, F. A. J. (2021). Resolving inflammation by TAM receptor activation. Pharmacol. Ther. 227, 107893. doi:10.1016/j.pharmthera.2021.107893
Vo, N. V., Hartman, R. A., Yurube, T., Jacobs, L. J., Sowa, G. A., and Kang, J. D. (2013). Expression and regulation of metalloproteinases and their inhibitors in intervertebral disc aging and degeneration. Spine J. 13 (3), 331–341. doi:10.1016/j.spinee.2012.02.027
Wang, D., Li, Z., Huang, W., Cao, S., Xie, L., Chen, Y., et al. (2023). Single-cell transcriptomics reveals heterogeneity and intercellular crosstalk in human intervertebral disc degeneration. iScience 26 (5), 106692. doi:10.1016/j.isci.2023.106692
Yang, S., Jing, S., Wang, S., and Jia, F. (2024). From drugs to biomaterials: a review of emerging therapeutic strategies for intervertebral disc inflammation. Front. Cell. Infect. Microbiol. 14, 1303645. doi:10.3389/fcimb.2024.1303645
Yang, W., Jia, C., Liu, L., Fu, Y., Wu, Y., Liu, Z., et al. (2023). Hypoxia-inducible factor-1α protects against intervertebral disc degeneration through antagonizing mitochondrial oxidative stress. Inflammation 46 (1), 270–284. doi:10.1007/s10753-022-01732-y
Ye, F., Lyu, F. J., Wang, H., and Zheng, Z. (2022). The involvement of immune system in intervertebral disc herniation and degeneration. JOR Spine 5 (1), e1196. doi:10.1002/jsp2.1196
Zhang, B., Pan, C., Feng, C., Yan, C., Yu, Y., Chen, Z., et al. (2022). Role of mitochondrial reactive oxygen species in homeostasis regulation. Redox Rep. 27 (1), 45–52. doi:10.1080/13510002.2022.2046423
Keywords: intervertebral disc degeneration, single-cell genomics, nucleus pulposus, hypoxia-inducible factor 1 alpha, dendritic cells
Citation: Liu P, Ren X, Zhang B, Guo S and Fu Q (2024) Investigating the characteristics of mild intervertebral disc degeneration at various age stages using single-cell genomics. Front. Cell Dev. Biol. 12:1409287. doi: 10.3389/fcell.2024.1409287
Received: 29 March 2024; Accepted: 03 June 2024;
Published: 02 July 2024.
Edited by:
Jun Zhang, Zhejiang Provincial People’s Hospital, ChinaReviewed by:
Joe Kodama, University of Maryland, United StatesMartha Elena Diaz Hernandez, Emory University, United States
Copyright © 2024 Liu, Ren, Zhang, Guo and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Fu, ZnVxaWFuZ3RqMjAyMUAxNjMuY29t
†These authors have contributed equally to this work