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

Identification and Validation of Circadian Rhythm-Related Genes Involved in Intervertebral Disc Degeneration and Analysis of Immune Cell Infiltration via Machine Learning

Background

Low back pain is a significant burden worldwide, and intervertebral disc degeneration (IVDD) is identified as the primary cause. Recent research has emphasized the significant role of circadian rhythms (CRs) and immunity in affecting intervertebral discs (IVD). However, the influence of circadian rhythms and immunity on the mechanism of IVDD remains unclear. This study aimed to identify and validate key rhythm-related genes in IVDD and analyze their correlation with immune cell infiltration.

Methods

Two gene expression profiles related to IVDD and rhythm-related genes were obtained from the Gene Expression Omnibus and GeneCards databases to identify differentially expressed rhythm-related genes (DERGs). Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA) were conducted to explore the biological functions of these genes. LASSO regression and SVM algorithms were employed to identify hub genes. We subsequently investigated the correlation between hub rhythm-related genes and immune cell infiltration. Finally, nucleus pulposus-derived mesenchymal stem cells (NPMSCs) were isolated from normal and degenerative human IVD tissues. Hub rhythm-related genes expression in NPMSCs was confirmed by real-time quantitative PCR (RT-qPCR).

Results

Six hub genes related to CRs (CCND1, FOXO1, FRMD8, NTRK2, PRRT1, and TFPI) were screened out. Immune infiltration analysis revealed that the IVDD group had significantly more M0 macrophages and significantly fewer follicular helper T cells than those of the control group. Specifically, M0 macrophages were significantly associated with FRMD8, PRRT1, and TFPI. T follicular helper cells were significantly associated with FRDM8, FOXO1, and CCND1. We further confirmed that CCND1, FRMD8, NTRK2, and TFPI were dysrhythmic within NPMSCs from degenerated IVD in vitro.

Conclusion

Six genes (CCND1, FOXO1, FRMD8, NTRK2, PRRT1 and TFPI) linked to circadian rhythms associated with IVDD progression, together with immunity. The identification of these DEGs may provide new insights for the diagnosis and treatment of IVDD.

Role of macrophage in intervertebral disc degeneration

Intervertebral disc degeneration is a degenerative disease where inflammation and immune responses play significant roles. Macrophages, as key immune cells, critically regulate inflammation through polarization into different phenotypes. In recent years, the role of macrophages in inflammation-related degenerative diseases, such as intervertebral disc degeneration, has been increasingly recognized. Macrophages construct the inflammatory microenvironment of the intervertebral disc and are involved in regulating intervertebral disc cell activities, extracellular matrix metabolism, intervertebral disc vascularization, and innervation, profoundly influencing the progression of disc degeneration. To gain a deeper understanding of the inflammatory microenvironment of intervertebral disc degeneration, this review will summarize the role of macrophages in the pathological process of intervertebral disc degeneration, analyze the regulatory mechanisms involving macrophages, and review therapeutic strategies targeting macrophage modulation for the treatment of intervertebral disc degeneration. These insights will be valuable for the treatment and research directions of intervertebral disc degeneration.

Single cell transcriptomics in blood of patients with chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide. Single-cell RNA sequencing (scRNA-seq) provides gene expression profiles at the single-cell level. Hence, we evaluated gene expression in the peripheral blood of patients with COPD.

METHODS

Peripheral blood samples from seven healthy controls and eight patients with COPD were obtained in this study. The 10X Genomics Chromium Instrument and cDNA synthesis kit were utilized to generate a barcoded cDNA library for single cell RNA-sequencing. We compared the scRNA-seq data between the COPD and control groups using computational analysis. Functional analyses were performed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses.

RESULTS

scRNA-seq was used to analyze the transcriptome of peripheral blood mononuclear cells from seven normal controls and eight patients with COPD. We found an increased number of monocyte/macrophages in the COPD group compared to the normal control group. Among the differentially expressed genes (DEGs) in monocyte/macrophages, we identified 15 upregulated genes (EGR1, NR4A1, CCL3, CXCL8, PTGS2, CD83, BCL2A1, SGK1, IL1B, BTG2, NFKBIZ, DUSP2, MAFB, PLAUR and CCL3L1) and 7 downregulated genes (FOLR3, RPS4Y1, HLA-DRB5, NAMPT, CD52, TMEM176A and TMEM176B) in the COPD group compared to the normal control group.

CONCLUSIONS

Using scRNA-seq, we found differences in cell type distribution, especially in monocyte/ macrophages. Several upregulated and downregulated genes were found in the monocyte/macrophages of the COPD group.

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.

"Dictionary of immune responses" reveals the critical role of monocytes and the core target IRF7 in intervertebral disc degeneration

Background

Intervertebral disc degeneration (IDD) is widely regarded as the primary contributor to low back pain(LBP). As an immune-privileged organ, upon the onset of IDD, various components of the nucleus pulposus (NP) are exposed to the host's immune system, accumulating cytokines. Cytokines facilitate intercellular communication within the immune system, induce immune cells polarisation, and exacerbate oxidative stress in IDD.

Methods

Machine learning was used to identify crucial immune cells. Subsequently, Immune Response Enrichment Analysis (IREA) was conducted on the key immune cells to determine their cytokine responses and polarisation states in IDD. "CellChat" package facilitated the analysis of cell-cell communication. Differential gene expression analysis, PPI network, GO and KEGG pathway enrichment analysis, GSVA, co-expressed gene analysis and key gene-related networks were also performed to explore hub genes and their associated functions. Lastly, the differential expression and functions of key genes were validated through in vitro and in vivo experiments.

Results

Through multiple machine learning methods, monocytes were identified as the crucial immune cells in IDD, exhibiting significant differentiation capacity. IREA revealed that monocytes in IDD polarize into an IFN-a1 and IFN-b enriched Mono-a state, potentially intensifying inflammation. Cell-cell communication analysis uncovered alteration in ANNEXIN pathway and a reduction in CXCL signaling between macrophages and monocytes, suggesting immune response dysregulation. Furthermore, ten algorithms identified three hub genes. Both experiments conducted in vitro and in vivo have conclusively shown that IRF7 serves as a crucial target for the treatment of IDD, and its knockdown alleviates IDD. Eight small-molecule drugs were predicted to have therapeutic potential for IDD.

Conclusion

These findings offer a multidimensional understanding of the pathogenesis of IDD, pinpointing monocytes and key genes as potential diagnostic and therapeutic targets. They provide novel insights into potential diagnostic and therapeutic targets for IDD.

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 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.