Chemo- and optogenetic activation of hypothalamic Foxb1-expressing neurons and their terminal endings in the rostral-dorsolateral PAG leads to tachypnea, bradycardia, and immobility
eLife assessment
This paper describes valuable results from studies investigating circuits in the brain that underlie behavioral responses in fearful situations. The authors identified a role for a class of neurons that are sufficient to cause these stereotyped behaviors including freezing behaviors. These solid studies increase our understanding of brain pathways regulating these types of behaviors.
https://doi.org/10.7554/eLife.86737.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Solid: Methods, data and analyses broadly support the claims with only minor weaknesses
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Foxb1 -expressing neurons occur in the dorsal premammillary nucleus (PMd) and further rostrally in the parvafox nucleus, a longitudinal cluster of neurons in the lateral hypothalamus of rodents. The descending projection of these Foxb1+ neurons end in the dorsolateral part of the periaqueductal gray (dlPAG). The functional role of the Foxb1+ neuronal subpopulation in the PMd and the parvafox nucleus remains elusive. In this study, the activity of the Foxb1+ neurons and of their terminal endings in the dlPAG in mice was selectively altered by employing chemo- and optogenetic tools. Our results show that in whole-body barometric plethysmography, hM3Dq-mediated, global Foxb1+ neuron excitation activates respiration. Time-resolved optogenetic gain-of-function manipulation of the terminal endings of Foxb1+ neurons in the rostral third of the dlPAG leads to abrupt immobility and bradycardia. Chemogenetic activation of Foxb1+ cell bodies and ChR2-mediated excitation of their axonal endings in the dlPAG led to a phenotypical presentation congruent with a ‘freezing-like’ situation during innate defensive behavior.
Introduction
Mapping the circuits that form the basis of specific behaviors is essential to understand the brain and is an essential task for the neuroscience of our time. The circuitries that underlie fear are important in the clinic and consequently have become a hotspot of neuroscientific research in the last few years (Tovote et al., 2015; Mobbs et al., 2020; Silva and McNaughton, 2019). Here, we have used chemogenetics and optogenetics to show that the activation of a specific, well delimited neuronal group in the hypothalamus leads to innate defensive behavior, including tachypnea, immobility, and bradycardia.
During hypothalamic development, two neuronal migratory streams expressing the transcription factor Foxb1 form the nuclei of the mamillary region as well as one nucleus of the lateral hypothalamus (LHA). One of these streams forms most of the mamillary body while the other forms the dorsal premammillary nucleus (PMd) and the Foxb1-expressing neurons of the parvafox nucleus (parvafoxFoxb1) in the LHA (Alvarez-Bolado et al., 2000). The ventrolateral part of the PMd (vlPMd) projects to the dorsolateral periaqueductal gray (dlPAG; Canteras and Swanson, 1992) while the parvafoxFoxb1 axon terminals occupy a region straddling the dlPAG and the lPAG (Bilella et al., 2016). Solid evidence indicates that the dlPAG is associated with innate defensive responses like ‘freezing’ and escape behavior (Canteras and Goto, 1999; Vianna et al., 2001; Bittencourt et al., 2004; Deng et al., 2016; Kunwar et al., 2015; Souza and Carobrez, 2016; Evans et al., 2018).
In addition to inputs from Foxb1+ cells of the hypothalamus, the dlPAG also receives afferences from neurons of the ventromedial hypothalamus (VMH) expressing steroidogenic factor 1 (SF1+). The optogenetic activation of their cell bodies in the VMH (Kunwar et al., 2015) and of their terminal endings in the dlPAG (Wang et al., 2015a) leads to immobility. As the stimulation intensity increases, immobility is followed by bursts of activity (running and jumping; Kunwar et al., 2015). This behavioral change is ascribed to the co-activation of another axon collateral of the SF1+ neurons, innervating the anterior hypothalamic nucleus (AHN; Wang et al., 2015a). As animals remain immobile, increasing the intensity of the stimulation of SF1+ cell bodies in the VMH leads to a change of the cardiovascular responses from tachycardia to bradycardia (Wang et al., 2015b).
The third afference to the dlPAG (Schenberg et al., 2005) stems from medially located glutamatergic neurons of the superior colliculus (mSC), which is responsive to innately aversive looming stimuli (Evans et al., 2018). mSC neurons integrate threat evidence and pass the information through a synaptic threshold to the glutamatergic dlPAG neurons that initiate escape (Evans et al., 2018).
We set out to establish the function of two distinct neuronal groups identifiable through Foxb1 expression, the parvafoxFoxb1 of the LHA and the PMdFoxb1 of the hypothalamic mamillary region, both projecting to specific regions of the PAG (Canteras and Swanson, 1992; Bilella et al., 2016). To that end, we used the specific expression of Cre recombinase in Foxb1tm1(cre)Gabo mice as a tool.
In our study (Figure 1), projection-unspecific chemogenetic activation of the Foxb1+ cell bodies of the PMd and of the parvafoxFoxb1 led to tachypnea, while selective optogenetic activation of theire terminal endings in the rostral part of the dlPAG provoked immobility, accompanied by bradycardia.
Results
Chemogenetic modulation of the Foxb1 neurons increases breaths per minute (BPM) without altering tidal volume (TV)
The activation of the hM3Dq (activating DREADD receptors) expressed in the parvafoxFoxb1 by intraperitoneal clozapine or CNO injection led to changes in several parameters of respiration consistent with an increased respiratory effort Figure 2; breaths per minute (BPM), F(2,26)=6.061, p=0.00691, g=; inspiratory time (IT), F(2,26) = 5.831, p=0.00809; total time (TT) F(2,26)=4.655, p=0.01973; minute volume adjusted for bodyweight (MVadjPerGram), F(2,26)=3.998, p=0.03265; peak inspiratory flow (PIFadj), F(2,26) = 4.284, p=0.03439; all reported results are obtained by a 3x3 mixed-design ANOVA for the ‘Condition’ factor, p values are Huynh-Feldt sphericity corrected. Namely, we observed a significant increase in three parameters within saline, clozapine and CNO injected hM3Dq animals: (1) BPM (Figure 2a–I; saline vs. clozapine: t(5) = –3.199, p=0.024, g = –1.380 (large); saline vs. CNO: t(5) = –2.916, p=0.033, g = –1.192(large); clozapine vs. CNO: t(5) = 3.069, p=0.028, g=0.226(small)), (2) MVadjPerGram: (Figure 2a–VI; saline vs. clozapine: t(5) = –2.790, p=0.038, g=–0.974(large); saline vs. CNO: t(5) = –2.600, p=0.048, g = –0.928(large)) and (3) PIFadj (Figure 2a–VII; saline vs. clozapine: t(5) = –2.726, p=0.041, g = –1.217(large); all reported results are post-hoc two-tailed paired student’s t-Tests; g values represent the effect size according Hedge’s g corrected for paired data, the magnitude of the effect size is provided in parentheses [i.e. negligible, small, medium, or large]). Remarkably, the tidal volume normalized for bodyweight (TVadjPerGram; Figure 2a-IV) was not altered by the intervention (3x3 mixed-design ANOVA: F(2,26) = 1.165, p=0.322). Furthermore, since minute volume equals the tidal volume multiplied by BPM, it follows that, the principal factor explaining the changes in MVadjPerGram was the increased BPM. The inspiratory time (IT, Figure 2a–II; ; saline vs. clozapine: t(5) = 3.296, p=0.022, g=1.170[large]) and the total respiratory time (TT, Figure 2a-IV; saline vs. clozapine:: t(5) = 2.683, p=0.043, g=1.140[large]) decreased, while however the expiratory time (Figure 2a–III) remained unaffected (ET; mixed-design ANOVA: F(2,26) = 1.656, p=0.217). Therefore, the increase in BPM was achieved by selectively shortening the IT through an increase in PIFadj, while expiratory parameters (e.g. PEFadj (Figure 2a–VIII)) did not achieve statistical significance level and therefore do not seem to have contributed significantly to the increase in BPM.
It is important to note that none of the parameters were altered to a level of statistical significance neither in the control group (DREADD_neg) nor in the inhibitory DREADD group (hM4Di). This indicates that the effect was caused by neuronal excitation of the parvafoxFoxb1 and not by an inherent effect of the injected substances (i.e. clozapine or CNO). This conclusion is further strengthened by the increase of c-Fos immunoreactivity evident in the parvafox nucleus of hM3Dq expressing mice, but not in those of hM4Di expressing mice nor in control mice (Figure 2—figure supplement 1).
To test for a potential bias of respiratory recordings by alterations in gross locomotor activity, an open-field test was performed (Figure 2c and d). Global hM3Dq activation of the parvafoxFoxb1 decreased the distance moved to a level of statistical significance, however, with only a small effect size (t(4) = 3.774, p=0.02, g=0.225(small)). The time spent in an immobile state was not affected by chemogenetic activation of the parvafoxFoxb1 (t(4) = 2.257, p=0.087). Activation of the inhibitory hM4Di receptors did not affect gross locomotor activity distance moved: t(4) = –0.697, p=0.524; time spent in immobile state: t(4) = 0.695, p=0.525. Injections of clozapine in control mice (DREADD_neg) did not alter distance moved (t(3) = –2.708, p=0.073), but lead to a statistically significant decrease of time spent in an immobile state, however, with a negligible effect size t(3) = 4.336, p=0.027, g=0.081 (negligible). Thus, despite reaching statistical significance for distance moved in hM3Dq animals and for time spent in an immobile state in DREADD_neg animals, the magnitude of the effect size of chemogenetic modulation of the parvafoxFoxb1 on gross locomotor activity is only small and negligible, respectively, and is unlikely to explain the large size of the effects observed in the alterations of respiratory patterns.
Optogenetic modulation of Foxb1 terminals in the dlPAG induces immobility
Respiratory results aside (see previous paragraph), the most prominent effect observed in the optogenetic experiments was the immobility (‘freezing behavior’) displayed by the group of ChR2-expressing mice during photoactivation of the hypothalamic Foxb1+ axonal projections to the rostral part of the dlPAG (Figure 1; Bregma –3.40/–4.04; Figure 3). These mice exhibited a short-latency immobility response immediately after onset of the photostimulation with blue light (7–15 mW, burst of 500ms duration with an intraburst frequency of 30 Hz (5ms pulse duration) and an interburst interval of 500ms). A decrease of locomotor activity was observed in 13 out of 16 ChR2-expressing mice (Figure 3a–c), of which eight mice showed almost complete absence of locomotion. Three of these mice were still capable of moving their heads as a sign of attentive behavior towards its surroundings (158/19, 160/19, and 35B-20), while all four limbs remained largely immobile (Video 1). Five other mice (162/19, 106/21–10, 34/21–7, 34/21–10, and 35E-20) displayed an even greater effect and remained completely immobile (including head activity) during almost the entire duration of the LEDon period (Video 2). No activity bursts were observed, neither during nor after the photostimulation period. Increasing the intensity of the stimulation from 70 to 222 mW (measured at patchcord) using a laser beam did not trigger motor activities.
To quantify the locomotive behavior induced by optogenetic modulation, we performed open field experiments, where mice were recorded for 3 min without photo-stimulation (BL) and 3 min with photo-stimulation (Stim).
Since we observed different degrees of immobility displayed by ChR2-expressing mice, we allocated each mouse to either an ‘OnTarget_antPAG’ or an ‘OffTarget’ group based on the histologically confirmed optic fiber placement. OnTarget_antPAG animals had the tip of the optic fiber implant located above the dlPAG at an anterior-posterior level AP-4.04mm (from bregma) or proximal. The OffTarget group contains animals with fiber tips located below (i.e. ventral to) the dlPAG and/or located more distal than AP –4.04 mm. Two-tailed paired Welch’s t-Tests revealed significant effects of photoactivation in the following parameters in OnTarget_antPAG mice with ChR2 expression in the ventral-posterior hypothalamus (Figure 3a–h): Distance moved (t(8) = 5.934, p<0.001, g=2.506 [large]), speed.moving: (t(8) = 6.003, p<0.001, g=2.612 [large]), time in immobile state: (t(8) = –4.704, p=0.002, g = –1.976 [large]), time in center (t(8) = 4.108, p=0.003, g=1.178 [large]) and time in periphery (t(8) = –2.841, p=0.022, g=–0.903 [large]). All these parameters had a trend in the same direction in the OffTarget group; however, none besides speed.moving were altered to a level of statistical significance (Figure 3d–h): Distance moved (t(7) = 1.725, p=0.128), speed.moving: (t(7) = 2.895, p=0.023, g=0.779 [medium]), time in immobile state: t(7)=–0.996, p=0.352, time in periphery (t(7) = –0.899, p=0.398). Besides speed.moving: (t(3) = 0.096, p=0.929), photoinhibition in ArchT3.0 OnTarget_antPAG mice consistently changed locomotor parameters in a direction opposite to the ChR2 mice (Figure 3a–h): Distance moved (t(3) = –3.176, p=0.050, g = –0.832 [large]), time in immobile state: (t(3) = 5.369, p=0.013, g = –1.712 [large]), time in center (t(3) = –4.590, p=0.019, g = –0.885 [large]) and time in periphery (t(3) = 1.367, p=0.265).
A comment about the specificity of the optogenetic results: the large volume of virus injected (200 nl) and the highly penetrant and strongly expressing viral serotype of the vector used (AAV5) led to expression of ChR2 not only (as intended) in the parvafoxFoxb1 neurons but also in other Foxb1-expressing hypothalamic neurons in the neighboring mammillary region.
However, of all the mammillary nuclei only the dorsal premammillary nuclei (PMd) send axons to the rostral dlPAG (Vertes, 1992; Meller and Dennis, 1986; Canteras and Swanson, 1992), where the tip of the glass fibers for the optogenetic activation were implanted. In our experiments, the terminals from the PMdFoxb1 were therefore co-activated with those of the parvafoxFoxb1. Projections to the rostral dlPAG from other mamillary neurons are very scarce (Allen and Hopkins, 1990; Shen, 1983; Beart et al., 1990; Cruce, 1977; Canteras et al., 1992).
In summary, optogenetic activation of Foxb1+ terminals originating from the hypothalamic parvafox nucleus and/or PMd and projecting to the anterior PAG induces a state of immobility and hypoactivity. Such an activation namely decreases the distance moved, the speed during locomotor periods (speed.moving), the time spent in the center region of the arena, and increases the time spent in an immobile state and the time spent in the periphery of the arena. Stimulation at more posterior positions along the PAG and/or ventral to the dlPAG columns are inefficient in inducing this behavior. The observed effect is bidirectional, in the sense that inhibition of the same terminals in the anterior PAG induces increased distance moved, and time spent in the center of the arena, while decreasing the time spent in an immobile state.
scRNA seq dataset reveals distinct Foxb1 expression in the PMd
Our results to this point indicated that a population of Foxb1-expressing neurons in the PMd induces immobility in mice. In contrast to this finding, a cholecystokinin (Cck)-expressing population of neurons in the PMd has been shown to induce escape behavior in mice (Wang et al., 2021). We hypothesized that the Cck-expressing and the Foxb1-expressing PMd neurons are distinct, separate neuronal groups regulating opposite behaviors. To validate this hypothesis, we performed a reanalysis of previously published single-cell RNA sequencing datasets focusing on the murine posterior ventral hypothalamus (Mickelsen et al., 2020; Figure 4).
After quality control, normalization, and integration of all four datasets (from two female and two male mice), k-nearest neighbor clustering resulted in 24 distinct cell types detected in the ventral-posterior hypothalamus (Figure 4a). Plotting the expression levels of a set of PMd markers identified by Mickelsen et al. (i.e. Cck, Foxb1, Synpr, Dlk1, Ebf3, and Stxbp6) onto the UMAP plot identified cluster 9 as the PMd (Figure 4b). The cluster identity was further confirmed by calculating differentially-expressed genes between the PMd cluster and the other Cck- and Foxb1-expressing clusters (i.e. clusters 5, 7, and 8) and qualitative comparison of expression patterns from in situ hybridization data provided by the Allen Brain Atlas (Figure 4c and Figure 4—figure supplement 1; Lein et al., 2007). Plotting gene expression with fixed scales across both genes reveals significantly higher expression levels of Cck than Foxb1 in the PMd cluster (column d1 in Figure 4d). Simultaneous visualization of both Cck and Foxb1 expression levels within the PMd cluster shows a biased expression pattern of Cck- and Foxb1-expressing cells towards opposing sides of the PMd cluster (column d2 in Figure 4d). There is a substantial number of single-positive cells for each of the two genes (column d3 in Figure 4d). There is indeed a group of PMd neurons co-expressing both markers, but very few of them show high expression levels of both genes.
In summary, although both Cck and Foxb1 are expressed throughout the entire PMd cluster, there are two well-defined, distinct, subpopulations of neurons expressing either Cck or Foxb1 but not both. We were hence able to introduce the presence of the PMdFoxb1 as a novel subdivision of the PMd with functional distinction from the PMdCck.
Hot plate experiments
Two studies had previously identified the parvalbumin+ subpopulation of the parvafox nucleus (parvafoxPvalb) to be involved in nociceptive behaviors (Roccaro-Waldmeyer et al., 2018; Siemian et al., 2019). We therefore extended the scope of our project to further investigate a possible reciprocal effect of the parvafoxFoxb1 on pain sensation.
In these experiments, we did not observe any significant difference between BL and Stim conditions in any of the three tested animal groups (Figure 5, ChR2-OnTarget_antPAG: Number of shakes before endpoint t(5) = –0.442, p=0.677 and latency until endpoint t(5) = 0.515, p=0.629; ChR-OffTarget: Number of shakes before endpoint t(10) = –1.076, p=0.307 and latency until endpoint t(10) = 1.716, p=0.117; ArchT3.0-OnTarget_antPAG: Number of shakes before endpoint t(5) = –0.344, p=0.753 and latency until endpoint z=1.150, p=0.875; all reported results are from Welch’s t tests, except the last result that was obtained by a Wilcoxon signed-rank test due to non-normality of the data).
These results reveal that the parvafoxFoxb1 does not excert a reciprocal effect on the parvafoxPvalb’s role in pain sensation. Interestingly though, the absence of altered latencies until endpoint behavior demonstrates that the immobility phenotype of ChR2-expressing Foxb1-Cre mice can be escaped, when thermal stimulus intensity reaches nociceptive threshold (Figure 5a).
Alterations of cardiovascular parameters upon optogenetic modulation of Foxb1 terminals in the dlPAG is observed in a small cohort of animals
The measurement of cardiovascular parameters in a group of 8 Foxb1-Cre mice revealed a sudden onset of bradycardia in three mice (106/21–10, 34/21–7, 34/21–10), immediately after starting the optogenetic activation (Figure 6a). The heart rate (HR) decreased abruptly by 30–60 beats/min but could reach 290 beats/min in some experiments (Figure 6b; 106/21–10). The light-induced cardiovascular responses had an on-kinetic of 1–10 s and at light offset, the HR returned to the baseline levels nearly instantaneously. Augmenting the light intensity of the stimulation increased the extent of HR deceleration. In these three mice, the glass fibers were located bilaterally over the rostral part of the dlPAG (Figure 7a and b). Glass fibers located over the Su3 or over the intermediate or caudal dlPAG or the PV2 nucleus did not trigger cardiovascular reactions. In the five Pvalb-Cre mice with similar injection of Chr2 in the parvafoxPvalb and glass fibers positioned bilaterally over the rostral dlPAG, no changes in cardiovascular parameters were measured.
In mouse 106-21/10 (Supplementary file 5), the HRV changed from 80ms during the baseline period to 120ms during the optogenetic stimulation the STDEV of cycle duration changed from 3.8 ms to 20.6ms [+445%].
The results of these cardiovascular experiments are to be considered as an observational report. Due to a limited number or animals, we cannot draw any statistically valid conclusions at a group level.
Discussion
Chemogenetic activation of Foxb1+ cell bodies in the PMd and parvafox nuclei of the hypothalamus, as well as optogenetic excitation of their axonal endings in the dlPAG both lead to a phenotypic presentation congruent with the reaction observed during innate defensive behavior: tachypnea, immobility, and bradycardia.
As both CNO and clozapine have been reported to activate DREADDs (with different temporal dynamics), we have designed the experimental outline of this study in a way that would allow us to also investigate the effect of both DREADD ligands against each other and against the effect of saline injections. We have shown that injection of both DREADD ligands reproducibly led to the same phenotype in DREADD-expressing animals but not in control animals. The activation of hM3Dq (activating DREADD receptors) expressed in the parvafoxFoxb1 led to changes in several respiratory parameters, particularly to an increase in respiratory frequency (measured in breaths per minute), consistent with an increased respiratory effort. We infer that this effect is at least partly mediated by the projections of the parvafoxFoxb1 neurons to the dlPAG. In accordance with this assumption, injections of excitatory aminoacids D,L-hopmocysteic acid (DLH) in the dlPAG raise respiratory activity by an increase in respiratory rate (Iigaya et al., 2010; Subramanian and Holstege, 2014; Dampney et al., 2013).
Immobility is the absence of movements associated with an increase in respiration and muscle tone (Fanselow, 1994; Kalin and Shelton, 1989) and is consistently correlated with heart rate deceleration (Vianna and Carrive, 2005; Walker and Carrive, 2003). The body immobility observed in the open field test perturbed the execution and therefore may have influenced the results of the hot-plate test.
In our optogenetic experiments, we cannot differentiate between the respective contributions of the Foxb1-terminals belonging to the parvafoxFoxb1 versus those belonging to the PMd to the behavioral (immobility) and autonomic (bradycardia) effects observed. Are these effects generated by a single common population of neurons, for example, all Foxb1+ neurons together (i.e. ‘command neurons’) or instead by the distinct sub-populations of neurons in the PMd or in the parvafoxFoxb1?
A subpopulation of neurons in both nuclei express Foxb1 and a distinct subpopulation consists of parvalbumin positive cells in the parvafox (Bilella et al., 2014) and Cck positive cells in the PMd (Wang et al., 2021). Ten percent of parvalbumin neurons in the parvafox express also Foxb1, whereas the situation for Cck in the PMd is yet unknown.
The PMd and the parvafoxFoxb1 originate from the same Foxb1+ cell lineage have differential connectivity pattern and as a result they may assume distinct functional tasks (Bilella et al., 2014; Wang et al., 2021). The PMd receives afferences from the infralimbic and prelimbic areas of the medial frontal cortex (Comoli et al., 2000), whereas the parvafoxFoxb1 neurons receives their inputs from the lateral and ventrolateral orbitofrontal cortex (Babalian et al., 2018). The neurons of the PMd have a bifurcated output, sending their axons rostrally to the anterior hypothalamic nucleus (AHN) and caudally to the dorsolateral sector of the PAG (Canteras and Swanson, 1992). Instead, the efferences of the parvafoxFoxb1 neurons have only minor rostral projections to the septal region and no rostral projections to other hypothalamic nuclei such the SF1+ neurons of the VHN (Kunwar et al., 2015). Instead, the parvafoxFoxb1 rather projects caudally to a wedge-shaped field straddling the lateral (lPAG) and dorsolateral (dlPAG) column (Bilella et al., 2016). A further columnar field of terminals is located in the Su3-region of the ventromedial PAG (Bilella et al., 2016).
The Foxb1+ neurons of the parvafox nucleus also use glutamate as a neurotransmitter (Bilella et al., 2016) and express the gene coding for the neuropeptide adenylate cyclase-activating polypeptide1 (Adcyap1; Girard et al., 2011), found to be involved among others in stress disorders (Ressler et al., 2011) and cardiorespiratory control (Shi et al., 2021; Barrett et al., 2019).
The PMd is predestined to play a pivotal functional role in the coordination of defensive behaviors, to both innate and conditioned threats (Wang et al., 2021). It consists of a predator activated ventrolateral part (vlPMd) which project to the dlPAG and a dorsomedial part (dmPMd) activated by encounter with aggressive conspecific which avoid projecting to the dlPAG. Ibotenic-acid lesions of the PMd eliminates escape and freezing responses (Canteras et al., 1997).
Recording Ca2+-transients in PMdCck neurons with fiber photometry, the group of Adhikari (Wang et al., 2021), revealed that these neurons are activated during escape but not during freezing and that their chemogenetic inhibition decreases escape speed from threats (Wang et al., 2021). The optogenetic activation of the PMdCck cells presumably leads to the activation of the glutamatergic neurons in the dlPAG, which then provoke the fleeing. In our experiments, however, activation of the larger group of Foxb1+ neurons in parvafoxFoxb1 and PMd led to immobility and not to escape. How can this paradoxical effect be explained? One potential explanation would be, that the parvafoxFoxb1 neurons, at least in the rostral dlPAG, innervate a different group of neurons, namely the abundant inhibitory GABA-cells (Barbaresi, 2005). It is to be assumed that these short-axon cells, with overlapping distribution in the dlPAG, would locally inhibit the activity of the glutamatergic neurons activated by the PMdCck neurons. Another subgroup of nitric oxide expressing neurons in the dlPAG (Onstott et al., 1993) also display a tonic inhibitory effect, via a potentiation of GABAergic synaptic inputs (Xing et al., 2008). The role of a fourth group of acetylcholinesterase positive neurons in the dlPAG (Illing, 1996) is unknown. The immobility of the mouse would then derive from the local inhibition of the excitatory elements of the dlPAG. Future work will be required to clarify this issue. A second potential explanation is rooted in the expression pattern of Foxb1 and Cck within the PMd. Our reanalysis of the scRNA seq data of the murine ventral-posterior hypothalamus has revealed, that there is a substantial number of Foxb1 and Cck single-positive cells in the PMd. Since our study and the one from Adhikari’s group both made use of single-gene Cre knock-in mouse lines, the optogenetic stimulation in both studies led to activation of both single-positive and double-positive neurons. We therefore propose that the observed effects in both studies may be exclusively mediated by single-positive neuronal populations within the PMd. Using an intersectional approach combining Cre- and Flp-recombination dependent expression of optogenetic tools, would allow to investigate whether the same behavioral effects can also be elicited by exclusively activating double-positive (Foxb1+/Cck+) neurons within the PMd.
Which other parts of the brain are recruited by the dlPAG neurons to initiate bradycardia and immobility? The dlPAG projects to the cuneiform nucleus (Meller and Dennis, 1991), which is involved in defensive locomotion (Jordan, 1998) through its connection with the lateral gigantocellular reticular nuclei (Tovote et al., 2016), implicated in high-speed locomotion (Caggiano et al., 2018). Additional targets of the cuneiform projections are the motor nucleus of the vagus and the nucleus tractus solitarius (Korte et al., 1992) and a few immunostained boutons were found around the nucleus ambiguus in the rostroventrolateral medullary nucleus (Figure 5K in Korte et al., 1992). These projections are probably involved in the bradycardic response to the optogenetic stimulation in the dlPAG. The projections of the cuneiform nucleus to the rostral ventrolateral medulla promote sympathetic vasomotor activity (Verberne, 1995). An alternative pathway for the expression of the behavioral and autonomic effects of dlPAG activation is through a link with the dorsomedial hypothalamus (DMH; Dampney et al., 2013).
Bradycardia starts immediately after the beginning of the optogenetic activation of the Foxb1+ terminals in the dlPAG. The increase in the HRV during the optogenetic stimulation indicates an intervention of the parasympathetic system on the sinoatrial node of the heart (Standish et al., 1995). The parasympathetic impact on the heart is mediated by acetylcholine neurotransmission (Wu et al., 2019) and has a very short latency of response, with peak effect at about 0.5 s and return to baseline within 1 s (Pumprla et al., 2002).
During ultrasonic vocalization tests in our laboratory (data not reported in this study), we observed that the Foxb1-Cre mice did not escape the ChR2-induced state of immobility, even when another genotype-matched and novel intruder mouse was placed into the ChR2-expressing mouse’s home cage. When approached by the intruder mouse, the mice were still able to move their heads under photo-stimulation, did sniff and interacted with the intruder mouse but did not actively follow the intruder. A mouse displaying complete immobility remained completely motionless even when the intruder was actively approaching it and even when it started digging into the cage bedding underneath the resident mouse.
In an interesting recent publication of Chen’s group (Liu et al., 2022), the projection of the parvafoxFoxb1 to the lateral PAG (lPAG) indeed has been found to drive social avoidance in mice. Inhibition of the parvafoxFoxb1 neurons by the GABAergic input from the lateral septum reversed the deficit in social novelty preference induced by chronic social defeat stress. The authors conclude that activation of the parvafoxFoxb1 leads to avoidance of general social situations.
Hence, our observations of a reduced social interaction phenotype agree with the results from Chen’s group. However, since activation of the parvafoxFoxb1 terminals in the dlPAG also leads to immobility of the mouse in the absence of any social stimulus, we argue that the outcome of the social avoidance test could have been influenced by the failing motor performance of the Foxb1-Cre mice and might not be of pure social nature. In light of these observations, it is further interesting to note, that the results from our hot plate experiments indicate that the immobile/hyporeactive phenotype of these mice can indeed be escaped, when an adequate stimulus is presented (e.g. thermal stimulus intensity reaching nociceptive levels).
The possibility of antidromic propagation of ChR2-mediated action potentials poses a major caveat to the interpretation of results obtained from experiments in which presynaptic terminals are optogenetically activated. This mechanism could lead to coincident activation and realease of neurotransmitters from collaterals originally not intended to be targeted by the intervention. The use of inhibitory optogenetic tools, however, circumvents such a mechanism of uncontrolled backpropagating action potentials (Rost et al., 2022). In our optogenetic experiments in which we have used the inhibitory opsin ArchT3.0, we detected statistically significant differences to the baseline conditions in multiple locomotive parameters. This observation of bidirectional modulation is an indication of at least some baseline activity of the Foxb1+ neurons in the hypothalamus in our experimental settings. Importantly, these results and observations additionally support the notion that the locomotor phenotype is due to the specific projection of Foxb1+ neurons onto the PAG. If the hypolocomotor phenotype in ChR2 expressing mice was driven by backpropagation-mediated activation of collaterals, the presynaptic ArchT3.0-mediated inhibition of the Foxb1+ projections to the PAG would be unlikely to cause an effect opposing the ChR2-mediated phenotype. In our ArchT3.0 experiments, we specifically employed a stimulation paradigm in which we had light illumination on for 10 s and then off for another 10 s. This was done to reduce pH-mediated side effects of ArchT3.0 axonal activation observed with sustained stimulation, like spontaneous release of neurotransmitters (Mahn et al., 2016). If such side effects had occured in our experiments, we would expect to see a phenotype that resembles the one observed in the ChR2-mediated activation. However, in our experiments, we saw no indication of such a behavior. Recently, two new optogenetic tools have been developed that allow efficient and targeted presynaptic silencing through GPCRs without the alteration in intracellular pH levels (Mahn et al., 2021). For future experiments, such tools might be more suitable for presyaptic silencing experiments.
In summary, our results demonstrate evidence for a role of the hypothalamic Foxb1+ neuronal population located in the parvafoxFoxb1, and the PMd in the expression of bradycardia and immobility, as well as in increasing respiratory rate. This effect is in accordance with distinctive coping patterns, the freezing-like behavior accompanied by bradycardia being associated with innate defensive behavior.
In view of the recent finding that PMd neurons are activated during escape (Wang et al., 2021), the immobility that we observe could be mediated by the terminal endings of the parvafoxFoxb1, synapsing on the inhibitory GABA-interneurons of the rostral third of the dlPAG or could be the result of the activation of the Foxb1 single-positive PMd neuronal population (i.e. Foxb1+/Cck-). Further investigations on the neuronal activity of these hypothalamic Foxb1+ populations during exposure to natural threats will allow for a better understanding of the differential roles of Foxb1+ and Cck+ PMd neurons. Such experiments could for example consist of fiber photometry calcium recordings or microendoscopic calcium imaging experiments during predator exposure in Foxb1 and/or Cck Cre-driver lines.
Materials and methods
Mice
A total of 42 mice of both sexes were used for the purpose of this study. Most animals were Foxb1tm1(cre)Gabo mice that express Cre-recombinase under the control of the promoter for Foxb1. For the cardiovascular experiments, five mice belonging to the Pvalb-Cre strain (129P2-Pvalbtm1(cre)Arbr/J) were used (Hippenmeyer et al., 2005). All animals were maintained at a constant temperature of 24 °C in state-of-the-art animal facilities with a 12hr-light/12hr-dark cycle and had ad libitum access to food and water.
All experiments were approved by the Swiss federal and cantonal committee for animal experimentation (2016_20E_FR and 2021–10-FR) and were conducted in accordance with the institutional guidelines of the University of Fribourg.
Intracerebral AAV injection
Request a detailed protocolMice in which we aimed at modulating neuronal activity were injected intracerebrally either with (i) Channelrhodopsin (AAV5-EF1α-DIO-hChR2(H134R)-eYFP), (ii) Archeorhodopsin (AAV2-EF1a-DIO-eArchT3.0-eYFP), (iii) activating DREADD (AAV2-hSyn-DIO-hM3D(Gq)-mCherry) or, (iv) inhibiting DREADD (AAV2-hSyn-DIO-hM4D(Gi)-mCherry). All optogenetic and chemogenetic agents injected into Foxb1-Cre and Pvalb-Cre mice were Cre-dependent and were bilaterally injected into the posterior portion of the parvafox nuclei of the LHA.
All intracerebral viral vector injections were conducted according to the following standard protocol:
The animal was weighed and was anesthetized with an intraperitoneal injection (i.p.) of a mixture of ketamine (40–60 mg.kg–1 of body weight) and xylazine (10–15 mg.kg–1 of body weight) diluted in physiological (0.9 %) saline. The fur covering the cranium was shaved and the mouse was head fixed into a stereotaxic apparatus (Kopf instruments, model 5000) equipped with a heating pad to maintain body temperature. Depth of anesthesia was assessed regularly throughout the entire duration of the surgery. If necessary, additional doses of the anesthetic agent were injected. Eye ointment was applied, and the eyes were protected from exposure to direct light. Once tail pinch and toe pinch reflexes were vanished, a sagittal skin incision above the midline of the cranium was performed and the cranial sutures were identified to locate bregma. Craniotomy was performed bilaterally above the site of injection with a dental steel bur. The viral vector was then aspirated into a 2.5 μl Hamilton syringe via a fine-bored 34-gauge needle (external diameter of 0.14 mm). The Hamilton syringe was mounted onto a manual microinjection unit and fixed to the stereotaxic frame. After identification of bregma, the needle tip was placed just above it and the anterior-posterior and medial-lateral coordinates for the injection sites were calculated relative to bregma (AP –1.3 mm, ML +/-1.3 mm). The needle was then placed on the brain surface above the target injection site and the dorsal-ventral coordinate was calculated relative to the brain surface (DV –5.5 mm). The needle was subsequently lowered into the brain until the desired depth was reached. 200 nl of the viral construct were injected into the parvafox nuclei bilaterally at a rate of 100 nl/min. Before retraction of the needle, the needle was left in place for 5 min to allow diffusion of the virus and to minimize backflush of the virus along the entry path of the needle. The same procedure was then repeated on the contralateral side. Once both viral injections were completed, the mouse was released from the stereotaxic frame and the skin incision was closed with one to two surgical stiches. The mouse was placed into a separate cage to recover from the surgery and was later put back into its home cage once consciousness was regained.
Fiber optic cannulas and implantation procedure
Request a detailed protocolFiber optic cannulas were custom made with ceramic ferrules (230 μm ID, 1.25 mm OD; Prizmatix Ltd., Israel) and optical fibers (200 µm OD, 0.66 NA; Prizmatix Ltd., Israel). The length of the protruding glass fiber end was set to 0.4 mm for the cannulas to be implanted above the dlPAG. The non-protruding end of the optical fiber was polished on a series of lapping sheets with decreasing grit size (5 μm, 3 μm, 1 μm, 0.3 μm). Each optical fiber implant was measured for its coupling efficiency before implantation. Implantation of the fiber optic cannulas was performed 2–3 weeks after virus injection.
The surgical procedure for implantation of the cannulas was as follows:
Initial preparation of the animal for access to the skull was performed as described above for viral vector injections. Once bregma was identified, small marks on the skull were made to identify the position of the bilateral bur holes for the cannulas and the three to four skull fixation screws. The bur holes were hand drilled with a dental steel bur and three to four fixation screws were screwed into the skull. The skull fixation screws were additionally locked to the skull by cyanoacrylate adhesive. A cannula was then mounted to the stereotaxic frame via a standard electrode holder (Kopf instruments) and was inserted into the brain at the desired coordinates. The bilateral cannulas implanted above the dlPAG were inserted at a 20° angle in the coronal plane (V-shaped arrangement of cannulas) to allow enough space for later connection to the patch cords. The coordinates for dlPAG cannulas insertion sites were AP –4.0 mm, ML +/-1.5 mm relative to bregma and insertion depth (at an angle of 20°) was –2.8 mm from the surface of the skull. Once the fibers were in place, self-curing acrylic (Palidur powder and liquid, Heraeus Kulzer GmbH, Hanau, Germany) was applied to the skull, the ceramic ferrule and the fixation screws. The mouse was then released from the stereotaxic apparatus and a 0.3 ml subcutaneous injection of physiological (0.9 %) saline into the neck scruff was made to support recovery from the surgery. Further, Carprofen (Rimadyl, 5 mg.kg–1 of body weight) was subcutaneously administered as an analgesic following surgery. After surgery, the mouse was placed in a separate cage under a warming lamp until consciousness was fully regained. Whenever possible, mice were put back into their home cages together with their cage mates. If fighting behavior was observed between male cage mates, the mice were kept in separate cages.
Optogenetic stimulation
Request a detailed protocolTime-resolved optogenetic gain-of-function manipulation of Foxb1+ axon terminals was performed bilaterally. A dual LED light source (MP-Nr: 34117; Prizmatix Ltd., Israel) with blue (peak λ=453 nm) and lime green (peak λ=536 nm) light emitting diodes (LEDs) were connected to a 1.5 mm optical fiber (NA 0.63; MP-Nr: 34131; Prizmatix Ltd.) which terminated into a rotary joint (MP-Nr: 34043, Prizmatix Ltd.), connecting to a dual fiber patch cord (2x500 μm, NA 0.63, MP-Nr: 34115; Prizmatix Ltd.). The two patch cord ends were attached to the cannulas on the head of the mouse by ceramic sleeves (ID = 1.25 mm, MP-Nr: 34071; Prizmatix Ltd.). ChR2 mice were stimulated with an intensity of 7–15 mW per fiber tip. The protocol for ChR2 activation consisted of bursts of 500ms duration with an intraburst frequency of 30 Hz, a pulse duration of 5ms and an interburst interval of 500ms and was based on the previously published firing properties of Foxb1+ neurons of the medial mammillary complex (Alonso and Llinás, 1992). ArchT3.0 stimulation protocol consisted of alternating 10 s windows with continuous LEDon and LEDoff, respectively (intensity of 7–15 mW). The optogenetic pulses were generated on a PulserPlus (MP-Nr.:34192, Prizmatix Ltd.) installed on a PC and were sent to the dual LED source and to an ACQ-7700 USB amplifier (PNM-P3P-7002SX, Data Sciences International [DSI], St. Paul, MN, USA).
A laser system of the type LRS-0473-PFO-00500–01 LabSpec (473 nm DPSS) from Laserglow technologies, North York, Ontario, Canada, was employed for augmenting the power output in the cardiovascular experiments (70–222 mW). The estimated power at the specimen was measured with a photodiode (Thorlabs).
Chemogenetic stimulation
Request a detailed protocolMice stereotaxically injected with activating or inhibiting DREADDs, as well as animals without expression of any DREADDs (i.e. DREADD_neg) were injected i.p. with CNO, clozapine or physiological saline in a given experimental block. CNO was administered at a dose of 1 mg∙kg–1 of bodyweight 30 min before the start of the experiment. Due to faster pharmacokinetics, clozapine was administered i.p. at a subthreshold dose of 0.1 mg∙kg–1 of bodyweight immediately before the start of the experiment. Attention was given to eventual backflush or unintentional subcutaneous administration. Dose calculations were adjusted individually to bodyweight before each injection.
Whole body barometric plethysmography (WBP)
Request a detailed protocolTo measure respiratory parameters, mice were placed in whole-body barometric plethysmography chambers (item Nr. 601-0001-011, DSI) that were connected to a Buxco Bias flow pump (item Nr. 601-2201-001, DSI, St. Paul, MN, USA) to avoid CO2 accumulation inside the chambers and to circulate ambient air through the chambers at a rate of 1 L∙min–1. Each chamber was equipped with a high sensitivity differential pressure transducer (Buxco TRD5700 Pressure, item Nr. 600-1114-002, DSI) as well as a temperature and humidity probe (item Nr. 600-2249-001, DSI). All sensors were connected to the ACQ-7700 amplifier for pre-processing of signals.
Acquisition and processing of respiratory signals as well as the TTL signal from the LED pulser was performed within the Ponemah Software (PNM-P3P-CFG, DSI). The respiratory flow signals were sampled at a rate of 500 Hz and were filtered with a 30 Hz low pass filter. Using the Unrestraint Plethysmography Analysis Module (URPM) within Ponemah, respiratory cycles were validated in the flow signal and several respiratory parameters were derived from them. For each animal, the attribute analysis settings were individually adjusted to obtain optimal validation marks of respiratory cycles. The derived respiratory parameters were then averaged into 5 s averages and exported as excel files for further statistical analysis and plotting in R and RStudio (RStudio, Inc, Boston, MA, USA).
On the experimental day, mice were brought into the experimental room and were allowed to habituate to the new environment for at least 45 min. Each animal had two habituation sessions of 90 min on two different days prior to the first baseline measurements. For these habituation sessions, DREADD mice were injected with 0.2 ml of physiological (0.9 %) saline i.p. 30 min before they were placed into the WBP chamber.
Each animal’s baseline (BL) was measured for 3x90 min on 3 consecutive days. BL condition for DREADD animals consisted in a 0.2 ml i.p. injection of physiological (0.9 %) saline 30 min before measurement. After BL recordings, DREADD animals were measured for 3x90 min under CNO stimulation on 3 consecutive days and for another 3x90 min under clozapine stimulation on another 3 consecutive days. Mice were given a break of at least 3 days between CNO and clozapine experimental blocks to allow for complete clearance of the substances and to reduce stress on the animal. WBP chambers were cleaned with soap and water after each recording to avoid olfactory stimulation of the next mouse. All WBP experiments for a given mouse were performed at the same time of the day to account for circadian variability.
Cardiovascular measurements with telemetry
Request a detailed protocolFor measuring cardiovascular functions, an implantable telemetry system from Datascience International (DSI) was used. The transmitter was implanted in 13 mice following manufacturer instructions and as previously described (Huetteman and Bogie, 2009; Pillai et al., 2018). Briefly, mice were anesthetized with the aid of isoflurane and implanted with a PA-C-10 transmitter (DSI). This small pressure sensing telemetry tool was implanted in the left carotid artery taking care to place the pressure-sensitive tip in the aortic arch. The radio transmitting device (RTD) was placed under the skin along the right flank of the mice. Mice were given analgesics for 3 days post-surgery. One week after telemetric sensor implantation, a test measurement was conducted to ensure that the catheter, the RTD and the receiver (RPC-1) were functioning well. The optogenetic experiment was performed after 2 more weeks.
Measurements were recorded with Dataquest ART (version 3.1) and RespiRate (DSI). Recordings were taken three times for 3 weeks post-surgery. Recordings were continuously taken for 30 s every 5 min during these sessions. Systolic, diastolic and mean blood pressure, heart rate and activity were analyzed. Pulse wave signals were used to measure the heart rate variability (HRV).
Open field test
Request a detailed protocolAs in all other experiments reported in this paper, mice were transported into the experimental room and were allowed to acclimatize to the new environment for at least 45 min. The experimental arena for the open field test consisted in a 40x40 cm cage with transparent plexiglas walls and a gray colored metal floor. After each recording of a mouse, the arena was thoroughly cleaned with 70% ethanol. Each mouse was recorded four times (2 x saline/LEDoff and 2 x clozapine/LEDon) spread across 2 days with at least 2 days between experimental days. In DREADD experiments, the saline condition was performed in the morning and the clozapine session in the afternoon. To account for potential circadian bias, optogenetic morning and afternoon sessions for each mouse were alternated, so that the two recordings for LEDoff and LEDon, respectively, were recorded once in the morning session and once in the afternoon session. Recordings consisted of a 5-min habituation period inside the experimental arena, uninterruptedly followed by a 5-min recording period. Mice for DREADD experiments were injected with clozapine (i.p. 0.1 mg∙kg–1 of bodyweight) 30 min before the recording time window. Optogenetic mice remained in their home cages until the start of the 5-min habituation session, shortly before which they were attached to the patch cords. In all optogenetic animals (including controls), a stop watch signal marked the end of the 5 min of habituation. The experimenter then initiated the LEDon period accompanied by another brief auditory signal.
Pose estimation of open field data
Request a detailed protocolFor body part tracking, DeepLabCut (version 2.2.2) (Mathis et al., 2018; Nath et al., 2019) was used. Specifically, 500 frames taken from 25 videos were labeled and 95% was used for training. A ResNet-50-based neural network with default parameters was used for 600,000 training iterations. We validated with 1 shuffle, and found the test error was: 2.03 pixels, train: 2.39 pixels (image size was 1280 by 720 pixels). We then used a p-cutoff of 0.95 to condition the X,Y coordinates for future analysis. This network was then used to analyze videos from similar experimental settings.
Hot plate test
Request a detailed protocolTo assess thermal nociceptive perception in optogenetic Foxb1-Cre mice, they were first connected to the patch cords and placed onto an insulation layer of cork and several layers of paper towels on the hot plate for 3 min to acclimatize to the new environment. The hot plate (Analgesia meter for rodents, IITC life sciences Inc, Woodland Hills, CA, USA) was maintained at 51+/-0.1 °C. After 3 min had passed, the insulation layer was removed and the latencies until hindlimb shaking, hindlimb licking and jumping were recorded as baseline. In case the mouse did not display any of the two endpoint behaviors (i.e. hind paw licking or jumping), the recording was terminated 50 s after the onset of the thermal stimulus to prevent tissue damage. After termination of the baseline recording, the insulation layer was again placed between the mouse and the hot plate and the optogenetic stimulation was initiated (see ‘Optogenetic stimulation’ for detailed parameters). After 3 min of optogenetic stimulation, the insulation was once again removed and the mouse was placed back onto the hot plate, while the stimulation continued. Just like in the baseline condition, the same latencies were recorded, or the experiment was terminated after maximally 50 s. As observed in pilot experiments, the ChR2-injected mice displayed reduced locomotor activity and seemed to be limited in their ability to lick their hind paw. We therefore also terminated the recording, when the mouse displayed an obvious attempt to lick its hind paw and the latency until the display of such an event was recorded. The same procedure was repeated one more time on another day to record two baseline and two LED condition recordings for each optogenetic animal.
Histology
Request a detailed protocolBefore mice were perfused, they went through the same routine as in a regular experimental condition (LED stimulation or CNO injection) for later detection of c-Fos immunofluorescence. Two hours after the start of the activation/inhibition experiments, the mice were deeply anesthetized with the same anesthetic agent used for surgical procedures (see above). Once the pain reflexes vanished, the thorax was fenestrated, and the mouse was transcardially perfused with physiological (0.9 %) saline for 3 mins and subsequently with 4% paraformaldehyde (PFA) in PBS 0.1 M (pH 7.4) for 5 min. Decapitation was performed and the head was placed in 4% PFA until extraction of the glass fiber implants. After cannulas extraction from the skull, the cannulas were stored for post-extraction coupling efficiency measurement. The brains were then harvested and further immersed in TBS 0.1 M+18% sucrose +0.02% Na-azide overnight for cryoprotection.
Brain tissue was cryo-sectioned into 40-μm-thick sections on a sliding microtome (SM 2010R, Leica) connected to a freezing unit (Microm KS 34, Thermo Fisher Scientific) and every sixth section was selected for incubation with antibodies.
The histological processing of free-floating sections was performed as follows:
Sections were washed for 3x5 min in TBS 0.1 M. The sections were then incubated for 2 days at 4 °C with the primary antibodies at their corresponding dilutions (see Supplementary file 6) in TBS 0.1 M+0.1% Triton X-100 and 10% bovine serum (BS). After incubation with the primary antibodies, sections were washed for 3x5 min in TBS 0.1 M and further incubated for 2 hr at room temperature with a biotinylated antibody diluted 1:200 in TBS 0.1 M and 10% BS. Subsequently, sections were washed 1x5 min in TBS 0.1 M and 2x5 min in Tris pH 8.2 before they were incubated with Cy2-, Cy3- and Cy5-conjucated secondary antibodies or streptavidin, which were diluted 1:200 in Tris pH 8.2 for 2 hr at room temperature. After incubation with Cy-conjugated secondary antibodies and Cy-conjugated streptavidin, the sections were washed 1x5 min in Tris pH 8.2 and 2x5 min in TBS 0.1 M and then incubated with DAPI 1:5000 in TBS 0.1 M for 5 min at room temperature.
Sections were mounted onto Superfrost +glass slides (Thermo Scientific) and were left to dry for 2 hr at 37 °C. Slides were then quickly washed in dH2O before standard cover slips were mounted to the slides with Hydromount mounting medium (National Diagnostics, Atlanta, GA, USA).
Genotyping
Request a detailed protocolGenotyping was performed before animal selection and again after perfusion to exclude any mix-up during the testing period. For the first genotyping, tissue samples from the toe clipping procedure were used. For the second genotyping, tissue samples were taken from the tail after the animal was deeply anesthetized and before perfusion with 4% PFA.
The primer sequences used for Foxb1-Cre genotyping were as follows:
-EGFP-f: 5’-CTC GGC ATG GAC GAG CTG TAC AAG-3’
-GAB20: 5’-CAC TGG GAT GGC GGG CAA CGT CTG-3’
-GAB22: 5’-CAT CGC TAG GGA GTA CAA GAT GCC-3’
Reanalysis of single-cell RNA sequencing data set
Request a detailed protocolThe raw data used for the reanalysis of single-cell RNA sequencing reads from mouse (postnatal day 30–34) ventral-posterior hypothalami was kindly made available to the public by the original authors (Mickelsen et al., 2020) through the gene expression omnibus under the accession number ‘GSE146692’.
The entire analysis workflow was executed as follows using the Seurat package V4.1.1 in R (Hao et al., 2021):
Datasets of two male and two female animals were imported into R and initialized as four separate Seurat objects. Subsequently, the four separate objects were merged into a common Seurat object and the percentage of mitochondrial transcripts as well as the number of hemoglobin transcripts per cell were calculated. Quality control was performed by analyzing number of features, number of counts, percentage of mitochondrial RNA, and number of hemoglobin gene transcripts. Data meeting the following criteria were kept for downstream analysis: nFeatures >200 & nFeatures <7500 & percent.mtRNA <15 & nHemoglobin_RNA <50. The filtered dataset was then normalized within each original identity (male1, male2, female1, female2) using the SCTransform function and 3000 integration features were selected. To correct for batch effects, integration anchors were detected based on the selected integration features prior to integration of the four datasets. Next, uniform manifold approximation and projection (UMAP) was performed for dimensional reduction before clusters were identified and plotted for inspection. We then plotted the expression levels of a set of candidate genes (features) to identify the cluster representing the PMd and to differentiate it from the lateral and medial premammillary nuclei (LM and MM). To confirm our cluster identification, we extracted markers for the identified PMd cluster as well as differential markers for the PMd cluster vs. MM and LM clusters. Qualitative comparison of these markers with in situ hybridization data from the Allen Mouse Brain Atlas confirmed PMd cluster identity. Within the PMd cluster, we then extracted and plotted the level of co-expression of Cck and Foxb1 within each cell of the cluster.
Data analysis and statistics
Request a detailed protocolStatistical analysis and data handling was performed with custom written codes in R/RStudio and with Python in Jupyter Notebooks.
WBP data was first plotted as line plots and violin plots with the ggplot2 package within R. To check whether data matched the test assumptions, data was first assessed with a Shapiro-Wilks test and with visual inspection of Q-Q plots.
Before the analysis of the WBP data, outliers above and below 1.5 x the interquartile range were removed and the data was averaged across each condition (saline, clozapine, CNO). WBP experiments were analyzed by a 3x3 mixed-design ANOVA with Huynh-Feldt Sphericity correction, followed by post-hoc two-tailed paired Welch’s t-Tests. Hedge’s g with a correction for paired data (Gibbons et al., 1993) was calculated for all statistically significant pairwise comparisons.
The hot plate and open field experiments were analyzed with two-tailed paired Welch’s t-Tests if assumptions of normality were met. Otherwise, Wilcoxon signed-rank tests were used. Statistical analysis was limited to the groups that contained at least three subjects. Hedge’s g with a correction for paired data (Gibbons et al., 1993) was calculated for all statistically significant pairwise comparisons. Pose estimation data from open field experiments were further analyzed and plotted in R. A set of functions for the analysis of DeepLabCut labeled data (Sturman et al., 2020) was used and modified where necessary according to the analysis requirements. The open field arena was partitioned into arena, periphery, center, and corner zones according to the default parameters. Animal movement and zone visits were calculated with a movement cut-off of 1 and an integration period of 13 (with a video frame rate of 25 frames/second). For optogenetic experiments, a 3-min bin immediately before (baseline) and 3 min immediately after (stim) the start of the optogenetic stimulation were used for statistical comparison and plotting. For chemogenetic experiments, 5 min bins (following a 5-min habituation) were used from each experiment (2 x BL, 2 x Clo for each animal) for statistical comparison and plotting. The telemetrically recorded cardiovascular and movement-related data were analyzed in Ponemah. To be able to identify the LEDon and LEDoff periods in these data, two video cameras were used in parallel: camera 1 was recording the animal (photostimulation light visible), and camera 2 was recording a computer screen displaying the telemetric recording in real-time. To quantify the effect of the optogenetic activation (respectively of the change from LEDoff of LEDon in control animals), data of the last 60 s of the LEDoff condition were used as baseline data and statistically compared to the data of the first 60 s of the LEDon condition. Multivariate analysis of variance (MANOVA) was used to test for the presence of overall differences between the LEDoff vs. LEDon conditions. If significant, pair-wise ordinary least squares (OLS) was used for post-hoc analysis. To correct for multiple testing (four dependent variables), p=0.5/4=0.125 was set as threshold level for post-hoc significance (Bonferroni correction).
To assess heart rate variability (HRV), the heart rate data were analyzed in full, that is 1-epoch resolution. The duration of each individual heart cycle was derived by calculating the difference between the timestamp of each pair of two subsequent heart cycles. The mean, median and SD of all heart cycle durations during the last baseline minute (LEDoff), respectively during the first minute of the LEDon condition, were calculated. The SD of all heart cycle durations during that minute was used as an indicator for HRV.
Data availability
Raw data generated in this study was deposited on Dryad at: https://doi.org//10.5061/dryad.rv15dv4fn.Custom written code generated for this publication are available on GitHub in the following repository: https://github.com/rbcola/Parvafox-eLife-2024, (copy archived at Cola, 2024).
-
Dryad Digital RepositoryChemo- and optogenetic activation of hypothalamic Foxb1-expressing neurons and their terminal endings in the rostral-dorsolateral PAG leads to tachypnea, bradycardia, and immobility.https://doi.org/10.5061/dryad.rv15dv4fn
-
NCBI Gene Expression OmnibusID GSE146692. Single cell RNA sequencing to classify molecularly distinct neuronal and non-neuronal cell types in the mouse ventral posterior hypothalamus.
References
-
Topography and synaptology of mamillary body projections to the mesencephalon and pons in the ratThe Journal of Comparative Neurology 301:214–231.https://doi.org/10.1002/cne.903010206
-
Electrophysiology of the mammillary complex in vitro. II. medial mammillary neuronsJournal of Neurophysiology 68:1321–1331.https://doi.org/10.1152/jn.1992.68.4.1321
-
Expression of Foxb1 reveals two strategies for the formation of nuclei in the developing ventral diencephalonDevelopmental Neuroscience 22:197–206.https://doi.org/10.1159/000017442
-
Impaired neonatal cardiorespiratory responses to hypoxia in mice lacking PAC1 or VPAC2 receptorsAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 316:R594–R606.https://doi.org/10.1152/ajpregu.00250.2018
-
The Foxb1-expressing neurons of the ventrolateral hypothalamic parvafox nucleus project to defensive circuitsThe Journal of Comparative Neurology 524:2955–2981.https://doi.org/10.1002/cne.24057
-
Projections of the ventral premammillary nucleusThe Journal of Comparative Neurology 324:195–212.https://doi.org/10.1002/cne.903240205
-
SoftwareParvafox-eLife-2024, version swh:1:rev:57072f2dcc7375e330f5f8fa99b5c003fe484cb5Software Heritage.
-
An autoradiographic study of the descending connections of the mammillary nuclei of the ratThe Journal of Comparative Neurology 176:631–644.https://doi.org/10.1002/cne.901760411
-
Periaqueductal gray neuronal activities underlie different aspects of defensive behaviorsThe Journal of Neuroscience 36:7580–7588.https://doi.org/10.1523/JNEUROSCI.4425-15.2016
-
Neural organization of the defensive behavior system responsible for fearPsychonomic Bulletin & Review 1:429–438.https://doi.org/10.3758/BF03210947
-
Estimation of effect size from a series of experiments involving paired comparisonsJournal of Educational Statistics 18:271–279.https://doi.org/10.3102/10769986018003271
-
Gene expression analysis in the parvalbumin-immunoreactive PV1 nucleus of the mouse lateral hypothalamusThe European Journal of Neuroscience 34:1934–1943.https://doi.org/10.1111/j.1460-9568.2011.07918.x
-
Topographical specificity of regulation of respiratory and renal sympathetic activity by the midbrain dorsolateral periaqueductal grayAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 299:R853–R861.https://doi.org/10.1152/ajpregu.00249.2010
-
Initiation of locomotion in mammalsAnnals of the New York Academy of Sciences 860:83–93.https://doi.org/10.1111/j.1749-6632.1998.tb09040.x
-
Mesencephalic cuneiform nucleus and its ascending and descending projections serve stress-related cardiovascular responses in the ratJournal of the Autonomic Nervous System 41:157–176.https://doi.org/10.1016/0165-1838(92)90137-6
-
Biophysical constraints of optogenetic inhibition at presynaptic terminalsNature Neuroscience 19:554–556.https://doi.org/10.1038/nn.4266
-
DeepLabCut: markerless pose estimation of user-defined body parts with deep learningNature Neuroscience 21:1281–1289.https://doi.org/10.1038/s41593-018-0209-y
-
Space, time, and fear: survival computations along defensive circuitsTrends in Cognitive Sciences 24:228–241.https://doi.org/10.1016/j.tics.2019.12.016
-
Kidney Mass Reduction Leads to l-Arginine Metabolism-Dependent Blood Pressure Increase in MiceJournal of the American Heart Association 7:5.https://doi.org/10.1161/JAHA.117.008025
-
Functional assessment of heart rate variability: physiological basis and practical applicationsInternational Journal of Cardiology 84:1–14.https://doi.org/10.1016/s0167-5273(02)00057-8
-
Optogenetics at the presynapseNature Neuroscience 25:984–998.https://doi.org/10.1038/s41593-022-01113-6
-
Functional specializations within the tectum defense systems of the ratNeuroscience and Biobehavioral Reviews 29:1279–1298.https://doi.org/10.1016/j.neubiorev.2005.05.006
-
Central neuronal circuit innervating the rat heart defined by transneuronal transport of pseudorabies virusThe Journal of Neuroscience 15:1998–2012.https://doi.org/10.1523/JNEUROSCI.15-03-01998.1995
-
Neuronal circuits for fear and anxietyNature Reviews Neuroscience 16:317–331.https://doi.org/10.1038/nrn3945
-
Cuneiform nucleus stimulation produces activation of medullary sympathoexcitatory neurons in ratsThe American Journal of Physiology 268:R752–R758.https://doi.org/10.1152/ajpregu.1995.268.3.R752
-
PHA-L analysis of projections from the supramammillary nucleus in the ratThe Journal of Comparative Neurology 326:595–622.https://doi.org/10.1002/cne.903260408
-
Changes in cutaneous and body temperature during and after conditioned fear to context in the ratThe European Journal of Neuroscience 21:2505–2512.https://doi.org/10.1111/j.1460-9568.2005.04073.x
-
How do amusement, anger and fear influence heart rate and heart rate variability?Frontiers in Neuroscience 13:1131.https://doi.org/10.3389/fnins.2019.01131
Article and author information
Author details
Funding
Swiss National Science Foundation (grant 31003A_160325)
- Marco R Celio
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Christiane Marti and Laurence Clément for their technical support. The first part of the plethysmography experiments was conducted by SN in the laboratory of Prof. Julian Paton, Physiology Dept., University of Bristol (UK), now at the University of Auckland, New Zealand. We also thank Prof. Dr. Matthias Hänggi, Intensive Care Unit, University Hospital Bern for lending us the telemetric receiver equipment. We are further grateful for the financial support provided by the Swiss National Science Foundation (SNSF grant 31003 A_160325).
Ethics
All experiments were approved by the Swiss federal and cantonal committee for animal experimentation (2016_20E_FR and 2021-10-FR) and were conducted in accordance with the institutional guidelines of the University of Fribourg.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.86737. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Cola 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,079
- views
-
- 51
- downloads
-
- 0
- 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
The relation between neural activity and behaviorally relevant variables is at the heart of neuroscience research. When strong, this relation is termed a neural representation. There is increasing evidence, however, for partial dissociations between activity in an area and relevant external variables. While many explanations have been proposed, a theoretical framework for the relationship between external and internal variables is lacking. Here, we utilize recurrent neural networks (RNNs) to explore the question of when and how neural dynamics and the network’s output are related from a geometrical point of view. We find that training RNNs can lead to two dynamical regimes: dynamics can either be aligned with the directions that generate output variables, or oblique to them. We show that the choice of readout weight magnitude before training can serve as a control knob between the regimes, similar to recent findings in feedforward networks. These regimes are functionally distinct. Oblique networks are more heterogeneous and suppress noise in their output directions. They are furthermore more robust to perturbations along the output directions. Crucially, the oblique regime is specific to recurrent (but not feedforward) networks, arising from dynamical stability considerations. Finally, we show that tendencies toward the aligned or the oblique regime can be dissociated in neural recordings. Altogether, our results open a new perspective for interpreting neural activity by relating network dynamics and their output.
-
- Neuroscience
The circadian clock, an internal time-keeping system orchestrates 24 hr rhythms in physiology and behavior by regulating rhythmic transcription in cells. Astrocytes, the most abundant glial cells, play crucial roles in CNS functions, but the impact of the circadian clock on astrocyte functions remains largely unexplored. In this study, we identified 412 circadian rhythmic transcripts in cultured mouse cortical astrocytes through RNA sequencing. Gene Ontology analysis indicated that genes involved in Ca2+ homeostasis are under circadian control. Notably, Herpud1 (Herp) exhibited robust circadian rhythmicity at both mRNA and protein levels, a rhythm disrupted in astrocytes lacking the circadian transcription factor, BMAL1. HERP regulated endoplasmic reticulum (ER) Ca2+ release by modulating the degradation of inositol 1,4,5-trisphosphate receptors (ITPRs). ATP-stimulated ER Ca2+ release varied with the circadian phase, being more pronounced at subjective night phase, likely due to the rhythmic expression of ITPR2. Correspondingly, ATP-stimulated cytosolic Ca2+ increases were heightened at the subjective night phase. This rhythmic ER Ca2+ response led to circadian phase-dependent variations in the phosphorylation of Connexin 43 (Ser368) and gap junctional communication. Given the role of gap junction channel (GJC) in propagating Ca2+ signals, we suggest that this circadian regulation of ER Ca2+ responses could affect astrocytic modulation of synaptic activity according to the time of day. Overall, our study enhances the understanding of how the circadian clock influences astrocyte function in the CNS, shedding light on their potential role in daily variations of brain activity and health.