# NMDA-receptor-dependent plasticity in the bed nucleus of the stria terminalis triggers long-term anxiolysis

Christelle Glangetas, Léma Massi, Giulia R. Fois, Marion Jalabert, Delphine Girard, Marco Diana, Keisuke Yonehara, Botond Roska, Chun Xu, Andreas Lüthi, et al.

► To cite this version:

Christelle Glangetas, Léma Massi, Giulia R. Fois, Marion Jalabert, Delphine Girard, et al., NMDA-receptor-dependent plasticity in the bed nucleus of the stria terminalis triggers long-term anxiolysis. Nature Communications, 2017, 8, pp.14456. 10.1038/ncomms14456. hal-01485729

HAL Id: hal-01485729

https://hal.science/hal-01485729

Submitted on 8 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

Received 2 Feb 2016 | Accepted 3 Jan 2017 | Published 20 Feb 2017

DOI: 10.1038/ncomms14456

OPEN

# NMDA-receptor-dependent plasticity in the bed nucleus of the stria terminalis triggers long-term anxiolysis

Christelle Glangetas<fnref n="1" />,²,*, Léma Massi³,*, Giulia R. Fois⁴,*, Marion Jalabert⁵,⁶, Delphine Girard<fnref n="1" />,²,⁴, Marco Diana⁷, Keisuke Yonehara³, Botond Roska³, Chun Xu³, Andreas Lüthi³, Stéphanie Caille⁸,⁹ & François Georges<fnref n="1" />,²,⁴

Anxiety is controlled by multiple neuronal circuits that share robust and reciprocal connections with the bed nucleus of the stria terminalis (BNST), a key structure controlling negative emotional states. However, it remains unknown how the BNST integrates diverse inputs to modulate anxiety. In this study, we evaluated the contribution of infralimbic cortex (ILCx) and ventral subiculum/CA1 (vSUB/CA1) inputs in regulating BNST activity at the single-cell level. Using trans-synaptic tracing from single-electroporated neurons and *in vivo* recordings, we show that vSUB/CA1 stimulation promotes opposite forms of *in vivo* plasticity at the single-cell level in the anteromedial part of the BNST (amBNST). We find that an NMDA-receptor-dependent homosynaptic long-term potentiation is instrumental for anxiolysis. These findings suggest that the vSUB/CA1-driven LTP in the amBNST is involved in eliciting an appropriate response to anxiogenic context and dysfunction of this compensatory mechanism may underlie pathologic anxiety states.

<footnote n="1">¹Université de Bordeaux, Interdisciplinary Institute for Neuroscience, UMR 5297, F-33076 Bordeaux, France. ²Centre National de la Recherche Scientifique, Interdisciplinary Institute for Neuroscience, UMR 5297, F-33076 Bordeaux, France. ³Friedrich Miescher Institute for Biomedical Research, Maulbeerstrasse 66, 4058 Basel, Switzerland. ⁴Centre National de la Recherche Scientifique, Neurodegeneratives Diseases Institute, UMR 5293, F-33076 Bordeaux, France. ⁵Université de la Méditerranée UMR S901, F-13009 Aix-Marseille 2, France. ⁶INMED, F-13009 Marseille, France. ⁷‘G. Minardi’ Cognitive Neuroscience Laboratory, Department of Chemistry and Pharmacy, University of Sassari, 07100 Sassari, Italy. ⁸Université de Bordeaux, Institut de Neurosciences Cognitives et Intégrative d'Aquitaine, BP31, F-33076 Bordeaux, France. ⁹Centre National de la Recherche Scientifique, UMR 5287-Institut de Neurosciences Cognitives et Intégrative d'Aquitaine, F-33076 Bordeaux, France. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to F.G. (email: francois.georges@u-bordeaux.fr).</footnote>
Anxiety is a physiological negative emotion that triggers a state of alert to possible threat and promotes survival. The bed nucleus of the stria terminalis (BNST) belongs to a neuronal network of interconnected limbic regions and exerts a pivotal role in the expression of anxiety both in humans¹,² and in rodent models³,⁴. Previous studies using electrical⁵, lesioning⁶,⁷, pharmacological⁸ or optogenetic³,⁴ manipulations targeting directly specific BNST nuclei or cell types have been informative but have not been able to determine precisely how different BNST inputs influence anxiety. It has been shown that the ventral subiculum/ CA1 (vSUB/CA1) and infralimbic cortex (ILCx) massively projects to the BNST⁹. However, the integrative properties of BNST neurons at the single-cell level have never been explored. One hypothesis is that cortical and hippocampal information is processed at the single-cell level in BNST to trigger anxiolysis.

One of the most likely neural mechanisms underlying persistent anxiety is long-lasting plasticity in the neuronal network¹⁰. Probably owing to their anatomically convergent and segregated excitatory inputs from the infralimbic Cortex (ILCx) and the vSUB/ CA1 and their high content in stress-related neuromodulators, BNST circuits display morphological or synaptic plastic adaptations in response to stress and anxiety¹¹⁻¹³. Ex vivo studies in slice have correlated plasticity at glutamate synapses within the BNST with alterations in anxiety levels¹¹, but it remains unclear how the BNST integrates diverse inputs to modulate anxiety-related behaviors. Synaptic homeostasis seems to be a crucial process to compensate a long-lasting enhancement in signal transmission and maintain the stability of neuronal activity¹⁴,¹⁵. Recent computational modeling of synaptic plasticity have shown that the homeostatic processes that control the network stability are supported by the interaction of homosynaptic plasticity with heterosynaptic plasticity¹⁶. To unravel the integrative properties of BNST neurons at the single-cell level and the mechanism of their plastic changes in response to specific input stimulation, we combined conventional tracing, trans-synaptic tracing from single-electroporated neurons, in vivo single-cell recordings, pharmacological and behavioral techniques. Here, we demonstrate: (1) that both vSUB/CA1 and ILCx converge on the same amBNST neurons; (2) that in vivo stimulation of the vSUB/CA1 promotes an N-methyl-D-aspartic acid or N-methyl-D- aspartate (NMDA)-dependent long term potentiation (LTP) at the vSUB/CA1-amBNST synapses (LTP<sub>vSUB/CA1</sub>), associated with an NMDA-independent long-term depression (LTD) at the ILCx- amBNST synapses (LTD<sub>ILCx</sub>); and (3) that the induction of in vivo NMDA-R-dependent plasticity in the amBNST triggers long-term changes of anxiety state.

## Results

### BNST neurons are connected to both vSUB/CA1 and ILCx inputs.

It has been shown that the BNST could integrate information from the vSUB/CA1 or the ILCx (refs 9,12,17–19). However, it was unknown whether these inputs were integrated separately by different BNST neurons or if the same BNST neuron integrates both inputs. To address this question, we first injected a retrograde tracer, the cholera-toxin-B subunit (CTb), into the BNST, and confirmed that the BNST receives strong innervations from the vSUB/CA1 and the ILCx (Supplementary Fig. 1a–c). We then injected two different anterograde tracers, Phaseolus-vulgaris-leucoagglutinin (PHAL) and biotinylated-dextran-amine (BDA) in the vSUB/CA1 and in the ILCx, respectively (Fig. 1a,b). We used the Fox3 protein as a neuronal marker in the BNST. Using confocal imaging, we showed strong convergences of ILCx and vSUB/CA1 terminal fibers on the same Fox3-positive amBNST neuron (Fig. 1c). By using *in vivo* electrophysiology in anesthetized rats, we investigated whether amBNST neurons were controlled by both vSUB/CA1 and ILCx projections at the single-

<figure><img src="image_3.png" /><figcaption>**Figure 1 | Characterization of vSUB/CA1-ILCx-amBNST pathways.**</figcaption></figure>

(a,b) Injection sites: anterograde tracers (PHAL, BDA) were respectively injected in the vSUB/CA1 (a, magenta) and the ILCx (b, green). Scale, 0.5 mm.

(c) Immunofluorescence confocal images showing strong convergences (white arrows) of vSUB/CA1 (magenta) and ILCx (green) terminals, on amBNST neuron (Fox3 labelling, blue). Scale bar, 10 µm. (d) Experimental protocol. (e) Percentage of amBNST neurons responding to the stimulation of the ILCx, vSUB/CA1 or both. (f) Typical PSTHs (+ raster) showing respective excitatory responses of a same amBNST neuron to ILCx (top) and vSUB/CA1 (bottom) stimulation. Stimulus: t0 (gray lines). PSTH with 500 ms scale and Bin width, 5 ms (left) and PSTH of the same neuron with 50 ms scale and bin width, 1 ms on the right. Insets: representative traces.

cell level. When we tested the evoked responses of all recorded amBNST neurons to the stimulation of the ILCx and the vSUB/ CA1, we observed that 85.2% of the recorded neurons responded
to the stimulation of one of these inputs (Fig. 1d). We quantified that 70% of these responding neurons responded to both ILCx and vSUB/CA1 stimulation, (Fig. 1e,f), whereas only 17% and 13% responded to either ILCx or vSUB/CA1 stimulation, respectively (Fig. 1e). The basal firing rate of the responding amBNST neurons was 0.7 ± 0.2 Hz. The onset of the evoked responses for the ILCx-amBNST pathway was 6.94 ± 0.59 ms and 9.75 ± 0.77 ms for the vSUB/CA1-amBNST pathway (*P* < 0.01). Durations of evoked excitatory responses did not differ between these two pathways (ILCx-amBNST: 12.06 ± 1.06 ms; vSUB/CA1-amBNST: 12.78 ± 1.16 ms). The intensity of stimulation used to trigger evoked responses in the amBNST was in the same range between the two inputs (0.2–1 mA). To further dissect connectivity of amBNST inputs onto single-amBNST neurons, we performed *in vivo* single-cell electroporation of amBNST neuron followed by retrograde monosynaptic tracing with pseudotyped rabies virus strategy (Fig. 2a,b)²⁰. Five days after electroporation, retrogradely infected neurons were detected in the ILCx and in the vSUB/CA1 (Fig. 2b). Unfortunately, starter cells were not detected in the amBNST, most probably because rabies virus infection will eventually induce cytotoxicity in infected neurons²¹. Together, our findings confirm that a large population of

amBNST neurons was synaptically controlled, at the single-cell level, by both vSUB/CA1 and ILCx inputs.

**HFS<sub>vSUB/CA1</sub> promotes *in vivo* input-specific LTD and LTP.** To unravel the integrative properties of BNST neurons at the single-cell level, we combined *in vivo* single-cell recordings and pharmacological approaches. Given that vSUB/CA1 neurons fire at hundreds of Hertz in basal condition (Supplementary Fig. 2a–c), we first assessed the impact of a high-frequency stimulation (HFS; Supplementary Fig. 2d,e) of the vSUB/CA1 (HFS<sub>vSUB/CA1</sub>), on the ability of the vSUB/CA1-amBNST and ILCx-amBNST synapses to undergo plasticity (Fig. 3a–d). This HFS<sub>vSUB/CA1</sub>, originally described by Abraham *et al.*²², induced an extremely long-lasting (> 8 days), robust and stable LTP in the projection from vSUB/CA1 to medial prefrontal cortex (mPFC)²³. Here, we established that HFS<sub>vSUB/CA1</sub> triggers input-specific neuroplastic changes in the BNST (Fig. 3). After 30–40 min, HFS<sub>vSUB/CA1</sub>, amBNST neurons exhibited a 42.1% increase in the vSUB/CA1-evoked spike probability (LTP<sub>vSUB/CA1</sub>) together with a 47.4% decrease in the ILCx-evoked spike probability (LTD<sub>ILCx</sub>; $F_{(1,15)} = 35.97$, $P < 0.0001$ for the HFS<sub>vSUB/CA1</sub> protocol × input-stimulated interaction; $F_{(1,15)} = 80.60$, $P < 0.0001$ for the effect of the input stimulated, Fig. 3c–e). No change of the basal firing rate of amBNST neurons was observed before and after the HFS<sub>vSUB/CA1</sub> (2.74 ± 1.03 Hz versus 1.91 ± 0.74 Hz; $P = 0.401$). We next delivered HFS<sub>vSUB/CA1</sub> in the presence of intra-BNST infusion of the NMDA-R antagonist AP5 (100 µM). Under these conditions, HFS<sub>vSUB/CA1</sub> failed to elicit LTP<sub>vSUB/CA1</sub> (41.79 ± 13.16% of baseline, $n = 10$, $P < 0.001$; $F_{(1,15)} = 33.89$, $P < 0.001$ for the effect of AP5 treatment × effect of HFS<sub>vSUB/CA1</sub> protocol interaction; $F_{(1,15)} = 28.60$, $P < 0.0001$ for the effect of the AP5 treatment; Fig. 3f). Under the same conditions, HFS<sub>vSUB/CA1</sub> still induced LTD<sub>ILCx</sub> (2.76 ± 2.76% of baseline, $n = 5$, $P < 0.001$; $F_{(1,13)} = 9.68$, $P < 0.01$ for the effect of AP5 treatment × effect of HFS<sub>vSUB/CA1</sub> protocol interaction; $F_{(1,13)} = 97.77$, $P < 0.0001$ for the effect of HFS<sub>vSUB/CA1</sub> protocol; $F_{(1,13)} = 6.90$, $P < 0.05$ for the effect of AP5 treatment, Fig. 3g). In addition, HFS protocol applied in the ILCx (HFS<sub>ILCx</sub>) was ineffective to trigger plasticity in the ILCx nor in the vSUB/CA1 (Supplementary Fig. 3).

**In vivo LTP in the amBNST triggers anxiolysis.** Finally, we conducted behavioral experiments to test our hypothesis that vSUB/CA1-driven NMDA-receptor-dependent LTP in the amBNST triggers anxiolysis ( timeline Fig. 4a), as both regions are implicated in anxiety-related behaviors²⁴,²⁵. To test whether HFS<sub>vSUB/CA1</sub> promotes anxiolytic effect in basal but also in anxiogenic situation, we used two complementary anxiety-tests based on the innate aversion of rodents to brightly illuminated areas, with the light–dark test (Fig. 4b–f) or to open spaces with the elevated plus maze (EPM; Fig. 4e,g,h)²⁶. In the light–dark test, rats exposed to HFS<sub>vSUB/CA1</sub> spent more time in the light compartment in basal situation compared with SHAM group (percentage of time in light; HFS<sub>vSUB/CA1</sub> group, 36.8 ± 3.8%; SHAM group, 27.5 ± 2.7%; $P < 0.05$, Fig. 4c). Intra-BNST AP5 infusion prevented the anxiolytic-like effect induced by HFS<sub>vSUB/CA1</sub>. AP5 injection into the BNST decreased the % of time spent in light only in the HFS<sub>vSUB/CA1</sub> group ($P < 0.05$ for the effect of AP5 injection, Fig. 4c,d). Neither novelty-induced locomotor activity (Supplementary Fig. 4a) nor circadian rhythms of general activity (Supplementary Fig. 4b) were altered after HFS<sub>vSUB/CA1</sub> or after manipulating the activity on the NMDA receptors in the BNST, thereby reinforcing the specific role of homeostatic plasticity on behavioral anxiety. There was no significant difference in the number of transitions between groups (light–dark test: SHAM/aCSF: 15.20 ± 0.60 transitions, SHAM/AP5: 16.50 ± 1.50 transitions,

<figure><img src="image_4.png" /><figcaption>**Figure 2 | Trans-synaptic tracing from single-electroporated amBNST neurons.** (a) Experimental protocol of *in vivo* single-cell electroporation of amBNST neuron coupled with a recombinant monosynaptic retrograde-pseudotyped rabies virus strategy. A single-amBNST neuron, responding to both ILCx and vSUB/CA1 stimulations, was electroporated *in vivo* with a plasmid coding for the rabies glycoprotein, the TVA receptor and a fluorescent marker (pAAV-EF1a-Rabies G-TVA-tdTomato). Two days later, a pseudotyped rabies (EnvA-SAD ○ G-GFP) was injected into the amBNST and will infect only the electroporated amBNST neuron. (b) Confocal images of vSUB/CA1 (top) and ILCx (bottom) GFP-positive neurons (scale bar, 10 µm).</figcaption></figure>
HFS<sub>vSUB/CA1</sub>/aCSF: 17.20 ± 2.20 transitions, HFS<sub>vSUB/CA1</sub>/AP5: 16.30 ± 1.50 transitions; *P* > 0.05) (Supplementary Fig. 4c). In the light–dark test or in the EPM, the light serves as an anxiogenic stimulus<sup>27</sup>. To create an anxiogenic situation, the lighting was increased from 560 Lux to 1,230 Lux, and rats were restrained for 5 min before the light–dark test and the lighting was set-up at 260 Lux in open arms for the EPM (Fig. 4e)<sup>28</sup>. In the light–dark test, as expected, when rats where submitted to anxiogenic situation, they spent less time in the light compared with the SHAM group in basal situation (SHAM group in basal situation, 27.5 ± 2.7%; SHAM group in anxiogenic situation, 10.16 ± 3.1%,

P < 0.005, Fig 4c,f). In addition, when rats are exposed to anxiogenic situation in the light–dark test or EPM (Fig. 4e–h), HFS<sub>vSUB/CA1</sub> has still an anxiolytic effect (light–dark test: SHAM group in anxiogenic situation, 10.16 ± 3.1%; HFS<sub>vSUB/CA1</sub> group in anxiogenic situation, 19.76 ± 4.2%; *P* < 0.05, Fig. 4f; EPM: percentage of time in the open arms, SHAM group: 14.95 ± 2.51%; HFS<sub>vSUB/CA1</sub> group: 23.38 ± 3.65%; *P* < 0.05, Fig. 4h). No differences were observed between groups for the number of transitions in the light–dark test or EPM (light–dark test: SHAM group: 9.90 ± 0.89 transitions; HFS<sub>vSUB/CA1</sub> group: 11.29 ± 1.60 transitions; *P* > 0.05, Supplementary Fig. 4d; EPM: SHAM group: 16.14 ± 1.47 transitions; HFS<sub>vSUB/CA1</sub> group: 19.29 ± 1.96 transitions; *P* > 0.05, Supplementary Fig. 4e). Taken together, these data indicate that (1) HFS<sub>vSUB/CA1</sub> decreases the steady-state anxiety level through a NMDA-receptor-dependent mechanism in the BNST (2) acute HFS<sub>vSUB/CA1</sub> diminishes the anxiety induced by an anxiogenic situation.

Discussion

Using tract-tracing, trans-synaptic tracing from single-electro-
porated neurons and *in vivo* recordings, we report that the
majority of neurons of the anteromedial BNST integrates
information from the vSUB/CA1 and the ILCx at the single-cell
level. Considering the differences in the straight-line distances
between the ILCx-BNST (≈4 mm) and vSUB/CA1-BNST
(≈6 mm), the 3 ms difference in the latency to the onset of the
stimulation response between the two inputs is probably not
supported by a difference in conduction velocity, but by the
length of the axonal projections. We next found that HFS<sub>vSUB/CA1</sub>
triggers *in vivo* an evoked spike potentiation (LTP<sub>vSUB/CA1</sub>),
which requires the activation of NMDA-Rs in the BNST. This
HFS protocol was efficient to trigger plasticity in the amBNST
when applied in the vSUB/CA1, but not in the ILCx
(Supplementary Fig. 3). This is probably due to the fact that
vSUB/CA1 is one of the few major output structures of the
hippocampal formation and transmits learning and memory-
related signals in a high-frequency bursting mode (Supple-
mentary Fig. 2) repeated at a low frequency (0.5–2 Hz)<sup>29</sup>.
A pioneering study in the hippocampus demonstrated that
changes in synaptic strengths affect network activity and
shape neuronal integration in an input-specific manner<sup>30</sup>. Here,
we report that in response to an excessive activity of the vSUB/
CA1-amBNST inputs (LTP<sub>vSUB/CA1</sub>), amBNST neurons at the
single-cell level down-regulate the efficacy of their ILCx-amBNST
inputs (LTD<sub>ILCx</sub>) and maintain their basal activity stable. This
interaction of homosynaptic plasticity (LTP<sub>vSUB/CA1</sub>) with
heterosynaptic plasticity (LTD<sub>ILCx</sub>) occurring at the single-cell
level at two amBNST excitatory synapses could be considered as
homeostatic plasticity<sup>16</sup>. In fact, these neuroplastic changes
correspond to a form of homeostasis necessary to maintain

<caption>Figure 3. | *In vivo* input-specific opposite plasticity in the amBNST.</caption>

(a) Experimental protocol. (b) Cartography of recording sites in the BNST. Neurons which have been tested for the ILCX are represented in green, those tested for the vSUB are in magenta and those tested for both inputs are in black. (c,d) Kinetic (c) and quantification (d) of the mean percentage change (± s.e.m.) in vSUB/CA1 (magenta) and ILCx (green) evoked spike probability, normalized to the baseline, after HFS<sub>vSUB/CA1</sub>. Rmag, excitatory response magnitude. (e) Typical PSTHs (+ rasters) illustrating responses of a same single-BNST neuron before (left) and after (right) HFS<sub>vSUB/CA1</sub> on the vSUB/CA1-amBNST pathway (top) and on the ILCx-amBNST pathway (bottom). Stimulus: t0 (gray lines). Bin width: 1 ms. Traces in insets. (f,g) Quantification of the mean percentage change (± s.e.m.) in the vSUB/CA1 spike probability (f) or in ILCx spike probability (g) in control condition or after AP5 infusion in the amBNST.

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4** | *In vivo* NMDA-R antagonist infusion in BNST neurons blocks anxiolytic effect induced by HFS<sub>vSUB/CA1</sub>. (a) Timeline of the behavioral experiments and experimental protocol. (b) Experimental protocol in basal situation. (c) Histograms showing the performances at the light-dark box test after AP5 (grey) or vehicle in the BNST (aCSF, white) followed by HFS<sub>vSUB/CA1</sub> or SHAM manipulation in basal situation. (d) Histological control. Scale bar, 1.0 mm. (e) Experimental protocol in anxiogenic situation. (f) Histograms showing the performances at the light-dark box test followed by HFS<sub>vSUB/CA1</sub> (purple) or SHAM manipulation (white) in anxiogenic situation. (g) Example trace of a SHAM rat (left) and a HFS<sub>vSUB/CA1</sub> treated rat (right) performances at the EPM (h) Histograms showing the percentage of time spent in the open arms in the EPM in the SHAM group and in HFS vSUB/CA1 group in anxiogenic situation.</figcaption></figure>

stability in the amBNST, in response to the strengthening of the vSUB/CA1-amBNST synapses (LTP<sub>vSUB/CA1</sub>) associated with a negative feedback at the ILCx-amBNST synapses (LTD<sub>ILCx</sub>). Regardless of the cellular explanation, the functional effect of LTD<sub>ILCx</sub> is complementary to that of LTP<sub>vSUB/CA1</sub>, that is it

optimizes the signal-to-noise ratio by reinforcing the functional weight of the recently potentiated synapses. Interestingly, in presence of AP5 in the amBNST, the direction of plasticity at the vSUB/CA1-amBNST inputs elicited by HFS<sub>vSUB/CA1</sub>, was switch from an LTP<sub>vSUB/CA1</sub> to an LTD<sub>vSUB/CA1</sub> (Fig. 3f). One possibility is that AP5, locally infused in the amBNST, only partially blocks the NMDA receptors, and this partial blockade reverses the direction of plasticity elicited HFS<sub>vSUB/CA1</sub> (refs 31,32). Another unexpected result was that the intra-amBNST blockade of NMDA-Rs potentiates LTD<sub>ILCx</sub> elicited by HFS<sub>vSUB/CA1</sub> (Fig. 3g). One possibility is that is that HFS<sub>vSUB/CA1</sub> triggers concomitant activation of NMDA and metabotropic glutamate receptors in amBNST neurons, leading to a more profound LTD in the presence of AP5 (ref. 33). Further experiments are necessary to determine the molecular mechanisms by which in the absence of NMDA-Rs stimulation, HFS<sub>vSUB/CA1</sub> triggers LTD<sub>vSUB/CA1</sub> and potentiates LTD<sub>ILCx</sub>. Finally, we can not exclude that HFS<sub>vSUB/CA1</sub> also triggers plasticity in the mPFC<sup>23</sup> or the basolateral amygdala<sup>34</sup>, but we provide behavioral evidence that vSUB/CA1-driven NMDA-R-dependent LTP in the amBNST triggers anxiolytic-like effects. This is in line with pioneer studies showing that changing the activity in the amBNST has a direct impact on the perception of aversive contextual stimuli<sup>35</sup> or production of stress hormones<sup>5</sup>. In fact, here we have demonstrated, using two different anxiety assays that HFS<sub>vSUB/CA1</sub> induced an anxiolytic effect in basal situation but also in anxiogenic situation (Fig. 4). Together, these data support the conclusion that the amBNST plays a crucial role in integrating and sending information related to anxiety<sup>9,36</sup>. Previous studies have shown that anxiety is controlled by multiple circuits in the brain, many of which share robust and reciprocal connections with the BNST<sup>4,37</sup>. These circuits include projections from the basolateral nucleus of the amygdala (BLA) to the ventral hippocampus<sup>38</sup>, from BLA to the central nucleus of the amygdala (CeA)<sup>39</sup>, from the ventral hippocampus to the medial prefrontal cortex (mPFC)<sup>40</sup>, from mPFC to BLA<sup>41</sup> and from BLA to BNST<sup>4</sup>. Our anatomical and functional characterization of the vSUB/CA1-amBNST projection on a circuit and synaptic level furthers the understanding of the role played by amBNST in the modulation of anxiety<sup>4,37</sup>. In conclusion, we show that in response to HFS<sub>vSUB/CA1</sub>, homeostasis in amBNST neurons is guaranteed at the single-cell level by an NMDA-R-dependent up-scaling of the vSUB/CA1-amBNST synapses associated with an NMDA-R independent down-regulation of the efficacy of its ILCx-amBNST inputs (LTD<sub>ILCx</sub>; Supplementary Fig. 5). Together these findings elucidate the molecular targets that contribute to changes in synaptic functions in the amBNST, and highlight important future directions where manipulation of inputs to the amBNST using opto- or chemogenetic tools may be critical for changing network output, physiological manifestations of anxiety and anxiety-associated disorders<sup>42</sup>.

## Methods

**Animals.** Male sprague Dawley rats (275–300 g; 10 weeks old; Elevage Janvier, France) were used. Rats were housed three or four per cage under controlled conditions (22–23 °C, 40% relative humidity, 12 h light/dark illumination cycle; lights on from 07:00 hours to 19:00 hours), were acclimatized to laboratory conditions 1 week before the experiment, with food and water *ad libidum*. All procedures were conducted in accordance with the European directive 2010-63-EU and with approval from Bordeaux University Animal Care and Use Committee (N° 50120205-A).

**Surgery.** *In vivo* electrophysiology. Stereotaxic surgeries for electrophysiology, tract-tracing and for light–dark test experiments were performed under 1.0–1.2% isoflurane (in 50% air/50% O₂; 11 min⁻¹) anesthesia<sup>43</sup>. Stimulation electrodes, recording pipettes or injection pipettes were, respectively, inserted into the ILCx (+3.0 mm/bregma, 0.5 mm/midline, 4.5 mm/brain surface), the vSUB/CA1 (−6.0 mm/bregma, 5.1 mm/midline, 7.1 mm/brain surface), the amBNST (0.0 mm/bregma, 1.3 mm/midline, 6–7.5 mm/brain surface).
**Electrical stimulation of the ILCx and the vSUB/CA1.** Bipolar electrical stimulation of the vSUB/CA1 and ILCx was conducted with a concentric electrode (Phymep, France) and a stimulus isolator (500 µs duration, 0.2–2 mA; Digitimer, UK). Baseline was recorded for 10 min (2 × 100 pulses; 0.5 Hz). To induce LTP, high-frequency stimulation protocol was performed in the vSUB/CA1 (HFS<sub>vSUB/CA1</sub>) at the same intensity used for baseline (0.2–1 mA). To avoid the confounding effect of epilepsy driven behavioral changes occurring in awake freely moving animals<sup>23</sup>, HFS was delivered in anesthetized animals. HFS<sub>vSUB/CA1</sub> consisted in 50 trains (500 pulses at 400 Hz, 250 µs duration pulse) presented as bursts of five trains. The frequency of the five trains was 1 Hz. Each burst of five trains was presented five times at 1 min interval (Supplementary Fig. 2).

**amBNST recordings.** A glass micropipette (tip diameter, 1–2 µm; 10–15 MΩ) filled with a 2% pontamine sky blue solution in 0.5 M sodium acetate was lowered into the amBNST. The extracellular potential was recorded with an Axoclamp-2B amplifier and filter (300 Hz/0.5 Hz). Spikes were collected online (CED 1401, SPIKE 2; Cambridge Electronic Design; UK). During electrical stimulation of the ILCx or vSUB/CA1, cumulative peristimulus time histograms (PSTH, 5 ms bin width) of amBNST activity were generated for each neuron recorded.

**Pharmacological treatment.** For local delivery of 100 µM AP5, double barrel pipettes were used<sup>43</sup>. For behaviour, a mixture of 180 nl of AP5 (100 µM) and 0.2% Fluorogold (to mark the injection site) was injected bilaterally in the BNST.

**Histology.** At the end of each recording experiment, the recording pipette placement was marked with an iontophoretic deposit of pontamine sky blue dye (–20 µA, 30 min). To mark electrical stimulation sites, +50 µA was passed through the stimulation electrode for 90 s. After, brains were removed and snap-frozen in a solution of isopentane stored at –80 °C.

**Plasmid solution.** A plasmid encoding for the rabies glycoprotein, the avian virus receptor TVA and a fluorescent marker (tdTomato) was used in this study. For recording followed by electroporation experiments, the electrode was filled with the plasmid (pAAV-EF1a-G-TVA-tdTomato, 17.5 ng µl⁻¹) diluted in standard intracellular solution.

**In vivo single-cell electroporation.** Single-cell electroporation was performed as described previously<sup>44</sup>. After recording amBNST neurons responding to both ILCx and vSUB/CA1 stimulations, they were electroporated with a solution containing a plasmid DNA (pAAV-EF1a-G-TVA-tdTomato). We applied +10 V square-pulses delivered at 50 Hz for 1 s. Only one cell per rat was electroporated. After 2 days, an EnVA pseudotyped G-deleted rabies virus (EnVA-SAD ○ G-GFP) was injected into the amBNST. After 5 days, the electroporation protocol, rats were killed for immunohistological experiments.

**Tract-tracing method.** Tracer injections were performed as described previously<sup>43</sup> with the following modifications. For retrograde tracing, 30 nl of 0.5% CTb (Sigma Aldrich, France) were infused by pressure into the amBNST (n = 5). Animals received a single iontophoretic injection of a 2.5% solution of an anterograde tracer PHAL (Vector Laboratories; UK) in the vSUB/CA1 (n = 5). Animals received a single iontophoretic injection of a 2.5% solution of an anterograde tracer BDA in the ILCx (n = 3). Animals were allowed to survive 7–14 days.

**Immunohistochemical methods.** Immunocytochemical detection was performed by standard light, confocal microscopy as previously described<sup>43</sup>. Rats were perfused transcardially (4% paraformaldehyde solution). Sections were incubated (overnight per 4 °C) with rabbit or goat anti-PHAL primary antibody (1/2,000, Clinisciences, France), rabbit anti-Fox3 primary antibody (1/1,000; Abcam, France), goat anti-CTb primary antibody (1/10,000, List biological laboratories, USA). For fluorescence microscopy, after washing sections were incubated overnight at 4 °C with a donkey anti-rabbit secondary antibody (labelling of PHAL, 1/1,000, Invitrogen, France, Alexa 568), donkey anti-goat secondary antibody (labelling of CTb, 1/1,000, Jacskon Immunoresearch, UK, Alexa 647), biotinylated antibody coupled with an Alexa 488 (labelling of BDA, Vector laboratories) or donkey anti-goat secondary antibody (labelling of PHAL, 1/1,000, Invitrogen, France, alexa 568), donkey anti-rabbit secondary antibody (labelling of Fox3, 1/1,000, Jacskon immunoresearch, UK, Alexa 647). Sections were washed and then mounted in Vectashield medium (Vector Laboratories), coverslipped and imaged on a laser scanning confocal microscope. The confocal microscope was a Leica TCS-SP5-STED on an inverted stand DMI6000 (Leica Microsystems, Germany), using objectives HCX Plan Apo CS × 10 dry NA 0.40, HCX Plan Apo CS × 20 multi-immersion NA 0.70, HCX Plan Apo CS × 63 oil NA 1.40. For conventional confocal microscopy the lasers used were Argon (458 nm, 476 nm, 488 nm, 496 nm, 514 nm) Green Helium-Neon 561 nm, Orange Helium-Neon 594 nm and Red Helium-Neon 633 nm. The scanning was done with a resonant scanner (8,000 Hz–16,000 Hz). The microscope was equipped with 5 internal photomultiplier tube (PMT), two of them are hybrid detectors, 2 NDD PMT,

2 APD and 1 PMT Trans. Photomicrographs were taken and displayed using Image J to perform z stack images.

**Elevated plus maze.** One week later SHAM and HFS stimulations, it has been performed an EPM test. The apparatus is in plexiglass, with four elevated arms arranged in a cross-like disposition, it consisted in two opposite enclosed arms and two open arms, having at their intersection a central square platform which gave access to any of the four arms. Each arm is 50 cm; the wall for the closed arm is 40 cm height. All floor surfaces were grey and the open arms were under an anxiogenic illumination of 260 Lux. In brief, the rats were placed individually in the central square of the EPM facing an open arm and then allowed to start exploring the maze freely during 5 min test. Video recordings were analysed offline using video tracking software (Videotrack from Viewpoint, Lyon, France). The following behaviours were measured: time spent in each compartments and number of transitions between compartment. The percentage of time spent in open arms was evaluated to assess anxiety.

**Light–dark test.** Eight days after the application of HFS in the vSUB/CA1 in a new group of anesthetized rats, we performed the light–dark test. The test lasted 5 min and was performed in a two-compartment box (40 × 40 × 35 cm) with two equal compartments that limit exploratory behavior. An aperture enabled the rats to pass from one compartment to the other. One was completely enclosed by black opaque plastic sides, with a lid of the same material, while the other was white, had no lid, and was brightly illuminated (560 lux). At the start of the experiment a rat was placed in the center of the lighted box with its head facing the wall opposite to the door. The latency for the first emergence from the dark to the light compartment, time spent in each compartment and frequency of explorations of the light compartment were recorded. Time in a zone was considered when the animal placed its four paws in that zone. Rats performed only one time the light–dark test.

A separate group of rats was exposed to anxiogenic situation in order to increase the anxious state of the rats. The intensity of the lighting was increases from 560 Lux to 1,230 lux, and animals were restrained for 5 min, and after 2 min of recovery were exposed to the light–dark test<sup>45</sup>. The same behavioral parameters were measured in anxiogenic situation as in normal situation.

**Data analysis.** For *in vivo* electrophysiological experiments, cumulative PSTHs of aBNST activity were generated during stimulation of ILCx or vSUB/CA1. Excitatory magnitudes ($R_{mag}$ values) were normalized for different levels of baseline impulse activity. Baseline activity was calculated on each PSTH, during the 500 ms preceding the stimulation. For each PSTH, $R_{mag}$ values for excitation were calculated according to: excitation $R_{mag}$ = (number of spikes in excitatory epoch)–(mean number of spikes per baseline bin × number of bins in excitatory epoch). The cortical or hippocampal excitation strength onto amBNST neurons was determined as the amount of current needed to obtain a 50% spike probability for ILCx-evoked responses or vSUB/CA1-evoked responses ($R_{mag}$ ranging from 30 to 60). Results are expressed as mean ± s.e.m. For statistic, two-group comparisons were achieved using Student's *t*-tests or Mann–Whitney when necessary. For multiple comparisons, values were subjected to a two-way Anova followed if significant by Bonferroni post hoc tests or to Kruskall–Wallis Anova for the behavioral part.

**Data availability.** The data that support the findings of this study are available from the corresponding author upon reasonable request.

## References

1. McMenamin, B. W., Langeslag, S. J., Sirbu, M., Padmala, S. & Pessoa, L. Network organization unfolds over time during periods of anxious anticipation. *J. Neurosci.* **34**, 11261–11273 (2014).

2. Somerville, L. H., Whalen, P. J. & Kelley, W. M. Human bed nucleus of the stria terminalis indexes hypervigilant threat monitoring. *Biol. Psychiatry* **68**, 416–424 (2010).

3. Jennings, J. H. *et al.* Distinct extended amygdala circuits for divergent motivational states. *Nature* **496**, 224–228 (2013).

4. Kim, S. Y. *et al.* Diverging neural pathways assemble a behavioural state from separable features in anxiety. *Nature* **496**, 219–223 (2013).

5. Dunn, J. D. Plasma corticosterone responses to electrical stimulation of the bed nucleus of the stria terminalis. *Brain Res.* **407**, 327–331 (1987).

6. Choi, D. C. *et al.* Bed nucleus of the stria terminalis subregions differentially regulate hypothalamic-pituitary-adrenal axis activity: implications for the integration of limbic inputs. *J. Neurosci.* **27**, 2025–2034 (2007).

7. Duvarci, S., Bauer, E. P. & Pare, D. The bed nucleus of the stria terminalis mediates inter-individual variations in anxiety and fear. *J. Neurosci.* **29**, 10357–10361 (2009).

8. Walker, D. L., Miles, L. A. & Davis, M. Selective participation of the bed nucleus of the stria terminalis and CRF in sustained anxiety-like versus phasic fear-like responses. *Prog. Neuropsychopharmacol. Biol. Psychiatry* **33**, 1291–1308 (2009).
9. Cullinan, W. E., Herman, J. P. & Watson, S. J. Ventral subiculum interaction with the hypothalamic paraventricular nucleus: evidence for a relay in the bed nucleus of the stria terminalis. *J. Comp. Neurol.* **332**, 1–20 (1993).

10. Maggio, N. & Segal, M. Persistent changes in ability to express long-term potentiation/depression in the rat hippocampus after juvenile/adult stress. *Biol. Psychiatry* **69**, 748–753 (2011).

11. Conrad, K. L., Louderback, K. M., Gessner, C. P. & Winder, D. G. Stress-induced alterations in anxiety-like behavior and adaptations in plasticity in the bed nucleus of the stria terminalis. *Physiol. Behav.* **104**, 248–256 (2011).

12. Glangetas, C. et al. Stress switches cannabinoid type-1 (CB1) receptor-dependent plasticity from LTD to LTP in the bed nucleus of the stria terminalis. *J. Neurosci.* **33**, 19657–19663 (2013).

13. Pego, J. M. et al. Dissociation of the morphological correlates of stress-induced anxiety and fear. *Eur. J. Neurosci.* **27**, 1503–1516 (2008).

14. O'Brien, R. J. et al. Activity-dependent modulation of synaptic AMPA receptor accumulation. *Neuron* **21**, 1067–1078 (1998).

15. Turrigiano, G. G., Leslie, K. R., Desai, N. S., Rutherford, L. C. & Nelson, S. B. Activity-dependent scaling of quantal amplitude in neocortical neurons. *Nature* **391**, 892–896 (1998).

16. Zenke, F., Agnes, E. J. & Gerstner, W. Diverse synaptic plasticity mechanisms orchestrated to form and retrieve memories in spiking neural networks. *Nat. Commun.* **6**, 6922 (2015).

17. Dong, H. W., Petrovich, G. D. & Swanson, L. W. Topography of projections from amygdala to bed nuclei of the stria terminalis. *Brain Res. Brain Res. Rev.* **38**, 192–246 (2001).

18. Massi, L. et al. Cannabinoid receptors in the bed nucleus of the stria terminalis control cortical excitation of midbrain dopamine cells in *vivo*. *J. Neurosci.* **28**, 10496–10508 (2008).

19. Vertes, R. P. Differential projections of the infralimbic and prelimbic cortex in the rat. *Synapse* **51**, 32–58 (2004).

20. Marshel, J. H., Mori, T., Nielsen, K. J. & Callaway, E. M. Targeting single neuronal networks for gene expression and cell labeling in vivo. *Neuron* **67**, 562–574 (2010).

21. Wickersham, I. R., Finke, S., Conzelmann, K. K. & Callaway, E. M. Retrograde neuronal tracing with a deletion-mutant rabies virus. *Nat. Methods* **4**, 47–49 (2007).

22. Abraham, W. C. et al. Correlations between immediate early gene induction and the persistence of long-term potentiation. *Neuroscience* **56**, 717–727 (1993).

23. Taylor, C. J., Ohline, S. M., Moss, T., Ulrich, K. & Abraham, W. C. The persistence of long-term potentiation in the projection from ventral hippocampus to medial prefrontal cortex in awake rats. *Eur. J. Neurosci.* **43**, 811–822 (2016).

24. Herman, J. P. & Mueller, N. K. Role of the ventral subiculum in stress integration. *Behav. Brain Res.* **174**, 215–224 (2006).

25. Crestani, C. C. et al. Mechanisms in the bed nucleus of the stria terminalis involved in control of autonomic and neuroendocrine functions: a review. *Curr. Neuropharmacol.* **11**, 141–159 (2013).

26. Crawley, J. & Goodwin, F. K. Preliminary report of a simple animal behavior model for the anxiolytic effects of benzodiazepines. *Pharmacol. Biochem. Behav.* **13**, 167–170 (1980).

27. Lister, R. G. Ethologically-based animal models of anxiety disorders. *Pharmacol. Ther.* **46**, 321–340 (1990).

28. Pare, D. & Gaudreau, H. Projection cells and interneurons of the lateral and basolateral amygdala: distinct firing patterns and differential relation to theta and delta rhythms in conscious cats. *J. Neurosci.* **16**, 3334–3350 (1996).

29. Sharp, P. E. & Green, C. Spatial correlates of firing patterns of single cells in the subiculum of the freely moving rat. *J. Neurosci.* **14**, 2339–2356 (1994).

30. Abraham, W. C. & Goddard, G. V. Asymmetric relationships between homosynaptic long-term potentiation and heterosynaptic long-term depression. *Nature* **305**, 717–719 (1983).

31. Bains, J. S., Longacher, J. M. & Staley, K. J. Reciprocal interactions between CA3 network activity and strength of recurrent collateral synapses. *Nat. Neurosci.* **2**, 720–726 (1999).

32. Cummings, J. A., Mulkey, R. M., Nicoll, R. A. & Malenka, R. C. $Ca^{2+}$ signaling requirements for long-term depression in the hippocampus. *Neuron* **16**, 825–833 (1996).

33. Huang, C. C. & Hsu, K. S. Activation of NMDA receptors reduces metabotropic glutamate receptor-induced long-term depression in the nucleus accumbens via a CaMKII-dependent mechanism. *Neuropharmacology* **63**, 1298–1307 (2012).

34. Sotres-Bayon, F., Sierra-Mercado, D., Pardilla-Delgado, E. & Quirk, G. J. Gating of fear in prelimbic cortex by hippocampal and amygdala inputs. *Neuron* **76**, 804–812 (2012).

35. Haufler, D., Nagy, F. Z. & Pare, D. Neuronal correlates of fear conditioning in the bed nucleus of the stria terminalis. *Learn Mem* **20**, 633–641 (2013).

36. Radley, J. J. & Sawchenko, P. E. A common substrate for prefrontal and hippocampal inhibition of the neuroendocrine stress response. *J. Neurosci.* **31**, 9683–9695 (2011).

37. Herman, J. P. Neural pathways of stress integration: relevance to alcohol abuse. *Alcohol Res.* **34**, 441–447 (2012).

38. Felix-Ortiz, A. C. et al. BLA to vHPC inputs modulate anxiety-related behaviors. *Neuron* **79**, 658–664 (2013).

39. Tye, K. M. et al. Amygdala circuitry mediating reversible and bidirectional control of anxiety. *Nature* **471**, 358–362 (2011).

40. Adhikari, A., Topiwala, M. A. & Gordon, J. A. Single units in the medial prefrontal cortex with anxiety-related firing patterns are preferentially influenced by ventral hippocampal activity. *Neuron* **71**, 898–910 (2011).

41. Likhtik, E., Stujenske, J. M., Topiwala, M. A., Harris, A. Z. & Gordon, J. A. Prefrontal entrainment of amygdala activity signals safety in learned fear and innate anxiety. *Nat. Neurosci.* **17**, 106–113 (2014).

42. Luyten, L., Hendrickx, S., Raymaekers, S., Gabriels, L. & Nuttin, B. Electrical stimulation in the bed nucleus of the stria terminalis alleviates severe obsessive-compulsive disorder. *Mol. Psychiatry* **21**, 1272–1280 (2015).

43. Georges, F. & Aston-Jones, G. Activation of ventral tegmental area cells by the bed nucleus of the stria terminalis: a novel excitatory amino acid input to midbrain dopamine neurons. *J. Neurosci.* **22**, 5173–5187 (2002).

44. Oyama, K. et al. Long-lasting single-neuron labeling by *in vivo* electroporation without microscopic guidance. *J. Neurosci.* Methods **218**, 139–147 (2013).

45. Haas, D. A. & George, S. R. Single or repeated mild stress increases synthesis and release of hypothalamic corticotropin-releasing factor. *Brain Res.* **461**, 230–237 (1988).

## Acknowledgements

The microscopy was done in the Bordeaux Imaging Center a service unit of the CNRS-INSERM and Bordeaux University, member of the national infrastructure France BioImaging. We thank Philippe Legros for technical support and precious advices for confocal images acquisition. We thank Francis Chaouloff and Mario Carta for scientific discussion. We also thank for Ed Callaway for generously sharing materials. This work was supported by grants from Centre National de la Recherche Scientifique (CNRS), University of Bordeaux, Agence Nationale de la Recherche (ANR-12-BSV4-0022 to F.G.), by LABEX BRAIN ANR-10-LABX-43 and Region Aquitaine and by a NARSAD Young Investigator Grant from the Brain & Behavior Research Foundation to Léma Massi.

## Author contributions

C.G., S.C. and F.G. designed the experiments, L.M., C.X. and A.L. designed and performed the *in vivo* single-cell electroporation experiments. K.Y. and B.R. made the plasmid used for electroporation. C.G., M.J., G.R.F., S.C., L.M. and D.G. collected and analyzed the data. C.G., M.D., S.C., L.M., A.L., G.R.F. and F.G. wrote the manuscript. All authors discussed the results and commented the manuscript.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Glangetas, C. et al. NMDA-receptor-dependent plasticity in the bed nucleus of the stria terminalis triggers long-term anxiolysis. *Nat. Commun.* **8**, 14456 doi: 10.1038/ncomms14456 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this

article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017