Abstract
From lamprey to monkeys, the organization of the descending control of locomotion is conserved across vertebrates. Reticulospinal neurons (RSNs) form a bottleneck for descending commands, receiving innervation from diencephalic and mesencephalic locomotor centers and providing locomotor drive to spinal motor circuits. Given their optical accessibility in early development, larval zebrafish offer a unique opportunity to study reticulospinal circuitry. In fish, RSNs are few, highly stereotyped, uniquely identifiable, large neurons spanning from the midbrain to the medulla. Classically labeled by tracer dye injections into the spinal cord, recent advances in genetic tools have facilitated the targeted expression of transgenes in diverse brainstem neurons of larval zebrafish. Here, we provide a comparative characterization of four existing and three newly established transgenic lines in larval zebrafish. We determine which identified neurons are consistently labeled and offer projection-specific genetic access to subpopulations of RSNs. We showcase transgenic lines that label most or all RSNs (nefma, adcyap1bccu96Et) or subsets of RSNs, including ipsilateral (vsx2, calcaccu75Et), contralateral (pcp4accu97Tg) or all (tiam2ay264Et) components of the Mauthner array, or midbrain-only RSNs (s1171tEt). In addition to RSNs, selected transgenic lines (nefma, s1171tEt, calcaccu75Et) labeled other potential neurons of interest in the brainstem. For those, we performed in situ hybridization to show expression patterns of several excitatory and inhibitory neurotransmitters at larval stages as well as glutamatergic expression patterns in juvenile fish. We provide an overview of transgene expression in the brainstem of larval zebrafish that serves to lay a foundation for future studies in the supraspinal control of locomotion.
Significance Statement
Genetic access to subpopulations of brainstem neurons greatly facilitates the dissection of supraspinal circuitry and function. Here, we present several new transgenic lines and rigorously describe existing ones, all in terms of their degree of overlap with the reticulospinal system, variability in transgenic labeling, and neurotransmitter identity. Having transgenic access to different subpopulations of reticulospinal neurons enables targeted functional calcium imaging, anatomical tracing, and optogenetic manipulations to decipher the role of individual reticulospinal neurons in movement production, ultimately clarifying existing understanding and facilitating future studies in the supraspinal control of locomotion.
Introduction
Animals are constantly faced with the need to produce flexible behavior in response to a rapidly changing environment. From foraging for food, finding a mate, to escaping from predators, movement is central to survival. As such, many regions of the central nervous system are dedicated to its production (Brownstone and Chopek, 2018). A key processing site for descending locomotor commands is the reticular formation (RF), situated in the brainstem and composed of several nuclei. The RF contains cholinergic (Josset et al., 2018; Le Ray et al., 2022), monoaminergic (McLean and Fetcho, 2004), gamma-aminobutyric acid (GABA)/glycinergic, and glutamatergic neurons (Higashijima et al., 2004), with glutamatergic reticulospinal neurons (RSNs) forming the key excitatory descending output (Perreault and Giorgi, 2019). RSNs receive descending input from the mesencephalic locomotor region (MLR) and diencephalic locomotor region and the cerebellum, as well as ascending signals from the spinal cord allowing for sensorimotor integration. Their main output is to central pattern generators composed of spinal inter- and motor neurons, responsible for movement production (for a comprehensive review on the descending control of locomotion, see Grillner and El Manira, 2020). A central question in neuroscience has long been whether different behaviors are generated by distinct or overlapping populations of RSNs (Siegel, 1979), and if so, whether activity is localized or distributed across the brainstem. The ability to label, activate, and manipulate subsets of RSNs is critical in addressing these questions.
Larval zebrafish are ideally suited for studying brainstem circuitry due to their small size, high fecundity, array of available genetic tools, optical accessibility, and well-characterized behavioral repertoire (Fetcho and Liu, 1998). For decades, researchers have labeled RSNs by spinal injections of tracer compounds, such as horseradish peroxidase or dextran-conjugated dyes (Kimmel et al., 1982; Metcalfe et al., 1986; Kimmel et al., 1990; Orger et al., 2008; Severi et al., 2014). While these have been instrumental in disentangling reticulospinal circuity, labeling can be variable depending on the skill of the experimenter, and it is challenging to consistently label descending neurons of small axon caliber. In addition, spinal injections are necessarily damaging spinal axon tracts and carry the risk of affecting subsequent behavioral studies. A solution comes from several recently established transgenic lines with expression in the brainstem (Kimura et al., 2006; Scott and Baier, 2009; Marquart et al., 2015; Eschstruth et al., 2020). Given the sometimes stochastic expression patterns in Gal4 lines (Akitake et al., 2011), it would be useful to understand precisely which RSNs are labeled by each transgenic line across a number of animals. Additionally, gaining genetic access to all as well as to selected subpopulations of RSNs, for instance, ipsi- versus contralaterally (CL) projecting RSNs, would offer an excellent opportunity when combined with optogenetic tools to decipher cell-specific functions in the control of locomotion.
Here, we characterized seven transgenic lines driving gene expression in the brainstem in larval zebrafish at 6 days-post-fertilisation: four existing transgenic lines (nefma, vsx2, s1171tEt, and tiam2ay264Et) and three newly established transgenic lines (calcaccu75Et, pcp4accu97Tg, and adcyap1bccu96Et). For each transgenic line, we performed retrograde labeling of RSNs followed by immunohistochemistry and quantified overlap across fish at the single-cell level. We showcase transgenic lines that label most or all RSNs (nefma, adcyap1bccu96Et), midbrain-only RSNs (s1171tEt), or subsets of RSNs, including ipsilateral (vsx2, calcaccu75Et), contralateral (pcp4accu97Tg), or all (tiam2ay264Et) components of the Mauthner array.
For select transgenic lines (nefma, s1171tEt, and calcaccu75Et), we performed in situ hybridization to show expression patterns of genes associated with several neurotransmitter phenotypes (vglut1, vglut2a, vglut2b, chata, gad1b, gad2, glyt1, and glyt2) at larval stages as well as glutamatergic expression in juvenile fish demonstrating consistency across development. By providing a comprehensive overview of transgene expression in the brainstem of larval zebrafish, we set a foundation for future studies in the supraspinal control of locomotion.
Materials and Methods
Fish husbandry
Adult fish were raised and bred at 28°C on a 14 h light/10 h dark cycle following standard husbandry methods as detailed in Martins et al. (2016). All fish colonies were maintained under meticulous plans involving importation of wild types every 1–3 years and line-specific breeding schemes designed to reduce inbreeding depression (Martins et al., 2016). Embryos were collected and larvae were raised at 28°C in E3 embryo medium (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, and 0.33 mM MgSO4, changed daily) at a density of 100 larvae per 200 ml until euthanasia at 6 days-post-fertilisation (dpf). A subset of fish from the calcaccu75Et line used for Extended data Figure 1-1 was raised in the presence of the pigmentation inhibitor 1-phenyl-2-thiourea (Sigma-Aldrich, P7629) at 0.2 mM from 1 dpf onwards. From 5 dpf onwards, approximately 10 ml of a live L-type rotifer poly-culture (containing 1,000–2,000 rotifers per ml) were added to each dish twice a day and larvae were allowed to feed freely. Zebrafish do not sexually differentiate until approximately 3 months of age, therefore the sex of the animals cannot be reported. All experimental procedures were approved by the Champalimaud Foundation Ethics Committee and the Portuguese Direcção Geral Veterinária, and were performed according to the European Directive 2010/63/EU.
Transgenic lines
The following transgenic lines were used in a nacre (mitfa −/ −) background (Table 1). With the exception of vsx2, which we received as an mRFP line, and s1171tEt, which was crossed with Tg[10xUAS:GCaMP6fccu1Tg] (Van Opbergen et al., 2018), the lines used in the experiments were created by crossing with Tg[10xUAS:GCaMP6fEF05ccu2Tg] (Félix et al., 2024).
Transgenic fish lines used in this study
Cloning
adcyap1b:Gal4FF
A 1,679 bp promoter region upstream of the adcyap1b start codon was cloned into pCR8/GW/TOPO (Invitrogen) using the following primers:
5′-GAACTGGAACACTTGGTGGCAGTATTG-3′
5′-GATCTGGCCAGGCTGTAAAGATACAAGAAAG-3′
The adcyap1b promoter (Entry Clone) was then recombined into an Gal4FF destination vector (Gateway LR recombination, Invitrogen), derived from the Tol2Kit (Kwan et al., 2007), so that the construct was bracketed by two Tol2 (Kawakami et al., 2000) inverted terminal repeats.
calca:Gal4FF
A 5,023 bp promoter region upstream of the calca start codon was cloned into pCR8/GW/TOPO (Invitrogen) using the following primers:
5′-GTGCCTGCTGAGGAGCATAAC-3′
5′-GGTCCCCTGTAGTAAAACATC-3′
The calca promoter (Entry Clone) was then recombined into the Tol2 Gal4FF destination vector (Gateway LR recombination, Invitrogen).
pcp4a:Gal4FF
The Tg[pcp4a:Gal4FFccu97Tg] line was generated using a bacterial artificial chromosome (BAC) recombineering following Suster et al. (2011). In short, the iTol2-amp cassette was introduced into BAC CH211-231M12. Positive clones were selected to further introduce the Gal4FF-pA-FRT-kan-FRT at the ATG site of pcp4a. To do so, homology arms, short DNA sequences of about 300 bp that flanked the ATG site and with sequence overlap to the Gal4FF-pA-FRT-kan-FRT, were amplified from the BAC:
pcp4a-HI-for CACACACCAATGCATACATCAAAGCG
pcp4a-HI-rev ACAGTAGCTTCATGGTGGCGCTGGATGAAGAGTATGAAGATGAAGGAA GAAG
pcp4a-HII-for CCAGCCTACACGCGGGTGAGCTTTCCTCCATACACATTGCACA
pcp4a-HII-rev GCGCACATACAATATCCTCCATCCCT
Similarly, the Gal4FF-pA-FRT-kan-FRT cassette was amplified using primers overlapping with the homology arms:
pcp4a-GFF-for CATCTTCATACTCTTCATCCAGCGCCACCATGAAGCTACTGTCTTCTATCGAAC
pcp4a-GFF-rev TGTGTATGGAGGAAAGCTCACCCGCGTGTAGGCTGGAGCTGCTTC
The three DNA fragments were fuzed using the Gibson Assembly Cloning Kit (New England Biolabs) and subcloned into the pCR2.1-TOPO vector (Invitrogen). The pcp4a-HI-for and pcp4a-HII-rev primers were used for further amplification of the cassette. Five hundred nanogram of it was used for recombination with the CH211-231M12/iTol2-amp BAC.
Line establishment
DNA constructs were injected together with Tol2 transposase mRNA and non-integrating UAS:GFP plasmid into 1–2 cell stage mitfa −/ − eggs. Embryos with GFP expression were raised and screened as adults for germ line transmission. Progeny of positive animals with stable expression pattern were selected as founders for the respective Gal4FF driver line.
Screening
Larvae were pre-screened at 3–4 dpf to select fish with positive expression of GCaMP or RFP. While this is a standard practice, we note that some lines appear to have less than Mendelian numbers of offspring with expression presumably due to silencing (Akitake et al., 2011). For instance, the pcp4accu97Tg line seems to be particularly prone to silencing, which can be remedied by setting multiple crosses and rigorous pre-screening.
Immunohistochemistry
The staining was performed in all seven transgenic lines (Table 1) according to a modified protocol (Randlett et al., 2015).
Dye injections
Larvae were raised in standard conditions (see above) until 5 dpf, fed with rotifers in the morning and injected in the afternoon. Larvae were injected while mounted sideways on a small petri dish containing a layer of agarose gel (SeaKem Low Electroendosmosis Agarose, #50004, Lonza), following previously described methods (O’Malley et al., 1996). The dye (10,000 MW at 50 mg/ml of Dextran, Alexa Fluor 647, Invitrogen by Life Technologies) was pressure-injected into the spinal cord near myomere 8 using a pulled capillary glass needle (GC100F-10, Harvard Apparatus) positioned with a micro-manipulator (MN-153, Narishige) and a stereo-microscope (Stereo Discovery V8, Zeiss). Pressure was applied using a pneumatic picopump (WPI, PV820). Following injection, larvae were allowed to recover overnight and checked for normal swimming behavior before euthanasia.
Immunohistochemistry
At 6 dpf, larvae were anesthetized in 15 mM tricaine (E10521, Sigma-Aldrich) for 10 min and fixed with 4% paraformaldehyde (PFA) for 2 h at room temperature while covered and with agitation. From this step onwards, larvae were kept in darkness. To stop fixation, larvae were rinsed and washed 2 × 5 min in phosphate-buffered saline with 0.25% Triton (PBT) with agitation. For epitope retrieval, larvae were rinsed twice, washed 1 × 5 min, and finally incubated in 150 mM Tris–HCl with pH 9.0 at 70°C in a water bath. Following incubation, samples were cooled on ice and washed 3 × 5 min in PBT. To permeabilise, larvae were incubated in 0.05% Trypsin–EDTA in phosphate-buffered saline (PBS) for 5 min on ice, followed by a rinse and 2 × 5 min washes in PBT. Next, larvae were incubated in blocking solution (PBS, bovine serum albumin, normal goat serum, dimethyl sulfoxide, Triton, azide, and sterilized H2O) at 4°C overnight. A primary antibody solution was prepared by diluting primary antibodies for anti-total extracellular signal-related kinase (anti-tERK) and anti-GFP or anti-mCherry (Table 2) in blocking solution at 1:500. Blocking solution was replaced by the primary antibody solution and larvae were allowed to incubate for at least three over-nights at 4°C under cover with agitation. To remove primary antibodies, larvae were rinsed three times and washed 3 × 30 min in PBT. A secondary antibody solution was prepared by diluting secondary antibodies in blocking solution at 1:500. Larvae were incubated in secondary antibody solution for at least three over-nights at 4°C with agitation. To remove secondary antibodies (Table 2), larvae were rinsed three times and washed 3 × 30 min in PBT. Samples were stored in PBT at 4°C in the dark until imaging.
Reagents used for immunohistochemistry experiments in this study
In situ hybridization chain reaction in larval zebrafish
In situ hybridization chain reaction (isHCR) reagents, including probes, hairpins, and buffers, were purchased from Molecular Instruments for detection of mRNAs. The staining was performed according to the “HCR v3.0 protocol for whole-mount zebrafish embryos and larvae” protocol provided by Molecular Instruments (Choi et al., 2018). The following transgenic lines were used in a nacre (mitfa −/ −) background: Tg[nefma:KalTa4, 10xUAS:GCaMP6fEF05], Tg[−5.0calca:Gal4FFccu75Et, 10xUAS:GCaMP6fEF05], and Tg[−0.6hsp70l:Gal4-VP16s1171tEt+, 10xUAS:GCaMP6f].
Tissue fixation
At 6 dpf, larvae were anesthetized in 15 mM tricaine (E10521, Sigma-Aldrich) and fixed with 4% PFA overnight at 4°C while covered and gently agitated. From this step onwards, larvae were continuously kept in darkness.
Tissue preparation
On the following day, larvae were washed 3 × 5 min in PBS to stop fixation, dehydrated and permeabilised with a series of 100% methanol (MeOH) washes, and stored at −20°C for up to 6 months before use. MeOH-fixed larvae were rehydrated progressively with a series of graded MeOH in PBS with Tween 0.1% (PBST) washes at room temperature with agitation for 5 min each. Next, larvae were treated with proteinase K (30 μg/ml) for 45 min at room temperature, rinsed twice with PBST, post-fixed with 4% PFA for 20 min at room temperature, and finally washed thoroughly 5 × 5 min with PBST with agitation.
Detection
Larvae were pre-hybridised with 500 μl of probe hybridization buffer for 30 min at 37°C. Probe solutions were prepared by adding 2 pmols of each probe set (Table 3) to 500 μl of probe hybridization buffer at 37°C. The pre-hybridization solution was replaced with the probe solution mix and larvae were incubated overnight at 37°C. The following day, excess probes were removed by washing larvae 4 × 15 min with 500 μl of probe wash buffer at 37°C, followed by 2 × 5 min washes with 5 × sodium chloride sodium citrate buffered with 0.1% Tween (5xSSCT) at room temperate with agitation. Note that in the isHCR 3.0 method, detection of an mRNA requires a mixture of primers, each of them including a part that is complementary to the target RNA and another part that is used for amplification of the signal (named B1–B5). Table 3 summarizes the probes indicating their target and the amplification reagent used.
Reagents used for in situ hybridization experiments in this study
Amplification
Pre-amplification was performed by incubating with 500 μl of amplification buffer for 30 min at room temperature. Next, 30 pmol of hairpin h1 and 30 pmol of hairpin h2 were prepared by snap-cooling 10 μl of 3 μM stock: hairpins were heated separately at 95°C for 90 s and allowed to cool down to room temperature in a dark drawer for 30 min. After cooling down, the hairpin solution was prepared by adding snap-cooled h1 and h2 hairpins to 500 μl of amplification buffer at room temperature. The pre-amplification solution was replaced with the hairpin solution and larvae were incubated overnight in the dark at room temperature. On the next day, excess hairpins were removed by washing the samples with 500 μl of 5×SSCT at room temperature with agitation for 2 × 5 min, 2 × 30 min, and 1 × 5 min. Samples were stored at 4°C in 5×SSCT in the dark until imaging.
isHCR in juvenile zebrafish
To perform isHCR in juvenile fish, the aforementioned isHCR 3.0 protocol for larval zebrafish from Molecular Instruments was adapted using methods from Kenney et al. (2021) as well as considering advice from Molecular Instruments on adapting the protocol for adult zebrafish (Dr. Chanpreet Singh, personal communication).
Tissue fixation
Juvenile zebrafish Tg[−0.6hsp70l:Gal4-VP16s1171tEt+] at 4 weeks post-fertilisation were anesthetized in 15 mM tricaine (E10521, Sigma-Aldrich) for 5 min and immersed in ice-cold fish facility water for 20 min. The absence of reflexes were assessed before fish were dissected caudal to cloaca using a razor blade and heads were placed in ice-cold PBS for 10 min to let blood drain. Samples were fixed in 4% PFA overnight at 4°C with gentle agitation.
Tissue preparation
The next day, samples were washed 2 × for 5 min in cold PBS, followed by careful dissection of brains into cold, sterile PBS and stored at 4°C until further processing. Following dissection, samples were split into groups of 4 and washed 3 × for 5 min in PBS at room temperature with gentle agitation. Then, samples were dehydrated using a series of MeOH/PBS mixtures (1 h each in 20%, 40%, 60%, 80% and 2 × 100% MeOH). Samples were washed 2 × for 5 min in 100% MeOH and incubated in 5% hydrogen peroxide in MeOH overnight at 4°C. The next day, samples were rehydrated using a series of MeOH/PBS mixtures (1 h each in 80%, 60%, 40%, and 20% MeOH). Samples were washed 2 × for 5 min in PBS and then 2 × for 1 h in PBS with 0.2% TritonX-100 (PTx.2). Samples were then washed overnight at 37°C in permeabilisation solution (PTx.2, 0.3 M glycine, and 20% DMSO).
Detection and amplification
The next day, samples were washed 2 × for 5 min in PTx.2. From here on, RNA detection was performed following the detection and amplification sections described above, with the following modifications: in the detection stage, probe solutions were prepared by adding 10 pmols of each probe set to 500 μl of probe hybridization buffer at 37°C. In the amplification stage, 45 pmol of hairpin h1 and 30 pmol of hairpin h2 were used by snap-cooling 15 μl of 3 μM stock.
Clearing
Samples were washed 2 × for 5 min in 50%SSC/50%PBS, dehydrated using a series of MeOH/PBS mixtures (1 h each in 20%, 40%, 60%, 80%, and 2 × 100% MeOH), and incubated in 100% MeOH overnight at 4°C. The following day, samples were incubated in 66% dichloromethane in MeOH for 3 h at room temperature. Samples were washed 2 × for 15 min in dichloromethane and then stored in dibenzyl ether until imaging.
Imaging and image processing
Imaging conditions for larval zebrafish
Samples were cut to remove most of the spinal cord for easier handling and mounted dorsal-side up in 1% low-melting point agarose (Ultrapure Low Melting Point Agarose, Cat#16520100, Invitrogen by Life Technologies) prepared in PBS or SSC. Imaging of samples was performed on an upright confocal laser-scanning microscope (Zeiss, LSM 980) with a 25 × multi-immersion objective (Zeiss, NA 0.8, Plan-Apochromat), using laser wavelengths 488, 594, and 650 nm. Images were acquired with a pixel size of 0.79 μm in x and y and sampled with a 1 μm interval in the z axis. Pinhole size was 35 μm, corresponding to 1.9 μm confocal section.
Imaging conditions for juvenile zebrafish
Samples were mounted dorsal-side up by attaching to a small petri dish lid using UV-cured glue (Bondic, UV Liquid Plastic Welder Starter Kit) and filling the lid with immersion oil (Gelest, PHENYLMETHYLSILOXANE OLIGOMER #PDM-7040) matching the refractive index of dibenzyl ether. Imaging of samples was performed on an upright confocal laser-scanning microscope (Zeiss, LSM 980) with a 25 × multi-immersion objective (Zeiss, NA 0.8, Plan-Apochromat), using laser wavelengths 594 and 650 nm. Images were acquired with a pixel size of 0.79 μm in x and y and sampled with a 2 μm interval in the z axis. Pinhole size was 35 μm, corresponding to 1.9 μm confocal section.
Image processing
Image analysis was performed using custom-written scripts in ImageJ processing package (Fiji). Briefly, where applicable, tiles were stitched and converted to a resolution that matches an in-house reference image. Using advanced normalization tools (Avants et al., 2009), each fish was registered to an in-house reference image of tERK, which labels all neurons, or a transgenic line-specific reference image. All image processing was performed on a computer with 64.0 GB RAM and a Intel(R) Core(TM) i7-5820K CPU at 3.30 GHz processor.
Data availability
Stacks of example fish for each transgenic line were uploaded in a data sharing repository (www.zenodo.org, https://doi.org/10.5281/zenodo.15102728). First channel (green) is the transgenic line, second channel (gray) is the tERK counter-stain, and third channel (magenta) is the reticulospinal backfill.
Results
In this study, we characterized four existing (nefma, Tg[nefma:KalTa4]; tiam2ay264Et, Tg[tiam2ay264Et(B)]; s1171tEt, Tg[−0.6hsp70l:Gal4-VP16s1171tEt+]; vsx2, TgBAC[vsx2: Gal4FFnns18Tg]) and three new (calcaccu75Et, Tg[−5.0calca: Gal4FFccu75Et]; pcp4accu97Tg, TgBAC[pcp4a:Gal4FFccu97Tg]; adcyap1bccu96Et, Tg[−1.7adcyap1b:Gal4FFccu96Et]) transgenic lines. Each line was crossed with Tg[10xUAS:GCaMP6fEF05] unless reported otherwise. We note that transgenic lines based on small promoter regions will have expression patterns that may reflect the locus of insertion in the genome. We have therefore denoted these lines as enhancer traps, indicating that they are not expected to be faithful reporters of the natural expression pattern of the genes from which the construct was derived. Here, we aim to characterise these transgenic lines as tools in themselves to lay a basis for dissecting the role of premotor brainstem neurons in movement production. We chose these transgenic lines as they label reticulospinal (RSNs) and other potential neurons of interest in the brainstem of larval zebrafish at 6 days-post-fertilisation (dpf) (Fig. 1A). Briefly, nefma, adcyap1bccu96Et, and tiam2ay264Et label cells across the brainstem, s1171tEt in the tegmentum, and pcp4accu97Tg, calcaccu75Et, and vsx2 in the hindbrain. Larval zebrafish of each transgenic line were injected with dextran-conjugated tracer dye to retrogradely label RSNs and compared to their transgene expression profiles following immunohistochemistry and confocal imaging.
Overview of transgenic lines and neurotransmitter-associated gene expression patterns. A, Graphical summary of methods for reticulospinal backfills paired with immunohistochemistry (left) or in situ hybridization (right), followed by confocal imaging, image registration to an average reference brain, and analysis. B, GCaMP expression in transgenic lines used in this study as indicated in each panel. C, Expression of genes associated with neurotransmitter phenotypes as indicated in each panel. For B and C, each panel shows a maximum intensity projection (MIP) from the dorsal and sagittal views, taken from average stack following image registration. Number of fish: nefma n = 7, calcaccu75Et n = 12, vsx2 n = 11, s1171tEt n = 5, pcp4accu97Tg n = 6, adcyap1bccu96Et n = 8, tiam2ay264Et n = 6, vglut1 n = 16, vglut2 n = 13, chata n = 8, gad1b/2 n = 4, glyt1 n = 16, and glyt2 n = 15. Scale bar 100 μm. For brainwide expression patterns in example fish from selected transgenic lines, see Extended data Figure 1-1.
Figure 1-1
Example fish to show brain-wide expression patterns in selected transgenic lines. The nefma line has expression in the tectum, pre-tectum, tegmentum and hindbrain, as well as labelling the anterior and posterior lateral line ganglia, the trigeminal ganglion and neuromasts. The tiam2ay264Et line has expression in the olfactory epithelium, optic chiasm, tectum, interpeduncular nucleus, tegmentum, hindbrain and cerebellum. The calcaccu75Et line mostly labels cells in the hindbrain, as well as a small number of cells in the mid- and forebrain, and outer retina. The adcyap1bccu96Et line has expression in the olfactory epithelium, olfactory bulb, tegmentum, hindbrain, anterior and posterior lateral line ganglia, as well as sparse labelling in the tectum. The pcp4accu97Tg line has sparse labelling of cells in the forebrain and tectum, expression in the habenula, hindbrain, and the anterior and posterior lateral line ganglia. Scale bar is 100μm. Download Figure 1-1, TIF file.
Here, we provide an overview of the expression patterns of each transgenic line (Extended data Fig. 1-1). The nefma line has expression in the tectum, pretectum, tegmentum, and hindbrain, as well as labeling the anterior and posterior lateral line ganglia, the trigeminal ganglion, and neuromasts. Expression in the vsx2 line is primarily in the hindbrain, whereas in the s1171tEt line, expression is restricted to the tegmentum. The tiam2ay264Et line has expression in several areas throughout the brain, including the olfactory epithelium, optic chiasm, tectum, interpeduncular nucleus, tegmentum, hindbrain, and cerebellum. The calcaccu75Et line mostly labels cells in the hindbrain, as well as a small number of cells in the mid- and forebrain, and outer retina. The adcyap1bccu96Et line has expression in the olfactory epithelium, olfactory bulb, tegmentum, hindbrain, anterior and posterior lateral line ganglia, as well as sparse labeling in the tectum. Finally, the pcp4accu97Tg line has sparse labeling of cells in the forebrain and tectum, expression in the habenula, hindbrain, and the anterior and posterior lateral line ganglia.
For transgenic lines (nefma, calcaccu75Et, and s1171tEt) that labeled cells beyond the classical RSNs, we performed isHCR to study the mRNA expression of different excitatory and inhibitory neurotransmitters: glutamate (vglut1, vglut2a, and vglut2b), acetyl choline transferase (chata), GABA (gad1b, gad2), and glycine (glyt1 and glyt2; Fig. 1B). The vesicular glutamate transporter vglut2 has two orthologs in zebrafish, vglut2a and vglut2b, which had largely overlapping brain-wide expression of mRNA, consistent with data from the zebrafish brain atlas mapzebrain (Kunst et al., 2019). Except where explicitly stated, we refer to vglut2 as a combination of vglut2a and vglut2b.
The nefma line labels all RSNs and cranial nerve nuclei
We examined the organization of RSNs in a knock-in reporter line in the nefma gene (Eschstruth et al., 2020) using dextran-conjugated retrograde labeling in larval zebrafish at 6 dpf (n = 10). The nefma line has expression in the tectum, pretectum, tegmentum, and hindbrain, as well as labeling the anterior and posterior lateral line ganglia, the trigeminal ganglion, and neuromasts (Extended data Fig. 1-1). In this line, the Mauthner cell is not present at this age because it degenerates between 3 and 5 dpf (data not shown; Dr. J. Bin and Dr. D. Lyons, The University of Edinburgh, personal communication). The Mauthner homologs (MiD2cm, MiD2cl, MiD3cm, and MiD3cl) remain intact. The back-labeled cells closely match the organization previously described (Kimmel et al., 1982), which we term the “classical” RSNs. All other classical RSNs that were back-labeled were present and expressing GCaMP (Fig. 2A,B). This includes the four named large cells (MeM1, MeLm, MeLr, and MeLc) in the nucleus of the medial longitudinal fasciculus (nMLF), RSNs in the hindbrain, including the rostral RoL1 and RoM1-3 cells, various medial cell groups (MiV1, MiV2, MiM1, MiR1, and MiR2), and the caudal MiD2, MiD3, and MiT cells (Extended data Fig. 6-1). Note that the rostro-lateral RoL2-RoL3 and the caudal CaD and CaV cells were not successfully back-labeled, most likely requiring a different injection strategy. However, it is highly likely that these cells are still present in the nefma line. We see numerous small GCaMP-positive cells lateral to the RoM3 cells, which could be the RoL2-3, and there are large GCaMP-positive cells situated caudally to the MiD3 cells most likely constituting the CaD/CaV cells (Kimmel et al., 1982; Orger et al., 2008).
The nefma line labels all RSNs as well as cranial nerve nuclei. A, MIPs from dorsal and sagittal views of an exemplary nefma fish (n = 10) at 6 dpf with RSNs labeled via retrograde dye injection. B, Close ups at several planes to illustrate overlap between nefma line and backfill, indicated by red triangles. Putative CaD and CaV cells indicated by blue triangle. C, MIPs from dorsal and sagittal views of a nefma fish (n = 24) at 6 dpf with vglut2 mRNA expression. D, Close ups at several planes to illustrate considerable overlap between nefma line and vglut2. E, MIPs from dorsal and sagittal views of a nefma fish (n = 16) at 6 dpf with chata mRNA expression. F, Close ups at several planes to illustrate cranial nerves III–VII, IX–X, and OLe are labeled in the nefma line. For A, C, and E, scale bar 100 μm, for B, D, and F, scale bar is 25 μm. For GABAergic (gad1b/2), glycinergic (glyt1 and glyt2), or glutamatergic (vglut1) expression in neurons labeled by the nefma line, see Extended data Figure 2-1.
Figure 2-1
No GABAergic (gad1b/2), glycinergic (glyt1, glyt2) or glutamatergic (vglut1) expression in neurons labelled by the nefma line. A,C,E,G) Maximum intensity projections from dorsal and sagittal view of an exemplary nefma fish at 6dpf with A) gad1b/2 (n = 8), C) glyt1 (n = 16), E) glyt2 (n = 15) or G) vglut1 (n = 14) mRNA expression. Scale bar is 100μm. Close ups at several planes to better illustrate no overlap between nefma line and B) gad1b/2, D) glyt1, F) glyt2 or H) vglut1 mRNA-expressing neurons. Scale bar is 25μm. Download Figure 2-1, TIF file.
The majority of cells labeled by the nefma line, including the classical RSNs, are glutamatergic (vglut2; Fig. 2C,D) as confirmed by isHCR (n = 24). In addition, the nefma line (n = 16) labels cholinergic cells, which could be identified as the cranial nuclei (CNIII, oculomotor; CNIV, trochlear; CNV, trigeminal; CNVI, abducens; CNVII, facial; CNIX, glossopharyngeal; CNX, vagus) and octavolateralis efferent neurons (OLe; Fig. 2E,F), with locations consistent with anatomical annotations in zebrafish atlases (Kunst et al., 2019). There was no vglut1, GABAergic (gad1b/2), or glycinergic (glyt1 and glyt2) mRNA expression in cells labeled by the nefma line (Extended data Fig. 2-1).
Transgenic lines offer access to subpopulations of RSNs
calcaccu75Et
We established a new Gal4-transgenic line using a 5kb fragment of the calca promoter. The calca and adcyap1b promoters were cloned as part of a larger screen for short promoter sequences with neural expression, based on two characteristics: (1) annotation of expression in the hindbrain, among other brain regions (www.zfin.org, Thisse et al., 2004), and (2) a compact genomic locus (around 5 kb), with a short first intron before the start codon that could be included within a short promoter sequence. The reason for this latter consideration is that many successful short promoter sequences used in zebrafish include the first intron, which has been associated with improved transgene expression (Higashijima et al., 1997). The calcaccu75Et line mostly labels cells in the hindbrain, as well as a small number of cells in the mid- and forebrain, and outer retina (Extended data Fig. 1-1). The organization of RSNs in the calcaccu75Et line was examined using dextran-conjugated retrograde labeling at 6 dpf (n = 12). Many of the classical RSNs that were visualized by the dye injection into the rostral spinal cord were also labeled by the calcaccu75Et line (Fig. 3A,B). This includes the rostral RoM2, RoM3, and RoV3 cells, various medial cell groups (MiV1, MiV2, MiM1, MiR1, and MiR2), and notably the caudal ipsilaterally projecting (IL) MiD2i and MiD3i cells, but not their CL projecting counterparts. Following the projections of MiD2i and MiD3i at higher spatial resolution confirmed their identity as no decussation at the midline occurred (data not shown). In this line, the Mauthner cell is present but not labeled. For detailed RSN labeling across multiple fish, see Extended data Figure 6-3.
The calcaccu75Et line labels only IL-RSNs and closely matches the vsx2 line. A, MIPs from dorsal and sagittal views of an exemplary calcaccu75Et fish (of n = 12) at 6 dpf with RSNs labeled via retrograde dye injection. B, Close ups at several planes to illustrate overlap between calcaccu75Et line and backfill. Note that only IL-RSNs are present in calcaccu75Et as indicated by red triangles. C, MIPs from dorsal and sagittal views of an exemplary calcaccu75Et fish (of n = 8) at 6 dpf with vglut2 mRNA expression. D, Close ups at several planes to illustrate overlap between calcaccu75Et line and vglut2, particularly of the RSNs as well as a rostro-caudal glutamatergic (vglut2) stripe in the dorsal hindbrain. E, Single planes from dorsal (105 μm) to ventral (180 μm) at 15 μm intervals in calcaccu75Et, vsx2, and composite. For both top and middle panels, images are from two single fish (of n = 12 each) that were registered to a common reference brain (tERK). For A, C, and E, scale bar 100 μm, for B and D, scale bar 25 μm. For cholinergic (chata), GABAergic (gad1b/2), or glycinergic (glyt1 and glyt2) expression in neurons labeled by the calcaccu75Et line, see Extended data Figure 3-1.
Figure 3-1
No cholinergic (chata), GABAergic (gad1b/2) or glycinergic (glyt1, glyt2) expression in neurons labelled by the calcaccu75Et line. A,C,E,G) Maximum intensity projections from dorsal and sagittal view of an exemplary calcaccu75Et fish at 6dpf with A) chata (n = 15), C) gad1b/2 (n = 7), E) glyt1 (n = 16), G) glyt2 (n = 16) mRNA expression. Scale bar is 100μm. B,D,F,H) Close ups at several planes to better illustrate no overlap between calcaccu75Et line and B) chata, D) gad1b/2, F) glyt1 or H) glyt2 or mRNA-expressing neurons. Scale bar is 25μm. Download Figure 3-1, TIF file.
There was no cholinergic (chata), GABAergic (gad1b/2), or glycinergic (glyt1 and glyt2) mRNA expression in cells labeled by the calcaccu75Et line (Extended data Fig. 3-1). In fact, most cells labeled in the calcaccu75Et line were glutamatergic (vglut2), with a particularly prominent glutamatergic stripe spanning the whole medial dorsal hindbrain (Fig. 3C,D). This strongly resembled the well-characterized vsx2 line that has a medial stripe of glutamatergic V2a neurons in the hindbrain (Kimura et al., 2006).
vsx2
We therefore compared the expression profiles of two calcaccu75Et and vsx2 fish registered to a shared reference brain. This revealed a high degree of overlap between the two transgenic lines, with the exception of additional dorso-lateral transgene expression close to the midbrain boundary in the calcaccu75Et line (Fig. 3E).
To illustrate the similarity between the calcaccu75Et and vsx2 lines at a cellular level, we also performed back-labeling of RSNs in vsx2 fish at 6 dpf (n = 12). Again, only IL projecting RSNs in the hindbrain were present in vsx2 (Fig. 4A,B). This includes the rostral RoM2, RoM3, and RoV3 cells, various medial cell groups (MiV1, MiV2, MiM1, MiR1, and MiR2), and the caudal IL projecting MiD2i and MiD3i cells. Following the projections of MiD2i and MiD3i at higher spatial resolution confirmed their identity with no decussation at the midline present (data not shown). Similarly to the calcaccu75Et line, the Mauthner cell was present but not labeled in vsx2. For detailed RSN labeling across multiple fish, see Extended data Figure 6-4.
Reticulospinal cell labeling in four different transgenic lines. A, C, E, and G, MIPs from dorsal and sagittal views of an exemplary fish from A, vsx2 (n = 12), C, pcp4accu97Tg (n = 12), E, adcyap1bccu96Et (n = 11), and G, tiam2ay264Et (n = 12) at 6 dpf with RSNs labeled via retrograde dye injection. Scale bar is 100 μm. B, D, F, and H, Close ups at several planes to illustrate overlap between B, vsx2, D, pcp4accu97Tg, F, adcyap1bccu96Et, or H, tiam2ay264Et and RS backfill. Blue triangle in F, rostral panel indicates putative RoL1. Blue triangle in F and H, caudal panel indicates putative CaD and CaV cells. Scale bar is 25 μm.
pcp4accu97Tg
We established a new transgenic line driving GAL4FF expression under the pcp4a promoter using a recombineered BAC. Pcp4a had been shown to be expressed in the telencephalon, habenula, pretectum, pre-glomerular complex, mammillary bodies, optic tectum, and a subset of RSNs (Mione et al., 2006). The pcp4accu97Tg line has sparse labeling of cells in the forebrain and tectum, expression in the habenula, hindbrain, and the anterior and posterior lateral line ganglia (Extended data Fig. 1-1). The reticulospinal back-labeling (n = 12) specifically and reliably includes the dorso-caudal CL projecting MiD2cm, MiD2cl, MiD3cm, and MiD3cl but not their ipsilateral counterparts (Fig. 4C,D). Labeling of other RSNs in this line was more variable, with different combinations of various medial cell groups (MiV1, MiV2, MiM1, MiR1, and MiR2) being labeled in 50% of fish (Extended data Fig. 6-6). To our knowledge, this is the first transgenic line specifically labeling the CL projecting MiD cells of the Mauthner array.
adcyap1bccu96Et
We established a new transgenic line in which GAL4FF expression is driven by a 1.7 kb promoter region upstream of the start codon of adenylate cyclase-activating polypeptide 1 b (adcyap1b) gene. The adcyap1b gene encodes pituitary adenylate cyclase-activating polypeptide 2 (PACAP2) and is evolutionarily conserved and expressed in the telencephalon, diencephalon, rhombencephalon, and dorsal spinal cord (Alexandre et al., 2011). PACAP2 plays a role in brain development (Wu et al., 2006) and the neuropeptide adcyap1b has been shown to enhance sensory responsiveness in larval zebrafish (Woods et al., 2014). adcyap1b was included in the previously mentioned screen based on the same characteristics described above for calca. In addition to labeling cells in the tegmentum and hindbrain, our adcyap1bccu96Et line has expression in the olfactory epithelium, olfactory bulb, anterior and posterior lateral line ganglia, as well as sparse labeling in the tectum (Extended data Fig. 1-1). Following back-filling with dextran-conjugated tracer dye (n = 11), we report overlap with GCaMP expression in the four identified cells of the nMLF (MeM1, MeLm, MeLr, and MeLc), the Mauthner cell, and most other RSNs including the vestibular cells (Fig. 4E,F). We observed labeling of most of the caudal CL projecting MiD cells (MiD2cm, MiD2cl, MiD2i, MiD3cm, and MiD3cl) but not the MiD3i and MiT cells. For detailed RSN labeling across multiple fish, see Extended data Figure 6-5. As mentioned previously, the CaD, CaV, RoL1, and RoL2 cells were not retrogradely labeled in our preparation. However, it is likely that the adcyap1bccu96Et line labels those cells well, as there are large cell bodies visible in their respective putative locations (Kimmel et al., 1982; Orger et al., 2008).
tiam2ay264Et
The previously established tiam2ay264Et line (Marquart et al., 2015) was of interest due to labeling the Mauthner cell and homologs. It has expression throughout the brain, including the olfactory epithelium, optic chiasm, tectum, interpeduncular nucleus, tegmentum, hindbrain, and cerebellum (Extended data Fig. 1-1). Backfills (n = 12) revealed that this line labeled MeM1 more reliably than the other identified nMLF neurons (MeLm, MeLr, and MeLc), as well as consistently labeling the rostral RoM2l, RoM3m, and RoM3l, the Mauthner cell, and both the CL- and IL-projecting MiD2 and MiD3 cells in the hindbrain (Fig. 4G,H). Labeling of other RSNs was more variable across fish (Extended data Fig. 6-7). Similarly to above, it is likely that the CaD and CaV cells are also labeled in the tiam2ay264Et line (Kimmel et al., 1982; Orger et al., 2008).
The s1171tEt line labels vglut2-expressing neurons in the nMLF
The organization of midbrain RSNs in the s1171tEt line was examined using dextran-conjugated retrograde labeling in larval zebrafish at 6 dpf (n = 12). The s1171tEt has a triangular expression profile centered around the nMLF in the tegmentum, rostral to the cholinergic oculomotor nucleus (Thiele et al., 2014). It includes the four canonical large nMLF cells (MeM1, MeLm, MeLr, and MeLc), as well as small mesencephalic cells (MeS; Fig. 5A,B). For detailed RSN labeling across multiple fish, see Extended data Figure 6-2.
The s1171tEt line labels RSNs in the nMLF and is glutamatergic (vglut2) across development. A, Dorsal and sagittal MIPs from an example s1171tEt fish (of n = 12) at 6 dpf with RSNs labeled via retrograde dye injection. Scale bar 100 μm. B, Close ups at several planes to illustrate overlap between s1171tEt line and backfill, note the four main RSNs as indicated by red triangles. Scale bar 25 μm. C, Close ups at several planes of exemplary fish to show considerable overlap between s1171tEt and vglut2 (n = 8), particularly in RSNs. No overlap between s1171tEt and vglut1 (n = 7), chata (n = 8), gad1b/2 (n = 8), or glyt2 (n = 7). Scale bar 25 μm. D, vglut1 and vglut2 expression in exemplary 4-week-old juvenile fish (n = 4). Left) Two planes in dorsal view, scale bar 100 μm. Right) Matching transverse view, scale bar 50 μm. Blue triangles indicate blood vessels, red triangles point toward region of nMLF. For more examples of glutamatergic expression patterns in 4-week-old juvenile fish, see Extended data Figure 5-1.
Figure 5-1
Glutamatergic expression patterns in three 4 week-old juvenile fish (of n = 4). A,C,E) Left panels show vglut1, vglut2, composite in two planes from the dorsal view, scale bar is 100μm. Right panels show vglut1, vglut2, composite from the transverse view, scale bar is 50μm. B,D,F) Composite images of several planes from dorsal to ventral. Note the presence of torus longitudinalis (*) and cerebellum (+) in vglut1 (green), with brain wide expression in vglut2 (magenta). Auto-fluorescence of blood vessels is seen in the vglut1 channel, as indicated by blue triangles. Scale bar 200μm. Download Figure 5-1, TIF file.
The majority of cells in the s1171tEt line were glutamatergic, specific for the vesicular glutamate transporter 2 (vglut2; Fig. 5B). We did not detect vglut1-expressing cells near the nMLF, as well as no cholinergic (chata), GABAergic (gad1b/2), or glycinergic (glyt1 and glyt2) mRNA expression in cells labeled by the s1171tEt line (Fig. 5C). A previous study distinguished lateral and medial regions of the nMLF in zebrafish at 4-weeks post-fertilisation based on different levels of expression observed in vglut1 and vglut2a transgenic lines (Berg et al., 2023). This difference with the isHCR results could reflect a change in the expression pattern over development, or alternatively could be due to expression in these transgenic lines that does not reflect the natural gene expression patterns.
To understand this, we asked whether glutamatergic expression profiles change across development. We performed a modified isHCR protocol on brains dissected from juvenile zebrafish (n = 4; 4-weeks post-fertilisation). Using landmarks identified with the adult zebrafish brain atlas (AZBA; Kenney et al., 2021)—such as the torus longitudinalis (TL), valvula cerebelli, and the diencephalic ventricle—allowed us to identify several brain regions. Glutamatergic expression appears stable across development, with vglut1 expression being restricted to the TL and the cerebellum, in contrast to brain-wide vglut2 expression, including the putative nMLF (Extended data Fig. 5-1). To identify the nMLF in the juvenile brain, we again utilized the AZBA. Matching the triangular shape of the nMLF on AZBA in the transverse view, situated between the diencephalic ventricle and a cluster of vglut2 cells most likely belonging to the red nucleus, is a cluster of large vglut2 cells characteristic of the nMLF (Fig. 5D). We did not observe vglut1 expression near this cluster of cells, only auto-fluorescent blood vessels.
Summary of RSN labeling across transgenic lines
The transgenic lines characterized in this study can be utilized in a complementary fashion to study reticulospinal circuits at larval stages (Fig. 6). The nefma line (n = 10) reliably labels all RSNs. However, this line cannot be used to study the Mauthner cell since it degenerates early in development. The adcyap1bccu96Et line (n = 11) labels most RSNs including the Mauthner cell. However, the MiD3i and MiT cells are not labeled, and we emphasize the need for rigorous pre-screening due to silencing. The most suitable line to study Mauthner and homolog activity is the tiam2ay264Et line (n = 12), with labeling of other RSNs being more variable. The calcaccu75Et/vsx2 and pcp4accu97Tg lines offer an interesting opportunity by giving genetic access to complementary projection patterns in the Mauthner array, with the calcaccu75Et line (n = 12) only labeling the IL projecting MiD cells (as well as other rostral and medial RSNs) and the pcp4accu97Tg line (n = 12) specifically labeling CL projecting MiD cells, while labeling of various medial RSNS is much more rare. In addition, we have demonstrated that the calcaccu75Et line closely matches the existing vsx2 line (n = 12), both in broad expression pattern and at the single-RSN-level—though it labels additional cells in the rostro-lateral hindbrain. For exact numbers across fish for all transgenic lines, see Extended data Figures 6-1–6-7. Finally, for wider nMLF-related studies, the s1171tEt line (n = 12) offers broad labeling in the tegmentum, including the four canonical nMLF-RSNs (MeM1, MeLr, MeLc, and MeLm). We further demonstrate that cells labeled in these transgenic lines are glutamatergic (vglut2), with the nefma line additionally labeling the cholinergic cranial neurons.
Graphical summary of the number of RSNs labeled in each transgenic line across multiple fish. nefma n = 10, calcaccu75Et n = 12, vsx2 n = 12, s1171tEt n = 12, pcp4accu97Tg n = 12, adcyap1bccu96Et n = 11, and tiam2ay264Et n = 12. For data on which reticulospinal cells are labeled in each fish per transgenic line, see Extended data Figures 6-1–6-7.
Figure 6-1
Reticulospinal cells labelled in each nefma fish at 6dpf. Download Figure 6-1, XLSX file.
Figure 6-2
Reticulospinal cells labelled in each s1171tEt fish at 6dpf. Download Figure 6-2, XLSX file.
Figure 6-3
Reticulospinal cells labelled in each calcaccu75Et fish at 6dpf. Download Figure 6-3, XLSX file.
Figure 6-4
Reticulospinal cells labelled in each vsx2 fish at 6dpf. Download Figure 6-4, XLSX file.
Figure 6-5
Reticulospinal cells labelled in each adcyap1bccu96Et fish at 6dpf. Download Figure 6-5, XLSX file.
Figure 6-6
Reticulospinal cells labelled in each pcp4accu97Tg fish at 6dpf. Download Figure 6-6, XLSX file.
Figure 6-7
Reticulospinal cells labelled in each tiam2ay264Et fish at 6dpf. Download Figure 6-7, XLSX file.
Discussion
The aim of this study was to provide a detailed account of seven transgenic lines with gene expression in the brainstem of larval zebrafish. We characterized four existing transgenic lines (nefma, vsx2, s1171tEt, and tiam2ay264Et) and presented three newly established transgenic lines (calcaccu75Et, pcp4accu97Tg, and adcyap1bccu96Et). For each transgenic line, we performed retrograde labeling of RSNs followed by immunostaining against GCaMP/RFP and the whole-brain neural marker total ERK for subsequent image registration purposes.
As previously reported (Thiele et al., 2014), the s1171tEt line reliably labels the four large cells (MeLr, MeLc, MeLm, and MeM1) in the nMLF as well as many other neurons in the surrounding area. We observed occasional overlap with other backfilled RSNs, most likely constituting the small mesencephalic cells (MeS), first characterized by Thiele et al. (2014). The nMLF in fish is thought to be homologous to the interstitial nucleus of Cajal (Yamamoto et al., 2017; Carbo-Tano et al., 2023) or a part of such a structure (Wullimann et al., 1996), due to its involvement in postural control and projections to the spinal cord during early development (Thiele et al., 2014; Liu and Bagnall, 2023; Sugioka et al., 2023).
Two transgenic lines labeled almost all RSNs, with notable exceptions: the nefma line reliably labels all RSNs except for the Mauthner cell, which degenerates earlier in development. We can speculate as to why the Mauthner cell degenerates in this transgenic line. The degeneration appears to be independent of the presence of a UAS reporter, consistent with reports that the KaLTA4 construct (Distel et al., 2009) can be cytotoxic at high concentrations (Zhang et al., 2019). An alternative explanation is that high axon caliber neurons like the Mauthner cell are more sensitive to perturbations of nefma gene expression that arise from the insertion of the transgene in that locus (Baraban et al., 2013). Another transgenic line has been made by insertion into the same locus (Liu et al., 2020), in which the Mauthner cell does not degenerate (Hamling et al., 2023). Multiple differences exist between the two lines, which could account for the discrepancy. First, in the line used here, KalTA4 is integrated in frame at the end of the nefma gene with a self-cleaving peptide sequence between the two, while in the other line the construct is integrated upstream of the nefma gene. Second, a different GAL4 construct was used, with the other line using GAL4FF rather than KalTA4. It is likely that the other line has a similar expression pattern, with most RSNs labeled, but this remains to be systematically documented. We were also not able to confirm the labeling of the rostral RoL2-3 and caudal CaD and CaV neurons, as these were not successfully labeled by back-filling in our preparation and may require a different injection strategy. However, there are GCaMP-expressing cells in the nefma line where we would expect the RoL2-3 and CaD and CaV cells to be based on previous studies (Kimmel et al., 1982; Orger et al., 2008). We thus conclude that the nefma line labels all RSNs except the Mauthner cell. The newly generated adcyap1bccu96Et line reliably labels most RSNs including the Mauthner cells and most homologs and the putative CaD and CaV cells, with a specific exception: we did not observe labeling of the MiT and IL projecting MiD3i cells in the adcyap1bccu96Et line.
To study the Mauthner array, three complementary transgenic lines can be used. The tiam2ay264Et line reliably labels the Mauthner cells and homologs (including the MiD2i and MiD3i); however, labeling of other RSNs is stochastic, as previously reported by Tabor et al. (2014). We presented two new lines that offer complementary genetic access to different projection patterns within the Mauthner array. The calcaccu75Et line largely overlaps with the well-known line vsx2, with the addition of gene expression in the rostral lateral hindbrain. The calcaccu75Et line reliably labels the RoM2, RoM3, and various medial RSNs as well as the IL projecting Mauthner homologs MiD2i and MiD3i. Conversely, the pcp4accu97Tg line reliably labels the CL projecting medial and lateral MiD2cm, MiD2cl, MiD3cm, and MiD3cl as well as more rarely various medial cells. To our knowledge, this is the first description of a transgenic line in zebrafish providing specific genetic access to the CL projecting Mauthner homologs.
There is ample evidence for the usefulness of transgenic lines labeling different subpopulations in the hindbrain to dissect premotor circuitry and function. For instance, using a Gal4-transgenic line under the nefma promoter allowed experimenters to reliably record convergent afferent signals in vestibulospinal neurons in vivo (Liu et al., 2020). It was later shown in the same line that vestibulospinal neuron activity increases as a function of ipsilateral tilt amplitude using tilt-in-place microscopy (Hamling et al., 2023).
Another transgenic line identified in a Gal4 enhancer trap screen, s1171tEt, with expression in the thalamus, cerebellum, and trunk musculature (Scott and Baier, 2009) was instrumental in showing that the four large (MeLr, MeLc, MeLm, and MeM1) as well as smaller (MeS) nMLF neurons are involved in steering movements (Thiele et al., 2014), tail bending (Dal Maschio et al., 2017), and relay commands from the pretectal AF7 to the hindbrain, potentially involved in hunting behavior (Semmelhack et al., 2014). In addition to their terminal arbors in the spinal cord, it was also demonstrated that nMLF cells labeled in the s1171tEt line have extensive axon collaterals in the hindbrain, suggesting they may be involved in coordination of different premotor regions (Thiele et al., 2014).
Information on RSNs involved in steering behaviors also comes from another study using a transgenic line that labels V2a neurons, showing they receive innervation from the MLR and are also involved in forward swimming (Carbo-Tano et al., 2023). Initially generated as a BAC line named chx10 after the mammalian homolog, it has been named alx and vsx2 in fish. It labels medial IL projecting cells in the hindbrain and spinal cord, shown to be glutamatergic descending interneurons providing excitatory drive to motoneurons in the spinal cord (Kimura et al., 2006). It was subsequently confirmed that they overlap with the most medial V2a glutamatergic stripe in the hindbrain (Kinkhabwala et al., 2011). Optogenetic activation of chx10-expressing neurons in the hindbrain using channelrhodopsin evoked swimming, while forced inactivation using Archearhodopsin3 or Halorhodopsin reliably stopped ongoing swimming (Kimura et al., 2013).
Finally, the tiam2ay264Et line labels the Mauthner array reliably and other RSNs more stochastically, as well as the anterior lateral line ganglia amongst others. Studies using the tiam2ay264Et line showed that direct activation of the Mauthner cell by electric field pulses could drive ultrarapid escape responses (Tabor et al., 2014), and that it receives pre-pulse inhibition from gsx2-glutamatergic neurons (Tabor et al., 2018).
Together, these studies illustrate the usefulness of transgenic lines labeling different sets of RSNs in order to decipher their role in movement control. Another option to record neural activity in this population is to retrogradely label the neurons by spinal injections of a dextran-conjugated fluorescent calcium indicator dye, such as Cal520 (Tada et al., 2014). However, spinal injections can be variable, require experimental skill, and are time consuming. Importantly, their invasive nature, which necessarily damages the axons of the cells of interest, has the potential to impact the fish’s behavior, altering the circuits that we wish to study (Gahtan and O’Malley, 2001). To complement functional imaging studies, the transgenic lines described here can be used to express optogenetic tools to perturb neuronal function, enabling optical manipulation of cell activity. Since all of the lines have expression in multiple targets, some additional spatial targeting of the incident light will be needed to achieve specificity. This has been achieved in practice in zebrafish through the use of optical fibers (Arrenberg et al., 2009), patterned projections (Wyart et al., 2009; Zhu et al., 2012), and 2-photon holography (Dal Maschio et al., 2017). Despite this caveat, restriction of expression to a small subset of possible neurons will minimize off-target effects and simplify experimental design. Other tools can benefit from similar intersections of optical and genetic specificity. For example, the photosensitizer KillerRed can be used to target cells for light-mediated ablation (Teh et al., 2010), while optical highlighter proteins can be used to reveal single neuron morphologies and projections with minimal background (Datta and Patterson, 2012).
In addition to reticulospinal backfills, we confirmed that cells labeled by the nefma, calcaccu75Et, and s1171tEt transgenic lines are glutamatergic (vglut2). The nefma line also labels the cholinergic cranial motor neurons III–VII and IX–X. A previous study had reported medial vglut2-positive and lateral vglut1-positive cells in the nMLF in zebrafish at 4–6 weeks post-fertilisation based on expression patterns of transgenic lines (Berg et al., 2023). We did not observe vglut1-positive cells in the nMLF in our study on larval zebrafish, instead, vglut1-expression was restricted to the torus longitudinalis (TL) and cerebellum. This finding is corroborated by isHCR data available in the larval zebrafish atlas mapzebrain (Kunst et al., 2019). We thus wanted to understand whether glutamatergic expression profiles change over development. We performed isHCR against vglut1 and vglut2 in brains dissected from 4-week-old zebrafish, an age consistent with the study by Berg et al. (2023). Similarly to the larvae, we only observed vglut1-expression in the TL and cerebellum, identified using the AZBA (Kenney et al., 2021). Vglut2-expression was more widespread, including expression in a region likely constituting the nMLF—identified using landmarks, such as the diencephalic ventricle, from AZBA. Again, we wondered whether glutamatergic expression profiles in the nMLF would change at later stages of development. However, this does not seem to be the case: a study in adult zebrafish (older than 90 days) using in situ hybridization reports vglut1-expression in the TL and cerebellum only (Bae et al., 2009). We therefore think that the distinction between nMLF subregions in Berg et al. (2023) reflects variation in transgenic expression patterns, rather than the natural expression of these genes. However, this hypothesis needs to be tested by directly comparing in situ hybridization data with expression in those transgenic lines.
In conclusion, we have provided a detailed account of the degree of RSN labeling in seven transgenic lines with glutamatergic (vglut2) expression in the brainstem in larval zebrafish, offering projection-specific genetic access to subpopulations of RSNs. This resource provides a useful basis for future research, including the possibility to combine with genetically encoded activity reporters and optogenetic tools, in the hope to uncover fundamental principles in the descending control of locomotion.
Footnotes
The authors declare no competing financial interests.
We thank Dr. S. Lackner (Champalimaud Foundation) for cloning of the calca promoter fragment for Tg[−5.0calca:Gal4FFccu75Et]. We thank I. Oliveira and J. Monteiro (Champalimaud Foundation) for early screening of the Tg[nefma:KalTa4] line. We thank the Champalimaud Foundation Advanced Bio-imaging and Bio-optics Experimental Platform for technical support with confocal imaging and the Champalimaud Foundation Fish Platform for meticulous fish care. We thank Dr. Chanpreet Singh (Molecular Instruments) for technical advice on adapting the in situ hybridization protocol to juvenile fish. We thank J. Bin and D. Lyons (University of Edinburgh, UK) for sharing unpublished data on the Mauthner cell degeneration in Tg[nefma:KalTa4]. We would like to thank the Giudicelli lab (Sorbonne University, France) for generating and sharing embryos from Tg[nefma:KalTa4], the Schuster lab (University of Bayreuth, Germany) for sending embryos from Tg[tiam2ay264Et(B)], and the Wyart lab (Paris Brain Institute, France) for sending embryos from TgBAC[vsx2:Gal4FFnns18Tg, UAS:mRFP]. We were supported by Champalimaud Foundation. Additionally, E.C. was funded by a grant from the Fundação para a Ciência e a Tecnologia, FCT (Bolsa SFRH/BD/147089/2019) and M.O. secured funding through an ERC Consolidator Grant (Neurofish-DLV-773012), the Volkswagen Stiftung Life? Initiative (A126151), and Champalimaud Foundation. Fish care was supported by the research infrastructure CONGENTO, co-financed by Lisboa Regional Operational Programme (Lisboa2020), under the PORTUGAL 2020 Partnership Agreement through the European Regional Development Fund (ERDF) and FCT under the project LISBOA-01- 0145-FEDER-022170.
This paper contains supplemental material available at: https://doi.org/10.1523/ENEURO.0581-24.2025
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.