Identification of Biomarkers, Pathways, Immune Properties of Mitophagy Genes, and Prediction Models for Intervertebral Disc Degeneration

• Published in: Journal of inflammation research 2024-05, Vol.17, p.2959-2975
• Main Authors: Huang, Yongxiong, Qiu, Xianshuai, Liu, Jinlian, Wan, Jiangtao, Yu, Cheng, Liu, Chun, Duan, Yang, Chen, Chong, Dai, Jingxing, Ouyang, Jun, Liu, Ming, Min, Shaoxiong, Qiu, Sujun
• Format: Article
• Language: English
• Subjects: Backache; Biological markers; diagnostic model; Genes; immune infiltration; intervertebral disc degeneration; Medical research; Medicine, Experimental; MicroRNA; mitophagy-related differentially expressed genes; nucleus pulposus cells; Original Research; Physiological aspects; sqstm1
• ISSN: 1178-7031; 1178-7031
• DOI: 10.2147/JIR.S461668

Intervertebral disc degeneration (IDD) is the leading cause of low back pain (LBP). The mechanism of IDD development and progression is not fully understood. Peripheral biomarkers are increasingly vital non-radioactive methods in early detection and diagnosis for IDD. Nevertheless, less attention has been paid to the role of mitophagy genes in the progress of IDD. This study aimed to identify the mitophagy disease-causing genes in the process of IDD and mitophagy diagnostic biomarkers for IDD. Mitophagy-related differentially expressed genes (MRDEGs) related to IDD were investigated by analyzing the microarray datasets of IDD cases from GEO, PathCards and Molecular Signatures Databases. We used R software, WGCNA, PPI, mRNA-miRNA, mRNA-TF, GO, KEGG, GSEA, GSVA and Cytoscape to analyze and visualize the data. We further used ssGSEA for immunoinfiltration analysis to obtain different immune cell infiltration. LASSO model was developed to screen for genes that met the diagnostic gene model requirements. Finally, qRT-PCR, Western blotting and HE were used to verify hub genes and their expression from clinical IDD samples. We identified 14 MRDEGs and 12 hub genes. GO, KEGG, GSEA and GSVA analyses demonstrated that hub genes were critical for the development of IDD. LASSO diagnostic model consisted of six hub genes, among which SQSTM1, ATG7 and OPTN were significantly different between the two IDD disease subtypes. At the same time, SQSTM1 also had a high correlation with immune characteristic subtypes. The results of qRT-PCR and Western blotting also indicated that these genes were significantly differentially expressed in nucleus pulposus cells (NPCs) of the IDD group. We explored an association between MRDEGs-associated signature in IDD and validated that hub genes like SQSTM1 might serve as biomarkers for diagnostic and therapeutic targets for IDD. Meanwhile, this study can provide new insights into the functional characteristics and mechanism of mitophagy in the development of IDD.