Guide for reproducing our BioRxiv Manuscript.

In this post, I would like to provide guidelines on how to reproduce results and figures presented in our recent manuscript in BioRxiv. I will show how to download the raw data, run the analysis routines,  and produce figures that are shown in the paper. The text below is also present in the READ.me section of the associated Matlab toolbox that is published as an open-source in GitHub.

About

With this repository you can reproduce all the analysis and figures presented in our paper:
Fear Generalization as Threat Prediction: Adaptive Changes in Facial Exploration Strategies revealed by Fixation-Pattern Similarity Analysis. Bioarxiv
In this paper, we investigate how humans explore faces and focus how these exploration patterns change when faces are associated with an aversive outcome.

Initial Setup

You can download the data and the associated repositories, including this one from the Open Science Framework here. Add these repositories to your Matlab path.
addpath('/home/onat/Documents/Code/Matlab/FPSA_FearGen/');
addpath('/home/onat/Documents/Code/Matlab/globalfunctions//');
addpath('/home/onat/Desktop/fancycarp');

Examples

Get the list of participants.

The data you downloaded from the OSF contains all recorded participants. However, few had to be excluded for the main analysis. It is important to include all participants as it provide the possibility to reproduce the selection criteria used in the report or test the results with another selection criteria. 'get_subjects' action returns all the selected participants, totalling to 74.
>> FPSA_FearGen('get_subjects')
ans =
  Columns 1 through 45
     1     2     3     4     5     6     7     8     9    10    11    12    14    15    16    17    18    19    20    21    22    23    24    25    26    27    28    29    30    31    32    33    35    36    37    39    40    41    42    43    44    45    46    47    48
  Columns 46 through 74
    49    50    51    53    54    56    57    59    60    61    62    63    64    65    66    67    68    69    70    71    72    73    74    77    78    79    80    81    82
>> 

Sanity check 1.

In eye tracking it is usual that few trials are excluded due to blinks or bad calibration. The following code plots the number of trials per condition across the participant pool. The figure shows that almost all the particpants have their trials recorded as expected. This sanity check also shows that few participants have less trials than others.
FPSA_FearGen('get_trialcount',4)
4 above refers to the generalization phase. (2: baseline; 3: conditioning; 4: test phase).

Get the data to workspace using the get_fixmat action.

Fixmat is a compact way of storing large amounts of eye movement recordings in the form of fixation points.
>> fixmat = FPSA_FearGen('get_fixmat')
ans = 
  Fixmat with properties:

     subject: [1×118188 uint32]
       phase: [1×118188 int32]
       start: [1×118188 int32]
        stop: [1×118188 int32]
           x: [1×118188 single]
           y: [1×118188 single]
         eye: [1×118188 uint8]
    deltacsp: [1×118188 int32]
        file: [1×118188 int32]
     oddball: [1×118188 int32]
     trialid: [1×118188 int32]
         ucs: [1×118188 int32]
         fix: [1×118188 int32]
       chain: [1×118188 double]
       isref: [1×118188 double]
It stores every fixation's attributes in the form of separate vectors. For example, information such as x and y coordinates, participant's index, the image, the condition are stored in the Fixmat.
However, Fixmat variable above is not a simple Matlab structure though, but rather an instance of a Fixmat object, as coded in the FancyCarp Toolbox. The benefit of having a Fixmat object is that there are useful methods built-in in the Fixmat object, such as for example visualizing fixation density maps as heatmaps.
For example, the following code can be used to plot a fixation density map (FDM) based on all subjects (stored in S) during the generalization phase:
S = FPSA_FearGen('get_subjects');
fixmat = FPSA_FearGen('get_fixmat');
v{1} = {'subject' S 'phase' 4};
fixmat.getmaps(v{:});
fixmat.plot
The Fixmat.getmaps method creates an FDM based on the cell array argument v. In the example above,v is used to select all fixations that belong to both phase 4 and subjects S. The method Fixmat.plot plots the computed FDM. In order to create a separate FDM for different conditions or participants we would create a separate cell array for each of the required filtering conditions.
v=[];
c=0;
for ns = S([11 3 8 23]);
  c=c+1;
  v{c} = {'subject' ns 'deltacsp' 0 'phase' 4};
end
fixmat.getmaps(v{:});
fixmat.plot
This will plot a separate FDM for the 4 different participants.
Note how different participants scan faces differently. This is the basis for recent reports that analyzed the scan-path idiosyncrasy during viewing of faces. And the major reason for us to come up with the FPSA as a methodology to investigate how eye movement strategies change with aversive learning during viewing of faces.

Get FDM for the Fixation-Pattern Similarity Analysis

FPSA analysis is conducted on single-participants at a time. The get_fixmap action can be used to gather the required FDMs. It is basically a wrapper around the Fixmat.getmaps method.
>> maps = FPSA_FearGen('get_fixmap',fixmat,{'subject' 2});;
>> size(maps)
ans =
      250000          16
As you can see in maps each FDM is represented as a column vector. As we have 2 phases and 8 different conditions (faces), this amounts to 16 different vectors. This representation is appropriate for the similarity analysis as it can be readily used as an argument to Matlab's pdist

Fixation-Pattern Similarity Analysis

The following command will run a correlation-based similarity analysis on the FDMs of single participants, before (phase = 2) and after (phase =4) aversive learning. {'fix',1:100} indicates that the analysis will include 1st to 100th fixations, that is all the fixations.
1:3 ensures that similarity analysis is ran separately for the 3 runs of the generalization phase (phase = 4). This is important because the baseline phase (phase = 2) entails only one single run. In order to have a valid comparison of the similarity values across the baseline and test phases, the FDMs should contain on average same number of fixations. Computing separate FPSA matrices for each run ensures this.
sim = FPSA_FearGen('get_fpsa_fair',{'fix',1:100},1:3);
>> sim
sim = 
  struct with fields:

    correlation: [74×120 double]
The resulting matrix sim.correlation contains in each row the similarity data from a given participant. The pair-wise similarity values between the 16 FDMs are stored across the columns in a non-redundant format. You can visualize the resulting 16x16 matrix after putting the similarity values back to their positive-definite matrix format with the squareform function. 'plot_fpsa' action is used to do this.
The first quadrant of this matrix shows the similarity relationships between the 8 fixation-density maps recorded during the baseline period, before aversive learning has taken place. The second quadrant shows the same after learning has taken place.

Multi-dimensional Similarity Analysis

Another extremely useful way of representing the similarity relationships between the fixation maps consists of using the multi-dimensional scaling method. Below, I ran an MDS analysis using two dimensions.
fixmat = FPSA_FearGen('get_mdscale',squareform(mean(sim.correlation)),2)
This places a set of dots in such a way that the distances between each dot corresponds to the dissimilarity between the FDMs. Therefore, it is nice visual tool that summarizes the complex disssimilarity matrice, which is sometimes hard to digest. 
With the following piece of code, the upper part of the figure 03 of the manuscript can be produced summing up bundling previous figures.
FPSA_FearGen('figure_03A',FPSA_FearGen('get_fpsa_fair',{'fix' 1:100},1:3))

Analysis of dissimilarity matrices

One major aim in our paper was to understand how the similarity relationships between the FDMs changed with aversive learning. We take a linear modelling approach:
Y = bX
where Y is a non-redundant dissimilarity matrix, X is one of the 3 models we are testing, b are the coefficients that are to be fitted.
>> FPSA_FearGen('FPSA_get_table',{})
ans = 
      FPSA_B        FPSA_G         circle        specific      unspecific     Gaussian    subject    phase
    __________    ___________    ___________    ___________    ___________    ________    _______    _____
      -0.52828      -0.014812        0.70711     4.3298e-17        0.70711    0.48254     1          1    
     -0.064526        -0.3752    -1.0147e-17           -0.5            0.5    0.32924     1          1    
      -0.21919        -0.3426       -0.70711       -0.70711     8.6596e-17    0.25385     1          1    
      -0.28268       -0.20612             -1           -0.5           -0.5    0.32924     1          1    
      -0.35528      -0.052057       -0.70711    -1.2989e-16       -0.70711    0.48254     1          1    
       0.10652     -0.0068816    -1.7938e-16            0.5           -0.5    0.50102     1          1    
       0.27702        0.12585        0.70711        0.70711    -1.7319e-16    0.49959     1          1    
      -0.05106       -0.13439        0.70711    -4.3298e-17        0.70711    0.56373     1          1    
      -0.10373      -0.067047     6.1232e-17    -6.1232e-17     1.2246e-16      0.523     1          1    
       -0.1175       -0.39994       -0.70711    -4.3298e-17       -0.70711    0.56373     1          1  
       ...
The first two columns are the observed dissimilarity values for the first subject in the first phase (see last columns). Columns with names circlespecificunspecific and Gaussian are the modelled dissimilarity values that will be fitted for each participant separately.
In Matlab's statistical toolbox one can fit a linear model (pretty much the same way as in R) as follows:
fitlm(T,'FPSA_B ~ 1 + circle');
To to this for all the participants and models separately, we can run
T = FPSA_FearGen('FPSA_get_table',{});
[A,C] = FPSA_FearGen('FPSA_model_singlesubject',T)  
>> A.model_01
ans = 
  struct with fields:

    w1: [74×2 double]
>> A.model_01
ans = 
  struct with fields:

    w1: [74×2 double]
>> A.model_02
ans = 
  struct with fields:

    w1: [74×2 double]
    w2: [74×2 double]
The results of this analysis can be plotted with the following code, reproducing bottom panels of the figure 03.
FPSA_FearGen('figure_03C',{'fix' 1:100})

Towards a new understanding of fear generalization and its neural origin

This is an opinion article examining potential research directions to make progress on the topic of fear generalization.

This is based on the grant proposal I have written a year ago. 


Introduction

One way of dealing with the ungraspable complexity of the environment consists of making generalizations (1,2). Previously learnt regularities of the environment can be useful when applied to novel situations. For example, a novel nutriment can be categorized as inedible based on past experiences with truly harmful ones.

This competence called fear generalization (FG) is a remarkably high-level cognitive ability that builds upon more basic skills such as object recognition and categorization, statistical learning, perceptual learning, memory, affective processing and conceptual learning. FG provides an important opportunity to study how basic cognitive abilities, which are typically studied in isolation, function collectively to generate adaptive behavior in a complex world.

Notably, dissonance between these abilities manifests as maladaptive behavior and may result in mental health disorders (3–8), such as specific phobia. These are characterized by an overgeneralization of previous harmful encounters, leading to the perception of truly safe situations as harmful. Therefore, understanding the neuronal and computational mechanisms of FG is crucial both for basic, as well as clinical neuroscience.


Figure 1 Faces form a circular similarity gradient (color: distance from the CS+, see color wheel). BOLD responses form a fear tuning profile, which can be parametrically characterized16. 

The study of human FG has benefited enormously from well-established experimental paradigms dating back to Pavlov (9–14). The rationale behind these paradigms consists of characterizing how learning generalizes to other events based on their perceptual similarity with a harmful item. During conditioning, humans learn the characteristics of truly harmful (CS+) and safe (CS–) events. The harmful quality of the CS+ is established by pairing it with an aversive outcome (UCS; e.g. mild electric shock on the hand) using well-established conditioning paradigms (15), where learning can be objectively monitored.

Empirically, FG is characterized by measuring fear-related responses to other stimuli organized to form a continuous similarity gradient (Fig. 1). Typically responses decay with decreasing similarity to the CS+ resulting in graded fear tuning profiles (1). The strength of this paradigm consists of parametric characterization of behavioral and neuronal fear tuning profiles based on their peak positions and widths (16) (Fig. 1). Hence, it provides a powerful paradigm to investigate neuronal mechanisms responsible for enacting adaptive and selective fear responses.

Roadmap of Objectives for Progress in Fear Generalization
Accounts of FG differ on how they attempt to explain graded fear tuning profiles. According to different perceptual models, graded responses are a mere reflection of perceptual similarity to the behaviorally relevant stimulus (17,18). However, for FG to be an adaptive process in a complex world, it must be flexible and be regulated independently from perceptual factors. Several findings suggest that this is indeed the case;

(1) Fear tuning profiles do not always peak on the objectively most harmful stimulus19. Such peak-shift are well-documented (20,21) and indicate that FG is prone to subjective biases. 

(2) Patients with anxiety disorders typically show wider fear-tuning than healthy controls (3,4,6–8,22,23), even though they have been presented with the same perceptual stimulus material. 

(3) It has been shown that participants readily generalize to semantically related objects that are part of the same category but which do not necessarily bear close physical resemblance (e.g. a hammer and a saw) (24). 

Despite these observations, there are up to date no theoretical frameworks to understand how flexibility and adaptivity emerge during FG. With the FearGen project, I aim to go beyond the current conception of FG as a sole result of physical similarity. Compelling empirical evidence and a unifying framework will deliver ground-breaking insights on how FG can be adaptively tailored to ensure survival in a complex world. 

To realize this conceptual shift, I propose the following roadmap to be incorporated in my research agenda:

1. Descriptive to Normative Transition (WP1): Two factors presumably contribute to fear tuning profiles. (1) My previous work provided evidence that uncertainty about the occurrence of harmful events is an important factor for FG (16). However, the precise role of uncertainty in tailoring FG strategies in an adaptive manner is still largely unknown. (2) There is evidence that prior belief about the harmfulness of different stimuli along the generalization gradient plays an important role that might lead to biases in fear tuning profiles (19). The term “bias“ implicitly suggests an error in detection performance. However, this behavior can be compatible with an agent behaving optimally while integrating different sources of knowledge with the aim of disambiguating the source of harmful event in uncertain conditions. Normative models can characterize optimality of behavior. I will therefore use modern theoretical framework of predictive coding (25–27), where Bayesian inference (28–30) takes a central role, as a powerful tool to advance our understanding of FG as an optimal, adaptive phenomenon.

2. Role of categorization during fear generalization (WP2): Categorization constitutes a fundamental mechanism for organizing and transferring knowledge (31). It can therefore support FG (24). Yet, our knowledge on its contribution to FG is scarce (32). The role of categorical knowledge can be investigated in two forms. (A) Categorical knowledge can be used to transfer knowledge in abstract ways between events independently from perceptual similarity (31,33). Hence, the emergence of categorical knowledge predicts a qualitative change in fear tuning (34). This transition can be investigated in humans with representational similarity analysis (35) of multivariate activity patterns recorded in EEG and fMRI. This can inform us where, when and how abstract aversive representations form during FG. (B) Categorical knowledge can also lead to ambiguity, especially when a given item simultaneously belongs to multiple categories. Understanding how ambiguity is resolved has clinical relevance (36–38). Using the correct category for FG among many competing ones can only be achieved by collecting statistical regularities about the occurrence of harmful events. I will investigate whether this ambiguity is resolved faster in healthy humans in comparison to groups suffering from anxiety disorders.

3. Establishing causation between neural activity and fear generalization (WP3): Understanding neuronal underpinnings of FG requires ultimately establishing causality between neuronal activity and fear tuning profiles. Since most of the research about neuronal mechanisms of FG is of correlational nature, studies establishing causality can provide crucial insights. In particular, parametrically organized stimulus gradients offer the possibility to quantify effects of causal interventions by biasing FG profiles. Knowing which brain structures shape FG is not only of high interest for the progress of field but also of utmost importance for research into the aetiology of anxiety disorders and the development of clinical applications.

4. Space and time of neuronal dynamics (WP2 & 3): Neuronal activity unfolds both in time and space. In the past, investigations of FG has benefited enormously from fMRI (4,6,16,18,22), but much less so from EEG (but see 39). Consequently, our understanding of fast temporal dynamics of FG during a single trial is scarce. For a thorough characterization of rapid neuronal mechanisms of FG I will recourse to EEG and iEEG methods in addition to fMRI. Representational similarity analysis (35,40) (RSA) of multivariate activity patterns is an appropriate method to investigate FG both with EEG41 and fMRI (42). This is because FG relies ultimately on a subjective metric that evaluates the relatedness of different stimuli to CS+. Hence, similarity of activity patterns between conditions can capture subjective strategies used during FG. In combination with EEG, RSA is a powerful method (42) to investigate temporal dynamics of FG.

5. Clinical relevance (WP1 & 2 & 3): Ultimately research in FG must elucidate why anxiety disorders are associated with wide fear tuning profiles. Therefore, experiments proposed here will also be conducted in groups with anxiety disorders. For this objective, points 1, 2 and 3 on the roadmap will provide key insights: In (1), I will bring a normative approach that will give a computational account of wider fear tuning profiles. In (2), I will probe patients for a compromised strategy of updating their internal hypothesis about the source of threat. And in (3) I will conduct an intervention study elucidating causal neuronal mechanisms responsible for tailoring fear tuning.

Work Package 1: Fear Generalization as Harm Prediction: A Bayesian Integration

This work package addresses points 1 (“Normative Framework”) and 5 (“Clinical Relevance”) of the roadmap. The goal of WP1 is to cast FG as cognitive ability used to predict future harm. To do so, I conceive FG as a Bayesian inference problem for recovering the cause of threat in an uncertain environment. Here, fear tuning reflects the degree of belief about the potential harm of different stimuli. In order to reduce uncertainty about the source of threat, humans integrate different sources of available knowledge: (1) Learnt threat likelihood and (2) prior beliefs. Likelihoods reflect the probability of different stimuli to predict objectively harmful outcomes; therefore it reflects the conditioning regime imposed by the experimenter. Prior beliefs reflect one’s previous opinion (i.e. before conditioning) on different stimuli to be harmful. The integration of these sources results in the observed fear tuning, which reflects the posterior, the integrated high-level FG.

1.1 A new Bayesian framework for fear generalization

As a first step I will establish a new experimental paradigm, which explicitly manipulates uncertainty levels associated with the prediction of harmful events. This framework predicts that humans will rely more heavily on their prior beliefs when the sensory evidence for the prediction of harmful events is less reliable. This requires controlling uncertainty and using stimulus material where humans can use prior knowledge.

Controlling uncertainty: To introduce uncertainty as an experimental factor I will control UCS administration with a probability distribution along the generalization gradient during the conditioning phase. Depending on the width of the probability distribution (min: only one item along the gradient predicts harmful outcome; max: all items equally predict harmful outcome), participants will be provided with more or less reliable sensory information about the source of threat. This allows us to parametrically control uncertainty associated with the delivery of harmful outcomes and provides an elegant extension of the classical conditioning paradigms.

Introducing prior knowledge: Prior beliefs consist of the accumulated experience reflecting regularities acquired across longer time-scales in real life. Therefore, to bring prior knowledge into the game we must use stimuli that have ecological validity. Faces are an excellent choice for this objective, as they are a rich source of information in social situations. Even though prior beliefs are not directly accessible, we can observe their influence by using social priors that people commonly associate with faces. I will use gender (43), emotional expression (19), pupil size (44), ethnicity (45) and gaze direction (46). Independent evidence from fear learning literature indicates how these features modulate fear learning (e.g. males perceived more dangerous than females (43)). Here, I will bring these disparate observations about social priors in an encompassing theoretical framework. These features can be (1) manipulated parametrically, and (2) used as facial elements without interfering with the identity of a face that was previously learnt to predict UCS. By parametrically introducing prior biases along the generalization gradient, I will therefore test the predictions of the Bayesian framework at different uncertainty levels. The availability of good software support to generate faces makes this a feasible goal.

Predictions: The optimal Bayesian integration makes two clear predictions on empirical fear tuning profiles depending on how threat-likelihood and prior knowledge are aligned with each other, and the uncertainty levels in the threat likelihood. (i) With increasing uncertainty in threat-likelihood, FG will rely more strongly on prior beliefs. Hence, stronger deviations away from the objectively most harmful stimulus (i.e. where the likelihood peaks) with increasing uncertainty are expected (Fig. 2). (ii) At a given uncertainty-level, fear tuning must become increasingly sharper with increasing alignment of prior beliefs and threat-likelihood (i.e. smaller separation between the peaks of likelihood and prior). This is because the evidence from threat-likelihood and prior beliefs will add up to produce highly selective fear tuning. Learning with fear-relevant stimuli47 can be taken as an illustration of this effect. When objects that are also dangerous in real life are used as stimuli stronger and more persistent learning can be established. The predictions will be statistically tested based on the parameters (e.g peak position, tuning width) describing the empirical fear tuning profiles.


Fig. 2. Optimal Bayesian integration for FG predicts larger deviations in fear tuning with increasing uncertainty levels in the likelihood. 

Experiment: During the conditioning phase, faces along the generalization gradient will be probabilistically paired with UCSs based on a Gaussian distribution (i.e. likelihood). While the width parameter of this Gaussian controls uncertainty, its peak position determines which face predicts best UCS, hence the alignment between prior and likelihood. I will measure autonomic nervous system activity in the form of skin-conductance responses to monitor the establishment of learning. I will complement this measure with explicit ratings of UCS likelihood, which will be gathered at the end of the conditioning and test phases. During the test phase, same faces will be presented, but additionally they will exhibit facial features of social priors.

1.2 Computational characterization of neuronal fear tuning profiles

In our previous work16, we found that fear tuning in prefrontal cortex is significantly wider than in insula, which is characterized by a sharp fear tuning. The framework I outlined above allows us to understand these descriptive observations in computational terms. Wide prefrontal fear tuning is compatible with a fear representation that results from the integration of different sources of information. On the other hand, sharp insular fear tuning can reflect stimulus-UCS contingencies before the integration of prior knowledge. Hence, it is compatible with a representation based on threat-likelihood.

Once the experimental paradigm is established, I will use fMRI to obtain neuronal fear tuning profiles based on BOLD responses. Hence, testing the predictions (i) and (ii) of the Bayesian framework is a feasible objective by building on my previous expertise. It will advance our understanding of how brain forms aversive representations and whether these can be understood in terms of optimal aversive representations.

1.3. Clinical study

The Bayesian framework will also provide a basis to advance our understanding of anxiety disorders. We will use this paradigm to investigate a widespread anxiety disorder. Patients with specific phobias (ICD10 - F40.2) exhibit intense fear responses that are triggered by very specific situations47,48 and exhibit overgeneralization behavior during FG7,49. This offers an optimal scenario to study predictions from the Bayesian framework. In this framework, this type of behavior can be obtained by an inability to form optimal threat representations, or alternatively by the presence of over-precise prior beliefs. By comparing the width of the threat-likelihood and the recovered priors with healthy individuals, we will be able to identify which of these effects cause intense fear responses in phobia patients.

Work Package 2: Role of Categorization during Fear Generalization

This work package addresses points 2 (“Categorical Knowledge”), 4 (“Neurodynamics”) and 5 (“Clinical Relevance”) of the road map. The first part is concerned with the emergence of categorical knowledge during FG and the associated changes in neuronal representations and dynamics (34). In the last part, I will use hierarchically organized categories of commonly-known objects (50) to induce ambiguity and investigate FG strategies in anxiety patients and healthy controls (36–38).

2.1 Two-stage model of fear generalization

FG has been mainly investigated with stimuli organized along a similarity gradient. FG based on perceptual similarity could simply constitute one specific form. My working hypothesis is that this type of similarity-based generalization is the predecessor of category-based generalization, which instead requires an abstraction from superficial perceptual aspects. However, as long as the organism has not yet experienced enough aversive events, it will not be possible to extract features that can abstractly describe these events. At this initial stage, grouping different stimuli based on their perceptual similarity could be the best available strategy. However, gradual learning leads to the emergence of harmful and safe categories. I predict that this “Aha!” moment will be marked by a transition from fuzzy generalization profiles across perceptually similar stimuli to a binary yes-or-no type generalization profile reflecting their category membership.

To capture this transition, I will establish a novel experimental paradigm where the CS+ and CS– will characterize two probabilistic category structures31,51,52 defined across two facial features (e.g. gender and age). Across interleaved conditioning and test phases, I will give participants the possibility to extract the underlying category structure53. Importantly, I will pit category membership of faces against perceptual similarities to investigate their independent contributions over the course of the experiment. To this end, all stimuli will be characterized both (1) by their similarity to previous harmful faces, and (2) by their category membership. By modeling fear-related responses (i.e. SCR, explicit ratings) with these two predictors I will quantify the contribution of perceptual and categorical factors. I predict that with the emergence of categorical knowledge the contribution of perceptual factors will diminish. This will therefore establish an important link between two cognitive abilities that were so far studied separately.

2.2 fMRI on category learning during fear generalization

Bringing this paradigm to fMRI, I will identify neuronal mechanisms responsible for the emergence of categorical knowledge during FG. This will allow me to address the extent to which perceptual and categorical FG shares common neuronal mechanisms. As category-based FG relies on the use of more abstract knowledge, it is possible that it depends on different neuronal mechanisms than perceptual FG. This echoes an important dichotomy in the categorization research regarding abstract vs. similarity based categorical learning54. RSA40 is powerful and sensitive method that can contribute to the elucidation of neuronal mechanisms. As it evaluates between-condition similarity of multi-voxel activity patterns, it predicts different similarity geometries depending on whether FG proceeds with categorical24,55 or perceptual factors. It is therefore the appropriate tool for the identification of neuronal mechanisms of FG when multiple factors are available.

2.3 EEG on category learning during fear generalization fear generalization

Temporal dynamics of neuronal activity during FG within a single trial are largely unknown. Using the same experimental paradigm in an EEG setting, I will investigate fast neuronal dynamics (e.g time-frequency analysis) of FG and characterize temporal unfolding of perceptual and categorical factors. I will use RSA across EEG sensors41, an analysis that is appropriate for the identification of perceptual and categorical factors. Therefore, this EEG experiment will bring a different but synergistic perspective to the insights gathered in 2.2 with fMRI. In particular, using RSA we will be able to merge insights gathered from EEG and fMRI modalities42.

2.4 Overgeneralization across hierarchically organized categories

Overgeneralization in individuals with anxiety disorders has been observed in FG paradigms using perceptual gradients3,6,8,23. There is evidence that this finding can be accounted, at least in part, by perceptual confusion22. Moreover, perceptual performance is influence by learning56,57. It is therefore crucial to test the finding of overgeneralization in situations that do not require fine perceptual discrimination.

In hierarchically organized taxonomies, a single item simultaneously belongs to multiple categories (e.g. sub-ordinate, basic-level). Therefore, during a FG experiment with such stimuli, it should be impossible to unambiguously assign a CS+ to a given hierarchical level. Hence, with such stimuli one can measure overgeneralization independent of perceptual factors. In this experiment we will test FG performance of patients with specific phobias. Overgeneralization, if true, predicts that patients will consistently tune in on higher levels in the hierarchy in comparison to healthy individuals for the prediction of harmful events.

Work Package 3: Intracranial Recordings during Fear Generalization

This work package addresses point 3 (“Establishing causation”) and 4 (“Neurodynamics”) of the roadmap. Presurgical epilepsy patients with implanted deep electrodes offer a unique possibility to directly investigate activity of neuronal populations during complex cognitive tasks58,59, such as FG. These recordings will provide a detailed picture of population dynamics during FG in the form of local field potentials at precisely known neuronal sites. Most importantly, through stimulation of neuronal activity with subthreshold electrical currents60,61, we will aim to establish causality between neuronal activity and fear tuning profiles.

Electrodes are targeted to specific neuronal sites depending on each patient’s requirements. Typical locations include amygdala, hippocampus, insula and temporal cortices. Hence, even though electrodes have customized locations for every patient, the distribution of localizations is suitable for an analysis of how limbic and perceptual areas work in concert during FG.

I recently developed a paradigm to investigate the dynamic emergence of fear tuning (Fig. 3A). Pairing the CS+ face with UCS at unpredictable moments resulted in the emergence of fear tuning that dynamically grew along the course of the experiment (Fig. 3A, shock symbols). This paradigm is well suited for introducing biases in fear tuning profiles during learning with subthreshold electrical stimulation. Hence, my know-how from this previous fMRI experiment will contribute collecting best quality data with an already tested paradigm.

3.1 Dynamic emergence of fear generalization profiles in affective brain structures

Fear tuning in human brain has been almost exclusively shown with fMRI4,16,18,19. Given the sluggish nature of BOLD responses, our knowledge on neuronal signatures responsible for encoding fear responses is largely unknown. Local field potentials capture population dynamics at high temporal resolution, and thus provide a great source of information. The parametric nature of the FG experiment makes it possible to identify fear tuning across different frequency channels. This will provide important insights for understanding neuronal mechanisms of FG.

3.2 Establishing causality

Understanding how neuronal activity causally contributes to the regulation of behavioral fear tuning profiles is crucial for an understanding of neuronal mechanisms implicated in anxiety disorders62. Causal intervention through electrical stimulation is a powerful method that can be used to investigate neuronal sites that can potentially bias fear tuning profiles61. The parametric nature of fear tuning profiles provides an objective and quantitative method to investigate these biases at the behavioral level.
Fig. 3 A. Emergence of FG across time and faces in Insula. (color: BOLD amplitude, shock symbols: UCS delivery with CS+ presentation). FG is shown only for trials where no UCS is administrated. B. Causal intervention with electric stimulation biases fear tuning in a reversible manner. C. Deep implanted electrodes targeting amygdala. 

To achieve this objective, I will aim to introduce biases on fear tuning profiles via subthreshold electrical stimulation. I will use the loudness of white noise auditory bursts as UCSs with presurgical patients. Therefore along the generalization gradient faces will be paired with increasing loudness levels. The CS+ face will be paired with the loudest UCS. I will aim to increment the aversive quality of faces closely neighboring the CS+ face with electrical stimulation in a reversible manner across two different runs (Fig. 3B). For stimulation, we will use electrode contacts that are functionally related to FG, which will be characterized previously. Using this methodology I will investigate the causal contribution of different neuronal sites to the production of fear tuning profiles.

References


1. Shepard RN. Toward a universal law of generalization for psychological science. Science. 1987;237(4820):1317-1323.
2. Tenenbaum JB, Griffiths TL. Generalization, similarity, and Bayesian inference. Behav Brain Sci. 2001;24(4):629–640. doi:10.1017/S0140525X01000061.
3. Lissek S, Rabin S, Heller RE, et al. Overgeneralization of conditioned fear as a pathogenic marker of panic disorder. Am J Psychiatry. 2009;167(1):47-55. doi:10.1176/appi.ajp.2009.09030410.
4. Greenberg T, Carlson JM, Cha J, Hajcak G, Mujica-Parodi LR. Ventromedial Prefrontal Cortex Reactivity Is Altered in Generalized Anxiety Disorder During Fear Generalization. Depress Anxiety. 2013;30(3):242-250. doi:10.1002/da.22016.
5. Lissek S, Kaczkurkin AN, Rabin S, Geraci M, Pine DS, Grillon C. Generalized Anxiety Disorder Is Associated With Overgeneralization of Classically Conditioned Fear. Biol Psychiatry. 2014;75(11):909-915. doi:10.1016/j.biopsych.2013.07.025.
6. Cha J, Greenberg T, Carlson JM, DeDora DJ, Hajcak G, Mujica-Parodi LR. Circuit-Wide Structural and Functional Measures Predict Ventromedial Prefrontal Cortex Fear Generalization: Implications for Generalized Anxiety Disorder. J Neurosci. 2014;34(11):4043-4053. doi:10.1523/JNEUROSCI.3372-13.2014.
7. Dymond S, Dunsmoor JE, Vervliet B, Roche B, Hermans D. Fear Generalization in Humans: Systematic Review and Implications for Anxiety Disorder Research. Behav Ther. http://www.sciencedirect.com/science/article/pii/S0005789414001348. Accessed March 2, 2015.
8. Jasnow AM, Lynch JF, Gilman TL, Riccio DC. Perspectives on fear generalization and its implications for emotional disorders. J Neurosci Res. 2017;95(3):821-835. doi:10.1002/jnr.23837.
9. Anrep GV. The irradiation of conditioned reflexes. Proc R Soc Lond Ser B Contain Pap Biol Character. 1923:404–426.
10. Bass MJ, Hull CL. The irradiation of a tactile conditioned reflex in man. J Comp Psychol. 1934;17(1):47-65.
11. Hovland CI. The Generalization of Conditioned Responses: I. The Sensory Generalization of Conditioned Responses with Varying Frequencies of Tone. J Gen Psychol. 1937;17(1):125-148. doi:10.1080/00221309.1937.9917977.
12. Guttman N, Kalish HI. Discriminability and stimulus generalization. J Exp Psychol. 1956;51(1):79-88.
13. Hull CL. The problem of primary stimulus generalization. Psychol Rev. 1947;54(3):120-134. doi:10.1037/h0061159.
14. Lashley KS, Wade M. The Pavlovian theory of generalization. Psychol Rev. 1946;53(2):72-87.
15. Lonsdorf TB, Menz MM, Andreatta M, et al. Don’t fear “fear conditioning”: Methodological considerations for the design and analysis of studies on human fear acquisition, extinction, and return of fear. Neurosci Biobehav Rev. 2017;77:247-285. doi:10.1016/j.neubiorev.2017.02.026.
16. Onat S, Büchel C. The neuronal basis of fear generalization in humans. Nat Neurosci. 2015;advance online publication. doi:10.1038/nn.4166.
17. Lissek S. Toward an account of clinical anxiety predicated on basic, neurally mapped mechanisms of Pavlovian fear-learning: the case for conditioned overgeneralization. Depress Anxiety. 2012;29(4):257–263. doi:10.1002/da.21922.
18. Lissek S, Bradford DE, Alvarez RP, et al. Neural substrates of classically conditioned fear-generalization in humans: a parametric fMRI study. Soc Cogn Affect Neurosci. 2014;9(8):1134-1142. doi:10.1093/scan/nst096.
19. Dunsmoor JE, Prince SE, Murty VP, Kragel PA, LaBar KS. Neurobehavioral mechanisms of human fear generalization. NeuroImage. 2011;55(4):1878-1888. doi:10.1016/j.neuroimage.2011.01.041.
20. Purtle RB. Peak shift: A review. Psychol Bull. 1973;80(5):408-421.
21. Thomas DR, Mood K, Morrison S, Wiertelak E. Peak shift revisited: a test of alternative interpretations. J Exp Psychol Anim Behav Process. 1991;17(2):130-140.
22. Laufer O, Israeli D, Paz R. Behavioral and Neural Mechanisms of Overgeneralization in Anxiety. Curr Biol. 2016;26(6):713-722. doi:10.1016/j.cub.2016.01.023.
23. Dunsmoor JE, Paz R. Fear Generalization and Anxiety: Behavioral and Neural Mechanisms. Biol Psychiatry. 2015;78(5):336-343. doi:10.1016/j.biopsych.2015.04.010.
24. Dunsmoor JE, Kragel PA, Martin A, LaBar KS. Aversive Learning Modulates Cortical Representations of Object Categories. Cereb Cortex. 2014;24(11):2859-2872. doi:10.1093/cercor/bht138.
25. Rao RP, Ballard DH. Predictive coding in the visual cortex: a functional interpretation of some extra-classical receptive-field effects. Nat Neurosci. 1999;2(1):79–87. doi:10.1038/4580.
26. Friston K. The free-energy principle: a unified brain theory? Nat Rev Neurosci. 2010;11(2):127-138. doi:10.1038/nrn2787.
27. Vilares I, Kording K. Bayesian models: the structure of the world, uncertainty, behavior, and the brain. Ann N Y Acad Sci. 2011;1224(1):22-39. doi:10.1111/j.1749-6632.2011.05965.x.
28. Knill DC, Pouget A. The Bayesian brain: the role of uncertainty in neural coding and computation. TRENDS Neurosci. 2004;27(12):712–719.
29. Ernst MO, Bulthoff HH. Merging the senses into a robust percept. Trends Cogn Sci. 2004;8(4):162–169. doi:10.1016/j.tics.2004.02.002.
30. Koerding KP, Wolpert DM. Bayesian integration in sensorimotor learning. Nature. 2004;427(6971):244–247. doi:10.1038/nature02169.
31. Ashby FG, Maddox WT. Human Category Learning. Annu Rev Psychol. 2005;56(1):149-178. doi:10.1146/annurev.psych.56.091103.070217.
32. Dunsmoor JE, Murphy GL. Categories, concepts, and conditioning: how humans generalize fear. Trends Cogn Sci. 2015;19(2):73-77. doi:10.1016/j.tics.2014.12.003.
33. Medin DL. Concepts and conceptual structure. Am Psychol. 1989;44(12):1469.
34. Ohl FW, Scheich H, Freeman WJ. Change in pattern of ongoing cortical activity with auditory category learning. Nature. 2001;412(6848):733-736. doi:10.1038/35089076.
35. Kriegeskorte N, Mur M, Ruff DA, et al. Matching Categorical Object Representations in Inferior Temporal Cortex of Man and Monkey. Neuron. 2008;60(6):1126-1141. doi:10.1016/j.neuron.2008.10.043.
36. Nader K, Balleine B. Ambiguity and anxiety: when a glass half full is empty. Nat Neurosci. 2007;10(7):807-808. doi:10.1038/nn0707-807.
37. Grupe DW, Nitschke JB. Uncertainty and anticipation in anxiety: an integrated neurobiological and psychological perspective. Nat Rev Neurosci. 2013;14(7):488-501. doi:10.1038/nrn3524.
38. Lake JI, LaBar KS. Unpredictability and Uncertainty in Anxiety: A New Direction for Emotional Timing Research. Front Integr Neurosci. 2011;5. doi:10.3389/fnint.2011.00055.
39. McTeague LM, Gruss LF, Keil A. Aversive learning shapes neuronal orientation tuning in human visual cortex. Nat Commun. 2015;6:7823. doi:10.1038/ncomms8823.
40. Kriegeskorte N, Mur M, Bandettini P. Representational Similarity Analysis – Connecting the Branches of Systems Neuroscience. Front Syst Neurosci. 2008;2. doi:10.3389/neuro.06.004.2008.
41. Kietzmann TC, Gert AL, Tong F, König P. Representational Dynamics of Facial Viewpoint Encoding. J Cogn Neurosci. 2017;29(4):637-651. doi:10.1162/jocn_a_01070.
42. Cichy RM, Pantazis D, Oliva A. Resolving human object recognition in space and time. Nat Neurosci. 2014;17(3):455-462. doi:10.1038/nn.3635.
43. Navarrete CD, Olsson A, Ho AK, Mendes WB, Thomsen L, Sidanius J. Fear Extinction to an Out-Group Face: The Role of Target Gender. Psychol Sci. 2009;20(2):155-158. doi:10.1111/j.1467-9280.2009.02273.x.
44. Demos KE, Kelley WM, Ryan SL, Davis FC, Whalen PJ. Human Amygdala Sensitivity to the Pupil Size of Others. Cereb Cortex. 2008;18(12):2729-2734. doi:10.1093/cercor/bhn034.
45. Olsson A, Ebert JP, Banaji MR, Phelps EA. The role of social groups in the persistence of learned fear. Science. 2005;309(5735):785-787. doi:10.1126/science.1113551.
46. N’Diaye K, Sander D, Vuilleumier P. Self-relevance processing in the human amygdala: gaze direction, facial expression, and emotion intensity. Emotion. 2009;9(6):798.
47. Öhman A, Mineka S. Fears, phobias, and preparedness: toward an evolved module of fear and fear learning. Psychol Rev. 2001;108(3):483.
48. Ipser JC, Singh L, Stein DJ. Meta-analysis of functional brain imaging in specific phobia. Psychiatry Clin Neurosci. 2013;67(5):311-322. doi:10.1111/pcn.12055.
49. Dymond S, Schlund MW, Roche B, Whelan R. The spread of fear: Symbolic generalization mediates graded threat-avoidance in specific phobia. Q J Exp Psychol. 2014;67(2):247-259. doi:10.1080/17470218.2013.800124.
50. Xu F, Tenenbaum JB. Word learning as Bayesian inference. Psychol Rev. 2007;114(2):245.
51. Nosofsky RM. Attention, similarity, and the identification–categorization relationship. J Exp Psychol Gen. 1986;115(1):39.
52. Sigala N, Logothetis NK. Visual categorization shapes feature selectivity in the primate temporal cortex. Nature. 2002;415(6869):318-320. doi:10.1038/415318a.
53. Kietzmann TC, König P. Perceptual learning of parametric face categories leads to the integration of high-level class-based information but not to high-level pop-out. J Vis. 2010;10(13):20-20. doi:10.1167/10.13.20.
54. Sloman SA. The empirical case for two systems of reasoning. Psychol Bull. 1996;119:3–22.
55. Visser RM, Scholte HS, Beemsterboer T, Kindt M. Neural pattern similarity predicts long-term fear memory. Nat Neurosci. 2013;16(4):388-390. doi:10.1038/nn.3345.
56. Resnik J, Sobel N, Paz R. Auditory aversive learning increases discrimination thresholds. Nat Neurosci. 2011;14(6):791-796. doi:10.1038/nn.2802.
57. Holt DJ, Boeke EA, Wolthusen RPF, Nasr S, Milad MR, Tootell RBH. A parametric study of fear generalization to faces and non-face objects: relationship to discrimination thresholds. Front Hum Neurosci. 2014;8:624. doi:10.3389/fnhum.2014.00624.
58. Lachaux J-P, Axmacher N, Mormann F, Halgren E, Crone NE. High-frequency neural activity and human cognition: past, present and possible future of intracranial EEG research. Prog Neurobiol. 2012;98(3):279-301. doi:10.1016/j.pneurobio.2012.06.008.
59. Vidal JR, Perrone-Bertolotti M, Kahane P, Lachaux J-P. Intracranial spectral amplitude dynamics of perceptual suppression in fronto-insular, occipito-temporal, and primary visual cortex. Front Psychol. 2015;5. doi:10.3389/fpsyg.2014.01545.
60. Cogan SF. Neural Stimulation and Recording Electrodes. Annu Rev Biomed Eng. 2008;10(1):275-309. doi:10.1146/annurev.bioeng.10.061807.160518.
61. Guillory SA, Bujarski KA. Exploring emotions using invasive methods: review of 60 years of human intracranial electrophysiology. Soc Cogn Affect Neurosci. 2014;9(12):1880-1889. doi:10.1093/scan/nsu002.
62. Etkin A. Neurobiology of Anxiety: From Neural Circuits to Novel Solutions? Depress Anxiety. 2012;29(5):355-358. doi:10.1002/da.21957.

Export SPM results table to spreadsheet for publication

When using SPM (Statistical Parametric Mapping) for the analysis of fMRI data, one soon or later needs to export the numbers in the "SPM Results Table" to a spreadsheet format to include in a publication. While the SPM Table is neatly organized within a Matlab Figure, it is not easy to export the numbers to a friendly format.

SPM results table.
Here is a small Matlab function that simplifies the process. It neatly dumps all values in an SPM table into a .txt file.
Same table in the text editor.

You can then easily import this text file to your spreadsheet engine (excel, gdoc, etc.) as a tab separated .txt file. With little cosmetic efforts, it is ready to be inserted in your next publication.

Great Night at the Night of Science in Hamburg 2017

During the Night of Science event in Hamburg, we (me, Lea Kampermann and Lukas Neugebauer) introduced the eye-tracking technique to our guests, and illustrated it with the classical change blindness experiment.

We explained the basics of the eye-tracking and illustrated it with a classical experiment in visual neurosciences, namely the phenomenon of change blindness. We received about 60 people, and recorded eye-movements from 10 volunteers. Below I prepared an animated GIF that shows both the images shown to volunteers and the location that are most fixated by one volunteers. Overall our efforts were rewarded as being the highest rated demo during the night in our department.

Nine images shown to ten volunteers during the Night of Science event in Hamburg. The animation above consists of 3 different images.  The first two consist of flickering images presented during the experiment to induce change blindness; and the last one is the semi-transparent fixation map showing locations that were most attended by all volunteers.

For more information on change blindness, there is probably no better source than Kevan O'Regan's webpage whose name is closely associated with this phenomenon. For preparing the demo, we actually used many of the images that were used in the original publication. Furthermore, his webpage provides also a rich and original source of information on vision and perception.

Just submitted the following short paper to the 1st Cognitive Computational Neuroscience meeting. I am looking forward to participate!

Model-Based Fixation-Pattern Similarity Analysis Reveals Adaptive Changes in Face-Viewing Strategies Following Aversive Learning

Lea Kampermann (l.kampermann@uke.de)
Department of Systems Neuroscience
University Medical Center Hamburg-Eppendorf

Niklas Wilming (n.wilming@uke.de)
Department of Neurophysiology and Pathophysiology
University Medical Center Hamburg-Eppendorf

Arjen Alink (a.alink@uke.de)
Department of Systems Neuroscience
University Medical Center Hamburg-Eppendorf

Christian Büchel (c.buechel@uke.de)
Department of Systems Neuroscience
University Medical Center Hamburg-Eppendorf

Selim Onat (sonat@uke.de)
Department of Systems Neuroscience
University Medical Center Hamburg-Eppendorf


Abstract:


Learning to associate an event with an aversive outcome typically leads to generalization when similar situations are encountered. In real-world situations, generalization must be based on the sensory evidence collected through active exploration. However, our knowledge on how exploration can be adaptively tailored during generalization is scarce. Here, we investigated learning-induced changes in eye movement patterns using a similarity-based multivariate fixation-pattern analysis. Humans learnt to associate an aversive outcome (a mild electric shock) with one face along a circular perceptual continuum, whereas the most dissimilar face on this continuum was kept neutral. Before learning, eye-movement patterns mirrored the similarity characteristics of the stimulus continuum, indicating that exploration was mainly guided by subtle physical differences between the faces. Aversive learning increased the dissimilarity of exploration patterns. In particular, this increase occurred specifically along the axis separating the shock predicting face from the neutral one. We suggest that this separation of patterns results from an internal categorization process for the newly learnt harmful and safe facial prototypes.
Keywords: Eye movements; Generalization; Categorization; Face Perception; Aversive Learning; Multivariate Pattern Analysis; Pattern Similarity

To avoid costly situations, animals must be able to rapidly predict future adversity based on actively harvested information from the environment. In humans, a central part of active exploration involves eye movements, which can rapidly determine what information is available in a scene. However, we currently do not know the extent to which eye movement strategies are flexible and can be adaptive following aversive learning.

We investigated how aversive learning influences exploration strategies during viewing of faces that were designed to form a circular perceptual continuum (Fig. 1A). One randomly chosen face along this continuum (CS+; Fig. 1, red, see colorwheel) was paired with a mild electric shock, which introduced an adversity gradient based on  physical  similarity  to  the


Figure 1: (A) 8 exploration patterns (FDMs, colored frames) from a representative individual overlaid on 8 face stimuli (numbered 1 to 8) calibrated to span a circular similarity continuum across two dimensions. A pair of maximally dissimilar faces was randomly selected as CS+ (red border) and CS– (cyan border; see color wheel for color code). The similarity relationships among the 8 faces and the resulting exploration patterns are depicted as two 8×8 matrices. (B-E) Multidimensional-scaling representations (top row) and the corresponding dissimilarity matrices (bottom row) depicting four possible scenarios on how learning could change the similarity geometry between the exploration maps (same color scheme; red: CS+; cyan: CS–). These matrices are decomposed onto covariate components (middle row) centered either on the CS+/CS– (specific component) or +90°/–90° faces (unspecific component). A third component is uniquely centered on the CS+ face (adversity component). These components were fitted to the observed dissimilarity matrices, and model selection procedure was carried out.

CS+ face. The most dissimilar face (CS–; Fig. 1, cyan) separated by 180° on the circular continuum was not reinforced and thus stayed neutral. Using this paradigm, we were able to investigate how exploration strategies were modified by both the physical similarity relationships between faces, and the adversity gradient introduced through aversive learning.

We used a variant of representational similarity analysis (Kriegeskorte, Mur, & Bandettini, 2008) that we term “fixation-pattern similarity analysis” (FPSA). FPSA considers exploration patterns as multivariate entities and assesses between-condition dissimilarity of fixation patterns for individual participants (Fig. 1A). We formulated 4 different hypotheses (Bottom-up saliency, increased arousal, adversity categorization, adversity tuning) based on how aversive learning might alter the similarity relationships between exploration patterns when one face on the continuum started to predict adversity (Fig. 1B-E).

Before learning, eye movement patterns mirrored the similarity characteristics of the stimulus continuum, indicating that exploration was mainly guided by subtle physical differences between the faces. Aversive learning resulted in a global increase in dissimilarity of eye movement patterns following learning. Model-based analysis of the similarity geometry indicated that this increase was specifically driven by a separation of patterns along the adversity gradient, in agreement with the adversity categorization model (Fig. 1D). These findings show that aversive learning can introduce substantial remodeling of exploration patterns in an adaptive manner during viewing of faces. In particular, we suggest that separation of patterns for harmful and safe prototypes results from an internal categorization process operating along the perceptual continuum following learning.

References

Kriegeskorte, N., Mur, M., & Bandettini, P. (2008). Representational Similarity Analysis – Connecting the Branches of Systems Neuroscience. Frontiers in Systems Neuroscience, 2. https://doi.org/10.3389/neuro.06.004.2008

 



Talk given at the EMHFC Conference


I gave this talk at the European Meeting on Human Fear Conditioning about "Temporal Dynamics of Aversive Generalization".

Abstract:
Temporal Dynamics of Aversive Learning and Generalization in Amygdala
The amygdala is thought to orchestrate coordinated bodily responses important for the survival
of the organism during threatening situations. However its contribution to the generalization of
previously learnt aversive associations is not well understood. As amygdala responses in the
context of fear conditioning are temporally phasic and adapt quickly, we designed a new
paradigm to investigate its temporal dynamics during fear generalization. We used faces that
formed a perceptual circular similarity continuum, allowing us to gather two-sided generalization
gradients. While one face predicted an aversive outcome (UCS), the most dissimilar face was
kept neutral. Importantly, participants were compelled to learn these associations throughout the
fMRI recording which they started naive, allowing us to collect temporally resolved
generalization gradients for BOLD and skin-conductance responses. Following fMRI, we
evaluated subjective likelihood for single faces to be associated with UCS, and complemented
these with behavioral measurements using eye-movement recordings to assess how the saliency
associated with faces were modified. Aversive generalization in the amygdala emerged late
during the task, and temporal dynamics were characterized by low learning rates. We observed
significant differences in amygdala responses for participants who exhibited a behavioral effect
in addition to verbal ratings of UCS likelihood. Amygdalar responses contrasted with temporal
dynamics in the insula where generalization gradients emerged earlier and gradually increased
with higher learning rates, similar to skin-conductance responses. Overall our results imply a
weak and late contribution of the amygdala to aversive generalization, in comparison to insular
responses that are stronger and contribute early during learning.