Lmx1b is required for the glutamatergic fates of a subset of spinal cord neurons

Transcriptional Regulators with Broad Expression in the Zebrafish Spinal Cord

Background

The spinal cord is a crucial part of the vertebrate CNS, controlling movements and receiving and processing sensory information from the trunk and limbs. However, there is much we do not know about how this essential organ develops. Here, we describe expression of 21 transcription factors and one transcriptional regulator in zebrafish spinal cord.

Results

We analyzed the expression of aurkb, foxb1a, foxb1b, her8a, homeza, ivns1abpb, mybl2b, myt1a, nr2f1b, onecut1, sall1a, sall3a, sall3b, sall4, sox2, sox19b, sp8b, tsc22d1, wdhd1, zfhx3b, znf804a, and znf1032 in wild-type and MIB E3 ubiquitin protein ligase 1 zebrafish embryos. While all of these genes are broadly expressed in spinal cord, they have distinct expression patterns from one another. Some are predominantly expressed in progenitor domains, and others in subsets of post-mitotic cells. Given the conservation of spinal cord development, and the transcription factors and transcriptional regulators that orchestrate it, we expect that these genes will have similar spinal cord expression patterns in other vertebrates, including mammals and humans.

Conclusions

Our data identify 22 different transcriptional regulators that are strong candidates for playing different roles in spinal cord development. For several of these genes, this is the first published description of their spinal cord expression.

Excitatory Spinal Lhx9-Derived Interneurons Modulate Locomotor Frequency in Mice

Locomotion allows us to move and interact with our surroundings. Spinal networks that control locomotion produce rhythm and left-right and flexor-extensor coordination. Several glutamatergic populations, Shox2 non-V2a, Hb9-derived interneurons, and, recently, spinocerebellar neurons have been proposed to be involved in the mouse rhythm generating networks. These cells make up only a smaller fraction of the excitatory cells in the ventral spinal cord. Here, we set out to identify additional populations of excitatory spinal neurons that may be involved in rhythm generation or other functions in the locomotor network. We use RNA sequencing from glutamatergic, non-glutamatergic, and Shox2 cells in the neonatal mice from both sexes followed by differential gene expression analyses. These analyses identified transcription factors that are highly expressed by glutamatergic spinal neurons and differentially expressed between Shox2 neurons and glutamatergic neurons. From this latter category, we identified the Lhx9-derived neurons as having a restricted spinal expression pattern with no Shox2 neuron overlap. They are purely glutamatergic and ipsilaterally projecting. Ablation of the glutamatergic transmission or acute inactivation of the neuronal activity of Lhx9-derived neurons leads to a decrease in the frequency of locomotor-like activity without change in coordination pattern. Optogenetic activation of Lhx9-derived neurons promotes locomotor-like activity and modulates the frequency of the locomotor activity. Calcium activities of Lhx9-derived neurons show strong left-right out-of-phase rhythmicity during locomotor-like activity. Our study identifies a distinct population of spinal excitatory neurons that regulates the frequency of locomotor output with a suggested role in rhythm-generation in the mouse alongside other spinal populations.

Molecular analyses of zebrafish V0v spinal interneurons and identification of transcriptional regulators downstream of Evx1 and Evx2 in these cells

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 wild-type and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

RESULTS

Our data reveal two molecularly distinct subtypes of zebrafish 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 interneuron 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 zebrafish 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 these cells. 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 and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

Molecular Analyses of V0v Spinal Interneurons and Identification of Transcriptional Regulators Downstream of Evx1 and Evx2 in these Cells

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 single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double 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 these cells. 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 and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

sox1a:eGFP transgenic line and single-cell transcriptomics reveal the origin of zebrafish intraspinal serotonergic neurons

Sox transcription factors are crucial for vertebrate nervous system development. In zebrafish embryo, sox1 genes are expressed in neural progenitor cells and neurons of ventral spinal cord. Our recent study revealed that the loss of sox1a and sox1b function results in a significant decline of V2 subtype neurons (V2s). Using single-cell RNA sequencing, we analyzed the transcriptome of sox1a lineage progenitors and neurons in the zebrafish spinal cord at four time points during embryonic development, employing the Tg(sox1a:eGFP) line. In addition to previously characterized sox1a-expressing neurons, we discovered the expression of sox1a in late-developing intraspinal serotonergic neurons (ISNs). Developmental trajectory analysis suggests that ISNs arise from lateral floor plate (LFP) progenitor cells. Pharmacological inhibition of the Notch signaling pathway revealed its role in negatively regulating LFP progenitor cell differentiation into ISNs. Our findings highlight the zebrafish LFP as a progenitor domain for ISNs, alongside known Kolmer-Agduhr (KA) and V3 interneurons.

Evolution of lbx spinal cord expression and function

Ladybird homeobox (Lbx) transcription factors have crucial functions in muscle and nervous system development in many animals. Amniotes have two Lbx genes, but only Lbx1 is expressed in spinal cord. In contrast, teleosts have three lbx genes and we show here that zebrafish lbx1a, lbx1b, and lbx2 are expressed by distinct spinal cell types, and that lbx1a is expressed in dI4, dI5, and dI6 interneurons, as in amniotes. Our data examining lbx expression in Scyliorhinus canicula and Xenopus tropicalis suggest that the spinal interneuron expression of zebrafish lbx1a is ancestral, whereas lbx1b has acquired a new expression pattern in spinal cord progenitor cells. lbx2 spinal expression was probably acquired in the ray-finned lineage, as this gene is not expressed in the spinal cords of either amniotes or S. canicula. We also show that the spinal function of zebrafish lbx1a is conserved with mouse Lbx1. In zebrafish lbx1a mutants, there is a reduction in the number of inhibitory spinal interneurons and an increase in the number of excitatory spinal interneurons, similar to mouse Lbx1 mutants. Interestingly, the number of inhibitory spinal interneurons is also reduced in lbx1b mutants, although in this case the number of excitatory interneurons is not increased. lbx1a;lbx1b double mutants have a similar spinal interneuron phenotype to lbx1a single mutants. Taken together these data suggest that lbx1b and lbx1a may be required in succession for correct specification of dI4 and dI6 spinal interneurons, although only lbx1a is required for suppression of excitatory fates in these cells.

Induced Torpor as a Countermeasure for Low Dose Radiation Exposure in a Zebrafish Model

The development of the Artemis programme with the goal of returning to the moon is spurring technology advances that will eventually take humans to Mars and herald a new era of interplanetary space travel. However, long-term space travel poses unique challenges including exposure to ionising radiation from galactic cosmic rays and potential solar particle events, exposure to microgravity and specific nutritional challenges arising from earth independent exploration. Ionising radiation is one of the major obstacles facing future space travel as it can generate oxidative stress and directly damage cellular structures such as DNA, in turn causing genomic instability, telomere shortening, extracellular-matrix remodelling and persistent inflammation. In the gastrointestinal tract (GIT) this can lead to leaky gut syndrome, perforations and motility issues, which impact GIT functionality and affect nutritional status. While current countermeasures such as shielding from the spacecraft can attenuate harmful biological effects, they produce harmful secondary particles that contribute to radiation exposure. We hypothesised that induction of a torpor-like state would confer a radioprotective effect given the evidence that hibernation extends survival times in irradiated squirrels compared to active controls. To test this hypothesis, a torpor-like state was induced in zebrafish using melatonin treatment and reduced temperature, and radiation exposure was administered twice over the course of 10 days. The protective effects of induced-torpor were assessed via RNA sequencing and qPCR of mRNA extracted from the GIT. Pathway and network analysis were performed on the transcriptomic data to characterise the genomic signatures in radiation, torpor and torpor + radiation groups. Phenotypic analyses revealed that melatonin and reduced temperature successfully induced a torpor-like state in zebrafish as shown by decreased metabolism and activity levels. Genomic analyses indicated that low dose radiation caused DNA damage and oxidative stress triggering a stress response, including steroidal signalling and changes to metabolism, and cell cycle arrest. Torpor attenuated the stress response through an increase in pro-survival signals, reduced oxidative stress via the oxygen effect and detection and removal of misfolded proteins. This proof-of-concept model provides compelling initial evidence for utilizing an induced torpor-like state as a potential countermeasure for radiation exposure.

Hmx3a Has Essential Functions in Zebrafish Spinal Cord, Ear and Lateral Line Development

Transcription factors that contain a homeodomain DNA-binding domain have crucial functions in most aspects of cellular function and embryonic development in both animals and plants. Hmx proteins are a subfamily of NK homeodomain-containing proteins that have fundamental roles in development of sensory structures such as the eye and the ear. However, Hmx functions in spinal cord development have not been analyzed. Here, we show that zebrafish (Danio rerio) hmx2 and hmx3a are coexpressed in spinal dI2 and V1 interneurons, whereas hmx3b, hmx1, and hmx4 are not expressed in spinal cord. Using mutational analyses, we demonstrate that, in addition to its previously reported role in ear development, hmx3a is required for correct specification of a subset of spinal interneuron neurotransmitter phenotypes, as well as correct lateral line progression and survival to adulthood. Surprisingly, despite similar expression patterns of hmx2 and hmx3a during embryonic development, zebrafish hmx2 mutants are viable and have no obviously abnormal phenotypes in sensory structures or neurons that require hmx3a In addition, embryos homozygous for deletions of both hmx2 and hmx3a have identical phenotypes to severe hmx3a single mutants. However, mutating hmx2 in hypomorphic hmx3a mutants that usually develop normally, results in abnormal ear and lateral line phenotypes. This suggests that while hmx2 cannot compensate for loss of hmx3a, it does function in these developmental processes, although to a much lesser extent than hmx3a More surprisingly, our mutational analyses suggest that Hmx3a may not require its homeodomain DNA-binding domain for its roles in viability or embryonic development.

Small RNA-Seq reveals novel miRNAs shaping the transcriptomic identity of rat brain structures

Small RNA-Seq of the rat central nervous system reveals known and novel miRNAs specifically regulated in brain structures and correlated with the expression of their predicted target genes, suggesting a critical role in the transcriptomic identity of brain structures.

In the central nervous system (CNS), miRNAs are involved in key functions, such as neurogenesis and synaptic plasticity. Moreover, they are essential to define specific transcriptomes in tissues and cells. However, few studies were performed to determine the miRNome of the different structures of the rat CNS, although a major model in neuroscience. Here, we determined by small RNA-Seq, the miRNome of the olfactory bulb, the hippocampus, the cortex, the striatum, and the spinal cord and showed the expression of 365 known miRNAs and 90 novel miRNAs. Differential expression analysis showed that several miRNAs were specifically enriched/depleted in these CNS structures. Transcriptome analysis by mRNA-Seq and correlation based on miRNA target predictions suggest that the specifically enriched/depleted miRNAs have a strong impact on the transcriptomic identity of the CNS structures. Altogether, these results suggest the critical role played by these enriched/depleted miRNAs, in particular the novel miRNAs, in the functional identities of CNS structures.

Tal1, Gata2a, and Gata3 Have Distinct Functions in the Development of V2b and Cerebrospinal Fluid-Contacting KA Spinal Neurons

Vertebrate locomotor circuitry contains distinct classes of ventral spinal cord neurons which each have particular functional properties. While we know some of the genes expressed by each of these cell types, we do not yet know how several of these neurons are specified. Here, we investigate the functions of Tal1, Gata2a, and Gata3 transcription factors in the development of two of these populations of neurons with important roles in locomotor circuitry: V2b neurons and cerebrospinal fluid-contacting Kolmer-Agduhr (KA) neurons (also called CSF-cNs). Our data provide the first demonstration, in any vertebrate, that Tal1 and Gata3 are required for correct development of KA and V2b neurons, respectively. We also uncover differences in the genetic regulation of V2b cell development in zebrafish compared to mouse. In addition, we demonstrate that Sox1a and Sox1b are expressed by KA and V2b neurons in zebrafish, which differs from mouse, where Sox1 is expressed by V2c neurons. KA neurons can be divided into ventral KA″ neurons and more dorsal KA′ neurons. Consistent with previous morpholino experiments, our mutant data suggest that Tal1 and Gata3 are required in KA′ but not KA″ cells, whereas Gata2a is required in KA″ but not KA′ cells, even though both of these cell types co-express all three of these transcription factors. In gata2a mutants, cells in the KA″ region of the spinal cord lose expression of most KA″ genes and there is an increase in the number of cells expressing V3 genes, suggesting that Gata2a is required to specify KA″ and repress V3 fates in cells that normally develop into KA″ neurons. On the other hand, our data suggest that Gata3 and Tal1 are both required for KA′ neurons to differentiate from progenitor cells. In the KA′ region of these mutants, cells no longer express KA′ markers and there is an increase in the number of mitotically-active cells. Finally, our data demonstrate that all three of these transcription factors are required for later stages of V2b neuron differentiation and that Gata2a and Tal1 have different functions in V2b development in zebrafish than in mouse.

Zebrafish transgenic constructs label specific neurons in Xenopus laevis spinal cord and identify frog V0v spinal neurons

A correctly functioning spinal cord is crucial for locomotion and communication between body and brain but there are fundamental gaps in our knowledge of how spinal neuronal circuitry is established and functions. To understand the genetic program that regulates specification and functions of this circuitry, we need to connect neuronal molecular phenotypes with physiological analyses. Studies using Xenopus laevis tadpoles have increased our understanding of spinal cord neuronal physiology and function, particularly in locomotor circuitry. However, the X. laevis tetraploid genome and long generation time make it difficult to investigate how neurons are specified. The opacity of X. laevis embryos also makes it hard to connect functional classes of neurons and the genes that they express. We demonstrate here that Tol2 transgenic constructs using zebrafish enhancers that drive expression in specific zebrafish spinal neurons label equivalent neurons in X. laevis and that the incorporation of a Gal4:UAS amplification cassette enables cells to be observed in live X. laevis tadpoles. This technique should enable the molecular phenotypes, morphologies and physiologies of distinct X. laevis spinal neurons to be examined together in vivo. We have used an islet1 enhancer to label Rohon-Beard sensory neurons and evx enhancers to identify V0v neurons, for the first time, in X. laevis spinal cord. Our work demonstrates the homology of spinal cord circuitry in zebrafish and X. laevis, suggesting that future work could combine their relative strengths to elucidate a more complete picture of how vertebrate spinal cord neurons are specified, and function to generate behavior. © 2017 Wiley Periodicals, Inc. Develop Neurobiol 77: 1007-1020, 2017.