Specialized neurons in the right habenula mediate response to aversive olfactory cues
Abstract
Hemispheric specializations are well studied at the functional level but less is known about the underlying neural mechanisms. We identified a small cluster of cholinergic neurons in the dorsal habenula (dHb) of zebrafish, defined by their expression of the lecithin retinol acyltransferase domain containing 2 a (lratd2a) gene and their efferent connections with a subregion of the ventral interpeduncular nucleus (vIPN). The lratd2a-expressing neurons in the right dHb are innervated by a subset of mitral cells from both the left and right olfactory bulb and are activated upon exposure to the odorant cadaverine that is repellent to adult zebrafish. Using an intersectional strategy to drive expression of the botulinum neurotoxin specifically in these neurons, we find that adults no longer show aversion to cadaverine. Mutants with left-isomerized dHb that lack these neurons are also less repelled by cadaverine and their behavioral response to alarm substance, a potent aversive cue, is diminished. However, mutants in which both dHb have right identity appear more reactive to alarm substance. The results implicate an asymmetric dHb-vIPN neural circuit in the processing of repulsive olfactory cues and in modulating the resultant behavioral response.
Editor's evaluation
This work presents a conceptual advance on our understanding of the habenula in vertebrate species, by revealing interesting functions to specific cell types within this region of the brain.
https://doi.org/10.7554/eLife.72345.sa0Introduction
Fish use the sense of smell to search for food, detect danger, navigate, and communicate social information by detecting chemical cues in their aquatic environment (Yoshihara, 2014). As with birds and mammals, perception of olfactory cues is lateralized and influences behavior (Siniscalchi, 2017). In zebrafish, nine glomerular clusters in the olfactory bulb (OB) receive olfactory information from sensory neurons in the olfactory epithelium and, in turn, transmit signals to four forebrain regions: the posterior zone of the dorsal telencephalon (Dp), the ventral nucleus of the ventral telencephalon (Vv), the posterior tuberculum (PT), and the dorsal habenular region (dHb) (Miyasaka et al., 2014; Yoshihara, 2014). In contrast to all other target regions, which are located on both sides of the forebrain, only the right nucleus of the dHb is innervated by mitral cells that emanate from medio-dorsal (mdG) and ventro-medial (vmG) glomerular clusters in both OBs (Miyasaka et al., 2014; Yoshihara, 2014). Moreover, calcium imaging experiments suggest that the right dHb shows a preferential response to odorants compared to the left dHb (Chen et al., 2019; Dreosti et al., 2014; Jetti et al., 2014; Krishnan et al., 2014). The identity of the post-synaptic neurons within the right dHb that receive olfactory input and the purpose of this asymmetric connection are unknown.
The habenulae are highly conserved structures in the vertebrate brain and, in teleosts such as zebrafish, consist of dorsal and ventral (vHb) nuclei, which are equivalent to the medial and lateral habenulae of mammals, respectively (Amo et al., 2010). The neurons of the dHb are largely glutamatergic and contain specialized subpopulations that also produce acetylcholine, substance P, or somatostatin (deCarvalho et al., 2014; Hsu et al., 2016; Lee et al., 2019). In zebrafish, the number of neurons within each subtype differs between the left and right dHb (deCarvalho et al., 2014). The dHb have been implicated in diverse states such as reward, fear, anxiety, sleep, and addiction (Duboué et al., 2017; Hikosaka, 2010; Lee et al., 2019; Okamoto et al., 2012). Accordingly, the right dHb was shown to respond to bile acid and involved in food-seeking behaviors (Chen et al., 2019; Krishnan et al., 2014), whereas the left dHb was found to be activated by light and attenuate fear responses (Dreosti et al., 2014; Duboué et al., 2017; Zhang et al., 2017). However, the properties of the dHb neurons implicated in these behaviors, such as their neurotransmitter identity and precise connectivity with their unpaired target, the midbrain interpeduncular nucleus (IPN), have yet to be determined.
Here, we describe a group of cholinergic neurons defined by their expression of the lecithin retinol acyltransferase domain containing 2 a (lratd2a) gene [formerly known as family with sequence similarity 84 member B (fam84b)], that are predominantly located in the right dHb where they are selectively innervated by the olfactory mitral cells that originate from both sides of the brain (Miyasaka et al., 2009), and form efferent connections with a restricted subregion of the ventral IPN (vIPN). Activity of the lratd2a-expressing neurons is increased following exposure to the aversive odorant cadaverine and, when these neurons are inactivated, adult zebrafish show a diminished repulsive response. Our findings provide further evidence for functional specialization of the left and right habenular nuclei and reveal the neuronal pathway that mediates a lateralized olfactory response.
Results
lratd2a-expressing neurons in the right dHb receive bilateral olfactory input
A subset of medio-dorsal mitral cells that are labeled by Tg(lhx2a:gap-YFP) were previously shown to project their axons bilaterally through the telencephalon and terminate in the right dHb (Miyasaka et al., 2009), in the vicinity of a small population of lratd2a-expressing neurons [(deCarvalho et al., 2013) and Figure 1A]. To characterize this subset of dHb neurons, we used CRISPR/Cas9-mediated targeted integration (Kimura et al., 2014) to introduce the QF2 transcription factor (Ghosh and Halpern, 2016; Subedi et al., 2014) under the control of lratd2a transcriptional regulation (Figure 1B). QF2 does not disrupt transcription at the lratd2a locus (Figure 1—figure supplement 1) and drives expression of QUAS-regulated genes encoding fluorescent protein reporters in a similar pattern to endogenous lratd2a gene expression in the nervous system (deCarvalho et al., 2013). In 5 days post fertilization (dpf) larval zebrafish, this includes a subset of neurons in the OB, the bilateral vHb, and asymmetrically distributed neurons in the dHb. There are 2.5 times more cells in the right dHb than the left (Figure 1C–E’’ and Figure 1—figure supplement 2) and those on the left consistently show weaker lratd2a expression and less GFP labeling than those on the right (Figure 1A and Figure 1—figure supplements 1–2). Double labeling confirmed that axons of the bilateral lhx2a positive olfactory mitral cells terminate precisely at the lratd2a-expressing neurons in the right dHb (Figure 1C–D’’, Figure 1—figure supplement 3B), co-localized with the pre-synaptic protein synaptophysin (Figure 1—figure supplement 3A). The lratd2a-expressing dHb neurons, in turn, project to a restricted region of the ventral IPN (Figure 1F–G’’).
Aversive olfactory cues activate lratd2a neurons in the right dHb
Several studies using transgenic expression of the genetically encoded calcium indicator GCaMP have demonstrated that dHb neurons are activated by olfactory cues in larval zebrafish (Jetti et al., 2014; Krishnan et al., 2014). We had previously examined the habenular response of adult zebrafish to cadaverine, a known aversive olfactory cue that is released from decaying fish (Hussain et al., 2013), and to chondroitin sulfate, a component of alarm substance (also known as Schreckstoff), which is released from the skin of injured fish (Mathuru et al., 2012). However, we did not have the necessary transgenic tools to specifically record the activity of the lratd2a cell population (deCarvalho et al., 2014).
To determine whether olfactants activate the lratd2a-positive neuronal cluster in the dHb, we measured calcium signaling in Tg(lratd2a:QF2)c644, Tg(QUAS:GCaMP6f)c587 larvae at 7 dpf. We monitored the responses of individual cells in the right dHb, as GCaMP6f labeling was weakly or not detected in neurons on the left (data not shown). The observed changes in fluorescence intensity after the addition of two applications of an aversive cue were compared to those observed after two additions of vehicle alone (i.e. deionized water), all delivered 1 min apart (Figure 2A and B). Cadaverine elicited a significant increase in GCaMP6f labeling in the same neurons following the first and second applications (16-fold increase averaged over both doses; Figure 2C). On average, a greater than fourfold increase in calcium signaling was measured in individual lratd2a-expressing neurons following addition of chondroitin sulfate compared to their response to vehicle alone (Figure 2D).
To examine whether the response to aversive odorants persists in the olfactory-dHb pathway of adult zebrafish, we used expression of the fos gene as an indicator of neuronal activation (deCarvalho et al., 2013; Hong et al., 2013). Consistent with previous findings (Dieris et al., 2017), cadaverine broadly activated OB mitral neurons in the dorsal glomerulus (dG), dorso-lateral glomerulus (dlG), medio-anterior glomerulus (maG), medio-dorsal glomerulus (mdG), lateral glomerulus (lG). In addition, we observed a threefold increase in the number of fos-expressing neurons in the right dHb following exposure of adult zebrafish to cadaverine relative to delivery of water alone (15.73 ± 1.25 vs 5.61 ± 1.07 cells, Figure 2E and F). The position of the fos-expressing cells in the right dHb corresponded to that of the lratd2a-expressing neurons (11.00% ± 1.07 of lratd2a-expressing region; Figure 2F’), supporting that the lratd2a subpopulation responds to cadaverine in both larvae and adults.
Exposure to alarm substance prepared from adult zebrafish increased the number of fos-expressing cells in the lateral glomerulus (lG) and dlG of the OB as would be expected (Mathuru et al., 2012; Yoshihara, 2014), but also in the dorso-lateral region of the dHb (Figure 2G and G’), where transcripts co-localized to lratd2a-expressing neurons. In contrast to cadaverine, alarm substance activated neurons equally in both the left and right dHb (20.72 ± 2.70 cells on the left and 20.31 ± 2.53 on the right; Figure 2H).
Synaptic inhibition of right dHb lratd2a neurons reduces aversive response to cadaverine
To confirm that lratd2a expressing dHb neurons play a role in processing of aversive olfactory cues, we inhibited synaptic transmission in these cells and tested adults for their reaction to cadaverine. The odorant was introduced at one end of a test tank and the time individuals spent within or outside of this region of the tank was measured. We used a preference index that is based on the position of an individual fish at a given time relative to the application site of the odorant (refer to Materials and methods).
The Tg(lratd2a:QF2) driver line is expected to inhibit lratd2a-expressing neurons in both the vHb as well as in the right dHb. We therefore devised an intersectional strategy to block the activity of neurons selectively in the right dHb, which combines Cre/lox-mediated recombination (Förster et al., 2017; Satou et al., 2013; Tabor et al., 2019) and the QF2/QUAS system (Ghosh and Halpern, 2016; Subedi et al., 2014). We produced transgenic fish expressing Cre recombinase under the control of the endogenous solute carrier family 5 member 7 a (slc5a7a) gene using CRISPR/Cas9 targeted integration (Kimura et al., 2014; Figure 3A–B). slc5a7a encodes a choline transporter involved in acetylcholine biosynthesis and, in zebrafish larvae, is strongly expressed in the right dHb and not in the vHb (Hong et al., 2013). Accordingly, in larvae bearing the three transgenes Tg(lratd2a:QF2)c601, Tg(slc5a7a:Cre)c662 and Tg(QUAS:loxP-mCherry-loxP-GFP-CAAX)c679, Cre-mediated recombination resulted in a switch in reporter labeling from red to green in right dHb neurons (Figure 3C). We followed a similar approach to inhibit synaptic transmission from lratd2a right dHb neurons using Botulinum neurotoxin (Lal et al., 2018; Sternberg et al., 2016; Zhang et al., 2017). A BoTxBLC-GFP fusion protein was placed downstream of a floxed mCherry reporter to generate Tg(QUAS:loxP-mCherry-loxP-BoTxBLC-GFP)c674. To validate the effectiveness of this transgenic line, a neuron specific promoter from a Xenopus neural-specific beta tubulin (Xla.Tubb2) gene was used to drive QF2 expression. Larvae bearing Tg(Xla.Tubb2:QF2;he1.1:mCherry)c663; Tg(slc5a7a:Cre)c662 and Tg(QUAS:loxP-mCherry-loxP-BoTxBLC-GFP)c674 showed a significantly reduced response to a touch stimulus, indicating that the neurotoxin was produced in the presence of Cre recombinase (Figure 3—figure supplement 2, Video 1).
BoTxBLC-GFP was selectively expressed in lratd2a/slc5a7a neurons of the right dHb (Figure 3D) in individuals bearing the three transgenes Tg(lratd2a:QF2)c601, Tg(slc5a7a:Cre)c662, and Tg(QUAS:loxP-mCherry-loxP-BoTxBLC-GFP)c674. Axons labeled by BoTxBLC-GFP terminated at the vIPN (Figure 3E), in the same location as those observed in Tg(lratd2a:QF2), Tg(QUAS:mApple-CAAX) fish (Figure 1G) suggesting that botulinum neurotoxin inhibits synaptic transmission within this restricted region of the vIPN. We confirmed that BoTxBLC-GFP labeling persisted in the right dHb neurons throughout development, although some variability was observed in the number of BoTxBLC-GFP-positive cells between individuals (Figure 3—figure supplement 3).
To determine whether the lratd2a neurons in the right dHb contributed to the aversive response to cadaverine, we monitored the behavior within and between groups of adults that had or did not have the BoTxBLC-GFP transgene. Fish lacking the transgene showed a significantly reduced preference for the side of the tank where cadaverine had been applied (Figure 3G).
When the response of individual fish within each group was compared over time (Figure 3—figure supplement 1), adults with or without BoTxBLC-GFP initially avoided the side of the test tank where cadaverine had been introduced. However, aversion was sustained for 4 min in control fish, but not in those expressing BoTxBLC-GFP in lratd2a neurons. These findings, from statistical tests on individuals both within and between groups, suggest that lratd2a neurons in the right dHb are required for a prolonged aversive response to cadaverine.
Disruption of synaptic transmission in lratd2a-expressing Hb neurons alone did not alter the response of zebrafish to alarm substance, which typically triggers erratic, rapid swimming and bottom dwelling, followed by freezing behavior (Diaz-Verdugo et al., 2019; Jesuthasan and Mathuru, 2008). Similar to controls, both juveniles and adults expressing BoTxBLC-GFP under the control of Tg(lratd2a:QF2)c601 showed rapid swimming/darting behavior within 22–25 s after delivery of alarm substance, first doubling their speed of swimming (Figure 3H–K and Figure 3—figure supplement 4), and then freezing for the duration of the 5 min recording period. Blocking the activity of lratd2a neurons in the right dHb is therefore insufficient to diminish the robust behavioral changes elicited by alarm substance (Figure 3H–K and Figure 3—figure supplement 4).
Zebrafish mutants with habenular defects show altered responses to aversive cues
We examined the response to aversive odorants by tcf7l2zf55 mutants that develop with symmetric left-isomerized dHb and lack the vHb, but are viable to adulthood (Hüsken et al., 2014; Muncan et al., 2007). In agreement with the transformation of dHb identity, projections from OB mitral cells do not terminate in the right dHb of homozygous mutants nor are lratd2a-expressing neurons or their efferents to the vIPN detected (Figure 4A–D).
Following application of cadaverine, tcf7l2zf55 homozygous adults failed to exhibit the characteristic aversive behavior of their wild-type siblings (Figure 4E and F and Figure 4—figure supplement 1). Exposure to alarm substance also did not elicit a significant increase in swimming speed from baseline (1.13 ± 0.22 cm/s before and 1.89 ± 0.56 cm/s after) relative to WT siblings (2.84 ± 0.48 cm/s before and 4.88 ± 0.63 cm/s after, Figure 4G, I and J). Homozygous tcf7l2zf55 mutants tended to swim more slowly (Figure 4G) and spend more time in the top half of a novel test tank than wild-type adults, although the latter behavior was suppressed in the presence of alarm substance (Figure 4H). tcf7l2zf55 mutants are thus partly responsive to this fearful cue.
To further assess the role of lratd2a-expressing neurons in aversive olfactory processing, we looked at homozygous mutants of the brain-specific homeobox (bsx) gene, which develop right-isomerized dHb [(Schredelseker and Driever, 2018) and Figure 5A] and are viable to adulthood (Schredelseker and Driever, 2018). As might be expected when both dHb have right identity, equivalent populations of lratd2a-expressing neurons were found on both sides of the brain (Figure 5C). Instead of innervating only the right dHb as in controls, the axons of lhx2a:gap-YFP-labeled olfactory mitral cells terminated in the left and right dHb [(Dreosti et al., 2014) and Figure 5B], where the clusters of lratd2a neurons are situated (data not shown). Projections from the lratd2a dHb neurons coursed bilaterally through the left and right fasciculus retroflexus (FR) and innervated the same limited region of the ventral IPN (Figure 5C).
To measure the reaction to cadaverine in bsxm1376/m1376 adults with bilaterally symmetric lratd2a neurons, we counted the number of cells expressing fos in the dHb and found an increase in the left nucleus compared to heterozygous siblings (Figure 5D–E). Despite the symmetric activation of dHb neurons, bsx homozygotes and heterozygotes both showed reduced responsiveness to cadaverine (Figure 5F and Figure 5—figure supplement 1). Overall, homozygous mutants were slower swimmers than heterozygotes (Figure 5G–J); however, after exposure to alarm substance, their swimming speed relative to baseline was twofold faster than that of their heterozygous siblings (Figure 5G), indicative of an enhanced response to this aversive cue.
Discussion
From worms to humans, stimuli including odors are differently perceived by left and right sensory organs to elicit distinct responses (Güntürkün and Ocklenburg, 2017; Güntürkün et al., 2020). Honeybees, for example, show an enhanced performance in olfactory learning when their right antenna is trained to odors (Guo et al., 2016; Letzkus et al., 2006; Rogers and Vallortigara, 2008). In mice, over one third of mitral/tufted cells were found to be interconnected between the ipsilateral and contralateral olfactory bulbs for sharing of odor information received separately from each nostril, and for coordinated perception (Grobman et al., 2018). The zebrafish provides a notable example of a lateralized olfactory pathway, with the discovery of a subset of bilateral mitral cells that project to the dorsal habenulae but terminate only at the right nucleus (Miyasaka et al., 2014; Miyasaka et al., 2009). This finding prompted us to ask what is different about the post-synaptic dHb neurons that receive this olfactory input and what function does this asymmetric pathway serve.
Aversive olfactory cues activate identified neurons in the right dHb
We previously showed that the olfactory mitral cells that express lhx2a and are located in medio-dorsal and ventro-medial bilateral glomerular clusters (Miyasaka et al., 2014; Miyasaka et al., 2009) and project their axons to a subregion of the right dHb where the lratd2a gene is transcribed (deCarvalho et al., 2013). From transgenic labeling with membrane-tagged fluorescent proteins and synaptophysin, we now confirm that the lhx2a olfactory neurons precisely terminate at a cluster of lratd2a/slc5a7a expressing cholinergic neurons present in the right dHb.
Through calcium imaging using a genetically encoded calcium indicator, we validated that the right dHb is responsive when larval zebrafish are exposed to aversive odors such as cadaverine or chondroitin sulfate (Jetti et al., 2014; Krishnan et al., 2014), a component of alarm substance (Mathuru et al., 2012), more specifically, that the lratd2a-expressing neurons of the right dHb significantly respond to these aversive olfactory cues above their response to water alone. As has also been observed by others (Jesuthasan et al., 2021), application of vehicle alone, even when introduced slowly into a testing chamber, is sufficient to elicit a change in GCaMP fluorescence.
In adults, we used fos expression as a measure of neuronal activation and showed that transcripts colocalized to lratd2a-expressing cells. Interestingly, cadaverine predominantly activated neurons in the right dHb in larvae and adults, whereas neurons responsive to alarm substance were detected in both the left and right dHb nuclei of adult zebrafish. Different types of olfactory cues activate distinct glomeruli in the OB (Friedrich and Korsching, 1997; Yoshihara, 2014), and consistent with the prior studies, we observed that, in adults, cadaverine significantly increased fos expression in the mdG and dG regions, the location of lhx2a neurons that project to the right dHb. By contrast, alarm substance predominantly activated neurons in the lG and dlG regions of the OB that innervate the telencephalon and posterior tuberculum (Miyasaka et al., 2014; Miyasaka et al., 2009), suggesting that both dHb receive input via this route rather than through direct olfactory connections. Indeed, we found that more neurons reacted to alarm substance than cadaverine throughout the brain, including in the Dp, Vv and thalamic areas (data not shown).
In previous experiments (deCarvalho et al., 2013), we did not detect activated neurons in the right dHb of adult zebrafish following exposure to cadaverine or alarm substance. Several factors could account for the difference from the earlier study: we now have the transgenic tools to examine lratd2a neurons directly, we used higher concentrations of cadaverine and alarm substance and, in contrast to delivering odorants to groups of zebrafish, we tested the neuronal response in individual adults.
It has been suggested that lateralized olfactory and visual functions of the dHb are more prominent early in development and less so at later stages (Fore et al., 2020). However, the presence of lratd2a-expressing neurons in the right dHb and their preferential response to cadaverine from larval to adult stages supports the persistence of lateralized activity and illustrates the value of examining defined neuronal populations.
Right dHb neurons mediate aversive behavioral responses
As a group, zebrafish in which the synaptic activity of lratd2a neurons in the right dHb was inhibited by BoTxBLC-GFP did not exhibit repulsion to cadaverine, although some individuals showed a mild aversive response that was not sustained relative to sibling controls. This subtle behavioral effect suggests that the entire population of lratd2a neurons may not be effectively inactivated in all animals. The described intersectional approach to suppress neuronal activity may be incomplete, as evidenced by the observed variability in BoTxBLC-GFP labeling between individual larvae. Additionally, compensatory mechanisms occurring between the onset of Cre expression in slc5a7a neurons (at 3 dpf) and adulthood could allow for partial recovery of the response to repulsive cues. It is also possible that the lratd2a neurons are a heterogeneous population, with different subsets controlling the magnitude of the aversive response or the duration. More finely tuned techniques for temporal or spatial regulation of neuronal inactivation would help resolve these issues.
Despite both being aversive cues (Hussain et al., 2013; Mathuru et al., 2012), cadaverine and alarm substance elicit different behavioral responses by adult zebrafish. Control fish show active repulsion to cadaverine for the first to 2–4 min of a 5 min testing period, whereas alarm substance triggers immediate erratic behavior such as rapid swimming and darting that is typically followed by prolonged freezing (Hussain et al., 2013; Mathuru et al., 2012). Therefore, it is not necessarily expected that the same neuronal populations will mediate the response to both substances.
Perturbation of the lratd2a-expressing right dHb neurons either selectively by BoTxBLC-mediated synaptic inactivation, or in tcf7l2zf55 homozygous mutants that completely lack them, reduced aversion to cadaverine, either in the length or degree of the response. In contrast to juveniles or adults with BoTxBLC inactivated neurons that displayed a similar response to alarm substance as controls, tcf7l2zf55 mutants, showed no difference in their swimming behavior before and after its addition. One explanation is that many regions throughout the brain are likely involved in directing the complex repertoire of behaviors elicited by alarm substance and inactivation of lratd2a neurons in the habenular region alone is insufficient to weaken the overall response. Furthermore, the tcf7l2zf55 mutation could disrupt other brain regions that regulate behaviors elicited by alarm substance since the tfc7l2 gene is expressed in neurons throughout the brain, including the anterior tectum, dorsal thalamus and the hindbrain (Young et al., 2002).
Similar to tcf7l2zf55, the bsxm1376 mutation is pleiotropic resulting in right-isomerization of the dHb due to the absence of the parapineal (Schredelseker and Driever, 2018), and also the loss of the terminal tuberal hypothalamus, mammillary hypothalamic regions and secondary prosencephalon (Schredelseker et al., 2020). We did not observe enhanced or prolonged aversion to cadaverine in bsxm1376 homozygotes relative to controls. However, although homozygous mutants did show a hyperactive response to alarm substance, we cannot discount the involvement of other affected brain regions. Albeit technically challenging in adults, a more selective test such as optogenetic activation of only the lratd2a dHb neurons in wild-type and mutant zebrafish could help resolve their contribution to the alarm response.
The identification of a subset of neurons in the right dHb that receive olfactory input and terminate their axons at a defined subregion of the ventral IPN lays the groundwork for tracing an entire pathway from olfactory receptors to the neurons directing the appropriate behavioral response. The midline IPN has been morphologically defined into subregions (deCarvalho et al., 2014 ; Lima et al., 2017; Quina et al., 2017), but their connectivity and functional properties are not well studied. Recent work has begun to assign different functions to given subregions, such as the role of the rostral IPN in nicotine aversion (Morton et al., 2018; Quina et al., 2017). Neurons in the ventral IPN project to the raphe nucleus (Agetsuma et al., 2010; Lima et al., 2017), but the precise identity of raphe neurons that are innervated by the lratd2a-expressing dHb neurons remains to be determined. Transcriptional profiling of the IPN should yield useful information on its diverse neuronal populations and likely lead to the identification of the relevant post-synaptic targets in the ventral IPN and their efferent connections. Elaboration of this pathway may also help explain the advantage of lateralization in the processing of aversive information. It has been argued, for instance, that the antennal specialization to aversive odors in bees is correlated with directed turning away from the stimulus and escape (Rogers and Vallortigara, 2019). Directional turning has also been observed in larval zebrafish (Horstick et al., 2020), but whether it is correlated with laterality of the Hb-IPN pathway is unclear. Beyond olfaction, left-right asymmetry appears to be a more general feature of stress-inducing, aversive responses as demonstrated for the rat ventral hippocampus (Sakaguchi and Sakurai, 2017) and human pre-frontal cortex, where heightened anxiety also activates more neurons on the right than on the left (Avram et al., 2010; Ocklenburg et al., 2016).
Materials and methods
Zebrafish
Zebrafish were maintained at 27 °C under a 14:10 hr light/dark cycle in a recirculating system with dechlorinated water (system water). The AB wild-type strain (Walker, 1998), transgenic lines Tg(lratd2a:QF2)c601, Tg(slc5a7a:Cre)c662, Tg(Xla.Tubb2:QF2;he1.1:mCherry)c663, Tg(QUAS:GCaMP6f)c587, Tg(QUAS:GFP)c403 (Subedi et al., 2014), Tg(QUAS:mApple-CAAX;he1.1:mCherry)c636, Tg(QUAS:loxP-mCherry-loxP-GFP-CAAX)c679, and Tg(QUAS:loxP-mCherry-loxP-BoTxBLC-GFP)c674, Tg(–10lhx2a:gap-EYFP)zf177 (formally known as Tg(lhx2a:gap-YFP)) (Miyasaka et al., 2009), Tg(lhx2a:syp-GFP)zf186 (Miyasaka et al., 2009), and mutant strains tcf7l2zf55 (Muncan et al., 2007) and bsxm1376 (Schredelseker and Driever, 2018) were used. For imaging, embryos and larvae were transferred to system water containing 0.003% phenylthiourea (PTU) to inhibit melanin pigmentation. All zebrafish protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of the Carnegie Institution for Science [Protocol #122] or Dartmouth College [Protocol #00002253(m3a)].
Generation of transgenic lines by Tol2 transgenesis
Request a detailed protocolThe MultiSite Gateway-based construction kit (Kwan et al., 2007) was used to create transgenic constructs for Tol2 transposition. A 16 bp QUAS sequence (Potter et al., 2010), was cloned into the 5’ entry vector (pDONRP4-P1R, #219 of Tol2kit v1.2) via a BP reaction (11789020, Thermo Fisher Scientific). Middle entry vectors (pDONR221, #218 of Tol2kit v1.2 Kwan et al., 2007) were generated for QF2, mApple-CAAX, loxP-mCherry-stop-loxP and GCaMP6f. Sequences corresponding to the SV40 poly A tail, the SV40 poly A tail followed by a secondary marker consisting of the zebrafish hatching enzyme one promoter (Xie et al., 2012) driving mCherry, or to BoTxBLC-GFP (Lal et al., 2018; Sternberg et al., 2016; Zhang et al., 2017) were cloned into the 3’ entry vector (pDONRP2R-P3, #220 of Tol2kit v1.2 Kwan et al., 2007). Final constructs were created using an LR reaction (11791020, Thermo Fisher Scientific) into a Tol2 destination vector (pDestTol2pA2, #394 of the Tol2kit v1.2 Kwan et al., 2007; Supplementary file 1).
To produce Tol2 transposase mRNA, pCS-zT2TP was digested by NotI and RNA synthesized using the mMESSAGE mMACHINE SP6 Transcription Kit (AM1340, Thermo Fisher Scientific). RNA was purified by phenol/chloroform-isoamyl extraction, followed by chloroform extraction and isopropanol precipitation (Suster et al., 2011). A solution containing QF2/QUAS plasmid DNA (~25 ng/μl), transposase mRNA (~25 ng/μl) and phenol red (0.5%) was microinjected into one-cell stage zebrafish embryos, which were raised to adulthood. To identify transgenic founders, F0 adult fish were outcrossed to AB and embryos were assessed for the presence of the secondary marker by screening for mCherry labeling of hatching gland cells after 24 hr post fertilization (hpf) and raised to adulthood.
Generation of transgenic lines by genome editing
Request a detailed protocolFor generating transgenic lines at targeted sites, we performed CRISPR/Cas9-mediated genome editing using the method of Kimura et al., 2014, which relies on homology-independent repair of double-strand breaks for integration of donor DNA. To construct the donor DNA, we combined GFP bait sequences (Gbait) and the hsp70 promoter fragment (Kimura et al., 2014), with a QF2 sequence, which contains the DNA binding and transcriptional activation domains of the QF transcription factor of Neurospora crassa (Ghosh and Halpern, 2016; Subedi et al., 2014). The Gbait-hsp70 sequence was amplified with forward 5’- GGCGAGGGCGATGCCACCTACGG-3’ and reverse 5’- CCGCGGCAAGAAACTGCAATAAAAAAAAC-3’ primers, using Gbait-hsp70:Gal4 donor DNA (Kimura et al., 2014). QF2 sequence was amplified with forward 5’- ACTAGTATGCCACCCAAGCGCAAAACGC-3’ and reverse 5’- CTGCAGCAACTATGTATAATAAAGTTGAAA-3’ primers, using pDEST:QF2 template DNA. Subsequently, the Gbait-hsp70 fragment and QF2 fragment were independently inserted into pGEM T-easy (A1360, Promega) and subsequently combined into one vector by SacII digestion and ligation (Addgene, plasmid #122563). The Cre sequence was amplified using pCR8GW-Cre-pA-FRT-kan-FRT as template DNA (Suster et al., 2011) (forward 5’-ACTAGTGCCACCATGGCCAATTTACTG-3’ and reverse 5’-CTGCAGGGACAAACCACAACTAGA-3’ primers), and inserted into pGEM T-easy. The Gbait-hsp70 fragment was subcloned into the Cre vector by SacII digestion and ligation (Addgene, plasmid #122562).
Production of sgRNAs and Cas9 RNA was performed as described previously (Hwang et al., 2013; Jao et al., 2013). Potential sgRNAs were designed using Zifit (Sander et al., 2010). Pairs of synthetic oligonucleotides (lratd2a sense, 5’-TAGGACTGGACACCGAAGAAGA-3’; lratd2a anti-sense, 5’-AAACTCTTCTTCGGTGTCCAGT-3’; slc5a7a sense, 5’-TAGGCTCTTTGTGCACTGTTGG-3’; slc5a7a anti-sense, 5’-AAACCCAACAGTGCACAAAGAG-3’), 5’-TAGG-N18-3’ and 5’-AAAC-N18-3’, were annealed and inserted at the BsaI site of the pDR274 vector (Addgene, plasmid #42250). To make sgRNA and Cas9 mRNA, template DNA for sgRNAs and pT3TS nCas9n (Addgene, plasmid #46757) were digested by DraI and XbaI, respectively. The MAXIscript T7 Transcription Kit (AM1312, Thermo Fisher Scientific) was used for synthesis of sgRNAs from linearized DNA template and the mMESSAGE mMACHINE T3 Transcription Kit (AM1348, Thermo Fisher Scientific) for synthesis of Cas9 RNA. RNA was purified by phenol/chloroform and precipitated by isopropanol.
A solution containing sgRNA for the targeted gene (~50 ng/μl), sgRNA (~50 ng/μl) to linearize donor plasmids at the Gbait site (Auer et al., 2014; Kimura et al., 2014), the Gbait-hsp70-QF2-pA and Gbait-hsp70-Cre-pA (~50 ng/μl) plasmids, Cas9 mRNA (~500 ng/μl), and phenol red (0.5%) was microinjected into one-cell stage embryos. To verify integration of donor DNA at the target locus, PCR was performed using primers that correspond to sequences flanking the integration site and within the donor plasmid (hsp70 reverse, 5’-TCAAGTCGCTTCTCTTCGGT-3’). (For lratd2a, the forward primer is 5’-CTGCTGAAGTGGCATTTATGGGC-3’ and the reverse primer is 5’-CCTGGAAGTCCCCGACATAC-3’; for slc5a7a the forward primer is 5’-CACATCTCTCTGACGTCCATC-3’ and the reverse is 5’-GTTGCTGCGCAGGACTTAAAA-3’). Sequence analysis of PCR products confirmed integration at the targeted sites.
RNA in situ hybridization
Request a detailed protocolWhole-mount in situ hybridization was performed as previously described (deCarvalho et al., 2014; Gamse et al., 2002). In brief, larvae and dissected brains were fixed in 4% paraformaldehyde (P6148, Sigma-Aldrich) in 1 X PBS (phosphate-buffered saline) at 4 °C overnight. To synthesize RNA probes, the following restriction enzymes and RNA polymerases were used: lratd2a (BamHI/T7), fos (NotI/SP6), slc5a7a (NotI/SP6), kctd12.1 (EcoRI/T7) (deCarvalho et al., 2013; Hong et al., 2013). Probes were labeled with UTP-digoxigenin (11093274910, Roche) and samples incubated at 70 °C in hybridization solution containing 50% formamide. Hybridized probes were detected using alkaline phosphatase-conjugated antibodies (Anti-Digoxigenin-AP, #11093274910, and Anti-Fluorescein-AP, #11426338910, Sigma-Aldrich) and visualized by staining with 4-nitro blue tetrazolium (NBT, #11383213001, Roche), 5-bromo-4-chloro-3-indolyl-phosphate (BCIP, #11383221001, Roche) and 2-(4-Iodophenyl)–3-(4-nitrophenyl)–5-phenyltetrazolium Chloride (INT, #I00671G, Fisher Scientific).
Preparation of odorants
Request a detailed protocolAlarm substance was freshly prepared on the day of testing. Adult zebrafish (6 female and six male) were anesthetized in 0.02% tricaine (E10521, Ethyl 3-aminobenzoate methanesulfonate; Sigma-Aldrich). Shallow lesions were made on the skin (10 on each side) using a fresh razor blade and single fish were consecutively immersed in a beaker containing distilled water (25 ml for fos experiment and 50 ml for behavioral analyses) for 30 s at 4 °C. The solution was filtered using a 0.2 μm filter (565–0020, ThermoFisher Scientific) and stored at 4 °C until used (Mathuru et al., 2012). The cadaverine (#33211, final concentration 100 μM) and chondroitin sulfate (#c4384, 100 μg/ml) were purchased from Sigma-Aldrich, and stock solutions prepared in distilled water.
Calcium imaging in larval zebrafish
Request a detailed protocolTo monitor activity within the same population of lratd2a neurons in response to odorants or to a vehicle control, calcium signaling was performed on seven dpf Tg(lratd2a:QF2)c644; Tg(QUAS:GCaMP6f)c587 larvae. Larvae were paralyzed by immersion in α-bungarotoxin (20 µl of 1 mg/ml solution in system water, B1601, ThermoFisher Scientific) Duboué et al., 2017; Severi et al., 2014 followed by washing in fresh system water. Individuals were embedded in 1.5% low melting agarose (SeaPlaque Agarose, Lonza) in a petri dish (60 mm) with a custom-designed mold. After solidification, the agarose around the nostrils was carefully removed with forceps for access to odorants, and the individual immersed in fresh system water. The dish was placed under a 20 X (NA = 0.5) water immersion objective on a Zeiss LSM 980 confocal microscope. Images were acquired in XYT acquisition mode with 512 × 200 pixel resolution at a rate of 2 Hz and digitized eight bit from single focal planes. The vehicle control (deionized water) followed by the odorant solution (0.2 ml) was slowly expelled through plastic tubing (Tygon R-3606; 0.8 mm ID, 2.4 mm OD) attached to a 1 ml syringe (BD 309659) and with the other end positioned in front of an individual’s face. Two sequential applications of each were spaced by 1 min intervals. To calculate fluorescence intensity, regions of interest (ROI) were manually drawn around each cell in the average focal plane with the polygon tool and ROI manage in Fiji (Schindelin et al., 2012). To normalize calcium activity for each neuron to baseline fluorescence, the fractional change in fluorescence (ΔF/F) was calculated before the application of odorants (average of 320 frames from each neuron), according to the formula F = (Fi-F0)/F0, where Fi is the fluorescence intensity at a single time point and F0 is the baseline fluorescence. All data and images were analyzed using custom programs in MATLAB (MathWorks, version 7.3), Excel software and GraphPad Prism9.
Assay of Fos expression in adult zebrafish
Request a detailed protocolIndividual adult zebrafish (7–9 months old) were placed in a tank with 1 L system water and acclimated for at least 1 hr prior to odorant exposure. Each odorant solution (1 ml) was gently pipetted into the tank water and the fish was kept there as the odorant diffused. After 30 min, the fish was sacrificed in an ice water slurry, and the brain dissected out and fixed in 4% paraformaldehyde in 1 X PBS overnight at 4 °C. Fixed brains were embedded in 4% low melting agarose in 1 X PBS and sectioned at 50 μm (for juvenile brains) or 70 μm (for adult brains) using a vibratome (VT1000S, Leica Biosystems, Inc). For more precise counting of fos expressing cells in adult brains, habenular sections were 35 μm thick. Sections were covered in 50% glycerol in 1 X PBS under coverslips. Bright-field images were captured with a Zeiss AxioCam HRc camera mounted on a Zeiss Axioskop. A Leica SP5 confocal microscope was used for fluorescent images. Data from fos RNA in situ hybridization experiments were quantified using ImageJ/Fiji software (Schindelin et al., 2012).
Behavioral assays
Request a detailed protocolBehavioral assays were performed using 5–7 week-old juvenile zebrafish and adults that were of 4–8 months of age. Responses to odorants were measured between 10:00 a.m. and 4:00 p.m. and fish were starved for 1 day prior to testing (Koide et al., 2009). Individual adults were placed in a 1.5 L test tank (Aquatic Habitats) in 1 L of system water and allowed to acclimate for at least 1 hr. For experiments with juveniles, individuals were acclimated to the behavior room for 1 hr, gently netted into the test tank (20 × 9 x 8.3 cm, 1.5 L mating cage) containing 0.6 L system water and maintained there for 5 min prior to testing. Swimming activity was recorded for 5 min (4 min for juveniles) before and after the application of odorants. Odorants (2 ml for adults, 1 ml for juveniles) were slowly expelled through plastic tubing (Tygon R-3606; 0.8 mm ID, 2.4 mm OD) attached on one end to a 3 ml syringe (BD 309657) and on the other positioned at one end of the test tank. A preference index was calculated using the formula: Preference to odorant = [(Total time spent in the half of the tank where odorant was delivered) − (Total time spent in the other half of the tank)]/Total time (Koide et al., 2009; Wakisaka et al., 2017). A preference index of –1 indicates that the position of an individual fish is maximally distanced from the application site of the odorant.
Quantification and statistical analyses
Request a detailed protocolData were collected using custom written scripts in MATLAB (MathWorks). Statistical analyses were performed using GraphPad Prism9 (GraphPad Software, Inc). The Wilcoxon signed-rank test was used to compare a set of matched samples. An unpaired t-test was used to compare two independent data sets and two-way ANOVA followed by Bonferroni’s post hoc test was used to compare datasets greater than two. All tests were two-tailed and results depicted as non-significant (ns, p> 0.05) or significant (*, p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001).
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting file. Source data files have been provided for Figures 1-5. Behavioral analyses were performed using custom written scripts in MATLAB and uploaded as Source code Files.
References
-
The habenula is crucial for experience-dependent modification of fear responses in zebrafishNature Neuroscience 13:1354–1356.https://doi.org/10.1038/nn.2654
-
Identification of the zebrafish ventral habenula as a homolog of the mammalian lateral habenulaThe Journal of Neuroscience 30:1566–1574.https://doi.org/10.1523/JNEUROSCI.3690-09.2010
-
Frontal EEG activation asymmetry reflects cognitive biases in anxiety: evidence from an emotional face Stroop taskApplied Psychophysiology and Biofeedback 35:285–292.https://doi.org/10.1007/s10484-010-9138-6
-
Mating Suppresses Alarm Response in ZebrafishCurrent Biology 29:2541–2546.https://doi.org/10.1016/j.cub.2019.06.047
-
Left Habenular Activity Attenuates Fear Responses in Larval ZebrafishCurrent Biology 27:2154–2162.https://doi.org/10.1016/j.cub.2017.06.017
-
BookTranscriptional regulation using the Q system in transgenic zebrafishIn: Detrich M, Leonard IZ, editors. Methods in Cell Biology. Academic Press. pp. 205–218.https://doi.org/10.1016/bs.mcb.2016.05.001
-
Brain Lateralization: A Comparative PerspectivePhysiological Reviews 100:1019–1063.https://doi.org/10.1152/physrev.00006.2019
-
The habenula: from stress evasion to value-based decision-makingNature Reviews. Neuroscience 11:503–513.https://doi.org/10.1038/nrn2866
-
Molecular and cellular determinants of motor asymmetry in zebrafishNature Communications 11:1170.https://doi.org/10.1038/s41467-020-14965-y
-
Efficient genome editing in zebrafish using a CRISPR-Cas systemNature Biotechnology 31:227–229.https://doi.org/10.1038/nbt.2501
-
The alarm response in zebrafish: innate fear in a vertebrate genetic modelJournal of Neurogenetics 22:211–228.https://doi.org/10.1080/01677060802298475
-
Neural correlates of state transitions elicited by a chemosensory danger cueProgress in Neuro-Psychopharmacology & Biological Psychiatry 111:110110.https://doi.org/10.1016/j.pnpbp.2020.110110
-
The right dorsal habenula limits attraction to an odor in zebrafishCurrent Biology 24:1167–1175.https://doi.org/10.1016/j.cub.2014.03.073
-
The Tol2kit: a multisite gateway-based construction kit for Tol2 transposon transgenesis constructsDevelopmental Dynamics 236:3088–3099.https://doi.org/10.1002/dvdy.21343
-
Lateralization of olfaction in the honeybee Apis melliferaCurrent Biology 16:1471–1476.https://doi.org/10.1016/j.cub.2006.05.060
-
Afferent and efferent connections of the interpeduncular nucleus with special reference to circuits involving the habenula and raphe nucleiThe Journal of Comparative Neurology 525:2411–2442.https://doi.org/10.1002/cne.24217
-
Stress and laterality - The comparative perspectivePhysiology & Behavior 164:321–329.https://doi.org/10.1016/j.physbeh.2016.06.020
-
Specific connections of the interpeduncular subnuclei reveal distinct components of the habenulopeduncular pathwayThe Journal of Comparative Neurology 525:2632–2656.https://doi.org/10.1002/cne.24221
-
Complementary Specializations of the Left and Right Sides of the Honeybee BrainFrontiers in Psychology 10:280.https://doi.org/10.3389/fpsyg.2019.00280
-
Left-right functional asymmetry of ventral hippocampus depends on aversiveness of situationsBehavioural Brain Research 325:25–33.https://doi.org/10.1016/j.bbr.2017.02.028
-
ZiFiT (Zinc Finger Targeter): an updated zinc finger engineering toolNucleic Acids Research 38:W462–W468.https://doi.org/10.1093/nar/gkq319
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Bsx controls pineal complex developmentDevelopment 145:dev163477.https://doi.org/10.1242/dev.163477
-
BookOlfactory LateralizationIn: Rogers LJ, Vallortigara G, editors. Lateralized Brain Functions: Methods in Human and Non-Human Species. Springer. pp. 103–120.https://doi.org/10.1007/978-1-4939-6725-4_4
-
Transposon-mediated BAC transgenesis in zebrafishNature Protocols 6:1998–2021.https://doi.org/10.1038/nprot.2011.416
-
An Adenosine Receptor for Olfaction in FishCurrent Biology 27:1437–1447.https://doi.org/10.1016/j.cub.2017.04.014
-
BookChapter 3 Haploid Screens and Gamma-Ray MutagenesisIn: Detrich HW, Westerfield M, Zon LI, editors. Methods in Cell Biology. Academic Press. pp. 43–70.https://doi.org/10.1016/S0091-679X(08)61893-2
-
BookZebrafish Olfactory SystemIn: Mori K, editors. The Olfactory System: From Odor Molecules to Motivational Behaviors. Springer. pp. 71–96.https://doi.org/10.1007/978-4-431-54376-3_5
-
Expression and splice variant analysis of the zebrafish tcf4 transcription factorMechanisms of Development 117:269–273.https://doi.org/10.1016/s0925-4773(02)00180-6
Article and author information
Author details
Funding
National Institutes of Health (R01HD078220)
- Marnie E Halpern
National Institutes of Health (R37HD091280)
- Marnie E Halpern
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to Bethany Malskis for animal care and to Dr. Krishan Ariyasiri and Dr. Eugene Demidenko for advice on statistical analyses. We thank Dr. Shin-ichi Higashijima for the Gbait-hsp70:Gal4 donor plasmid, Dr. Wenbiao Chen for pT3TS nCas9n plasmid, Dr. Koichi Kawakami for the UAS:zBoTXBLC-GFP construct, Paul Krieg for the Xla.Tubb promoter construct, Dr. Claire Wyart for GCaMP6f plasmid and Dr. Wolfgang Driever and Dr. Tatjana Piotrowski for providing bsxm1376 and tcf7l2zf55 mutant zebrafish, respectively. This work was supported by NIH grants R37HD091280 and R01HD078220 to MEH.
Ethics
All zebrafish protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of the Carnegie Institution for Science [Protocol #122] or Dartmouth College [Protocol #00002253(m3a)].
Copyright
© 2021, Choi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,937
- views
-
- 253
- downloads
-
- 24
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Sour taste, which is elicited by low pH, may serve to help animals distinguish appetitive from potentially harmful food sources. In all species studied to date, the attractiveness of oral acids is contingent on concentration. Many carboxylic acids are attractive at ecologically relevant concentrations but become aversive beyond some maximal concentration. Recent work found that Drosophila ionotropic receptors IR25a and IR76b expressed by sweet-responsive gustatory receptor neurons (GRNs) in the labellum, a peripheral gustatory organ, mediate appetitive feeding behaviors toward dilute carboxylic acids. Here, we disclose the existence of pharyngeal sensors in Drosophila melanogaster that detect ingested carboxylic acids and are also involved in the appetitive responses to carboxylic acids. These pharyngeal sensors rely on IR51b, IR94a, and IR94h, together with IR25a and IR76b, to drive responses to carboxylic acids. We then demonstrate that optogenetic activation of either Ir94a+ or Ir94h+ GRNs promotes an appetitive feeding response, confirming their contributions to appetitive feeding behavior. Our discovery of internal pharyngeal sour taste receptors opens up new avenues for investigating the internal sensation of tastants in insects.
-
- Neuroscience
Time estimation is an essential prerequisite underlying various cognitive functions. Previous studies identified ‘sequential firing’ and ‘activity ramps’ as the primary neuron activity patterns in the medial frontal cortex (mPFC) that could convey information regarding time. However, the relationship between these patterns and the timing behavior has not been fully understood. In this study, we utilized in vivo calcium imaging of mPFC in rats performing a timing task. We observed cells that showed selective activation at trial start, end, or during the timing interval. By aligning long-term time-lapse datasets, we discovered that sequential patterns of time coding were stable over weeks, while cells coding for trial start or end showed constant dynamism. Furthermore, with a novel behavior design that allowed the animal to determine individual trial interval, we were able to demonstrate that real-time adjustment in the sequence procession speed closely tracked the trial-to-trial interval variations. And errors in the rats’ timing behavior can be primarily attributed to the premature ending of the time sequence. Together, our data suggest that sequential activity maybe a stable neural substrate that represents time under physiological conditions. Furthermore, our results imply the existence of a unique cell type in the mPFC that participates in the time-related sequences. Future characterization of this cell type could provide important insights in the neural mechanism of timing and related cognitive functions.