Derivation and comprehensive analysis of ageing-related genes in intervertebral disc degeneration for prediction and immunology

Study on molecular mechanism of intervertebral disc degeneration by single cell hdWGCNA combined with transcriptome sequencing

Background

Intervertebral disc degeneration (IVDD) is one of the important causes of lower back pain, seriously affecting people's health and quality of life. This research employs single-cell analysis to identify the specific cellular subtypes and key regulatory genes associated with IVDD.

Methods

We analyzed the single-cell data and screened cells that closely associated with the development of IVDD. The differential expression of feature genes between IVDD and control groups was analyzed. Additionally, drugs and regulatory transcription factors that interact with feature genes were predicted and clinically validated by reverse transcription quantitative real-time PCR (RT-qPCR), immunohistochemistry (IHC), and enzyme-linked immunosorbent assay (ELISA).

Results

Our study identified the Chond2 cell subtype associated with IVDD and selected four feature genes influencing the development of IVDD, namely IGFBP3, ACAN, VAPA and TMEM45A, through the high-dimensional weighted gene co-expression network analysis (hdWGCNA) analysis, least absolute shrinkage and selection operator (LASSO), and random forest (RF). Besides, compared to the MDD group, IGFBP3 and TMEM45A were significantly upregulated in the SDD group, while ACAN and VAPA showed no significant difference between the two groups. ELISA testing revealed a positive correlation between IGFBP3 concentration and the grading of IVDD. Furthermore, Celecoxib may be used to treat IVDD by inhibiting IGFBP3.

Conclusion

Our study identified the Chond2 cell subtype associated with IVDD and selected four feature genes influencing the development of IVDD, namely IGFBP3, ACAN, VAPA and TMEM45A. Our findings establish a robust theoretical foundation for the clinical diagnosis and treatment of IVDD patients.

A novel proteomic prognostic signature characterizes the immune landscape and predicts nasopharyngeal carcinoma prognosis

Background

Nasopharyngeal carcinoma (NPC) is a highly diverse and aggressive cancer type, leading to varying prognoses and responses to immunotherapy. This study aims to develop a protein-based signature that provides new insights into assessing the prognosis and immunotherapeutic response in NPC patients.

Methods and Results

We obtained transcriptomic and proteomic data for NPC from TCGA and CPTAC databases, respectively. Differentially expressed proteins with prognostic significance were identified using limma combined with uniCox analysis. A prognostic protein signature was created utilizing the LASSO algorithm. Receiver operating characteristic (ROC) curve analysis along with Kaplan-Meier survival analysis was conducted to assess the predictive accuracy of this signature. To evaluate immune infiltration levels among patients categorized by high or low risk scores (RPscores), we employed ssGSEA and ESTIMATE methods, while TIDE was used to forecast responses to immunotherapy. Our research pinpointed four critical prognostic proteins: CdSTA, AGR3, DUSP14, and LRRC17, allowing us to compute risk scores (RPscores). Kaplan-Meier curves demonstrated that individuals in the low-risk category exhibited better survival rates. Furthermore, RPscore effectively predicted overall survival across both training and testing cohorts. The ssGSEA results indicated that RPscore is linked with an immune-suppressive microenvironment correlating with diminished immune responses. Notably, DUSP14 showed significant upregulation in NPC cases; its role in promoting cell invasion and metastasis was confirmed through in vitro studies.

Conclusion

We have established a robust protein-related signature capable of accurately forecasting prognosis as well as immunotherapy outcomes for NPC patients. Moreover, DUSP14 emerged as a potential therapeutic target due to its strong association with patient prognosis in nasopharyngeal carcinoma.

Effective delivery of miR-150-5p with nucleus pulposus cell-specific nanoparticles attenuates intervertebral disc degeneration

BACKGROUND

The use of gene therapy to deliver microRNAs (miRNAs) has gradually translated to preclinical application for the treatment of intervertebral disc degeneration (IDD). However, the effects of miRNAs are hindered by the short half-life time and the poor cellular uptake, owing to the lack of efficient delivery systems. Here, we investigated nucleus pulposus cell (NPC) specific aptamer-decorated polymeric nanoparticles that can load miR-150-5p for IDD treatment.

METHODS

The role of miR-150-5p during disc development and degeneration was examined by miR-150-5p knockout (KO) mice. Histological analysis was undertaken in disc specimens. The functional mechanism of miR-150-5p in IDD development was investigated by qRT-PCR assay, Western blot, coimmunoprecipitation and immunofluorescence. NPC specific aptamer-decorated nanoparticles was designed, and its penetration, stability and safety were evaluated. IDD progression was assessed by radiological analysis including X-ray and MRI, after the annulus fibrosus needle puncture surgery with miR-150-5p manipulation by intradiscal injection of nanoparticles. The investigations into the interaction between aptamer and receptor were conducted using mass spectrometry, molecular docking and molecular dynamics simulations.

RESULTS

We investigated NPC-specific aptamer-decorated polymeric nanoparticles that can bind to miR-150-5p for IDD treatment. Furthermore, we detected that nanoparticle-loaded miR-150-5p inhibitors alleviated NPC senescence in vitro, and the effects of the nanoparticles were sustained for more than 3 months in vivo. The microenvironment of NPCs improves the endo/lysosomal escape of miRNAs, greatly inhibiting the secretion of senescence-associated factors and the subsequent degeneration of NPCs. Importantly, nanoparticles delivering miR-150-5p inhibitors attenuated needle puncture-induced IDD in mouse models by targeting FBXW11 and inhibiting TAK1 ubiquitination, resulting in the downregulation of NF-kB signaling pathway activity.

CONCLUSIONS

NPC-targeting nanoparticles delivering miR-150-5p show favorable therapeutic efficacy and safety and may constitute a promising treatment for IDD.

ATG9A as a potential diagnostic marker of intervertebral disc degeneration: Inferences from experiments and bioinformatics analysis incorporating sc-RNA-seq data

BACKGROUND

Disfunctional autophagy plays a pivotal role in Intervertebral Disc Degeneration (IDD) progression. however, the connection between Autophagy-related gene 9A (ATG9A) and IDD has not been reported.

METHODS

Firstly, transcriptome datasets from the GEO and Autophagy-related genes (ARGs) from GeneCards were carried out using R. Following this, IDD-specific signature genes were identified through methods such as least absolute shrinkage and selection operator (LASSO), random forest (RF), and support vector machine (SVM) analyses. Validation of these findings proceeded through in vitro experiments, evaluation of independent datasets, and analysis of receiver operating characteristic (ROC) curves. Subsequent steps incorporated co-expression analysis, Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, Gene Set Enrichment Analysis (GSEA), and construction of competing endogenous RNA (ceRNA) network. The final section established the correlation between immune cell infiltration, ATG9A, and IDD utilizing the CIBERSORT algorithm and single-cell RNA (scRNA) sequencing data.

RESULTS

Research identified 87 differentially expressed genes, with only ATG9A noted as an IDD signature gene. Analysis of in vitro experiments and independent datasets uncovered a decrease in ATG9A expression within the degeneration group. The area under the curve (AUC) of ATG9A exceeded 0.8 following ROC analysis. Furthermore, immune cell infiltration and scRNA sequencing data analysis elucidated the substantial role of immune cells in IDD progression. A ceRNA network was constructed, centered around ATG9A, included 4 miRNAs and 22 lncRNAs.

CONCLUSION

ATG9A was identified as a diagnostic gene for IDD, indicating its viability as a effective target for therapy disease.

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.