- Research article
- Open Access
Lmx1b is required for the glutamatergic fates of a subset of spinal cord neurons
- William C. Hilinski1, 2,
- Jonathan R. Bostrom†3,
- Samantha J. England†1,
- José L. Juárez-Morales†1,
- Sarah de Jager4,
- Olivier Armant5,
- Jessica Legradi5,
- Uwe Strähle5,
- Brian A. Link3 and
- Katharine E. Lewis1Email author
- Received: 20 April 2016
- Accepted: 8 August 2016
- Published: 23 August 2016
Abstract
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.
Keywords
- Spinal cord
- Interneuron
- Zebrafish
- Lmx1b
- Excitatory
- Neurotransmitter
- CNS
- Transcription factor
- V0v
- dI5
Background
Neurons in the central nervous system (CNS) must specify and maintain several properties in order to integrate and function properly within neuronal circuitry [1]. One crucial neuronal characteristic that must be specified correctly and usually must be maintained (for some exceptions see [2]) is the neurotransmitter phenotype [1]. Failure to correctly specify or maintain neurotransmitter phenotypes can result in incorrect levels of excitatory or inhibitory neurotransmitter release and lead to diseases such as epilepsy, autism spectrum disorder, and Alzheimer’s [3–6].
Neurotransmitter phenotypes, like many other neuronal properties, are initially specified by transcription factors that individual neurons express as they start to differentiate [7–12]. These neurotransmitter phenotypes are then maintained either by these same transcription factors or by additional ones [7, 13–17]. However, in many cell types the transcription factors that specify and/or maintain neurotransmitter phenotypes are still unknown. This is a critical gap in our knowledge and one that we need to address in order to potentially develop better treatments for some of the aforementioned diseases and disorders.
In this paper, we investigate the functions of Lmx1b transcription factors in the zebrafish spinal cord. Lmx1b has been implicated in a variety of functions in different regions of the vertebrate CNS including cell migration, cell survival, as well as correct specification and/or maintenance of cell identity, neuronal connectivity and neurotransmitter phenotypes [18–25]. However, it remains unclear if Lmx1b is required for neurotransmitter specification and/or maintenance in the spinal cord.
Zebrafish have two Lmx1b ohnologs, lmx1ba and lmx1bb, that we show are probably expressed in overlapping spinal cord domains. Consistent with previous analyses in mouse, we show that lmx1bb is expressed by dI5 neurons, and for the first time in any animal, we show that V0v neurons (cells that form in the ventral part of the V0 domain [11, 12, 26–31]) also express lmx1bb. Both dI5 and V0v cells are glutamatergic [8, 11, 16, 31, 32] and consistent with this we demonstrate that the vast majority of lmx1bb-expressing cells are glutamatergic.
We also show in zebrafish lmx1bb homozygous mutants that glutamatergic neurons are correctly specified during early development but are reduced in number at later developmental time points. Interestingly, we see the same phenotype in lmx1ba homozygous mutants, lmx1ba;lmx1bb double mutants and lmx1ba;lmx1bb double heterozygous embryos suggesting that lmx1ba and lmx1bb act at least partially redundantly in a dose-dependent manner and that three functional lmx1b alleles are required for the specification or maintenance of correct numbers of spinal cord glutamatergic cells at later developmental stages. In contrast to the reduction in the number of glutamatergic neurons, there is no reduction in the numbers of V0v or dI5 cells in lmx1bb homozygous mutants and there is no increase in cell death. This suggests that lmx1b-expressing spinal neurons are still present in normal numbers at these later stages of development, but that fewer of them are glutamatergic. Interestingly, there is no increase in the number of inhibitory neurons, suggesting that the cells that are no longer excitatory do not become inhibitory. Finally, we demonstrate that lmx1ba and lmx1bb expression in V0v cells requires Evx1 and Evx2. In combination with a previous study that showed that Evx1 and Evx2 are required for V0v cells to become glutamatergic [11], this suggests that Lmx1ba and Lmx1bb act downstream of Evx1 and Evx2 either to maintain V0v glutamatergic fates or to specify the glutamatergic fates of a later-forming subset of V0v cells.
Methods
Zebrafish husbandry and fish lines
Zebrafish (Danio rerio) were maintained on a 14-h light/10-h dark cycle at 28.5 °C. Embryos were obtained from natural paired and/or grouped spawnings of wild-type (WT) (AB, TL or AB/TL hybrid) fish or identified heterozygous lmx1bb jj410 , lmx1ba mw80 , evx1 i232 ;evx2 sa140 or smoothened b641 mutant fish or Tg(slc17a6:EGFP) [formerly called Tg(vGlut2a:EGFP)] [33] or Tg(evx1:EGFP) SU1 [11] transgenic fish or lmx1bb jj410 crossed into the background of either Tg(slc17a6b(vglut2a):loxP-DsRed-loxP-GFP) nns14 [41, 42] or Tg(evx1:EGFP) SU1 fish respectively. Embryos were reared at 28.5 °C and staged by hours post fertilization (h) and/or days post fertilization (dpf). Most embryos were treated with 0.2 mM 1-phenyl 2-thiourea (PTU) at 24 h to inhibit melanogenesis [34–36].
The evx1 i232 , evx2 sa140 and lmx1bb jj410 mutants have been previously described [11, 37–39]. All three of these mutations are single base pair changes that lead to premature stop codons before the homeobox. Therefore, if any of these RNAs are not degraded by nonsense mediated decay, the resulting proteins will lack the DNA binding domain. lmx1ba mw80 mutant zebrafish were generated using TALENs constructs that target the sequences TCAAGTAGACATGCTGGACG and TCCGCTCCTGTCCTGAACTG within the first exon of lmx1ba. Constructs were made using steps 1–38 outlined in [40]. To generate mRNA encoding the TALENs, approximately 5 μg of plasmid DNA was digested with ApoI and purified via the Invitrogen PureLink PCR Purification Kit (ThermoFisher, K310001). RNA was synthesized using the Ambion mMessage mMachine T7 kit (ThermoFisher, AM1344) with a poly(A) tail added from the Poly(A) Tailing Kit (Ambion, AM1350) and purified with the Megaclear Kit (Ambion, AM1908). 100 pg of RNA for each TALEN was co-injected into 1-cell WT embryos. The lmx1ba mw80 allele was recovered and identified as a single base pair deletion 20 bp into the coding sequence. This results in a frameshift after the first six amino acids and a premature stop codon 11 amino acids later. This stop codon is upstream of both the Lim and homeobox domains, suggesting that this allele is likely to be a complete loss of function.
Genotyping
DNA for genotyping was isolated from both anesthetized adults and fixed embryos via fin biopsy or head dissections respectively. Fin biopsy and evx1 and evx2 genotyping of adults were performed as previously described [11, 37]. KASP assays, designed by LGC Genomics LLC, using DNA extracted from head dissections, were used to identify embryos carrying the evx1 i232 and evx2 sa140 mutations. These assays use allele-specific PCR primers which differentially bind fluorescent dyes that we quantified with a BioRad CFX96 real-time PCR machine to distinguish genotypes. The proprietary primers used are: Evx1_y32_i232 and Evx2_sa140.
Heads of fixed embryos were dissected in 80 % glycerol/20 % phosphate-buffered saline (PBS) with insect pins. Embryo trunks were stored in 70 % glycerol/30 % PBS at 4 °C for later analysis. DNA was extracted via the HotSHOT method [41] using 20 μL of NaOH and 2 μL of Tris-HCl (pH-7.5).
The lmx1ba mw80 and lmx1bb jj410 alleles were identified by restriction enzyme digestion assays as both of these mutations disrupt endogenous restriction enzyme sites. For lmx1ba mw80 , a 540 bp amplicon encompassing the mutation site was generated with the following primers:
Forward GATCCTCAAGAGGAGCTCATACACA and Reverse CATGCACATTTAACTATGATCTGAGCCGTG.
This amplicon was digested with MluCI to yield 311 bp and 142 bp and 87 bp (WT), 453 bp and 87 bp (homozygous mutant), or 453 bp and 311 bp and 142 bp and 87 bp (heterozygous mutant) products. Similarly, for lmx1bb jj410 , a 264 bp amplicon encompassing the mutation site was generated with the following primers: Forward GAAGGCTCGTCTCTGCTGTGTGGTG and Reverse CGTTATGGATGCGCTGAGACTGAATACC. This amplicon was digested with BfaI to yield 211 bp and 53 bp (WT), 264 bp (homozygous mutant), or 264 bp and 211 bp and 53 bp (heterozygous mutant) products.
Expression profiling V0v neurons & microarray design
To identify transcription factors expressed by V0v neurons, V0v spinal neurons, all spinal cord neurons and all cells within the trunk were isolated from live transgenic zebrafish embryos at 27 h using fluorescence activated cell-sorting (FACS). Prior to FACS, embryos were prim-staged, de-yolked, dissected and dissociated as in [42, 43]. Heads and tails were removed from all samples to ensure that only trunk or spinal cord cells were collected. 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 [44]. V0v neurons are EGFP-positive cells from Tg(evx1:EGFP) SU1 trunks [11]. Total RNA was extracted using an RNeasy Micro Kit (Qiagen, 74004). The quality of RNA was assessed via an Agilent 2100 Bioanalyser (RNA 6000 Pico Kit, Agilent, 5067–1513) before being converted to fluorescently-labeled cDNA (Ovation Pico WTA System V2, Pico, 3302) and hybridized to a custom-designed Agilent microarray (Agilent #027382). Data pre-processing and normalization was performed using Bioconductor software (https://www.bioconductor.org/). A three-class ANOVA analysis was performed using GEPAS software [45, 46]. Relative expression levels were subjected to a Z-transformation normalization and are presented as Z scores where mean = 0 and standard deviation = +1 (red) to -1 (blue) [47–49]. All reported statistics were corrected for multiple testing [50].
To generate the custom-designed Agilent microarray (Agilent #027382) we first performed comprehensive bioinformatic searches for proteins that contain at least one of the 483 InterPro domains identified in [51] as being specific to transcriptional regulators. These domains comprise three functional classes: DNA binding, chromatin remodeling and general transcription machinery. We identified 3192 potential transcription factors. 2644 of these proteins were identified in Zv8 (Ensembl release 54) of the zebrafish genome and a further 548 non-overlapping transcription factors were identified in the zebrafish Unigene dataset (release 117). Our custom arrays contain 33784 probes corresponding to eight distinct 60-mer probes for each of the transcripts associated with these 3192 proteins. We also included 170 housekeeping genes (five copies of eight probes each), 23 positive controls, such as neurotransmitter markers (two copies of eight probes each) and 49 negative controls (Arabidopsis sequences; multiple copies of eight probes each) on the arrays. Four biological replicates were performed per sample type. Microarray data are deposited at NCBI GEO entry number GSE83723.
in situ hybridization
Embryos were fixed in 4 % paraformaldehyde and single and double in situ hybridization experiments were performed as previously described [52, 53]. Probes for in situ hybridization experiments were prepared using the following templates: evx1 [30], evx2 [29], lbx1a [54] and lmx1ba [24]. A probe for lmx1bb was generated from cDNA as previously described [11, 43] with the following primers: forward CTGGATATCAAGCCGGAGAA; reverse AATTAACCCTCACTAAAGGGATCCGAACATCACATTTCAACA. The lmx1bb probe sequence was selected to avoid cross-hybridization with lmx1ba and other lmx1 family members.
To try and improve signal strength of the lmx1ba probe, we also hydrolyzed the full length lmx1ba probe described above [24] to approximately 200 bp fragments as outlined in [55] and tested two additional lmx1ba probes. The second probe was synthesized from a plasmid containing the last 584 bp of the coding sequence of lmx1ba. The third probe, which recognizes the 3’ coding sequence and UTR of lmx1ba, was generated from cDNA, as previously described [11, 43], with the following primers: forward CGCATGCGTTGGTATCTATG; reverse AATTAACCCTCACTAAAGGGAAAGCATCCTCCACAATGTCC. As these probes did not improve the signal quality when compared to the first probe described above [24], results from these in situ hybridization experiments are not included in this paper.
To determine neurotransmitter phenotypes, we used in situ probes for genes that function as transporters of neurotransmitters or that synthesize specific neurotransmitters as these are some of the most specific molecular markers of these cell fates (e.g. see [56] and references therein). A mixture of probes to slc17a6a and slc17a6b, which encode glutamate transporters, was used to label glutamatergic neurons [56, 57]. To label inhibitory cells we used slc32a1, which encodes a vesicular inhibitory amino acid transporter [33]. To label glycinergic cells a mixture of probes (glyt2a and glyt2b) for the gene slc6a5 were used [56, 57]. The slc6a5 gene encodes for a glycine transporter necessary for glycine reuptake and transport across the plasma membrane. GABAergic neurons were labeled by a mixture of probes to gad1b and gad2 genes (probes previously called gad67a, gad67b and gad65) [56, 57]. The gad1b and gad2 genes encode for glutamic acid decarboxylases, which are necessary for the synthesis of GABA from glutamate.
Immunohistochemistry
Embryos were fixed in 4 % paraformaldehyde and stored in PBS with 0.1 % tween20. To permeabilize embryos they were treated with acetone at -20 °C for 30 min (36 h or younger), 1 h (48 h) or 3 h (7 dpf) and then processed as previously described [11]. Primary antibodies were: mouse anti-GFP (Roche Applied Science, 11814460001, 1:500), rabbit anti-DsRed (Clontech, 632496, 1:200) or rabbit anti-activated caspase-3 (Fisher Scientific/BD, BDB559565, 1:500). Secondary antibodies were: Alexa Fluor 568 goat anti-rabbit (Molecular Probes, A11036, 1:500), Alexa Fluor 488 goat anti-mouse (Molecular Probes, A11029, 1:500) or Alexa Fluor 488 goat anti-rabbit (Molecular Probes, A11034, 1:500).
Double stains
Both double in situ hybridization and immunohistochemistry plus in situ hybridization double labeling experiments were performed as in [52].
Acridine orange treatment
A stock acridine orange base (Sigma-Aldrich, 235474) solution of 2.5 mg/mL in dimethyl sulfoxide (DMSO) was made and stored at -20 °C until used. At 24 h, 36 h and 48 h acridine orange stock solution was added to embryo medium (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2 · 2H2O and 0.33 mM MgSO4 · 7H2O in water) to make a final concentration of 5 μg/mL. Embryos were bathed in the acridine orange / embryo medium solution in the dark at 28.5 °C for 28 min. Embryos were then washed five times in embryo medium for 5 min each and analyzed using fluorescent microscopy on a Zeiss Axio Imager M1 compound microscope and Olympus SZX16 dissecting microscope.
Imaging
Embryos were mounted in 70 % glycerol, 30 % PBS and differential interference contrast (DIC) pictures were taken using an AxioCam MRc5 camera mounted on a Zeiss Axio Imager M1 compound microscope. Embryos from acridine orange experiments and anti-activated caspase-3 experiments were mounted in 2 % 1,4-diazabicyclo[2.2.2] octane (DABCO) and imaged in the same way. Zeiss LSM 710 and LSM 780 confocal microscopes were used to image embryos mounted in DABCO from fluorescent in situ hybridization and immunohistochemistry experiments. Images were processed using Adobe Photoshop software (Adobe, Inc), GNU Image Manipulations Program (GIMP 2.6.10, http://gimp.org) and Image J software (Abramoff et al. [58]). 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
For acridine orange staining and activated caspase-3 immunohistochemistry experiments, cells were counted along both sides of the entire rostral-caudal axis of the spinal cord. For all other experiments, we identified somites 6–10 in each embryo and counted the number of labeled cells in that stretch of the spinal cord. In all cases, embryos were mounted laterally with the somite boundaries on each side of the embryo exactly aligned and the apex of the somite over the middle of the notochord. This ensures that the spinal cord is straight along its dorsal-ventral axis and that cells in the same dorsal/ventral position on opposite sides of the spinal cord will be directly above and below each other. Cell counts for fluorescently-labeled cells were performed by analyzing all focal planes in a confocal stack of the appropriate region(s) of the spinal cord. Labeled cells in embryos analyzed by DIC were counted while examining embryos on the Zeiss Axio Imager M1 compound microscope. We adjusted the focal plane as we examined the embryo to count cells at all medial/lateral positions (both sides of the spinal cord; also see [7, 11, 52, 59]). Values are reported as the mean +/- the standard error of the mean. Results were analyzed using the student’s t-test.
Results
lmx1ba and lmx1bb are expressed by zebrafish dI5 and V0v neurons
lmx1b expression in zebrafish spinal cord. a Three-class ANOVA comparison of V0v cells (class 3), trunk cells (class 1) and all post-mitotic neurons (class 2). p values test hypothesis that there is no differential expression among the 3 classes. Columns represent individual microarray experiments. Rows indicate relative expression levels as normalized data, subjected to a Z-transformation where mean = 0 and standard deviation = 1, where red = normalized expression value of +1 and blue = normalized expression value of -1 (see Methods for more details). lmx1ba and lmx1bb are expressed by V0v neurons. Positive control evx1 is also expressed by V0v neurons. Negative controls eng1b and myod1 are expressed by other neurons (V1 cells) and trunk cells respectively. βactin is a housekeeping gene that is expressed by all populations. b-l Lateral views of zebrafish spinal cord at 27 h (b and c), 30 h (h-l), 36 h (d and e) and 48 h (f and g). Anterior left, dorsal top. in situ hybridization for lmx1ba (b, d and f) and lmx1bb (c, e and g). Black dashed line (b-g) is just below ventral limit of spinal cord, floor plate is right above this, in the most ventral part of the spinal cord, roof plate is the most dorsal part of the spinal cord. Double in situ hybridization for lmx1bb (red) and lbx1a (green) in WT embryo, merged view (h) and magnified single confocal plane of white dotted box region (h’-h’”). in situ hybridization for lmx1bb (red) and EGFP immunohistochemistry (green) in Tg(evx1:EGFP) SU1 embryo, merged view (i) and magnified single confocal plane of white dotted box region (i’-i’”). Double in situ hybridization for lmx1bb (red) and slc17a6 (green) in WT embryo, merged view (j) and magnified single confocal plane of white dotted box region (j’-j’”). in situ hybridization for lmx1bb (red) and EGFP immunohistochemistry (green) in Tg(slc17a6:EGFP) embryo, merged image (k) and magnified single confocal plane of white dotted box region (k’-k’”). White dashed line (k) marks the dorsal limit of the spinal cord. Red staining above the dashed line is outside the spinal cord. Double in situ hybridization for lmx1bb (red) and slc32a1 (green) in WT embryo, merged image (l) and magnified single confocal plane of white dotted box region (l’-l’”). In all cases (h-l) * indicates co-labeled cell, x indicates single labeled lmx1bb-expressing cell. In all cases, at least two independent double-labeling experiments were conducted (h-l). Results were similar for each replicate. Numbers of single and double-labeled cells and number of embryos counted are provided in Tables 1 and 2. Scale bar = 50 μm (b-g), 70 μm (h-l) and 20 μm (h’-l’”)
lmx1ba and lmx1bb are expressed in zebrafish spinal cord
lmx1ba-expressing cells | lmx1bb-expressing cells | |||||||
|---|---|---|---|---|---|---|---|---|
27 h | 30 h | 36 h | 48 h | 27 h | 30 h | 36 h | 48 h | |
Mean | 3.5 | 8.6 | 11.8 | 22.5 | 31.1 | 37.6 | 57 | 80.4 |
SEM | 1 | 1.8 | 0.5 | 2.4 | 1.4 | 1.7 | 1.3 | 1.9 |
n | 4 | 5 | 4 | 4 | 11 | 15 | 4 | 5 |
In contrast, at 27 h, lmx1bb spinal cord expression already extends along the entire rostral-caudal axis in a narrow dorsal-ventral domain (Fig. 1c). Like lmx1ba, lmx1bb is also expressed in the roof plate and floor plate at this stage. As development progresses, more spinal cord neurons express lmx1bb and roof plate expression becomes more prominent while floor plate expression is lost by 36 h (Fig. 1c, e and g; Table 1). By 48 h, lmx1ba and lmx1bb are expressed in presumably overlapping domains, although, as all lmx1ba in situ probes tested produced very weak staining (see methods), it was not possible to confirm this with co-labeling experiments.
Co-expression of other genes with lmx1bb
lmx1bb + Tg(slc17a6:EGFP) double labeling experiments | |||
30 h | lmx1bb | Tg(slc17a6:EGFP) | co-labeled |
Mean | 30 | 105.7 | 21 |
SEM | 3 | 9.7 | 2.2 |
n | 7 | 7 | 7 |
% | 70 % | 20 % | n/a |
lmx1bb + slc17a6 double labeling experiments | |||
30 h | lmx1bb | slc17a6 | co-labeled |
Mean | 32.5 | 99.8 | 25.8 |
SEM | 1.1 | 4.4 | 1.3 |
n | 4 | 4 | 4 |
% | 79 % | 26 % | n/a |
lmx1bb + slc32a1 double labeling experiments | |||
30 h | lmx1bb | slc32a1 | co-labeled |
Mean | 28.3 | 142.7 | 3 |
SEM | 1.2 | 2.4 | 0.6 |
n | 6 | 6 | 6 |
% | 10 % | 2 % | n/a |
lmx1bb + Tg(evx1:EGFP) double labeling experiments | |||
30 h | lmx1bb | Tg(evx1:EGFP) | co-labeled |
Mean | 36 | 70.5 | 13.5 |
SEM | 2.1 | 2.4 | 1.7 |
n | 6 | 6 | 6 |
% | 38 % | 19 % | n/a |
lmx1bb + lbx1a double labeling experiments | |||
30 h | lmx1bb | lbx1a | co-labeled |
Mean | 29.4 | 40 | 13.3 |
SEM | 1.7 | 2.7 | 1 |
n | 7 | 7 | 7 |
% | 45 % | 33 % | n/a |
As mentioned above, our expression profiling of V0v neurons suggested that zebrafish lmx1b genes may also be expressed by these cells (Fig. 1a). To confirm these results we performed EGFP immunohistochemistry and lmx1bb in situ hybridization in Tg(evx1:EGFP) SU1 embryos that express EGFP in V0v neurons [11]. These experiments showed that at 30 h at least 38 % of lmx1bb-expressing neurons are V0v neurons (Fig. 1i; Table 2).
Both V0v cells and dI5 cells are glutamatergic [8, 11, 16, 33, 34]. Moreover, Lmx1b-expressing neurons are glutamatergic in the amniote spinal cord [8, 16, 32]. Therefore, to further confirm the identity of zebrafish lmx1bb-expressing spinal neurons we performed double-labeling experiments. Double in situ hybridization for lmx1bb and glutamatergic markers slc17a6a + slc17a6b (a mixture of probes for both genes, referred to here as slc17a6; see methods), showed that at 30 h at least 79 % of lmx1bb-expressing cells co-express slc17a6 (Fig. 1j; Table 2). To further confirm that most lmx1bb-expressing neurons are glutamatergic, we also performed double staining for EGFP and lmx1bb in 30 h Tg(slc17a6:EGFP) embryos in which many glutamatergic neurons express EGFP [33, 65–67]. In these embryos, we found that approximately 70 % of lmx1bb-expressing neurons also express EGFP (Fig. 1k; Table 2). In contrast, double in situ hybridization with lmx1bb and slc32a1, which labels all spinal cord inhibitory neurons [33, 68], revealed that only 10 % of lmx1bb neurons are inhibitory (Fig. 1l; Table 2). Taken together, these data suggest that the vast majority of zebrafish lmx1bb-expressing cells are glutamatergic and that these glutamatergic cells correspond to dI5 and V0v neurons.
lmx1bb is required for glutamatergic neurotransmitter phenotypes at later developmental stages but does not repress inhibitory neurotransmitter phenotypes
Expression of lmx1b RNAs in lmx1b mutants. Lateral view of zebrafish spinal cord at 48 h (a-f). Anterior left, dorsal top. in situ hybridization of lmx1ba (a-c) or lmx1bb (d-f) in WT (a and d), lmx1ba mutant (b and e) and lmx1bb mutant (c and f). Lower magnification insert in (f) shows expression remaining in hindbrain region. The rest of the head was removed for genotyping. One in situ hybridization of at least 40 embryos was conducted for each of b and e. Two independent in situ hybridizations of at least 50 embryos each were conducted for a, c, d and f. In these cases, results were the same for each replicate experiment. At least three genotyped mutant and wild-type embryos were analyzed in detail for each experiment. Scale bar = 50 μm
lmx1bb is required for glutamatergic phenotypes at later developmental stages but does not repress inhibitory phenotypes. Lateral view of zebrafish spinal cord at 27 h (a, b, h and i), 36 h (c, d, j and k) and 48 h (e-f’, l, m and o-t), anterior left, dorsal top. in situ hybridization for slc17a6a + slc17a6b (slc17a6) (a-f’), slc32a1 (h-m), gad1b + gad2 (GAD) (o, p), slc6a5 (q and r) and pax2a (s and t). (e’ and f’) are magnified views of black dashed box region in (e and f) respectively. Mean number of cells (y-axis) expressing markers slc17a6 (g), slc32a1 (n) and GAD, slc6a5 or pax2a at 48 h (u) in spinal cord region adjacent to somites 6–10 in WT embryos (white) and lmx1bb homozygous mutants (grey) (x-axis). Statistically significant (p < 0.05) comparisons are indicated with square brackets and stars. Error bars indicate standard error of the mean. Two independent experiments were conducted for all slc17a6 and slc32a1 experiments (a-m). Cells count results were similar for each replicate. One experiment was conducted for (o-t). Cell count data presented here (g, n and u) are average values for 4 to 17 embryos from the same in situ hybridization experiment. Precise numbers of embryos counted and p values are provided in Tables 3 and 4. Scale bar = 50 μm (a-f, h-m and o-t) and 25 μm (e’ and f’)
Lmx1bb is required for excitatory and not inhibitory neurotransmitter phenotypes
27 h | 36 h | 48 h | |||||
|---|---|---|---|---|---|---|---|
Marker | WT | lmx1bb -/- | WT | lmx1bb -/- | WT | lmx1bb -/- | |
slc17a6 | Mean | 121.6 | 127.2 | 137.2 | 123.6 | 211 | 175 |
SEM | 5.3 | 4.1 | 2.5 | 2.4 | 5.5 | 2.9 | |
n | 10 | 10 | 13 | 17 | 10 | 13 | |
p value | 0.411 | <0.001 | <0.001 | ||||
slc32a1 | Mean | 149.7 | 148 | 173.2 | 169.7 | 210.5 | 202 |
SEM | 2.7 | 5.8 | 10.9 | 17.1 | 4 | 10.8 | |
n | 10 | 6 | 14 | 7 | 12 | 6 | |
p value | 0.77 | 0.85 | 0.48 | ||||
Expression of genes in WT and lmx1bb mutant embryos
Marker | 48 h | WT | lmx1bb -/- |
|---|---|---|---|
GAD | Mean | 114.2 | 109.5 |
SEM | 1.2 | 4.5 | |
n | 5 | 4 | |
p value | 0.54 | ||
slc6a5 | Mean | 208.3 | 201.2 |
SEM | 4.9 | 6.2 | |
n | 6 | 6 | |
p value | 0.38 | ||
pax2a | Mean | 138 | 140 |
SEM | 3.4 | 3.4 | |
n | 7 | 5 | |
p value | 0.7 | ||
lbx1a | Mean | 96.5 | 96.1 |
SEM | 3.6 | 2.4 | |
n | 8 | 8 | |
p value | 0.93 | ||
evx | Mean | 90.9 | 92.1 |
SEM | 1.9 | 3.9 | |
n | 10 | 6 | |
p value | 0.78 | ||
Tg(evx1:EGFP) SU1 | Mean | 94.9 | 96.3 |
SEM | 2.6 | 0.94 | |
n | 8 | 8 | |
p value | 0.63 | ||
lmx1ba and lmx1bb single mutants and lmx1ba;lmx1bb double mutants have the same spinal cord phenotype
As shown above (Fig. 1d-g), lmx1ba and lmx1bb are expressed in potentially overlapping domains within the zebrafish spinal cord during the developmental time points that we detected neurotransmitter phenotypes in lmx1bb mutants. This suggested that these two ohnologs might function redundantly in the spinal cord. Therefore, we examined spinal cord neurotransmitter phenotypes in lmx1ba single and lmx1ba;lmx1bb double mutants.
Three functional lmx1b alleles are required for correct numbers of glutamatergic cells at later developmental stages. Lateral view of zebrafish spinal cord at 48 h (a-g’, i and j), anterior left, dorsal top. in situ hybridization for slc17a6a + slc17a6b (slc17a6) (a-g’) and slc32a1 (i and j). (a’-g’) are magnified views of black dashed box regions in panels (a-g). Columns on left indicate lmx1ba and lmx1bb genotype. Mean number of cells (y-axis) expressing slc17a6 (h) and slc32a1 (k) in spinal cord region adjacent to somites 6–10 at 48 h (x-axis). Square brackets and star in (h) indicates that each of the first three columns is statistically significantly different from each of the last four columns (p < 0.05). Embryo genotype is indicated below graph. Error bars indicate standard error of the mean. Two independent experiments were conducted for (a-g). Cell count results were similar in each replicate. One experiment was conducted for (i and j). Cell count data presented here (h and k) are average values of 4–13 embryos. Precise numbers of embryos counted and p values are provided in Table 5. Scale bar (g) = 50 μm (a-g) and 20 μm (a’-g’) and scale bar (j) = 50 μm (i, j)
lmx1ba and lmx1bb mutant alleles act redundantly in a dose-dependent manner
48 h | lmx1ba +/+ lmx1bb +/+ | lmx1ba +/- lmx1bb +/+ | lmx1ba +/+ lmx1bb +/- | lmx1ba +/- lmx1bb +/- | lmx1ba -/- lmx1bb +/+ | lmx1ba +/+ lmx1bb -/- | lmx1ba -/- lmx1bb -/- | |
|---|---|---|---|---|---|---|---|---|
slc17a6 | Mean | 246.3 | 250.2 | 233.8 | 196.4 | 199.4 | 204.5 | 193.1 |
SEM | 7.4 | 5.9 | 4.5 | 5.5 | 10.4 | 10 | 6.2 | |
n | 8 | 5 | 4 | 8 | 11 | 9 | 13 | |
p value 1 | n/a | 0.72 | 0.3 | 0.001 | 0.001 | 0.003 | 0.001 | |
p value 2 | 0.3 | 0.08 | n/a | 0.005 | 0.03 | 0.05 | 0.004 | |
p value 3 | 0.001 | 0.001 | 0.005 | n/a | 0.66 | 0.38 | 0.78 | |
p value 4 | 0.003 | 0.004 | 0.05 | 0.38 | 0.7 | n/a | 0.45 | |
p value 5 | 0.001 | 0.001 | 0.03 | 0.66 | n/a | 0.7 | 0.78 | |
slc32a1 | Mean | 179.8 | 188.5 | |||||
SEM | 12.9 | 14.1 | ||||||
n | 6 | 6 | ||||||
p value 1 | n/a | 0.94 |
Given the similarity of the phenotypes in lmx1ba and lmx1bb single and double mutants, we tested whether lmx1bb is required for lmx1ba spinal cord expression or vice versa. However, when we analyzed expression of lmx1ba in lmx1bb mutants and expression of lmx1bb in lmx1ba mutants we saw no obvious differences between WT and mutant embryos (Fig. 2c and e). This suggests that the phenotypic similarities between the mutants are not due to cross-regulation of these two lmx1b genes.
Together, these results suggest that lmx1ba and lmx1bb function partially redundantly in the spinal cord and that the presence of two or more mutant alleles (regardless of whether the mutation is in lmx1ba or lmx1bb) is sufficient to cause a reduction in the number of glutamatergic cells in the spinal cord. To test this, we examined the number of glutamatergic spinal neurons in lmx1ba;lmx1bb double heterozygous embryos and both lmx1ba and lmx1bb single heterozygous embryos. Consistent with our hypothesis, the reduction in the number of glutamatergic neurons in lmx1ba;lmx1bb double heterozygous embryos was not statistically significantly different from the reduction in lmx1ba mutants (p = 0.66), lmx1bb mutants (p = 0.38) or lmx1ba;lmx1bb double mutants (p = 0.78; Fig. 4d, d’ and h; Table 5). In contrast, neither lmx1ba nor lmx1bb single heterozygous embryos had a statistically significant reduction in the number of glutamatergic neurons when compared to WT siblings (p = 0.72 and p = 0.3 respectively; Fig. 4b-c’ and h; Table 5).
To test the possibility that lmx1ba might compensate for the loss of lmx1bb in the repression of inhibitory neurotransmitter phenotypes, we also analyzed the expression of slc32a1 in lmx1ba;lmx1bb double mutants. However, like the lmx1bb single mutant results, the lmx1ba;lmx1bb double mutants had no statistically significant change (p = 0.94) in the number of spinal inhibitory neurons (Fig. 4i-k; Table 5). These data suggest that lmx1ba and lmx1bb are not required to repress (or specify) inhibitory neurotransmitter phenotypes and that the reduction in spinal cord glutamatergic cells in these mutants does not correlate with an increase in inhibitory cells.
The reduction in spinal glutamatergic cells is not due to cell death
To test whether the reduction in glutamatergic neurons might be an indirect consequence of increased cell death, we used both acridine orange (AO) and an activated caspase-3 antibody [76, 77]. As the glutamatergic phenotype is comparable among lmx1ba;lmx1bb double mutants and both single mutants, we used lmx1bb single mutants for these and all subsequent experiments.
There is no increase in apoptosis in lmx1bb mutants between 36 h and 72 h. Lateral view of zebrafish spinal cord at 27 h (a-b’), 36 h (c-d’ and g-h’) and 48 h (e-f’ and i-j’), anterior left, dorsal top. Acridine orange (AO) treatment (a-f’) and activated caspase-3 immunohistochemistry as anterior-posterior montages (g-j’). Sib. in (a) is a sibling embryo to smoothened mutant in (b). (a’-j’) are magnified view of corresponding boxed region. Mean number of cells (y-axis) with activated caspase-3 staining in WT embryos (white) and lmx1bb homozygous mutants (grey) (x-axis) at indicated developmental times. Error bars indicate standard error of the mean. Two independent experiments were conducted for (c-j). Expression (a-j) and cell count data (k) were similar in each replicate. Data presented in (k) are average values of 5-8 embryos from the same experiment. Precise number of embryos counted and p values are provided in Table 6. Scale bar = 100 μm (a-j) and 80 μm (a’-j’)
Activated caspase-3 is not increased in zebrafish lmx1bb mutants during the first 72 h of development
36 h | 48 h | 72 h | |||||
|---|---|---|---|---|---|---|---|
Marker | WT | lmx1bb -/- | WT | lmx1bb -/- | WT | lmx1bb -/- | |
Activated Caspase-3 | Mean | 1.1 | 1.8 | 2.8 | 1.8 | 0.4 | 1 |
SEM | 0.77 | 0.94 | 0.95 | 0.58 | 0.24 | 0.68 | |
n | 8 | 7 | 6 | 5 | 6 | 5 | |
p value | 0.63 | 0.4 | 0.46 | ||||
V0v and dI5 cells form in normal numbers in lmx1bb mutants
lmx1bb is expressed by dI5 and V0v neurons and these cells form in normal numbers in lmx1bb mutants. Lateral view of zebrafish spinal cord at 48 h (a-f), anterior left, dorsal top. in situ hybridization for lbx1a in WT (a) and lmx1bb homozygous mutant (b) embryos. in situ hybridization for evx1 + evx2 (evx) in WT (c) and lmx1bb homozygous mutant (d) embryos. Immunohistochemistry for EGFP in Tg(evx1:EGFP) SU1 WT (e) and lmx1bb homozygous mutant (f) embryos. g Mean number of cells (y-axis) in WT embryos (white) and lmx1bb homozygous mutants (grey) expressing lbx1a, evx or Tg(evx1:EGFP) SU1 in spinal cord region adjacent to somites 6–10 (x-axis). Two independent experiments were conducted for (a-f). Cell count results were similar in each replicate. Data shown here (g) are average values of 6–10 embryos from the same experiment. Precise number of embryos counted and p values are provided in Table 4. Error bars indicate standard error of the mean. Scale bar = 50 μm (a-f)
lmx1bb is required for the glutamatergic phenotype of at least a subset of V0v interneurons
lmx1bb is required for V0v interneuron glutamatergic fates at later stages of development. Lateral view of zebrafish spinal cord at 48 h (a and b) and 7 dpf (c and d), anterior left, dorsal top. in situ hybridization slc17a6 (red) and EGFP immunohistochemistry (green) in WT (a) and lmx1bb mutant (b) Tg(evx1:EGFP) SU1 embryos. Single magnified confocal plane from white dashed box region (a’-a’” and b’-b’”). Immunohistochemistry for EGFP (green) and DsRed (red) in WT (c) and lmx1bb mutant (d) Tg(slc17a6b(vglut2a):loxP-DsRed-loxP-GFP) nns14 ;Tg(evx1:EGFP) SU1 embryos. Single magnified confocal plane from white dashed box region (c’-c’” and d’-d’”). * indicates co-labeled cell, x indicates single labeled EGFP-expressing (V0v) cell. (e and f) Mean number of cells (y-axis) expressing slc17a6 or DsRed (glutamatergic), EGFP (V0v) and slc17a6 or DsRed + EGFP (co-labeled) (x-axis) in WT (white) and lmx1bb homozygous mutants (grey). Error bars indicate standard error of the mean. Three independent experiments were conducted for (c and d). Cell count results were similar in each replicate. One experiment was conducted for (a and b). Data shown here (e and f) are average values of 5–12 embryos. Precise number of embryos counted and p values are provided in Table 7. The glutamatergic and V0v numbers include co-labeled cells. Statistically significant (p < 0.05) comparisons are indicated with square brackets and stars. Scale bar = 30 μm (a-d) and 25 μm (a’-d’”)
Lmx1b regulates the glutamatergic phenotype of a subset of V0v neurons
slc17a6 (glutamatergic) | EGFP (V0v) | co-labeled | ||||
48 h | WT | lmx1bb -/- | WT | lmx1bb -/- | WT | lmx1bb -/- |
Mean | 155 | 142.2 | 92.3 | 93.4 | 44.3 | 32.4 |
SEM | 1.54 | 2.22 | 2.72 | 2.58 | 0.9 | 1.9 |
n | 6 | 5 | 6 | 5 | 6 | 5 |
p value | <0.001 | 0.76 | <0.001 | |||
DsRed (glutamatergic) | EGFP (V0v) | co-labeled | ||||
7 dpf | WT | lmx1bb -/- | WT | lmx1bb -/- | WT | lmx1bb -/- |
Mean | 132.6 | 105.3 | 100.1 | 91.3 | 42.9 | 21.8 |
SEM | 5.5 | 5.2 | 1.4 | 2.1 | 4.9 | 2.7 |
n | 9 | 12 | 9 | 12 | 9 | 12 |
p value | <0.005 | 0.002 | <0.001 | |||
We were also interested in establishing if the reduction of glutamatergic cells in general and/or the reduction in the number of glutamatergic V0v cells persists at later stages of development. As our slc17a6 RNA probe does not label cells effectively in double staining experiments at later stages of development, we created fish transgenic for both Tg(slc17a6b:loxP-DsRed-loxP-GFP) nns14 and Tg(evx1:EGFP) SU1 and heterozygous for the lmx1bb mutation. Embryos from these parents express DsRed in glutamatergic neurons and EGFP in V0v neurons [11, 85, 86]. We then counted single and double-labeled cells in WT and lmx1bb mutant embryos at 7 dpf to determine if the number of excitatory cells and excitatory V0v neurons in particular were reduced in the lmx1bb mutants. We again observed a statistically significant (p < 0.005) reduction in the total number of glutamatergic, DsRed-labeled neurons (Fig. 7d and f; Table 7) as well as a statistically significant reduction in the number of glutamatergic DsRed-labeled V0v neurons (p < 0.001; Fig. 7d and f; Table 7) in lmx1bb mutant embryos compared to WT siblings. More surprisingly, we also observed a very small but statistically significant (p = 0.002) reduction in the number of V0v (EGFP-labeled) neurons at 7 dpf in lmx1bb mutants (Fig. 7d’” and f; Table 7). However, this slight reduction was substantially less than the reduction in double-labeled glutamatergic V0v neurons, demonstrating that the second result cannot be explained by the first. Taken together, these results suggest that lmx1bb is required for the glutamatergic neurotransmitter phenotype of at least a subset of V0v neurons at later developmental stages.
lmx1ba and lmx1bb expression requires evx1 and evx2
Evx1 and Evx2 are required for lmx1b expression. Lateral view of zebrafish spinal cord at 30 h (a-d), anterior, left dorsal top. in situ hybridization for lmx1ba (a, b) and lmx1bb (c, d). Brown coloration in (a-d) is pigment from melanocytes. e Mean number of cells (y-axis) that express lmx1ba and lmx1bb in WT (white) and evx1;evx2 double mutants (grey) (x-axis). Error bars indicate standard error of the mean. Statistically significant (p < 0.05) comparisons are indicated with square brackets and stars. One experiment was conducted for (a and b). Two independent experiments were conducted for (c and d) and cell count results were similar in each replicate. Data presented here (e) are average values for 5–12 embryos. Precise number of embryos counted and p values are provided in Table 8. f. Proposed mechanism for how the excitatory (glutamatergic) neurotransmitter phenotype of at least a subset of V0v neurons is specified and/or maintained. We previously demonstrated that Evx1 & Evx2 specifies the excitatory (glutamatergic) neurotransmitter phenotype and represses inhibitory (glycinergic) phenotypes in V0v cells [11]. The current study demonstrates that Evx1 & Evx2 are also required for lmx1ba and lmx1bb expression. Furthermore, we show that Lmx1bb is required at later developmental stages either to maintain the excitatory (glutamatergic) neurotransmitter phenotype for at least a subset of V0v neurons or to specify the glutamatergic phenotype of a later-forming subset of V0v cells. Scale bar = 50 μm
Evx1 and Evx2 are required for lmx1b expression
Marker | 30 h | WT | evx1 -/- ;evx2 -/- |
|---|---|---|---|
lmx1ba | Mean | 8.6 | 4.8 |
SEM | 0.81 | 1.1 | |
n | 5 | 6 | |
p value | 0.029 | ||
lmx1bb | Mean | 35.9 | 19.7 |
SEM | 6.4 | 5.3 | |
n | 12 | 6 | |
p value | <0.001 | ||
Discussion
In this paper we demonstrate that lmx1ba and lmx1bb (zebrafish ohnologs of Lmx1b) are both expressed by spinal cord interneurons during development (Figs. 1, 2 and 8). This is consistent with a previous report where lmx1bb was shown to be expressed by anterior spinal neurons [24], but in contrast, this earlier study suggested that lmx1ba was not expressed in the zebrafish spinal cord [24]. Given that we have only detected very weak spinal cord staining with our in situ hybridizations for lmx1ba, despite trying three different RNA in situ probes (see methods), we think that the spinal cord staining of lmx1ba was too weak to be easily detected in these previous experiments.
Consistent with results in mouse [16, 18, 32, 62–64, 87], we demonstrate here that in the zebrafish spinal cord lmx1bb is predominantly expressed by glutamatergic neurons (Fig. 1j and k). Our results demonstrate that at least 79 % of lmx1bb-expressing neurons are glutamatergic (Fig. 1j and k; Table 2). We consistently see fewer labeled cells with the Tg(slc17a6:EGFP) and Tg(slc17a6b:loxP-DsRed-loxP-GFP) nns14 lines than we do with single in situ hybridization for slc17a6 (Fig. 1j and k; Table 2; data not shown), suggesting these transgenic lines either do not label all spinal cord glutamatergic cells at the stages examined or there is a delay in the expression of the fluorescent proteins compared to slc17a6 RNA. slc17a6 is also not a strong probe in double-labeling experiments and as a result it labels slightly fewer cells in double in situ hybridizations than in single in situ hybridizations. This suggests that some of the lmx1bb-positive, slc17a6-negative cells in our double labels may also be glutamatergic. Therefore, it is possible that more than 79 % of lmx1bb-expressing neurons are glutamatergic, especially as we also show that only about 10 % of lmx1bb-expressing cells are inhibitory (Fig. 1l; Table 2).
Also consistent with results in amniotes [16, 18, 32, 62–64, 87], our analyses suggest that a subset of lmx1b-expressing spinal cord cells are dI5 cells. dI5 cells constitute about a third of all Lbx1-expressing spinal cord cells and they are also the only excitatory, Lbx1-expressing spinal cells [8–10, 16]. We find that at least 45 % of lmx1bb-expressing spinal cells co-express lbx1a and that these co-expressing cells constitute about a third of the lbx1a-expressing cells (Fig. 1h; Table 2). Together with the fact that most lmx1bb-expressing cells are glutamatergic, this suggests that at least most of the cells that co-express lmx1bb and lbx1a are dI5 cells.
However, in contrast to previous reports in amniotes, we also find that a substantial proportion of lmx1bb-expressing spinal cells (at least 38 %) are V0v neurons (Fig. 1i; Table 2). This is the first time that lmx1bb expression has been described in this cell type in any animal. However, a small subset of Lmx1b cells are located in the ventral spinal cord of E12.5 mice [32] in a region similar to where Evx1, a V0v marker, is expressed [12]. Therefore, it is possible that some mouse V0v neurons may also express Lmx1b. Interestingly, at the stages that we examined, only a subset of zebrafish V0v neurons express lmx1bb (Fig. 1i; Table 2). This suggests that lmx1bb may be expressed by a specific subset of V0v interneurons. Interestingly, Satou and colleagues have shown that V0v neurons can be divided into three subsets based on their morphology [31]. Alternatively, it is possible that all V0v cells express lmx1bb, but either only transiently, or only at later stages of their development, resulting in only a subset being co-labeled at any particular time.
While we were unable to successfully perform double-labeling experiments with lmx1ba and lmx1bb, due to very weak staining with our lmx1ba RNA probes, our results suggest that these two genes are expressed by the same spinal cord neurons. Firstly, their spinal cord expression patterns are very similar, although lmx1ba is expressed later than lmx1bb in most spinal cord domains (Figs. 1d-g, 2a-f and 8a and c). Secondly, lmx1ba single mutant, lmx1bb single mutant, lmx1ba;lmx1bb double mutant and lmx1ba;lmx1bb double heterozygous embryos all have the exact same spinal cord phenotype (the same reduction in the number of glutamatergic neurons) suggesting that the two ohnologs have redundant functions in the spinal cord, and must, therefore, be co-expressed in at least some cells (Fig. 4; Table 5). Interestingly, other studies using zebrafish to examine lmx1ba and lmx1bb functions in the isthmus, diencephalon and eye have also found that these two genes have overlapping expression and function redundantly in those tissues [21, 23, 24]. Finally, in evx1;evx2 double mutants the number of cells expressing either lmx1ba or lmx1bb is reduced (Fig. 8; Table 8). Given that evx1 and evx2 are only expressed in V0v neurons in the spinal cord, this strongly suggests that both of the lmx1b ohnologs are expressed in V0v cells.
In this paper, we demonstrate that Lmx1b transcription factors are required in zebrafish for correct numbers of spinal glutamatergic cells at later stages of development. At 27 h there is no change in the number of spinal glutamatergic cells in the spinal cord in lmx1bb mutants, but by 36 h there is a reduction in the number of glutamatergic cells and this reduction becomes more severe by 48 h. We also demonstrate that this phenotype persists until at least 7 days. In contrast, lmx1b-expressing dI5 and V0v cells are still present in normal numbers at 48 h and we observe no increase in cell death at either 36 h or 48 h. These results suggest that the reduction of glutamatergic cells is not a consequence of cell death, changes in cell proliferation or global changes in cell fate.
Prior to this study, Lmx1b had been implicated in correct neuronal migration, connectivity and viability [18, 19, 25, 62, 88, 89]. However, while data from previous studies suggested that Lmx1b may have a role in the development of a subset of spinal cord glutamatergic cells (e.g. [8, 9, 18, 62]), a precise function in glutamatergic fate specification or maintenance had not been identified. Interestingly, when Lmx1b was conditionally ablated specifically in the mouse spinal cord, at E18.5 there was also a reduction in the number of glutamatergic neurons and no change in the number of inhibitory neurons, which is similar to our results in zebrafish (Figs. 3 and 4) [18]. However, the authors of this study attributed this reduction in glutamatergic neurons to cell death because they also observed a reduction in the total number of cells in the dorsal horn and an increase in caspase-3-positive neurons. Despite this, the authors speculated that Lmx1b may function in the maintenance of the neurotransmitter phenotype prior to the death of these cells [18].
Our results differ from this mouse study because we do not see any evidence of an increase in cell death in the spinal cord of zebrafish lmx1bb mutants, at least between 36 and 72 h, even though we see statistically significant reductions in spinal glutamatergic cells at these stages. It is possible that the slight reduction in EGFP-labeled V0v cells that we see at 7 dpf might be due to cell death, but if this is the case this is likely to be a second, later phenotype as it occurs much later than the reduction in glutamatergic cells and affects a much smaller number of cells than the glutamatergic phenotype. Future studies could perform cell death assays to test whether V0v cells die at later stages, but this would not be trivial as this death could occur any time between 72 h and 7 dpf and the number of cells lost at 7 dpf is very small.
The lack of cell death in zebrafish lmx1bb mutants at earlier developmental stages where we see neurotransmitter phenotypes might seem surprising given the evidence that mis-programmed cells die in the mammalian spinal cord [18, 88, 89]. This result could possibly be the consequence of different developmental strategies being utilized in mouse and zebrafish spinal cords. In zebrafish embryos, with the exception of Rohon Beard cells, there is very little apoptosis in the spinal cord during development [90], compared to substantially more programed cell death in the mouse spinal cord [91–93]. It is possible that the fast development and/or smaller size of zebrafish embryos makes it more difficult to utilize a strategy of creating and pruning excess neurons. In this case, mice might be better equipped to eliminate neurons with incorrect functional characteristics than zebrafish. Consistent with this hypothesis, in an earlier study when we removed Pax2 and Pax8 function in zebrafish, several subsets of spinal interneurons lost their inhibitory neurotransmitter phenotypes but they were still present in normal numbers and had normal morphologies and axonal trajectories [7], suggesting that their viability was not affected. Similarly, in zebrafish evx1;evx2 double mutants, V0v spinal neurons switch their neurotransmitter phenotypes from excitatory to inhibitory but they retain normal axon projections until at least 48 h, again suggesting that their viability is not affected [11].
Regardless of the reason, we do not see any evidence of increased cell death in the spinal cord of zebrafish lmx1bb mutants between 36 and 72 h. We also do not see any change in the numbers of V0v or dI5 cells, suggesting that the cells that usually express lmx1bb still form and are present in normal numbers (Fig. 6; Table 4). However, in contrast and as discussed above, our data suggest that Lmx1b transcription factors are required either to maintain the glutamatergic neurotransmitter phenotype of a subset of excitatory spinal neurons or to specify the glutamatergic phenotype of a later-forming subset of spinal neurons (Figs. 3 and 4). If Lmx1bb is required to maintain a subset of glutamatergic fates in the spinal cord, this would be consistent with Lmx1b function in some other regions of the CNS, specifically the mouse raphe nucleus and trigeminal brainstem complex, where Lmx1b is required to maintain specific neurotransmitter phenotypes [19, 25, 94, 95]. However, in the case of the raphe nucleus it is a serotonergic phenotype rather than a glutamatergic phenotype that Lmx1b maintains.
As lmx1bb has not previously been shown to be expressed by V0v neurons, we were particularly interested in testing whether lmx1bb is specifically required for the glutamatergic phenotypes of these cells. We found that there is a statistically significant reduction in the number of glutamatergic V0v cells in lmx1bb mutants at both 48 h and 7dpf (Fig. 7; Table 7). Despite the fact that previous reports have shown that all V0v cells are excitatory [11, 31], not all WT V0v neurons were co-labeled with slc17a6 or Tg(slc17a6b(vglut2a):loxP-DsRed-loxP-GFP) nns14 in these experiments. This is probably because, as discussed above, slc17a6 is a weak probe in double labelling experiments and the Tg(slc17a6b(vglut2a):loxP-DsRed-loxP-GFP) nns14 transgenic line only labels a subset of glutamatergic spinal cord cells in our hands. In contrast, we observed the same number of evx-expressing cells by in situ hybridization as EGFP-positive cells in Tg(evx1:EGFP) SU1 embryos (Fig. 6c-g; Table 4), suggesting that this transgenic line labels all V0v cells, even at later stages of development, as has previously been shown for earlier stages [11].
In addition to showing that lmx1bb is expressed by V0v cells, our results also demonstrate that lmx1ba and lmx1bb expression in these cells requires Evx1 and Evx2 activity. When we examine lmx1ba and lmx1bb expression in evx1;evx2 double mutant embryos we see a statistically significant reduction in the number of lmx1ba and lmx1bb-expressing neurons (Fig. 8a-e; Table 8). Our previous work demonstrated that Evx1 and Evx2 function partially redundantly to specify the glutamatergic neurotransmitter phenotype of V0v neurons [11]. Combined with these earlier results, the data in this paper start to elucidate a pathway of neurotransmitter fate specification and maintenance for V0v cells, with Evx1 and Evx2 specifying the glutamatergic neurotransmitter phenotype as well as lmx1ba and lmx1bb expression. The lmx1b genes then function downstream of Evx1 and Evx2, either to maintain the glutamatergic neurotransmitter phenotype of at least a subset of V0v neurons or to specify the glutamatergic neurotransmitter phenotype of a late-forming subset of V0v cells (Fig. 8f).
Interestingly, our results also show that correct Lmx1b function requires three functional lmx1b alleles in zebrafish, but it does not seem to matter which lmx1b alleles these are (Fig. 4). This suggests that lmx1ba and lmx1bb are at least partially redundant. Given that lmx1ba and lmx1bb are ohnologs that presumably arose in the teleost specific genome duplication event [96, 97], this requirement for three functional lmx1b alleles must have arisen in the teleost lineage. Interestingly, in humans, just one mutant allele of Lmx1b causes the autosomal dominant disorder nail-patella syndrome (NPS), suggesting that gene dosage is also important in mammals [98, 99]. Moreover, a quarter of NPS patients experience peripheral neurological symptoms which may be the result of improper specification of spinal cord neurons [18, 62, 100], suggesting that our results may have direct relevance to this human disorder.
During these studies we also discovered that slc32a1, previously believed to be expressed by all inhibitory neurons, does not label all inhibitory spinal neurons at 48 h in zebrafish [33, 68]. At 48 h slc32a1 expression is restricted to a band of neurons in the middle of the dorsal-ventral axis of the spinal cord (Fig. 3l), whereas the GADs (markers of GABAergic cells) are also expressed in more ventral regions at this stage (Fig. 3l and o). At 27 h slc32a1 seems to be expressed by all inhibitory neurons, as previously reported [33, 68]. However, at 36 h slc32a1 expression starts to be restricted to more dorsal inhibitory populations, although some ventral slc32a1-expressing cells are still detected. It is this sporadic ventral slc32a1 expression that likely causes the larger variations in the number of slc32a1-expressing cells detected at 36 h when compared to 27 h and 48 h (Fig. 3j and k; Table 3).
Conclusions
In conclusion, we demonstrate that lmx1ba and lmx1bb are expressed by V0v and dI5 spinal interneurons. These genes are required, partially redundantly, in a dose-dependent manner, for the glutamatergic neurotransmitter phenotype of at least a subset of these neurons at later developmental stages. However, lmx1ba and lmx1bb are not required to repress (or specify) inhibitory neurotransmitter phenotypes as there is no statistically significant change in the number of inhibitory cells in either lmx1ba or lmx1bb single or double mutants. We also show that lmx1ba and lmx1bb require Evx1 and Evx2 for their expression in V0v neurons, suggesting that lmx1ba and lmx1bb act downstream of Evx1 and Evx2 in specifying or maintaining the glutamatergic neurotransmitter phenotype of at least a subset of V0v neurons. Taken together, our results provide new and powerful insights into the mechanisms required for excitatory neurotransmitter phenotypes within the spinal cord.
Abbreviations
- AO:
-
Acridine Orange
- CNS:
-
Central Nervous System
- DABCO:
-
1,4-diazabicyclo[2.2.2]octane
- DIC:
-
Differential Interference Contrast
- DMSO:
-
Dimethyl Sulfoxide
- dpf:
-
Days Post Fertilization
- FACS:
-
Fluorescent Activated Cell Sorting
- GADs:
-
Glutamic Acid Decarboxylases
- h:
-
Hours Post Fertilization
- IACUC:
-
Institutional Animal Care and Use Committee
- NPS:
-
Nail-patella Syndrome
- PBS:
-
Phosphate-buffered Saline
- PTU:
-
1-phenyl 2-thiourea
- WT:
-
Wild-type
Declarations
Acknowledgments
We would like to thank Sean Megason for providing us with the lmx1bb jj410 line, 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, labeling, and hybridization, Ian Swinburne, William Haws, Matthew Allen and Christiane Voufo for help with genotyping, Henry Putz, Jessica Bouchard, Annika Swanson, Paul Campbell and several SU undergraduate fish husbandry workers for help with maintaining zebrafish lines and Santanu Banerjee for helpful input during the writing of this manuscript.
Funding
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 and Syracuse University start-up funds awarded to KL as well as NIH R01EY016060 awarded to BL and European IP ZF-Health funds awarded to U.S.
Authors’ contributions
WH performed most of the experiments in the paper, including most of the single and double mutant experiments, many of the double-labeling experiments, cell death experiments, and all of the cell counts; JM created the Tg(evx1:EGFP) SU1 line, created and maintained the evx1;evx2 double mutant line and performed some of the double-labeling experiments; JB and BL generated the lmx1ba mw80 mutant and maintained the lmx1ba;lmx1bb double mutants; SE and SdJ performed cell dissociations and RNA extractions; SE performed the FACS analyses, data pre-processing and initial analyses in GEPAS, curation of data for NCBI and some in situ hybridization experiments; OA and JL performed the bioinformatic searches to identify potential transcription factor genes for the microarrays under the direction and mentorship of US; OA, JL, US, and KL designed the microarrays; KL designed and directed the study; WH and KL wrote most of the paper with input from the other authors. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Ethics approval and consent to participate
All zebrafish experiments in this research were approved either by the UK Home Office or by the Syracuse University IACUC committee.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Deneris ES, Hobert O. Maintenance of postmitotic neuronal cell identity. Nat Neurosci. 2014;17:899–907.View ArticlePubMedPubMed CentralGoogle Scholar
- Spitzer NC. Activity-dependent neurotransmitter respecification. Nat Rev Neurosci. 2012;13:94–106.PubMedPubMed CentralGoogle Scholar
- 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:510–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Pittenger C, Bloch MH, Williams K. Glutamate abnormalities in obsessive compulsive disorder: neurobiology, pathophysiology, and treatment. Pharmacol Ther. 2011;132:314–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Xue M, Atallah BV, Scanziani M. Equalizing excitation–inhibition ratios across visual cortical neurons. Nature. 2014;511:596–600.View ArticlePubMedPubMed CentralGoogle Scholar
- Zadori D, Veres G, Szalardy L, Klivenyi P, Toldi J, Vecsei L. Glutamatergic dysfunctioning in Alzheimer’s disease and related therapeutic targets. J Alzheimer’s Dis. 2014;42:177–87.Google Scholar
- Batista MF, Lewis KE. Pax2/8 act redundantly to specify glycinergic and GABAergic fates of multiple spinal interneurons. Dev Biol. 2008;323:88–97.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:1510–5.View ArticlePubMedGoogle Scholar
- 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:770–8.View ArticlePubMedGoogle Scholar
- Glasgow SM, Henke RM, Macdonald RJ, Wright CVE, Johnson JE. Ptf1a determines GABAergic over glutamatergic neuronal cell fate in the spinal cord dorsal horn. Development. 2005;132:5461–9.View ArticlePubMedGoogle Scholar
- Juárez-Morales JL, Schulte CJ, Pezoa SA, Vallejo GK, Hilinski WC, England SJ, et al. Evx1 and Evx2 specify excitatory neurotransmitter fates and suppress inhibitory fates through a Pax2-independent mechanism. Neural Dev. 2016;11:5.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:385–99.View ArticlePubMedGoogle Scholar
- 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:357–66.View ArticlePubMedGoogle Scholar
- Simon HH, Saueressig H, Wurst W, Goulding MD, O’Leary DD. Fate of midbrain dopaminergic neurons controlled by the engrailed genes. J Neurosci. 2001;21:3126–34.PubMedGoogle Scholar
- Albéri L, Sgadò P, Simon HH. Engrailed genes are cell-autonomously required to prevent apoptosis in mesencephalic dopaminergic neurons. Development. 2004;131:3229–36.View ArticlePubMedGoogle Scholar
- 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:510–7.View ArticlePubMedGoogle Scholar
- 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:394–405.View ArticlePubMedGoogle Scholar
- Szabo NE, da Silva RV, Sotocinal SG, Zeilhofer HU, Mogil JS, Kania A. Hoxb8 intersection defines a role for lmx1b in excitatory dorsal horn neuron development, spinofugal connectivity, and nociception. J. Neurosci. 2015;35:5233–46.View ArticleGoogle Scholar
- Zhao Z-Q, Scott M, Chiechio S, Wang J-S, Renner KJ, Gereau RW, et al. Lmx1b is required for maintenance of central serotonergic neurons and mice lacking central serotonergic system exhibit normal locomotor activity. J Neurosci. 2006;26:12781–8.View ArticlePubMedGoogle Scholar
- Ding Y-Q, Marklund U, Yuan W, Yin J, Wegman L, Ericson J, et al. Lmx1b is essential for the development of serotonergic neurons. Nat Neurosci. 2003;6:933–8.View ArticlePubMedGoogle Scholar
- Filippi A, Dürr K, Ryu S, Willaredt M, Holzschuh J, Driever W. Expression and function of nr4a2, lmx1b, and pitx3 in zebrafish dopaminergic and noradrenergic neuronal development. Dev Biol. 2007;7:135.Google Scholar
- Smidt MP, Asbreuk CH, Cox JJ, Chen H, Johnson RL, Burbach JP. A second independent pathway for development of mesencephalic dopaminergic neurons requires Lmx1b. Nat Neurosci. 2000;3:337–41.View ArticlePubMedGoogle Scholar
- O’Hara FP, Beck E, Barr LK, Wong LL, Kessler DS, Riddle RD. Zebrafish Lmx1b.1 and Lmx1b.2 are required for maintenance of the isthmic organizer. Development. 2005;132:3163–73.View ArticlePubMedPubMed CentralGoogle Scholar
- McMahon C, Gestri G, Wilson SW, Link B a. Lmx1b is essential for survival of periocular mesenchymal cells and influences Fgf-mediated retinal patterning in zebrafish. Dev Biol. 2009;332:287–98.View ArticlePubMedPubMed CentralGoogle Scholar
- Laguna A, Schintu N, Nobre A, Alvarsson A, Volakakis N, Jacobsen JK, et al. Dopaminergic control of autophagic-lysosomal function implicates Lmx1b in Parkinson’s disease. Nat Neurosci. 2015;18:826–35.View ArticlePubMedGoogle Scholar
- Avaron F, Thaëron-Antono C, Beck CW, Borday-Birraux V, Géraudie 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:145–56.View ArticlePubMedGoogle Scholar
- Dollé P, Fraulob V, Duboule D. Developmental expression of the mouse Evx-2 gene: relationship with the evolution of the HOM/Hox complex. Dev Supp. 1994;143–53. http://www.ncbi.nlm.n2ih.gov/pubmed/7579515.
- Griener A, Zhang W, Kao H, Wagner C, Gosgnach S. Probing diversity within subpopulations of locomotor-related V0 interneurons. Dev Neurobiol. 2015;75:1189–203.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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:167–72.View ArticlePubMedGoogle Scholar
- Satou C, Kimura Y, Higashijima S-I. Generation of multiple classes of v0 neurons in zebrafish spinal cord: progenitor heterogeneity and temporal control of neuronal diversity. J Neurosci. 2012;32:1771–83.View ArticlePubMedGoogle Scholar
- Dai JX, Hu ZL, Shi M, Guo C, Ding YQ. Postnatal ontogeny of the transcription factor Lmx1b in the mouse central nervous system. J Comp Neurol. 2008;509:341–55.View ArticlePubMedGoogle Scholar
- 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:5684–97.View ArticlePubMedGoogle Scholar
- Karlsson J, von Hofsten J, Olsson PE. Generating transparent zebrafish: a refined method to improve detection of gene expression during embryonic development. Mar Biotechnol. 2001;3:522–7.View ArticlePubMedGoogle Scholar
- Whittaker JR. An analysis of melanogenesis in differentiating pigment cells of ascidian embryos. Dev Biol. 1966;14:1–39.View ArticlePubMedGoogle Scholar
- Epping JJ. Melanogenesis in amphibians.I.A study of the fine structure of the normal and phenylthiourea treated pigmented epithelium in Rana pipiens tadpole eyes. Cell Tissue Res. 1970;103:238–46.Google Scholar
- Schulte CJ, Allen C, England SJ, Juárez-Morales JL, Lewis KE. Evx1 is required for joint formation in zebrafish fin dermoskeleton. Dev Dyn. 2011;240:1240–8.View ArticlePubMedGoogle Scholar
- Schibler A, Malicki J. A screen for genetic defects of the zebrafish ear. Mech Dev. 2007;124:592–604.View ArticlePubMedGoogle Scholar
- Obholzer N, Swinburne IA, Schwab E, Nechiporuk AV, Nicolson T, Megason SG. Rapid positional cloning of zebrafish mutations by linkage and homozygosity mapping using whole-genome sequencing. Development. 2012;139:4280–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanjana NE, Cong L, Zhou Y, Cunniff MM, Feng G, Zhang F. A transcription activator-like effector toolbox for genome engineering. Nat Protoc. 2012;7:171–92.View ArticlePubMedGoogle Scholar
- Truett G, Heeger P, Mynatt R, Truett A. Preparation of PCR-quality mouse genomic DNA with hot sodium hydroxide and tris (HotSHOT). Biotechniques. 2000;29:29–30.Google Scholar
- Cerda GA, Hargrave M, Lewis KE. RNA profiling of FAC-sorted neurons from the developing zebrafish spinal cord. Dev Dyn. 2009;238:150–61.View ArticlePubMedGoogle Scholar
- England SJ, Hilinski WC, de Jager S, Andrzejczuk L, Campbell P, Chowdhury T, et al. Identifying transcription factors expressed by ventral spinal cord interneurons [Internet]. ZFIN on-line Publ. 2014. Available from: 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:279–93.View ArticlePubMedGoogle Scholar
- Montaner D, Tárraga J, Huerta-Cepas J, Burguet J, Vaquerizas JM, Conde L, et al. Next station in microarray data analysis: GEPAS. Nucleic Acids Res. 2006;34:W486–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Tárraga J, Medina I, Carbonell J, Huerta-Cepas J, Minguez P, Alloza E, et al. GEPAS, a web-based tool for microarray data analysis and interpretation. Nucleic Acids Res. 2008;36:308–14.View ArticleGoogle Scholar
- Tan PK, Downey TJ, Spitznagel EL, Xu P, Fu D, Dimitrov DS, et al. Evaluation of gene expression measurements from commercial microarray platforms. Nucleic Acids Res. 2003;31:5676–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheadle C, Vawter MP, Freed WJ, Becker KG. Analysis of microarray data using Z score transformation. J Mol Diagn. 2003;5:73–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Sehgal M, Gupta R, Moussa A, Singh TR. An integrative approach for mapping differentially expressed genes and network components using novel parameters to elucidate key regulatory genes in colorectal cancer. PLoS One. 2015;10:1–18.Google Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser. 1995;57:289–300.Google Scholar
- Armant O, März M, Schmidt R, Ferg M, Diotel N, Ertzer R, et al. Genome-wide, whole mount in situ analysis of transcriptional regulators in zebrafish embryos. Dev Biol. 2013;380:351–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Batista MF, Jacobstein J, Lewis KE. Zebrafish V2 cells develop into excitatory CiD and Notch signalling dependent inhibitory VeLD interneurons. Dev Biol. 2008;322:263–75.View ArticlePubMedGoogle Scholar
- 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:2835–46.PubMedGoogle Scholar
- Ochi H, Westerfield M. Lbx2 regulates formation of myofibrils. Dev Biol. 2009;9:13.Google Scholar
- Cox KH, DeLeon DV, Angerer LM, Angerer RC. Detection of mRNAs in sea urchin embryos by in situ hybridization using asymmetric RNA probes. Dev Biol. 1984;101:485–502.View ArticlePubMedGoogle Scholar
- Higashijima S-I, Mandel G, Fetcho JR. Distribution of prospective glutamatergic, glycinergic, and GABAergic neurons in embryonic and larval zebrafish. J Comp Neurol. 2004;480:1–18.View ArticlePubMedGoogle Scholar
- Higashijima S-I, Schaefer M, Fetcho JR. Neurotransmitter properties of spinal interneurons in embryonic and larval zebrafish. J Comp Neurol. 2004;480:19–37.View ArticlePubMedGoogle Scholar
- Abràmoff MD, Magalhães PJ, Ram SJ. Image processing with ImageJ. Biophotonics Int. 2005;11:36–43.Google Scholar
- 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:5121–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Z-F, Rebelo S, White F, Malmberg AB, Baba H, Lima D, et al. The paired homeodomain protein DRG11 is required for the projection of cutaneous sensory afferent fibers to the dorsal spinal cord. Neuron. 2001;31:59–73.View ArticlePubMedGoogle Scholar
- Chizhikov VV, Millen KJ. Control of roof plate development and signaling by Lmx1b in the caudal vertebrate CNS. J Neurosci. 2004;24:5694–703.View ArticlePubMedGoogle Scholar
- Ding Y-Q, Yin J, Kania A, Zhao Z-Q, Johnson RL, Chen Z-F. Lmx1b controls the differentiation and migration of the superficial dorsal horn neurons of the spinal cord. Development. 2004;131:3693–703.View ArticlePubMedGoogle Scholar
- Müller T, Brohmann H, Pierani A, Heppenstall PA, Lewin GR, Jessell TM, et al. The homeodomain factor lbx1 distinguishes two major programs of neuronal differentiation in the dorsal spinal cord. Neuron. 2002;34:551–62.View ArticlePubMedGoogle Scholar
- Gross MK, Dottori M, Goulding M. Lbx1 specifies somatosensory association interneurons in the dorsal spinal cord. Neuron. 2002;34:535–49.View ArticlePubMedGoogle Scholar
- Bae Y-K, 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:406–26.View ArticlePubMedGoogle Scholar
- Kani S, Bae YK, Shimizu T, Tanabe K, Satou C, Parsons MJ, et al. Proneural gene-linked neurogenesis in zebrafish cerebellum. Dev Biol. 2010;343:1–17.View ArticlePubMedGoogle Scholar
- Kinkhabwala A, Riley M, Koyama M, Monen J, Satou C, Kimura Y. A structural and functional ground plan for neurons in the hindbrain of zebrafish. PNAS. 2010;108:1164–9.View ArticleGoogle Scholar
- Sagne C, El Mestikawy S, Isambert M-F, Hamon M, Henry J-P, Giros B, et al. Cloning of a functional vesicular GABA and glycine transporter by screening of genome databases. FEBS Lett. 1997;417:177–83.View ArticlePubMedGoogle Scholar
- Erlander MG, Tillakaratne NJK, Feldblum S, Patel N, Tobin AJ. Two genes encode distinct gluamate decarboxylases. Neuron. 1991;7:91–100.View ArticlePubMedGoogle Scholar
- Kaufman DL, Houser CR, Tobin AJ. Two forms of the y-aminobutyric acid synthetic enzyme glutamate decarboxylase have distinct intraneuronal distributions and cofactors interactions. J Neurochem. 1991;56:720–3.View ArticlePubMedGoogle Scholar
- Martin SC, Heinrich G, Sandell JH. Sequence and expression of glutamic acid decarboxylase isoforms in the developing zebrafish. J Comp Neurol. 1998;396:253–66.View ArticlePubMedGoogle Scholar
- Jursky F, Nelson N. Localization of glycine neurotransmitter trasnsporter (GLYT2) reveals correlation with the distribution of glycine receptor. J Neurochem. 1995;66:589–98.Google Scholar
- Luque JM, Nelson N, Richards JG. Cellular expression of glycine transporter 2 messenger RNA exclusively in rat hindbrain and spinal cord. Neuroscience. 1995;64:525–35.View ArticlePubMedGoogle Scholar
- Zafra F, Aragón C, Olivares L, Danbolt NC, Giménez C, Storm-Mathisen J. Glycine transporters are differentially expressed among CNS cells. J Neurosci. 1995;15:3952–69.PubMedGoogle Scholar
- Zafra F, Gomeza J, Olivares L, Aragon C, Gimenez C. Regional distribution and developmental variation of the glycine transporters GLYT1 and GLYT2 in the rat CNS. Eur J Neurosci. 1995;7:1342–52.View ArticlePubMedGoogle Scholar
- Chen W, Burgess S, Hopkins N. Analysis of the zebrafish smoothened mutant reveals conserved and divergent functions of hedgehog activity. Development. 2001;128:2385–96.PubMedGoogle Scholar
- Jette CA, Flanagan A, Ryan J, Pyati UJ, Carbonneau S, Stewart RA, et al. BIM and other BCL-2 family proteins exhibit cross-species conservation of function between zebrafish and mammals. Cell Death Differ. 2008;15:1063–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Barrallo-Gimeno A, Holzschuh J, Driever W, Knapik EW. Neural crest survival and differentiation in zebrafish depends on mont blanc/tfap2a gene function. Development. 2004;131:1463–77.View ArticlePubMedGoogle Scholar
- Berry FB, Skarie JM, Mirzayans F, Fortin Y, Hudson TJ, Raymond V, et al. FOXC1 is required for cell viability and resistance to oxidative stress in the eye through the transcriptional regulation of FOXO1A. Hum Mol Genet. 2008;17:490–505.View ArticleGoogle Scholar
- Chen HL, Yuh CH, Wu KK. Nestin is essential for zebrafish brain and eye development through control of progenitor cell apoptosis. PLoS One. 2010;5:e9318.View ArticlePubMedPubMed CentralGoogle Scholar
- Yeh LK, Liu CY, Chien CL, Converse RL, Kao WWY, Chen MS, et al. Molecular analysis and characterization of zebrafish Keratocan (zKera) gene. J Biol Chem. 2008;283:506–17.View ArticlePubMedGoogle Scholar
- Sidi S, Sanda T, Kennedy RD, Hagen AT, Jette CA, Hoffmans R, et al. Chk1 suppresses a caspase-2 apoptotic response to DNA damage that bypasses p53, Bcl-2, and caspase-3. Cell. 2008;133:864–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorrells S, Carbonneau S, Harrington E, Chen AT, Hast B, Milash B, et al. Ccdc94 protects cells from ionizing radiation by inhibiting the xxpression of p53. PLoS Genet. 2012;8:e1002922.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorrells S, Toruno C, Stewart RA, Jette C. Analysis of apoptosis in zebrafish embryos by whole-mount immunofluorescence to detect activated Caspase 3. J Vis Exp. 2013;20:e51060.Google Scholar
- Miyasaka N, Morimoto K, Tsubokawa T, Higashijima S, Okamoto H, Yoshihara Y. From the olfactory bulb to higher brain centers: genetic visualization of secondary olfactory pathways in zebrafish. J Neurosci. 2009;29:4756–67.View ArticlePubMedGoogle Scholar
- Koyama M, Kinkhabwala A, Satou C, Higashijima S, Fetcho J. Mapping a sensory-motor network onto a structural and functional ground plan in the hindbrain. PNAS. 2011;108:1–6.View ArticleGoogle Scholar
- 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:367–84.View ArticlePubMedGoogle Scholar
- Tiret L, Le Mouellic H, Maury M, Brûlet P. Increased apoptosis of motoneurons and altered somatotopic maps in the brachial spinal cord of Hoxc-8-deficient mice. Development. 1998;125:279–91.PubMedGoogle Scholar
- Ross SE, Mardinly AR, McCord AE, Zurawski J, Cohen S, Jung C, et al. Loss of inhibitory interneurons in the dorsal spinal cord and elevated Itch in Bhlhb5 mutant mice. Neuron. 2010;65:886–98.View ArticlePubMedPubMed CentralGoogle Scholar
- van Ham TJ, Mapes J, Kokel D, Peterson RT. Live imaging of apoptotic cells in zebrafish. FASEB J. 2010;24:4336–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamamoto Y, Henderson CE. Patterns of programmed cell death in populations of developing spinal motoneurons in chicken, mouse, and rat. Dev Biol. 1999;214:60–71.View ArticlePubMedGoogle Scholar
- White FA, Keller-Peck CR, Knudson CM, Korsmeyer SJ, Snider WD. Widespread elimination of naturally occurring neuronal death in Bax-deficient mice. J Neurosci. 1998;18:1428–39.PubMedGoogle Scholar
- Vanderhaeghen P, Cheng H-J. Guidance molecules in axon pruning and cell death. Cold Spring Harb Perspect Biol. 2010;2:a001859.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiang C-X, Zhang K-H, Johnson RL, Jacquin MF, Chen Z-F. The transcription factor, Lmx1b, promotes a neuronal glutamate phenotype and suppresses a GABA one in the embryonic trigeminal brainstem complex. Somatosens Mot Res. 2012;29:1–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Song N-N, Xiu J-B, Huang Y, Chen J-Y, Zhang L, Gutknecht L, et al. Adult raphe-specific deletion of Lmx1b leads to central serotonin deficiency. PLoS One. 2011;6:e15998.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolfe K. Robustness - it’s not where you think it is. Nat Genet. 2000;25:3–4.View ArticlePubMedGoogle Scholar
- Postlethwait J. The zebrafish genome in context: ohnologs gone missing. J Exp Zool. 2007;308:563–77.View ArticleGoogle Scholar
- Turner J. An hereditary arthrodysplasia associated with hereditary dystrophy of the nails. JAMA. 1933;100:882–4.View ArticleGoogle Scholar
- Lucas L, Opitz JM. The nail-patella syndrome. J Pediatr. 1966;68:273–88.View ArticleGoogle Scholar
- Sweeney E, Fryer A, Mountford R, Green A, McIntosh I. Nail patella syndrome: a review of the phenotype aided by developmental biology. J Med Genet. 2003;40:153–62.View ArticlePubMedGoogle Scholar
