Baseline regulatory programs in larval and adult neural progenitors converge towards an injury-induced state after spinal cord injury

Abstract

Regeneration after spinal cord injury requires progenitor cells to convert injury-associated signals into coordinated remodeling of gene regulatory programs. Mammalian spinal progenitors show limited neurogenic output after injury, whereas zebrafish regenerate spinal neurons and recover motor function. To investigate the regulatory changes that allow ependymo-radial glia (ERG) cells, the progenitor cells of the zebrafish spinal cord, to generate new neurons, we combined single-nucleus gene expression and chromatin accessibility profiling across embryonic, larval, and adult stages with topic-based gene regulatory network (GRN) inference. We found that larval and adult ERGs enter the injury response from distinct regulatory baselines: larval progenitors are characterized by a gliogenic program, whereas adult progenitors maintain a comparatively quiescent state. Following injury, both populations gradually change their baseline programs and shift towards a lesion-associated module marked by stress-responsive and chromatin-associated regulators, including jun, hmga1a, hmga2, ybx1, and foxj1a. The shift away from homeostatic states is supported by decreased expression of the Notch-associated regulators nuclear factor I A (nfia) and hey1 in larvae, while in adults, downregulation of the same nuclear factor and other TFs such as bhlhe41 is associated with quiescence exit. Pathway analysis showed stage-specific alterations after injury, characterized predominantly by extracellular signaling and cytoskeletal reorganization in larvae and by metabolic and translational remodeling in adults. Despite divergence from the homeostatic states, injury-induced larval and adult GRNs remain distinct from embryonic hERG regulatory programs. Thus, larval and adult progenitors follow different trajectories from their baselines towards a related lesion-reactive state, in which shared regeneration-associated features are acquired within respective contexts.