Identification of Core Genes and Screening of Potential Targets in Intervertebral Disc Degeneration Using Integrated Bioinformatics Analysis

Overview

Journal Front Genet

Date 2022 Jun 17

PMID 35711934

Authors

Jianjun Li

Cheng Yu

Songjia Ni

Yang Duan

Affiliations

Soon will be listed here.

Abstract

Intervertebral disc degeneration (IDD), characterized by diverse pathological changes, causes low back pain (LBP). However, prophylactic and delaying treatments for IDD are limited. The aim of our study was to investigate the gene network and biomarkers of IDD and suggest potential therapeutic targets. Differentially expressed genes (DEGs) associated with IDD were identified by analyzing the mRNA, miRNA, and lncRNA expression profiles of IDD cases from the Gene Expression Omnibus (GEO). The protein-protein interaction (PPI) network, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis as well as miRNA-lncRNA-mRNA networks were conducted. Moreover, we obtained 71 hub genes and performed a comprehensive analysis including GO, KEGG, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), Disease Ontology (DO), methylation analysis, receiver operating characteristic (ROC) curve analysis, immune infiltration analysis, and potential drug identification. We finally used qRT-PCR to verify 13 significant DEGs in normal and degenerative nucleus pulposus cells (NPCs). We identified 305 DEGs closely related to IDD. The GO and KEGG analyses indicated that changes in IDD are significantly associated with enrichment of the inflammatory and immune response. GSEA analysis suggested that cell activation involved in the inflammatory immune response amide biosynthetic process was the key for the development of IDD. The GSVA suggested that DNA repair, oxidative phosphorylation, peroxisome, IL-6-JAK-STAT3 signaling, and apoptosis were crucial in the development of IDD. Among the 71 hub genes, the methylation levels of 11 genes were increased in IDD. A total of twenty genes showed a high functional similarity and diagnostic value in IDD. The result of the immune cell infiltration analysis indicated that seven genes were closely related to active natural killer cells. The most relevant targeted hub genes for potential drug or molecular compounds were MET and PIK3CD. Also, qRT-PCR results showed that ARHGAP27, C15orf39, DEPDC1, DHRSX, MGAM, SLC11A1, SMC4, and LINC00887 were significantly downregulated in degenerative NPCs; H19, LINC00685, mir-185-5p, and mir-4306 were upregulated in degenerative NPCs; and the expression level of mir-663a did not change significantly in normal and degenerative NPCs. Our findings may provide new insights into the functional characteristics and mechanism of IDD and aid the development of IDD therapeutics.

Citing Articles

Integrated bulk and single-cell RNA sequencing to identify potential biomarkers in intervertebral disc degeneration.

Fan C, Xu W, Li X, Wang J, He W, Shen M Eur J Med Res. 2025; 30(1):102.

PMID: 39953636 PMC: 11827443. DOI: 10.1186/s40001-025-02346-4.

Assessment of XCI skewing and demonstration of XCI escape region based on single-cell RNA sequencing: comparison between female Grave's disease and control.

Baek I, Sim S, Suh B, Kim T, Cho W BMC Mol Cell Biol. 2025; 26(1):8.

PMID: 39891056 PMC: 11786500. DOI: 10.1186/s12860-025-00533-z.

Deciphering SPP1-related macrophage signaling in the pathogenesis of intervertebral disc degeneration.

Yu X, Zou P, Li T, Bai X, Wang S, Guan J Cell Biol Toxicol. 2025; 41(1):33.

PMID: 39825191 PMC: 11748470. DOI: 10.1007/s10565-024-09948-4.

A diagnostic signatures for intervertebral disc degeneration using TNFAIP6 and COL6A2 based on single-cell RNA-seq and bulk RNA-seq analyses.

Song H, Li Z, Huang Q, Wu C, Li M, Shen J Ann Med. 2024; 57(1):2443568.

PMID: 39704340 PMC: 11703038. DOI: 10.1080/07853890.2024.2443568.

Experimental confirmation and bioinformatics reveal biomarkers of immune system infiltration and hypertrophy ligamentum flavum.

Liu F, Zhong M, Yang L, Song C, Chen C, Xu Z JOR Spine. 2024; 7(3):e1354.

PMID: 39071860 PMC: 11272949. DOI: 10.1002/jsp2.1354.

References

Stich S, Jagielski M, Fleischmann A, Meier C, Bussmann P, Kohl B . Degeneration of Lumbar Intervertebral Discs: Characterization of Anulus Fibrosus Tissue and Cells of Different Degeneration Grades. Int J Mol Sci. 2020; 21(6). PMC: 7139657. DOI: 10.3390/ijms21062165. View

Liu X, Che L, Xie Y, Hu Q, Ma C, Pei Y . Noncoding RNAs in human intervertebral disc degeneration: An integrated microarray study. Genom Data. 2015; 5:80-1. PMC: 4583642. DOI: 10.1016/j.gdata.2015.05.027. View

Wang Y, Dai G, Jiang L, Liao S, Xia J . Microarray analysis reveals an inflammatory transcriptomic signature in peripheral blood for sciatica. BMC Neurol. 2021; 21(1):50. PMC: 7856817. DOI: 10.1186/s12883-021-02078-y. View

Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov J, Tamayo P . The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2016; 1(6):417-425. PMC: 4707969. DOI: 10.1016/j.cels.2015.12.004. View

Szklarczyk D, Gable A, Lyon D, Junge A, Wyder S, Huerta-Cepas J . STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2018; 47(D1):D607-D613. PMC: 6323986. DOI: 10.1093/nar/gky1131. View

Yang Y, Zhang B, Nakata M, Nakae J, Mori M, Yada T . Islet β-cell-produced NUCB2/nesfatin-1 maintains insulin secretion and glycemia along with suppressing UCP-2 in β-cells. J Physiol Sci. 2019; 69(5):733-739. PMC: 10717817. DOI: 10.1007/s12576-019-00689-2. View

Wang J, Hu J, Chen X, Huang C, Lin J, Shao Z . BRD4 inhibition regulates MAPK, NF-κB signals, and autophagy to suppress MMP-13 expression in diabetic intervertebral disc degeneration. FASEB J. 2019; 33(10):11555-11566. DOI: 10.1096/fj.201900703R. View

Newman A, Steen C, Liu C, Gentles A, Chaudhuri A, Scherer F . Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019; 37(7):773-782. PMC: 6610714. DOI: 10.1038/s41587-019-0114-2. View

Liu L, He J, Liu C, Yang M, Fu J, Yi J . Cartilage intermediate layer protein affects the progression of intervertebral disc degeneration by regulating the extracellular microenvironment (Review). Int J Mol Med. 2021; 47(2):475-484. PMC: 7797476. DOI: 10.3892/ijmm.2020.4832. View

10.

Bakula D, Muller A, Zuleger T, Takacs Z, Franz-Wachtel M, Thost A . WIPI3 and WIPI4 β-propellers are scaffolds for LKB1-AMPK-TSC signalling circuits in the control of autophagy. Nat Commun. 2017; 8:15637. PMC: 5460038. DOI: 10.1038/ncomms15637. View

11.

Gallardo-Montejano V, Yang C, Hahner L, McAfee J, Johnson J, Holland W . Perilipin 5 links mitochondrial uncoupled respiration in brown fat to healthy white fat remodeling and systemic glucose tolerance. Nat Commun. 2021; 12(1):3320. PMC: 8175597. DOI: 10.1038/s41467-021-23601-2. View

12.

Ren H, Yu X, Shen G, Zhang Z, Shang Q, Zhao W . miRNA-seq analysis of human vertebrae provides insight into the mechanism underlying GIOP. Bone. 2018; 120:371-386. DOI: 10.1016/j.bone.2018.11.013. View

13.

Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A . ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009; 25(8):1091-3. PMC: 2666812. DOI: 10.1093/bioinformatics/btp101. View

14.

Pang W, Hu F . Cellular and physiological functions of C9ORF72 and implications for ALS/FTD. J Neurochem. 2020; 157(3):334-350. PMC: 8842544. DOI: 10.1111/jnc.15255. View

15.

Fan Y, Xia J . miRNet-Functional Analysis and Visual Exploration of miRNA-Target Interactions in a Network Context. Methods Mol Biol. 2018; 1819:215-233. DOI: 10.1007/978-1-4939-8618-7_10. View

16.

Han S, Zhang Y, Zhang X, Zhang H, Meng S, Kong M . Single-Cell RNA Sequencing of the Nucleus Pulposus Reveals Chondrocyte Differentiation and Regulation in Intervertebral Disc Degeneration. Front Cell Dev Biol. 2022; 10:824771. PMC: 8899542. DOI: 10.3389/fcell.2022.824771. View

17.

Wang H, Yu X, Liu Z, Cheng X, Samartzis D, Jia L . Deregulated miR-155 promotes Fas-mediated apoptosis in human intervertebral disc degeneration by targeting FADD and caspase-3. J Pathol. 2011; 225(2):232-42. DOI: 10.1002/path.2931. View

18.

Wang H, Liu W, Yu B, Yu X, Chen B . Identification of Key Modules and Hub Genes of Annulus Fibrosus in Intervertebral Disc Degeneration. Front Genet. 2021; 11:596174. PMC: 7875098. DOI: 10.3389/fgene.2020.596174. View

19.

Li C, Ma X, Ni C, Xu J, Xie Y, Kan J . LncRNA NEAT1 promotes nucleus pulposus cell matrix degradation through regulating Nrf2/ARE axis. Eur J Med Res. 2021; 26(1):11. PMC: 7818737. DOI: 10.1186/s40001-021-00481-2. View

20.

Tonomura H, Nagae M, Takatori R, Ishibashi H, Itsuji T, Takahashi K . The Potential Role of Hepatocyte Growth Factor in Degenerative Disorders of the Synovial Joint and Spine. Int J Mol Sci. 2020; 21(22). PMC: 7698933. DOI: 10.3390/ijms21228717. View