A conjoint analysis of bulk RNA-seq and single-nucleus RNA-seq for revealing the role of ferroptosis and iron metabolism in ALS

Abstract

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by progressive and selective degeneration of motor neurons in the motor cortex of brain and spinal cord. Ferroptosis is a newly discovered form of cell death and reported to mediate selective motor neuron death in the mouse model of ALS. The growing awareness of ferroptosis and iron metabolism dysfunction in ALS prompted us to investigate the expression pattern of ferroptosis and iron metabolism-related genes (FIRGs) in ALS. Here, we performed a conjoint analysis of bulk-RNA sequence and single-nucleus RNA sequence data using the datasets from Gene Expression Omnibus (GEO) to reveal the role of FIRGs in ALS, especially in selective motor neuron death of ALS. We first investigated the differentially expressed genes (DEGs) between ALS and non-neurological controls. Weighted gene co-expression network analysis constructed the gene co-expression network and identified three modules closely associated with ALS. Fifteen FIRGs was identified as target genes based on least absolute shrinkage and selection operator regression analysis as follows: ACSL4, ANO6, ATP6V0E1, B2M, CD44, CHMP5, CYBB, CYBRD1, HIF1A, MOSPD1, NCF2, SDCBP, STEAP2, TMEM14C, ULK1. These genes could differentiate ALS patients from non-neurological controls (p < 2.2e−16) and had a valid value in predicting and diagnosing ALS (AUC = 0.881 in primary dataset and AUC = 0.768 in validation dataset). Then we performed the functional enrichment analysis of DEGs between ALS cases, the most significantly influenced by target genes, and non-neurological controls. The result indicated that the most significantly influenced functions in ALS pathogenesis by these identified FIRGs are synapse pathways, calcium signaling pathway, cAMP signaling pathway, and phagosome and several immune pathways. At last, the analysis of single- nuclear seq found that CHMP5, one of the 15 FIRGs identified by bulk single-nucleus RNA-seq data, was expressed significantly higher in ALS than pathologically normal (PN), specifically in excitatory neuron populations with layer 2 and layer 3 markers (Ex L2_L3), layer 3 and layer 5 markers (Ex L3_L5). Taken together, our study indicates the positive correlation between FIRGs and ALS, presents potential markers for ALS diagnosis and provides new research directions of CHMP5 function in selective motor neuron death in ALS.