Integrative Bioinformatics Analysis Revealed Mitochondrial Dysfunction-Related Genes Underlying Intervertebral Disc Degeneration

Abstract

Objective. Mitochondrial dysfunction plays an important role in intervertebral disc degeneration (IDD). We aim to explore the pathways and key genes that cause mitochondrial dysfunction during IDD and to further reveal the pathogenesis of IDD based on bioinformatic analyses. Methods. Datasets GSE70362 and GSE124272 were downloaded from the Gene Expression Omnibus. Differentially expressed genes (DEGs) of mitochondrial dysfunction between IDD patients and healthy controls were screened by package limma package. Critical genes were identified by adopting gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) pathways, and protein-protein interaction (PPI) networks. We collected both degenerated and normal disc tissues obtained surgically, and we performed western blot and qPCR to verify the key DEGs identified in intervertebral disc tissues. Results. In total, 40 cases of IDD and 24 healthy controls were included. We identified 152 DEGs, including 67 upregulated genes and 85 downregulated genes. Four genes related to mitochondrial dysfunction (SOX9, FLVCR1, NR5A1 and UCHL1) were screened out. Of them, SOX9, FLVCR1, and UCHL1 were down-regulated in peripheral blood and intervertebral disc tissues of IDD patients, while NR5A1 was up-regulated. The analysis of immune infiltration showed the concentrations of mast cells activated were significantly the highest in IDD patients. Compared with the control group, the level of T cells CD4 memory resting was the lowest in the patients. In addition, 24 cases of IDD tissues and 12 cases of normal disc tissues were obtained to verify the results of bioinformatics analysis. Both western blot and qPCR results were consistent with the results of bioinformatics analysis. Conclusion. We identified four genes (SOX9, FLVCR1, NR5A1 and UCHL1) associated with mitochondrial dysfunction that play an important role in the progress of disc degeneration. The identification of these differential genes may provide new insights for the diagnosis and treatment of IDD.