In order to better understand what aspects of defense are encoded in VMHvl neuron firing, we used in vivo microendoscopic calcium imaging to measure neuronal response properties in mice subjected to social defeat (Figure 1a–e). Mice were habituated to a home chamber from which they were given access to a corridor and far chamber for a brief period each day. On the social defeat day, mice were closed into the far chamber and an aggressive mouse was introduced. Following social defeat, the far chamber door was opened to allow the mouse to escape and exhibit approach-avoidance behavior toward the aggressor who remained restricted to the far chamber. Many neurons (100/246, 40%, Social+) showed a tonic increase in activity during the social defeat phase that returned to baseline levels.

Figure 1 with 2 supplements see all

Download asset Open asset

Video 1

Download asset

Next, we examined neural activity patterns during the post-defeat, approach-avoidance phase. The mice repeatedly advanced in a cautious manner toward the aggressor, exhibiting frequent risk assessment behaviors in which they stretched their body in the direction of the far chamber (stretch-attend or stretch-approach, Blanchard et al., 2011). Once close to the far chamber mice often turned and fled back to the home chamber (flight, Video 1). Calcium imaging identified many cells that were robustly activated during risk assessment either close to or far away from the far chamber (160/326, 49%, Assessment+). Smaller sets of cells were either deactivated during risk assessment (63/326, 19%, Assessment-) or activated or deactivated during flight (30/198, 15%, Flight+; 60/198, 30%, Flight-; Figure 1l–q). Notably, Assessment+ cells were activated during risk assessment events that occurred at a distance from the far chamber, particularly at the junction of the home cage and corridor (Figure 1o; Video 2). A comparison of neuronal response properties revealed that a majority of Assessment+ cells (81/152, 53% vs chance 17%, p=4.23×10⁻¹¹) overlapped with Defeat+ cells, suggesting that they may encode a generalized internal state associated with both direct threat as well as the assessment of threat even when this does not involve close social contact (Figure 1r). Moreover, many Assessment+ cells (37/100, 37% vs chance 14%, p=0.00031) were also Flight- cells, showing activation as the mouse approached the far chamber and then turning off abruptly when the animal fled (Figure 1—figure supplement 1s; Video 1). A similar firing response pattern has been reported for neurons in VMHdm during approach toward a predator (Masferrer et al., 2018). Overall, Defeat+, Assessment+, Flight- cells showed a high degree of overlap (26% vs chance 5%, p=2.0×10⁻⁶) and made up the largest fraction of responsive neurons (Figure 1t). These findings suggest that VMHvl neurons may encode a generalized social threat state rather than particular behaviors associated with threat response.

Video 2

Download asset

Our observation that a large fraction of neurons in VMHvl were activated during risk assessment both close and far from the social stimulus could be explained by the multi-sensory inputs that VMHvl receives (Canteras et al., 1995; Garfield et al., 2014; Lo et al., 2019; Wong et al., 2016). In particular, our earlier observation that VMHdm is required for the expression of defensive behaviors in a context previously associated with a predator suggested that the medial hypothalamus can be recruited by stimulus-associated cues (Silva et al., 2016a) and is consistent with evidence for activation of VMHvl during nose poke in anticipation of aggression (Falkner et al., 2016). To test for the recruitment of VMHvl by purely contextual cues, we performed in vivo calcium endoscopy in animals who experienced social defeat and were re-exposed 24 hr later to the defeat context in the absence of the aggressor (Figure 2a). Analysis of neuronal response properties revealed a large fraction of neurons with robust activation in the defeat chamber (93/343, 27%, Far chamber+, Figure 2b–d). Unexpectedly, a second group of neurons showed marked activation in the home chamber (78/343, 23%, Home+, Figure 2e–g). Notably, Home+ and Far chamber+ cells abruptly turned off when the animal passed from the chamber to the corridor, suggesting that they specifically responded to the far and home chamber contexts, rather than to a gradient of social threat cues, for example. A comparison of neuronal response properties between the context, approach-avoidance, and defeat phases of the test showed that Far chamber+ neurons strongly overlapped with Social+ neurons (45/67, 67% vs chance 11%, p=9.597×10⁻¹¹) and Home+ neurons with Social- neurons (37/60, 62% vs chance 7%, p=1.065×10⁻¹⁰, Figure 2i,j and Figure 2—figure supplement 1) and that Far chamber+ neurons also overlapped with Defeat+ neurons (35/84, 42% vs chance 10%, p=2.0×10⁻⁶; Figure 2k) and, in a three-way comparison, a larger fraction of Far chamber+ than Home+ neurons were Assessment+ or Defeat+ (Figure 2l,o). These data suggest that VMHvl encodes a generalized state of social threat that can be reactivated by contextual cues (Social+, Defeat+, Far chamber+) as well as features of social territory (Home+).

Figure 2 with 3 supplements see all

Download asset Open asset

To understand whether the recruitment by context was linked to the past experience of the animal, we compared context encoding in animals before and after social defeat. Only a small number of Far chamber+ and Home+ cells could be identified during the pre-defeat habituation phase even though the animal explored the apparatus to a similar extent in both phases (Figure 2—figure supplement 2c). To quantify changes in neuron activity before and after defeat, we employed linear discriminant analysis (LDA) of neuron ensemble activity. Before defeat, neuron ensemble activity in the home, corridor, and far chambers was overlapping as visualized using principal LDA discriminants (Figure 2m,n and Figure 2—figure supplement 2a). Following defeat, however, home, corridor, and far chamber ensemble activity became significantly more separated (Figure 2p,q and Figure 2—figure supplement 2b; Home vs. Far chamber LDA distance, p=0.0075) confirming the observation that social defeat induced a marked enhancement of the encoding of spatial context and suggesting that VMHvl may dynamically encode features of territory.

Having established that VMHvl encodes features of social threat and context we examined the hypothesis that social threat neurons would also be active during aggression. First, we performed serial cFos labeling using TRAP-tagging (cFos::CreERT2; Rosa26::LSL-tomato; Guenthner et al., 2013) to determine the extent of overlap in recruitment during social defeat and resident-intruder aggression. Naive animals were subjected to two resident-intruder tests at one week interval in which they were either the resident or intruder, in all possible combinations (Aggression-Aggression, Defense-Defense, Aggression-Defense, Defense-Aggression). Immediately following the first test mice were treated with 4-OHT to induce persistent labeling of cFos+ cells, followed by immunolabeling of cFos+ cells after the second test. Animals subjected to the same behavioral experience (Aggression-Aggression or Defense-Defense) showed 50 ± 15% overlap of cFos+ cells while those with different experiences (Aggression-Defense or Defense-Aggression) showed 18 ± 5% overlap and those under control conditions 10 ± 2% overlap (ANOVA: F = 20.89, p=0.0001; Figure 3a–d). These data are consistent with experiments using a viral cFos-tagging strategy (Sakurai et al., 2016) and suggest that approximately one quarter of the cells in VMHvl activated during aggression and defense are recruited by both experiences and that this population of cells may support a common function during male-male social interaction.

Figure 3

Download asset Open asset

Next, we used in vivo calcium endoscopy to understand whether VMHvl neurons recruited during aggression might overlap with the social threat and context neurons identified following defeat. Robust aggression was elicited by confining the mouse to its home cage and introducing a subordinate male mouse. During aggression phasic modulation of neuron activity occurred during bouts in which the resident actively investigated (56/310, 18% Sniff+; Figure 3e–g) or attacked (38/266, 14% Attack+; Figure 3h–j) the intruder. Comparison of neuronal response properties revealed that most Attack+ neurons (19/37, 51% vs chance 3%, p=2×10⁻⁶) were Sniff+ neurons (Figure 3k), consistent with a strong overlap between these populations in previous single unit recording studies (Lin et al., 2011; Falkner et al., 2014). Comparison of neuronal response properties across social defeat and aggression showed that a large fraction of Defeat+ neurons (24/68, 35% vs chance 6%, p=0.00003) were reactivated during active investigation of the subordinate (Sniff+), indicating that these may encode a common social threat state in both defenders and attackers (Figure 3l). A smaller fraction of Defeat+ neurons (12/60, 20% vs chance 5%, p=0.02458) were Attack+ neurons (Figure 3m) and a three-way comparison revealed that 42% of Defeat+ cells overlap with either Sniff+ or Attack+ (Figure 3n), a finding that is consistent with our cFos data (Figure 3a–d) and confirms the recruitment of both common and unique neuron ensembles during defense and aggression. Comparison of aggression responsive neurons with those activated by the defeat context (Far chamber+) revealed an overlap between Attack+, Defeat+, and Far chamber+ neurons (Figure 3o). Unexpectedly, Attack+ neurons showed more overlap with Far chamber+ (14/34, 41% vs chance 4%, p=0.00022) than Home+ neurons (4/34, 12% vs chance 3%, p=0.35591; Figure 3p) despite the fact that the attack occurred in the home cage, reinforcing the idea that aggression elicits features of social threat, rather than safety.

Following on our observation that social defeat elicited changes in contextual encoding in VMHvl (Figure 2) and previous experiments showing that sexual experience can remodel VMHvl ensemble activity elicited by exposure to a social stimulus (Remedios et al., 2017), we investigated the possibility that social defeat might also induce changes in neural activity elicited during social encounter. Animals were exposed to repeated social defeat experiences in the far chamber of a dual chamber apparatus as described earlier (Figure 1e and Figure 1—figure supplement 1). The experimental animals showed similar frequencies of behaviors across the two days of social defeat (Figure 4—figure supplement 1). Nevertheless, to minimize any confounding effect of potential changes in behavior across experimental sessions, we restricted our ensemble analysis to data extracted from bouts of close social encounter (defeat, upright, orientating, sniffing, following, see Materials and methods). Linear discriminant analysis (LDA) was used to track and quantify shifts in neuron ensemble encoding across sessions. Neuron ensemble activity elicited during the first and second social defeat experiences were significantly non-overlapping (Figure 4ab and Figure 4—figure supplement 2a). As a control, we used LDA to track and quantify neuron ensemble activity across sessions of resident-intruder aggression carried out as previously described (Figure 4—figure supplement 2c). Unlike social defeat and consistent with previous studies (Remedios et al., 2017) repeated aggression was not associated with changes in neuron ensemble activity (Figure 4ab and Figure 4—figure supplement 2a) demonstrating that defeat has a unique capacity to alter the encoding of social cues in VMHvl. Notably, this change was found only when defeats preceded, but not when they followed aggression experiences (Figure 4—figure supplement 2bc). Notably, the changes in ensemble activity across defeat days were not associated with a change in defensive behaviors, although, as expected, upright postures were significantly lower in the group that experienced aggression first (Figure 4—figure supplement 1). This finding demonstrates that changes in neuron ensemble encoding seen across days are likely due to experience-dependent plasticity rather than spontaneous fluctuations in firing patterns, and suggests that repeated exposure to aggression and winning may impart resistance to the transforming effects of social defeat on social encoding in VMHvl.

Figure 4 with 3 supplements see all

Download asset Open asset

A previous in vivo calcium endoscopy study found activation of estrogen receptor α expressing (Esr1+) neurons in VMHvl during social interaction (Remedios et al., 2017) and a recent bulk calcium imaging study (Wang et al., 2019) as well as single-cell RNA-Seq study (Kim et al., 2019) reported activation of Esr1+ neurons during both social aggression and defeat suggesting that this subpopulation of VMHvl neurons may contribute to the encoding of social threat we observed in the general population of VMHvl neurons. To test this hypothesis, we carried out in vivo calcium endoscopy in mice expressing GCaMP6 selectively in VMHvl Esr1+ neurons. Single unit neural activity was assessed from mice during social defeat, post-defeat approach/avoidance, and during re-exposure to the defeat chamber as described earlier (Figure 2a). During the social defeat phase of the test, many Esr1+ neurons showed tonic increases (47/102, 46%, Social+) or decreases (18/102, 18%, Social-) that returned to baseline levels during the subsequent approach-avoidance phase (Figure 4c). A subset of Esr1+ neurons showed phasic responses that were time-locked to individual defeat events (32/102, 31%, Defeat+; 11/102,11%, Defeat-; Figure 4d). During the subsequent approach/avoidance phase, many Esr1+ cells were robustly activated (63/98, 64%, Assessment+) or deactivated (7/98, 7%, Assessment-) during risk assessment and flight behavior (16/98, 16%, Flight+; 15/98, 15%, Flight-; Figure 4e–f). Overall, Defeat+ and Assessment+ cells showed a moderate, but not significant degree of overlap (32% vs chance 20%, p=0.44031) and made up a major fraction of responsive neurons (Figure 4g). These findings demonstrate that VMHvl Esr1+ neurons encode a generalized social threat state in a manner largely similar to the general VMHvl neuron population.

Analysis of neuronal response properties during context re-exposure revealed a large fraction of neurons with robust activation in the far chamber (37/106, 35%, Far chamber+, Figure 4h) or home (19/106, 18%, Home+, Figure 4i) chambers. A comparison of neuronal response properties between the context and approach-avoidance phases of the test showed that a larger fraction of Far chamber+ than Home+ neurons were Assessment+ (Figure 4j) and linear discriminant analysis of ensemble neural encoding in Esr1+ cells showed a clear separation of firing patterns in Home and Far chambers (Figure 4k). These data confirm that VMHvl Esr+ neurons encode features of social territory in a similar way to the general VMHvl population.

Finally, to understand whether social defeat could remodel neuronal ensemble activity in VMHvl Esr1+ neurons animals were exposed to two consecutive rounds of social defeat as described earlier (Figure 1e). Linear discriminant analysis (LDA) of neuron ensemble activity extracted from bouts of close social encounter (defeat, upright, orientating, sniffing, following, see Materials and methods) revealed that neuron ensemble activity elicited during the first and second social defeat experiences were significantly non-overlapping (Figure 4l and Figure 4—figure supplement 3c; Kolmogorov-Smirnov test, p=0.0001 for each mouse). Machine learning using Esr1+ cells data achieved similar levels of accuracy in distinguishing defeat days as using the data from all cells (Figure 4—figure supplement 3d). These data confirm that Esr1+ ensemble activity patterns are remodeled by social defeat in a manner indistinguishable from the general VMHvl population.

Changes in neuron ensemble encoding in VMHvl following social defeat could reflect plasticity in upstream brain regions that provide afferent inputs to VMHvl or they could be due to plasticity within VMHvl. These possibilities can be distinguished using a gain-of-function approach. Artificial activation of VMHvl before and after social defeat should result in different behavioral outputs if plasticity occurs within or downstream of VMHvl, but not if it is restricted to afferent inputs. Previous work showed that optogenetic activation of VMHvl Esr1+ neurons can elicit aggression against females or castrated males (Lee et al., 2014; Wang et al., 2019). However, in some cases, such animals showed brief avoidance responses during the initial stimulation trials (Lee et al., 2014) suggesting that VMHvl Esr1+ neurons are capable of promoting both aggression and avoidance, but that this may occur in in a manner that depends on context or past social experience.

To test whether social defeat might induce changes in VMHvl Esr1+ function, we expressed the blue light-activated cationic channel, channelrhodopsin 2 (ChR2), in Esr1+ VMHvl neurons (AAV-Ef1a::FLEX-ChR2, Esr1::Cre mice) and subjected the animals to two sessions of social defeat (Figure 5a–c). Several hours after social defeat, a subordinate intruder mouse was introduced into the home cage of the animal and following a period of habituation, light pulses were delivered to activate Esr1+ neurons. Following social defeat optogenetic activation of Esr1+ neurons elicited rapid and robust defensive behaviors, including freezing, flight, a reduction in social interaction, and immobility in all animals as well as in a subset of mice (N = 3) that had not been previously stimulated (Figure 5d). In one particularly striking case, optogenetic activation of Esr1+ neurons occurring while the resident was attacking the intruder caused an abrupt cessation of aggression and evoked sudden flights away from the subordinate animal (Video 3). Notably, no significant increase in defensive behavior was elicited by optogenetic activation in ChR2-expressing mice that were stimulated during a resident-intruder test on the day before the social defeat (Figure 5e), nor in mice expressing a control virus (AAV-Ef1a::FLEX-YFP, Esr1::Cre). These data demonstrate that VMHvl neurons promote defense in a manner that depends on social experience and argue that at least part of the plasticity in neuron ensemble recruitment during social interaction (Figure 4ab) occurs within or downstream of VMHvl.

Figure 5

Download asset Open asset

Video 3

Download asset

All experimental procedures involving the use of animals were carried out in accordance with EU Directive 2010/63/EU and under approval of the EMBL Animal Use Committee and Italian Ministry of Health License 541/2015-PR to C.G. Animals were maintained in a temperature and humidity controlled environment with food and water provided ad libitum and 12 hr/12 hr light-dark cycle with lights on at 7:00. Experimental male C57BL/6J wild type (Charles River) and Esr1::Cre (Stock No. 017913, Jackson Laboratory) mice were switched to reverse dark-light cycle (lights off at 9:00) and singly housed in the home cage of the behavioral apparatus at least 1 week before initiating the experimental procedures and tested during the dark period under red lighting (two 1W LED lights). The custom Plexiglas behavioral apparatus consisted of three detachable parts: 1) home cage with dimensions 25 × 25×25 cm with a Y-shaped slit, 4 cm wide at the bottom serving as an entrance closed by sliding doors, 2) stimulus chamber identical to home cage, and 3) a corridor 46 × 12×30 cm connecting home cage and stimulus chamber (modified from Silva et al., 2013). Aggressor mice were CD-1 adult retired breeders (Charles River) screened for robust aggression and singly housed (Franklin et al., 2017). Subordinate mice were 9–15 week old BALB/c males bred at EMBL. Mice cages were changed weekly.

Mice were anesthetised before surgery with 3% isoflurane (Provet) in oxygen and placed in stereotaxic frame (Kopf). Anesthesia was maintained with continuous 1–2% isoflurane administration in breathing air enriched with oxygen. Body temperature was maintained with a heating pad. During surgery the skull was exposed, aligned, and cleaned with 0.3% hydrogen peroxide solution. For optogenetic activation experiments, 0.1–0.2 µl of AAV5-Ef1a::DIO-hChR2(E123T/T159C)-EYFP or AAV5-Ef1a::DIO-EYFP virus (UNC Vector Core) was infused bilaterally into VMHvl (from Bregma L:+ / - 0.67, A/P: −0.98, D/V: −5,75). After 5–10 min, the glass capillary was retracted and custom-made optic fibre connectors were implanted (0.66 NA, 200 µm core fibre and ceramic ferule with 230 µm/1250 µm internal/external diameter; from Bregma L: + / - 0.67, A/P: −0.98, D/V: −5.55 and L: + / - 1.14, A/P: −0.98, D/V: −5.6, at 5° angle). For in vivo calcium endoscopy 0.2–0.3 µl of AAV5-hSyn::GCaMP6s or AAV5-CAG::Flex-GCaMP6s (Penn Vector Core) virus was injected unilaterally into VMHvl and the endoscope lens (Snap-imaging cannula model L type E, Doric Lenses) was implanted at a very slow rate (from Bregma L: + / - 0.67, A/P: −0.98, D/V: −5.7). All implants were secured to the skull using miniature screws (RWD) and dental cement (Duralay). The wound was cleaned and skin was stitched around the implant. After the surgery mice received intra-peritoneal injection of 0.4 ml saline and were placed into heated cages with drinking water containing paracetamol for ~1 week. Mice were maintained in isolation for ~4 weeks before experimentation.

On the social defeat day, the mouse was allowed to explore the apparatus freely for 5 min, after which the animal was closed in the stimulus chamber and an aggressor was introduced for 10 min. In a few cases, the animal was allowed to escape earlier to avoid excessive defeat. The mouse was released from the stimulus chamber and the door closed to confine the aggressor to the stimulus chamber. After social defeat the mouse was allowed to explore the apparatus freely for at least 5 min. The memory test was conducted on the day following after social defeat, during which the mouse could explore the apparatus freely. The apparatus was washed with detergent and 50% alcohol between mice to avoid any remaining smell that could influence behavior of the test subject. CD-1 aggressors were screened for aggression for 3 days. Every day an intruder was placed in the aggressor’s cage for 3 min. Only mice that attacked the intruder on every occasion were selected.

On the aggression day a BALB/c intruder was introduced for 10 min after which the intruder was removed. BALB/c intruders used for this test were 9–15 weeks old and housed 3–5 mice per cage. For cFos experiments, animal were allowed to explore the entire apparatus for 5 min before and after introduction of the intruder in order to match exploration in the social defeat test.

Behavior was recorded at 40–50 Hz from above with up to two cameras (acA1300-60gmHIR, Basler) with GigE connections using Pylon software. Frame-by-frame behavioral annotation was carried out manually using Observer XT11 (Noldus) and Solomon Coder software. The experimenter was blind to genotype or calcium trace of recorded neurons when scoring behavior. The following behaviors were scored: defeat – biting attack toward the experimental animal in which the animal exhibited avoidance behavior (whole body movement away from intruder); assessment – stretch attend or stretch approach behavior in which the animal extended its body in the direction of the threat or threat chamber from an immobile or slowly moving position; flight – rapid movement away from the threat or threat chamber; attack – biting attack or vigorous anogenital sniffing toward the intruder animal; sniffing – close contact of the nose of the experimental animal with the intruder; upright – animal rises on back paws with head up, keeping front paws stretch out; put down – keeping other animal down usually with two front paws and staying on top of it; follow – experimental animal closely follows intruder; cornering – staying in the corner of the apparatus (animal has body contact with two walls); locomotion – free movement and exploration of experimental apparatus; freeze – no body movement; head/body orientation – turning head or whole body towards another conspecific. Any animals with misplaced viral infections or optic fiber implants were excluded from the analysis.

All animals were deeply anesthetized and transcardially perfused with PBS (Invitrogen) followed by 4% PFA (Sigma) in 0.1M PB solution and then post-fixed in 4% PFA at 4°C for 24 hr and cut using a vibratome (Leica VT 1000 s) in PBS (80 µm) for injection or implant location check, or staining (50 µm) procedure. If not used immediately, sections were stored in PBS with 0.1% sodium azide. For anti-cFos staining (SC-52G, Lot FO215, Santa Cruz) sections were washed three times in PBS for 10 min, blocked with 10% normal donkey serum, 0.2% Triton-X in PBS for at least 1 hr, and incubated with primary antibody (1:500) containing 5% normal donkey serum, 0.2% Triton-X overnight at 4°C. Sections were washed three times in PBS, incubated with secondary antibody solution (1:1000) containing 10% normal donkey serum for 2 hr at room temperature, washed twice with PBS, and stained with DAPI for 15 min, before washing twice with PBS and mounting on SuperFrost Plus slides (ThermoFischer) with Moviol.

Double heterozygous cFos::CreERT2;RC::LSL-tdTomato mice (Stock No. 021882 and 007914, Jackson Laboratory) were isolated for 7 days and subjected to 3 days of habituation to the apparatus before social defeat. A single dose (50 mg/kg) of 4-hydroxytamoxifen (4-OHT, 70% z-isomer, Sigma) was injected intra-peritoneally <2 min after social defeat to induce fluorophore expression in neurons recruited by defeat. One week later, mice were again habituated for 3 days to the apparatus and subjected to social aggression, and 1.5 hr later trans-cardially perfused and processed for cFos immunofluorescence. In a second group, the order of social defeat and social aggression were swapped to produce five experimental groups: defense-defense (same), aggression-aggression (same), defense-aggression (different), aggression-defense (different), defense-control (control). Control indicates mice exposed only to context prior to labelling.

Heterozygous Esr1::Cre mice were injected with virus and allowed to recover for 2–3 weeks, housed under reverse light-cycle for 2 weeks, and then handled and habituated to the optical cables for at least 2 days prior to testing. Initially, after a 2-min free exploration period, each animal was stimulated with light (30 s, 20 Hz) every 1–2 min with increasing power (0.5, 1, 3, 6, 10 mW) to identify the optimal intensity to elicit a behavioral response (immobility or locomotion) and this was subsequently used for further stimulation (3 mW for most animals, with some responding at 1 mW and 6 mW). Control animals received 6 mW stimulation. For stimulation in the presence of a subordinate intruder-free exploration was allowed for ~3 min during which a single light stimulus was delivered (30 s, 20 Hz) after which a BALB/c intruder was introduced and the animal was further stimulated (3–5 times, 30 s, 20 Hz) during periods of social interaction. Mice were then subjected to 2 days of social defeat (5 min) by a CD-1intruder in a novel plexiglas cage following 1 min free exploration. Several hours later, stimulation in the presence of a subordinate intruder was repeated as above. To control for possible repeated optogenetic activation effects a subset of animals from the ChR2 group (N = 3) was not stimulated before defeat episodes. Optical stimulation of ChR2 was achieved using a 465 nm LED (Plexbright, Plexon) attached to a manual rotatory joint with 1 m patch cables (Plexbright High Performance, Plexon). Power at the end of the patch cables was measured before each experiment with a portable optical power meter (Thor Labs) and stimulation trains were generated using V2.2 Radiant software (Plexon). Animals with mistargeted viral injections or optic fibers were excluded from the analysis.

Following GCaMP6 virus infection, GRIN lens imaging canulae (Model L, Doric Lenses) were stereotaxically implanted over the VMHvl and mice allowed to recover for 2–3 weeks followed by 2 weeks isolation under reverse light-cycle and at least 3 days habituation to the microscope plugging and unplugging procedure. Mice were then habituated to the behavioral apparatus for 3 days, and for the second and third days the microscope body was attached. The following week mice (N = 4) were subjected to two consecutive social defeats tests (day 1: social defeat, day 2: memory) over 4 days, and 2–3 days later, three consecutive aggression tests over 3 days. In a second group of mice (N = 3), the testing order as reversed. Esr1 mice (N = 3) were subjected to two consecutive social defeats tests (day 1: social defeat, day 2: memory) over 4 days. Recordings were done with 15–25% LED intensity using 50 ms or 100 ms exposure times. Calcium imaging data was successfully collected from 19/32 animals that underwent surgery. Animals with mistargeted endoscope placements, very few (<15) recorded neurons, recordings exhibiting excessive movement artifacts, or who showed insufficient instances of relevant behaviors were excluded from the analysis (9/19).

Image processing was done using Fiji software. Briefly, videos were first loaded as a stack of images in. tiff format. The stack was duplicated and for each frame a background frame was generated using a band pass filter (lower band: 100, higher band: 10,000). Next, each frame of the recording was divided by the corresponding background frame and the resulting background-filtered stack was aligned using TurboReg via the translation batch algorithm. Neuronal ROIs were manually selected from the maximum intensity projection image, with detection aided by inspecting recordings at increased speed to discern dim neurons with slow dynamics. Finally, the mean intensity ROI traces were extracted and ΔF/F was calculated where F is the mean intensity over the entire recording period. To track ROIs over different recordings and days, an ROI mask was projected onto each new recording and translated if necessary to account for possible field of view movement between recordings (Figure 1—figure supplement 1). ROIs that could not be assigned to the mask were treated as new ROIs.

For the analysis of neuronal response properties, receiver operating characteristic (ROC) curves were calculated using custom Python scripts (scikit-learn) (Krzywkowski and Penna, 2020). Briefly, all frames in a selected calcium imaging recording were scored as either positive of negative for a particular behavior and a ROC curve was generated for each neuron by plotting the true positive and false positive rates across the distribution. The area under the ROC curve (auROC) was calculated for each recording and averaged across days to give a measure of the responsiveness of a neuron to a given behavior. Neurons with auROC values greater than 0.65 or less than 0.35 were considered to respond to the behavior. Significance of neuronal response auROC values was estimated by calculating the mean and standard deviation (SD) of the distribution of auROC values obtained by shuffling true and false positive labels 1000 times. Neuronal response auROC values ≥ 3 SD away from the mean were considered significant. Neuron needed to fulfill both of the above criteria to be classified as positive or negative responding during any given behavior (e.g. Assessment+, Defeat-). For the analysis of auROC values for flight behavior 2 s before and after the initiation of each flight were labeled as true negative and true positive, respectively. Only recordings with at least two instances of relevant behavior were selected for auROC analysis. For example, during pre-defeat context recordings only four animals out of seven transitioned between all the compartments of the apparatus two or more times and the remaining three animals were excluded from analysis. Heat maps were generated by extracting calcium imaging data (meaned z-score over all trials) for individual neurons from before and after the onset or offset of a given behavior. Significant change across behavioral onset or offset was estimated for each neuron by performing a Wilcoxon signed rank-sum test (Python, scikit-learn) on the distribution of values that resulted from averaging the z-scores of frames from the before and after periods for each trial. For linear discriminant analysis (LDA algorithm, Python, scikit-learn), data was normalized within days by calculating z-scores for each recording separately. Each frame was then expressed as a vector containing calcium values for all neurons. The distance between clusters was quantified by calculating the average distance between data points of a given cluster and all other data points from other clusters. For Venn diagrams, we included only neurons responding to the behaviors (e.g. Defeat+) that were recorded across all relevant behaviors, resulting in different numbers of neurons across Venn analyses (Supplementary file 1). Probability of chance overlap was computed by multiplying the probabilities of responsive neurons. Significance was assessed using Fisher’s test.

Prism Graphpad five software or custom scripts in Python were used to generate graphs and perform statistical analysis. For calcium imaging Wilcoxon rank-sum test (Python, scikit-learn package) was performed. For cFos experiments, one-way ANOVA with Tukey’s post-hoc test was used. For optogenetic experiments repeated measures ANOVA with Tukey’s post-hoc test was used. For assessing significance of overlap between different populations of neurons Fisher’s test was used. For assessing significance differences in LDA representations the one-way ANOVA and T test were used.

Custom code written for this study is made available on GitHub platform https://github.com/GrossLab/VMHvl-Calcium-Imaging-and-Visualisation (copy archived at https://github.com/elifesciences-publications/VMHvl-Calcium-Imaging-and-Visualisation). Behavioral and imaging data will be made available upon reasonable request.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We thank Dominic Evans and Tiago Branco for sharing expertise in calcium endoscopy and Daniel Rossier, Irene Ayuso, Senthil Deivasigamani, Hiroki Asari and Santiago Rompani for helpful comments on the manuscript. The work was supported by EMBL and the European Research Council (ERC) Advanced Grant COREFEAR to CTG.

Animal experimentation: All experimental procedures involving the use of animals were carried out in accordance with EU Directive 2010/63/EU and under approval of the EMBL Animal Use Committee and Italian Ministry of Health License 541/2015-PR to CG.

1. Received: March 23, 2020
2. Accepted: August 12, 2020
3. Version of Record published: September 21, 2020 (version 1)

© 2020, Krzywkowski 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.