Single unit recording of aIC neurons during conditioned fear expression.

(A) Schematic of the auditory fear conditioning and single unit recording procedures. Neuronal activities and behaviors were simultaneously recorded during the expression of conditioned fear. The control group underwent the same procedures except for fear conditioning.

(B) Histological sample showing the tip of a recording electrode, indicated with a red arrow. Brain atlas maps indicating the locations of aIC recording sites used for analysis. The black dots represent the recording site. GI: granular insular cortex, DI: dysgranular insular cortex, AID: dorsal agranular insular cortex, AIV: ventral agranular insular cortex.

(C) Sample cluster cutting of neuronal recording acquired from a tetrode.

(D) Juxtaposed action potential shape of all recorded neurons (left). Putative pyramidal neurons and interneurons were distinguished based on the shape of the action potential. The plot on the right shows action potential duration and repolarization. Putative pyramidal neurons are represented by green dots, while putative interneurons are represented by gray dots.

(E) Freezing (fear) behavior of conditioned or control mice during five tone presentations (shaded blue bars).

(F) Freezing of conditioned (left) or control (right) mice during five tone-on and tone-off sessions.

(G) Normalized (z-scored) firing rate of aIC pyramidal neurons simultaneously recorded with freezing (shaded blue bars represent tone presentations).

(H) Firing rate of conditioned (left) or control (right) mice during five tone-on and tone-off sessions. Paired sample t-test was used to compare means between tone-on and tone-off sessions in each group (E-H). *** P<0.001

(I) Spectrogram of aIC local field potential recording of the conditioned and control groups before and during tone presentation.

(J) Power spectrum of aIC local field potential of the conditioned and control groups.

(K) Theta power analysis before tone (left) and during tone (right) of the conditioned and control groups. Cond.: conditioned group, Cont.: control group. The Mann–Whitney U test was used to assess the statistical difference between the conditioned and control groups in each session. *** P<0.001

The activity of aIC pyramidal neurons bidirectionally correlates with freezing (fear) behaviors.

(A) Example of aIC neurons that decreased (inhibited-cell, left), increased (excited-cell, center), or did not change (no-change cell, right) their firing rate over five CS presentations (CS1-5). Top row: raster plot, each tick mark indicates when an action potential occurred. Middle row: peri-stimulus time histogram of each cell type over five CS presentations. Bottom row: plot of neuronal activity changes and freezing behavior of each cell type. Blue line: activity of inhibited cell, Red line: activity of excited cell, Gray line: activity of no-change cell, Purple line: freezing behavior.

(B) Distribution of Pearson’s correlation coefficient values between freezing behavior and aIC neuronal activities of the fear conditioned group. The red bars represent cells with significant positive correlation (freezing-excited cells), blue bars represent cells with significant negative correlation (freezing-inhibited cells), and gray bars represent cells with no significant correlation (non-responsive) with freezing behavior. Significance was determined at P<0.05. Relative ratio of the non-responsive and responsive (freezing-excited or freezing-inhibited) cell types recorded in the fear conditioned group (right pie chart).

(C) Locations of cell-types recorded. Red dot: freezing-excited or freezing-excited/non-responsive cells, blue dot: freezing-inhibited or freezing-inhibited/non-responsive cell, red and blue gradation dot: combination of freezing-excited/freezing-inhibited/non-responsive cells, gray dot: non-responsive cells.

(D) Distribution of all recorded aIC neuronal activities (z-scored) before, during, and after CS (tone) presentation. Each line represents the mean response of one cell to five CSs, aligned from the lowest to the highest firing rate change in response to CS. The blue box represents cells classified as freezing-inhibited cells, and the red box represents cells classified as freezing-excited cells.

(E) Baseline firing rates of freezing-excited cells, freezing-inhibited cells, non-responsive cells, and control cells. The Mann-Whitney U test was used to assess statistical differences between responsive types (** P<0.01, *** P<0.001).

(F) Pearson correlation analysis for each neuronal response type and freezing behavior.

(G) Firing rate changes of freezing-inhibited and freezing-excited cells analyzed in 1 s intervals (left). The box plot shows the difference in time that each cell type takes to reach peak changes. The Mann–Whitney U test was used to assess the statistical difference between the time that freezing-inhibited and freezing-excited cells take to reach peak firing rate change. * P<0.05

(H) Tone onset and offset analysis of the conditioned and the control groups. Blue line: freezing-inhibited, red line: freezing-excited, gray line: non-responsive cells. Mean (± standard deviation) spikes per second.

(I) Schematic drawing of the conditioned fear extinction procedure.

(J) Freezing (fear) behavior of mice during conditioned fear extinction (ten tone presentations, N=7 mice).

(K) Distribution of Pearson’s correlation coefficient values between freezing behavior and aIC neuronal activities. The red bars represent cells with significant positive correlation (freezing-excited cells), blue bars represent cells with significant negative correlation (freezing-inhibited cells), and gray bars represent cells with no significant correlation with freezing behavior. Significance was determined at P<0.05.

(L) Pearson correlation analysis for each neuronal response type and freezing behavior.

(M) Distribution of Pearson’s correlation coefficient values between tone and aIC neuronal activities.

Optogenetic regulation of aIC pyramidal neuron output modulates conditioned fear expression.

(A) Schematic drawing of the experimental protocol. Experimental groups consisted of viral vectors injected for excitation (AAV-CaMKII-ChR2-eYFP), inhibition (AAV-CaMKII-NpHR3-eYFP), and control (AAV-CaMKII-eYFP).

(B) Histological sample showing viral expression pattern and optic fiber placement. Dashed white line depicts location of optic fiber placement. Yellow dotted square indicates the 40x magnified region shown below.

(C) Velocity before (baseline), during (light on), and after (light off) optostimulation of the three groups in an open field test.

(D) Simultaneous single unit and local field potential (LFP) recording during the 20 Hz blue-light activation of ChR2-expressing neurons. Juxtaposed unit activity (each tick indicates an action potential), peri-stimulation time histogram, and LFP (left). Local brain rhythm (frequency) changes induced by blue light stimulation (right).

(E) Simultaneous single unit and local field potential (LFP) recording during the continuous (cont.) yellow-light activation of NpHR-expressing neurons. Juxtaposed unit activity (each tick indicates an action potential), peri-stimulation time histogram, and LFP (left). Local brain rhythm (frequency) changes induced by yellow light stimulation (right).

(F) Mean (± SEM) freezing behavior of mice during the five CS (tone) and optostimulation delivery.

(G) Freezing behavior of the five tone-on with optostimulation and tone-off periods (Kruskal–Wallis test; **P<0.01, ***P<0.001).

(H) Comparison between the first and the last tone+optostimulation and tone-off session freezing behavior among groups (Wilcoxon signed ranks test; **P<0.01).

(I) Mean (± SEM) freezing behavior of mice during five CS-only presentations without any optostimulation.

(J) Freezing behavior of the five tone-on and tone-off periods.

(K) Comparison between the first and the last tone and tone-off session freezing behavior among groups.

(F-K) Blue line represents optical excitation (ChR2, n=13), yellow line represents optical inhibition (NpHR3, n=7), and gray line represents control (eYFP, n=13).

(F, G, L, J) The Kruskal–Wallis test was used to test statistical significance of among groups (H and K) The Wilcoxon signed ranks test was used to compare means of freezing behavior between first and last sessions in each group.

Tracing of aIC neurons that projects to the thalamus or the amygdala.

(A) Schematic drawing of cholera toxin B (CTB) retrograde tracers injected into the amygdala (green) and thalamus (red).

(B) Image showing the injection sites.

(C) High-resolution of the panel b image. sm: stria medullaris of the thalamus, LHb: lateral habenular nucleus, MHb: medial habenular nucleus, DG: dentate gyrus, PV: paraventricular thalamic nucleus, IMD: intermediodorsal thalamic nucleus, MDM: mediodorsal thalamic nucleus, medial part, MDC: mediodorsal thalamic nucleus, central part, MDL: mediodorsal thalamic nucleus, lateral part, CM: central medial thalamic nucleus, LaDL: lateral amygdaloid nucleus, dorsolateral part, and BLA: basolateral amygdaloid nucleus, anterior part.

(D) Sample image of neurons in the aIC that either project to the thalamus (red) or the amygdala (green). A high-resolution image of the thalamus-projecting aIC neurons is shown on the left, and that of the amygdala-projecting aIC neurons is shown on the right. M2: secondary motor cortex, M1: primary motor cortex, S1: primary somatosensory cortex, GI: granular insular cortex, DI: dysgranular insular cortex, AID: agranular insular cortex, dorsal part, AIV: agranular insular cortex, ventral part, Pir: piriform cortex, Den: dorsal endopiriform nucleus, CPu: caudate putamen (striatum), and AcbC: accumbens nucleus, core.

(E) The number of aIC neurons that were labeled to either project to the thalamus (Alexa fluor 594, n=3 mice) or the amygdala (Alexa fluor 488, n=3 mice). Data are presented as mean ± SEM. The Mann-Whitney U test was used to assess statistical difference (* P<0.05).

Specific optogenetic activation of aIC projection to the medial thalamus or the amygdala bidirectionally regulate conditioned fear expression.

(A) Schematic drawing of the adeno-associated virus (AAV) injection and optostimulation sites. Viral vectors for opsin expression (AAV-CaMKII-ChR2-eYFP) or control (AAV-CaMKII-eYFP) were bilaterally injected in the aIC, and bilateral optostimulation was delivered in the mid-thalamic region.

(B) Example expression pattern of aIC projection terminals in the thalamic region and the location of bilateral optic fiber placement outlined by dashed lines. LHb: lateral habenular nucleus, MHb: medial habenular nucleus, PV: paraventricular thalamic nucleus, IMD: intermediodorsal thalamic nucleus, MDM: mediodorsal thalamic nucleus, medial part, MDC: mediodorsal thalamic nucleus, central part, MDL: mediodorsal thalamic nucleus, lateral part, CM: central medial thalamic nucleus.

(C) Mean (±SEM) freezing behavior of mice expressing ChR2 (n=12) or control (eYFP, n=8) induced by optostimulation of aIC projection terminating in the thalamus to repeated CS (tone). Blue bars indicate when tone and optostimulation were delivered. Last point indicates freezing behavior when tone was delivered without optostimulation 24h after the extinction protocol (retrieval).

(D) Freezing behavior during the “Tone on+optostim” and “Tone off” of the aIC thalamus ChR2 and control groups.

(E) Schematic drawing of AAV injection and optostimulation sites. Viral vectors for opsin expression (AAV-CaMKII-ChR2-eYFP) or control (AAV-CaMKII-eYFP) were bilaterally injected in the aIC, and bilateral optostimulation was delivered in the anterior-amygdala region.

(F) Example expression pattern of aIC projection terminals in the amygdala and location of bilateral optic fiber placement outlined by dashed lines. CeA: central nucleus of the amygdala, BLA: basolateral amygdaloid nucleus, anterior part.

(G) Mean (± SEM) freezing behaviors of mice expressing ChR2 (n=6) or control (eYFP, n=6) induced by optostimulation of aIC projection terminating in the amygdala to repeated CS (tone). Blue bars indicate when tone and optostimulation were delivered. Last point indicates freezing behavior when tone was delivered without optostimulation 24h after the extinction protocol (retrieval).

(H) Freezing behavior during the “Tone on+optostim” and “Tone off” of the aIC-amygdala ChR2 and control groups.

(C and G) Two-way repeated measure ANOVAs were performed to compare freezing behavior between groups (C and G).

(D and H) Paired t-test was used to compare means of “Tone on+optostim” and “Tone off” in each group. (*P<0.05, **P<0.01, ***P<0.001).