Molecular Analyses of V0v Spinal Interneurons and Identification of Transcriptional Regulators Downstream of Evx1 and Evx2 in These Cells. [scRNA-Seq]
Background: V0v spinal interneurons are highly-conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells. Methods: To identify candidate members of V0v gene regulatory networks, we FAC-sorted WT and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and scRNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 mutants and wild-type siblings.  Results: Our data reveal two molecularly-distinct subtypes of V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuronal expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression. Conclusions: This study identifies two molecularly-distinct subsets of V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of V0v interneurons. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type. Overall design: Two samples were analysed in total, both at 48 hours post fertilization (hpf). These two samples represent two biological replicates of all V0v spinal interneurons (isolated from embryos from an incross of heterozygous evx1i232/+;evx2sa140/+;Tg(evx1:EGFP)SU2 parents.