Abstract:Background: 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. Methods: 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). Results: 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. Conclusion: Our findings may provide new insights into the functional characteristics and mechanism of IDD and aid the development of IDD therapeutics.

Screening of differentially expressed microRNAs in degenerative intervertebral disc and investigation about the involvement of JNK pathway in IVDD

Identification of Cervical Disc Degeneration and Genes Associated with Immune Infiltration

Identification of microRNAs and target genes associated with the development of human intervertebral disc degeneration by bioinformatic analysis and verification

Screening and identification of microRNAs mediating cartilage endplate degeneration in human intervertebral disks

Comprehensive evaluation of differential lncRNA and gene expression in patients with intervertebral disc degeneration.

Screening of miR-15a-5p as a potential biomarker for intervertebral disc degeneration through RNA-sequencing

RETRACTED: Inhibited microRNA-494-5p promotes proliferation and suppresses senescence of nucleus pulposus cells in mice with intervertebral disc degeneration by elevating TIMP3 (Retracted Article)

The Role of Mir-31-5p in the Development of Intervertebral Disc Degeneration and Its Therapeutic Potential.

Characterization of Microrna Expression Profiles in Patients with Intervertebral Disc Degeneration

Identifying critical modules and biomarkers of intervertebral disc degeneration by using weighted gene co-expression network

[MicroRNAs: a type of novel regulative factor for intervertebral disc degeneration].

The Role of the Mir-4306/Pak6 Axis in Degenerative Nucleus Pulposus Cells in Human Intervertebral Disc Degeneration

Identification of the Potential Molecular Targets for Human Intervertebral Disc Degeneration Based on Bioinformatic Methods.

Bioinformatics Analysis of Molecular Mechanisms Involved in Intervertebral Disc Degeneration Induced by TNF-α and IL-1β.

Gene expression profile identifies potential biomarkers for human intervertebral disc degeneration.

Identification of novel gene signatures and immune cell infiltration in intervertebral disc degeneration using bioinformatics analysis

Nucleus pulposus related lncRNA and mRNA expression profiles in intervertebral disc degeneration

Identification of miRNA-mRNA Pairs in Relation to TNF- α/IL-1 β Induced Inflammatory Response in Intervertebral Disc Degeneration

[Gene Expression Profiles of Long Non-Coding RNAs in Human Degenerated Intervertebral Disc Tissue].

Potential Role of lncRNAs in Contributing to Pathogenesis of Intervertebral Disc Degeneration Based on Microarray Data.

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