Skip to main content
eLife logoLink to eLife
. 2024 Feb 16;13:e86581. doi: 10.7554/eLife.86581

Functionally refined encoding of threat memory by distinct populations of basal forebrain cholinergic projection neurons

Prithviraj Rajebhosale 1,, Mala R Ananth 1,, Ronald Kim 1, Richard Crouse 2,§, Li Jiang 1, Gretchen López-Hernández 3, Chongbo Zhong 1, Christian Arty 4,#, Shaohua Wang 5, Alice Jone 6,, Niraj S Desai 1, Yulong Li 7,8,9, Marina R Picciotto 2,10, Lorna W Role 1,‡,, David A Talmage 1,‡,
Editors: Joshua Johansen11, Laura L Colgin12
PMCID: PMC10928508  PMID: 38363713

Abstract

Neurons of the basal forebrain nucleus basalis and posterior substantia innominata (NBM/SIp) comprise the major source of cholinergic input to the basolateral amygdala (BLA). Using a genetically encoded acetylcholine (ACh) sensor in mice, we demonstrate that BLA-projecting cholinergic neurons can ‘learn’ the association between a naive tone and a foot shock (training) and release ACh in the BLA in response to the conditioned tone 24 hr later (recall). In the NBM/SIp cholinergic neurons express the immediate early gene, Fos following both training and memory recall. Cholinergic neurons that express Fos following memory recall display increased intrinsic excitability. Chemogenetic silencing of these learning-activated cholinergic neurons prevents expression of the defensive behavior to the tone. In contrast, we show that NBM/SIp cholinergic neurons are not activated by an innately threatening stimulus (predator odor). Instead, VP/SIa cholinergic neurons are activated and contribute to defensive behaviors in response to predator odor, an innately threatening stimulus. Taken together, we find that distinct populations of cholinergic neurons are recruited to signal distinct aversive stimuli, demonstrating functionally refined organization of specific types of memory within the cholinergic basal forebrain of mice.

Research organism: Mouse

Introduction

Acetylcholine (ACh) is critical for cognition. Basal forebrain cholinergic neurons (BFCNs), neurons that synthesize and release ACh that are sparsely distributed throughout the base of the forebrain, provide the primary source of ACh to the cortex, hippocampus, and amygdala. Disruptions to normal cholinergic transmission are thought to contribute to several neuropsychiatric disorders (Sarter et al., 1999; Higley and Picciotto, 2014) as well as to cognition (Ananth et al., 2023) and salience-related behaviors (Jiang et al., 2016; Hersman et al., 2017; Crouse et al., 2020). BFCNs are anatomically divided into several clusters: the medial septum/diagonal band complex (MS/DB), the ventral pallidum (VP), the substantia innominata (SI), and the nucleus basalis (NBM). Between and within these anatomical groupings, BFCNs comprise heterogenous subclusters (Zaborszky and Gyengesi, 2012). How this heterogeneity contributes to the significant control that cholinergic signaling exerts over large, behaviorally relevant circuits is unclear (Zaborszky et al., 2015; Gielow and Zaborszky, 2017).

ACh plays an important role in modulating emotionally salient memories (Luchicchi et al., 2014; Ballinger et al., 2016; Knox, 2016; Ananth et al., 2023). We and others have found that cholinergic signaling in the basolateral amygdala (BLA) is important for generating defensive behaviors in response to both learned and innate threats (Power and McGaugh, 2002; Jiang et al., 2016; Wilson and Fadel, 2017). Optogenetic manipulation of endogenous ACh release in the BLA during learning modulates the expression of threat response behaviors in mice upon recall of a conditioned stimulus (Jiang et al., 2016). Stimulating release of ACh increases activity of BLA principal neurons, in part by increasing the release probability of glutamatergic inputs to these neurons, and is sufficient to induce long-term potentiation (LTP) when paired with minimal (non-LTP generating) stimulation of glutamatergic input to the BLA (Unal et al., 2015; Jiang et al., 2016). Memory formation and retrieval are associated with fast synaptic mechanisms that are modulated by ACh, that are in turn necessary for the proper learning and expression of threat response behaviors (Nonaka et al., 2014). Given the broad distribution of cholinergic input across the BLA, and the well-established role of ACh in modulating BLA plasticity, the basal forebrain cholinergic system is well positioned to serve an important role in the encoding of threat memories and generation of threat response behaviors (Ananth et al., 2023).

The BLA receives dense cholinergic input from neurons located in various regions within the basal forebrain (such as the VP, SI, and NBM). In this study, we asked how these distinct populations of BLA-projecting BFCNs contribute to threat responses. Using a genetically encoded ACh sensor, activity-dependent genetic tagging, chemogenetic manipulations, and electrophysiological recordings, we identify a population of BFCNs in the NBM/SIp (SIp defined as the portion of the sublenticular SI posterior to bregma −0.4 mm) that are required for learned threat responsiveness. We find that NBM/SIp cholinergic neurons are necessary for freezing behavior following cue-conditioned threat learning while freezing behavior elicited by an innately threatening stimulus activates cholinergic neurons in the VP/SIa (VP/SIa; SIa defined as the portion of the SI ventral to the anterior commissure located anterior to bregma −0.4 mm).

Results

Animals recognize varied sensory stimuli and categorize them as either threatening or non-threatening. Recognition of threatening stimuli can be innate or acquired, for example, by association of an aversive experience with an innocuous, co-occurring sensory input. In this study, we sought to understand if the basal forebrain cholinergic system participates in the encoding of associative threat or in response to innate threat.

ACh is released in the basal lateral amygdala in response to threat

The BLA plays a central role in associative threat learning and in the generation of threat responses. We have previously demonstrated that silencing cholinergic input to the BLA during cue-conditioned threat learning (pairing a naive tone with a foot shock) blunts learned freezing in response to the conditioned stimulus (tone) (Jiang et al., 2016). Given this, the first question we asked was whether ACh was released in the BLA during associative threat learning (Figure 1, Figure 1—figure supplements 14). To monitor acute changes in extracellular ACh levels during the cue conditioned threat-learning task, we expressed a genetically encoded ACh sensor, GRABACh3.0 (Jing et al., 2018; Jing et al., 2020) in BLA neurons and visualized fluorescence using fiber photometry (Figure 1A). Our associative threat-learning protocol involved placing mice in a novel chamber and exposing them to an 80 dB tone for 30 s. During the final 2 s of the tone the mice received a foot shock (0.7 mA). The tone–shock pairing was repeated twice (for a total of three pairings). Twenty-four hours later, mice were placed in a different chamber (with different tactile, visual, and olfactory cues to the training chamber) and exposed to tone alone.

Figure 1. Acetylcholine (ACh) is released in basolateral amygdala (BLA) during threat recall (see also Figure 1—figure supplements 14).

(A) Left: Schematic of strategy for targeting the genetically encoded ACh sensor (GRABACh3.0) to BLA. Right: Image of ACh sensor expression (green). White dotted oval delineates ACh sensor-labeled BLA. White dotted box denotes prior location of optical fiber. Scale bar = 100 µm. Please refer to Figure 1—figure supplement 4 for fiber placement mapping in all mice. (B) Schematic of the associative threat-learning protocol employed consisting of three tone + shock pairings during the training period and tone alone during the recall session. (C) Average traces of ACh release in response to tone; shading represents standard error of the mean (SEM): naive tone (Tone 1 during training) in black, recall tone in deep blue (tone onset indicated by bar above; n = 11). BL = baseline. (D) Quantification of ACh release (area under the curve (AUC)) during baseline period (1 s prior to tone onset) and in response to the first (naive) tone and in response to the recall tone (1 s following tone onset). Naive tone did not induce significant increase in ACh release in the BLA (Wilcoxon matched-pairs signed rank test, BL vs. Tone 1, p = 0.8311, W = −6). Recall tone induced a significant increase in ACh release in the BLA (paired t-test, BL vs. Recall tone, ** p = 0.0039 (two-tailed), t = 3.732, df = 10). See also Figure 1—figure supplements 14.

Figure 1—source data 1. Processed ACh release fiber photometry data for Figure 1.

Figure 1.

Figure 1—figure supplement 1. Exposure to tones during tone–shock pairings does not significantly alter acetylcholine (ACh) release in the basolateral amygdala (BLA) during training.

Figure 1—figure supplement 1.

(A) Left: Schematic of the associative threat-learning protocol employed consisting of three tone + shock pairings during the training period. Right: Quantification of typical freezing behavior throughout the training session (30 s time bins; n = 9). Shaded blue bars indicate periods of tone (Inline graphic)–shocks (Inline graphic) pairings. (B) Sample traces of raw GRABACh3.0 fluorescence in the BLA (green) and reference channel (black). (C) Average traces of ACh release for the three tone + shock pairings (30 s) with 5 s prior to tone onset and 5 s after tone (n = 11); shading represents standard error of the mean (SEM) (Tone 1: black, Tone 2: deep sky blue, Tone 3: navy blue). (D) ACh release in response to tone (Inline graphic) during training (n = 11). Average traces of ACh release in response to tone; shading represents SEM. Quantification of ACh release (area under the curve (AUC)) during BL (1 s prior to onset of Tone 1) and during Tones 1, 2, and 3 (1 s after tone onset) (Tone 1: black, Tone 2: deep sky blue, Tone 3: navy blue) for each tone + shock pairing. There were no significant changes in ACh release in response to tone during training (Friedman test, p = 0.1381).
Figure 1—figure supplement 2. Repeat tone exposures without foot shock during either a training session or a recall session fail to induce significant changes in acetylcholine (ACh) release in the basolateral amygdala (BLA).

Figure 1—figure supplement 2.

(A) Schematic of tone only training protocol consisting of three 30 s tones on day 1 without shocks followed by presentation of the same tone 24 hr later, on day 2. (B) Average traces of ACh release in response to tone; shading represents standard error of the mean (SEM): naive tone (Tone 1 on day 1) in black, 24 hr tone in deep blue (tone onset indicated by bar above; n = 5). (C) Quantification of ACh release (area under the curve (AUC)) during baseline period (1 s prior to tone onset) and in response to the first (naive) tone and in response to the 24 hr tone (1 s following tone onset). Insets display traces from individual mice in response to the respective tones; bar above trace indicates tone onset. Neither the naive tone nor the repeat tone presentation (24 hr tone) significantly increased ACh release relative to baseline (paired t-test. BL vs. Tone 1, p = 0.8437, t = 0.2103, df = 4. BL vs. 24 hr tone, p = 0.3152, t = 1.147, df = 4). (D) Average traces of ACh release for the three tones on day 1 (30 s) with 5 s prior to tone onset (n = 5). Shading represents SEM (Tone 1: black, Tone 2: deep sky blue, Tone 3: navy blue). Red inset showing 1 s before and after onset of the tones is expanded on the right.
Figure 1—figure supplement 2—source data 1. Processed ACh release fiber photometry data for Figure 1—figure supplement 2.
Figure 1—figure supplement 3. Shock alone does not sensitize cholinergic tone responses in the basolateral amygdala (BLA).

Figure 1—figure supplement 3.

(A) Top: Schematic of shock only training and recall protocol consisting of three 2 s shocks during the training session (day 1) and a single 30-s tone during the recall session (day 2). Acetylcholine (ACh) release in response to shocks (Inline graphic) during training and to tone alone (Inline graphic) 24 hr following exposure to three-foot shocks (shock without tone) (n = 5). (B) Left: Average trace of ACh release in response to shocks (Inline graphic) shows robust ACh release in response to each shock. Right: Average trace of ACh release in response to tone (Inline graphic). Shading represents standard error of the mean (SEM). (C) Freezing behavior during the ‘recall’ (tone alone 24 hr after shock only training) session. Blue shaded bar represents the ‘24 hr tone’. There was no significant increase in freezing in response to the tone when training was conducted with shock alone (p = 0.2418, repeated measures (RM) one-way analysis of variance (ANOVA), F(2.381, 9.523) = 1.657). (D) Quantification of ACh release (area under the curve (AUC)) during the 24 hr tone session during the 1 s following 24 hr tone onset compared to baseline (1 s prior to tone onset). There was no change in ACh release in response to the tone when training with shock alone (no tone–shock pairing) (AUC, p = 0.7272 (two-tailed), t = 0.3742, df = 4, paired t-test).
Figure 1—figure supplement 3—source data 1. Processed ACh release data for Figure 1—figure supplement 3.
Figure 1—figure supplement 4. Remapping of optic fiber placements.

Figure 1—figure supplement 4.

(A) Relocalization of fiber optic placement mapped across the anterior–posterior axis of the basolateral amygdala (BLA). Top: Atlas images showing fiber tip locations (red squares). (B) Left: Average traces of acetylcholine (ACh) release for the 2 s prior to tone onset through 28 s of recall tone during the naive tone (Tone 1) presentation (black line) from mice that underwent the training (tone + shock) + recall (tone alone) paradigm. Right: Average traces of ACh release for 2 s prior to tone onset through 28 s of the recall tone (deep blue line) from mice that underwent the training (tone + shock) + recall (tone alone) paradigm; shading represents standard error of the mean (SEM; n = 11). Tone onset indicated by black bar above trace. (C) Expanded traces of ACh release from individual mice for the naive tone (Tone 1, left) and recall tone (right) over 2 s including 1 s prior to tone onset and 1 s after tone onset from mice that underwent the training (tone + shock) followed by recall (1× tone) paradigm (n = 11) (Tone 1, left: black lines; Recall tone, right: deep blue lines).

Foot shock, either alone or paired with tone, increased ACh release in the BLA whereas the naive tone that is the first tone before shock presentation (Tone 1), did not (Figure 1—figure supplement 1C, Figure 1C, D, left; baseline (BL) vs. Tone 1, p = 0.8311). In contrast to Tone 1, the recall tone, presented 24 hr after the three tone–shock pairings, resulted in significant increase in ACh release in the BLA (Figure 1C, D, right; p = 0.0039). The change in tone-associated ACh release required pairing with foot shock: naive tone (Figure 1—figure supplement 2C left, p = 0.8437), three consecutive tones alone (without shock), or a subsequent repeat tone presentation after 24 hr (not previously paired with shock) (Figure 1—figure supplement 2C, right; p = 0.3152), did not induce significant changes in ACh release in the BLA (Figure 1—figure supplement 2).

To verify that the increases in ACh release were indeed specific to the tone–shock association and not due to generalization from prior shock exposure, we also subjected mice to three shocks (day 1) followed by a tone presentation 24 hr later (day 2) (Figure 1—figure supplement 3A). While mice demonstrated freezing behavior during the session on day 2, there was no significant increase in freezing behavior to the 24 hr tone presentation (Figure 1—figure supplement 3C, p = 0.2418). There was no increase in ACh in response to the tone when it was not explicitly paired with a shock, confirming that the changes in ACh release were indeed associative (Figure 1—figure supplement 3D; baseline (pre-tone, day 2)–24 hr tone (tone presentation, day 2): p = 0.7272). Therefore, after repeated tone–shock pairings, BLA-projecting cholinergic neurons acquire enhanced tone responsiveness.

NBM/SIp cholinergic neurons are activated by threat learning and reactivated during threat memory recall

Following associative threat-learning, cholinergic neurons exhibited increased ACh release in the BLA in response to a previously innocuous auditory stimulus; this increase occurred exclusively following pairing of the tone with a shock. Using a two-color labeling system, we asked whether NBM/SIp cholinergic neurons were activated during the training session and reactivated during the recall session. To do this, we injected the offspring of a cross of Chat-IRES-Cre x Fos-tTA:Fos-shGFP with a viral vector, AAV9-TRE-DIO-mCherry-P2A-tTAH100Y, resulting in activity (tTA) dependent, Cre dependent (aka ADCD) mCherry expression (see methods and Figure 2—figure supplement 1). These mice carry three transgenes: one encoding Cre recombinase in cholinergic neurons, a second doxycycline (Dox) repressible, tetracycline transactivator (tTA) expressed following activation of the fos promoter, and a third destabilized green fluorescent protein (short half-life GFP) also under transcriptional regulation of the fos promoter. tTA and shGFP are transiently expressed in activated neurons. In the absence of Dox (delivered via chow diet), activation of Cre-expressing cholinergic neurons leads to tTA expression and expression of the virally transduced mCherry along with a mutant tTA, which is insensitive to Dox. Thus, after closure of the labeling window by re-administration of Dox, cholinergic neurons activated during the Dox off period maintain mCherry expression permanently driven by the mutant tTA. When ADCD labeling is coupled with the transient expression of Fos−shGFP, we can label and visualize participation of cholinergic neurons in two separate behavioral sessions (mCherry+ = session 1 activated cells and GFP+ = session 2 activated cells) (Figure 2—figure supplement 1B).

Two to three weeks following injection with the ADCD virus, mice were either (1) kept in home cage throughout, (2) exposed to tone without foot shock (tone alone), or (3) put through the standard threat-learning paradigm (tone + shock). Twenty-four hours prior to the training session (session 1) mice were switched from Dox-containing to Dox-free chow to allow function of tTA. Immediately following tone–shock pairings, mice were placed back on Dox-containing chow (Figure 2A). This switch from Dox on → Dox off → Dox on was also performed for mice that remained in their home cages and for those that were exposed to tones without shock. Recall was performed 72 hr later (tone alone in a new context), and mice were sacrificed ~2.5 hr following recall (the peak of the Fos−shGFP expression). We quantified the number of mCherry+/GFP+ (double positive) neurons following session 2 (e.g. white arrow, Figure 2B). Significantly more double positive cholinergic neurons were seen following the complete associative threat-learning paradigm (tone + shock followed by tone recall) compared to mice that underwent session 1 without shocks (Figure 2C, p = 0.0249). To further ensure that the reactivation of these cholinergic neurons was not due to a generalized increase in responsiveness of these neurons following shock exposure, we quantified reactivated neurons in mice exposed to shock alone during session 1 followed by tone alone during session 2 (shock alone (session 1) → tone alone (session 2)) along with shock alone (session 1) → home cage (session 2), and home cage controls (Figure 2D). All three conditions showed few reactivated neurons and no differences between groups (p = 0.9471). Thus, associative threat-learning results in activation of NBM/SIp cholinergic neurons which are reactivated during subsequent cue-induced memory recall.

Figure 2. Nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic neurons are activated by threat learning and reactivated during threat memory recall (see Figure 2—figure supplement 1).

(A) Strategy for labeling activated NBM/SIp cholinergic neurons during both training and recall. Chat-IRES-Cre × Fos-tTA/shGFP mice (n = 14) were injected in the NBM/SIp with ADCD-mCherry virus (AAV9: TRE-DIO-mCherry-P2A-tTAH100Y). During session 1 (off Dox) mice either remained in their home cage, were exposed to three tones (Tone alone), or were exposed to three tone–shock pairings (training session). During session 2, mice remained in home cage or were exposed to a single tone (recall session). Cholinergic neurons activated during training express ADCD-mCherry stably after training (red during training), and neurons activated during recall transiently express GFP (green during recall). (B) Image of the NBM/SIp showing cholinergic neurons activated during training (red arrow) or by both training and recall (reactivated – white arrow; image taken at A/P ~ −0.8 from Bregma; scale bar = 50 µm). (C) Quantification of the number of cholinergic neurons activated during session 1 (ADCD-mCherry+) that were reactivated during session 2 (both mCherry and GFP positive (activated both during session 1 (training) and during session 2 (recall))). Home cage (n = 7 sections from 4 mice), tone only (n = 4 sections from 2 mice) and tone + shock (n = 17 sections from 8 mice) conditions. Significantly more cholinergic neurons were reactivated by tone following tone–shock pairings (Kruskal–Wallis p = 0.0249). Tone–shock compared to tone only (*p = 0.0464, Dunn’s corrected). (D) Quantification of number of reactivated cholinergic neurons (activated both during session 1 and during session 2 vs. the total number of cholinergic neurons activated during session 1). Home cage → home cage (n = 9 sections from 3 mice), shock → home cage (n = 17 sections from 4 mice), and shock → tone (n = 11 sections from 3 mice). Kruskal–Wallis p = 0.9471, KW = 0.1219.

© 2024, BioRender Inc

Figure 2 was created using BioRender, and is published under a CC BY-NC-ND license. Further reproductions must adhere to the terms of this license.

Figure 2.

Figure 2—figure supplement 1. Construction of viral vectors for tagging activated cholinergic neurons.

Figure 2—figure supplement 1.

(A) Top: Diagrams of 1. ADCD and 2. ADCD hM4Di viral constructs. Activity dependence is conferred by the Tet response element (TRE – 7 repeats of the tetO sequence followed by a minimal promoter). Cre dependence is conferred by pairs of loxP and lox2272 sites flanking the ‘cargo’ in an antisense orientation. Cargos: permanent labeling ADCD includes an oChIEF-mCherry fusion followed by a P2A element and a doxycycline insensitive tTA (tTA*). ADCD-hM4Di cargo includes an hM4Di-mCherry fusion protein without the dox-insensitive tTA. Bottom: Images showing lack of ADCD expression in Chat-Cre and Fos-tTA mice, and ADCD expression in Chat-IRES-Cre × Fos-tTA/shGFP mice. Scale bar = 50 µm. Inset shows magnified image of two cells expressing ADCD-mCherry (white arrows). Scale bar = 10 µm. (B) Schematic representing transient expression of the shGFP under control of the Fos promoter. Activated neurons express shGFP (green line), which peaks in expression ~2.5 hr following activation. This GFP expression is not influenced by the presence or absence of doxycycline and is transient. The virally delivered ADCD construct can only be expressed in (cre expressing) cholinergic neurons, in the absence of doxycycline (DOX OFF). Once expression initiates, it is maintained even in the presence of Dox by a dox-insensitive tTA (red line), thereby resulting in permanent labeling. However, labeling in cells that were not activated during the DOX OFF period is prevented. Intercepting the peak of a shGFP signal (for example following a second behavioral assay) can allow for identification of cells that express mCherry (activated during DOX OFF) and cells that express GFP (activated 2.5 hr prior to euthanasia). (C) Test to determine the minimal time off dox to allow ADCD expression. Mice were kept on dox food for 2 weeks (starting 2 days prior to ADCD +AAV9-camk2a-Cre-IRES-GFP injection into the basolateral amygdala (BLA)) (Dox ON). They were then shifted to regular chow either 6, 12, or 24 hr prior to training (Dox OFF). They were returned to dox-chow immediately following training (ON). ADCD expression was quantified 72 hr later as % of mCherry+ cells out of total GFP/Cre+ cells in the BLA. (Note the background level of expression in the presence of Dox. See Methods for details.)

Reactivation of cholinergic neurons activated by training is required for learned behavioral responses

BLA-projecting cholinergic neurons acquire tone responsiveness following associative threat learning (Figure 1) and a population of NBM/SIp cholinergic neurons are activated during tone–shock pairing and reactivated during the recall session (Figure 2). If these cholinergic neurons are indeed part of a threat memory engram, then their reactivation would be required for generation of learned threat responses. To block reactivation of cholinergic neurons in response to tone, we expressed the inhibitory, designer receptor hM4Di, in an activity dependent, Cre-dependent manner in NBM/SIp cholinergic neurons (ADCD-hM4Di; Figure 3A) and subjected these mice to the threat-learning paradigm (Figure 3A). Mice were taken off Dox-chow 24 hr prior to the training session, immediately placed back on Dox-chow after training, and then tested for tone recall after 72 hr. ADCD-hM4Di and sham operated control mice were injected with clozapine (CLZ; 0.1 mg/kg; injected intraperitoneally (i.p.)) 10 min prior to the recall session to selectively silence the population of NBM/SIp cholinergic neurons that were previously activated during training (Figure 3A). Freezing behavior was quantified during both the training and recall sessions. Freezing was compared between the ‘Pre-Tone’ period and ‘Recall Tone Response’ (defined as freezing occurring from the onset of the recall tone through the end of the recall session) (Figure 3—figure supplement 1B). Both groups of mice showed the same freezing behavior during the training session (Figure 3C, p = 0.6482, Figure 3—figure supplement 1A). In the recall session, sham mice displayed typical freezing behavior in response to tone (Figure 3D, gray boxes; Pre-Tone vs. Recall Tone Response, p = 0.0001). In contrast, ADCD-hM4Di mice did not show increased freezing in response to the tone (Figure 3D, red boxes; BL vs. tone response, p = 0.8451). Overall ADCD-hM4Di mice showed lower freezing behavior compared to sham controls (Figure 3D; sham – gray, hM4Di – red: p = 0.0052), indicating that reactivation of training-activated NBM/SIp cholinergic neurons during the recall session was required for the expression of learned threat response behavior.

Figure 3. Re-activation of a subset of nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic neurons is required for threat memory retrieval (see Figure 3—figure supplement 1).

(A) ADCD-hM4Di (AAV9: TRE-DIO-hM4Di.mCherry) was injected into the NBM/SIp of Chat-IRES-Cre × Fos-tTA/shGFP mice. Two weeks later mice underwent training on regular chow (Dox-chow removed 24 hr prior to training session) to allow hM4Di.mCherry to be selectively expressed in training-activated cholinergic neurons. Three days later, recall was tested in Dox on conditions. Clozapine (CLZ) was injected 10 min before the recall session to activate the inhibitory DREADD, hM4Di specifically in previously activated cholinergic neurons. (B) Representative image taken at Bregma −0.8 mm of mCherry (ADCD-hM4Di.mCherry) expressing cells. Inset shows higher magnification images of ADCD expression (red arrows). Scale bar = 50 µm. (C) Freezing behavior during training in sham (gray, n = 8 mice) and ADCD-hM4Di injected (red, n = 7 mice) for each 30 s bin during tone presentation (Tone + Shock 1, 2, 3). There were no significant differences between the groups during the training session (repeated measures (RM) two-way analysis of variance (ANOVA) Time × Group p = 0.6482; Group p = 0.7311). (D) Freezing behavior during recall following selective hM4Di-mediated inhibition of training-activated cholinergic neurons in the NBM/SIp. Sham (gray, n = 8 mice) and hM4Di (red, n = 7 mice) groups. There were significant differences between pre-tone vs. tone-related freezing for sham (****Pre-Tone vs. Recall Tone Response, p = 0.0001, Bonferroni corrected), response to tone between sham and hM4Di (**p = 0.0026, Bonferroni corrected) and a significant main effect of Time × Group interaction (RM two-way ANOVA (GLM) Time × Group, p = 0.0052). (See Figure 3—figure supplement 1 for details on time periods comprising Pre-Tone and Recall Tone Response periods.)

© 2024, BioRender Inc

Figure 3 was created using BioRender, and is published under a CC BY-NC-ND license. Further reproductions must adhere to the terms of this license.

Figure 3.

Figure 3—figure supplement 1. Time-resolved freezing plot from recall following silencing of training-activated NBM/SI basal forebrain cholinergic neurons (BFCNs).

Figure 3—figure supplement 1.

(A) Freezing behavior during training in sham (black, n = 8 mice) and ADCD-hM4Di injected (red, n = 7 mice). Error bars represent standard error of the mean (SEM). Tone 1 (Inline graphic), Tone 2 (Inline graphic), Tone 3 (Inline graphic), and shock (Inline graphic) presentation denoted with symbols above traces. (B) Freezing behavior during recall in sham (black, n = 6 mice) and ADCD-hM4Di injected (red, n = 5 mice). Error bars represent SEM. Recall tone (Inline graphic) denoted by symbol above traces.

BLA-projecting NBM/SIp cholinergic neurons are reactivated during threat memory recall

To investigate whether NBM/SIp cholinergic neurons that are reactivated during recall are BLA projecting, we injected Chat-IRES-Cre × Fos-tTA:Fos-shGFP mice with ADCD-mCherry in the NBM/SIp, and simultaneously delivered the retrograde tracer Fast Blue (FB) into the BLA (Figure 4A). The mice were taken off doxycycline containing chow during the training period, returned to dox-chow for 72 hr and then exposed to the tone alone. We then quantified BLA-projecting cholinergic neurons that were reactivated by tone (ChAT immunoreactive, FB labeled and ADCD-mCherry+/Fos-shGFP+; Figure 4C). We found that ~20% of NBM/SIp cholinergic neurons in both the home cage and the threat-learning + recall paradigm group (at Bregma −0.8 mm) were labeled with FB, with no significant differences in the percentage of cholinergic neurons with retrograde label between groups (Figure 4D; p = 0.5192). Next, we quantified the percentage of BLA-projecting NBM/SIp cholinergic neurons that were active during session 1 and reactivated during session 2. We found that, on average, ~21% of BLA-projecting cholinergic neurons were reactivated during recall (Figure 4E). This reactivation of BLA-projecting BFCNs was significantly higher in mice that underwent training + recall compared to mice that remained in their home cage but still underwent the Dox on → Dox off → Dox on protocol (Figure 4E; p = 0.0183). Based on these data we conclude that BLA-projecting BFCNs are activated by associative threat learning and reactivated by threat recall.

Figure 4. Basolateral amygdala (BLA)-projecting nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic neurons are reactivated by the conditioned tone stimulus.

Figure 4.

(A) Left: Strategy for labeling activated NBM/SIp cholinergic during both training and recall along with mapping of BLA-projecting neurons. Chat-IRES-Cre × Fos-tTA/shGFP mice (n = 7) were injected in the NBM/SIp with ADCD-mCherry virus and in the BLA with Fast Blue dye. During session 1 (off Dox) mice either remained in their home cage or were exposed to three tone–shock pairings. During session 2 (recall session), mice remained in home cage or were exposed to a single tone. Cholinergic neurons activated during training express GFP transiently and express ADCD-mCherry stably after training (red during training), and neurons activated during recall transiently express GFP (green during recall). Neurons projecting to the BLA were labeled by Fast Blue (blue). Cholinergic neurons were identified by ChAT staining (magenta). (B) Image of a Fast Blue injection site in the BLA; Inset: Mapping of injection sites for all Fast Blue experiments. (C) Representative image showing (clockwise), ChAT+ neurons in the NBM/SIp at bregma −0.8 mm (magenta), BLA-projecting neurons (blue, blue arrowheads), training-activated cells (ADCD) (red), and recall-activated neurons (green, green arrow heads). BLA-projecting basal forebrain cholinergic neurons (BFCNs) activated by training and recall are denoted by double arrowheads (blue and green). Scale bar = 50 µm. Inset scale bar = 10 µm. (D) Quantification of percentage of ChAT+ neurons that were labeled by Fast Blue in mice from the home cage group (gray) and mice from the training + recall group (black) from bregma −0.8 mm. No significant differences were found between groups (Welch’s t-test, p = 0.5192). (E) Quantification of percentage of BLA-projecting BFCNs (ChAT+/Fast Blue+) at bregma −0.8 mm that were reactivated during session 2 (ADCD + GFP) in mice from the home cage group (n = 2) (gray) and mice from the training + recall group (n = 5) (black). Mice that underwent training and recall had significantly higher numbers of engram-enrolled BLA-projecting BFCNs (Welch’s t-test, *p = 0.0183).

© 2024, BioRender Inc

Figure 4 was created using BioRender, and is published under a CC BY-NC-NDlicense. Further reproductions must adhere to the terms of this license.

Silencing BLA-projecting BFCNs during training or recall prevents activation of BLA neurons and conditioned freezing behavior

To determine whether chemogenetic silencing of BLA-projecting cholinergic neurons during training or during recall interfered with the activation of BLA neurons, we injected the BLA of Chat-IRES-Cre mice with CAV2-DIO-hM4Di.mCherry and AAV9-camk2a-GCaMP (cav.hM4DiBLA mice) or AAV9-camk2a-GCaMP alone (sham mice) (Figure 5A, Figure 5—figure supplement 1A; GFP fluorescence from GCaMP was used to mark the injection sites). We found mCherry was expressed in cholinergic neurons predominantly in the NBM/SIp, followed by the VP/SIa, with a small contribution from the horizontal limb of the diagonal band of Broca (hDB) (Figure 5A, right). These data support previous findings (Zaborszky and Gyengesi, 2012) that NBM/SIp cholinergic neurons provide a major input to the BLA.

Figure 5. Basolateral amygdala (BLA)-projecting cholinergic neuronal activity is required both during training and during recall for learned threat processing (see Figure 5—figure supplements 1 and 2).

(A) Left: Strategy for retrograde targeting of hM4Di DREADD to BLA-projecting cholinergic neurons. Middle: Re-localization of BLA injection sites (using AAV9-camk2a-GCaMP6f to mark the injection site), and identification of retrogradely labeled cholinergic neurons within the nucleus basalis and posterior substantia innominata (NBM/SIp; scale bar = 50 µm). Right: Quantification of hM4Di-expressing cholinergic populations (mCherry+) across the basal forebrain (n = 8 mice, 56–80 sections) (Bregma +0.6 mm to −1.5 mm). (B) BLA-projecting cholinergic neurons were silenced by injecting mice with clozapine (CLZ) 10 min prior to training. Percent time freezing during the recall session including the pre-tone (baseline) period and in response to the conditioned tone. CLZ was only administered during the training session (repeated measures (RM) two-way analysis of variance (ANOVA), Time × Group p = 0.0047; Group p = 0.0007). Sham vs. DREADD (Tone response, ****p < 0.0001, Bonferroni corrected). (b’) DREADD-induced silencing of BLA-projecting cholinergic neurons during training reduced BLA Fos immunoreactivity following recall. Representative BLA images from sham injected and CAV2-DIO-hM4Di mice fixed and stained with anti-Fos antibodies (white) at 45–60 min following recall. Dotted line outlines the BLA (scale bar = 100 µm). (b’’) The density of recall-activated BLA neurons under sham injected conditions vs. following selective inhibition of the BLA-projecting cholinergic neurons using CAV2-hM4Di (Fos+). Fos+ cell density in BLA sham injected (black) vs. CAV2-DIO-hM4Di.mcherry (red) (n = 7 mice/group, averaged from 22 sections sham vs. 28 sections hM4Di). Mann–Whitney test: *p = 0.0286. Lines represent median for each group. (C) BLA-projecting cholinergic neurons were silenced during recall (clozapine given ONLY 10 min prior to the recall). Freezing differed significantly between pre-tone vs. recall tone response for sham and DREADD groups (RM two-way analysis of variance (ANOVA), pre-tone vs. recall tone response, sham ****p < 0.0001; DREADD ***p = 0.0003). There was a significant effect of group, p = 0.0312. sham and DREADD groups were significantly different in their response to the recall tone, p = 0.0279. All multiple comparisons were Bonferroni corrected. (c’) hM4Di-induced silencing of BLA-projecting cholinergic neurons during recall reduced BLA Fos immunoreactivity following recall. BLA images following Fos immunostaining (scale bar = 100 µm). (c’’) Fos+ cell density in BLA between sham injected (black) vs. CAV2-DIO-hM4Di.mcherry (red) injected mice (n = 4–5 mice/group, averaged from 30 sections sham vs. 32 sections hM4Di). Mann–Whitney test: *p = 0.0317.

Figure 5.

Figure 5—figure supplement 1. DREADD-induced silencing of basolateral amygdala (BLA)-projecting cholinergic neurons reduces threat induced activation of anterior BLA and CeC neurons.

Figure 5—figure supplement 1.

(A) Top: Strategy for retrograde targeting of hM4Di DREADD to BLA-projecting cholinergic neurons. BLA-projecting cholinergic neurons were silenced by injecting mice with clozapine (CLZ) prior to training. Recall was conducted 24 hr later and mice were euthanized 45–60 min after recall and probed for Fos immunoreactivity in the BLA. Bottom: Representative images from a mouse showing expression of GFP (AAV9-camk2a-GCaMP6f co-injected with the CAV2-DIO-hM4Di.mCherry virus) across the A–P extent of the BLA. Scale bar = 50 µm. (B) Representative images from home cage control, sham and cav.hM4DiBLA mice at bregma locations from −0.8 to −1.4 mm showing Fos immunoreactivity (grayscale). Silencing BLA-projecting basal forebrain cholinergic neurons (BFCNs) during training reduced Fos expression in the BLA during recall. BLA is outlined by white dashed lines. Scale bar = 50 µm. (C) Silencing BLA-projecting cholinergic neurons during Recall reduced Fos immunoreactivity in the capsular portion of the central amygdala (CeAc/CeC) when profiled following recall. Images of CeC Fos immunoreactivity in sham vs. CAV2-DIO-hM4Di.mcherry (image taken at A/P = −0.8 from Bregma− low-magnification images shown above with atlas overlay). Fos+ cell density in CeC of sham vs. CAV2-DIO-hM4Di.mcherry injected mice (n = 3 mice/group, averaged from 6 sections control vs. 7 sections hM4Di). Welch’s t-test: **p = 0.0091.
Figure 5—figure supplement 2. Role of nucleus basalis and posterior substantia innominata (NBM/SIp) in associative threat memory recall.

Figure 5—figure supplement 2.

(A) Top: NBM/SIp of Chat-IRES-Cre mice was directly targeted with an AAV9-DIO-hM4Di.mCherry (DREADD) or AAV9-DIO-eCFP (sham) and cholinergic neurons were silenced during recall. Bottom: Freezing behavior during recall (pre-tone period, recall tone response) (repeated measures (RM) two-way analysis of variance (ANOVA; GLM) Time × Group p = 0.0017); Sham mice displayed significant increase in freezing to the recall tone (pre-tone vs. recall tone response, ***p = 0.0004). Silencing the NBM/SIp significantly reduced freezing to the tone (p = **0.0058) (sham, n = 5, hM4Di n = 4).

We injected cav.hM4DiBLA or sham control mice with CLZ 10 min prior to initiating cue-conditioned threat learning (Figure 5B) or 10 min prior to the memory recall session (Figure 5C). In both experiments, mice were sacrificed 45–60 min following recall and assessed for Fos immunoreactivity (IR) in the BLA. We found that DREADD-mediated silencing of BLA-projecting cholinergic neurons during training alone blunted recall-induced freezing behavior and activation of BLA neurons (Figure 5B: freezing behavior, sham vs. cav.hM4DiBLA (Recall Tone Response), p < 0.0001, Figure 5b’, b’’: Fos density, sham vs. cav.hM4DiBLA p = 0.0286). Similarly, DREADD-mediated silencing of BLA-projecting cholinergic neurons during recall alone also reduced recall-induced freezing and activation of BLA neurons (Figure 5C: freezing behavior, sham vs. cav.hM4DiBLA (Recall Tone Response) p = 0.0279, Figure 5c’, c’’: Fos density, sham vs. cav.hM4DiBLA p = 0.0317). Mice in both sham groups showed equivalent freezing behavior (Figure 5B, C, gray boxes; comparing sham groups, p = 0.8155) and density of Fos-IR cells (Figure 5b’, b’’, black circles; comparing sham groups, p = 0.5273) indicating that 0.1 mg/kg CLZ alone (in the absence of DREADD expression) did not alter Fos expression or expression of the learned threat response behavior. Thus, activity of BLA-projecting cholinergic neurons is required during both training and recall for recall induction of Fos expression in BLA neurons and freezing behavior. Preventing cholinergic neuron activity during either training or recall significantly reduced the density of Fos+ BLA neurons and tone-induced freezing.

Differences in recall-induced Fos expression between sham and cav.hM4DiBLA mice were maximal in rostral portions of the BLA (between bregma −0.8 and −1.4 mm) (Figure 5—figure supplement 1B). This region of the rostral BLA has been shown to contain genetically distinguishable neurons that are activated by aversive stimuli and preferentially project to the capsular portion of the central amygdala (CeC), a region known to drive freezing behavior (Kim et al., 2016; Kim et al., 2017). We examined the CeC of mice in which BLA-projecting BFCNs were silenced during recall and found significantly reduced Fos+ cell density in cav.hM4DiBLA mice compared to control mice (Figure 5—figure supplement 1C sham vs. cav.hM4DiBLA p = 0.0091). Thus, silencing cholinergic input to the BLA altered activation of BLA circuits involved in the execution of defensive behaviors.

Mapping BLA-projecting BFCNs infected by CAV2-DIO-hM4Di revealed that the majority of the cholinergic input to the BLA originates in the NBM/SIp (Figure 5A). As such, we delivered AAV9-DIO-hM4Di.mCherry or AAV9-DIO-eCFP (sham mice) into the NBM/SIp of Chat-IRES-Cre mice (Figure 5—figure supplement 2). Both hM4Di and eCFP animals were injected with CLZ 10 min prior to the recall session. Animals in which NBM/SIp cholinergic neurons were silenced during the recall session did not show increased freezing in response to tone (Figure 5—figure supplement 2A, sham, gray boxes: Pre-Tone to Recall Tone Response, p = 0.0004; cav.hM4DiNBM, red boxes: Pre-Tone to Recall Tone Response, p > 0.9999). Thus, silencing NBM/SIp BFCNs was sufficient to block expression of the learned threat response behavior.

Recall-induced activation of NBM/SIp cholinergic neurons correlates with the degree of threat response behavior

During recall, we observed variability in individual freezing responses to the conditioned tone. Based on their responsiveness, we stratified the mice into two groups – high and low responders. ‘High Responders’ were defined as mice who showed a >10 percentage points increase in time spent freezing in response to the tone compared to the pre-tone period (see methods for further details). Mice with <10 percentage points increase in time spent freezing in response to the tone compared to the pre-tone period were defined as ‘Low Responders.’ When stratified as high or low responders according to this criterion, only High Responders showed a statistically significant increase in freezing during the recall tone compared to the pre-tone period (Figure 6A; Pre-tone vs. tone: High Responders, p = 0.0016; Low Responders, p > 0.9999). High Responders showed more freezing compared to Low Responders specifically during the recall tone presentation (High vs. Low responders: recall tone blue shading, p = 0. 0454). ‘High Responders’ spent more time freezing in response to the tone compared to the pre-tone period (Figure 6B).

Figure 6. The extent of cholinergic neuronal activation in the anterior nucleus basalis and posterior substantia innominata (NBM/SIp) co-varies with the behavioral performance during threat memory recall (see Figure 6—figure supplements 1 and 2).

(A) Behavioral performance (freezing) from recall session showing High (black, n = 9) and Low (gray, n = 9) responding mice. High Responders show significantly higher freezing to recall tone whereas low responders do not (two-way repeated measures (RM) analysis of variance (ANOVA)). Interaction effect (Time × Group classification, p = 0.0001; Time, p = 0.0042. High vs. Low Responder *p = 0.0454; Pre-tone vs. Tone: High Responder, p = 0.0016; Low Responder, p > 0.9999). All multiple comparisons were Bonferroni corrected. (B) Quantification of change in freezing responses during recall session in Low and High Responders (pre-tone to tone). Dotted line delineates 10% points change in freezing, which was set as criteria for separating the two populations (see Methods for rationale on stratification criteria, n = 9 Low Responder, n = 9 High Responder). (C) Mice injected in the NBM/SIp with ADCD-mCherry underwent training on Dox and recall off Dox to label recall-activated NBM/SIp cholinergic neurons (n = 11). (D) Quantification of change in number of cholinergic neurons activated (ADCD+) in low or High Responders relative to the home cage. The number of ADCD+ neurons differed significantly between Low and High Responders (Mann–Whitney test, **p = 0.01) (n = 3 Low Responder, n = 8 High Responder). Gray shading represents the range of fold-change in ADCD+ cells in individual home cage mice relative to the average of all home cage mice (n = 5) (Mann–Whitney test, home cage vs. Low Responder, p > 0.9999; home cage vs. High Responder, p = 0.0121). (E) Schematic showing anatomical distribution of ADCD-labeled NBM/SIp basal forebrain cholinergic neurons (BFCNs) activated during recall across the anterior (bregma ~−0.8 mm) to posterior (bregma ~−1.3 mm) extent of the NBM/SIp in High Responders (Top) vs. Low Responders (Bottom). Red circles highlight region of notable difference between High and Low responding mice.

© 2024, BioRender Inc

Figure 6 was created using BioRender, and is published under a CC BY-NC-ND license. Further reproductions must adhere to the terms of this license.

Figure 6.

Figure 6—figure supplement 1. Reactivation of training-activated cholinergic neurons scales with associative threat learning and with behavioral performance during memory recall.

Figure 6—figure supplement 1.

(A) Left: Strategy for labeling activated nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic during training and recall. Chat-IRES-Cre × Fos-tTA/shGFP mice (n = 14) were injected in the NBM/SIp with ADCD-mCherry virus. During Session 1 (off Dox) mice either remained in their home cage, were exposed to three tones (Tone alone), or were exposed to three tone–shock pairings. During session 2 (recall session, on Dox, conducted 72 hr after training), mice remained in home cage or were exposed to a single tone. Cholinergic neurons activated during training express mCherry (red during training), and cells activated during recall transiently express GFP (green during recall). Non-activated neurons are shown in gray. Right: Quantification of percent reactivated cholinergic neurons (activated both during session 1 and during session 2 vs. the total number of cholinergic neurons activated during session 1). Home cage (n = 4 mice), tone only (n = 2 mice) and tone + shock (n = 11 mice – High Responders (based on threshold of 10%-point change in freezing to recall tone) (n = 2 mice), Low responders (n = 6 mice)) conditions. Lines in scatter plots represent medians. (B) Left: Strategy for retrograde targeting of hM4Di DREADD to basolateral amygdala (BLA)-projecting cholinergic neurons. Right: Silencing BLA-projecting cholinergic neurons increases the proportion of low responding mice. Control and hM4Di expressing mice to which clozapine was injected 10 min prior to training (left) or recall (right) were stratified into Low Responders (black fill) and High Responders (white fill). Silencing BLA-projecting cholinergic neurons during training resulted in 100% Low Responders during recall (sham 1 of 5 Low Responders; hM4Di 7 of 7 Low Responders), whereas silencing BLA-projecting cholinergic neurons during recall shifted the proportion of Low Responders from ~11% in the sham (1 of 9) to 50% (3 of 6).
Figure 6—figure supplement 2. Recall-activated basal forebrain cholinergic neurons (BFCNs) in nucleus basalis and posterior substantia innominata (NBM/SIp) of ‘High Responder’ mice are located at a basolateral amygdala (BLA)-projecting locus in the cholinergic basal forebrain.

Figure 6—figure supplement 2.

(A) Representative images of the NBM/SIp at bregma locations −0.8, 1.0, and 1.3 mm from a Fos-tTA/shGFP High Responder (top) and Low Responder (bottom) mouse sacrificed 2.5 hr after recall. Brain sections were stained for ChAT (magenta), and GFP to amplify Fos signal (green). ChAT+ and Fos+ co-labeled cells are marked by yellow arrowheads. Scale bar = 50 µm. Insets show magnified images of ChAT+ and Fos+ neurons. Scale bar = 10 µm. (B) Representative images from Chat-IRES-Cre mice injected with CAV2-DIO-hM4Di.mCherry. mCherry-expressing (red) represent BLA-projecting cholinergic neurons retrogradely labeled. Scale bar = 100 µm. (C) Quantification of BLA-projecting BFCNs (red dots, left y-axis) and Recall-activated BFCNs (black dots, right y-axis) in the NBM/SIp stratified by bregma locations (n = 4 mice).

We next examined whether there was a relationship between the extent of freezing and the engagement of the cholinergic neurons. Since the majority of training-activated cholinergic neurons were reactivated during recall (in high responding mice −~82%, Figure 6—figure supplement 1A), we labeled cholinergic neurons activated during the recall session with ADCD-mCherry (on dox during training, off dox during recall; Figure 6C). Next, we quantified the fold change in the number of mCherry+ neurons in each group relative to corresponding home cage control mice (Figure 6D). While there was no difference in mCherry expression in Low Responders compared to the home cage group (fold change ~1, p > 0.9999), High Responders displayed a threefold increase (p = 0.0121) in mCherry+ cells (High Responders vs. Low Responders, p = 0.0121, Figure 6D).

Mapping of recall-activated NBM/SIp cholinergic neurons revealed that activated BFCNs in ‘High Responder’ mice were in anatomically distinct regions from those in ‘Low Responder’ mice (Figure 6E). In a different cohort of ‘wild-type’ mice, we assessed Fos and ChAT expression following recall and found that in the Low Responders, few ChAT and Fos co-labeled neurons were found. These co-labeled cells were located in caudal regions of the NBM/SIp (~Bregma −1.3; Figure 6—figure supplement 2A – bottom row). In High Responders an additional population of activated cholinergic neurons in more rostral portions of the NBM/SIp was found (~Bregma −0.8; Figure 6—figure supplement 2A – top row). Thus, a discrete population of activated cholinergic neurons in the rostral NBM/SIp is present in mice that respond to the learned threat. When comparing retrograde mapping of BLA-projecting cholinergic neurons using CAV2-DIO-hM4Di.mCherry (Figure 5) to the distribution of ADCD-mCherry-labeled activated neurons (Figure 6), we find a similar distribution along the rostro-caudal axis of the NBM/SIp (Figure 6—figure supplement 2B, C).

Finally, we examined the proportion of high and low responding mice in our experiments where we silenced BLA-projecting cholinergic neurons either during training or during recall (Figure 5B, C). We found that under sham conditions (no cholinergic silencing), 80–90% of the mice were ‘High Responders’. Silencing BLA-projecting cholinergic neurons during training shifted the proportion such that 100% of the mice were ‘Low Responders’ (Figure 6—figure supplement 1B sham vs. cav.hM4DiBLA inhibition during training). Silencing BLA-projecting cholinergic neurons during recall resulted in ~50% of the mice being ‘Low Responders’ (Figure 6—figure supplement 1B, sham vs. cav.hM4DiBLA inhibition during recall). Thus, silencing BLA-projecting cholinergic neurons only during recall resulted in an all-or-none behavioral phenotype (50:50 chance of becoming a High or Low Responder).

Cholinergic neurons activated during threat memory recall have altered intrinsic excitability

Changes in excitability of neurons have been consistently associated with the threat memory engram (Zhang and Linden, 2003; Zhou et al., 2009; Cai et al., 2016; Rashid et al., 2016; Pignatelli et al., 2019). We asked whether cholinergic neurons activated during memory recall differed in their intrinsic excitability compared to non-activated cholinergic neurons. To do this, we prepared acute brain slices from Fos-tTA/shGFP mice for electrophysiological recording of activated (Fos−GFP+) and non-activated (Fos−GFP−) NBM/SIp neurons two and a half hours after the recall session or from mice that remained in their home cage. Cholinergic identity was verified post-recording by single-cell reverse transcriptase polymerase chain reaction (scRT-PCR) of each recorded cell (Figure 7A).

Figure 7. Nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic neurons show increased intrinsic excitability following threat memory recall (see Figure 7—figure supplement 1).

(A) Schematic of electrophysiological profiling of activated (Fos−shGFP+) vs. non-activated (Fos−shGFP−) neurons from mice following recall or in home cage mice, with post hoc identification of cholinergic identity by single-cell RT-PCR and evaluation of chat expression. (B) Representative traces following injection of current into a Fos−shGFP− NBM/SIp cholinergic neuron (ChAT+/Fos−). Red line denotes action potential (AP) half-width measurement. (C) Representative traces following step current injection in Fos-shGFP+NBM/SIp cholinergic neuron (ChAT+/Fos+). Red line denotes AP half-width measurement. (D–F) Population data (dot plot + line at median) for the electrophysiological properties of post hoc identified cholinergic neurons. Analyses assess passive and active membrane properties including AP (D) half-width, (E) rheobase, and (F) maximal firing rate in response to 200–500 ms depolarization from rest potential (−60 mV), from home cage (HC; n = 10–11 ChAT+ neurons from 10 to 11 mice) and following recall to tone alone (n = 10 ChAT+Fos−shGFP− neurons from 5 mice vs. n = 11 ChAT+Fos−shGFP+ neurons from 6 mice). (D) Kruskal–Wallis tests; AP half-width: p = 0.0054 (Dunn’s Corrected p-values: HC vs. Fos−shGFP−: p = 0.8971, HC vs. Fos−shGFP+: ***p = 0.0006, Fos−shGFP− vs. Fos−shGFP+: *p = 0.0206). (E) Rheobase: KW = p = 0.05 (Dunn’s Corrected p-values: HC vs. Fos−shGFP−: p = 0.6153, HC vs. Fos−shGFP+: p = 0.0938, Fos−shGFP− vs. Fos−shGFP+: *p = 0.0228). (F) Max firing rate: p = 0.0032 (Dunn’s Corrected p-values: HC vs. Fos−shGFP−: p = 0.3206, HC vs. Fos−shGFP+: **p = 0.003, Fos−shGFP− vs. Fos−shGFP+: **p = 0.0034).

Figure 7.

Figure 7—figure supplement 1. Latency and afterhyperpolarization (AHP) amplitudes significantly differed in recall-activated nucleus basalis and posterior substantia innominata (NBM/SIp) cholinergic neurons.

Figure 7—figure supplement 1.

(A–F) Additional plots of population data on NBM/SIp cholinergic neurons presented as dot plot + median for the electrophysiological properties of post hoc identified cholinergic neurons from home cage (HC, n = 10–11 ChAT+ neurons from 10 to 11 mice) compared with those of Fos− and Fos+ cholinergic neurons following recall to tone alone (Fos−shGFP−: n = 10 ChAT+ neurons from 5 mice; Fos−shGFP+: n = 11 ChAT+ neurons from 6 mice). Kruskal–Wallis tests (A) resting membrane potential, p = 0.5272 (RMP); (B) action potential amplitude p = 0.1457; (C) threshold, p = 0.9635; (D) AHP half-width, p = 0.3235; (E) latency, p = 0.012 (Dunn’s Corrected p-values: HC vs. Fos−shGFP−: *p = 0.0163, HC vs. Fos−shGFP+: **p = 0.0062, Fos−shGFP− vs. Fos−shGFP+: p = 0.7394); and (F) AHP amplitude, p = 0.0174 (Dunn’s Corrected p-values: HC vs. Fos−shGFP−: p = 0.3702, HC vs. Fos−shGFP+: **p = 0.0041, Fos−shGFP− vs. Fos−shGFP+: p = 0.0952). (G) Population data presented as dot plot + median for max firing rate of post hoc identified cholinergic neurons from home cage, mice sac’d 2.5 hr following training, sac’d 2.5 hr post-recall (D0 – same as Figure 7F. Replotted for comparison here), and ADCD+ NBM neurons tagged during recall and identified 3 days post-recall (D3) and 5 days post-recall (D5) from Chat-IRES-Cre × Fos-tTA/shGFP mice injected with ADCD-mCherry with Dox off during recall. Kruskal–Wallis test, p = 0.0011 (D0 vs. all other groups *p < 0.05 – refer to stats table for details on all comparisons).

Cholinergic neurons that were Fos+ following the recall session differed significantly from Fos− cholinergic neurons (Figure 7B, C) and from cholinergic neurons from home cage mice. Properties that showed significant differences included: action potential (AP) half-width, rheobase, and maximum firing rate (Figure 7D; half-width: HC vs. Fos−shGFP+ p = 0.0006, Fos−shGFP− vs. Fos−shGFP+ p = 0.021; Figure 7E; rheobase: Fos−shGFP− vs. Fos−shGFP+ p = 0.023; Figure 7F; max firing rate: HC vs. Fos−shGFP+ p = 0.003, Fos−shGFP− vs. Fos−shGFP+ p = 0.0034) as well as latency to fire (Figure 7—figure supplement 1E; latency: HC vs. Fos−shGFP+ p = 0.0062) and afterhyperpolarization (AHP) amplitude (Figure 7—figure supplement 1F, HC vs. Fos−shGFP+ p = 0.0041). Resting membrane potential, AP amplitude, AP threshold, and AHP half-width did not differ (Figure 7—figure supplement 1A–D).

We also compared the firing rate of cholinergic neurons in home cage mice with those expressing Fos two and a half hours after training or at longer intervals following recall (measured 2.5 hr (Fos−shGFP) and at 3 and 5 days (ADCD labeling during recall) after the recall session, Figure 7—figure supplement 1G). We found no differences in firing rate between home cage cholinergic neurons and cholinergic neurons that expressed Fos after training: that is, the change in firing rate was only seen in cholinergic neurons activated during recall. This increase in maximal firing rate seen after recall returned to baseline within 3–5 days (compared to recall D0, p < 0.05 for all).

Distinct subsets of BLA-projecting cholinergic neurons differentially contribute to learned vs. innate threat processing

Given the importance of BFCNs in a learned threat paradigm, we next asked whether these cells participate in innate threat responses as well. We stimulated an innate threat response by exposing Fos-tTA/shGFP mice to predator odor (mountain lion urine; Figure 8A; Blanchard and Blanchard, 1990). Exposed mice increased active and passive defensive behaviors compared to mice exposed to a saline wetted pad, including freezing (Figure 8A, p = 0.028), avoidance (Figure 8—figure supplement 1B, left, p = 0.0012), and defensive digging (Figure 8—figure supplement 1B, right, p = 0.023).

Figure 8. Distinct population of basolateral amygdala (BLA)-projecting cholinergic neurons contribute to innate threat processing (see Figure 8—figure supplement 1).

(A) Fos-tTA/Fos–shGFP mice were placed in chambers containing a gauze pad spotted with either saline or with mountain lion urine (predator odor). Defensive behaviors were monitored for 5 min. Mice froze significantly more in the presence of predator odor than saline (Mann–Whitney, *p = 0.028). (B) Basal forebrain sections from the ventral pallidum (VP/SIa) of Fos-tTA/shGFP mice were immunostained for ChAT and GFP 45 min following odor exposure. Predator odor-activated cholinergic neurons (GFP+/ChAT+) were quantified. Predator odor exposure increased the number of activated cholinergic neurons in the VP/SIa (Mann–Whitney: **p = 0.0023), n = 7 control and n = 6 odor exposed mice. (C) Fast Blue was injected into the BLA to retrogradely label BLA-projecting neurons 6 days prior to odor exposure. After exposure to predator odor sections from the basal forebrain were immunostained with antibodies recognizing ChAT and Fos and the numbers of activated cholinergic neurons were counted (ChAT+Fos+/total ChAT+). In the VP/SIa over 90% of BLA-projecting cholinergic neurons were activated (ChAT+ in red, Fos+ in green, Fast Blue in blue, n = 3 mice). Dotted line indicates % of Fos+ cholinergic neurons in the home cage group in this experiment. (D) Chat-IRES-Cre mice injected in the BLA with a control virus (AAV9-camk2a-GFP) alone (sham) or in combination with CAV2-DIO-hM4Di were exposed to predator odor following injection with clozapine (CLZ). Freezing behavior was measured during a 5-min exposure (scatter plot, bar indicates mean; sham – black, hM4Di – red). Silencing BLA-projecting cholinergic neurons significantly blunted the freezing response (Mann–Whitney: *p = 0.019; sham: n = 6; hM4Di: n = 4 mice).

© 2024, BioRender Inc

Figure 8 was created using BioRender, and is published under a CC BY-NC-ND license. Further reproductions must adhere to the terms of this license.

Figure 8.

Figure 8—figure supplement 1. Predator odor exposure activates VP/SIa cholinergic neurons.

Figure 8—figure supplement 1.

(A) Left: Fos-tTA/shGFP mice were exposed to predator odor or saline and sacrificed 2.5 hr later. Brain sections were processed for ChAT immunohistochemistry. Representative images of hDB, VP/SIa, and nucleus basalis and posterior substantia innominata (NBM/SIp) from saline exposed and predator odor exposed mice showing cholinergic neurons (red) expressing Fos-shGFP (green). Scale bar = 100 µm. Insets show magnified images of ChAT+Fos−shGFP+ neurons. Scale bar = 10 µm. Right: Quantification of % activated cholinergic neurons in the respective regions. Cholinergic neurons in the hDB (top, p = 0.2465) or NBM/SIp (bottom, p = 0.4441) were not significantly activated compared to saline controls. n = 7 control and n = 6 odor exposed mice. VP/SIa data (middle, **p=0.0023) from Figure 8B replotted for comparison. (B) Top: Schematic of behavior paradigm. Bottom: Number of contacts with odor pad (left) and number of digging bouts (right) during a 5-min saline or predator odor exposure (saline, gray: n = 9 mice and predator odor, black: n = 6 mice). Mice displayed defensive behaviors, including avoiding the odor pad and increased digging bouts in response to the predator odor. Mann–Whitney test: avoidance ***p = 0.001; digging: *p = 0.023. (C) Top: Schematic of targeting hM4Di to basolateral amygdala (BLA)-projecting cholinergic neurons and subsequent behavioral paradigm. Clozapine was administered to all mice 10 min prior to predator odor exposure. Bottom: Number of contacts with odor pad (left) and digging bouts (right) during a 5-min predator odor exposure in clozapine injected mice expressing hM4Di in BLA-projecting cholinergic neurons (sham n = 6 mice and hM4Di n = 5 mice). Silencing BLA-projecting cholinergic neurons did not affect avoidance behavior (# odor pad contacts) or digging. Mann–Whitney test: avoidance p = 0.8485; digging: p = 0.0714.

We quantified the number of cholinergic neurons expressing Fos (Fos−shGFP+) after saline or predator odor exposure (Figure 8B/Figure 8—figure supplement 1A; Fos−shGFP+/ChAT+). The number of Fos−shGFP expressing cholinergic neurons was significantly elevated in the predator odor exposed group in the VP/SIa (Figure 8B, Figure 8—figure supplement 1A – middle row, p = 0.0023), but not NBM/SIp (Figure 8—figure supplement 1A – bottom row, p = 0.4441), or the hDB (Figure 8—figure supplement 1A – top row, p = 0.2465).

VP/SIa cholinergic neurons formed the second largest source of cholinergic input to the BLA in our retrograde mapping experiments (Figure 5A). Since VP/SIa cholinergic neurons were found to be activated during predator odor exposure, rather than NBM/SIp or hDB cholinergic neurons, we asked if the BLA-projecting pool of VP/SIa cholinergic neurons was activated by predator odor exposure. We injected the retrograde tracer FB into the BLA of Fos-tTA/shGFP mice and then exposed them to either saline (control) or predator odor (Figure 8C, left). FB labeled approximately 30% of ChAT-IR neurons located in the VP/SIa (data not shown). Nearly, the entire subset of BLA-projecting VP/SIa cholinergic neurons (median 94% ± Std.dev 12.5) were also GFP+ (Figure 8C, right).

To determine whether activity of these BLA-projecting cholinergic neurons was necessary for mice to freeze in response to predator odor, we used CAV2-DIO-hM4Di to silence BLA-projecting cholinergic neurons. Silencing during predator odor exposure resulted in significantly less freezing compared to sham mice (Figure 8D, sham vs. cav.hM4DiBLA p = 0.019). Other measures of active avoidance of the predator odor were not significantly altered by silencing BLA-projecting cholinergic neurons (Figure 8—figure supplement 1C; avoidance p = 0.8485; defensive digging p = 0.0714). These data support the conclusion that activity of BLA-projecting cholinergic neurons is critical for normal freezing behavior in response to innate threat. Taken together, we find that distinct populations of BLA-projecting BFCNs are involved in associative threat learning and the response to innately threatening stimuli.

Discussion

A small number of sparsely distributed cholinergic neurons in the basal forebrain provide extensive innervation to most of the brain. These cholinergic neurons and their network of axonal terminal fields play a critical role in modulating cognitive processes (Ballinger et al., 2016; Záborszky et al., 2018; Ananth et al., 2023).

To begin addressing whether the cholinergic system encodes stimulus-specific information, or whether it is generally recruited with salient experiences we monitored ACh release in the BLA during threat learning and retrieval. We anatomically mapped and electrophysiologically characterized behaviorally relevant BFCNs, and then investigated the contribution of different subsets of BFCNs to threat response behaviors. Taken together, our results demonstrate distinct populations of cholinergic neurons that are an integral part of encoding a learned threat memory or contribute to innate threat responses.

Cholinergic modulation of associative threat learning

In the BLA, several molecular changes occur in response to learning conditioned stimulus-unconditioned stimulus (CS-US) associations, including new gene expression and protein synthesis (Sears et al., 2014). We used chemogenetics for projection-specific, cell-type-specific silencing of cholinergic neurons. We used CLZ activation of hM4Di, acting at cholinergic cell bodies and/or cholinergic terminals (Krashes et al., 2011; Ray et al., 2011; Ferguson et al., 2013; Stachniak et al., 2014; Zhang et al., 2017; Jin et al., 2019; Nishioka et al., 2020; O’Neal et al., 2020) to silence BLA-projecting BFCNs during training or during recall. This resulted in loss of freezing behavior as well as significantly reduced density of Fos expressing neurons in the BLA following recall (Figure 5). This reduction of Fos expression in the BLA indicates that cholinergic signaling in the BLA contributes to appropriate BLA engagement during the acquisition and recall of threat memory. While our experiments did not directly measure the BLA engram per se (i.e. activation–reactivation of the same neurons within the BLA), our data support the hypothesis that BLA-projecting cholinergic neurons play a critical role in the formation and/or activation of the BLA engram.

We have previously demonstrated that activation of presynaptic ACh receptors can induce sustained potentiation of glutamate release (McGehee et al., 1995; Zhong et al., 2008; Jiang et al., 2013; Zhong et al., 2013; Zhong et al., 2015; Jiang et al., 2016; Zhong et al., 2017). BLA neurons recruited during memory recall exhibit increased presynaptic glutamatergic activity (Nonaka et al., 2014). We further demonstrated that the increased glutamatergic transmission in BLA was dependent on presynaptic nicotinic acetylcholine receptors (nAChRs) located on glutamatergic terminals in the BLA, and that nAChR activation in the BLA was necessary for acquisition of conditioned threat memories (Jiang et al., 2016). Based on these findings, we propose that chemogenetic silencing of BLA-projecting cholinergic neurons during threat learning or during recall results in loss of Fos expression due to alterations in presynaptic glutamatergic transmission resulting in disruption to the formation and/or recruitment of the BLA engram.

BFCNs ‘learn’ to respond to the conditioned stimulus

In this study, we used a genetically encoded ACh sensor (GRABACH3.0) to monitor endogenous ACh release in the BLA during threat learning and recall. First, we found that foot shock rapidly and reliably evoked ACh release, in line with previous observations (Hangya et al., 2015; Jing et al., 2020). When we examined responses to the tone (CS, Figure 1—figure supplement 1), we did not detect a significant increase in ACh in the BLA in response to a naive, unexpected tone. However, following conditioning, when mice were exposed to the conditioned tone in a novel environment 24 hr later, we observed robust ACh release in the BLA compared with the naive tone (Figure 1D). This enhancement of ACh release supports the notion that BLA-projecting BFCNs undergo physiological changes which allow robust responsiveness to previously naive sensory stimuli. When mice were exposed to tones in the absence of foot shocks and then exposed to the same tone 24 hr later, we did not detect increased ACh release in the BLA (Figure 1—figure supplement 2). Thus, plasticity of ACh release in the BLA in response to the tone requires pairing of the tone with a salient stimulus such as a foot shock.

Changes in excitability of Fos+ cholinergic neurons

It has been proposed that alterations to synaptic weights and changes in ionic conductance resulting from learning-induced transcriptional programs allow for increased response fidelity during memory retrieval (Yap and Greenberg, 2018). To assess whether such changes occurred in recruited cholinergic neurons following memory retrieval, we recorded properties of neuronal excitability from activated NBM/SIp BFCNs (Fos+) and compared them with Fos− BFCNs recorded in the same brain slices (Figure 7). Recall-activated NBM/SIp cholinergic neurons showed increased excitability which lasted for at least several hours following threat memory retrieval, returning to baseline within days. This finding is in line with previous reports of learning-associated changes in electrical properties, which are found shortly after recall, but disappear at later time points despite the persistence of the learned behavior (Moyer et al., 1996; Pignatelli et al., 2019). Observed changes in the electrophysiological properties were not see in Fos-shGFP+ cholinergic neurons immediately following training. Thus, many of the changes in electrical properties we observed were specific to recall-activated cholinergic neurons. Within recall-activated cholinergic neurons we find several changes consistent with an increased excitability such as decreased AP half-width, decreased rheobase, and an increase in maximum firing rate. Common features of activated neurons previously reported include similar increases in firing rate, with reductions in adaptation, decreased duration of post-burst AHP, decreased AHP amplitude, and synaptic alterations (Whitaker and Hope, 2018).

Differential contribution of distinct BLA-projecting BFCNs in learned vs. innate threat processing

Amygdala microcircuits play an important role in the regulation of active vs. passive avoidance behaviors (Rickenbacher et al., 2017; Terburg et al., 2018). Our finding that silencing cholinergic input to the BLA resulted in a selective loss of threat-motivated freezing behavior supports potential specificity of cholinergic modulation within BLA microcircuits for freezing, but not active, defensive behaviors. We found that BLA-projecting cholinergic neurons were necessary for freezing in response to a learned threat-associated cue (Figure 5), and for freezing in response to the innately threatening predator odor (Figure 8). Direct silencing of NBM/SIp cholinergic neurons attenuated learned threat induced freezing. Instead, predator odor activated BLA-projecting VP/SIp cholinergic neurons and resulted in a freezing response. Based on these data, we propose that distinct populations of BLA-projecting BFCNs control freezing in response to fundamentally distinct threatening situations (learned vs. innate). Additionally, we note that while silencing BLA-projecting BFCNs did reduce freezing in response to predator odor exposure, it did not alter avoidance of the odor pad indicating that threat detection was still intact in these mice.

Memory encoding in neuromodulatory systems

Our study joins a growing literature demonstrating stimulus-encoding and rapid stimulus-contingent responses in various neuromodulatory neurons indicating that plasticity within subcortical modulatory circuits might represent a critical component of normal learning and memory recall. The BLA receives various modulatory inputs including dopamine (DA) from the ventral tegmental area (VTA) (Tang et al., 2020), noradrenaline (NA) from the locus coeruleus (LC) (Uematsu et al., 2017), and ACh from the basal forebrain. Including our present study, all three of these modulatory systems have been shown to be engaged during associative threat learning and retrieval. Each modulatory system seems to respond rapidly and robustly to aversive stimuli like mild electrical shocks, and activity within these systems during conditioning (i.e. during CS–US pairing) is critical for generation of freezing behavior during memory recall (Uematsu et al., 2017; Tang et al., 2020). VTA dopaminergic neurons have also been shown to display plasticity in tone responsiveness such that a naive tone does not result in significant firing of DA neurons (Tang et al., 2020). However, following three pairings of the tone with shocks, VTA DA neurons begin responding to tone presentations with millisecond latencies, a response that is sustained the following day during memory retrieval. A majority of the shock-responsive DA neurons were also found to acquire tone responsiveness following pairing, a finding replicated within the cholinergic system in our study. While shock rapidly activates LC NA neurons, conditioned tone-related responses in these neurons seems to be slow, occurring on average several seconds following tone presentation (Uematsu et al., 2017). How signaling by these different modulators interacts in the BLA and informs plasticity of BLA neurons is an intriguing question.

In addition to these modulators, peptides such as oxytocin have also been shown to participate in threat memory formation. A recent study demonstrated presence of a threat memory engram within the hypothalamic oxytocinergic projection to the amygdala (Hasan et al., 2019). Interestingly, upon conditioning these neurons demonstrate a transmitter preference switch, releasing glutamate in the amygdala. Thus, subcortical neuromodulatory and peptidergic systems might display unique mechanisms of engram-related biophysical changes that have not been found in traditionally studied systems.

We demonstrate at least two populations of BLA-projecting cholinergic neurons that are engaged in learned vs. innate threat responses. Differences in function of other BLA-projecting BFCNs (NBM vs. HDB) in threat memory formation vs. extinction were recently demonstrated (Hasan et al., 2019; Crimmins et al., 2023), further highlighting that effects of ACh release in the BLA are highly specific to which axons release the ACh, despite the dense overlapping terminal fields from different BFCN populations within the BLA. Similar heterogeneity of responses has also been found in the dopaminergic and noradrenergic systems (Azcorra et al., 2023). It is possible that single-cell transcriptomic analyses of the cholinergic basal forebrain may provide insight into the functional heterogeneity observed in our study.

Is there a cholinergic component in the associative threat memory engram?

Studies examining mechanisms of learning and memory in recent years have revived Semon’s theory on memory engrams: learning must result in lasting biophysical changes that form the substrate for retrieval of the learned experience (Semon, 1921; Tonegawa et al., 2015). Josselyn and Tonegawa have recently updated the definition of engram cells, requiring that these be activated by learning, modified by learning, and reactivated by subsequent presentation of the recall-inducing stimuli, resulting in memory retrieval (Josselyn and Tonegawa, 2020). NBM/SIp BFCNs investigated in this study indeed fulfil these criteria as they are activated by learning, show induction of Fos and altered physiological properties with recall, are reactivated by recall, and the reactivation of previously, training-activated BFCNs was necessary for recall behavior.

Multiple studies have used threat and reward learning paradigms in rodents to examine allocation of neurons to memory engrams. These studies have looked for these engram cells in regions such as cortex, amygdala, and hippocampus focusing on glutamatergic pyramidal neurons (Josselyn et al., 2015). However, recent work has demonstrated that memory engrams are distributed across brain-wide networks, and that reactivation of a multi-region engram more closely recapitulates natural recall behavior (Roy et al., 2022).

In addition to the BLA, cholinergic neurons in the NBM/SIp region project to various limbic and sensory regions such as the lateral orbital cortex, cingulate cortex, somatosensory cortex, and mediodorsal thalamus (Ananth et al., 2023). This raises the interesting possibility that the cholinergic signaling modulates various nodes of the threat memory engram circuit in conjunction with the amygdala, allowing for coordinated retrieval of engrams across distributed networks. Such coordinated activation of distributed engrams has been recently demonstrated to more closely recapitulate natural memory retrieval (Roy et al., 2022). Furthermore, functionally related regions have been shown to receive their cholinergic input from the same cholinergic nucleus (Zaborszky et al., 2015). We propose that engram-enrolled cholinergic neurons bind distributed engrams to encode stimulus-convergent, efficient memory retrieval.

Materials and methods

Resource availibility

Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by Lead Contact, Dr. David Talmage ([email protected]).

Materials availability

Plasmids generated in this study have been deposited to Addgene and will be available upon publication under Talmage Lab.

Experimental model and subject details

Adult (3–6 month) male and female Chat-IRES-Cre (B6;129S6-Chattm2(cre)Lowl/J), Jax stock number: 006410 (Rossi et al., 2011), Fos-tTA, Fos−shGFP (TetTag, Jax stock number: 018306, referred to as Fos-tTA/shGFP or Fos-shGFP), and Chat-IRES-Cre × Fos-tTA/shGFP mice were used. Mice within each cage were randomly assigned to experimental and control conditions. In all electrophysiology experiments, hemizygous Fos-tTA/shGFP mice on a C57BL/6 background were used. Mice were housed in a 12-hr light/dark cycle environment that was both temperature and humidity controlled. Mice had free access to food and water. All animal care and experimental procedures were approved by the Animal Care and Use Committees (ACUC) of the National Institute of Neurological Disorders & Stroke (NINDS) (Protocol #1531), SUNY Research Foundation at Stony Brook University (Protocol #1618), and Yale University (Protocol #2019-07895).

Method details

Viral construct

Construction of the ADCD probe

All cloning unless otherwise specified was performed using In-Fusion HD (Clontech). ‘mCherry-P2A’ was amplified using Phusion High-Fidelity DNA Polymerase (NEB) from pV2SGE (obtained as a gift from Dr. Shaoyu Ge Stony Brook University). ‘oChIEF-LoxP-Lox2272’ was amplified from pV2.2 (synthesized gene block from IDT). The two fragments were cloned into pAAV-WPRE linearized by BamHI. The resulting plasmid was linearized by Pml I. ‘7xTetO-LoxP-Lox2272-tTAH100Y.SV40’ was amplified from pV2.1 (synthesized gene block from IDT) and cloned into the Pml I site. The final plasmid was packaged into AAV9 viral particles. Viral packaging was performed by the University of Pennsylvania Vector Core.

Note re: ADCD expression in BLA neurons in the presence of doxycycline

As shown in Figure 2—figure supplement 1C, we noted ‘leaky’ expression of ADCD-mCherry in the presence of doxycycline, in the BLA of Fos-tTA mice when co-injected with a Cre expression vector expressed from a camk2a promoter. Co-injection of camk2a-Cre and ADCD-mCherry into cortex and hippocampus of wild-type (C57) mice was also found to result in ‘leaky’ expression despite the absence of genetically encoded tTA. Injection of ADCD-mCherry in hippocampus of PV-Cre mice did not result in expression similar to injection in Chat-IRES-Cre mice (Figure 2—figure supplement 1A, bottom). These findings underscore the importance of performing the appropriate controls when using these vectors in vivo.

Construction of the ADCD-DREADD probe

‘BglII-hM4Di.mCherry-AscI’ was amplified using CloneAmpTM HiFi PCR Premix (Takara) from pAAV-hSyn-DIO-hM4D(Gi)-mCherry (Krashes et al., 2011) (gift from Dr.Bryan Roth; Addgene plasmid # 44362; http://n2t.net/addgene:44362; RRID:Addgene_44362). A backbone with TRE and Lox sites was ligated with ‘BglII-hM4Di.mCherry-AscI’ using T4 DNA Ligase (NEB). The final plasmid was packaged into AAV9 viral particles. Viral packaging was performed by the University of North Carolina Vector Core.

Stereotaxic surgery and viral delivery

Three- to four-month-old ChAT-IRES-Cre mice were anesthetized and stereotaxically injected bilaterally. Coordinates were calculated based on the Paxinos Mouse Brain Atlas (Franklin, K & Franklin and Paxinos, 1997): BLA (−1.4 mm A/P, ±3.5 mm M/L, −4.8 mm D/V) and NBM (−0.7 mm A/P, ±1.7 mm M/L, −4 mm D/V).

Tracers

3% wt/vol solution of FB (17740-1, Polysciences Inc) was prepared in sterile milliQ water. ~0.2 µl of 3% FB was injected into the BLA bilaterally of Fos-GFP or Chat-IRES-Cre × Fos-tTA/shGFP mice. Mice were euthanized 7 days following injection.

Behavioral testing and analysis

Threat conditioning

All training and assessments were completed with experimenter blind to condition. Both training and recall sessions were analyzed using FreezeFrame v.3 (see below).

Habituation

All mice were handled for a minimum of 5 min daily for three consecutive days before behavioral training began. For DREADD experiments, all mice were additionally habituated to restraint and injection with 100 µl saline administered i.p. daily.

Training

On training day, all chambers were cleaned with 70% ethanol. Mice were placed into the behavioral chamber for a 10-min session which consisted of 3 min of habituation, followed by three tone–shock pairings (30 s 80 dB, 5 kHz tone, co-terminated with a 2-s 0.7-mA foot shock with a 1.5-min interval between each pairing), and finally 2 min of exploration. For DREADD experiments, mice were given 0.1 mg/kg CLZ (administered i.p.) (Sigma-Aldrich) 10 min prior to being placed in the chamber.

Recall

Recall session took place 24–72 hr after completion of the training. To specifically test the response to tone-cued recall, the contextual features of the chambers were altered including texture of the floor, color of the walls, and scent of cleaner (mild lemongrass citrus-based solution). Mice were placed in the behavioral chamber for another 5 min session during which a single tone was delivered (30 s 80 dB 5 kHz tone) 2 min after being placed in the chamber. No shock was administered.

Analysis

Percent time spent freezing was quantified using FreezeFrame v.3 (Actimetrics). Bout duration (defined as minimum required duration when animal is frozen) was set to 1 s, and threshold was manually defined as highest motion index with no movement other than breathing. Percent time spent freezing (defined as periods of no movement) was quantified across the 10 min session in bins of 30 s. The following periods were defined for statistical analysis: Baseline (average of all bins prior to tone onset) and Tone response (average of all bins following tone onset).

High Responders were defined as those mice that exhibited at least a 10 percentage point increase in % time spent freezing in the 30 s bin during the tone from the average of the pre-tone period (e.g. Pre-tone freezing 10% to tone-induced freezing of ≥20%). All other mice were considered Low responders. Prior to any behavioral manipulation, mice showed up to 10% (of total time in given time bin) freezing indicating this level of freezing to be non-associative (potentially related to novelty or generalized fear). This criterion was found to give statistically significant difference between pre-tone vs. tone only for high responders and not for low responders providing further validity to the delineation of the Low and High Responder groups.

Analysis of population composition of High and Low responders (Figure 6—figure supplement 1) was performed within experiment. Cross-experiment comparisons for population composition of High and Low responders were not possible due to differences in conditions and variability within and between cohorts.

Engram labeling

Mice were placed on doxycycline hyclate-containing chow (Cat# TD.08541 Envigo) at least 2 days prior to injection of activity-dependent viral markers. Threat conditioning was performed as mentioned above. During doxycycline withdrawal, mice were transferred to a clean cage to prevent mice from eating dox food that was dragged into the cage or buried in the bedding. To minimize stress, some bedding containing fecal pellets and urine, and nest from the old cage were transferred to the new cage.

Predator odor exposure

Habituation

All mice were habituated to restraint and injection with 100 μl saline administered i.p. daily for 3 days prior to behavioral testing for DREADD experiments. On exposure day, mice were transported to the lab several hours prior to exposure and habituated to the room and ambient sounds.

Exposure

For exposure to predator odors, a vented mouse cage (L 13in × W 7.5in × H 5.5in) with corncob bedding (EnviroDri) was placed in a designated location in a laminar flow hood with overhead fluorescent lighting. Mt. Lion Pee (Maine outdoor solutions LLC) was obtained from https://predatorpeestore.com/ and stored at 4°C. 200 µl of urine was pipetted onto a 3in × 3in 12 ply gauze pad (Cat#6312, Dukal corp.) placed in a polystyrene Petri dish (VWR) at the vented end of the cage. Mice were placed into the cage in the end away from the odor and the cage was covered using a clear plexiglass barrier. Mice were exposed for 5 min and the session was filmed using an overhead digital camcorder (Sony). Following exposure, mice were returned to their home cage or a holding cage in the case of multiple housed mice to prevent any odor transfer. Control mice were exposed to 0.9% saline. For DREADD experiments, mice were given 0.1 mg/kg CLZ (administered i.p.; Sigma-Aldrich) 15 min prior to being placed in the chamber.

Analysis

behavior was manually scored using Jwatcher (v0.9). Defensive digging was defined as vigorous digging performed by the mice using their snout, flinging bedding up and away from the animal. Freezing was defined as immobility without any obvious motion besides breathing. Cloth contacts were defined as front paw touches to the odor pad.

Fiber photometry

Acquisition

Fiber photometry recordings were made using a Doric Lenses 1-site Fiber Photometry System. Signal was recorded using Doric Neuroscience Studio (V 5.3.3.4) via the Lock-In demodulation mode with sampling rate of 12.0 kS/s. Data were downsampled by a factor of 10 and saved as a comma-separated file. For details on connection of the setup refer to Crouse et al., 2020.

Analysis

Preprocessing of the raw data was performed using a MATLAB script provided by Doric. The baseline fluorescence (F0) was calculated using a least mean squares regression over the duration of the recording session. The change in fluorescence for a given timepoint (ΔF) was calculated as the difference between it and F0, divided by F0, and multiplied by 100 to yield % ΔF/F0. The % ΔF/F0 was calculated independently for both the signal (465 nm) and reference (405 nm) channels and a final ‘corrected % ΔF/F0’ was obtained by subtracting the reference % ΔF/F0 from the signal % ΔF/F0 at each timepoint. The corrected % ΔF/F0 was z-scored to give the final ‘Z % ΔF/F0’ reported. Area under the curve was calculated for 1-s duration before (baseline) and 1 s after tone onset. The average of all the baseline periods within each analysis was used as the baseline reading for the AUC analysis.

Electrophysiology

Brain slice preparation

For slice physiology, mice were anesthetized and transcardially perfused with cutting solution (sucrose 248 mM, KCl 2 mM, MgSO4 3 mM, KH2PO4 1.25 mM, NaHCO3 26 mM, glucose 10 mM, sodium ascorbate 0.4 mM and sodium pyruvate 1 mM, bubbled with 95% O2 and 5% CO2) at 40°C. The brain was then rapidly removed and sliced, coronally, at 300 µM in oxygenated cutting solution at 40°C. Prior to recording, slices were incubated in oxygenated incubation solution (sucrose 110 mM, NaCl 60 mM, KCl 2.5 mM, MgCl2 7 mM, NaH2PO4 1.25 mM, NaHCO3 25 mM, CaCl2 0.5 mM, MgCl2 2 mM, glucose 25 mM, sodium ascorbate 1.3 mM, and sodium pyruvate 0.6 mM) at room temperature.

Electrophysiological recording

During recording, slices were superfused with oxygenated artificial cerebral spinal fluid (Jiang et al., 2016). Fos+ neurons were identified by GFP expression. Signals were recording using patch electrodes between 4 and 6 MΩ, a MultiClamp 700B amplifier, and pClamp10 software. Pipette internal solution was as follows: 125 mM K-gluconate, 3 mM KCl, 1 mM MgCl2, 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), 0.2 mM CaCl2, 0.1 mM ethylene glycol tetraacetic acid (EGTA), 2 mM MgATP, and 0.2 mM NaGTP (pH = 7.3). Following recording, cytoplasm was harvested via aspiration for cell-type identification using single-cell RT-PCR. Ten to twelve basic electrical properties were determined and defined as previously described (López-Hernández et al., 2017). Recordings were excluded if they did not meet the following criteria: (1) membrane potential less than or equal to −50 mV, (2) input resistance between 100 and 300 MΩ, (3) series resistance <10 MΩ that was unchanged throughout the recording, and (4) firing a 45 mV AP at rheobase.

Single-cell reverse transcription-PCR

Single-cell samples were pressure ejected into a fresh RT buffer prep (Applied biosystems). Samples were sonicated in a total volume of 20 µl at 40°C for 10 min before addition of RT enzyme mix (Applied Biosystem). Tubes were incubated at 37°C for 60 min and then 95°C for 5 min. Two rounds of amplification (30 cycles each) were done for the detection of Chat transcripts. For the first round of amplification (reaction volume 25 µl) included 2× mastermix, sterile water, 0.2 mM of each primer, 1 ml of cDNA sample. For the second amplification, the reaction included 1 µl of the previous (first-round) PCR product, 2× mastermix, sterile water, and 0.2 mM of each primer. Whole brain cDNA was run in parallel with the single-cell samples. After amplification, the PCR products (159 bp) were analyzed on 3% gels.

Immunohistochemistry

Following perfusion, brains were fixed overnight at 4°C in 4% Paraformaldehyde (PFA) (in 1× phosphate-buffered saline (PBS)) and were then transferred to a 30% sucrose solution (in 1× PBS). Brains were flash frozen in OCT Compound (Tissue Tek) and stored at −80°C until cryosectioning. 50 µm cryosections were mounted onto Superfrost slides (Fisher Scientific) in sets of 3 and allowed to dry overnight at room temperature. Sections were blocked overnight at 4°C in a PBS solution containing 0.3% Triton X-100 and 3% normal donkey serum and then incubated with primary antibody in a PBS-T solution (0.1% Triton X-100 and 1% normal donkey serum), overnight (24 hr at 4°C). The next day, sections were rinsed in PBS-T and incubated in secondary antibody for 2 hr at room temperature in PBS-T along with NeuroTrace-435 (Invitrogen). Sections were treated with an autofluorescence eliminator reagent (EMD Millipore) according to the manufacturer’s guidelines and mounted in Fluoromount-G (Southern Biotech). Details regarding antibodies can be found in the Key Resources Table (KRT).

Quantification and statistical analysis

Imaging and analysis

All imaging was conducted on an Olympus wide-field slide-scanner microscope at 20× magnification (VS-120 and VS-200 systems, Z-step = 3 µm). Images were processed using the cell counter plugin on ImageJ. For Fos+ cell counts in the amygdala, only neurons (Nissl/ Neurotrace positive) with nuclear Fos stain were counted. The amygdala was identified, and a region of interest (ROI) defined using ROI manager in ImageJ. Total area of the ROI was measured and noted. Fluorescence threshold was set to eliminate background fluorescence in ImageJ (defined as hazy background signal detected in space between neurons and white matter). This eliminated non-specific fluorescence and out of focus signals. Fos+ nuclei were then counted using the cell counter plugin.

For ADCD cell counts, mCherry+ neurons at the NBM/SI injection site were counted. NBM was consistently identified as the cluster of cholinergic cell bodies at the base of the internal capsule in the Globus Pallidus and the SI as the area located directly ventral to the GP as denoted by the Paxinos Mouse Brain Atlas (3rd Edition). 100% of the analyzed area of every third brain section was counted (~150 µm apart). Since the NBM/SI regions lack defined boundaries, we present the data as cell counts as opposed to cell density.

For Fos analysis in the BLA, Fos+ cells were counted in the area enclosed within the external and amygdalar capsules. Since the shape of the BLA changes along the anterior–posterior axis, Fos+ cell counts were normalized to the area enclosed within the external and amygdalar capsules and presented as density of Fos+ cells.

Statistical analysis

Statistical analyses were done using GraphPad Prism (GraphPad Software Inc, San Diego, CA, USA), Sigmaplot 12.5 (Systat Software, Inc, San Jose, CA, USA), and OriginPro 9.1 (Origin Lab Corporation, Northampton, MA, USA). Normality of the data was assessed using Shapiro–Wilk and Smirnov–Kolmogorov tests. Data that were not normally distributed according to both normality tests, were analyzed using appropriate non-parametric tests. Detailed information on statistical tests used, p-values, and sample sizes, and other descriptive statistics can be found in the text (figure legends) and/or in the statistical reporting table (Supplementary file 1). Sample sizes for behavior experiments were determined using a power calculation based on effect sizes in pilot experiments with power set to 0.8.

Parametric tests used: Repeated measures (RM) one-way analysis of variance (ANOVA), RM two-way ANOVA, Welch’s ANOVA, paired t-test (two-tailed), Welch’s t-test.

Non-parametric tests used: Mann–Whitney test, Wilcoxon matched-pairs signed rank test, Kruskal–Wallis test, Friedman test.

p-value criteria: *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001.

Acknowledgements

This work was supported by the Intramural Research Program of NINDS. This work was also supported by DA14241, DA037566, and MH077681 to MRP, and early phases by NS022061 and MH109104 to LWR and DAT. RBC was supported by a NINDS Training Grant (T32) NS007224. We thank Dr. Shaoyu Ge (Stony Brook University, NY) for providing reagents and insightful discussions aiding in the conceptualization of the project. We also thank Drs Josh Dubnau and Qiaojie Xiong (Stony Brook University, NY) for providing feedback and discussions on experiments presented in this manuscript. We thank Wendy Akmentin, Dr. Li Bai, and Taylor Muir for expert technical assistance in data curation. Figure schematics were created with BioRender.com.

Appendix 1

Appendix 1—key resources table.

Reagent type (species) or resource Designation Source or reference Identifiers Additional information
Strain, strain background (Mus musculus) Chat-IRES-Cre The Jackson Laboratory B6;129S6-Chattm2(cre)Lowl/J Stock number: 006410
Strain, strain background (Mus musculus) Fos-tTA, Fos−shGFP The Jackson Laboratory TetTag Stock number: 018306
Strain, strain background (Escherichia coli) Stellar Competent Cells, HST08 Takara Cat#636766
Genetic reagent (AAV9) AAV9-camk2a-GCaMP6f-WPRE-SV40 Penn Vector Core
Genetic reagent (AAV9) AAV9-DIO-eCFP This paper, Vector Biolabs Custom made
Genetic reagent (AAV9) AAV9-hSyn-GACh4.3 Vigene Biosciences Inc
Genetic reagent (AAV8) AAV8-DIO-hM4Di-mCherry Addgene Cat#44362
Genetic reagent (AAV9) AAV9-TRE-DIO-oChIEF-mCherry-P2A-tTAH100Y.SV40 This paper plasmid DNA See Methods and
Figure 2—figure supplement 1, can be obtained from Talmage lab.
Genetic reagent (AAV9) AAV9-TRE-DIO-hM4Di-mCherry This paper Cat#169415 Deposited to Addgene, see methods and Figure 2—figure supplement 1, can be obtained from Talmage lab.
Genetic reagent (CAV2) CAV2-DIO-hM4Di-mCherry Dr. EJ Kremer, Institut de Génétique Moléculaire de Montpellier, France
Antibody anti-ChAT (Goat polyclonal) Millipore Cat# AB144P; RRID:AB_2079751 IHC (1:500)
Antibody anti-GFP (Rabbit polyclonal) Thermo Fisher Scientific Cat#: A-11122; RRID:AB_221569 IHC (1:1000)
Antibody anti-GFP (Rabbit polyclonal) Abcam Cat#: ab13970; RRID:AB_300798 IHC (1:500)
Antibody anti-mCherry (Mouse monoclonal) Takara Cat#: 632543; RRID:AB_2307319 IHC (1:500)
Antibody anti-DsRed (Rabbit polyclonal) Takara Cat#: 632496; RRID:AB_10013483 IHC (1:500)
Antibody anti-c-Fos (Rabbit polyclonal) Synaptic Systems Cat#: 226003; RRID:AB_2231974 IHC (1:500)
Antibody anti-Rabbit IgG (H+L)-AlexaFluor 488 (Donkey polyclonal) Thermo Fisher Cat#: A32790; RRID:AB_2762833 IHC (1:1000)
Antibody anti-Rabbit IgG (H+L)-Rhodamine Red-X (Donkey polyclonal) Jackson Immunoresearch Cat#: 711-295-152; RRID:AB_2340613 IHC (1:1000)
Antibody anti-Goat IgG (H+L)- AlexaFluor 594 (Donkey polyclonal) Thermo Fisher Cat#: A-11058; RRID:AB_142540 IHC (1:1000)
Antibody anti-Chicken IgY-Cy2 (Donkey polyclonal) Gift from Dr.Shaoyu Ge, Stony Brook University NY IHC (1:1000)
Chemical compound NeuroTrace 435/455 Blue Fluorescent Nissl Stain Thermo Fisher Cat#: N21479 IHC (1:500)
Recombinant DNA reagent pAAV-hSyn-DIO-hM4D(Gi)-mCherry (plasmid) Addgene Cat#44362
Recombinant DNA reagent pV2SGE (plasmid) This paper Gift from Dr.Shaoyu Ge, Stony Brook University NY Used in the construction of reagent #25
Recombinant DNA reagent pAAV-TRE-DIO-oChIEF-mCherry-P2A-tTAH100Y.SV40 (plasmid) This paper Deposited to Addgene Addgene Cat# 169414
Recombinant DNA reagent pAAV-TRE-DIO-hM4Di-mCherry (plasmid) This paper Deposited to Addgene Addgene Cat# 169415
Sequence-based reagent chat_F IDT PCR primers TCTGGCAACTTCGTCGGA
Sequence-based reagent chat_R IDT PCR primers CTCCTGGGCTGTTACGCAC
Sequence-based reagent pV2.1 – Gene Block 7xTetO-LoxP-Lox2272-tTAH100Y.SV40 IDT Gene block, custom
Sequence-based reagent pV2.2 – Gene Block oChIEF-LoxP-Lox2272 IDT Gene block, custom
Commercial assay or kit In-Fusion HD Cloning Plus Takara/Clontech Cat#: 638920
Commercial assay or kit High-Capacity cDNA Reverse Transcription Kit Applied Biosystems Cat#: 4368814
Peptide, recombinant protein T4 DNA Ligase NEB Cat#M0202S
Peptide, recombinant protein BglII NEB Cat#R0144S
Peptide, recombinant protein AscI NEB Cat#R0558S
Peptide, recombinant protein BamHI-HF NEB Cat#R3136S
Peptide, recombinant protein PmlI NEB Cat#R0532S
Peptide, recombinant protein Phusion High-Fidelity DNA Polymerase NEB Cat#M0530S
Chemical compound, drug Clozapine Sigma-Aldrich Cat#C6305-25MG
Chemical compound, drug Fast Blue Polysciences Inc Cat#17740-1
Other Mt.Lion Urine Maine outdoor solutions LLC Obtained from https://predatorpeestore.com/.
Software, algorithm Prism GraphPad Software Inc RRID:SCR_002798
Software, algorithm Sigmaplot 12.5 Systat Software Inc RRID:SCR_003210
Software, algorithm OriginPro 9.1 Origin Lab Corporation RRID:SCR_014212
Software, algorithm Fiji is just imagej Fiji RRID:SCR_002285
Software, algorithm FreezeFrame v3 Actimetrics RRID:SCR_014429
Software, algorithm MATLAB Mathworks RRID:SCR_001622
Software, algorithm Pre-processing analysis MATLAB Script for FiberPhotometry Doric
Software, algorithm ACh sensor analysis MATLAB script Crouse, Richard B., et al. Elife 9 (2020): e57335

Funding Statement

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

Contributor Information

Lorna W Role, Email: [email protected].

David A Talmage, Email: [email protected].

Joshua Johansen, RIKEN Center for Brain Science, Japan.

Laura L Colgin, University of Texas at Austin, United States.

Funding Information

This paper was supported by the following grants:

  • National Institute of Neurological Disorders and Stroke 1ZIANS009424 to David A Talmage.

  • National Institute of Neurological Disorders and Stroke 1ZIANS009416 to Lorna W Role.

  • National Institute of Neurological Disorders and Stroke NS22061 to David A Talmage.

  • National Institute of Mental Health U01-MH109104 to David A Talmage.

  • National Institute of Mental Health MH077681 to Marina R Picciotto.

  • National Institute on Drug Abuse DA14241 to Marina R Picciotto.

  • National Institute of Neurological Disorders and Stroke NS007224 to Richard Crouse.

  • National Institute of Neurological Disorders and Stroke 1ZIANS009422 to Lorna W Role.

  • National Institute on Drug Abuse DA037566 to Marina R Picciotto.

Additional information

Competing interests

No competing interests declared.

Author contributions

Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing.

Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing.

Data curation, Investigation, Methodology.

Data curation, Software, Formal analysis, Validation, Investigation, Methodology, Writing – review and editing.

Investigation.

Investigation, Methodology.

Investigation.

Investigation.

Investigation.

Investigation.

Investigation.

Resources, Methodology.

Supervision, Funding acquisition, Methodology, Writing – review and editing.

Conceptualization, Supervision, Funding acquisition, Methodology, Project administration, Writing – review and editing.

Conceptualization, Supervision, Funding acquisition, Methodology, Project administration, Writing – review and editing.

Ethics

All animal care and experimental procedures were approved by the Animal Care and Use Committees (ACUC) of the National Institute of Neurological Disorders & Stroke (NINDS) (Protocol #1531), SUNY Research Foundation at Stony Brook University (Protocol #1618), and Yale University (Protocol #2019-07895).

Additional files

Supplementary file 1. Statistical reporting table.

This file contains exact sample sizes for each group, group median, 95% confidence interval of the median, actual confidence interval, statistical test used, p-values, and test statistic for each of the reported plots in the manuscript.

elife-86581-supp1.xlsx (22.3KB, xlsx)
MDAR checklist

Data availability

Source data for the fiber photometry experiments presented in Figure 1 and supplements are provided as individual source data files. Code for fiber photometry data was previously published in Crouse et al., 2020.

References

  1. Ananth MR, Rajebhosale P, Kim R, Talmage DA, Role LW. Basal forebrain cholinergic signalling: development, connectivity and roles in cognition. Nature Reviews. Neuroscience. 2023;24:233–251. doi: 10.1038/s41583-023-00677-x. [DOI] [PMC free article] [PubMed] [Google Scholar]
  2. Azcorra M, Gaertner Z, Davidson C, He Q, Kim H, Nagappan S, Hayes CK, Ramakrishnan C, Fenno L, Kim YS, Deisseroth K, Longnecker R, Awatramani R, Dombeck DA. Unique functional responses differentially map onto genetic subtypes of dopamine neurons. Nature Neuroscience. 2023;26:1762–1774. doi: 10.1038/s41593-023-01401-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  3. Ballinger EC, Ananth M, Talmage DA, Role LW. Basal forebrain cholinergic circuits and signaling in cognition and cognitive decline. Neuron. 2016;91:1199–1218. doi: 10.1016/j.neuron.2016.09.006. [DOI] [PMC free article] [PubMed] [Google Scholar]
  4. Blanchard RJ, Blanchard DC. An Ethoexperimental Analysis of Defense, Fear, and Anxiety. University of Otago Press; 1990. [Google Scholar]
  5. Cai DJ, Aharoni D, Shuman T, Shobe J, Biane J, Song W, Wei B, Veshkini M, La-Vu M, Lou J, Flores SE, Kim I, Sano Y, Zhou M, Baumgaertel K, Lavi A, Kamata M, Tuszynski M, Mayford M, Golshani P, Silva AJ. A shared neural ensemble links distinct contextual memories encoded close in time. Nature. 2016;534:115–118. doi: 10.1038/nature17955. [DOI] [PMC free article] [PubMed] [Google Scholar]
  6. Crimmins BE, Lingawi NW, Chieng BC, Leung BK, Maren S, Laurent V. Basal forebrain cholinergic signaling in the basolateral amygdala promotes strength and durability of fear memories. Neuropsychopharmacology. 2023;48:605–614. doi: 10.1038/s41386-022-01427-w. [DOI] [PMC free article] [PubMed] [Google Scholar]
  7. Crouse RB, Kim K, Batchelor HM, Girardi EM, Kamaletdinova R, Chan J, Rajebhosale P, Pittenger ST, Role LW, Talmage DA, Jing M, Li Y, Gao XB, Mineur YS, Picciotto MR. Acetylcholine is released in the basolateral amygdala in response to predictors of reward and enhances the learning of cue-reward contingency. eLife. 2020;9:e57335. doi: 10.7554/eLife.57335. [DOI] [PMC free article] [PubMed] [Google Scholar]
  8. Ferguson SM, Phillips PEM, Roth BL, Wess J, Neumaier JF. Direct-pathway striatal neurons regulate the retention of decision-making strategies. The Journal of Neuroscience. 2013;33:11668–11676. doi: 10.1523/JNEUROSCI.4783-12.2013. [DOI] [PMC free article] [PubMed] [Google Scholar]
  9. Franklin K, Paxinos G. The Mouse Brain in Stereotaxic Coordinates. San Diego: Academic Press; 1997. [Google Scholar]
  10. Gielow MR, Zaborszky L. The input-output relationship of the cholinergic basal forebrain. Cell Reports. 2017;18:1817–1830. doi: 10.1016/j.celrep.2017.01.060. [DOI] [PMC free article] [PubMed] [Google Scholar]
  11. Hangya B, Ranade SP, Lorenc M, Kepecs A. Central cholinergic neurons are rapidly recruited by reinforcement feedback. Cell. 2015;162:1155–1168. doi: 10.1016/j.cell.2015.07.057. [DOI] [PMC free article] [PubMed] [Google Scholar]
  12. Hasan MT, Althammer F, Silva da Gouveia M, Goyon S, Eliava M, Lefevre A, Kerspern D, Schimmer J, Raftogianni A, Wahis J, Knobloch-Bollmann HS, Tang Y, Liu X, Jain A, Chavant V, Goumon Y, Weislogel J-M, Hurlemann R, Herpertz SC, Pitzer C, Darbon P, Dogbevia GK, Bertocchi I, Larkum ME, Sprengel R, Bading H, Charlet A, Grinevich V. A fear memory engram and its plasticity in the hypothalamic oxytocin system. Neuron. 2019;103:133–146. doi: 10.1016/j.neuron.2019.04.029. [DOI] [PubMed] [Google Scholar]
  13. Hersman S, Cushman J, Lemelson N, Wassum K, Lotfipour S, Fanselow MS. Optogenetic excitation of cholinergic inputs to hippocampus primes future contextual fear associations. Scientific Reports. 2017;7:2333. doi: 10.1038/s41598-017-02542-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
  14. Higley MJ, Picciotto MR. Neuromodulation by acetylcholine: examples from schizophrenia and depression. Current Opinion in Neurobiology. 2014;29:88–95. doi: 10.1016/j.conb.2014.06.004. [DOI] [PMC free article] [PubMed] [Google Scholar]
  15. Jiang L, Emmetsberger J, Talmage DA, Role LW. Type III neuregulin 1 is required for multiple forms of excitatory synaptic plasticity of mouse cortico-amygdala circuits. The Journal of Neuroscience. 2013;33:9655–9666. doi: 10.1523/JNEUROSCI.2888-12.2013. [DOI] [PMC free article] [PubMed] [Google Scholar]
  16. Jiang L, Kundu S, Lederman JD, López-Hernández GY, Ballinger EC, Wang S, Talmage DA, Role LW. Cholinergic signaling controls conditioned fear behaviors and enhances plasticity of cortical-amygdala circuits. Neuron. 2016;90:1057–1070. doi: 10.1016/j.neuron.2016.04.028. [DOI] [PMC free article] [PubMed] [Google Scholar]
  17. Jin J, Cheng J, Lee KW, Amreen B, McCabe KA, Pitcher C, Liebmann T, Greengard P, Flajolet M. Cholinergic neurons of the medial septum are crucial for sensorimotor gating. The Journal of Neuroscience. 2019;39:5234–5242. doi: 10.1523/JNEUROSCI.0950-18.2019. [DOI] [PMC free article] [PubMed] [Google Scholar]
  18. Jing M, Zhang P, Wang G, Feng J, Mesik L, Zeng J, Jiang H, Wang S, Looby JC, Guagliardo NA, Langma LW, Lu J, Zuo Y, Talmage DA, Role LW, Barrett PQ, Zhang LI, Luo M, Song Y, Zhu JJ, Li Y. A genetically encoded fluorescent acetylcholine indicator for in vitro and in vivo studies. Nature Biotechnology. 2018;36:726–737. doi: 10.1038/nbt.4184. [DOI] [PMC free article] [PubMed] [Google Scholar]
  19. Jing M, Li Y, Zeng J, Huang P, Skirzewski M, Kljakic O, Peng W, Qian T, Tan K, Zou J, Trinh S, Wu R, Zhang S, Pan S, Hires SA, Xu M, Li H, Saksida LM, Prado VF, Bussey TJ, Prado MAM, Chen L, Cheng H, Li Y. An optimized acetylcholine sensor for monitoring in vivo cholinergic activity. Nature Methods. 2020;17:1139–1146. doi: 10.1038/s41592-020-0953-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
  20. Josselyn SA, Köhler S, Frankland PW. Finding the engram. Nature Reviews. Neuroscience. 2015;16:521–534. doi: 10.1038/nrn4000. [DOI] [PubMed] [Google Scholar]
  21. Josselyn SA, Tonegawa S. Memory engrams: Recalling the past and imagining the future. Science. 2020;367:eaaw4325. doi: 10.1126/science.aaw4325. [DOI] [PMC free article] [PubMed] [Google Scholar]
  22. Kim J, Pignatelli M, Xu S, Itohara S, Tonegawa S. Antagonistic negative and positive neurons of the basolateral amygdala. Nature Neuroscience. 2016;19:1636–1646. doi: 10.1038/nn.4414. [DOI] [PMC free article] [PubMed] [Google Scholar]
  23. Kim J, Zhang X, Muralidhar S, LeBlanc SA, Tonegawa S. Basolateral to central amygdala neural circuits for appetitive behaviors. Neuron. 2017;93:1464–1479. doi: 10.1016/j.neuron.2017.02.034. [DOI] [PMC free article] [PubMed] [Google Scholar]
  24. Knox D. The role of basal forebrain cholinergic neurons in fear and extinction memory. Neurobiology of Learning and Memory. 2016;133:39–52. doi: 10.1016/j.nlm.2016.06.001. [DOI] [PMC free article] [PubMed] [Google Scholar]
  25. Krashes MJ, Koda S, Ye C, Rogan SC, Adams AC, Cusher DS, Maratos-Flier E, Roth BL, Lowell BB. Rapid, reversible activation of AgRP neurons drives feeding behavior in mice. The Journal of Clinical Investigation. 2011;121:1424–1428. doi: 10.1172/JCI46229. [DOI] [PMC free article] [PubMed] [Google Scholar]
  26. López-Hernández GY, Ananth M, Jiang L, Ballinger EC, Talmage DA, Role LW. Electrophysiological properties of basal forebrain cholinergic neurons identified by genetic and optogenetic tagging. Journal of Neurochemistry. 2017;142 Suppl 2:103–110. doi: 10.1111/jnc.14073. [DOI] [PMC free article] [PubMed] [Google Scholar]
  27. Luchicchi A, Bloem B, Viaña JNM, Mansvelder HD, Role LW. Illuminating the role of cholinergic signaling in circuits of attention and emotionally salient behaviors. Frontiers in Synaptic Neuroscience. 2014;6:24. doi: 10.3389/fnsyn.2014.00024. [DOI] [PMC free article] [PubMed] [Google Scholar]
  28. McGehee DS, Heath MJ, Gelber S, Devay P, Role LW. Nicotine enhancement of fast excitatory synaptic transmission in CNS by presynaptic receptors. Science. 1995;269:1692–1696. doi: 10.1126/science.7569895. [DOI] [PubMed] [Google Scholar]
  29. Moyer JR, Thompson LT, Disterhoft JF. Trace eyeblink conditioning increases CA1 excitability in a transient and learning-specific manner. The Journal of Neuroscience. 1996;16:5536–5546. doi: 10.1523/JNEUROSCI.16-17-05536.1996. [DOI] [PMC free article] [PubMed] [Google Scholar]
  30. Nishioka T, Hamaguchi K, Yawata S, Hikida T, Watanabe D. Chemogenetic suppression of the subthalamic nucleus induces attentional deficits and impulsive action in a five-choice serial reaction time task in mice. Frontiers in Systems Neuroscience. 2020;14:38. doi: 10.3389/fnsys.2020.00038. [DOI] [PMC free article] [PubMed] [Google Scholar]
  31. Nonaka A, Toyoda T, Miura Y, Hitora-Imamura N, Naka M, Eguchi M, Yamaguchi S, Ikegaya Y, Matsuki N, Nomura H. Synaptic plasticity associated with a memory engram in the basolateral amygdala. The Journal of Neuroscience. 2014;34:9305–9309. doi: 10.1523/JNEUROSCI.4233-13.2014. [DOI] [PMC free article] [PubMed] [Google Scholar]
  32. O’Neal TJ, Nooney MN, Thien K, Ferguson SM. Chemogenetic modulation of accumbens direct or indirect pathways bidirectionally alters reinstatement of heroin-seeking in high- but not low-risk rats. Neuropsychopharmacology. 2020;45:1251–1262. doi: 10.1038/s41386-019-0571-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  33. Pignatelli M, Ryan TJ, Roy DS, Lovett C, Smith LM, Muralidhar S, Tonegawa S. Engram cell excitability state determines the efficacy of memory retrieval. Neuron. 2019;101:274–284. doi: 10.1016/j.neuron.2018.11.029. [DOI] [PubMed] [Google Scholar]
  34. Power AE, McGaugh JL. Cholinergic activation of the basolateral amygdala regulates unlearned freezing behavior in rats. Behavioural Brain Research. 2002;134:307–315. doi: 10.1016/s0166-4328(02)00046-3. [DOI] [PubMed] [Google Scholar]
  35. Rashid AJ, Yan C, Mercaldo V, Hsiang H-LL, Park S, Cole CJ, De Cristofaro A, Yu J, Ramakrishnan C, Lee SY, Deisseroth K, Frankland PW, Josselyn SA. Competition between engrams influences fear memory formation and recall. Science. 2016;353:383–387. doi: 10.1126/science.aaf0594. [DOI] [PMC free article] [PubMed] [Google Scholar]
  36. Ray RS, Corcoran AE, Brust RD, Kim JC, Richerson GB, Nattie E, Dymecki SM. Impaired respiratory and body temperature control upon acute serotonergic neuron inhibition. Science. 2011;333:637–642. doi: 10.1126/science.1205295. [DOI] [PMC free article] [PubMed] [Google Scholar]
  37. Rickenbacher E, Perry RE, Sullivan RM, Moita MA. Freezing suppression by oxytocin in central amygdala allows alternate defensive behaviours and mother-pup interactions. eLife. 2017;6:e24080. doi: 10.7554/eLife.24080. [DOI] [PMC free article] [PubMed] [Google Scholar]
  38. Rossi J, Balthasar N, Olson D, Scott M, Berglund E, Lee CE, Choi MJ, Lauzon D, Lowell BB, Elmquist JK. Melanocortin-4 receptors expressed by cholinergic neurons regulate energy balance and glucose homeostasis. Cell Metabolism. 2011;13:195–204. doi: 10.1016/j.cmet.2011.01.010. [DOI] [PMC free article] [PubMed] [Google Scholar]
  39. Roy DS, Park YG, Kim ME, Zhang Y, Ogawa SK, DiNapoli N, Gu X, Cho JH, Choi H, Kamentsky L, Martin J, Mosto O, Aida T, Chung K, Tonegawa S. Brain-wide mapping reveals that engrams for a single memory are distributed across multiple brain regions. Nature Communications. 2022;13:1799. doi: 10.1038/s41467-022-29384-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  40. Sarter M, Bruno JP, Turchi J. Basal forebrain afferent projections modulating cortical acetylcholine, attention, and implications for neuropsychiatric disorders. Annals of the New York Academy of Sciences. 1999;877:368–382. doi: 10.1111/j.1749-6632.1999.tb09277.x. [DOI] [PubMed] [Google Scholar]
  41. Sears RM, Schiff HC, LeDoux JE. Molecular mechanisms of threat learning in the lateral nucleus of the amygdala. Progress in Molecular Biology and Translational Science. 2014;122:263–304. doi: 10.1016/B978-0-12-420170-5.00010-6. [DOI] [PubMed] [Google Scholar]
  42. Semon R. The Mneme. English Translation by Louis Simon. Macmillan Company; 1921. [Google Scholar]
  43. Stachniak TJ, Ghosh A, Sternson SM. Chemogenetic synaptic silencing of neural circuits localizes a hypothalamus→midbrain pathway for feeding behavior. Neuron. 2014;82:797–808. doi: 10.1016/j.neuron.2014.04.008. [DOI] [PMC free article] [PubMed] [Google Scholar]
  44. Tang W, Kochubey O, Kintscher M, Schneggenburger R. A VTA to basal amygdala dopamine projection contributes to signal salient somatosensory events during fear learning. The Journal of Neuroscience. 2020;40:3969–3980. doi: 10.1523/JNEUROSCI.1796-19.2020. [DOI] [PMC free article] [PubMed] [Google Scholar]
  45. Terburg D, Scheggia D, Triana Del Rio R, Klumpers F, Ciobanu AC, Morgan B, Montoya ER, Bos PA, Giobellina G, van den Burg EH, de Gelder B, Stein DJ, Stoop R, van Honk J. The basolateral amygdala is essential for rapid escape: a human and rodent study. Cell. 2018;175:723–735. doi: 10.1016/j.cell.2018.09.028. [DOI] [PMC free article] [PubMed] [Google Scholar]
  46. Tonegawa S, Liu X, Ramirez S, Redondo R. Memory engram cells have come of age. Neuron. 2015;87:918–931. doi: 10.1016/j.neuron.2015.08.002. [DOI] [PubMed] [Google Scholar]
  47. Uematsu A, Tan BZ, Ycu EA, Cuevas JS, Koivumaa J, Junyent F, Kremer EJ, Witten IB, Deisseroth K, Johansen JP. Modular organization of the brainstem noradrenaline system coordinates opposing learning states. Nature Neuroscience. 2017;20:1602–1611. doi: 10.1038/nn.4642. [DOI] [PubMed] [Google Scholar]
  48. Unal CT, Pare D, Zaborszky L. Impact of basal forebrain cholinergic inputs on basolateral amygdala neurons. The Journal of Neuroscience. 2015;35:853–863. doi: 10.1523/JNEUROSCI.2706-14.2015. [DOI] [PMC free article] [PubMed] [Google Scholar]
  49. Whitaker LR, Hope BT. Chasing the addicted engram: identifying functional alterations in Fos-expressing neuronal ensembles that mediate drug-related learned behavior. Learning & Memory. 2018;25:455–460. doi: 10.1101/lm.046698.117. [DOI] [PMC free article] [PubMed] [Google Scholar]
  50. Wilson MA, Fadel JR. Cholinergic regulation of fear learning and extinction. Journal of Neuroscience Research. 2017;95:836–852. doi: 10.1002/jnr.23840. [DOI] [PMC free article] [PubMed] [Google Scholar]
  51. Yap EL, Greenberg ME. Activity-regulated transcription: bridging the gap between neural activity and behavior. Neuron. 2018;100:330–348. doi: 10.1016/j.neuron.2018.10.013. [DOI] [PMC free article] [PubMed] [Google Scholar]
  52. Zaborszky L, Gyengesi E. The Mouse Nervous System. Academic Press; 2012. [Google Scholar]
  53. Zaborszky L, Csordas A, Mosca K, Kim J, Gielow MR, Vadasz C, Nadasdy Z. Neurons in the basal forebrain project to the cortex in a complex topographic organization that reflects corticocortical connectivity patterns: an experimental study based on retrograde tracing and 3D reconstruction. Cerebral Cortex. 2015;25:118–137. doi: 10.1093/cercor/bht210. [DOI] [PMC free article] [PubMed] [Google Scholar]
  54. Záborszky L, Gombkoto P, Varsanyi P, Gielow MR, Poe G, Role LW, Ananth M, Rajebhosale P, Talmage DA, Hasselmo ME, Dannenberg H, Minces VH, Chiba AA. Specific basal forebrain-cortical cholinergic circuits coordinate cognitive operations. The Journal of Neuroscience. 2018;38:9446–9458. doi: 10.1523/JNEUROSCI.1676-18.2018. [DOI] [PMC free article] [PubMed] [Google Scholar]
  55. Zhang W, Linden DJ. The other side of the engram: experience-driven changes in neuronal intrinsic excitability. Nature Reviews. Neuroscience. 2003;4:885–900. doi: 10.1038/nrn1248. [DOI] [PubMed] [Google Scholar]
  56. Zhang Y, Jiang YY, Shao S, Zhang C, Liu FY, Wan Y, Yi M. Inhibiting medial septal cholinergic neurons with DREADD alleviated anxiety-like behaviors in mice. Neuroscience Letters. 2017;638:139–144. doi: 10.1016/j.neulet.2016.12.010. [DOI] [PubMed] [Google Scholar]
  57. Zhong C, Du C, Hancock M, Mertz M, Talmage DA, Role LW. Presynaptic type III neuregulin 1 is required for sustained enhancement of hippocampal transmission by nicotine and for axonal targeting of alpha7 nicotinic acetylcholine receptors. The Journal of Neuroscience. 2008;28:9111–9116. doi: 10.1523/JNEUROSCI.0381-08.2008. [DOI] [PMC free article] [PubMed] [Google Scholar]
  58. Zhong C, Talmage DA, Role LW. Nicotine elicits prolonged calcium signaling along ventral hippocampal axons. PLOS ONE. 2013;8:e82719. doi: 10.1371/journal.pone.0082719. [DOI] [PMC free article] [PubMed] [Google Scholar]
  59. Zhong C, Talmage DA, Role LW. Live Imaging of Nicotine Induced Calcium Signaling and Neurotransmitter Release Along Ventral Hippocampal Axons. Journal of Visualized Experiments. 2015;1:e52730. doi: 10.3791/52730-v. [DOI] [PMC free article] [PubMed] [Google Scholar]
  60. Zhong C, Akmentin W, Du C, Role LW, Talmage DA. Axonal Type III Nrg1 Controls Glutamate Synapse Formation and GluA2 Trafficking in Hippocampal-Accumbens Connections. eNeuro. 2017;4:ENEURO.0232-16.2017. doi: 10.1523/ENEURO.0232-16.2017. [DOI] [PMC free article] [PubMed] [Google Scholar]
  61. Zhou Y, Won J, Karlsson MG, Zhou M, Rogerson T, Balaji J, Neve R, Poirazi P, Silva AJ. CREB regulates excitability and the allocation of memory to subsets of neurons in the amygdala. Nature Neuroscience. 2009;12:1438–1443. doi: 10.1038/nn.2405. [DOI] [PMC free article] [PubMed] [Google Scholar]

Editor's evaluation

Joshua Johansen 1

This important study examines the existence of a specific population of memory encoding acetylcholine neurons of the basal forebrain which regulate the amygdala for fear expression. Using a combination of techniques including genetic access to c-Fos expressing neurons, in-vivo chemogenetics, and optical detection of acetylcholine (ACh), the authors present convincing evidence that posteriorly located amygdala projecting basal forebrain cholinergic neurons participate in cue-specific threat learning and memory. This paper will be of interest to those studying circuit-level mechanisms of learning and emotion regulation.

Decision letter

Editor: Joshua Johansen1
Reviewed by: Ekaterina Likhtik2

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Basal forebrain cholinergic neurons are part of the threat memory engram" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Joshua Johansen as the Reviewing Editor and Reviewer #3, and the evaluation has been overseen by Laura Colgin as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Ekaterina Likhtik (Reviewer #2).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The authors use hM4Di to "silence" Fos-tagged neurons in the basal forebrain, but they have not validated the efficiency or the possible various effects of this reagent.

It is possible that hM4Di actually has a relatively small effect on suppressing the AP activity of neurons. Nevertheless, hM4Di might still be an effective manipulation, because it was shown to additionally reduce transmitter release at the nerve terminal (see e.g. Stachniak et al. (Sternson) 2014, Neuron). Thus, the authors should evaluate in control experiments whether hM4Di expression plus CNO actually electrically silences the AP-firing of ChAT neurons in the BF as they seem to suggest, and/or if it reduces ACh release at the terminals. For example, one experiment to test the latter would be to perfuse CNO locally in the BLA; after expressing hM4Di in the cholinergic neurons of the BF. At the very least, the assumed action of hM4Di, and the possible caveats in the interpretation of these results should be discussed in the paper. See Reviewer 1, Point 1.

2) Throughout the paper, the authors use comparisons of cell activity between groups to address questions about projection-specific and cue-specific cell activation and reactivation. However, statistical comparisons are sometimes done between biological replicates (animal groups e.g. Figure 5A), whereas a lot of them are done between technical replicates (data points/slices, e.g. Figure 2B, 5B, 7B). Adding statistics that compare biological replicates would help increase confidence in the results.

3) To demonstrate engram-like specificity, in figure 4C the authors show fold change in cholinergic reactivation in low and high responders (animals that show low and high defensive freezing upon cue presentation) as normalized by cell activity while sitting in the home cage. However, the authors also collected a better control for this comparison, which is shown in figure S4, where the animals were exposed to an unconditioned tone cue. Comparing fold change to this tone-alone condition would provide stronger evidence for the authors' point, as this would directly compare the specificity of cholinergic reactivation to a conditioned vs an unconditioned cue. A discussion of the same comparison is relevant for figure 2 (and is shown in figure S4) but is not mentioned in the text.

4) The significant correlation between cue-evoked percent change in defensive freezing from pretone and fold change in cholinergic cell activity relative to the home cage that is shown in figure 4D is somewhat confusing. Is the correlation considering all the points shown (high and low responders as depicted by black and grey points)? It's first reported as one correlation but then is discussed as two populations that have different results. Further, is the average amount of reactivation for the home-cage controls used here the same denominator for each reported animal? Similarly to the point above, a correlation looking at fold change from tone-alone would also be helpful to determine the degree to which cholinergic reactivation is specific to threat-association learning versus the more general attentional component that this system is known for.

5) Some important controls are missing from the 'engram' causal manipulation study (Figure 3) which makes it difficult to definitively support the authors' contention that the captured cells are in fact 'engram' cells. The established definition of engram cells refers to those which store memory and are activated during learning and memory expression. While the authors do show conditioning-induced changes in CS activation in NBM-ACh cells, they don't have the necessary controls for their causal manipulations to rule out other non-engram factors which could account for their results. For example, it is possible that the captured cells are simply activated by auditory cues (even prior to learning, the immuno approach in Figure 2 may not be sensitive enough to detect this) and that they are transmitting CS-activated ACh signals to BLA which are necessary for memory expression. It is also possible that inhibiting a small number of NBM-ACh neurons (not restricted to 'engram' cells) reduces memory expression (early engram studies by the Josselyn lab controlled for this). The authors could address this with further experiments such as including 1) a control group in which cells are captured (off-dox) during a habituation session with CS presentation prior to on-dox learning and 2) using a dilute virus to manipulate a similar number of cells as in the engram group. Alternatively, they could present a more careful discussion of this issue and potentially avoid the use of the term 'engram'.

6) The authors suggest that the ACh 'engram' cells participate in the expression of learned defensive responses through projections to the BLA. However, they have not demonstrated this definitively. To do this would require testing the effects of the inactivation of ACh-engram neurons which project to BLA on learned freezing responses. This may be technically challenging, but one approach could be to use the ADCD technique they've developed and inactivate terminals of captured cells in the BLA (chemogenetically or optogenetically). At a minimum, the authors should examine whether retrogradely labeled ACh cells projecting to BLA are Fos activated during a recall session (compared with appropriate controls).

7) For most experiments, significantly more raw data should be shown (e.g. raw example traces for GRAB-ACh3.0), and also brain section images for almost all experiments. Significantly more raw data should also be shown in the Main Figures. Please see Reviewer # 1, "Specific points relating to data presentation and statistical analysis", for details and specific examples."

Reviewer #1 (Recommendations for the authors):

Specific points relating to data presentation and statistical analysis:

1) The reporting of statistical data (N's, p- values, and the test used) is a bit cryptic. In general, the authors seem to report these values in the Figure legends. This, however, has the advantage that the reader has to constantly jump between the results text and legends text if he/she wishes to check the values for a given experiment. I would therefore suggest that the authors report all values of statistical tests in the results text (N's, p- values, and the test used).

Also, in some instances, N-values seem to be missing from the legend texts. Please check all legends / Results text accordingly.

2) Figure 1:

– Traces of GRAB-ACh3.0 should be shown more extensively. In Figure 1B, only the onset to the CS is shown – the response along the entire tone (30 s) should also be shown.

– Furthermore, the trace in Figure 1B seems to be an average over all mice. In addition, traces from individual mice should be shown, to allow an estimation of the S/N ratio of these measurements.

– In the quantification (Figure 1B, right), it would be easily possible to show the individual data points (there are 11 mice). I did not understand what the leftmost quantification tells us ("Pre-tone"). Please explain how this was calculated, and inasmuch it is relevant.

– What is the unit "Z%deltaF/F0"? In my view, this should either be Z-scored or else, deltaF/F0 – please clarify.

3) Figure 2:

– the authors should extensively show images of brain sections. It is curious to see that the entire paper doesn't show a single image at low magnification, which shows the investigated basal forebrain areas, and stained neurons in there, in their anatomical context. This would allow the reader to put the images into context, and to roughly validate whether the analyses by the authors (e.g. counted cells) make sense. For example, for Figure 2A, please show an additional low mag image of the corresponding brain area, that shows this brain area (NBM/SIp) in the wider context. The actual image can then be shown as an insert.

– In the plot of Figure 2B, a tick at the value of 0 is missing – please add this (this also applies to Figure 2C, 4B, 4C, 5A, right, 5B, right, 5C, right, 7A, 7D, S6B, S6C, S7A-C – please add all ticks at zero). Adding the tick at zero will help the reader eyeball the values of the data points.

– In some plots, the y-axis doesn't start at 0, but at an artificially chosen higher value. This increases "graphically" the effect size; however, it doesn't allow the reader to look at the data in the most unbiased way (indeed, some Journals do not allow this practice) (see Figure 7C, Figure S6A, right). Please newly plot the data from zero.

– The authors did not report in the Statistics part, which criteria they used to add "stars" onto the plots. For example, Figure 2B has *** (3 stars), but the reported p-value is p = 0.0017. According to the criteria used by many authors, p < 0.05 (1 star), p < 0.01 (two stars), and p < 0.001 (three stars). Thus, Figure 2B should have 2 stars, but not 3. Please explain the criteria you used, and please validate all star symbols according to the criteria you used.

– What is the unit of the cell numbers plotted in Figure 2B, C? Were these indeed absolute numbers as the axis labelling seems to suggest, or was there a normalization for the analyzed area? Please comment.

4) Figure 3:

– While the freezing data is shown at a quite high time resolution for the training day (Figure 3B), for recall (Figure 3C), only excessively binned data are shown. It would be advantageous for the reader to see the data at the same time resolution, in a quasi-continuous trace, to evaluate the effects also of silencing.

5) Figure 4:

– There are several issues with the definition of high and low freezers, and with the analysis in this Figure.

– Hig+low freezers: The Methods state "High responders were defined as those mice that exhibited at least a 10 percentage point increase in % time spent freezing in the bin during the tone from the bin immediately pre-tone (e.g. Pre-tone freezing 10% to tone-induced freezing of {greater than or equal to}20%)." However, the Results text speaks about a "3-fold increase" (l. 235) – these are different criteria, please explain.

– Related to the above, the y-axis labeling in Figure 4B is unclear – is this the difference in freezing (as the Methods would suggest) or a fold-change (like the Results would suggest)? Also, the 3-star symbol for p = 0.005 seems not justified, should be 2 stars according to often-used criteria (see also above).

– l. 238 states: "There was no difference in mCherry expression in 'Low Responders' compared to the home cage group (fold change ~ 1),"

– this should be tested statistically.

– I don't see that the freezing data in Figure 4D (y-axis) shows a "positive correlation" with the number of Fos-labelled cells (x-axis). Indeed, looking at the high responders only, the data falls very close to a flat line. Please explain how the correlation was analyzed – this was probably done by fitting all data, i.e. also including the "low" responders (?). The obvious absence of a correlation for the "high responders" alone (black data points) is worrisome.

– In Figure 4D, the location of the data points (y-axis, freezing; about -10 – 0%) does not seem to correspond to what is plotted in Figure 4B (~ 2%; please explain).

– Figure 4F, and Figure S5: The result about the location of labelled cells along the a-p axis of the brain should be documented with a series of original images, both at low and high mag (see general comments, above).

6) Figure 5:

– The location of back-labelled cells should again be amply documented with several low and high mag images from all analyzed areas (hDB, VP/SIa, NBM/SIp).

– 3 stars in Figure 5B might be 2 stars.

7) Figure S4:

– was the data in Figure S4C a further analysis of data in Figure 2B? If yes, why are the animal numbers seemingly different between the two datasets?

8) Figure S5:

– The legend doesn't state the N-numbers.

– The box plot should be explained (meaning of the boxes). Even better, individual data points should be shown instead, if this is feasible.

9) Figure S6:

– For the data in A, again example images at low and high magnification should be shown …

– For the images in B, please indicate lateral-medial, and dorsal-ventral directions in these images (assuming that these are coronal sections?).

– The data in C again suggest that there might also be an effect on contextual fear learning (see also above).

10) Figure S7

– The text describing the data in A states "did not freeze to the tone" (l. 305). This is not a scientifically warranted statement, since there seems to be residual freezing. Rather, the data should be tested statistically.

– The N-numbers for the data in A, and B should be stated in the legends.

– For these and possibly similar experiments, the authors need to show in detailed histological images, for each individual mouse, that the expression of hM4Di has occurred in the expected set of cholinergic BF neurons.

11) Figure S8:

– I wonder why the firing rate data of the "Training" group in panel H is extremely small (around 5 Hz or less), and smaller than in all other groups. Is there another type of modulation of intrinsic excitability, a down-modulation, after training? Please explain.

12) Figure S9:

– In panel C, could you show the freezing data at a higher time resolution, e.g. the 30s bins also used in Figure 3B?

Reviewer #2 (Recommendations for the authors):

– For the flow of the story, I recommend switching the order of presentation where the lack of effect on behavior after silencing the CeA-projecting NBM/SIp neurons (lines 286-301) comes after the data showing the effect of silencing the BLA-projecting NBM/SIp.

– I would also recommend moving figures S7A and S7B into the main figures as they provide important evidence for the author's point.

– Given that the authors don't see an effect of clozapine on memory when inactivating during the consolidation period (Figure S7C), it might be worthwhile to discuss why this might be the case, and the effects of cholinergic input starting cascades of neural activation in the BLA during training.

– Overall, there are missing Bregmas and mapping of placements in several places (e.g. Figure 1A – Bregma and mapping of placements/spread of the virus across animals, it appears that the recordings are relatively anterior in the BLA; Bregma in Figures 5, S3C, 6SB, S6C – not clear how much of a CeA there is at this location, 7D; a zoomed out section of where the example data in figure 7C is coming from).

– It would also be helpful to indicate the area where the signal fluorescence is being calculated in traces shown in figures 1 and S1.

Reviewer #3 (Recommendations for the authors):

1) For the 'engram' data in Figure 2, the logic of the double labeling is a bit confusing as described and it's not clear, from a simple reading of the paper, how the authors achieve double labeling of encoding (training activated) and recall activated cells. After digging around in the Methods, Figs/legends I was able to figure it out, but it was difficult initially. The authors should state this in a simple way in the main text along with a description of how this works at the level of transgene expression at different time points.

2) Some important controls are missing from the 'engram' causal manipulation study (Figure 3) which makes it difficult to definitively support the authors' contention that the captured cells are in fact 'engram' cells. The established definition of engram cells refers to those which store memory and are activated during learning and memory expression. While the authors do show conditioning-induced changes in CS activation in NBM-ACh cells, they don't have the necessary controls for their causal manipulations to rule out other non-engram factors which could account for their results. For example, it is possible that the captured cells are simply activated by auditory cues (even prior to learning, the immuno approach in Figure 2 may not be sensitive enough to detect this) and that they are transmitting CS-activated ACh signals to BLA which are necessary for memory expression. It is also possible that inhibiting a small number of NBM-ACh neurons (not restricted to 'engram' cells) reduces memory expression (early engram studies by the Josselyn lab controlled for this). The authors could address this with further experiments such as including 1) a control group in which cells are captured (off-dox) during a habituation session with CS presentation prior to on-dox learning and 2) using a dilute virus to manipulate a similar number of cells as in the engram group. Alternatively, they could present a more careful discussion of this issue and potentially avoid the use of the term 'engram'.

3) The authors suggest that the ACh 'engram' cells participate in the expression of learned defensive responses through projections to the BLA. However, they have not shown this definitively. To do this would require testing the effects of the inactivation of ACh-engram neurons which project to BLA on learned freezing responses. This may be technically challenging, but one approach could be to use the ADCD technique they've developed and inactivate terminals of captured cells in the BLA (chemogenetically or optogenetically). At a minimum, the authors should examine whether retrogradely labeled ACh cells projecting to BLA are Fos activated during a recall session (compared with appropriate controls).

4) In Figure 5, the authors report a reduction in BLA Fos-activated neurons with chemogenetic inhibition of BLA-projecting ACh neurons. However, this experiment requires a control group without training (or, ideally, with unpaired CS/US training) to demonstrate that the inhibition is affect learning-induced activation of BLA neurons. While others have reported this using other paradigms previously, it is still important to show here using their own preparation.

[Editors’ note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Basal forebrain cholinergic neurons are part of the threat memory engram" for further consideration by eLife. Your revised article has been evaluated by Laura Colgin (Senior Editor) and Josh Johansen (Reviewing Editor).

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

We thank the authors for their response to the reviewers' previous suggestions. Many of the findings have been strengthened through the addition of new experiments and analyses. However, there are remaining concerns about the statistical approaches used and low population size in some experiments which should be addressed. Related to sample sizes, power analyses or clear demonstrations of existing standards in the field should be used to justify animal numbers. If necessary, further animals should be added to experimental groups (e.g. the differential role of NBM/SIp and VP/SIa in learned vs. innate threats) or, as suggested by reviewer 1, data with low n could be removed if it does not seriously undermine the conceptual advance of the paper. In addition, the strong focus on the activated basal forebrain cholinergic cells as being part of the threat memory engram is not well supported. The authors should moderate the claims related to this aspect of the paper.

Furthermore, the title and/or abstract should provide a clear indication of the biological system under investigation (i.e., species name or broader taxonomic group, if appropriate). Please revise your title and/or abstract with this advice in mind.

Reviewer #1 (Recommendations for the authors):

For this first revision, improvements in the data display have been implemented, at least concerning the traces for GRAB-ACh measurements. In addition, new data has been added; thus, the ms has been reorganized quite substantially. Nevertheless, there are still a lot of questions regarding the statistical analyses, and the anatomical analysis. In one instance, it seems that a control data set was used twice (Figure 5C and Figure 5SD; see Specific point 16 below).

Since my list of points is long, the question arises what might be the best strategy to go forward with the paper? Most of my points are about data display, statistical analyses, and similar, and it should be feasible to address these points. On the other hand, I also feel that especially the anatomical analysis relating to Figures 5 and 6 is still not clear and that some of the points the authors want to make related to these Figures, especially about the a-p distribution of (re-)activated Ach neurons (differential role of NMB/SIp versus VPSIa neurons in learned fear), differential role of these pools in "high" and "low" responders, and differential role of the same in learned versus innate fear, are unclear, and might go beyond the scope of the present paper. In my view, the authors should therefore consider the possibility to remove some of the controversial data, and related claims (see my detailed comments – these data are also often plagued by low n-numbers and incomplete anatomical analyses), and concentrate on the remaining main findings of the paper.

Specific points (in order of appearance)

1) l. 156/157 "when it was not explicitly paired with a shock" (data in Figure 1S3D).

– Please explain in the legend the exact temporal sequence between shock(s) and tone(s) used for this "non-explicitly paired" protocol (or does this refer to the "shock-only" protocol?).

2) Figure 1 S3C, D: The legend to this Figure states that here, the results from a "shock-only" training are displayed for the subsequent "recall" session.

First, the authors should also show the freezing response for the standard group undergoing standard fear learning (3x tone-shock pairings); this is important for a comparison to the "shock-only" training group.

Second, the Results text states that Figure 1S3D shows the ACh response when "it was not explicitly paired with a shock".

– Is this the same as the "shock-only" protocol described in the legend? Please clarify, and make changes in the manuscript as appropriate.

3) Figure 1 S4 was not called out from the Results text; this should be changed i.e. what was done in this Figure, should be (briefly) mentioned in the Results.

4) The use of star symbols in Figure 2C, in conjunction with what is stated in the corresponding Results text, is misleading and tends to over-interpret the statistical significance of this data set.

Specifically, l. 202 states: "Significantly more double-positive cholinergic neurons were seen following the complete (tone + shock followed by tone recall) associative threat learning paradigm compared to mice that underwent session 1 without shocks (Figure 2C, p=0.0249)."

– This is misleading. The Figure legend and statistics table state that the value "p = 0 0.0249" refers to the "overall" significance of the Kruskal-Wallis (KW) test used here, but NOT to the specific comparison cited in the above sentence. For that specific comparison, the p-value was higher ("p=0.0464, Dunn's corrected", see Figure legend).

– In response, the authors should spell out the "overall" significance of this dataset (effect of group; p = 0.0249; KW) in the Results section, followed by the pertinent p-value for important specific comparisons (in this case, p = 0.046 for the comparison between the "tone-tone" and tone+shock – tone" groups). The n-values for these pertinent comparisons should also be given in the Results text.

Furthermore, in this dataset (Figure 2C), the post-hoc test for the comparisons between the home cage (HC) and standard fear conditioning group is NOT significant (p = 0.72). This value should be indicated by a "bracket" symbol in Figure 2C, and the meaning of this non-significance should be discussed in the manuscript.

5) The p-values, n-values, and the statistical tests applied, should be reported together either in the figure legends or in the Results text (please choose one section to use for reporting these statistics, in addition to the statistical table) throughout the entire manuscript (see point 4, and see also my previous point "Data Recommendation 1").

6) Figure 2D:

l. 210 states:

"All three conditions showed few reactivated neurons and no differences between groups (p=0.9471)"

– The finding that all these 3 control groups are not significantly different from each other, is not surprising and does not go to the point of the control experiment. Rather, it should be tested statistically whether the number of re-activated neurons in these three conditions is significantly different from the one after standard fear conditioning (the data in Figure 2C, right column).

7) I think the title of the Sub-header:

"Reactivation of cholinergic neurons activated by training is required for learned behavioral responses" (l. 214) is wrong, it should read "Activation" instead (and then replace "activated by" with "during").

To my understanding, the chemogenetic silencing experiments in Figure 3 test the role of activated cholinergic neurons but do not suppress the "re-activation" of neurons (during recall…), since ADCD expression of hM4Di is driven during the training session.

8) In Figure 3 S1B, the average trace of time-resolved freezing is shown, for the same data set as the one in Figure 3D. However, the number of mice included in the time-resolved display was lower (n = 5 and 6) than in the display of Figure 3D (n = 7 and 8). This must be homogenized, i.e. all data used for the statistical analysis in Figure 3D (n = 8 and 7 mice) should also be shown in a time-resolved fashion for Figure 3S1B.

9) l. 273 states:

" We injected cav.hM4DiBLA or sham control mice with clozapine (CLZ) 10 min prior to initiating cue-conditioned threat learning"

– Please add here that these were IP (?) injections, for a quick reference for the reader. Indeed, this information should possibly also be added to the Methods (I did not find it upon a rapid check).

10) The authors should control for the possibility of leaky expression of hM4Di, following presumed Cre-dependent retrograde expression from a CAV2 virus after injections in the BLA (Figure 5). Because, as I assume, CLZ was not injected into the basal forebrain (BF) but rather IP (see point 9), a leaky expression of hM4Di beyond the cholinergic BF has the potential to strongly change the interpretation of this experiment. Thus, the authors have to show, first, whether hM4Di-mCherry is expressed at the injection site in the BLA itself (this would lead to a direct inhibition of BLA neurons). This should be done by showing example images of hM4di-mCherry fluorescence on the level of the BLA, and by analyzing the number of hM4Di-mCherry+ neurons in the BLA over the entire dataset. Second, the authors should show extensive raw example images of the BF at different a-p locations (and of further brain areas in a zoomed-out fashion), in an attempt to exclude the possibility that hM4Di might also be expressed in non-cholinergic neurons (ChAT could be visualized by IHC). Only a non-leaky expression of hM4Di-mCherry, essentially limited to ChAT-positive neurons in the BF, would support the conclusions of the authors.

11) On l. 298, the results of a control experiment is reported, which is not documented ("data not shown"). Some journals do not allow the statement "data not shown". The authors should document the control data set, or else remove this statement from the manuscript. See also l. 441, another "data not shown", which should be similarly documented or else removed.

12) Data in Figure 5C.

When checking the p-values from the legends/statistics table, I was puzzled by the fact the comparison of the two "grey" data sets has a reported p-value of < 0.0001, although many of the data points between the two groups overlap. On the other hand, the comparison between the "sham" group and hM4Di group at the "pre-tone" time had a p-value of 0.096, although the overlap of the data points is only partial. Please comment, and validate the underlying tests.

13) Related to point 12, there is a strong tendency that freezing before the tone is lower in the hM4Di group as compared to the "sham" group. Although the authors reported a p-value of 0.096 (and therefore, "not significant") for this comparison (Figure 5C), this p-value should nevertheless be verified. Because the trend is seen in this data set, as well as in several other ones (Figure 5B and Figure 3D), the authors should discuss the possibility that chemogenetic suppression of ACh release might, in addition to the effects studied here, also reduce contextual fear learning (see also my previous "Specific comment 2").

14) Still related to the data in Figure 5C:

The time-resolved display of tone-evoked freezing during the recall day should be shown both for this data set (Figure 5C), and for the data set in Figure 5B, similarly as was done above for the data in Figure 3D (see panel Figure 3 S1B). This would also be helpful to gain insight into the p-values resulting from ANOVA / KW testing (see point 12), which were probably performed on similar time-resolved data sets (because the effect of "time" was always tested).

15) l. 299 states:

"Differences in recall-induced Fos expression between sham and cav.hM4DiBLA mice were maximal in rostral portions of the BLA (between bregma -0.8mm to -1.4mm) (Figure 5-Supplement 1B)"

– First, the analysis in Figure 5 S1B seems to have been done in only N = 1 mouse, which is problematic (please confirm the n-number, and state it in the Results text and in the legend to that Figure).

Second, I actually fail to see an a-p gradient as claimed by the authors (higher in anterior than posterior). If one looks at the number of Fos+ cells in the hM4Di group (red data points), then it seems that hM4Di suppresses the number of Fos+ cells, in relative terms, stronger at more posterior levels (e.g. reduction by > 80% at -1.2 mm) than at anterior levels (e.g. reduction by < 40% at -0.8 mm). Wouldn't this suggest that opposite from the author's interpretation, MORE cholinergic neurons are specifically (in a hM4Di-dependent manner) activated posteriorly than anteriorly? Please comment.

16) Visual inspection of the data in Figure 5 S1D strongly suggests that the data in the control ("sham") group (n = 9, grey data) are actually the same as the data for the sham group in Figure 5C. In Figure 5C, retrograde expression from the BLA was studied, whereas in Figure 5 S1D, the effects of retrograde expression of hM4Di from the CeA were studied. If the experiment was originally set up for 3 groups of mice (thus, sham control – CAV2-hM4Di in BLA – CAV2-hM4Di in CeA), then this should be reported as such, the data should be shown in one panel, and the statistical comparisons should be done accordingly (i.e. the "multiple comparison" correction would need to be done for 3 groups). If, on the other hand, the CeA experiments were done without proper controls, and were simply compared to the control group from another experiment, then the CeA data should be removed. Please comment, and make appropriate changes in the manuscript.

17) Figure 5 S2 claims that neurons in the "NBM/SIp" versus neurons in the "VP/SIa" have been targeted selectively with AAV9:DIO:hM4Di-mCherry. This should be documented extensively, by showing expression of hM4Di-mCherry throughout several brain slices in a de-zoomed fashion for an example mouse, and by showing average results from the quantification of the localizations in a suitable manner (for example, by plotting cell density/number along the a-p axis), for all investigated mice. Otherwise, it remains unclear whether the two areas have indeed been targeted differentially.

18) In Figure 5 S2A, the "triple star" symbol comparing the two control data points might be removed. It is not surprising that the tone recall causes increased freezing in the "sham" group.

19) The data set reported in Figure 5 S2B looks underpowered (only n = 3 and 4 mice). I highly doubt that with such a small sample size, significant conclusions can be drawn in a behavioral experiment. For some reason, the data in this experiment seem to scatter much less than in comparable experiments even with larger n's (compare to the neighboring panel A, bottom) – please comment. Also, the use of star symbols (3 stars?), and the use of the bracket symbols is misleading (the brackets suggest that the data was pooled amongst all "pre-tone" and "recall tone" groups; the rationale of this analysis approach remained cryptic). Importantly, there was no group effect in the ANOVA (p = 0.1574 was reported). Thus, no further post-hoc group comparisons should be done. Overall, I think the authors need to first, show that they can indeed differentially target NBM/SIp versus VP/SIa (see point 17). If this can be shown, then they would have to obtain significantly more data in both the control and effect groups for the VP/SIa targeting, to drive the point home that targeting these two areas has differential behavioral effects.

On the other hand, the authors might simply remove the data in Figure 5 S2A and B from the manuscript, and remove the claim about the differential effects of the NBM/SIp and VP/SIa. These claims are currently little supported by experimental evidence, and I feel they are not important for the overall conclusions of the paper.

20) Figure 6A shows time-resolved freezing responses on the recall day for n = 9 control mice, and n = 9 effect mice. However, it remained unclear from which of the previous experimental groups these data were taken; whether a sub-selection was made, and whether animals were pooled from various manipulations. This should be explained in detail in the Results text and legend.

21) The "fold-change" analysis of expression analyzed in Figure 6D remained unclear to me (sorry if I didn't get it). Please explain this briefly in the Results text, and make sure it is explained in more detail in the Methods. Also, please give details about the viral construct used in Figure 6C (currently it only states "ADCD").

22) The picture shown in Figure 6E is simply a scheme of the idea that neurons activated in "high responders", in your view, should be localized more anteriorly, but it is in no way a documentation of an anatomical observation. This scheme should be replaced by experimental evidence.

23) The data that might relate to the above issue (localization of ChAT+, activated neurons along the a-p axis in "high" and "low" responders), might be contained in Figure 6 S2. However, the data seems to be only partially analyzed, and it is confusing to fully understand how the plot in panel C was derived. First, it was unclear from which group of mice the images in Figure 6 S2A were taken. Please explain in the legend from which experiment these images were taken, which mouse genotype was used (Fos-tTA/shGFP with or without ChAT-Cre?), and which viral constructs were injected in these mice.

For panel C, I assume that the red data point represents an analysis of the images in panel B (correct?), but I am unsure from which data set the black data points were taken (# of recall-activated BFCNs (ADCD)") – please mention this in detail in the legend of the Figure. Were the black data points related to the images in panel A? were they taken from "high" or "low" responder mice?

24) The images shown in Figure 6 S2C, should be shown in a de-zoomed fashion for several representative brain images on the level of the BF, and beyond. Such a documentation of the data would represent the control for "leaky" expression mentioned above (best after co-staining against ChAT; see point 10).

25) The 2-star symbol in Figure 7D above the comparison "HC" – Fos+" should probably be 3-stars, not 2-stars (p = 0.0006).

26) The authors should use the less "pushy" and more conventional criteria for using star symbols for statistical significance (i.e. 1 star, p < 0.05; 2 stars, < 0.01; 3 stars < 0.001), and omit the "4 star" symbol and definition throughout. The 1-3 star criteria should be stated in the Methods (I didn't find it). The paper in its current version has only a single "4 star" symbol (Figure 3D), which should be changed to "3 star" according to the above definition, or removed (see also point 18) (see also my previous point "Data Recommendations 4").

27) The bracket symbol in Figure 7F is probably wrong. On the left side, it should only point to the "HC" group, but not to both the "HC" and "Fos-" group as it stands now.

28) The authors claim that "latency to fire an AP" and AHP might significantly change in re-activated ChAT neurons (Figure 7 S1E, F). However, these effects are somewhat more subtle than the other effects shown in the main Figure 7. For example, the effect of "latency" seems to be driven merely by some outlier neurons in the home cage group with very long latencies (the authors would need to show example recordings), whereas the data from the Fos- and Fos+ groups virtually overlay and the latter comparison is the more pertinent one. Thus, I think it is more cautious to drop this claim.

29) l. 427 states:

"Next, we quantified the number of cholinergic neurons expressing Fos (Fos shGFP+) after saline or predator odor exposure"

– Please explain how long the mice were exposed to the odors.

30) In Figure 8, the authors seem to analyze the results of a triple-labeling anatomical experiment; i.e. there are neurons that are Fos+; neurons that are back-labelled by fast Blue; neurons that are ChAT+, and all combinations thereof. Thus, to give the reader a better view of the various neuron pools, the authors should report the numbers of all pools of neurons (i.e. the "green", "red", "blue" neurons and all combinations thereof). This can be either done numerically or in the form of a Venn diagram. As it stands, only part of the information is given (Figure 8C, right; which seems to correspond to the percentage of back-labelled, ChAT+ neurons which were also Fos+). However, this is only part of a bigger picture. For example, the images in Figure 8C, and in Figure 8 S1A (middle, right) show also many non-ChAT+ neurons, that were re-activated in the VP/SIa area.

31) What does the dashed line label "Home cage" in Figure 8C, right panel, mean, and how was it derived?

32) l. 442:

" Nearly the entire subset of BLA-projecting VP/SIa cholinergic neurons were also Fos-shGFP+ (Figure 8C right)."

– Another non-scientific statement. Rather, the average plus/minus S.D. of the quantification should be reported. In addition, information about all labelled neuronal populations should be given (see above, point 30).

33) The retrograde expression experiment of hM4Di from the BLA in Figure 8D is expected to drive hM4Di expression in all cholinergic BLA-projectors, not selectively in the VP/SIa. This caveat in the interpretation should be mentioned more carefully when the results of this section are summarized (lines 451 – 454).

Reviewer #2 (Recommendations for the authors):

I thank the authors for addressing my concerns, and I'm happy to recommend the manuscript for publication in its current form. There are a few figures that would benefit from additional labeling (outlined below).

Please add anatomical references for basal forebrain in figures F1S4, F2S1a, F4c, F5S1a, F5S1b , F6S2a.

Reviewer #3 (Recommendations for the authors):

The authors have addressed many of my previous concerns and now nicely dissociate the role of NBM/SIp and VP/SIp ACh neurons in learned threat memory processing and show that NBM/SIp ACh cells are activated by memory recall and necessary for threat memory expression. However, I still have a lingering concern with their strong emphasis on the role of NBM/SIp ACh neuron (and their projections to BLA) as being memory 'engram' cells. While some of their data are consistent with this idea, they do not test key criteria for defining engram neurons, and the new data supporting a role of NBM/SIp ACh cells which project to the BLA in the threat memory 'engram' is only correlative.

Specifically, while there is evidence in support of the NBM/SIp ACh cell population (independent of where they project) in memory engrams (activated by memory cues, necessary for memory expression). There are several criteria for defining engram neurons (see Josselyn and Tonegawa Science 2020) which haven't been met including sufficiency (stimulation should produce memory in the absence of recall cues), mimicry (creation of an artificial engram) and long-lasting cellular or synaptic changes in engram neurons (here they only see excitability differences at 2.5 hr, not at long term memory timepoints). It is possible (even likely) that activity in these neurons modulates/facilitates threat memory formation and expression and are not storing the memory trace themselves.

There is even less evidence that NBM/SIp ACh cells which project to BLA being engram cells. The authors report an increase in ACh release in BLA (which isn't necessarily coming only from NBM/SIp) and added new data showing that NBM/SIp ACh cells projecting to BLA are fos activated.

While there is data supporting a potential role for these cells in engram processing, it is inappropriate to so strongly insist (in the title, abstract, and throughout) that NBM/SIp ACh neurons, and particularly those which project to BLA, are engram cells. The authors should moderate this aspect of the paper and perhaps broaden the discussion of the potential role of these cells in threat memory processing as they've already nicely done in the 'Cholinergic Modulation of Associative Threat Learning' section of the Discussion.

eLife. 2024 Feb 16;13:e86581. doi: 10.7554/eLife.86581.sa2

Author response


Essential revisions:

1) The authors use hM4Di to "silence" Fos-tagged neurons in the basal forebrain, but they have not validated the efficiency or the possible various effects of this reagent.

It is possible that hM4Di actually has a relatively small effect on suppressing the AP activity of neurons. Nevertheless, hM4Di might still be an effective manipulation, because it was shown to additionally reduce transmitter release at the nerve terminal (see e.g. Stachniak et al. (Sternson) 2014, Neuron). Thus, the authors should evaluate in control experiments whether hM4Di expression plus CNO actually electrically silences the AP-firing of ChAT neurons in the BF as they seem to suggest, and/or if it reduces ACh release at the terminals. For example, one experiment to test the latter would be to perfuse CNO locally in the BLA; after expressing hM4Di in the cholinergic neurons of the BF. At the very least, the assumed action of hM4Di, and the possible caveats in the interpretation of these results should be discussed in the paper. See Reviewer 1, Point 1.

We clearly demonstrate that activation of hM4Di expressed in basal forebrain cholinergic neurons with clozapine significantly alters behavior and/or activation of neurons in projection targets of BFCNs (see Figure 3, Figure 3-supplement 1, Figure 5, Figure 5 Supplement 1, Figure 5-Supplement 2, Figure 6-Supplement 1 and Figure 8). As discussed on Page 17 Lines 488-493 of the revised manuscript, whether the effects we see on behavior result from inhibition of action potential firing per se, or from decreasing ACh release from terminals (or both) is not clear. However, the effect on behavior is striking and consistent with prior demonstrations that activating hM4Di (or hM3Dq) in cholinergic neurons changes membrane potential, blocks action potential firing and reduces firing rate (Zhang et al. 2017 and Jin et al. 2019) demonstrating DREADD-mediated excitation and inhibition of basal forebrain cholinergic neurons can alter electrical activity and behavior.

2) Throughout the paper, the authors use comparisons of cell activity between groups to address questions about projection-specific and cue-specific cell activation and reactivation. However, statistical comparisons are sometimes done between biological replicates (animal groups e.g. Figure 5A), whereas a lot of them are done between technical replicates (data points/slices, e.g. Figure 2B, 5B, 7B). Adding statistics that compare biological replicates would help increase confidence in the results.

We have replotted our data as scatter plots showing values from individual animals

(biological replicates) in New Figures 1 – 8 and in Figure 1-Supplements1, 2, and 3,

Figure 5-Supplements 1 and 2, Figure 6-Supplements 1 and 2, Figure 7-Supplement 1,

and Figure 8-Supplement 1. All statistical analyses have been calculated using biological replicates. To further clarify these analyses, we have included a statistical reporting table that details sample size, group medians, confidence intervals, statistical tests used (and when appropriate, post-hoc corrections). This can be found as a separate file uploaded with the manuscript associated documents.

3) To demonstrate engram-like specificity, in figure 4C the authors show fold change in cholinergic reactivation in low and high responders (animals that show low and high defensive freezing upon cue presentation) as normalized by cell activity while sitting in the home cage. However, the authors also collected a better control for this comparison, which is shown in figure S4, where the animals were exposed to an unconditioned tone cue. Comparing fold change to this tone-alone condition would provide stronger evidence for the authors' point, as this would directly compare the specificity of cholinergic reactivation to a conditioned vs an unconditioned cue. A discussion of the same comparison is relevant for figure 2 (and is shown in figure S4) but is not mentioned in the text.

We compared cholinergic neuron response to tone under three conditions: (1) in animals that were exposed to three tone-shock pairings, and then to a single tone; (2) in animals exposed to three tones (no shock) and subsequently to a single tone; and (3) in animals exposed to three shocks (no tone) and subsequently to a single tone. These responses were compared to responses from animals that remained in the home cage. There was no significant increase in ACh release in the BLA in response to tone presentation unless the tone had previously been paired with foot shock (condition 1) (new Figure 1 and Figure 1Supplements 1-4), and no significant increase in tone associated re- activation of cholinergic neurons (using IEG expression as a read-out) relative to home cage unless the tone had been previously paired with foot shock (new Figures 2, 4), and no significant difference between responses to conditions 2, 3, or 4 (Figure 2C- home cage-home cage vs. tone-tone comparison p=0.5012). Based on these analyses, for the subsequent experiments, we have used the home cage conditions as the control group for comparison.

4) The significant correlation between cue-evoked percent change in defensive freezing from pretone and fold change in cholinergic cell activity relative to the home cage that is shown in figure 4D is somewhat confusing. Is the correlation considering all the points shown (high and low responders as depicted by black and grey points)? It's first reported as one correlation but then is discussed as two populations that have different results. Further, is the average amount of reactivation for the home-cage controls used here the same denominator for each reported animal? Similarly to the point above, a correlation looking at fold change from tone-alone would also be helpful to determine the degree to which cholinergic reactivation is specific to threat-association learning versus the more general attentional component that this system is known for.

We agree that our prior presentation was confusing, and we have removed the correlation plots (see new Figure 6 and new Figure 6-Supplement 1). We now use data from a distinct cohort of mice (distinct from the ADCD labeling experiments) that underwent behavioral testing to demonstrate high vs. low responder behavioral performance (new Figure 6A/B). To the reviewer’s point regarding the denominator, counts of ADCD+ neurons from animals that underwent the recall session were normalized to the average of number of ADCD+ neurons in home cage animals in their respective cohorts.

5) Some important controls are missing from the 'engram' causal manipulation study (Figure 3) which makes it difficult to definitively support the authors' contention that the captured cells are in fact 'engram' cells. The established definition of engram cells refers to those which store memory and are activated during learning and memory expression. While the authors do show conditioning-induced changes in CS activation in NBM-ACh cells, they don't have the necessary controls for their causal manipulations to rule out other non-engram factors which could account for their results. For example, it is possible that the captured cells are simply activated by auditory cues (even prior to learning, the immuno approach in Figure 2 may not be sensitive enough to detect this) and that they are transmitting CS-activated ACh signals to BLA which are necessary for memory expression. It is also possible that inhibiting a small number of NBM-ACh neurons (not restricted to 'engram' cells) reduces memory expression (early engram studies by the Josselyn lab controlled for this). The authors could address this with further experiments such as including 1) a control group in which cells are captured (off-dox) during a habituation session with CS presentation prior to on-dox learning and 2) using a dilute virus to manipulate a similar number of cells as in the engram group. Alternatively, they could present a more careful discussion of this issue and potentially avoid the use of the term 'engram'.

Point 5.1: We have measured NBM/SIp cholinergic response to the CS tone using two distinct approaches. First, we have used GRABACh3.0 to quantify changes in ACh levels in the BLA in response to the CS tone (New Figures 1, and Figure 1-Supplements 1-4). Second, we use an immediate early gene transcriptional response to identify cholinergic neurons that were activated in response to various conditions, (home cage, CS alone, USCS followed by CS). These latter experiments have used a viral approach (ADCD), immunostaining of Fos-shGFP transgene expression and immunostaining of endogenous Fos protein expression. Within the limits of sensitivity of each approach, we have not seen a significant increase in NBM/SIp activation in response to the tone (CS) unless the tone had been paired with foot shock. For example, we did not see increases in ACh release in the BLA in response to naïve tone (first tone presentation during training, Figure 1C and Figure 1-Supplement 1), to three repeat tones followed 24 hr later by tone alone (Figure 1Supplement 2) or in response to tone 24 hr after three foot shocks (Figure 1-Supplement 3). Likewise, we did not see activation or re-activation of NBM/SIp cholinergic neurons in response to tone (above levels seen in home cage animals), unless the tone had been paired with foot shock (Figures 2, 3 and Figure 6-Supplement 1).

Point 5.2: The BLA receives cholinergic input from both the NBM/SIp and the VP/SIa (roughly in a 2 or 3:1 ratio based on the number of cholinergic neurons retrogradely labeled, see for example Figure 5). We have compared silencing of each of these populations on freezing in response to the CS (using AAV-DIO-hM4Di injected specifically into either the NBM or the VP of Chat-IRES-Cre animals, new Figure 5-Supplement 2). Silencing the NBM cholinergic neurons eliminated a significant freezing response to the CS whereas inhibiting the VP cholinergic neurons did not affect freezing in response to the CS. The lack of effect on freezing in response to the CS following VP cholinergic silencing is not consistent with the possibility that generally reducing ACh in the BLA reduces memory expression.

6) The authors suggest that the ACh 'engram' cells participate in the expression of learned defensive responses through projections to the BLA. However, they have not demonstrated this definitively. To do this would require testing the effects of the inactivation of ACh-engram neurons which project to BLA on learned freezing responses. This may be technically challenging, but one approach could be to use the ADCD technique they've developed and inactivate terminals of captured cells in the BLA (chemogenetically or optogenetically). At a minimum, the authors should examine whether retrogradely labeled ACh cells projecting to BLA are Fos activated during a recall session (compared with appropriate controls).

We retrogradely labeled cholinergic neurons projecting to BLA using the retrograde tracer FastBlue in [chat-cre X fos-tTA/shGFP] mice. We simultaneously injected ADCD into the NBM/SIp. Animals were trained off-dox. Mice were sacrificed following recall and we quantified BLA-projecting (Fast Blue+) cholinergic neurons that were activated during training (ADCD+) and reactivated during recall (GFP+). ~20% of BLA-projecting cholinergic neurons were reactivated (engram) in the Training + Recall group, while only ~4% of BLAprojecting cholinergic neurons are reactivated in the home cage group (see new Figure 4). These data, combined with the results of silencing BLA-projecting cholinergic neurons, show that BLA-projecting cholinergic neurons are indeed part of the cholinergic threat memory engram.

7) For most experiments, significantly more raw data should be shown (e.g. raw example traces for GRAB-ACh3.0), and also brain section images for almost all experiments. Significantly more raw data should also be shown in the Main Figures. Please see Reviewer # 1, "Specific points relating to data presentation and statistical analysis", for details and specific examples."

We have revised all figures and added additional data to both main figures and supplemental figures to include more raw data to accompany the quantified data. In addition, we have provided additional detailed information on our statistical analyses (both in Figures, legends and in the new Statistical reporting table). Almost all plots now show individual data points. We have extended the GRABACh3.0 traces to show the full recording intervals (shown as average plots with SEMs and as individual traces). Specific changes include:

Figure1/Figure 1-Supplements1-3- we have changed the box plots to scatter plots.

Figure 1-Supplement1C/ Figure 1-Supplement 2D/ Figure 1-Supplement 3B- we have added the averaged fiber photometry trace over the full duration of the session.

Figure 1-Supplement 4C- we provide traces from individual mice showing responses to naïve tone vs. the conditioned tone at tone onset. And in panel B we show the averaged response for the full 28s (last 2s of naïve tone have the shock response and have been removed from the trace to preserve the scale of the y-axis to see tone responses).

Figure 1-Supplement 4A- we provide both raw images and atlas images showing mapping of fiber tip placements and GRABACh3.0 expression in the BLA.

Figure2B/3B- we show a schematic that delineates the region imaged and show both low and higher magnification images.

Figure 3-Supplement 1- we show a time-resolved freezing plot for ADCD-DREADD recall experiment and clearly delineate the time bins used for analysis presented in the main figures.

Figure 5B/C- we show a schematic delineating anatomical orientations.

Figure 5-Supplement 1A- we show images of GFP expression across the A-P extent of the BLA showing spread of the viral injection.

Figure 5-Supplement 1B- we show images of Fos IHC following recall from home cage, sham and cav-injected mice with clozapine delivered during training across the A-P extent of the BLA.

Figure 5-Supplement 1C- we show an atlas overlay with a low mag image of the amygdala with Fos IHC from control and cav-injected mice following recall with silencing during recall.

Figure 6-Supplement 2A- we provide low and high magnification images of ChAT+ neurons expressing Fos across the A-P axis of the NBM/SIp between high and low responders.

Figure 6-Supplement 2B- we show hM4Di-mCherry expression across the basal forebrain following injection of CAV-DIO-hM4Di-mCherry in the BLA of Chat-IRES-Cre mice.

Figure 8-Supplement 1A- we show low and high magnification images of ChAT and FosGFP labeled cells in the hDB, VP/SIa, and NBM/SIp following predator odor exposure.

Reviewer #1 (Recommendations for the authors):

Specific points relating to data presentation and statistical analysis:

1) The reporting of statistical data (N's, p- values, and the test used) is a bit cryptic. In general, the authors seem to report these values in the Figure legends. This, however, has the advantage that the reader has to constantly jump between the results text and legends text if he/she wishes to check the values for a given experiment. I would therefore suggest that the authors report all values of statistical tests in the results text (N's, p- values, and the test used).

Also, in some instances, N-values seem to be missing from the legend texts. Please check all legends / Results text accordingly.

We have now included a separate statistical reporting table which details sample size, statistical test, and p-value for each comparison reported in the manuscript. In addition, these details have been included in the figure legends for reference along with the corresponding figures. We have also included p-values of the statistical comparisons directly in the Results section.

2) Figure 1:

– Traces of GRAB-ACh3.0 should be shown more extensively. In Figure 1B, only the onset to the CS is shown – the response along the entire tone (30 s) should also be shown.

– Furthermore, the trace in Figure 1B seems to be an average over all mice. In addition, traces from individual mice should be shown, to allow an estimation of the S/N ratio of these measurements.

– In the quantification (Figure 1B, right), it would be easily possible to show the individual data points (there are 11 mice). I did not understand what the leftmost quantification tells us ("Pre-tone"). Please explain how this was calculated, and inasmuch it is relevant.

– What is the unit "Z%deltaF/F0"? In my view, this should either be Z-scored or else, deltaF/F0 – please clarify.

Please find our response to each individual point below:

1. We have now added the full GRABACh3.0 traces showing the average +/- SEM throughout the tone (Figure 1-Supplement 2D and Figure 1-Supplement 4B) and tone-shock pairings (Figure 1-Supplement 1C).

2. To aid in estimation of the S/N we have now included individual traces of data shown in Figure 1D in Figure 1-Supplement 4C. In addition we have included individual responses for the tone only controls (3X tone presentation on Day 1 without shock followed by tone presentation on Day 2) during naïve tone (Tone 1) presentation and 24h tone presentation in Figure 1-Supplement 2C.

3. We have now included individual data points (biological replicates) for all figures. We define the period prior to stimulus onset as the “pre-tone period.” This can better be referred to as baseline and has been amended accordingly in all figures and in the text. In Figure 1C you will note the baseline and tone-response period are indicated by a black bar above. We use 1s prior to tone onset for baseline measurements and the 1s following tone onset for our tone measurements. To combine our ‘slope’ and ‘max’ metrics and better represent the data, we have now reanalyzed our ACh release data to compute the area under the curve (AUC) at baseline or in response to the tone. Please note the AUC metric has now been incorporated throughout Figure 1 and Figure 1-Supplements 1-3.

4. Our GRABACh3.0 data is displayed as a z-score of the change in fluorescence over baseline fluorescence (deltaF/F0) over time. We have labeled the axis of each figure accordingly as “z-score X time”.

3) Figure 2:

– the authors should extensively show images of brain sections. It is curious to see that the entire paper doesn't show a single image at low magnification, which shows the investigated basal forebrain areas, and stained neurons in there, in their anatomical context. This would allow the reader to put the images into context, and to roughly validate whether the analyses by the authors (e.g. counted cells) make sense. For example, for Figure 2A, please show an additional low mag image of the corresponding brain area, that shows this brain area (NBM/SIp) in the wider context. The actual image can then be shown as an insert.

– In the plot of Figure 2B, a tick at the value of 0 is missing – please add this (this also applies to Figure 2C, 4B, 4C, 5A, right, 5B, right, 5C, right, 7A, 7D, S6B, S6C, S7A-C – please add all ticks at zero). Adding the tick at zero will help the reader eyeball the values of the data points.

– In some plots, the y-axis doesn't start at 0, but at an artificially chosen higher value. This increases "graphically" the effect size; however, it doesn't allow the reader to look at the data in the most unbiased way (indeed, some Journals do not allow this practice) (see Figure 7C, Figure S6A, right). Please newly plot the data from zero.

– The authors did not report in the Statistics part, which criteria they used to add "stars" onto the plots. For example, Figure 2B has *** (3 stars), but the reported p-value is p = 0.0017. According to the criteria used by many authors, p < 0.05 (1 star), p < 0.01 (two stars), and p < 0.001 (three stars). Thus, Figure 2B should have 2 stars, but not 3. Please explain the criteria you used, and please validate all star symbols according to the criteria you used.

– What is the unit of the cell numbers plotted in Figure 2B, C? Were these indeed absolute numbers as the axis labelling seems to suggest, or was there a normalization for the analyzed area? Please comment.

1. We have included additional supplemental figures to accompany Figures 1, 5, 6, and 8 that include raw images of IHC data for added transparency. We have also added schematics showing the exact location where analyses were conducted for Figures 2 and 3.

2. We have added tick marks at 0 for all plots.

3. We have updated all our plots to start at 0.

4. The criteria used for significance is as follows: * ≤0.05, **≤0.01, ***≤0.001, **** ≤ 0.0001. All plots have been checked to ensure consistency.

5. In Figure 2, we report the average number of double labeled (mCherry+ and GFP+) cells within each condition from slices at bregma -0.7/-0.8mm. This value is not normalized to area since the NBM/SIp lacks easily defined boundaries.

4) Figure 3:

– While the freezing data is shown at a quite high time resolution for the training day (Figure 3B), for recall (Figure 3C), only excessively binned data are shown. It would be advantageous for the reader to see the data at the same time resolution, in a quasi-continuous trace, to evaluate the effects also of silencing.

We have provided ‘quasi-continuous’ traces for both the training and recall sessions from Figure 3 in new Figure 3-Supplement 1A/B. Above the tracers we have indicated tone/shock onset, and the time-windows that were used for analysis of the recall data (pretone vs. recall tone response period). The binned comparisons of these data during the training and recall sessions are shown in Figure 3C/D.

5) Figure 4:

– There are several issues with the definition of high and low freezers, and with the analysis in this Figure.

– Hig+low freezers: The Methods state "High responders were defined as those mice that exhibited at least a 10 percentage point increase in % time spent freezing in the bin during the tone from the bin immediately pre-tone (e.g. Pre-tone freezing 10% to tone-induced freezing of {greater than or equal to}20%)." However, the Results text speaks about a "3-fold increase" (l. 235) – these are different criteria, please explain.

– Related to the above, the y-axis labeling in Figure 4B is unclear – is this the difference in freezing (as the Methods would suggest) or a fold-change (like the Results would suggest)? Also, the 3-star symbol for p = 0.005 seems not justified, should be 2 stars according to often-used criteria (see also above).

– l. 238 states: "There was no difference in mCherry expression in 'Low Responders' compared to the home cage group (fold change ~ 1),"

– this should be tested statistically.

– I don't see that the freezing data in Figure 4D (y-axis) shows a "positive correlation" with the number of Fos-labelled cells (x-axis). Indeed, looking at the high responders only, the data falls very close to a flat line. Please explain how the correlation was analyzed – this was probably done by fitting all data, i.e. also including the "low" responders (?). The obvious absence of a correlation for the "high responders" alone (black data points) is worrisome.

– In Figure 4D, the location of the data points (y-axis, freezing; about -10 – 0%) does not seem to correspond to what is plotted in Figure 4B (~ 2%; please explain).

– Figure 4F, and Figure S5: The result about the location of labelled cells along the a-p axis of the brain should be documented with a series of original images, both at low and high mag (see general comments, above).

We have amended this figure significantly to state our point more clearly. Please note: this is now new Figure 6.

1. We have added two panels of data into Figure 6 to better explain our responder criteria. In Figure 6A/B, we stratified the freezing response of animals into two groups: High vs. Low responders. High Responders were defined as mice who showed a >10 percentage points increase in time spent freezing in response to the tone compared to the pre-tone period. In Figure 6A, the full time-course of the recall trace is shown to highlight the difference in response to the tone between High and Low responders. In Figure 6B, the dotted line at 10% marks the criteria used to stratify the two groups. The methods section has now been edited to clearly state the criteria to match the figure (Page 43 Line 1231-1239). This is also described along with the associated figures in the Results section on (Page 11-12 Lines 329-343).

2. We have now presented this as the difference between pre-tone and tone period according to the criteria in Figure 6A. The statistical testing is now done and shown on the quasi-continuous trace of recall. The details of post-hoc tests comparing each time bin to the other can be found in the statistical reporting table.

3. We have added a grey shaded area to Figure 6D to indicate the range of ADCD+ neurons captured in home cage animals (ADCD+ neuron counts in each home cage animal normalized to average of all home cage animals). We have included a formal statistical comparison between home cage and low (home cage vs. low, fold change ~ 1, p>0.99) or high responders (home cage vs. high, fold change ~3, p=0.0121) in the associated figure legend and in the Results section on Page 12 Lines 350-353.

4. We have removed the correlation plots and analyses from this figure. Instead, we believe our data are consistent with a threshold effect, whereby mice with ~2 fold increase in activation of cholinergic neurons results in the freezing behavior distinguishable from baseline freezing. Additionally, we find that this “critical mass” of activated cholinergic neurons is localized to a specific region within the anterior NBM/SIp of the basal forebrain.

5. We have added raw images of ChAT and Fos IHC following recall in Figure 5Supplement 2A from a High and Low responder mouse to supplement the schematic shown in Figure 6E. Quantification of counts of ADCD+ cells captured during recall across the A-P extent can be found in Figure 5-Supplement 2C.

6) Figure 5:

– The location of back-labelled cells should again be amply documented with several low and high mag images from all analyzed areas (hDB, VP/SIa, NBM/SIp).

– 3 stars in Figure 5B might be 2 stars.

1. We have provided injection site mapping for data presented in Figure 5A in Figure 5-Supplement 1A. Images of CAV back-labeled cells are provided in Figure 6Supplement 2B across the hDB, SIa, and NBM/SIp.

2. We have replotted Figure 5b’’ (previously 5B) data as biological replicates rather than technical replicate. The significance value has updated accordingly in the figure as well in the text (p=0.0286, *).

7) Figure S4:

– was the data in Figure S4C a further analysis of data in Figure 2B? If yes, why are the animal numbers seemingly different between the two datasets?

We have now replotted and analyzed all our data by biological replicates and have rectified this discrepancy. Figure S4C can now be found as new Figure 6-Supplement 1A.

8) Figure S5:

– The legend doesn't state the N-numbers.

– The box plot should be explained (meaning of the boxes). Even better, individual data points should be shown instead, if this is feasible.

1. We have now included sample sizes in the figure legends and in the statistical reporting table.

2. This is now Figure 6-Supplement 2C. We have converted this plot to individual data points as suggested. The data is described in the Results section on Page 12 Lines 363-366.

9) Figure S6:

– For the data in A, again example images at low and high magnification should be shown …

– For the images in B, please indicate lateral-medial, and dorsal-ventral directions in these images (assuming that these are coronal sections?).

– The data in C again suggest that there might also be an effect on contextual fear learning (see also above).

1. Newly added Figure 5-Supplement 1A/B shows a detailed A/P mapping of the injection sites and Fos labeling for experiments presented in Figure 5.

2. Figure 5-Supplement 1 images have been additionally labeled with medial/lateral (M/L) and dorsal/ventral (D/V) designations for ease.

3. We find no statistically significant differences in pre-tone freezing between groups (Figure 5B, Sham vs. hM4Di, Pre-tone p=0.0679; Figure 5C, Sham vs. hM4Di, Pretone p=0.0966).

10) Figure S7

– The text describing the data in A states "did not freeze to the tone" (l. 305). This is not a scientifically warranted statement, since there seems to be residual freezing. Rather, the data should be tested statistically.

– The N-numbers for the data in A, and B should be stated in the legends.

– For these and possibly similar experiments, the authors need to show in detailed histological images, for each individual mouse, that the expression of hM4Di has occurred in the expected set of cholinergic BF neurons.

1. We have amended our text in the Results section to note that silencing BLA projecting cholinergic neurons during training (Figure 5B) blunted the recall induced freezing behavior (p<0.0001) and subsequent activation of BLA neurons (p<0.0286). Similarly, silencing BLA-projecting cholinergic neurons during recall (Figure 5C) reduced recall-induced freezing behavior (p=0.0279) and the activation of BLA neurons (p=0.0317). Amended text in the Results section can be found on Page 9-10 Lines 271-282. Further details of the statistical comparisons can be found in the Statistical Reporting Table.

2. Sample size has been included in the figure legends and in the statistical reporting table.

3. We have included additional raw images to support the data. An example injection site is now included as part of Figure 5A, center panel, along with an example of ChAT+ cells that also have mCherry expression. The distribution of these hM4diexpressing cholinergic neurons is quantified in Figure 5A, left panel. Additionally, images across the A-P axis of the basal forebrain of hM4Di-mCherry expressing cells have been provided in Figure 6-Supplement 2B. We visually inspected each animal in these experiments for bilateral GFP expression in the BLA and presence of mCherry+ cells in the basal forebrain prior to including (or excluding) mice in the analysis. We have provided representative images showing the A/P extent of the GFP expression (injection area) from one animal in Figure 5-Supplement 1A.

11) Figure S8:

– I wonder why the firing rate data of the "Training" group in panel H is extremely small (around 5 Hz or less), and smaller than in all other groups. Is there another type of modulation of intrinsic excitability, a down-modulation, after training? Please explain.

We have also been intrigued by this observation. However, we reserve any conclusions regarding this decrease as it is not statistically significant and overlaps with the firing rate from neurons recorded from the home-cage condition. Further experiments will be necessary to investigate whether there is a bonafide decrease in excitability in Fos+ BFCNs following training.

12) Figure S9:

– In panel C, could you show the freezing data at a higher time resolution, e.g. the 30s bins also used in Figure 3B?

Please find the requested plot in Author response image 1. The associated summary plot for the recall session is now included in new Figure 5-Supplement 2A.

Author response image 1.

Author response image 1.

Reviewer #2 (Recommendations for the authors):

– For the flow of the story, I recommend switching the order of presentation where the lack of effect on behavior after silencing the CeA-projecting NBM/SIp neurons (lines 286-301) comes after the data showing the effect of silencing the BLA-projecting NBM/SIp.

We thank the reviewer for this recommendation and have rearranged the results and figures accordingly to address this. The amended results can be found on Page 10-11 Lines 297-314.

– I would also recommend moving figures S7A and S7B into the main figures as they provide important evidence for the author's point.

We agree that this set of experiments is critical to the main point of the manuscript. We have taken this suggestion and have now included these experiments into a significantly modified main figure, Figure 5 (Panels B-C).

– Given that the authors don't see an effect of clozapine on memory when inactivating during the consolidation period (Figure S7C), it might be worthwhile to discuss why this might be the case, and the effects of cholinergic input starting cascades of neural activation in the BLA during training.

We agree the possible function of ACh during consolidation is of interest, but feel that this requires a more complete experimental investigation that is beyond the scope of the present study, and as such we have removed prior Figure S7c as too preliminary.

– Overall, there are missing Bregmas and mapping of placements in several places (e.g. Figure 1A – Bregma and mapping of placements/spread of the virus across animals, it appears that the recordings are relatively anterior in the BLA; Bregma in Figures 5, S3C, 6SB, S6C – not clear how much of a CeA there is at this location, 7D; a zoomed out section of where the example data in figure 7C is coming from).

We have now included raw images and schematics to provide transparency to the data. These have been added to Figures 2, 3, 4, and Figure 1-Supplement 4, Figure 5-Supplement 1, Figure 6-Supplement 2, and Figure 8-Supplement 1. In new Figure 1-Supplement 4A we have included a schematic of the relocalized fiber placements from fiber photometry experiments presented in Figure 1. In Figures 2 and 3, we have now included low magnification and high magnification images of ADCD labeled neurons within the NBM/SIp. An A/P extent of the Fast Blue injection locations for the newly added Figure 4 is included as panel within Figure 4A. We have now included raw images detailing the AP extent of the evaluated BLA as related to Figure 5 data, shown in Figure 5-Supplement 1A/B. In Figure 5-Supplement 1C, we have included low and high magnification images to localize the CeC within the image. In Figure 6-Supplement 2 we provide images across the A/P extent of the basal forebrain to show where activated cholinergic neurons are clustered in High vs. Low responders. Finally, in Figure 8-Supplement 1, we show raw images of across the basal forebrain of activated cholinergic neurons in response to predator odor. We believe these additions have substantially strengthened the manuscript.

– It would also be helpful to indicate the area where the signal fluorescence is being calculated in traces shown in figures 1 and S1.

We now indicate with a black bar above our traces in Figure 1 the period that is being quantified for baseline (1s before stimulus onset) and tone response (1s following tone onset). We have also included additional text in the methods section to clarify this on Page 45 Lines 1288-1291.

Reviewer #3 (Recommendations for the authors):

1) For the 'engram' data in Figure 2, the logic of the double labeling is a bit confusing as described and it's not clear, from a simple reading of the paper, how the authors achieve double labeling of encoding (training activated) and recall activated cells. After digging around in the Methods, Figs/legends I was able to figure it out, but it was difficult initially. The authors should state this in a simple way in the main text along with a description of how this works at the level of transgene expression at different time points.

We have updated the text to clarify our double-labeling strategy (Page 6 and 7 Lines 170188). We have also added a schematic in Figure 2-Supplement 1B accompanied by detailed text in the figure legend (Page 28 Lines 807-817) detailing how the ADCD-engram labeling works. In brief, we deliver the ADCD construct in the basal forebrain of Chat-Cre X FostTA-shGFP mice. In this manner, the virally delivered ADCD construct can only be expressed in cholinergic neurons, in the absence of doxycycline (DOX OFF). Once expression initiates, it is maintained even in the presence of Dox by a dox-insensitive tTa, thereby resulting in permanent labeling. Intercepting the peak of a shGFP signal (for example 2.5h following a second behavioral assay) allows for identification of cells that express mCherry (activated during DOX OFF) and cells that express GFP.

2) Some important controls are missing from the 'engram' causal manipulation study (Figure 3) which makes it difficult to definitively support the authors' contention that the captured cells are in fact 'engram' cells. The established definition of engram cells refers to those which store memory and are activated during learning and memory expression. While the authors do show conditioning-induced changes in CS activation in NBM-ACh cells, they don't have the necessary controls for their causal manipulations to rule out other non-engram factors which could account for their results. For example, it is possible that the captured cells are simply activated by auditory cues (even prior to learning, the immuno approach in Figure 2 may not be sensitive enough to detect this) and that they are transmitting CS-activated ACh signals to BLA which are necessary for memory expression. It is also possible that inhibiting a small number of NBM-ACh neurons (not restricted to 'engram' cells) reduces memory expression (early engram studies by the Josselyn lab controlled for this). The authors could address this with further experiments such as including 1) a control group in which cells are captured (off-dox) during a habituation session with CS presentation prior to on-dox learning and 2) using a dilute virus to manipulate a similar number of cells as in the engram group. Alternatively, they could present a more careful discussion of this issue and potentially avoid the use of the term 'engram'.

Throughout the manuscript, we have evaluated the cholinergic response to the tone using 1) GRABACh 3.0 to quantify changes in ACh levels in the BLA and 2) using an IEG response to tag and identify cholinergic neurons that were activated under specific behavioral conditions. Using either approach, we have not seen a significant increase in activation or ACh release in response to the unconditioned tone (Figure 1C-left, 1D-left, 2C/D, Figure 1-Supplement 2, Figure 1-Supplement 3, Figure 6-Supplement 1A) unless it has been paired with a footshock (see Figure 1C-right, 2C, 3D).

The BLA receives a majority of cholinergic input from the NBM/SIp with additional cholinergic input from the VP/SIa. We evaluated the contribution of each of these populations to the freezing response to the conditioned tone. We find that silencing the NBM/SIp results in a decrease in freezing to the conditioned tone (Figure 5-Supplement 2A), while silencing the VP/SIa population does not affect the freezing response (Figure 5Supplement 2B). The lack of effect with VP inhibition suggests that silencing a random, small number of BFCNs, or decreasing ACh release in the BLA alone is not sufficient to affect memory recall.

3) The authors suggest that the ACh 'engram' cells participate in the expression of learned defensive responses through projections to the BLA. However, they have not shown this definitively. To do this would require testing the effects of the inactivation of ACh-engram neurons which project to BLA on learned freezing responses. This may be technically challenging, but one approach could be to use the ADCD technique they've developed and inactivate terminals of captured cells in the BLA (chemogenetically or optogenetically). At a minimum, the authors should examine whether retrogradely labeled ACh cells projecting to BLA are Fos activated during a recall session (compared with appropriate controls).

In new Figure 4, we retrogradely label cholinergic neurons that project BLA using the retrograde tracer FastBlue in [Chat-IRES-Cre X Fos-tTA/shGFP] mice. Simultaneously, we injected our ADCD virus into the NBM/SI. Cells were captured off-dox during training. Mice were sacrificed following recall to capture GFP expression in recall-activated cells. We quantified BLA-projecting cholinergic neurons (Fast Blue+) that were activated during training (ADCD+) and became reactivated during recall (GFP+). We find that ~20% of BLAprojecting cholinergic neurons are reactivated (engram) in the Training + Recall group, while only ~4% of BLA-projecting cholinergic neurons are reactivated in the home cage group. These data show that BLA-projecting cholinergic neurons are indeed part of the cholinergic threat memory engram.

4) In Figure 5, the authors report a reduction in BLA Fos-activated neurons with chemogenetic inhibition of BLA-projecting ACh neurons. However, this experiment requires a control group without training (or, ideally, with unpaired CS/US training) to demonstrate that the inhibition is affect learning-induced activation of BLA neurons. While others have reported this using other paradigms previously, it is still important to show here using their own preparation.

We have now included images across the A/P axis in the BLA for Fos expression in home cage animals alongside Fos expression following recall in animals that had a sham surgery or those that were injected with CAV2-DIO-hM4Di with cholinergic inhibition during training (Figure 5-Supplement 1B).

[Editors’ note: what follows is the authors’ response to the second round of review.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

We thank the authors for their response to the reviewers' previous suggestions. Many of the findings have been strengthened through the addition of new experiments and analyses. However, there are remaining concerns about the statistical approaches used and low population size in some experiments which should be addressed. We suggest using the 'eLife Transparent Reporting Form' as a guide. Related to sample sizes, power analyses or clear demonstrations of existing standards in the field should be used to justify animal numbers. If necessary, further animals should be added to experimental groups (e.g. the differential role of NBM/SIp and VP/SIa in learned vs. innate threats) or, as suggested by reviewer 1, data with low n could be removed if it does not seriously undermine the conceptual advance of the paper. In addition, the strong focus on the activated basal forebrain cholinergic cells as being part of the threat memory engram is not well supported. The authors should moderate the claims related to this aspect of the paper.

Furthermore, the title and/or abstract should provide a clear indication of the biological system under investigation (i.e., species name or broader taxonomic group, if appropriate). Please revise your title and/or abstract with this advice in mind.

We have addressed these concerns in the following manner:

1. We excluded the data related to the lack of engagement of the VP in cue associated-learned threat from the manuscript, as suggested, as this does not compromise the main claims of the revised manuscript.

2. We removed the use of term “engram” from the title, abstract and results text. We do discuss the definition of “engram” in the Discussion section.

3. We tempered claims regarding our demonstration of a cholinergic engram throughout the text and now consider arguments pro and con in the context of the literature more carefully in the Discussion section.

Re: clarifying species used in title and/or abstract:

– Line 42 (Abstract) now reads: “Using a genetically-encoded acetylcholine (ACh) sensor in mice…”

– Lines 54-55 (Abstract) is also modified to read: “…specific types of memory within the cholinergic basal forebrain of mice.”

Reviewer #1 (Recommendations for the authors):

Specific points (in order of appearance)

1) l. 156/157 "when it was not explicitly paired with a shock" (data in Figure 1S3D).

– Please explain in the legend the exact temporal sequence between shock(s) and tone(s) used for this "non-explicitly paired" protocol (or does this refer to the "shock-only" protocol?).

2) Figure 1 S3C, D: The legend to this Figure states that here, the results from a "shock-only" training are displayed for the subsequent "recall" session.

First, the authors should also show the freezing response for the standard group undergoing standard fear learning (3x tone-shock pairings); this is important for a comparison to the "shock-only" training group.

Second, the Results text states that Figure 1S3D shows the ACh response when "it was not explicitly paired with a shock".

– Is this the same as the "shock-only" protocol described in the legend? Please clarify, and make changes in the manuscript as appropriate.

We provided schematics as well as text to explicitly delineate the behavioral protocol and to describe each result obtained. Please note that the full context of the phrase that was excerpted by the reviewer "it was not explicitly paired with a shock” reads as follows:

“To verify that the increases in ACh release were indeed specific to the tone-shock association and not due to generalization from prior shock exposure, we also subjected mice to 3 shocks (day 1) followed by a tone presentation 24 hr later (day 2) (Figure 1- Supplement 3A). While mice demonstrated freezing behavior during the session on day 2, there was no significant increase in freezing behavior to the 24h tone presentation (Figure 1-Supplement 3C, p=0.2418). There was no increase in ACh in response to the tone when it was not explicitly paired with a shock, confirming that the changes in ACh release were indeed associative (Figure 1-Supplement 3D; baseline (pre-tone, day 2) to 24h tone (tone presentation, day2): p=0.7272).”

When read within the context of the entire paragraph it seems clear that we are referring to the shock alone training, followed by a tone alone recall. The grouping of these sentences into a single paragraph, focused on discussing a single figure (that starts with a schematic of the experimental design) is as explicit a presentation of the protocol and the data obtained as we could put together. However, we welcome further clarification from the reviewer.

3) Figure 1 S4 was not called out from the Results text; this should be changed i.e. what was done in this Figure, should be (briefly) mentioned in the Results.

We had included these additional plots of the data and schematics in response to reviewer request during the last review. As stated, these plots refer to the data that has been described in the text (Figure 1). We believe that the figure legend accompanying Figure 1 S4 clarifies the contents of this figure.

See lines 125-126 where we reference Figure 1 S4:

“Given this, the first question we asked was whether acetylcholine was released in the BLA during associative threat learning (Figure 1 and Figure 1-Supplements, S1-S4).”

4) The use of star symbols in Figure 2C, in conjunction with what is stated in the corresponding Results text, is misleading and tends to over-interpret the statistical significance of this data set.

We use a single star to denote p≤0.05 (as defined in the manuscript; line 1284, as quoted below and as standard usage in some journals and statistical analysis algorithims). The identified major issue is the use of a single star in Figure 2C to identify a p value = 0.0464, (i.e. ≤0.05.) The use of the star and the accompanying text is just a statement of the p value and is not intended to be misleading or an interpretation.

Also please see Line 1293:

“p-value criteria: * p≤0.05, ** p≤0.01, *** p≤0.001, **** p≤0.0001

Specifically, l. 202 states: "Significantly more double-positive cholinergic neurons were seen following the complete (tone + shock followed by tone recall) associative threat learning paradigm compared to mice that underwent session 1 without shocks (Figure 2C, p=0.0249)."

– This is misleading. The Figure legend and statistics table state that the value "p = 0 0.0249" refers to the "overall" significance of the Kruskal-Wallis (KW) test used here, but NOT to the specific comparison cited in the above sentence. For that specific comparison, the p-value was higher ("p=0.0464, Dunn's corrected", see Figure legend).

– In response, the authors should spell out the "overall" significance of this dataset (effect of group; p = 0.0249; KW) in the Results section, followed by the pertinent p-value for important specific comparisons (in this case, p = 0.046 for the comparison between the "tone-tone" and tone+shock – tone" groups). The n-values for these pertinent comparisons should also be given in the Results text.

Furthermore, in this dataset (Figure 2C), the post-hoc test for the comparisons between the home cage (HC) and standard fear conditioning group is NOT significant (p = 0.72). This value should be indicated by a "bracket" symbol in Figure 2C, and the meaning of this non-significance should be discussed in the manuscript.

5) The p-values, n-values, and the statistical tests applied, should be reported together either in the figure legends or in the Results text (please choose one section to use for reporting these statistics, in addition to the statistical table) throughout the entire manuscript (see above point 4, and see also my previous point "Data Recommendation 1").

As suggested, we had included this information in the figure legends for which quantitative analyses was assessed.

As the reviewer directs, these statistical data should be reported “either in the figure legends or in the Results text”. We had complied with this request, presenting the data in the figure legends and provided a detailed table for reporting statistics for full transparency.

6) Figure 2D:

l. 210 states:

"All three conditions showed few reactivated neurons and no differences between groups

(p=0.9471)"

– The finding that all these 3 control groups are not significantly different from each other, is not surprising and does not go to the point of the control experiment. Rather, it should be tested statistically whether the number of re-activated neurons in these three conditions is significantly different from the one after standard fear conditioning (the data in Figure 2C, right column).

We have performed the requested analysis and the updated statistical metrics have been provided in the statistical reporting table.

7) I think the title of the Sub-header:

"Reactivation of cholinergic neurons activated by training is required for learned behavioral responses" (l. 214) is wrong, it should read "Activation" instead (and then replace "activated by" with "during").

To my understanding, the chemogenetic silencing experiments in Figure 3 test the role of activated cholinergic neurons but do not suppress the "re-activation" of neurons (during recall…), since ADCD expression of hM4Di is driven during the training session.

The design of this experiment, as illustrated in the schematic, was to explicitly assay the effect of selective inhibition on neurons that were activated in training and then subsequently re-activated by recall. We have used a classic experimental paradigm to do so (Josselyn and Tonegawa, 2020. Science). That is, we expressed the inhibitory hM4Di construct in cholinergic neurons that had been activated during training (Activity dependent Cre dependent virus, dox off during training), We then use the DREADD ligand, clozapine, right before and during recall to selectively prevent these previously activated neurons (activated during training) from being activated again (hence, re-activated) during the recall session. The schematic in this figure and the figure legend attempt to clarify this further as follows:

“A. ADCD-hM4Di was injected into the NBM/SIp of Chat-IRES-Cre x Fos-tTA/shGFP mice. Two weeks later animals underwent training on regular chow (Dox chow removed 24 hr prior to training session) to allow hM4Di.mCherry to be selectively expressed in training activated cholinergic neurons. Three days later, recall was tested in Dox on conditions. Clozapine (CLZ) was injected 10min before the recall session to activate the inhibitory DREADD, hM4Di specifically in previously activated cholinergic neurons.”

8) In Figure 3 S1B, the average trace of time-resolved freezing is shown, for the same data set as the one in Figure 3D. However, the number of mice included in the time-resolved display was lower (n = 5 and 6) than in the display of Figure 3D (n = 7 and 8). This must be homogenized, i.e. all data used for the statistical analysis in Figure 3D (n = 8 and 7 mice) should also be shown in a time-resolved fashion for Figure 3S1B.

The reviewer is correct in noting that two mice are missing from each group in the population data, time resolved plot. The reason for this omission is that those mice were part of a cohort, where the video aspects of the data were temporarily stored on a drive that was damaged in a major flood that affected multiple labs within building 35 (NIH; Bethesda) last December. As such, we could not go back to generate the specific time-resolved plots for those 4 mice. We provided the data for the remainder of the mice to aid in clarity, as requested. We also reanalyzed our data by removing the missing mice from the summary data plots and this did not alter the statistical validity of our findings as shown in Author response image 2.

Author response image 2.

Author response image 2.

9) l. 273 states:

" We injected cav.hM4DiBLA or sham control mice with clozapine (CLZ) 10 min prior to initiating cue-conditioned threat learning"

– Please add here that these were IP (?) injections, for a quick reference for the reader. Indeed, this information should possibly also be added to the Methods (I did not find it upon a rapid check).

A search of the text reveals at least 4 instances in which we have detailed that the mice were injected i.p. Lines: 1126, 1132, 1169, 1183. In addition, we have now added another note regarding this procedural detail in the Results section at:

Line 215 “(0.1 mg/kg; injected intraperitoneally (i.p.)).”

10) The authors should control for the possibility of leaky expression of hM4Di, following presumed Cre-dependent retrograde expression from a CAV2 virus after injections in the BLA (Figure 5). Because, as I assume, CLZ was not injected into the basal forebrain (BF) but rather IP (see point 9 above), a leaky expression of hM4Di beyond the cholinergic BF has the potential to strongly change the interpretation of this experiment. Thus, the authors have to show, first, whether hM4Di-mCherry is expressed at the injection site in the BLA itself (this would lead to a direct inhibition of BLA neurons). This should be done by showing example images of hM4di-mCherry fluorescence on the level of the BLA, and by analyzing the number of hM4Di-mCherry+ neurons in the BLA over the entire dataset. Second, the authors should show extensive raw example images of the BF at different a-p locations (and of further brain areas in a zoomed-out fashion), in an attempt to exclude the possibility that hM4Di might also be expressed in non-cholinergic neurons (ChAT could be visualized by IHC). Only a non-leaky expression of hM4Di-mCherry, essentially limited to ChAT-positive neurons in the BF, would support the conclusions of the authors.

Figure 5A shows that the expression of hM4Di is wholly restricted to ChAT+ neurons. We have no reason to suspect that the Chat-IRES-Cre mice used in this study display leaky cre expression. We refer the reviewer to the Jax Labs page for this mouse line, with particular reference to the details posted there relevant to the ways in which the line we have used was created to address a previous cre leak issue - https://www.jax.org/strain/031661.

Figure 6 Supplement 2 shows several zoomed out images of the BFCNs in response to this same request made previously by Reviewers.

11) On l. 298, the results of a control experiment is reported, which is not documented ("data not shown"). Some journals do not allow the statement "data not shown". The authors should document the control data set, or else remove this statement from the manuscript. See also l. 441, another "data not shown", which should be similarly documented or else removed.

We have now removed this phrase from the manuscript. We apologize for the prior inclusion of “data not shown” as we weren’t aware that eLife prohibits the inclusion of such references.

12) Data in Figure 5C.

When checking the p-values from the legends/statistics table, I was puzzled by the fact the comparison of the two "grey" data sets has a reported p-value of < 0.0001, although many of the data points between the two groups overlap. On the other hand, the comparison between the "sham" group and hM4Di group at the "pre-tone" time had a p-value of 0.096, although the overlap of the data points is only partial. Please comment, and validate the underlying tests.

These data were analyzed, as is appropriate, using a repeated measures two-way ANOVA (as noted in the figure legends and in the statistical reporting table). Given the design of this experiment and the exigencies of the statistical tests used, we think that this analysis is more quantitatively rigorous than the readers’ visual inspection of overlap between populations in a scatter plot. Please see the plot Author response image 3.

Author response image 3.

Author response image 3.

13) Related to point 12 above, there is a strong tendency that freezing before the tone is lower in the hM4Di group as compared to the "sham" group. Although the authors reported a p-value of 0.096 (and therefore, "not significant") for this comparison (Figure 5C), this p-value should nevertheless be verified. Because the trend is seen in this data set, as well as in several other ones (Figure 5B and Figure 3D), the authors should discuss the possibility that chemogenetic suppression of ACh release might, in addition to the effects studied here, also reduce contextual fear learning (see also my previous "Specific comment 2").

We have used standard, and in our opinions correct, statistical tests to evaluate the quantitative data in this study and therefore feel that it would be unscientific to discuss mechanisms underlying results that do not differ significantly as detected by these tests.

14) Still related to the data in Figure 5C:

The time-resolved display of tone-evoked freezing during the recall day should be shown both for this data set (Figure 5C), and for the data set in Figure 5B, similarly as was done above for the data in Figure 3D (see panel Figure 3 S1B). This would also be helpful to gain insight into the p-values resulting from ANOVA / KW testing (see above, point 12), which were probably performed on similar time-resolved data sets (because the effect of "time" was always tested).

We can provide these plots additionally as a supplemental figure. However, our summary plots in the main figures are the commonplace reporting standard for these types of data. We feel that adding still more supplemental figures and analyses that do not change the conclusions of the study would unnecessarily extend this manuscript length. The manuscript has already seen a significant “ballooning” of supplemental figures to show each individual data set obtained, in response to prior requests of the reviewer.

In case it is the consensus of the editors and other reviewers that presentation of each individual time course points is required, we provide them in Author response image 4 for inclusion at your discretion.

Author response image 4.

Author response image 4.

15) l. 299 states:

"Differences in recall-induced Fos expression between sham and cav.hM4DiBLA mice were maximal in rostral portions of the BLA (between bregma -0.8mm to -1.4mm) (Figure 5-Supplement 1B)"

– First, the analysis in Figure 5 S1B seems to have been done in only N = 1 mouse, which is problematic (please confirm the n-number, and state it in the Results text and in the legend to that Figure).

Second, I actually fail to see an a-p gradient as claimed by the authors (higher in anterior than posterior). If one looks at the number of Fos+ cells in the hM4Di group (red data points), then it seems that hM4Di suppresses the number of Fos+ cells, in relative terms, stronger at more posterior levels (e.g. reduction by > 80% at -1.2 mm) than at anterior levels (e.g. reduction by < 40% at -0.8 mm). Wouldn't this suggest that opposite from the author's interpretation, MORE cholinergic neurons are specifically (in a hM4Di-dependent manner) activated posteriorly than anteriorly? Please comment.

We have removed the plot shown in Figure 5 S1B as it is not critical to the point of the analysis. To address the reviewer’s point, anterior BLA is located from bregma -0.8 to ~-2.0 after which posterior BLA begins to predominate (these are referred to in the Paxinos mouse brain atlas as BLA and BLP respectively). In our studies, DREADD mediated suppression of Fos activity was seen in the regions analyzed and presented but not beyond those bounds. This observation was consistent with a previously published report showing the segregation of negative valence processing BLA neurons, and as such we had included that plot in this manuscript.

We would like to note that we make no comments or claims regarding the relationship between BLA anterior to posterior anatomy and the basal forebrain anterior to posterior anatomy. We are unsure what led the reviewer to this conclusion but would be willing to address it explicitly to aid in clarity if indeed this is deemed a major issue as noted by the reviewer.

16) Visual inspection of the data in Figure 5 S1D strongly suggests that the data in the control ("sham") group (n = 9, grey data) are actually the same as the data for the sham group in Figure 5C. In Figure 5C, retrograde expression from the BLA was studied, whereas in Figure 5 S1D, the effects of retrograde expression of hM4Di from the CeA were studied. If the experiment was originally set up for 3 groups of mice (thus, sham control – CAV2-hM4Di in BLA – CAV2-hM4Di in CeA), then this should be reported as such, the data should be shown in one panel, and the statistical comparisons should be done accordingly (i.e. the "multiple comparison" correction would need to be done for 3 groups). If, on the other hand, the CeA experiments were done without proper controls, and were simply compared to the control group from another experiment, then the CeA data should be removed. Please comment, and make appropriate changes in the manuscript.

Yes, in fact we in fact explicitly stated this in the figure legends:

“Note, sham group is the same as the sham group shown in Figure 5C as this experiment was conducted in the same cohort.”

The reviewer is correct in noting that this was a post-hoc identification of CeA injection sites in subset of mice. We are removing this data from the manuscript as the point made is not critical.

17) Figure 5 S2 claims that neurons in the "NBM/SIp" versus neurons in the "VP/SIa" have been targeted selectively with AAV9:DIO:hM4Di-mCherry. This should be documented extensively, by showing expression of hM4Di-mCherry throughout several brain slices in a de-zoomed fashion for an example mouse, and by showing average results from the quantification of the localizations in a suitable manner (for example, by plotting cell density/number along the a-p axis), for all investigated mice. Otherwise, it remains unclear whether the two areas have indeed been targeted differentially.

In our hands, the two injection sites used to target VP vs. NBM result in expression restricted to each region.

Nevertheless, based on comments and concerns raised in “point 19”, we have decided to remove the VP/SIa DREADD data from this manuscript.

18) In Figure 5 S2A, the "triple star" symbol comparing the two control data points might be removed. It is not surprising that the tone recall causes increased freezing in the "sham" group.

We respectfully disagree with the reviewer. While it might not be surprising that tone recall causes increased freezing in the “sham” group, it is important to validate this group as an appropriate control in this experiment. In addition, please note that it was the explicit recommendation of Reviewer 1 in the previous round of reviews that all comparisons are tested for statistical significance.

21) The "fold-change" analysis of expression analyzed in Figure 6D remained unclear to me (sorry if I didn't get it). Please explain this briefly in the Results text, and make sure it is explained in more detail in the Methods. Also, please give details about the viral construct used in Figure 6C (currently it only states "ADCD").

We clarified this in the “response to reviewers letter” submitted alongside the previously revised manuscript. Here is our prior reply:

Reviewer1 Data Recommendation 5:

We agree that our prior presentation was confusing, and we have removed the correlation plots (see new Figure 6 and new Figure 6-Supplement 1). We now use data from a distinct cohort of mice (distinct from the ADCD labeling experiments) that underwent behavioral testing to demonstrate high vs. low responder behavioral performance (new Figure 6A/B). To the reviewer’s point regarding the denominator, counts of ADCD+ neurons from animals that underwent the recall session were normalized to the average of number of ADCD+ neurons in home cage animals in their respective cohorts.

And please note that the figure legend states:

“D. Quantification of change in number of cholinergic neurons activated (ADCD+) in low or High Responders relative to the home cage. The number of ADCD+ neurons differed significantly between Low and High Responders (Mann-Whitney test, p=0.01) (n=3 Low Responder, n=8 High Responder). Grey shading represents the range of fold-change in ADCD+ cells in individual home cage animals relative to the average of all home cage animals (n=5). (Mann-Whitney test, home cage v. Low Responder, p>0.9999; home cage v. High Responder, p=0.0121).”

Also see Lines 326-327: Next, we quantified the fold change in the number of mCherry+ neurons in each group relative to corresponding home cage control mice (Figure 6D).

We believe that this is a clear description, but welcome specific suggestions as to what we need to do to make our procedures and analysis optimally transparent to the reviewer.

The viral construct used has been extensively documented in the text upon its introduction, in a supplemental figure dedicated to it, and in the methods section.

Additionally, Lines 324-325: we labeled cholinergic neurons activated during the recall session with ADCD-mCherry (on dox during training, off dox during recall; Figure 6C).

22) The picture shown in Figure 6E is simply a scheme of the idea that neurons activated in "high responders", in your view, should be localized more anteriorly, but it is in no way a documentation of an anatomical observation. This scheme should be replaced by experimental evidence.

This is indeed, as we mention in the figure legends, just a schematic. We would like to clarify, however, that this is not a hypothesized or speculated distribution but as previously requested, a documentation of our data. We include images showing this A-P difference in activation and quantification of a subset of high responding mice which were provided in Figure 6 Supp2.

23) The data that might relate to the above issue (localization of ChAT+, activated neurons along the a-p axis in "high" and "low" responders), might be contained in Figure 6 S2. However, the data seems to be only partially analyzed, and it is confusing to fully understand how the plot in panel C was derived. First, it was unclear from which group of mice the images in Figure 6 S2A were taken. Please explain in the legend from which experiment these images were taken, which mouse genotype was used (Fos-tTA/shGFP with or without ChAT-Cre ?), and which viral constructs were injected in these mice.

Once again, we urge the reviewer to read the figure legends and text associated with this figure. We provide ample description of how these data were gathered. Specifically:

“A. Representative images of the NBM/SIp at bregma locations -0.8, 1.0, and 1.3mm from a Fos-tTA/shGFP High Responder (top) and Low Responder (bottom) mouse sacrificed 2.5h after recall. Brain sections were stained for ChAT (magenta), and GFP to amplify Fos signal (green). ChAT+ and Fos+ co-labeled cells are marked by yellow arrowheads. Scale bar = 50µm. Insets show magnified images of ChAT+ and Fos+ neurons. Scale bar = 10µm.”

24) The images shown in Figure 6 S2C, should be shown in a de-zoomed fashion for several representative brain images on the level of the BF, and beyond. Such a documentation of the data would represent the control for "leaky" expression mentioned above (best after co-staining against ChAT; see point 10 above).

We refer the reviewer to point #10 above.

25) The 2-star symbol in Figure 7D above the comparison "HC" – Fos+" should probably be 3-stars, not 2-stars (p = 0.0006).

We have now corrected this.

26) The authors should use the less "pushy" and more conventional criteria for using star symbols for statistical significance (i.e. 1 star, p < 0.05; 2 stars, < 0.01; 3 stars < 0.001), and omit the "4 star" symbol and definition throughout. The 1-3 star criteria should be stated in the Methods (I didn't find it). The paper in its current version has only a single "4 star" symbol (Figure 3D), which should be changed to "3 star" according to the above definition, or removed (see also point 18) (see also my previous point "Data Recommendations 4").

We certainly do not intend to be pushy – In fact the statistical analysis software used (Graphpad Prism) provides “4 star” symbols. The convention followed is included in the text as delineated above.

Our provision of the detailed statistical reporting table should circumvent any confusion in this regard.

27) The bracket symbol in Figure 7F is probably wrong. On the left side, it should only point to the "HC" group, but not to both the "HC" and "Fos-" group as it stands now.

No, actually, the symbol in Figure 7 F is not wrong. The Fos+ group differs from both the HC and Fos- group. Please see the excerpt from the statistical reporting table Author response table 1.

Author response table 1.

Max Firing Rate 0.0032
HC vs. Fos- 0.3206
HC vs. Fos+ 0.003
Fos- vs. Fos+ 0.0034

28) The authors claim that "latency to fire an AP" and AHP might significantly change in re-activated ChAT neurons (Figure 7 S1E, F). However, these effects are somewhat more subtle than the other effects shown in the main Figure 7. For example, the effect of "latency" seems to be driven merely by some outlier neurons in the home cage group with very long latencies (the authors would need to show example recordings), whereas the data from the Fos- and Fos+ groups virtually overlay and the latter comparison is the more pertinent one. Thus, I think it is more cautious to drop this claim.

Indeed, as in standard practice, we show all the data and all of the results of the statistical analyses for all biophysical measures made whether “subtle” or not. These analyses consider the full range and distribution of the data without bias as to normality. Our commitment to presentation of all the features for all the recorded neurons, whether statistically significantly different or not, uses both the main and supplemental figures to provide full transparency. This is considered best practice for reporting biophysical data in general and the results of these single cell recording + post hoc PCR, labor-intensive experiments.

29) l. 427 states:

" Next, we quantified the number of cholinergic neurons expressing Fos (Fos shGFP+) after saline or predator odor exposure"

– Please explain how long the mice were exposed to the odors.

The figure legend states this:

“A. Fos-tTA/Fos-shGFP mice were placed in chambers containing a gauze pad spotted with either saline or with mountain lion urine (predator odor). Defensive behaviors were monitored for 5 min. Animals froze significantly more in the presence of predator odor than saline (Mann-Whitney, p=0.028).”

30) In Figure 8, the authors seem to analyze the results of a triple-labeling anatomical experiment; i.e. there are neurons that are Fos+; neurons that are back-labelled by fast Blue; neurons that are ChAT+, and all combinations thereof. Thus, to give the reader a better view of the various neuron pools, the authors should report the numbers of all pools of neurons (i.e. the "green", "red", "blue" neurons and all combinations thereof). This can be either done numerically or in the form of a Venn diagram. As it stands, only part of the information is given (Figure 8C, right; which seems to correspond to the percentage of back-labelled, ChAT+ neurons which were also Fos+). However, this is only part of a bigger picture. For example, the images in Figure 8C, and in Figure 8 S1A (middle, right) show also many non-ChAT+ neurons, that were re-activated in the VP/SIa area.

The purpose of this experiment, and indeed the entire study, is to investigate the contribution of different populations of basal forebrain cholinergic neurons to threat responses. At no time do we say or even imply that other neuronal types are not involved as well, but these “non-ChAT+” cells are not our focus. The requested analyses are outside the scope of this study and would not add or take away from our conclusions.

31) What does the dashed line label "Home cage" in Figure 8C, right panel, mean, and how was it derived?

The meaning of the dashed line (which we referred to as dotted) was delineated in the figure legend:

“Dotted line indicates % of Fos+ cholinergic neurons in the home cage group in this experiment.”

32) l. 442:

" Nearly the entire subset of BLA-projecting VP/SIa cholinergic neurons were also Fos-shGFP+ (Figure 8C right)."

– Another non-scientific statement. Rather, the average plus/minus S.D. of the quantification should be reported. In addition, information about all labelled neuronal populations should be given (see point 30).

Line 409 now adds: “(median 94% ± Std.dev 12.5)”.

33) The retrograde expression experiment of hM4Di from the BLA in Figure 8D is expected to drive hM4Di expression in all cholinergic BLA-projectors, not selectively in the VP/SIa. This caveat in the interpretation should be mentioned more carefully when the results of this section are summarized (lines 451 – 454).

The summary statement referred to is a conclusion based on the entirety of the mapping experiments in the manuscript. It does not refer specifically to a single figure.

Reviewer #2 (Recommendations for the authors):

I thank the authors for addressing my concerns, and I'm happy to recommend the manuscript for publication in its current form. There are a few figures that would benefit from additional labeling (outlined below).

Please add anatomical references for basal forebrain in figures F1S4, F2S1a, F4c, F5S1a, F5S1b , F6S2a.

We thank the reviewer for their positive recommendation and have added the necessary bregma annotations.

Reviewer #3 (Recommendations for the authors):

The authors have addressed many of my previous concerns and now nicely dissociate the role of NBM/SIp and VP/SIp ACh neurons in learned threat memory processing and show that NBM/SIp ACh cells are activated by memory recall and necessary for threat memory expression. However, I still have a lingering concern with their strong emphasis on the role of NBM/SIp ACh neuron (and their projections to BLA) as being memory 'engram' cells. While some of their data are consistent with this idea, they do not test key criteria for defining engram neurons, and the new data supporting a role of NBM/SIp ACh cells which project to the BLA in the threat memory 'engram' is only correlative.

Specifically, while there is evidence in support of the NBM/SIp ACh cell population (independent of where they project) in memory engrams (activated by memory cues, necessary for memory expression). There are several criteria for defining engram neurons (see Josselyn and Tonegawa Science 2020) which haven't been met including sufficiency (stimulation should produce memory in the absence of recall cues), mimicry (creation of an artificial engram) and long-lasting cellular or synaptic changes in engram neurons (here they only see excitability differences at 2.5 hr, not at long term memory timepoints). It is possible (even likely) that activity in these neurons modulates/facilitates threat memory formation and expression and are not storing the memory trace themselves.

There is even less evidence that NBM/SIp ACh cells which project to BLA being engram cells. The authors report an increase in ACh release in BLA (which isn't necessarily coming only from NBM/SIp) and added new data showing that NBM/SIp ACh cells projecting to BLA are fos activated.

While there is data supporting a potential role for these cells in engram processing, it is inappropriate to so strongly insist (in the title, abstract, and throughout) that NBM/SIp ACh neurons, and particularly those which project to BLA, are engram cells. The authors should moderate this aspect of the paper and perhaps broaden the discussion of the potential role of these cells in threat memory processing as they've already nicely done in the 'Cholinergic Modulation of Associative Threat Learning' section of the Discussion.

We thank the reviewer for their positive comments and thoughtful recommendations. We agree with the reviewer that the criteria as established by the study of glutamatergic projection neurons in the BLA and hippocampus are not completely demonstrated in our study. However, our findings do satisfy the criteria defining engram cells as mentioned in Josselyn and Tonegawa (Science, 2020). Among the 4 methods to evaluate engram cells, we have used 2 to demonstrate observation and necessity. Additionally, a similar study of oxytocinergic cells being part of a threat memory engram was recently published (Hasan MT., et al. Neuron, 2019) wherein they satisfy the same 2 criteria as our study.

Nonetheless, we have now moderated our claims regarding the presence of a cholinergic engram throughout the manuscript.

We have changed the title of the manuscript to reflect this, and modified figure titles and legends, as well as their mention in the Results section.

Finally, we have added some text to the Discussion section to reflect our thoughts on the matter:

Lines 567-595:

“Is there a Cholinergic Component in the Associative Threat Memory Engram?

Studies examining mechanisms of learning and memory in recent years have revived Semon’s theory on memory engrams: learning must result in lasting biophysical changes that form the substrate for retrieval of the learned experience (Semon 1921, Tonegawa, Liu et al. 2015). Josselyn and Tonegawa have recently updated the definition of engram cells, requiring that these be activated by learning, modified by learning, and reactivated by subsequent presentation of the recall-inducing stimuli, resulting in memory retrieval (Josselyn and Tonegawa 2020). NBM/SIp BFCNs investigated in this study indeed fulfil these criteria as they are activated by learning, show induction of Fos and altered physiological properties with recall, are reactivated by recall, and the reactivation of previously, training-activated BFCNs was necessary for recall behavior.

Multiple studies have used threat and reward learning paradigms in rodents to examine allocation of neurons to memory engrams. These studies have looked for these engram cells in regions such as cortex, amygdala and hippocampus focusing on glutamatergic pyramidal neurons (Josselyn, Köhler and Frankland 2015). However, recent work has demonstrated that memory engrams are distributed across brain-wide networks, and that reactivation of a multi-region engram more closely recapitulates natural recall behavior (Roy, Park et al. 2022).

In addition to the BLA, cholinergic neurons in the NBM/SIp region project to various limbic and sensory regions such as the lateral orbital cortex, cingulate cortex, somatosensory cortex, and mediodorsal thalamus (Ananth, Rajebhosale et al. 2023). This raises the interesting possibility that the cholinergic signaling modulates various nodes of the threat memory engram circuit in conjunction with the amygdala, allowing for coordinated retrieval of engrams across distributed networks. Such coordinated activation of distributed engrams has been recently demonstrated to more closely recapitulate natural memory retrieval (Roy, Park et al. 2022). Furthermore, functionally related regions have been shown to receive their cholinergic input from the same cholinergic nucleus (Zaborszky, Csordas et al. 2015). We propose that engram-enrolled cholinergic neurons bind distributed engrams to encode stimulus-convergent, efficient memory retrieval.”

We hope these changes will prove to be satisfactory.

Associated Data

    This section collects any data citations, data availability statements, or supplementary materials included in this article.

    Supplementary Materials

    Figure 1—source data 1. Processed ACh release fiber photometry data for Figure 1.
    Figure 1—figure supplement 2—source data 1. Processed ACh release fiber photometry data for Figure 1—figure supplement 2.
    Figure 1—figure supplement 3—source data 1. Processed ACh release data for Figure 1—figure supplement 3.
    Supplementary file 1. Statistical reporting table.

    This file contains exact sample sizes for each group, group median, 95% confidence interval of the median, actual confidence interval, statistical test used, p-values, and test statistic for each of the reported plots in the manuscript.

    elife-86581-supp1.xlsx (22.3KB, xlsx)
    MDAR checklist

    Data Availability Statement

    Source data for the fiber photometry experiments presented in Figure 1 and supplements are provided as individual source data files. Code for fiber photometry data was previously published in Crouse et al., 2020.


    Articles from eLife are provided here courtesy of eLife Sciences Publications, Ltd

    RESOURCES