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.

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.

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

Background

Alterations in neurotransmitter phenotypes of specific neurons can cause imbalances in excitation and inhibition in the central nervous system (CNS), leading to diseases. Therefore, the correct specification and maintenance of neurotransmitter phenotypes is vital. As with other neuronal properties, neurotransmitter phenotypes are often specified and maintained by particular transcription factors. However, the specific molecular mechanisms and transcription factors that regulate neurotransmitter phenotypes remain largely unknown.

Methods

In this paper we use single mutant, double mutant and transgenic zebrafish embryos to elucidate the functions of Lmx1ba and Lmx1bb in the regulation of spinal cord interneuron neurotransmitter phenotypes.

Results

We demonstrate that lmx1ba and lmx1bb are both expressed in zebrafish spinal cord and that lmx1bb is expressed by both V0v cells and dI5 cells. Our functional analyses demonstrate that these transcription factors are not required for neurotransmitter fate specification at early stages of development, but that in embryos with at least two lmx1ba and/or lmx1bb mutant alleles there is a reduced number of excitatory (glutamatergic) spinal interneurons at later stages of development. In contrast, there is no change in the numbers of V0v or dI5 cells. These data suggest that lmx1b-expressing spinal neurons still form normally, but at least a subset of them lose, or do not form, their normal excitatory fates. As the reduction in glutamatergic cells is only seen at later stages of development, Lmx1b is probably required either for the maintenance of glutamatergic fates or to specify glutamatergic phenotypes of a subset of later forming neurons. Using double labeling experiments, we also show that at least some of the cells that lose their normal glutamatergic phenotype are V0v cells. Finally, we also establish that Evx1 and Evx2, two transcription factors that are required for V0v cells to acquire their excitatory neurotransmitter phenotype, are also required for lmx1ba and lmx1bb expression in these cells, suggesting that Lmx1ba and Lmx1bb act downstream of Evx1 and Evx2 in V0v cells.

Conclusions

Lmx1ba and Lmx1bb function at least partially redundantly in the spinal cord and three functional lmx1b alleles are required in zebrafish for correct numbers of excitatory spinal interneurons at later developmental stages. Taken together, our data significantly enhance our understanding of how spinal cord neurotransmitter fates are regulated.

Evx1 and Evx2 specify excitatory neurotransmitter fates and suppress inhibitory fates through a Pax2-independent mechanism

BACKGROUND

For neurons to function correctly in neuronal circuitry they must utilize appropriate neurotransmitters. However, even though neurotransmitter specificity is one of the most important and defining properties of a neuron we still do not fully understand how neurotransmitter fates are specified during development. Most neuronal properties are determined by the transcription factors that neurons express as they start to differentiate. While we know a few transcription factors that specify the neurotransmitter fates of particular neurons, there are still many spinal neurons for which the transcription factors specifying this critical phenotype are unknown. Strikingly, all of the transcription factors that have been identified so far as specifying inhibitory fates in the spinal cord act through Pax2. Even Tlx1 and Tlx3, which specify the excitatory fates of dI3 and dI5 spinal neurons work at least in part by down-regulating Pax2.

METHODS

In this paper we use single and double mutant zebrafish embryos to identify the spinal cord functions of Evx1 and Evx2.

RESULTS

We demonstrate that Evx1 and Evx2 are expressed by spinal cord V0v cells and we show that these cells develop into excitatory (glutamatergic) Commissural Ascending (CoSA) interneurons. In the absence of both Evx1 and Evx2, V0v cells still form and develop a CoSA morphology. However, they lose their excitatory fate and instead express markers of a glycinergic fate. Interestingly, they do not express Pax2, suggesting that they are acquiring their inhibitory fate through a novel Pax2-independent mechanism.

CONCLUSIONS

Evx1 and Evx2 are required, partially redundantly, for spinal cord V0v cells to become excitatory (glutamatergic) interneurons. These results significantly increase our understanding of the mechanisms of neuronal specification and the genetic networks involved in these processes.

Evx1 is required for joint formation in zebrafish fin dermoskeleton

The transcription factor Evx1 is expressed in the joints between individual lepidotrichia (bony ray) segments and at the distal tips of the lepidotrichia in developing zebrafish fins. It is also expressed in the apical growth zone in regenerating fins. However, so far there is no functional evidence that addresses whether Evx1 is required for any aspect of fin development or regeneration. In this study, we use a novel mutation in evx1 to address this. We find that Evx1 is not required for either fin outgrowth or regeneration. All of the fins form normally in evx1 mutants, and there are no significant changes in fin length. In contrast, Evx1 is required for lepidotrichia joint formation during both fin development and regeneration. This is a very specific phenotype as both lepidotrichia hemisegment separations and lepidotrichia bifurcations still form normally in evx1 mutant fins, as do joints in the more proximal endoskeletal radials.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Analysis of HFE-codon 63/282 (H63D/C282Y) gene variants in mexican mestizos. Blood donors and patients with hereditary hemochromatosis

BACKGROUND

The prevalence of hereditary hemochromatosis (HH) (H63D/C282Y) gene variants in Mexico is unknown.

METHODS

Using amplification refractory mutation system polymerase chain reaction, an analysis of HFE-codon 63/282 (H63D/C282Y) gene variants was performed in a group of 153 Mexican mestizo blood donors and six individuals with familial iron overload.

RESULTS

In normal blood donors, three heterozygotes for the C282Y mutation (2.0%) were found, whereas 18 heterozygotes and one homozygote for the H63D mutation (11.8% and 0.6%, respectively) were identified; there was one compound heterozygote for the C282Y/H63D mutation. These data resulted in allele frequencies of 0.013 (+/-0. 2%, alpha = 0.05) and 0.062 (+/-0.9%, alpha = 0.05), respectively, for these two mutations, results similar to those found in whites. In the six patients with the HH phenotype, two were found to be heterozygous for C282Y and one heterozygous for H63D; three individuals with HH had no gene mutations. Two heterozygous HH individuals were found to have iron overload associated with other conditions: one heterozygous for C282Y infected with HIV, and another heterozygous for H63D with heterozygous beta-thalassemia.

CONCLUSIONS

The prevalence of C282Y and H63D HFE gene mutations in Mexican mestizos is similar to that found in other populations. In addition, other gene mutations responsible for HH in the Mexican mestizo population should be investigated, because, in three of six individuals with the HH phenotype, neither of the two mutations was recorded.