A microarray data analysis investigating the pathogenesis and potential biomarkers of autophagy and ferroptosis in intervertebral disc degeneration.

Background: Intervertebral disc degeneration (IDD) entails complex pathological changes and causes lower back pain (LBP). However, there is still a lack of understanding of the mechanisms involved in IDD, particularly regarding the roles of autophagy and ferroptosis. The current study used microarray data to investigate the pathogenesis of IDD and potential biomarkers related to autophagy and ferroptosis in IDD. Methods: Differentially expressed genes (DEGs) were identified by analyzing the mRNA and miRNA expression profiles of IDD patients from the Gene Expression Omnibus (GEO). The protein-protein interaction network, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and gene set enrichment analysis (GSEA) were utilized. The Human Autophagy Database (HADb) and Ferroptosis Database were used in conjunction with hub genes to identify autophagy- and ferroptosis-related genes. The Transcription Factor -hub gene-miRNA network was constructed. Lastly, the expression of DEGs in normal and degenerated nucleus pulposus cells (NPCs) was investigated via the quantitative reverse transcription polymerase chain reaction (qRT-PCR). Results: A total of 362 DEGs associated with IDD were identified. GO and KEGG analyses indicated that oxidative stress, extracellular matrix, PI3K-AKT signaling pathway, and ferroptosis were key factors in IDD occurrence. GSEA indicated that IDD was associated with changes in autophagy, iron ion homeostasis, extracellular matrix, and oxidative stress. Eighty-nine hub genes were obtained, including five that were autophagy-related and three that were ferroptosis-related. Of these, TP53 and SESN2 were the intersections of autophagy- and ferroptosis-related genes. In qRT-PCR analysis, CANX, SLC38A1, and TP53 were downregulated in degenerative NPCs, whereas GNAI3, SESN2, and VAMP3 were upregulated. Conclusion: The current study revealed aspects of autophagy- and ferroptosis-related genes involved in IDD pathogenesis, warranting further investigation.

Kuang W ，Jiang C ，Yu C ，Hu J ，Duan Y ，Chen Z ... -  《Frontiers in Genetics》

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

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.

Li J ，Yu C ，Ni S ，Duan Y ... -  《Frontiers in Genetics》

MAPK8 and CAPN1 as potential biomarkers of intervertebral disc degeneration overlapping immune infiltration, autophagy, and ceRNA.

Intervertebral disc degeneration (IDD) is one of the most common health problems in the elderly and a major causative factor in low back pain (LBP). An increasing number of studies have shown that IDD is closely associated with autophagy and immune dysregulation. Therefore, the aim of this study was to identify autophagy-related biomarkers and gene regulatory networks in IDD and potential therapeutic targets. We obtained the gene expression profiles of IDD by downloading the datasets GSE176205 and GSE167931 from the Gene Expression Omnibus (GEO) public database. Subsequently, differentially expressed genes (DEGs) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, gene ontology (GO), and gene set enrichment analysis (GSEA) were performed to explore the biological functions of DEGs. Differentially expressed autophagy-related genes (DE-ARGs) were then crossed with the autophagy gene database. The hub genes were screened using the DE-ARGs protein-protein interaction (PPI) network. The correlation between the hub genes and immune infiltration and the construction of the gene regulatory network of the hub genes were confirmed. Finally, quantitative PCR (qPCR) was used to validate the correlation of hub genes in a rat IDD model. We obtained 636 DEGs enriched in the autophagy pathway. Our analysis revealed 30 DE-ARGs, of which six hub genes (MAPK8, CTSB, PRKCD, SNCA, CAPN1, and EGFR) were identified using the MCODE plugin. Immune cell infiltration analysis revealed that there was an increased proportion of CD8+ T cells and M0 macrophages in IDD, whereas CD4+ memory T cells, neutrophils, resting dendritic cells, follicular helper T cells, and monocytes were much less abundant. Subsequently, the competitive endogenous RNA (ceRNA) network was constructed using 15 long non-coding RNAs (lncRNAs) and 21 microRNAs (miRNAs). In quantitative PCR (qPCR) validation, two hub genes, MAPK8 and CAPN1, were shown to be consistent with the bioinformatic analysis results. Our study identified MAPK8 and CAPN1 as key biomarkers of IDD. These key hub genes may be potential therapeutic targets for IDD.

Zhang Y ，Zhang J ，Sun Z ，Wang H ，Ning R ，Xu L ，Zhao Y ，Yang K ，Xi X ，Tian J ... -  《Frontiers in Immunology》

Bioinformatics Research and qRT-PCR Verify Hub Genes and a Transcription Factor-MicroRNA Feedback Network in Intervertebral Disc Degeneration.

The present study explores the potentials of bioinformatics analysis to identify hub genes linked to intervertebral disc degeneration (IDD) and explored the potential molecular mechanism of transcription factor-microRNA regulatory network. Furthermore, the hub genes were identified through quantitative reverse transcriptase PCR (qRT-PCR). GEO database expression profile datasets for candidate genes (GSE124272) were downloaded. Genes that were differentially expressed (DEGs) were detected utilizing limma technique in the R programming language. Search Tool for the Retrieval of Interacting Genes/Proteins and NetworkAnalyst software identified hub genes. The Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis as well as Gene Ontology annotation of the DEGs were performed using Metascape. Using Bioinformatics data from the TRRUST, StarBase, and TransmiR databases, a TF-miRNA-hub genes network was constructed. qRT-PCR was utilized to confirm the result. As compared to healthy persons, 521 DEGs, comprising 203 down-regulated and 318 up-regulated genes, as well as 7 core genes, were found in people with IDD. Analysis revealed that all seven essential genes were under-expressed. qRT-PCR further confirmed the low expression of these seven important genes. Based on the TRRUST database, 16 TFs that could target five junction genes were then predicted. According to the StarBase database, four miRNAs were linked to crucial genes, while the TransmiR database predicted regulatory connections between four miRNAs and five TFs. The expression of the TP53-(hsa-miR-183-5p)-CCNB1 TF-miRNA-mRNA interaction network was discovered to be correlated with IDD. Throughout this investigation, a network of TF-miRNA-mRNA connections was built for investigation of the probable molecular mechanisms responsible for IDD. The identification of hub genes associated with IDD may reveal promising IDD treatment strategies.

Liu M ，Li H ，Huo Z ，Chen H ，Kang X ，Xu B ... -  《-》

Identification and validation of ferroptosis-related gene signature in intervertebral disc degeneration.

Lower back pain (LBP) is a leading cause of disability in the elderly and intervertebral disc degeneration (IDD) is the major contributor to LBP. Ferroptosis is a newly discovered programmed cell death, characterized by iron-dependent lethal lipid peroxidation. Growing evidence has shown that ferroptosis plays important roles in various human diseases. However, the underlying mechanism of ferroptosis in IDD remains elusive. This study is aimed to uncover the key roles of ferroptosis in the pathogenesis and progression of IDD comprehensively. To investigate the ferroptosis related differentially expressed genes (FRDEGs) in IDD, we analyzed the microarray data from the Gene Expression Omnibus (GEO) database. Then we performed functional enrichment analysis and protein-protein interaction (PPI) network analysis, and screened out the hub FRDEGs. To further evaluate the predictive value of these hub FRDEGs, we performed ROC analysis based on the GSE124272 dataset. A total of 80 FRDEGs were identified, including 20 downregulated and 60 upregulated FRDEGs. The FRDEGs were primarily involved in the biological processes of response to chemical, and response to stress. KEGG pathway enrichment analysis showed that the FRDEGs were mainly involved in ferroptosis, TNF signaling pathway, HIF-1 signaling pathway, NOD-like receptor signaling pathway, and IL-17 signaling pathway. Ten hub OSRDEGs were obtained according to the PPI analysis, including HMOX1, KEAP1, MAPK1, HSPA5, TXNRD1, IL6, PPARA, JUN, HIF1A, DUSP1. The ROC analysis and RT-qPCR validation results suggested that most of the hub FRDEGs might be potential signature genes for IDD. This study reveals that ferroptosis might provide promising strategy for the diagnosis and treatment of IDD.

Xiang Q ，Zhao Y ，Li W 《Frontiers in Endocrinology》