Research article | Open | Published:
Evx1 and Evx2 specify excitatory neurotransmitter fates and suppress inhibitory fates through a Pax2-independent mechanism
Neural Developmentvolume 11, Article number: 5 (2016)
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.
In this paper we use single and double mutant zebrafish embryos to identify the spinal cord functions of Evx1 and Evx2.
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.
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.
Hundreds of millions of people across the world are affected by neurological diseases and injuries. Understanding how functional neuronal circuitry is established in the vertebrate central nervous system (CNS) is essential for developing better treatments for these conditions. How neuronal circuitry develops is also a fundamental question in developmental neuroscience. To answer this question, we need to identify how the functional properties of distinct neurons are specified; since these properties determine which circuits the neurons participate in, the functional roles that the neurons have within those circuits and the resulting outputs of the circuitry. The spinal cord is a powerful system for establishing fundamental principles of neuronal fate specification, function and circuit assembly, as it is relatively simple and experimentally tractable compared to the brain. This has enabled considerable progress in establishing the functions of different ventral spinal cord interneurons in locomotor circuitry (e.g. [1–9]). However, we still know relatively little about how the functional properties of these cells are determined.
For neurons to function correctly they must synthesize and utilize correct neurotransmitters. Within neuronal circuitry, if they inhibit rather than excite their synaptic partners, or vice versa, then the behaviors and functional outputs of those circuits will be dramatically disturbed, and may give rise to pathological conditions. For example, disruptions in the balance of excitatory and inhibitory neurons in the CNS have been implicated in epilepsy, autism, Alzheimer’s and many other neurological disorders (e.g. [10–13]). 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.
Many neuronal properties are determined by the transcription factors that cells express as they start to differentiate. We already know a few transcription factors (e.g. Ptf1a, Lhx1, Lhx5, Lbx1, Pax2) that specify the inhibitory (GABAergic and/or glycinergic) fates of several subsets of spinal interneurons [14–18]. Strikingly, most of these transcription factors function in dorsal spinal neurons and all of them act through Pax2 [14–21]. In contrast, we only know two transcription factors, Tlx1 and Tlx3, that are required for the specification of excitatory (glutamatergic) fates and these are only expressed in dorsal dI3, dI5 and DILB cells [15, 16, 22]. Interestingly, Tlx1 and Tlx3 determine the glutamatergic fates of dI3 and dI5 cells at least in part by down-regulating Pax2 . These results suggest that Pax2 is a crucial player in neurotransmitter fate specification with its presence being required for inhibitory fates and its absence required for excitatory fates. However, we still do not know which transcription factors regulate the neurotransmitter fates of many excitatory spinal neurons, including those in the ventral spinal cord, whose correct functional specification is essential for locomotion.
In this paper we identify two transcription factors, Evx1 and Evx2, which are required for a subset of excitatory fates in the ventral spinal cord. In mammals, the spinal cord expression of Evx1 and Evx2 is restricted to a population of cells located in an intermediate dorso-ventral position corresponding to V0 cells (e.g. [23–28]). V0 cells are post-mitotic cells that form from the p0 (Dbx1-positive, Nkx6.2-negative) progenitor domain [23, 27–29]. These cells develop into interneurons that are important components of locomotor circuitry and they can be subdivided into an Evx1-positive sub-population called V0v cells and an Evx1-negative sub-population called V0D cells. These names reflect the fact that V0v cells form more ventrally than V0D cells (e.g. [23–28, 30–34]). Evx2 is expressed in a similar pattern to Evx1 in the mouse CNS, suggesting that it may also be expressed by V0v cells. This is consistent with the observation that Evx2 spinal cord expression is lost in mouse Evx1 mutants . However, co-expression of Evx1 and Evx2 in the mouse spinal cord has not yet been demonstrated .
In mammals, both V0v and V0D interneurons are crucial for correct left-right alternation during locomotion, with V0v cells in particular being required for hindlimb left-right alternation during fast locomotion [9, 34]. While the functions of V0 cells in specific behaviors have so far only been assayed in mouse, these cells have highly conserved commissural axon trajectories in all animals examined so far ([23–28, 32, 33, 35, 36]; this paper), suggesting that their functional properties are likely to be highly conserved across the vertebrate lineage. However, when we started this work the neurotransmitter phenotype of V0v cells had not been identified.
In zebrafish, evx1 and evx2 are expressed in a similar intermediate dorsal-ventral spinal cord position to that observed in other vertebrates [26, 32, 33], although again, co-expression of these two genes has not previously been demonstrated. In this paper, we confirm that evx1 and evx2 are co-expressed by V0v cells and we show that V0v cells are glutamatergic and have a Commissural Ascending (Comissural Secondary Ascending or CoSA) morphology. We also provide the first analysis of evx1;evx2 double mutants in any vertebrate and the first analysis of the spinal cord phenotype of evx2 mutants. Significantly, we demonstrate that Evx1 and Evx2 are required, partially redundantly, to specify the glutamatergic fates of V0v cells. Given that we know so little about how excitatory fates are specified in the spinal cord and particularly the ventral spinal cord, these findings add considerably to our understanding of CNS circuit development.
In the absence of both Evx1 and Evx2, V0v cells lose their glutamatergic fates but other functional characteristics like soma/cell body morphology and axon trajectory are unchanged. In addition, and in contrast to a previously described mouse Evx1 mutant , these cells do not express markers of neighboring cell types. This suggests that V0v cells are not transfating into a different class of neuron; they have just changed some of their functional properties. Strikingly, in evx1;evx2 double mutants V0v cells become inhibitory, but they do not express Pax2, suggesting that they are acquiring their inhibitory fates through a novel Pax2-independent mechanism.
All zebrafish experiments in this research were approved either by the UK Home Office or by the Syracuse University IACUC committee.
Zebrafish husbandry and fish lines
Zebrafish (Danio rerio) were maintained on a 14 h light/10 h dark cycle at 28.5 °C and embryos were obtained from natural, paired and/or grouped spawnings of wild-type (WT) adults (AB, TL or AB/TL hybrids), identified heterozygous or homozygous Tg(slc17a6:EGFP) (used to be called Tg(vGlut2a:EGFP); [36, 37]), Tg(evx1:EGFP) SU1 or Tg(evx1:EGFP) SU2 adults, double heterozygous evx1 i232 ;evx2 sa140 mutants or double heterozygous evx1 i232 ;evx2 sa140 mutants that also carried one of the Tg(evx1:EGFP) lines (see below). Embryos were reared at 28.5 °C and staged by hours post fertilization (h), days post or prim staging/or prim staging .
The evx1 i232 mutation has been described before . The evx2 sa140 mutant was received from the Wellcome Trust Sanger Centre, (https://www.sanger.ac.uk/sanger/Zebrafish_Zmpbrowse). Both mutations produce a single base pair change (a C to a T in the case of evx2 sa140) that results in a premature stop codon before the homeobox (; Fig. 1a). Therefore, the truncated proteins, if formed, will lack DNA binding domains. The evx2 sa140 mutation creates a BfaI recognition site that enables us to genotype individual fish and embryos (see below; Fig. 1b). In this paper we demonstrate that the evx2 sa140 allele does not make Evx2 protein (Fig. 1c), strongly suggesting that it is a null allele. The evx1 i232 mutant is also probably a null allele . However, in contrast to the evx1 i232 mutant which is viable, the evx2 sa140 mutant is embryonic lethal (see Results and Additional file 1: Results).
Genotyping of mutant alleles was performed on both live adults and fixed embryos using DNA extracted from fin clips and dissected heads respectively. Fin clipping and evx1 genotyping were performed as in . To extract DNA from embryos, heads were removed in 80 % glycerol / 20 % PBS with insect pins. Embryonic trunks were stored in 70 % glycerol / 30 % PBS at 4 °C for later analysis. Heads were incubated in 50 μL of Proteinase K buffer solution (1 M Tris-HCl, pH 8.2; 0.5 M EDTA; 1 M NaCl; 20 % SDS; 10 mg/ml Proteinase K) for 2 h at 55 °C. Proteinase K was heat inactivated at 100 °C for 10 min and tubes were centrifuged for 20 min at 13,000 rpm. DNA was precipitated with 100 % ethanol at -20 °C overnight, centrifuged to pellet the DNA and re-suspended in 20 μL of water. 2 μL of DNA was used for each PCR.
The evx2 sa140 mutation creates a BfaI recognition site. A genomic region flanking the mutation site was PCR amplified using the following conditions: 94 °C for 60 s, followed by 5 cycles of 92 °C for 30 s, 54 °C for 30 s, 72 °C for 60 s; followed by 40 cycles of 92 °C for 20 s, 52 °C for 30 s, and 72 °C for 60 s, followed by a final extension at 72 °C for 5 min. Forward primer: GTAATGCGATCCCAAAACG. Reverse primer: TTATTTTAGATTTGGCAATGG.
PCR products were digested with BfaI and analysed on a 1 % agarose gel. The WT product is 454 bp, whereas the mutant product is cut into 218 bp and 236 bp fragments. These fragments are close enough in size that they are usually detected as one band on an agarose electrophoresis gel (Fig. 1b).
Creation of Tg(evx1:EGFP) lines
Potential evx1 enhancer regions were identified by multispecies sequence comparisons using the global alignment program Shuffle-LAGAN  and visualized using VISTA . Zebrafish (Danio rerio) evx1 genomic sequence (ENSDARG00000099365) and orthologous sequences from human (ENSG00000106038) and mouse (ENSMUSG00000005503) were obtained from Ensembl (http://www.ensembl.org). Zebrafish sequence was used as the baseline and annotated using exon/intron information from Ensembl. The alignment was performed using a 100 bp window and a cutoff score of 70 % identity. A multi-species comparison of approximately 23Kb of Danio rerio genomic DNA sequence containing evx1 and extending 5 Kb into flanking regions revealed high conservation in both coding and non-coding sequences among compared species. We identified three Conserved Non-coding Elements (CNEs) located 5’ and 3’ to evx1. The first is located 79 bp upstream of zebrafish evx1 coding sequence and extends over 100 bp. The other two are located 3' to evx1. One is 2354 bp downstream of the stop codon and is 184 bp long whereas the other is 2979 bp downstream of the stop codon and extends over 140 bp (Fig. 1d). One amplicon encompassing these two 3’ CNEs and the intervening sequence (1.34 Kb) was PCR-amplified from genomic DNA. Forward primer: AAGATTGGAATGGAATGTCT. Reverse primer: GCATTTTCGCCTTTGCATCA.
The Tg(evx1:EGFP) SU1 line was generated by cloning this 3’ CNE amplicon into the pDONRTMP4-P1R vector from Invitrogen using Gateway technology [42, 43]. The final construct was assembled using the pENTRbasegfp plasmid and the pCSDest2 vector . This resulted in a vector containing Tol2:1.3Kb 3’ zfish evx1:ßcarp minimal promoter:EGFP:Tol2.
The same 3’ CNE amplicon was used to generate the Tg(evx1:EGFP) SU2 line. In this case, the GAL4VP16;UAS-EGFP cassette was taken from the pBGAL4VP16;UAS-EGFP plasmid  and cloned into a middle entry vector from Invitrogen [42, 43]. An oligonucleotide containing the cfos minimal promoter sequence  plus 17 bp of the 5’ arm of GAL4 was synthesized and used with a RVeGFPAttb2 primer to PCR amplify the GAL4VP16;UAS-EGFP cassette.
FWcfosGAL4VP16 primer: CACTCATTCATAAAACGCTTGTTATAAAAGCAGTGGCTGCGGCGCCTCGTACTCCAACCGCATCTGCAGCGAGCAACTGAGAAGCCAAGACTGAGCCGGCGGCCTTTGTACAAAAAAGCAG
RVeGFPAttb2 primer: GGGGACCACTTTGTACAAGAAAGCTGGGTTTACTTGTACAGCTCGTCCA.
This PCR product was used to generate a second PCR product using the primer FWattB1cfos: GGGGACAAGTTTGTACAAAAAAGCAGGCTCACTCATTCATAAAATCGCTT and the RVeGFPAttb2 primer.
The final amplicon was cloned into pDONR™221 using gateway technology. This middle entry vector was used to generate a final vector containing Tol2:1.3Kb 3’ zfish evx1:cfos minimal promoter:GAL4VP16;UAS-EGFP:Tol2.
Each of these plasmids was separately co-injected with transposase mRNA into 1-2 cell embryos as described by . Embryos were raised to adulthood and out-crossed to identify founders. In each case one stable transgenic line was generated. The Tg(evx1:EGFP) SU2 or Tg(1.3 kb evx1:cfos:GAL4-UAS:EGFP) line has the advantage that it contains a GAL4-UAS cassette to amplify EGFP expression. This facilitates visualization of axons. However, this line has a slightly more variegated expression than the Tg(evx1:EGFP) SU1 or Tg(1.3Kb evx1:ßcarp:EGFP) line, presumably because of stochastic silencing of the construct due to the GAL4-UAS sequences . In contrast the Tg(evx1:EGFP) SU1 line labels all V0v cells more consistently, but the EGFP expression is slightly weaker and we were never able to obtain evx1 -/- ;evx2 -/- double mutant embryos that contained this transgene, even though we could obtain evx1 -/- ;evx2 +/+ and evx1 -/- ;evx2 +/- embryos. This suggests that the Tg(1.3Kb evx1:ßcarp:EGFP) construct integrated close to the WT evx2 allele.
Approximately 5 nl of a 1:1 combination of two Evx2 ATG Morpholino antisense oligonucleotides (MOs) at 1.25 mg/ml each were injected into 1-2 cell embryos (evx2-1 MO: TTCTTTTCTTATCCTCTCCATCATG; evx2-2 MO: AATCCAAAGTCCCAGGGCTGGTGCT). In all cases, we confirmed that the MOs had completely knocked down Evx2 using immunohistochemistry for zebrafish Evx2.
Expression profiling V0v cells
To determine which neurotransmitters V0v cells express and to identify additional transcription factors expressed by these cells, different combinations of spinal cord and trunk cells were extracted from live transgenic zebrafish embryos at 27 h using fluorescence activated cell-sorting (FACS). Prior to FACS, embryos were prim-staged, deyolked, dissected and dissociated as in [49, 50]. In all cases, the heads were removed to ensure that only trunk or spinal cord cells were collected. Pure populations of cells were obtained using combinations of the following transgenic lines: Tg(elav13:EGFP), Tg(evx1:EGFP) SU1, Tg(pax2a:GFP), Tg(Xla.Tubb:DsRed (formerly Tg(NBT:DsRed)), Tg(vsx2:DsRed) and Tg(gata1:GFP) [8, 51–54]. Trunk samples correspond to FAC-sorted trunk cells (spinal cord and other tissues). All neuron samples are EGFP-positive cells from Tg(elav13:EGFP) trunks. V0v neurons are EGFP-positive cells from Tg(evx1:EGFP) SU1 trunks. V1 neurons are double-positive EGFP-positive, DsRed-positive cells from Tg(pax2a:GFP);Tg(Xla.Tubb:DsRed) trunks. V2a neurons are double-positive DsRed-positive, EGFP-positive cells from Tg(vsx2:DsRed);Tg(elavl3:EGFP) trunks. V2b + KA neurons are double-positive EGFP-positive, DsRed-positive cells from Tg(gata1:GFP);Tg(Xla.Tubb:DsRed) trunks. Total RNA was extracted using an RNeasy Micro Kit (Qiagen, 74004). RNA quality and quantity was assayed on an Agilent 2100 Bioanalyser (RNA 6000 Pico Kit, Agilent, 5067-1513), before converting to fluorescently-labelled cDNA (Ovation Pico WTA System V2, Pico, 3302) and hybridizing to a custom-designed Agilent microarray (EMBL Genomics Core, Heidelberg). Details of this microarray will be described elsewhere, along with the characterization of additional genes identified from these analyses. Data pre-processing and normalization was performed using Bioconductor software (https://www.bioconductor.org/). Two-class eBayes and three-class ANOVA analyses were performed using GEPAS software (Tárraga, (2008)). All reported statistics were corrected for multiple testing (Benjamini and Hochberg (1995)).
in situ hybridization
Embryos were fixed in 4 % paraformaldehyde and single and double in situ hybridizations were performed as previously described [55, 56]. RNA probes were prepared using the following templates, dbx1a and dbx1b , evx1 , evx2 , eve1 , pax2a, pax2b, pax8  and eng1b . To determine neurotransmitter phenotypes we used probes for genes that encode proteins that transport or synthesize specific neurotransmitters. A mixture of two probes (glyt2a and glyt2b) for slc6a5 (previously called glyt2) was used to label glycinergic cells [61, 62]. slc6a5 encodes for a glycine transporter necessary for glycine reuptake and transport across the plasma membrane. A mixture of two probes to gad1b (previously called gad67, probes used to be called gad67a and gad67b) and one probe to gad2 (previously called gad65) was used to label GABAergic cells [61, 62]. gad1b and gad2 encode for glutamic acid decarboxylases, necessary for the synthesis of GABA from glutamate. A mixture of slc17a6b (formerly called vglut2.1) and slc17a6a (formerly called vglut 2.2) probes was used to label glutamatergic cells [61, 62]. These genes encode proteins responsible for transporting glutamate to the synapse. In all of these cases, a mix of equal concentrations of the relevant probes was used [61, 62]. We also used slc32a1 (formerly called viaat), which encodes for a vesicular inhibitory amino acid transporter, to label all inhibitory cells .
The DNA template for the skor2 (ZDB-GENE-060825-57) probe was generated by PCR-amplifying the 3’ region of skor2 from cDNA using a reverse primer containing a T3 promoter sequence at the 5’ end (indicated in italics below). Total RNA was extracted by homogenizing 50-100 mg of 27hpf wild-type zebrafish embryos in 1 mL of TRIzol reagent (Ambion, 15596-026). cDNA was synthesized using Bio-Rad iScript Reverse Transcription Supermix kit (Bio-Rad, 170-8891). A 50 μL PCR was assembled containing 5 μL cDNA and 1 unit of Phusion High-Fidelity DNA Polymerase (NEB, M0530L). PCR conditions were: 94 °C for 3 min followed by 35 cycles of 94 °C for 30 s, 56.5 °C for 30 s, 72 °C for 1.5 min and then a final extension step of 72 °C for 10 min. PCR product was purified by phenol:chloroform extraction. Forward primer: CGCAAGACGCTTTTTATCC
Reverse primer: AATTAACCCTCACTAAAGGGAAAATGGAGAGCTGCCTTTCAG.
ZFIN Identification numbers are provided for all genes in Additional file 1: Table S2.
Primary antibodies used were rabbit anti-Evx2 (a kind gift from Dr Higashijima, described in Satou et al., 2012, raised against the first 168 amino acids of zebrafish Evx2, a region with no significant homology to zebrafish Evx1, 1:300), mouse anti-GFP (Roche Applied Science, 11814460001, 1:500), rabbit anti-GFP (Molecular Probes A6465, 1:500) and mouse anti-Pax2 (Covance PRB-276P 1:300). The Pax2 antibody recognizes both Pax2a and Pax2b in zebrafish . Antibodies used for fluorescent in situ hybridization were mouse anti-Dig (Jackson ImmunoResearch 200-002-156, 1:5000) and rabbit anti-Flu (Invitrogen A889, 1:2500). These were detected with Invitrogen Tyramide kits #12 and #5. Secondary antibodies used were Alexa Fluor 568 goat anti-rabbit (Molecular Probes A11036, 1:500), Alexa Fluor 488 goat anti-rabbit (Molecular Probes A11034, 1:500) and Alexa Fluor 488 goat anti-mouse (Molecular Probes A11029, 1:500).
Embryos for immunohistochemistry were treated with acetone for 15 min (24 h embryos) or 20 min (30 h embryos) to permeabilize them, washed for 5 min in distilled water, then washed 2 x 10 min in PBS. Embryos were treated with Image-iT Signal Enhancer (Invitrogen, I36933) for 30 min, then incubated in block solution (2 % goat serum, 1 % BSA, 10 % DMSO and 0.5 % Triton) for 1 h at room temperature followed by incubation in primary antibody in fresh block solution at 4 °C overnight. Embryos were washed with PBT (PBS + 0.1 % Triton) for 2 h at room temperature and incubated with secondary antibody in block solution at 4 °C overnight. Embryos were then washed with PBT for at 2 h at room temperature and stored in 2 % DABCO (Acros Organics, AC11247-1000).
For 3,3’- diaminobenzidine (DAB) staining, after incubation with primary antibody, samples were incubated in fresh blocking solution with goat anti-rabbit IgG (Covance SMI-5030C, 1:200) at 4 °C overnight. Embryos were then washed with PBT for 2 h and incubated with rabbit PAP (Covance SMI-4010 L, 1:200) in block solution at 4 °C overnight. Embryos were then washed in PBT for 2 h. Staining was performed using Sigma Fast 3,3’- diaminobenzidine tablets (Sigma, D4293).
Embryos were mounted in 70 % glycerol, 30 % PBS and DIC pictures were taken using an AxioCam MRc5 camera mounted on a Zeiss Axio Imager M1 compound microscope. Fluorescent images were taken on a Zeiss LSM 710 confocal microscope. Images were processed using Adobe Photoshop software (Adobe, Inc) and Image J software (Abramoff et al., 2004). In some cases different focal planes were merged to show labeled cells at different medial lateral positions in the spinal cord.
Cell counts and statistics
In all cases, cells counts are for both sides of a five-somite length of the spinal cord adjacent to somites 6-10. Most values are an average of at least 5 embryos. Exceptions are the skor2 + Tg(evx1:EGFP) double-labeling experiments, the skor2 + Tg(slc17a6:EGFP) double-labeling experiments and the pax2a and eng1b in situ hybridization results. In all of these cases 4 embryos were counted. Results were analyzed using the student’s t-test; Error bars indicate standard deviation.
Zebrafish V0v cells express evx1 and evx2
In mouse, Evx1 is expressed in V0v cells and while double-labeling experiments have not yet been performed, the data suggest that Evx2 is probably co-expressed by these same cells [23, 24, 27, 28]. Previous reports described evx1 and evx2 expression in a similar region of zebrafish spinal cord [32, 33] but didn’t determine whether these genes are co-expressed or the specific cell types that express them.
To address these questions, we performed single and double in situ hybridization experiments and found that zebrafish evx1 and evx2 are co-expressed in an intermediate region of the dorso-ventral axis of the spinal cord (Figs. 1g and 2a & i). We further confirmed that evx1 and evx2 are co-expressed in zebrafish spinal cord using an EGFP line, Tg(evx1:EGFP) SU1. We constructed this line using enhancer sequences identified downstream of evx1 (see methods & Fig. 1d). We confirmed that the stable line recapitulates endogenous evx1 expression (Fig. 1f) and also shows co-expression of EGFP and Evx2 protein (Fig. 1e). Interestingly, at 27 h, all of the cells that express Evx2 also express EGFP (and hence evx1) (Fig. 1e), however, a few cells express EGFP but not Evx2. Similarly, at 24 h, a few cells express evx1 but not evx2 (Fig. 1g). This is consistent with earlier reports that suggest that evx1 may be expressed in the spinal cord slightly earlier than evx2 [32, 33], although we cannot rule out the possibility that there is a very small subset of evx1-expressing V0v cells that do not express evx2.
V0 cells develop from the p0 progenitor domain, which expresses dbx1 [28, 63]. Therefore, to confirm that evx1/2-expressing cells are V0v cells we performed EGFP immunohistochemistry and in situ hybridization for dbx1a and dbx1b in Tg(evx1:EGFP) SU1 embryos. We found that zebrafish evx genes are expressed lateral to cells expressing both of these dbx1 genes, as would be predicted for cells developing from the p0 domain. In addition, dbx1a and dbx1b expression persists in some EGFP-positive cells (Fig. 1h), suggesting that these genes continue to be expressed by V0v cells for a short while after they become post-mitotic.
Zebrafish also have a third evx gene, called eve1, but earlier studies suggested that this gene is not expressed in the spinal cord [26, 32, 59, 64]. We confirmed this by examining eve1 expression at multiple stages of spinal cord development (every two somites from 2-somites - 24 h; Fig. 3 and data not shown). In all cases we never saw any spinal cord expression, only expression in the developing tailbud. To check whether eve1 expression is altered in the absence of Evx1 and/or Evx2 we also examined expression in evx1;evx2 double mutants. However, we saw no change in eve1 expression in these double mutants (Fig. 3f). Therefore, this gene is not considered further in this paper. We also confirmed that no additional evx genes exist in zebrafish (Additional file 1: Results and Figure S1).
Zebrafish V0v cells develop into commissural ascending interneurons
In mouse, V0v cells develop commissural axons that ascend (grow rostrally) for one to four somite segments . We used both the Tg(evx1:EGFP) SU1 transgenic line and an additional transgenic line that has stronger expression in V0v cell axons, Tg(evx1:EGFP) SU2 (see methods), to examine the morphology of zebrafish V0v cells. We found that by 27-30 h, almost all of the cells have extended their axons ventrally and have at least started to cross the midline to the other side of the spinal cord (Fig 4a & b). By 48 h, most of the cells have reached the other side of the spinal cord and have turned towards the head, giving them a clear commissural ascending, or CoSA [65, 66], morphology (Fig. 4c & f).
V0v cells are glutamatergic
In wild-type mouse spinal cords, V0 cells develop into both inhibitory (glycinergic or GABAergic) and excitatory (glutamatergic) interneurons [27, 28, 36], but when we started this project the neurotransmitter phenotype of V0v cells had not been established. Zebrafish have both excitatory and inhibitory CoSA interneurons [61, 62], so we could not infer the neurotransmitter properties of V0v cells from their morphology alone. Therefore, we performed double-labeling experiments to establish the neurotransmitter fates of these cells. Double immunohistochemistry for Evx2 and EGFP in the Tg(slc17a6:EGFP) line which labels glutamatergic interneurons in the zebrafish spinal cord [37, 67] revealed that all of the Evx2-positive cells are glutamatergic (Fig. 5a). As expected, not all of the glutamatergic cells express Evx2, since there are several excitatory cell types in the zebrafish spinal cord and only V0v cells express Evx2. In addition, in situ hybridization for slc17a6 (vglut) genes, markers of glutamatergic cells (see methods for more details) combined with EGFP immunohistochemistry in Tg(evx1:EGFP) SU1 embryos also confirmed that evx1-expressing cells are glutamatergic (Fig. 5b). We also examined if evx1 is co-expressed with markers of any other spinal cord neurotransmitter fates. Using double in situ hybridization we found no co-expression between evx1 and markers for glycinergic or GABAergic markers (data not shown). We also observed no double-labeled cells when we performed EGFP immunohistochemistry and in situ hybridization for slc32a1 (formerly called viaat) which labels all inhibitory neurons [8, 68, 69] in Tg(evx1:EGFP) SU2 embryos (Fig. 5c).
Consistent with these analyses, when we FAC-sorted and expression profiled EGFP-labeled V0v cells using the Tg(evx1:EGFP) SU1 line we found that these cells express markers of glutamatergic fates (slc17a6a (formerly called vglut2.2) and slc17a6b (formerly called vglut2.1)) and do not express either glycinergic markers (slc6a9 (formerly called glyt1) or slc6a5 (formerly called glyt2)) or GABAergic markers (gad1b or gad2) (Fig. 5d & e).
Zebrafish V0v cells do not express Pax2
As discussed above, Pax2 is an important regulator of inhibitory spinal cord fates [14–21]. Given that we had shown that V0v cells are excitatory (glutamatergic) and we had not observed any inhibitory V0v cells (Figs 5a-e), we would not predict that V0v cells would express Pax2. However, the literature contains contradictory evidence as to whether evx1 and evx2 are co-expressed with pax2 in the spinal cord [35, 58, 70, 71]. To resolve this issue, we performed double-labeling experiments for evx1 and pax2 using several complementary approaches. These included in situ hybridization for evx1 and immunohistochemistry for Pax2, EGFP immunohistochemistry and in situ hybridization for pax2a, pax2b and pax8 (three highly-related pax genes that are co-expressed in zebrafish spinal cord cells ) in Tg(evx1:EGFP) SU2 embryos and double immunohistochemistry for EGFP and Pax2 in Tg(evx1:EGFP) SU2 embryos. In all cases, we observed no double-labeled cells, suggesting that evx1 and evx2 are not co-expressed with pax2/ pax8 genes in zebrafish spinal cord (Fig. 6g and data not shown). These analyses complement those of Satou and colleagues  who recently reported that inhibitory V0 cells, which presumably correspond to V0D cells, express Pax2, but Evx2-expressing cells do not. Taken together, these data strongly suggest that V0v cells do not express pax2 genes.
evx2 sa140 is a null allele
To identify the functions of evx1 and evx2 in zebrafish V0v cells we used evx1 and evx2 mutants (see Methods). Our previous analyses suggest that evx1 i232 is a null allele . The evx2 sa140 mutation introduces a premature stop codon just before the homeodomain, suggesting that if a truncated protein is synthesized it will have no DNA binding activity. However, it is possible that a truncated protein might retain some function in the embryo. To determine if any Evx2 protein is made in mutant embryos, we used an Evx2 antibody that was made against the first 168 amino acids of zebrafish Evx2 (which corresponds exactly to the region upstream of the premature stop codon in the evx2 sa140 mutant allele). We found that all of the WT (21/51) and heterozygous embryos (20/51) had Evx2 antibody staining but all of the homozygous mutants (10/51) did not (Fig. 1c). This strongly suggests that the evx2 mutant allele does not produce any protein and is a null allele.
evx2 sa140 homozygous mutants are not viable
Unlike the zebrafish evx1 mutant, which is homozygous viable , we never identified an adult evx2 homozygous mutant (n = 262 fish from incrosses of heterozygous evx2 fish, P < 0.0001 using chi-squared test). However, we did obtain fish homozygous for evx1 and heterozygous for evx2 (17 fish identified from a total of 191 fish from incrosses of heterozygous double mutants; P = 0.13 using chi-squared test). Our analyses of incrosses from identified evx2 heterozygous fish suggest that evx2 mutant embryos have no obvious morphological defects for the first few days of development but that most of them die by larval stages (see Additional file 1: Results).
V0v cells still form in evx1;evx2 double mutants
In mouse Evx1 mutants, expression of Evx2 is lost and there is an increase in the number of cells expressing the V1 marker En1. In addition, many of the cells that would normally have expressed Evx1 develop axon trajectories similar to V1 cells. Most strikingly their axons change from being commissural to ipsilateral . This suggests that in the absence of Evx1, most mouse V0v cells transfate to V1 cells.
In contrast, we found that in zebrafish evx1 mutants there is only a small reduction in the number of cells expressing evx2 RNA and Evx2 protein in the spinal cord (approximately a 22 % reduction for evx2 at 24 h and a 32 % reduction for Evx2 at 30 h; Fig. 2c, l & n; Table 1), although expression of evx2 RNA is lost completely in evx1;evx2 double mutants (Fig. 2b & n; Table 1). In contrast, there is no difference in evx1 expression in evx1 or evx2 single mutants when compared to WT siblings, although approximately 47 % of V0v cells lose expression of evx1 in double mutants (Fig. 2j & m; Table 1). Consistent with the down-regulation of evx1 in double mutants, we observe a 30 % reduction in the number of EGFP-labeled V0v cells in double mutant Tg(evx1:EGFP) SU2 embryos (22.9 +/- 3.4 cells in WTs; 16.0 +/- 5.5 cells in double mutants). However, strikingly, most V0v cells are still labeled with EGFP and these cells have a normal commissural ascending CoSA morphology (Fig. 4d & g). Unfortunately, we were never able to identify double mutant embryos that carried the Tg(evx1:EGFP) SU1 transgene, suggesting that this transgene probably integrated in the vicinity of the WT evx2 allele (see methods). However, consistent with our results using the Tg(evx1:EGFP) SU2 transgenic line, when we injected evx2 ATG morpholinos, at a concentration that eliminates Evx2 protein, into evx1 single mutant Tg(evx1:EGFP) SU1 embryos, we also observed EGFP-labeled V0v cells with normal CoSA axon trajectories (Fig. 4e).
Consistent with this persistence of V0v cells, we also found no change in the number of eng1b-expressing spinal cord cells in evx1 and evx2 single or double mutants compared to WT embryos at 24 h. If anything, we observed a slight reduction in the number of eng1b cells in the double mutants, although this was not statistically significant (Fig. 2f & o, Table 1). To further confirm that V0v cells were not adopting a V1 fate and turning on eng1b expression, we repeated this experiment at 30 h. We still found no change in eng1b expression in either single or double mutants when compared to WT embryos (Fig. 2h & p, Table 1).
Taken together, these results suggest that V0v cells are not transfating into V1 cells in zebrafish, even in the absence of both Evx1 and Evx2. Instead at least most of these cells are maintaining their V0v identities. In addition, these data suggest that Evx1 and Evx2 act partially redundantly to maintain each other’s expression, although only evx2 expression requires Evx1/Evx2 activity as more than half of V0v cells still express evx1 in evx1;evx2 double mutants.
Evx1 and Evx2 are required for skor2 expression in V0v cells
To identify additional transcription factors that might be required for specification of V0v functional characteristics, we expression-profiled FAC-sorted V0v cells (see Methods; ). From these analyses, we identified skor2 as a transcription factor gene potentially expressed by V0v neurons. Our subsequent in situ hybridization experiments demonstrated that skor2 has two clear domains of spinal cord expression, a ventral domain and a more dorsal domain (Fig. 7c). Double labeling experiments show that in the ventral domain, at least most of the skor2-expressing cells are V0v cells (Fig. 7a). On average, 97 % of ventral skor2-expressing cells co-express EGFP in Tg(evx1:EGFP) SU1 embryos (73/75 cells counted in 4 embryos). Given this high number of double positive cells and the fact that there is usually a delay in EGFP expression it is possible that all of the ventral skor2-expressing cells are V0v cells. Interestingly, double labeling experiments with skor2 and Tg(slc17a6:EGFP) demonstrated that both the ventral and dorsal skor2-expressing cells are excitatory cells (Fig. 7b), suggesting that Skor2 may play in role in specifying excitatory fates. As skor2 is expressed by V0v cells, we tested whether it is regulated by Evx1 and Evx2. We found that the number of cells expressing skor2 in the ventral spinal cord is reduced in evx1 and evx2 single mutants compared to WT embryos. More strikingly, ventral skor2 expression is completely abolished in double evx1;evx2 mutants (Figs 7d & e; Table 1), demonstrating that Evx1 and Evx2 are required, partially redundantly for skor2 expression in V0v cells. In contrast, there was no change in the dorsal expression of skor2 (Fig. 7d & e; Table 1).
Evx1 and Evx2 are required to specify the glutamatergic fates of V0v cells
Given that V0v cells still form in the absence of Evx1 and Evx2 function and their axon trajectories appear to be unaffected, but expression of a novel excitatory cell marker skor2 is lost, we decided to test if V0v cell neurotransmitter phenotypes were changed. When we examined the expression of slc17a6 (vglut) genes in embryos from a cross of evx1;evx2 heterozygous parents, we saw a significant reduction of glutamatergic cells in both of the single mutants when compared to WT embryos and an even more severe reduction in double mutants (Fig. 8a-d & o; Table 1). These data indicate that Evx1 and Evx2 act partially redundantly to specify the glutamatergic phenotype of V0v cells. Strikingly, the number of glutamatergic cells lost in the double mutant (approximately 28 cells in the spinal cord region adjacent to somites 6-10) is equivalent to the number of V0v cells in that region of the spinal cord (approximately 29 cells express evx2 and 33 cells express evx1 in this region of the WT spinal cord at this stage; see Table 1), suggesting that probably all of the V0v cells have lost their glutamatergic phenotype.
V0v cells become inhibitory in evx1;evx double mutant embryos
Given that V0v cells lose their excitatory phenotype in evx1;evx2 double mutants, we asked whether they are acquiring an inhibitory neurotransmitter fate instead. When we examined expression of slc32a1, which is expressed by all inhibitory neurons [8, 68, 69], there was no significant difference in the number of cells expressing this gene between either of the single mutants and WT embryos (Fig. 8e, g, h & p, Table 1). However, interestingly, there was a significant increase (approximately 21 cells) in the number of slc32a1-expressing cells in evx1;evx2 double mutants (Fig. 8f & p, Table 1). This suggests that Evx1 and Evx2 act redundantly to repress the inhibitory fate in V0v cells.
To further confirm that V0v cells are switching to an inhibitory fate, we performed in situ hybridization for slc32a1 plus immunohistochemistry for EGFP in embryos from a cross of double heterozygous parents that carry the Tg(evx1:EGFP) SU2 transgene. In WT embryos we see no co-expression of slc32a1 and EGFP (Fig. 5c). However, in double mutant embryos most V0v cells express slc32a1 (77 % of EGFP-positive V0v cells (30/39 cells counted in 2 embryos); Fig. 5f).
To determine whether V0v cells are becoming GABAergic and/or glycinergic we examined expression of markers of these two fates. We see no significant difference in the number of cells expressing GABAergic markers in single or double mutant embryos (Fig. 8m, n & r). In contrast, there is an increase in the number of cells expressing glycinergic markers. Interestingly, and in contrast to the slc32a1 (viaat) result, we see a slight increase in both single mutants as well as a more pronounced increase in double mutants (Fig. 8i-l & q, Table 1). The increase in double mutants (approximately 28 cells) suggests that all V0v cells are becoming glycinergic.
V0v cells become glycinergic through a novel Pax2-independent mechanism
All of the transcription factors that have been identified so far as specifying inhibitory spinal fates act through Pax2 [14–21]. In addition, Tlx1 and Tlx3, the only other transcription factors that have been identified as specifying excitatory spinal cord fates [4, 5, 8], work at least in part by down-regulating Pax2 . Therefore, we decided to test if V0v cells turn on Pax2 expression in evx1;evx2 double mutants. However, when we analyzed pax2a expression there was no significant difference in the number of cells expressing this gene in either the single or double mutants compared to WT embryos (Fig. 6b & e; Table 1). To further confirm this result, we performed immunohistochemistry using a Pax2 antibody that recognizes Pax2a and Pax2b . Again, we saw no significant change in the number of cells expressing Pax2 protein in single or double mutants (Fig. 6d & f; Table 1). Finally, we also performed double-labeling experiments for Pax2 and EGFP in WT and mutant embryos that carried the Tg(evx1:EGFP) SU2 transgene, using either the Pax2 antibody or in situ hybridization with a mix of pax2a, pax2b and pax8 probes. In each case we examined at least two WT embryos and two double homozygous mutants and we did not observe any double-labeled cells (Fig. 6h and data not shown).
Taken together these results show that V0v cells are becoming glycinergic through a Pax2-independent mechanism. This is the first time that a Pax2-independent mechanism of glycinergic specification has been identified in spinal cord neurons.
Evx genes are found in a wide range of animals ranging from corals to humans . They encode transcription factors that contain both DNA-binding homeobox and C-terminal repressor domains [73–77]. Amniotes have two Evx genes (Evx1 and Evx2) and teleosts, including zebrafish, have three (evx1, evx2 and eve1), although only evx1 and evx2 are expressed in the spinal cord ([26, 32, 33, 59, 64]; this paper Figs. 2, 3, Additional file 1: Results and Figure S1). Interestingly, the genomic positions of these evx genes (adjacent to specific Hox clusters) along with phylogenetic analyses strongly suggest that the third evx gene in teleosts (eve1) is not the result of the extra genome duplication in the teleost lineage ([26, 78, 79], Additional file 1). Instead, it is likely that all three of these genes originated from the two rounds of whole genome duplication that occurred early in the vertebrate lineage  and eve1 was later lost in the tetrapod lineage ([26, 78], Additional file 1).
Here we provide the first comprehensive analysis of the functions of Evx1 and Evx2 in spinal cord interneuron development in any vertebrate. We demonstrate that, within the spinal cord, both of these transcription factors are expressed exclusively by V0v cells. We also show that V0v cells are glutamatergic. These findings complement and extend those of Satou and colleagues , who reported that Evx2-expressing cells that develop from the dbx1b progenitor domain in zebrafish express the glutamateric maker slc17a6b and Talpalar and colleagues , who showed that mouse V0v cells express the glutamatergic marker slc17a6 (vglut2). In addition, we confirm and extend previous reports that suggested that V0v neurons extend commissural axons [23, 27, 35, 36] and we identify these neurons as CoSA interneurons. Interestingly, while Satou and colleagues  observed similar V0v cell morphologies at the stages that we have examined, at later stages of development they also saw descending and bifurcating commissural excitatory V0 cells, suggesting that V0v cells may diversify morphologically at later developmental stages .
In mouse Evx1 mutants, most V0v cells completely change their fate and acquire characteristics of V1 cells, the cell-type that normally forms ventral to V0v cells. Cells that would have formed V0v interneurons lose expression of Evx1 and Evx2 and instead express the V1 marker Engrailed1 (En1) and develop axon trajectories and migration patterns characteristic of V1 interneurons . Most notably, their axon trajectories change from being contralateral to ipsilateral . In addition, experiments in chick embryos revealed that ectopic Evx is sufficient to suppress Engrailed expression and therefore presumably V1 cell fate . The role of Evx2 in mouse V0v cells is less well understood. Evx2 expression is dependent on Evx1, suggesting that it may be involved in the specification events described above , but spinal cord phenotypes of Evx2 mutants have not been described in mouse and before this study, Evx1;Evx2 double mutants had not been described in any vertebrate.
These amniote data suggest that Evx1 is required to inhibit the V1 fate in post-mitotic V0v cells. This global cell fate change is unusual for a transcription factor expressed in post-mitotic cells: it is more commonly seen with transcription factors expressed in spinal progenitor domains (e.g. [27–29, 81, 82]). For example, Nkx2.2 is a transcription factor expressed in the p3 progenitor domain and in Nkx2.2 mutant mice, cells that would have formed V3 interneurons change their fates (transfate) and become motoneurons instead . Similarly, Dbx1 is expressed in the progenitor domain (p0) from which V0 cells develop and in Dbx1 mutant mice, cells that would have become V0v cells assume the characteristics of V1 cells [27, 28].
Interestingly, our results are different and yield novel insights into Evx1 and Evx2 function in the spinal cord. We see no evidence of V0v cells transfating to V1 cells in zebrafish evx1 and evx2 single or double mutants. Notably, there is no increase in the number of cells expressing eng1b, which is specifically expressed in V1 cells, or Pax2, which is expressed by V1 and V0D cells. Instead most V0v cells continue to express evx1 mRNA and Tg(evx1:EGFP) and these EGFP-labelled V0v neurons have what appear to be normal CoSA axon morphologies. These data strongly suggest that V0v cells still form in zebrafish evx1;evx2 double mutants and that they do not become a different class of neuron. One possible explanation for the differences between our results and the previously reported analyses in mouse might be evolutionary changes in the functions of Evx1 and Evx2. However, it is also possible that the consequences of removing these transcription factors are different in mouse and zebrafish because of variations in the expression of Dbx1 and/or the timing of V0v development. Interestingly, in this paper, we have shown that expression of dbx1a and dbx1b persists in at least some V0v cells in zebrafish. Therefore, it is possible that in zebrafish evx1;evx2 mutants, Dbx might be able to inhibit post-mitotic V0v cells from becoming V1 cells (Fig. 9). Given the speed of zebrafish spinal cord development it is also possible that V0v cells become committed to their fate faster than in mouse and that, therefore, the window of time during which the V1 fate needs to be inhibited in V0v cells is much shorter in zebrafish than in mammals (see  for a different example of how the fast speed of zebrafish development can produce changes in spinal cord development). Regardless of how V0v global cell fate specification has evolved, the fact that V0v cells still form in zebrafish lacking Evx1 and Evx2, has provided us with a unique opportunity to identify Evx functions in V0v cells, independent of any role that these transcription factors may also have in repressing the V1 cell fate.
Our results show that in zebrafish, Evx1 and Evx2 act partially redundantly to specify the glutamatergic fate of V0v cells and inhibit an alternative glycinergic fate in these cells. Given that the only spinal cord cells that express evx1 and evx2 are V0v cells and that the only other trunk tissue that expresses either of these genes is the posterior gut, which expresses evx1, we consider that this requirement for Evx1 and Evx2 function is likely to be cell-autonomous. Interestingly, while there is a reduction of glutamatergic cells in both single and double mutants, expression of the inhibitory marker slc32a1 is only increased in double mutants, suggesting that the specification of glutamatergic fates and the inhibition of glycinergic fates may be independent processes which require different levels of Evx activity. However, in contrast to slc32a1, the number of cells expressing the glycinergic marker slc6a5 was slightly increased in single mutants, which suggests that the expression of different neurotransmitter transporter proteins is regulated independently and by distinct levels of Evx activity. These results are intriguing as they suggest that the regulation of neurotransmitter transporters and enzymes might be complex, with different components being regulated by distinct mechanisms.
V0v interneurons are a crucial part of locomotor circuitry as they are required for hindlimb left-right alternation during fast locomotion [9, 27, 34]. Therefore, changing the neurotransmitter fate of these cells might be expected to impair fast movements. Unfortunately, as evx2 mutants die by larval stages, we were not able to assess whether evx2 single mutants or evx1;evx2 double mutants have locomotion defects. In addition, evx1 single mutants lack joints in their fins , making it impossible to evaluate if any difference in evx1 single mutant behavior is due to this fin phenotype or a locomotive defect. Interestingly, we did not observe any obvious changes in V0v cell morphology or axon trajectory in evx1;evx2 double mutants. Given the changed neurotransmitter phenotype of V0v cells in these animals this might be considered surprising, although it is consistent with our previous analysis of V1 cells, that maintain their ipsilateral axon trajectories even when they lose their inhibitory fates in the absence of Pax2 and Pax8 . It is still possible though that there are subtle changes in V0v cell wiring and/or changes in V0v cell connectivity in evx1;evx2 double mutants as a result of their change in neurotransmitter fate. As there are fewer GFP-labelled V0v cells in evx1;evx2 double mutants it is also possible that V0v cells with inappropriate neurotransmitter fates eventually die, although alternatively this reduction in the number of GFP-positive cells may just reflect autoregulation of Evx expression.
In this paper, we also describe the expression of a different transcription factor gene expressed by V0v cells, skor2. Skor2 expression has also been reported in the mouse spinal cord but the cells that express it were not identified . Our results show that skor2 is expressed by a subset of V0v cells as well as at least one population of more dorsal excitatory spinal cord cells. We also demonstrate that expression of skor2 in V0v cells requires Evx1 and Evx2 activity. Given that skor2 is predominantly expressed by excitatory cells, it is possible that it acts downstream of Evx1 and Evx2 in V0v cells in either the specification of glutamatergic fates and/or the inhibition of glycinergic fates and that it might also have this function in other cells. However it is also possible that Skor2 acts downstream of Evx1 and Evx2 in some other as-yet-unidentified aspect of V0v cell specification. These alternatives can be tested by future loss-of-function analyses of Skor2.
Excitingly, in addition to demonstrating the roles of Evx1 and Evx2 in neurotransmitter specification, our data also show that these transcription factors function independently of Pax2 in specifying glutamatergic fates and inhibiting glycinergic fates. This is the first time that a Pax2-independent mechanism of inhibitory fate specification has been identified in the spinal cord. While several transcription factors have been identified that specify the inhibitory fates of particular spinal cord neurons, so far all of these act upstream of Pax2 [14–18]. In addition, as mentioned before, the only other transcription factors that have been identified as specifying excitatory spinal cord fates, Tlx1 and Tlx3 [4, 5, 8], work at least in part by down-regulating Pax2 . Therefore, our finding that V0v cells become inhibitory in evx1;evx2 double mutants but do not express Pax2 is a significant one as it demonstrates that there must be an additional Pax2-independent mechanism for specifying inhibitory neurons in the spinal cord.
In conclusion, in this paper we demonstrate that zebrafish V0v cells express evx1 and evx2 and develop into excitatory (glutamatergic) CoSA interneurons. We also show that Evx1 and Evx2 are required, partially redundantly for expression of skor2 and glutamatergic markers and inhibition of glycinergic markers in V0v cells and that in the absence of Evx1 and Evx2 function V0v cells become glycinergic through a novel Pax2-independent mechanism. Taken together, our data significantly increase our understanding of how neurotransmitter fates are specified and the genetic networks involved in these processes.
conserved non-coding element
central nervous system
commisural secondary ascending
days post fertilization
fluorescent activated cell sorting
glutamic acid decarboxylases
hours post fertilization
left hand side
right hand side
Goulding M. Circuits controlling vertebrate locomotion: moving in a new direction. Nat Rev Neurosci. 2009;10(7):507–18.
Lewis KE. How do genes regulate simple behaviours? Understanding how different neurons in the vertebrate spinal cord are genetically specified. Philos Trans R Soc Lond B Biol Sci. 2006;361(1465):45–66.
Gosgnach S, Lanuza GM, Butt SJ, Saueressig H, Zhang Y, Velasquez T, et al. V1 spinal neurons regulate the speed of vertebrate locomotor outputs. Nature. 2006;440(7081):215–9.
Li WC, Higashijima S, Parry DM, Roberts A, Soffe SR. Primitive roles for inhibitory interneurons in developing frog spinal cord. J Neurosci. 2004;24(25):5840–8.
Higashijima S, Masino M, Mandel G, Fetcho JR. Engrailed-1 Expression Marks a Primitive Class of Inhibitory Spinal Interneuron. J Neurosci. 2004;24(25):5827–39.
Crone SA, Zhong G, Harris-Warrick R, Sharma K. In mice lacking V2a interneurons, gait depends on speed of locomotion. J Neurosci. 2009;29(21):7098–109.
Crone SA, Quinlan KA, Zagoraiou L, Droho S, Restrepo CE, Lundfald L, et al. Genetic ablation of V2a ipsilateral interneurons disrupts left-right locomotor coordination in mammalian spinal cord. Neuron. 2008;60(1):70–83.
Kimura Y, Okamura Y, Higashijima S. alx, a zebrafish homolog of Chx10, marks ipsilateral descending excitatory interneurons that participate in the regulation of spinal locomotor circuits. J Neurosci. 2006;26(21):5684–97.
Talpalar AE, Bouvier J, Borgius L, Fortin G, Pierani A, Kiehn O. Dual-mode operation of neuronal networks involved in left-right alternation. Nature. 2013;500(7460):85–8.
Xue M, Atallah BV, Scanziani M. Equalizing excitation-inhibition ratios across visual cortical neurons. Nature. 2014;511(7511):596–600.
Zadori D, Veres G, Szalardy L, Klivenyi P, Toldi J, Vecsei L. Glutamatergic dysfunctioning in Alzheimer's disease and related therapeutic targets. J Alzheimers Dis. 2014;42 Suppl 3:S177–187.
Bateup HS, Johnson CA, Denefrio CL, Saulnier JL, Kornacker K, Sabatini BL. Excitatory/inhibitory synaptic imbalance leads to hippocampal hyperexcitability in mouse models of tuberous sclerosis. Neuron. 2013;78(3):510–22.
Pittenger C, Bloch MH, Williams K. Glutamate abnormalities in obsessive compulsive disorder: neurobiology, pathophysiology, and treatment. Pharmacol Ther. 2011;132(3):314–32.
Batista MF, Lewis KE. Pax2/8 act redundantly to specify glycinergic and GABAergic fates of multiple spinal interneurons. Dev Biol. 2008;323(1):88–97.
Cheng L, Samad OA, Xu Y, Mizuguchi R, Luo P, Shirasawa S, et al. Lbx1 and Tlx3 are opposing switches in determining GABAergic versus glutamatergic transmitter phenotypes. Nat Neurosci. 2005;8(11):1510–5.
Mizuguchi R, Kriks S, Cordes R, Gossler A, Ma Q, Goulding M. Ascl1 and Gsh1/2 control inhibitory and excitatory cell fate in spinal sensory interneurons. Nat Neurosci. 2006;9(6):770–8.
Glasgow SM, Henke RM, Macdonald RJ, Wright CV, Johnson JE. Ptf1a determines GABAergic over glutamatergic neuronal cell fate in the spinal cord dorsal horn. Development. 2005;132(24):5461–9.
Pillai A, Mansouri A, Behringer R, Westphal H, Goulding M. Lhx1 and Lhx5 maintain the inhibitory-neurotransmitter status of interneurons in the dorsal spinal cord. Development. 2007;134(2):357–66.
Guo Z, Zhao C, Huang M, Huang T, Fan M, Xie Z, et al. Tlx1/3 and Ptf1a control the expression of distinct sets of transmitter and peptide receptor genes in the developing dorsal spinal cord. J Neurosci. 2012;32(25):8509–20.
Huang M, Huang T, Xiang Y, Xie Z, Chen Y, Yan R, et al. Ptf1a, Lbx1 and Pax2 coordinate glycinergic and peptidergic transmitter phenotypes in dorsal spinal inhibitory neurons. Dev Biol. 2008;322(2):394–405.
Hori K, Cholewa-Waclaw J, Nakada Y, Glasgow SM, Masui T, Henke RM, et al. A nonclassical bHLH Rbpj transcription factor complex is required for specification of GABAergic neurons independent of Notch signaling. Gene Dev. 2008;22(2):166–78.
Cheng L, Arata A, Mizuguchi R, Qian Y, Karunaratne A, Gray PA, et al. Tlx3 and Tlx1 are post-mitotic selector genes determining glutamatergic over GABAergic cell fates. Nat Neurosci. 2004;7(5):510–7.
Moran-Rivard L, Kagawa T, Saueressig H, Gross MK, Burrill J, Goulding M. Evx1 is a postmitotic determinant of V0 interneuron identity in the spinal cord. Neuron. 2001;29(2):385–99.
Dolle P, Fraulob V, Duboule D. Developmental expression of the mouse Evx-2 gene: relationship with the evolution of the HOM/Hox complex. Dev Suppl. 1994;143–153.
Bastian H, Gruss P. A murine even-skipped homologue, Evx 1, is expressed during early embryogenesis and neurogenesis in a biphasic manner. EMBO J. 1990;9(6):1839–52.
Avaron F, Thaeron-Antono C, Beck CW, Borday-Birraux V, Geraudie J, Casane D, et al. Comparison of even-skipped related gene expression pattern in vertebrates shows an association between expression domain loss and modification of selective constraints on sequences. Evol Dev. 2003;5(2):145–56.
Lanuza G, Gosgnach S, Pierani A, Jessel T, Goulding M. Genetic Identification of Spinal Interneurons that Coordinate Left-Right Locomotor Activity Necessary for Walking Movements. Neuron. 2004;42:375–86.
Pierani A, Moran-Rivard L, Sunshine MJ, Littman DR, Goulding M, Jessell TM. Control of interneuron fate in the developing spinal cord by the progenitor homeodomain protein Dbx1. Neuron. 2001;29(2):367–84.
Briscoe J, Pierani A, Jessell TM, Ericson J. A homeodomain protein code specifies progenitor cell identity and neuronal fate in the ventral neural tube. Cell. 2000;101(4):435–45.
Ruiz i Altaba A, Melton DA. Bimodal and graded expression of the Xenopus homeobox gene Xhox3 during embryonic development. Development. 1989;106(1):173–83.
Ruiz i Altaba A. Neural expression of the Xenopus homeobox gene Xhox3: evidence for a patterning neural signal that spreads through the ectoderm. Development. 1990;108(4):595–604.
Sordino P, Duboule D, Kondo T. Zebrafish Hoxa and Evx-2 genes: cloning, developmental expression and implications for the functional evolution of posterior Hox genes. Mech Dev. 1996;59:165–75.
Thaeron C, Avaron F, Casane D, Borday V, Thisse B, Thisse C, et al. Zebrafish evx1 is dynamically expressed during embryogenesis in subsets of interneurones, posterior gut and urogenital system. Mech Dev. 2000;99:167–72.
Griener A, Zhang W, Kao H, Wagner C, Gosgnach S. Probing diversity within subpopulations of locomotor-related V0 interneurons. Developmental neurobiology. 2015.
Suster ML, Kania A, Liao M, Asakawa K, Charron F, Kawakami K, et al. A novel conserved evx1 enhancer links spinal interneuron morphology and cis-regulation from fish to mammals. Dev Biol. 2009;325(2):422–33.
Satou C, Kimura Y, Higashijima S. Generation of multiple classes of V0 neurons in zebrafish spinal cord: progenitor heterogeneity and temporal control of neuronal diversity. J Neurosci. 2012;32(5):1771–83.
Bae YK, Kani S, Shimizu T, Tanabe K, Nojima H, Kimura Y, et al. Anatomy of zebrafish cerebellum and screen for mutations affecting its development. Dev Biol. 2009;330(2):406–26.
Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF. Stages of embryonic development of the zebrafish. Dev Dyn. 1995;203(3):253–310.
Schulte CJ, Allen C, England SJ, Juarez-Morales JL, Lewis KE. Evx1 is required for joint formation in zebrafish fin dermoskeleton. Dev Dyn. 2011;240(5):1240–8.
Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, et al. LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003;13(4):721–31.
Mayor C, Brudno M, Schwartz JR, Poliakov A, Rubin EM, Frazer KA, et al. VISTA : visualizing global DNA sequence alignments of arbitrary length. Bioinformatics. 2000;16(11):1046–7.
Sasaki Y, Sone T, Yoshida S, Yahata K, Hotta J, Chesnut JD, et al. Evidence for high specificity and efficiency of multiple recombination signals in mixed DNA cloning by the Multisite Gateway system. J Biotechnol. 2004;107(3):233–43.
Suzuki Y, Kagawa N, Fujino T, Sumiya T, Andoh T, Ishikawa K, et al. A novel high-throughput (HTP) cloning strategy for site-directed designed chimeragenesis and mutation using the Gateway cloning system. Nucleic Acids Res. 2005;33(12):1–6.
Villefranc JA, Amigo J, Lawson ND. Gateway compatible vectors for analysis of gene function in the zebrafish. Dev Dyn. 2007;236(11):3077–87.
Köster RW, Fraser SE. Tracing transgene expression in living zebrafish embryos. Dev Biol. 2001;233(2):329–46.
Wang Y, Shen J, Arenzana N, Tirasophon W, Kaufman RJ, Prywes R. Activation of ATF6 and an ATF6 DNA binding site by the endoplasmic reticulum stress response. J Biol Chem. 2000;275(35):27013–20.
Kawakami K, Takeda H, Kawakami N, Kobayashi M, Matsuda N, Mishina M. A transposon-mediated gene trap approach identifies developmentally regulated genes in zebrafish. Dev Cell. 2004;7(1):133–44.
Halpern ME, Rhee J, Goll MG, Akitake CM, Parsons M, Leach SD. Gal4/UAS transgenic tools and their application to zebrafish. Zebrafish. 2008;5(2):97–110.
Cerda GA, Hargrave M, Lewis KE. RNA profiling of FAC-sorted neurons from the developing zebrafish spinal cord. Dev Dyn. 2009;238(1):150–61.
England S, Hilinski W, de Jager S, Andrzejczuk L, Campbell P, Chowdhury T, Demby C, Fancher W, Gong Y, Lin C et al: Identifying Transcription Factors expressed by Ventral Spinal Cord Interneurons. ZFIN on-line publication 2014, http://zfin.org/ZDB-PUB-140822-10.
Park HC, Kim CH, Bae YK, Yeo SY, Kim SH, Hong SK, et al. Analysis of upstream elements in the HuC promoter leads to the establishment of transgenic zebrafish with fluorescent neurons. Dev Biol. 2000;227(2):279–93.
Picker A, Scholpp S, Bohli H, Takeda H, Brand M. A novel positive transcriptional feedback loop in midbrain-hindbrain boundary development is revealed through analysis of the zebrafish pax2.1 promoter in transgenic lines. Development. 2002;129:3227–39.
Peri F, Nusslein-Volhard C. Live imaging of neuronal degradation by microglia reveals a role for v0-ATPase a1 in phagosomal fusion in vivo. Cell. 2008;133(5):916–27.
Kobayashi M, Nishikawa K, Yamamoto M. Hematopoietic regulatory domain of gata1 gene is positively regulated by GATA1 protein in zebrafish embryos. Development. 2001;128(12):2341–50.
Batista MF, Jacobstein J, Lewis KE. Zebrafish V2 cells develop into excitatory CiD and Notch signalling dependent inhibitory VeLD interneurons. Dev Biol. 2008;322(2):263–75.
Concordet JP, Lewis KE, Moore JW, Goodrich LV, Johnson RL, Scott MP, et al. Spatial regulation of a zebrafish patched homologue reflects the roles of sonic hedgehog and protein kinase A in neural tube and somite patterning. Development. 1996;122(9):2835–46.
Gribble SL, Nikolaus OB, Dorsky RI. Regulation and function of Dbx genes in the zebrafish spinal cord. Dev Dyn. 2007;236(12):3472–83.
Thaëron C, Avaron F, Casane D, Borday V, Thisse B, Thisse C, et al. Zebrafish evx1 is dynamically expressed during embryogenesis in subsets of interneurones, posterior gut and urogenital system. Mech Dev. 2000;99(1-2):167–72.
Joly JS, Joly C, Schulte-Merker S, Boulekbache H, Condamine H. The ventral and posterior expression of the zebrafish homeobox gene eve1 is perturbed in dorsalized and mutant embryos. Development. 1993;119(4):1261–75.
Pfeffer PL, Gerster T, Lun K, Brand M, Busslinger M. Characterization of three novel members of the zebrafish Pax2/5/8 family: dependency of Pax5 and Pax8 expression on the Pax2.1 (noi) function. Development. 1998;125(16):3063–74.
Higashijima S, Mandel G, Fetcho JR. Distribution of prospective glutamatergic, glycinergic, and GABAergic neurons in embryonic and larval zebrafish. J Comp Neurol. 2004;480(1):1–18.
Higashijima S, Schaefer M, Fetcho JR. Neurotransmitter properties of spinal interneurons in embryonic and larval zebrafish. J Comp Neurol. 2004;480(1):19–37.
Pierani A, Brenner-Morton S, Chiang C, Jessell TM. A Sonic Hedgehog Independent Retinoid-Activated Pathway of Neurogenesis in the Ventral Spinal Cord. Cell. 1999;97:903–15.
Cruz C, Maegawa S, Weinberg ES, Wilson SW, Dawid IB, Kudoh T. Induction and patterning of trunk and tail neural ectoderm by the homeobox gene eve1 in zebrafish embryos. Proc Natl Acad Sci U S A. 2010;107(8):3564–9.
Lewis KE, Eisen JS. From cells to circuits: development of the zebrafish spinal cord. Prog Neurobiol. 2003;69(6):419–49.
Bernhardt RR, Chitnis AB, Lindamer L, Kuwada JY. Identification of spinal neurons in the embryonic and larval zebrafish. J Comp Neurol. 1990;302:603–16.
Satou C, Kimura Y, Hirata H, Suster ML, Kawakami K, Higashijima S. Transgenic tools to characterize neuronal properties of discrete populations of zebrafish neurons. Development. 2013;140(18):3927–31.
McIntire SL, Reimer RJ, Schuske K, Edwards RH, Jorgensen EM. Identification and characterization of the vesicular GABA transporter. Nature. 1997;389(6653):870–6.
Sagne C, El Mestikawy S, Isambert MF, Hamon M, Henry JP, Giros B, et al. Cloning of a functional vesicular GABA and glycine transporter by screening of genome databases. FEBS letters. 1997;417(2):177–83.
Rabe N, Gezelius H, Vallstedt A, Memic F, Kullander K. Netrin-1-dependent spinal interneuron subtypes are required for the formation of left-right alternating locomotor circuitry. J Neurosci. 2009;29(50):15642–9.
Burrill JD, Moran L, Goulding MD, Saueressig H. PAX2 is expressed in multiple spinal cord interneurons, including a population of EN1+ interneurons that require PAX6 for their development. Development. 1997;124(22):4493–503.
Ahringer J. Posterior patterning by the Caenorhabditis elegans even-skipped homolog vab-7. Gene Dev. 1996;10(9):1120–30.
Briata P, Ilengo C, Van DeWerken R, Corte G. Mapping of a potent transcriptional repression region of the human homeodomain protein EVX1. FEBS letters. 1997;402(2-3):131–5.
Biggin MD, Tjian R. A purified Drosophila homeodomain protein represses transcription in vitro. Cell. 1989;58(3):433–40.
TenHarmsel A, Austin RJ, Savenelli N, Biggin MD. Cooperative binding at a distance by even-skipped protein correlates with repression and suggests a mechanism of silencing. Mol Cell Biol. 1993;13(5):2742–52.
Li C, Manley JL. Even-skipped represses transcription by binding TATA binding protein and blocking the TFIID-TATA box interaction. Mol Cell Biol. 1998;18(7):3771–81.
Han K, Manley JL. Functional domains of the Drosophila Engrailed protein. EMBO J. 1993;12:2723–33.
Amores A, Force A, Yan YL, Joly L, Amemiya C, Fritz A, et al. Zebrafish hox clusters and vertebrate genome evolution. Science. 1998;282(5394):1711–4.
Taylor JS, Braasch I, Frickey T, Meyer A, Van de Peer Y. Genome duplication, a trait shared by 22000 species of ray-finned fish. Genome Res. 2003;13(3):382–90.
Holland PW, Garcia-Fernandez J, Williams NA, Sidow A. Gene duplications and the origins of vertebrate development. Dev Suppl. 1994;125–133.
Ericson J, Rashbass P, Schedl A, Brenner-Morton S, Kawakami A, van Heyningen V, et al. Pax6 controls progenitor cell identity and neuronal fate in response to graded Shh signaling. Cell. 1997;90:169–80.
Sander M, Paydar S, Ericson J, Briscoe J, Berber E, German M, et al. Ventral neural patterning by Nkx homeobox genes: Nkx6.1 controls somatic motor neuron and ventral interneuron fates. Gene Dev. 2000;14(17):2134–9.
Briscoe J, Sussel L, Serup P, Hartigan-O'Connor D, Jessell TM, Rubenstein JL, et al. Homeobox gene Nkx2.2 and specification of neuronal identity by graded Sonic hedgehog signalling. Nature. 1999;398(6728):622–7.
England S, Batista MF, Mich JK, Chen JK, Lewis KE. Roles of Hedgehog pathway components and retinoic acid signalling in specifying zebrafish ventral spinal cord neurons. Development. 2011;138(23):5121–34.
Arndt S, Poser I, Schubert T, Moser M, Bosserhoff AK. Cloning and functional characterization of a new Ski homolog, Fussel-18, specifically expressed in neuronal tissues. Lab Investig. 2005;85(11):1330–41.
We thank Dr Shinichi Higashijima for kindly providing us with Evx2 antibody, Dr Derek Stemple and everyone working on the zebrafish mutation project at the Wellcome Trust Sanger Centre for the evx2 mutant allele, ZFIN for providing information on nomenclature and other essential zebrafish resources, Nigel Miller at the Flow Cytometry Facility, Department of Pathology, University of Cambridge, Cambridge, for his expert FAC-sorting, Tomi Ivacevic and Vladimir Benes at the Genomics Core Facility at EMBL, Heidelberg, for RNA amplification, labelling, and hybridization, Uwe Strähle, Olivier Armant and Jasmin Lampert for help with designing microarrays, Nicole Santos for help with genotyping, Henry Putz, Jessica Bouchard, Annika Swanson and several SU undergraduate fish husbandry workers for help with maintaining zebrafish lines and Dr Santanu Banerjee for helpful comments on a previous version of this manuscript.
This work was supported by NSF IOS-1257583, NIH NINDS R21NS073979, the Spinal Cord Injury Trust Fund through New York State Department of Health Contract #C030177, a Royal Society University Research Fellowship and Syracuse University start-up funds awarded to K.E.L; BBSRC, Cambridge Trust, DAAD and Daimler-Benz-Foundation funding awarded to C.J.S; NSF HRD-0703452 LSAMP and Syracuse University Ruth Meyer funding awarded to S.A.P and NSF HRD-1202480 LSAMP and Syracuse University Ruth Meyer funding awarded to G.K.V.
The authors have no competing interests.
JM created the Tg(evx:EGFP) lines and performed most of the experiments in the paper, including most of the single and double labeling experiments in WT, single mutant and double mutant embryos and all of the cell counts except those for the skor2 double labeling experiments; CS performed initial experiments with the evx1 mutant and evx1 and evx2 morpholinos and helped to formulate many of the hypotheses that this study tests; CS also performed many of the WT expression analyses of evx1, evx2 and eve1, including some of the initial double labeling experiments; SE performed the synteny and phylogenetic analyses; SP and GV genotyped embryos and adult fish and helped to establish stable Tg(evx:EGFP) WT and mutant lines; SP and GV also examined Evx2 expression in WT and evx2 single mutant embryos and performed the evx2 mutant survival experiments; SP also performed the double staining for Evx2 and EGFP in Tg(slc17a6:EGFP) embryos; WH performed the initial analysis of skor2 and the double labeling experiments between skor2 and other markers; SE and SdJ performed the FACS and expression profiling analyses; KL designed and directed the study; KL and JM wrote most of the paper with input from the other authors.
Supplementary Data. (PDF 2.67 mb)