Identification of Potential Ferroptosis Key Genes in the Pathogenesis of Lumbosacral Spinal Root Avulsion by RNA Sequencing and Bioinformatics Analysis

Objective: Ferroptosis is a type of cell death involved in various human diseases, including nerve injury. However, the role of ferroptosis in lumbosacral spinal root avulsion (LSRA) remains unknown. This study aims to investigate whether ferroptosis is induced after LSRA and the key ferroptosis-related genes and their potential function in LSRA. Methods: The biochemical and morphological changes of ferroptosis were determined by detection of iron accumulation and by transmission electron microscopy in a rat LSRA model. The transcriptional expression profile following LSRA was investigated by RNA sequencing and ferroptosis-related genes were downloaded from FerrDb and used to identify ferroptosis differentially expressed genes (DEGs). The differential expressions of ferroptosis DEGs were confirmed by qRT-PCR analysis. The potential functions of ferroptosis DEGs were revealed by DAVID 6.8 and WebGestalt. A protein–protein interaction (PPI) network and gene–miRNA interaction network were further constructed to identify key modules in ferroptosis DEGs, and the results were verified by qRT-PCR and western blot analysis. Results: LSRA was followed by ferroptosis-specific changes, such as shrunken mitochondria and increased iron accumulation, that can be alleviated by ferroptosis inhibitor deferoxamine (DFO). A total of 2,446 DEGs and 46 ferroptosis DEGs were identified after LSRA, and over 90% of the ferroptosis DEGs were confirmed to be differentially expressed following LSRA, which can also be eliminated by DFO treatment. Functional analysis demonstrated significant enrichment of the ferroptosis DEGs in pathways related to the oxidative stress response, the HIF-1 signaling pathway, and the tumor necrosis factor signaling pathway. PPI network analysis demonstrated that a set of key modules in ferroptosis DEGs were related to the HIF-1 signaling pathway: Il6, Nos2, Stat3, Hif1a, Vegfa, Cdkn1a, and Rela. Construction of a gene–miRNA network predicted miRNAs targeting four key ferroptosis DEGs—Stat3, Hif1a, Vegfa, and Rela, and further western blot analysis confirmed their upregulation after LSRA, which can be alleviated by DFO pretreatment. Conclusion: The data revealed the induction of ferroptosis in a rat LSRA model and identified possible regulatory roles for ferroptosis-related genes in the molecular mechanisms of LSRA, which provides new insights into the pathogenesis and helps to find new molecular targets for the treatment of LSRA.