Deciphering the sequential changes of monocytes/macrophages in the progression of IDD with longitudinal approach using single-cell transcriptome

M1 macrophage-derived exosomes promote intervertebral disc degeneration by enhancing nucleus pulposus cell senescence through LCN2/NF-κB signaling axis

Intervertebral disc degeneration (IVDD) is the primary factor contributing to low back pain (LBP). Unlike elderly patients, many young IVDD patients usually have a history of trauma or long-term abnormal stress, which may lead to local inflammatory reaction causing by immune cells, and ultimately accelerates degeneration. Research has shown the significance of M1-type macrophages in IVDD; nevertheless, the precise mechanism and the route by which it influences the function of nucleus pulposus cell (NPC) remain unknown. Utilizing a rat acupuncture IVDD model and an NPC degeneration model induced by lipopolysaccharide (LPS), we investigated the function of M1 macrophage-derived exosomes (M1-Exos) in IVDD both in vivo and in vitro in this study. We found that M1-Exos enhanced LPS-induced NPC senescence, increased the number of SA-β-gal-positive cells, blocked the cell cycle, and promoted the activation of P21 and P53. M1-Exos derived from supernatant pretreated with the exosome inhibitor GW4869 reversed this result in vivo and in vitro. RNA-seq showed that Lipocalin2 (LCN2) was enriched in M1-Exos and targeted the NF-κB pathway. The quantity of SA-β-gal-positive cells was significantly reduced with the inhibition of LCN2, and the expression of P21 and P53 in NPCs was decreased. The same results were obtained in the acupuncture-induced IVDD model. In addition, inhibition of LCN2 promotes the expression of type II collagen (Col-2) and inhibits the expression of matrix metalloproteinase 13 (MMP13), thereby restoring the equilibrium of metabolism inside the extracellular matrix (ECM) in vitro and in vivo. In addition, the NF-κB pathway is crucial for regulating M1-Exo-mediated NPC senescence. After the addition of M1-Exos to LPS-treated NPCs, p-p65 activity was significantly activated, while si-LCN2 treatment significantly inhibited p-p65 activity. Therefore, this paper demonstrates that M1 macrophage-derived exosomes have the ability to deliver LCN2, which activates the NF-κB signaling pathway, and exacerbates IVDD by accelerating NPC senescence. This may shed new light on the mechanism of IVDD and bring a fresh approach to IVDD therapy.

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

SSR1 and CKAP4 as potential biomarkers for intervertebral disc degeneration based on integrated bioinformatics analysis

Background

Intervertebral disc degeneration (IDD) is a significant cause of low back pain and poses a significant public health concern. Genetic factors play a crucial role in IDD, highlighting the need for a better understanding of the underlying mechanisms.

Aim

The aim of this study was to identify potential IDD-related biomarkers using a comprehensive bioinformatics approach and validate them in vitro.

Materials and Methods

In this study, we employed several analytical approaches to identify the key genes involved in IDD. We utilized weighted gene coexpression network analysis (WGCNA), MCODE, LASSO algorithms, and ROC curves to identify the key genes. Additionally, immune infiltrating analysis and a single-cell sequencing dataset were utilized to further explore the characteristics of the key genes. Finally, we conducted in vitro experiments on human disc tissues to validate the significance of these key genes in IDD.

Results

we obtained gene expression profiles from the GEO database (GSE23130 and GSE15227) and identified 1015 DEGs associated with IDD. Using WGCNA, we identified the blue module as significantly related to IDD. Among the DEGs, we identified 47 hub genes that overlapped with the genes in the blue module, based on criteria of |logFC| ≥ 2.0 and p.adj <0.05. Further analysis using both MCODE and LASSO algorithms enabled us to identify five key genes, of which CKAP4 and SSR1 were validated by GSE70362, demonstrating significant diagnostic value for IDD. Additionally, immune infiltrating analysis revealed that monocytes were significantly correlated with the two key genes. We also analyzed a single-cell sequencing dataset, GSE199866, which showed that both CKAP4 and SSR1 were highly expressed in fibrocartilage chondrocytes. Finally, we validated our findings in vitro by performing real time polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) on 30 human disc samples. Our results showed that CKAP4 and SSR1 were upregulated in degenerated disc samples. Taken together, our findings suggest that CKAP4 and SSR1 have the potential to serve as disease biomarkers for IDD.

Comprehensive analysis of senescence-related genes and immune infiltration in intervertebral disc degeneration: a meta-data approach utilizing bulk and single-cell RNA sequencing data

Objectives: This study aims to identify the key senescence genes and potential regulatory mechanisms that contribute to the etiology of intervertebral disc degeneration (IDD). Method: We analyzed GSE34095 and GSE70362 datasets, identifying key senescence-related differentially expressed genes (DEGs) in IDD using lasso regression. Risk scores classified patients into high- and low-risk groups. We compared pathways, functions, and immune infiltration between these groups. Diagnostic ability was assessed using ROC curves and a nomogram predicted IDD incidence. In single-cell dataset GSE165722, we evaluated expression of key senescence-related DEGs. Results: We identified 12 key senescence-related DEGs distinguishing high- and low-risk IDD patients. Enrichment analysis revealed cellular stress response, apoptotic signaling pathway, and protein kinase activation differences. Immune cell analysis showed elevated eosinophils in low-risk group and increased effector memory CD8 T, central memory CD4 T, myeloid-derived suppressor, natural killer, monocyte, Type 1 T helper, plasmacytoid dendritic, and natural killer T cells in high-risk group. A nomogram using AUC >0.75 genes (CXCL8, MAP4K4, MINK1, and TNIK) predicted IDD incidence with good diagnostic power. High senescence scores were observed in neutrophils. Conclusion: Our diagnostic model, based on key senescence-related DEGs and immune cell infiltration, offers new insights into IDD pathogenesis and immunotherapy strategies.

Analysis of global research hotspots and trends in immune cells in intervertebral disc degeneration: A bibliometric study

Intervertebral disc degeneration is an important pathological basis for spinal degenerative diseases. The imbalance of the immune microenvironment and the involvement of immune cells has been shown to lead to nucleus pulposus cells death. This article presents a bibliometric analysis of studies on immune cells in IDD in order to clarify the current status and hotspots. We searched the WOSCC, Scopus and PubMed databases from 01/01/2001 to 08/03/2023. We analyzed and visualized the content using software such as Citespace, Vosviewer and the bibliometrix. This study found that the number of annual publications is increasing year on year. The journal study found that Spine had the highest number of articles and citations. The country/regions analysis showed that China had the highest number of publications, the USA had the highest number of citations and total link strength. The institutional analysis found that Shanghai Jiao Tong University and Huazhong University of Science Technology had the highest number of publications, Tokai University had the highest citations, and the University of Bern had the highest total link strength. Sakai D and Risbud MV had the highest number of publications. Sakai D had the highest total link strength, and Risbud MV had the highest number of citations. The results of the keyword analysis suggested that the current research hotspots and future directions continue to be the study of the mechanisms of immune cells in IDD, the therapeutic role of immune cells in IDD and the role of immune cells in tissue engineering for IDD.

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.

Integrated analysis of single-cell and bulk RNA sequencing data identifies the characteristics of ferroptosis in lumbar disc herniation

Lumbar disc herniation (LDH) is the most common condition associated with low back pain, and it adversely impacts individuals' health. Ferroptosis has recently emerged as a novel factor in the pathogenesis of LDH; however, the specific impacts of ferroptosis on LDH have not been fully elucidated. Ferroptosis-related differentially expressed genes (FRDEGs) were identified from the transcriptomic datasets of LDH. Gene set enrichment analysis (GSEA) was conducted to identify biological mechanism and related pathways. LASSO and SVM-RFE algorithms were applied to detect signature genes. Function of the signature gene was confirmed by RT-qPCR. The CIBERSORT algorithm was used to compare immune infiltration between LDH and normal samples. Correlation analysis between MYB and immune cells was analyzed using the Pearson method. Additionally, we used scRNA-seq to dissect cell clusters and cellular interactions. AUCell scoring was used to analyze the ferroptosis scores of different cell types. We found that MYB, a highly expressed ferroptosis-related gene, was associated with LDH By leveraging bioinformatics analysis. In immune infiltration analysis, the abundance of monocytes and macrophages varied significantly between the LDH group and disc spondylolisthesis (DS) group. MYB was correlated with most immune cells. GSEA revealed MYB was significantly enriched in immune-related pathways. Furthermore, scRNA-seq analysis revealed the presence of eight distinct cell types. AUCell analysis showed that macrophages had a high ferroptosis score. Cell trajectory analysis revealed that chondrocyte 1 was at the beginning of the trajectory, while calcification inhibiting chondrocytes and fibrochondrocytes accumulated along the middle and tail end of the trajectory, respectively. Cell-cell communication analysis identified chondrocyte 1 had an extensive communication network with other clusters and interacted with nucleus pulposus through collagen signaling pathway. Our analysis demonstrated that MYB may be a potential therapeutic target for LDH. This study provides a resource for the orthopedics community that will facilitate additional discoveries directedly toward understanding the pathogenesis process of LDH.