Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Behavioral/Cognitive

Hippocampal Representation of Threat Features and Behavior in a Human Approach–Avoidance Conflict Anxiety Task

Aslan Abivardi, Saurabh Khemka and Dominik R. Bach
Journal of Neuroscience 26 August 2020, 40 (35) 6748-6758; DOI: https://doi.org/10.1523/JNEUROSCI.2732-19.2020
Aslan Abivardi
1Computational Psychiatry Research, Department of Psychiatry Psychotherapy and Psychosomatics, Psychiatric Hospital, University of Zurich, Zurich, 8032, Switzerland
2Zurich, Neuroscience Center Zurich, University of Zurich, Zurich, 8057, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Aslan Abivardi
Saurabh Khemka
1Computational Psychiatry Research, Department of Psychiatry Psychotherapy and Psychosomatics, Psychiatric Hospital, University of Zurich, Zurich, 8032, Switzerland
2Zurich, Neuroscience Center Zurich, University of Zurich, Zurich, 8057, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dominik R. Bach
1Computational Psychiatry Research, Department of Psychiatry Psychotherapy and Psychosomatics, Psychiatric Hospital, University of Zurich, Zurich, 8032, Switzerland
2Zurich, Neuroscience Center Zurich, University of Zurich, Zurich, 8057, Switzerland
3Wellcome Centre for Human Neuroimaging and Max Planck UCL Centre for Computational Psychiatry and Ageing Research, University College London, London, WC1N 3BG, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Dominik R. Bach
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Decisions under threat are crucial to survival and require integration of distinct situational features, such as threat probability and magnitude. Recent evidence from human lesion and neuroimaging studies implicated anterior hippocampus (aHC) and amygdala in approach–avoidance decisions under threat, and linked their integrity to cautious behavior. Here we sought to elucidate how threat dimensions and behavior are represented in these structures. Twenty human participants (11 female) completed an approach–avoidance conflict task during high-resolution fMRI. Participants could gather tokens under threat of capture by a virtual predator, which would lead to token loss. Threat probability (predator wake-up rate) and magnitude (amount of token loss) varied on each trial. To disentangle effects of threat features, and ensuing behavior, we performed a multifold parametric analysis. We found that high threat probability and magnitude related to BOLD signal in left aHC/entorhinal cortex. However, BOLD signal in this region was better explained by avoidance behavior than by these threat features. A priori ROI analysis confirmed the relation of aHC BOLD response with avoidance. Exploratory subfield analysis revealed that this relation was specific to anterior CA2/3 but not CA1. Left lateral amygdala responded to low and high, but not intermediate, threat probability. Our results suggest that aHC BOLD signal is better explained by avoidance behavior than by threat features in approach–avoidance conflict. Rather than representing threat features in a monotonic manner, it appears that aHC may compute approach–avoidance decisions based on integration of situational threat features represented in other neural structures.

SIGNIFICANCE STATEMENT An effective threat anticipation system is crucial to survival across species. Natural threats, however, are diverse and have distinct features. To be able to adapt to different modes of danger, the brain needs to recognize these features, integrate them, and use them to modify behavior. Our results disclose the human anterior hippocampus as a likely arbiter of approach–avoidance decisions harnessing compound environmental information while partially replicating previous findings and blending into recent efforts to illuminate the neural basis of approach–avoidance conflict in humans.

  • anterior hippocampus
  • decision-making under predation
  • high-resolution fMRI
  • hippocampal subfields
  • lateral amygdala
  • operant conflict test

Introduction

Integrating divergent situational demands is critical to survival, in particular, when predatory or metabolic threat is involved (Korn and Bach, 2015, 2018, 2019). A standard laboratory model of this situation is provided by approach–avoidance conflict (AAC) tests, for example, open-field test and elevated plus-maze (Calhoon and Tye, 2015), which are thought to reflect aspects of human clinical anxiety disorders (Aupperle and Paulus, 2010). Situational threat features are manifold and distinct in these tests, and even more so in biological scenarios (Evans et al., 2019). For a human during wintertime, there is a low probability of being attacked when encountering a hibernating bear and a higher probability when coming across wolves, who are short on food. The metabolic loss incurred by a bear chase, however, may be much higher than when being charged by a single wolf. How the neural system represents and integrates such different threat dimensions, and how they influence behavior (e.g., the decision to approach food under threat or passively avoid threat), remains unknown.

In rodent AAC tests, cautious (“anxiety-like”) behavior is consistently reduced by anxiolytic drugs, such as benzodiazepines (Gray and McNaughton, 2000). Ventral hippocampus lesions have a similar impact (Kjelstrup et al., 2002; Bannerman et al., 2003; McHugh et al., 2004; Ito and Lee, 2016; Kirlic et al., 2017), and it has been suggested that behavioral control requires interplay of hippocampal subfields (Schumacher et al., 2018). Theta oscillations of hippocampal local field potential (Gordon et al., 2005), and synchronization with prefrontal cortex (PFC) (Adhikari et al., 2010; Padilla-Coreano et al., 2016), are increased in AAC, while area-specific circuits influence decisions (Wallis et al., 2019). In a human computer game resembling open-field test, benzodiazepines (Korn et al., 2017) and other anxiolytics (Bach et al., 2018) reduced cautious behavior similar to hippocampus (Bach et al., 2014) and amygdala (Korn et al., 2017) lesions in humans and nonhuman primates (Chudasama et al., 2008; Machado et al., 2009). Amygdala contribution is inconsistently reported in rodents (Kirlic et al., 2017); in humans, it appears to be specifically required for retreat from threat after reward collection, rather than for the decision to approach (Bach et al., 2019).

While this suggests involvement of hippocampus and amygdala in behavioral control, it remains elusive how different threat features, ultimately determining behavior, are represented and integrated. Features, such as magnitude and probability of threat, are not experimentally controlled in many tests that build on innate anxiety or that are extended in time. For example, we have shown using fMRI that neural mass activity of anterior hippocampus (aHC) increases with threat probability in continuous-time AAC (Bach et al., 2014). However, fMRI studies with more abstract AAC tests not requiringimmediate behavior have yielded conflicting results, somesuggesting the same relation of aHC activity with threat probability (Korn and Bach, 2019); others, a relation of aHC activity (Loh et al., 2017) or multivoxel patterns (O'Neil et al., 2015) with behavior.

Operant conflict tests provide the opportunity to moreprecisely control threat features as demonstrated in rodents (Evenden et al., 2009; Oberrauch et al., 2019) and humans (Bach, 2015, 2017; Bach et al., 2019). Here, we capitalized on this latter operant AAC test to disambiguate representation of attack probability, its metabolic cost, and behavior, in aHC and amygdala. We previously used the same task to show that putative hippocampal gamma oscillations, and hippocampal theta synchronization with PFC, increased with threat probability (Khemka et al., 2017). Presently, we gained from the superior spatial resolution of fMRI collecting 1.5 mm isotropic BOLD images focused on amygdala and hippocampus while participants played the game. On each trial, they could either collect, or forego, a monetary token under threat of capture by a predator. Threat probability was defined by the predator wake-up rate and learned by experience; threat magnitude by potential token loss and explicitly signaled.

Materials and Methods

Participants

Twenty participants were recruited from general and student population in Zurich (mean age ± SD, 23.10 ± 3.34 years; 11 female). Participants had no prior history of neurologic or psychiatric disease and reported normal or corrected-to-normal vision. One participant was excluded from fMRI analysis because of a technical fault in MRI recordings, but included in behavioral analysis. Behavioral results remained consistent after removal of this participant. All participants gave their written informed consent before participation. The study protocol was in full accordance with the Declaration of Helsinki and approved by the governmental ethics committee (Kantonale Ethikkommission Zürich).

Experimental procedure

Participants performed an AAC computer game as previously used by Khemka et al. (2017), which was modified from Bach (2015). At the beginning of each trial, the human player was located in a “safe place” in the bottom block of a 2 × 2 diamond grid (Fig. 1) opposite of a sleeping predator, and was given the opportunity to collect a monetary token that would appear in the left or right grid block. Red diamonds underneath the grid explicitly signaled the number of tokens that would be lost (0-5) if captured by the predator. Threat probability was implicitly signaled through frame color (blue, pink, and orange). Threat probability was implemented by setting the wakeup rate per time unit to result in catch probabilities of 0.1, 0.2, or 0.3 per 100 ms spent outside of the safe grid block. These probabilities were learned from experience during 36 preceding training trials without token loss that did not count toward ultimate earnings.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

A, FOV focused on amygdala/hippocampus. Image represents the EPI coverage across participants (thresholded at p = 0.5), overlaid on a mean T1 image in MNI space. B, AAC task. In each trial, the human participant (green triangle) started out in a safe (dark gray) grid block opposite a sleeping predator (gray circle) and was presented with a reward token (yellow rhombus) on the left or right side. Threat probability was signaled by frame color (blue/magenta/orange). The player then had the choice to collect the token using left/right keys to move out of, and return to, the safe place. If caught while outside, the amount of tokens signaled in red below the frame (here two) was lost, thus constituting the magnitude of threat.

A token appeared after a random time interval drawn from a truncated γ distribution (k = 2, θ = 1; mean = 2 s, t ≤ 6 s). If the player chose not to collect the token, it would disappear after another time interval drawn from the same distribution, and the trial would end 1 s later. If the player went to acquire the token and successfully returned to the safe place, the trial would proceed until the same predetermined end time. Finally, if the predator caught the player, the predator changed its color from gray to red and remained on the screen until the predetermined end time of the trial. After a random intertrial interval also drawn from a γ distribution truncated at t ≤ 4 s, during which a blank screen was presented, the next trial would start. Participants completed 648 trials in random order, balanced for each combination of experimental factors (i.e., threat magnitude and threat probability). Participants were instructed beforehand that their payment depended on performance in six trials randomly drawn from the experiment excluding training trials. The experiment was programmed in Cogent (version 2000v1.25; www.vislab.ucl.ac.uk/Cogent) and MATLAB (version 7.14; The MathWorks).

Acquisition of MRI data

Data were recorded in a 3.0 Tesla MRI scanner (Phillips Achieva, Phillips Medical Systems) using a 32-channel head coil. Anatomical images were acquired using a 0.76 mm isotropic resolution T1-weighted scan (TR = 7.37 ms, TE = 3.29 ms, flip angle = 8°, FOV = 255 × 255 × 180 mm, matrix = 336 × 336, thickness = 0.76 mm, in-plane resolution = 0.76 × 0.76 mm2, slice tilt = 0°, 237 sagittal slices) and a 1.0 × 0.5 × 0.5 mm resolution T2-weighted scan centered on hippocampus (TR = 3200 ms, TE = 353 ms, flip angle = 90°, FOV = 200 × 52 × 200 mm, matrix = 400 × 400, thickness = 1 mm, in-plane resolution = 0.5 × 0.5 mm2, slice tilt = 22°, 104 transverse slices). B0 field maps were acquired with a double-echo fast gradient echo sequence (TR = 698.22 ms, TE = 4.10 and 7.10, flip angle = 44°, FOV = 240 × 224 × 240 mm, matrix = 80 × 80, thickness = 3 mm, in-plane resolution = 3 × 3 mm2, slice tilt = 0°, 2 × 64 sagittal slices). Functional images during the approach–avoidance paradigm were recorded with 1.5 mm isotropic resolution T2*-weighted EPI sequence (TR = 2800 ms, TE = 30 ms, flip angle = 85°, in-plane resolution = 1.5 × 1.5 mm2, FOV = 216 × 54 × 216 mm, matrix = 144 × 144; 36 transverse slices with thickness = 1.5 mm; slice order = interleaved ascending; slice tilt = −40°). FOV was centered on amygdala/hippocampus, but also encompassed striatum, thalamus, prefrontal cortices with exclusion of orbitofrontal cortex and cranio-posterior segments of frontal lobe, greater parts of temporal lobes and cerebellum, as well as complete coverage of insular cortices and brainstem (Fig. 1).

Preprocessing of MRI data

Preprocessing of functional images was performed using a standard pipeline in SPM12 (Statistical Parametric Mapping; Wellcome Center for Human Neuroimaging, London; http://www.fil.ion.ucl.ac.uk/spm/software/spm12). In a first step, slice time correction was performed to account for differences in acquisition time of individual brain slices (Sladky et al., 2011). Geometric distortions because of susceptibility-induced field inhomogeneities were addressed using a combined approach, which takes static distortions as well as changes in distortion because of head motion into account (Andersson et al., 2001; Hutton et al., 2002). Static distortions were derived for each subject individually from a B0 field map using the FieldMap toolbox in SPM12. EPIs were subsequently realigned and unwarped integrating the measured static distortion and the estimation of distortion caused by head motion, as well as head motion itself. EPI images as well as T2w images were then coregistered to the individual T1w whole-brain image using a 12-parameter affine transformation. Finally, EPI images were normalized into MNI space and smoothed using an isotropic 8 mm FWHM Gaussian kernel for primary mass-univariate analysis, and a 4 mm FWHM Gaussian kernel for a secondary analysis to improve localization of effects in amygdala and hippocampus. We note that the smoothing kernel must strike a balance between anatomic intersubject variability and regional specificity (Mikl et al., 2008). Thus, the larger smoothing kernel is expected to be more sensitive in detecting activations, but the smaller kernel can provide additional information on the localization of clusters. Unsmoothed EPI images in native space were used for ROI analysis.

fMRI analysis (focused brain coverage)

In a primary analysis (P1), we defined a GLM consisting of a δ function at token appearance (consistent with Khemka et al., 2017), convolved with a canonical HRF. Parametric modulators for linear and quadratic effect of threat probability (1-3), linear and quadratic effect of threat magnitude (0-5), and linear interaction effect of threat probability × magnitude were also convolved with the HRF. All parametric modulators were serially orthogonalized. Motion correction parameters were included as six additional regressors of no interest.

To distinguish effects of behavior from threat features, we ran a second parametric analysis (P2) with approach or avoidance as a first parametric modulator, followed by linear and quadratic effect of threat probability, linear and quadratic effect of threat magnitude, linear combination of threat probability × magnitude, linear combinations of approach × probability, approach × magnitude, and finally approach × probability × magnitude.

To extricate effects of threat probability and magnitude that were specific to ensuing behavior, we computed a third GLM with separate trial regressors for approach trials and avoidance trials (P3), each with parametric modulators for linear and quadratic effects of threat probability and magnitude as well as for effect of linear combination of probability × magnitude. To assess a potential relation of neural activity with response latencies, we defined three further models in an analogous manner with parametric regressors for linear and quadratic effects of approach and withdrawal latency during approach trials. We controlled for threat features in all models using serial orthogonalization. Since the player was often caught by the predator during attempts to obtain a token, data for withdrawal latency were available on fewer trials than for approach latency. Thus, we defined one model for approach latency over all trials without control for withdrawal latency, and two models over trials without capture, where approach and withdrawal latencies were orthogonalized in respect to each other.

ROI definition

Subcortical and cortical structures, including hippocampus and amygdala, were identified in native subject space using the “recon-all” pipeline in FreeSurfer version 6.0 (http://surfer.nmr.mgh.harvard.edu/) (Dale et al., 1999; Fischl et al., 1999a,b, 2002, 2008; Segonne et al., 2005; Desikan et al., 2006). Individual voxels were assigned neuroanatomical labels in an automated volumetric subcortical parcellation based on a probabilistic atlas from a manual training set (Fischl et al., 2002). The hippocampus segmentation was then further parcellated into anterior and mid-to-posterior hippocampus by automatically splitting the mask at one-third length along the anterior-posterior axis of the image in MATLAB (Strange et al., 2014). For exploratory purposes, CA1 and CA2/3 subfields as well as a mask for dentate gyrus of the were obtained from the higher resolution T2w images with FreeSurfer 6.0, which uses a statistical atlas based on ultra-high-resolution ex vivo data images for segmentation. (Iglesias et al., 2015). CA1, CA2/3, and dentate gyrus images were then multiplied with the binary aHC mask to focus only on the anterior segments.

For small-volume correction (SVC) of group-level analysis, a group-level bilateral hippocampus mask was generated by warping the individual bilateral hippocampus masks into MNI space using the deformation fields acquired during normalization of whole-brain T1w images in SPM12. These were then averaged, thresholded at 0.1, and binarized using the SPM12 function ImCalc. For visualization, group-level masks in MNI space for all significant clusters were extracted using SPM12 Results.

ROI fMRI analysis

For analysis of estimated condition × condition BOLD response, averaged within ROI, we defined a first-level GLM with separate regressors for 36 possible distinct combinations of threat probability (1-3), magnitude (0-5), and behavioral response (0/1). We extracted estimated condition × condition BOLD response for aHC, anterior subfields CA1 and CA2/3, anterior dentate gyrus, entire amygdala, centrocortical and basolateral amygdala subnucleus groups, and, for visualization, for significant clusters from focused brain analysis.

Statistical analysis

Image-based statistical tests for fMRI analysis were performed with SPM group-level analysis using cluster-level familywise error (FWE) correction for multiple comparisons at a voxel-inclusion threshold of p < 0.001 (correction for whole FOV, or small volume corrected for hippocampus) and applying a random-field theory based approach as implemented in SPM (Worsley et al., 1992).

For a priori ROIs amygdala and hippocampus, we implemented a mixed-effects analysis in R 3.4.3 (www.r-project.org) using function lmer (lme4 package) with the following fixed effects that followed the definition of the voxelwise analysis while adding a hemispheric difference: linear and quadratic effects of threat probability and magnitude, behavioral response and hemisphere, and ensuing interactions. We added a random intercept for subject. This resulted in the R formula (where all predictors are numerical rather than factors): Embedded Image

Exploratory analysis in aHC subfields (CA1, CA2/3) and amygdala subnuclei groups (basolateral and centrocortical) was then performed using the same formula. Significance level α was adjusted for multiple comparisons across two ROIs for a priori tests, and four ROI for exploratory analysis using the Holm-Bonferroni method (Holm, 1979). To further differentiate for region- and subfield-specific effects in an exploratory analysis, ROI was included as a fixed effect in one combined model for amygdala versus aHC and another for anterior CA1 versus CA2/3. Last, post hoc ROI analysis was performed in anterior dentate gyrus using the initial model without ROI as factor. Statistical analysis of behavioral data were likewise performed in R using a linear mixed-effects model (lme4 package), which can deal with the inherently unbalanced data (for details, see Bach, 2015; Khemka et al., 2017), using Satterthwaite approximation to degrees of freedom to appropriately control the false positive rate (Luke, 2017).

Data availability

A repository of unthresholded SPM activation maps for parametric analyses P1-P3 (group level; 4 and 8 mm kernel)is publicly available at https://github.com/a-abivardi/neural-threat-behavior-AAC-fMRI (Abivardi et al., 2020).

Results

Behavioral results

We first interrogated whether behavior was comparable to previous findings. Passive avoidance (i.e., the proportion of avoidance over approach decisions) increased with higher threat probability and magnitude. Behavioral inhibition, measured as approach latency, increased with higher threat probability and magnitude, whereas the opposite pattern was observed for withdrawal latency (Fig. 2; Table 1). These results replicate previous reports (Bach, 2015, 2017; Khemka et al., 2017; Bach et al., 2019) and are in concordance with known behavior from rodent studies.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

A, Proportion of approach–avoidance decisions ± SEM defined as SD of generalized linear mixed-effects model residuals divided by square root of number of data points. B, Approach and withdrawal latency, estimated from linear mixed-effects model ± SEM (defined as SD of model residuals divided by square root of number of data points).

View this table:
  • View inline
  • View popup
Table 1.

Linear and omnibus effects of threat features on behavioral responses

Mass-univariate fMRI results

As threat features and approach–avoidance behavior are strongly related, we chose a threefold parametric design (P1–P3) to disentangle distinct effects using serial orthogonalization as implemented in SPM12. In a primary analysis (P1; Table 2), we analyzed how BOLD signal related to linear and quadratic components of the two threat dimensions and their interactions, by including them as parametric modulators. A second analysis (P2; Table 3) prepended these modulators by behavioral response (approach–avoidance), making use of serial orthogonalization in SPM12, and further examined interactions between threat dimensions and behavior. Last (P3; Table 4), approach and avoidance trials were analyzed separately to account for behavior-specific effects of threat dimensions on brain activation. All results were corrected for FWE within the FOV. For bilateral hippocampus, additional FWE SVC was performed using a group-level bihemispheric mask, as we had strong a priori hypotheses for this region. Mass-univariate results are reported for images smoothed with an 8 mm Gaussian kernel, unless otherwise specified. Cluster anatomy was visually compared to the schematic brain atlas by Mai et al. (2016). Peak activation coordinates were labeled using the Automated Anatomical Labeling (AAL) atlas (Tzourio-Mazoyer et al., 2002) (Tables 2–4).

View this table:
  • View inline
  • View popup
Table 2.

Parametric analysis (P1): effects of threat features on brain activation

View this table:
  • View inline
  • View popup
Table 3.

Parametric analysis (P2): effects of approach–avoidance behavior and serially orthogonalized threat features on brain activation

View this table:
  • View inline
  • View popup
Table 4.

Parametric analysis (P3): effects of threat features on brain activation, separately for approach and for avoidance trials

In analysis P1, we observed higher BOLD signal with a combination of higher threat probability and higher threat magnitude in left aHC (specifically subiculum) and entorhinal cortex (linear × linear interaction, FOV-corrected corrected; Fig. 3; Table 2). This effect was not reproduced in a secondary analysis with a narrower smoothing kernel size of 4 mm. There were no mass-univariate effects in the amygdala. Exploratory analysis of the remaining brain coverage (Table 2) revealed higher BOLD signal with lower threat probability (linear negative effect of threat probability) in left dorsolateral PFC (dlPFC), a cluster extending into left putamen and anterior insula, and in the posterior lobe of the right cerebellum. Low threat magnitude was related to higher BOLD signal (linear negative effect of threat magnitude) in left internal capsule/putamen, posterior short gyrus of left insula, ventrolateral PFC (vlPFC), left inferior temporal gyrus, and multiple clusters in bilateral cerebellum and vermis.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Cluster-level significant aHC and amygdala clusters from parametric analysis with, for purposes of illustration, extracted estimated condition × condition BOLD response ± SEM as defined by SD of BOLD response amplitude estimates divided by square root of number of data points. Red represents primary analysis clusters using 8 mm FWHM smoothing kernel. Blue represents secondary analysis cluster (B) using 4 mm kernel. A, Left anterior subiculum-entorhinal cortex cluster modulated by combined threat probability and magnitude (linear-positive interaction effect, analysis P1). B, Left anterior subiculum-entorhinal cortex area relating to avoidance (P2; small-volume corrected). Secondary analysis localized this cluster to the left anterior CA3/dentate gyrus area. BOLD estimates ± SEM are displayed for the 4 mm cluster. C, Left lateral amygdala cluster quadratically modulated by threat probability (P1). All results are FEW-corrected at cluster level (p < 0.05; voxel inclusion threshold: p < 0.001).

In P2, there were no significant hippocampus or amygdala clusters at FOV correction. After SVC in bilateral hippocampus, we observed a cluster in which avoidance behavior related to higher BOLD activity. This cluster in left aHC (subiculum) and entorhinal cortex was located slightly posterior to the area modulated by combined threat magnitude and probability in P1 (Fig. 3; Table 3). This result was replicated in a secondary analysis using 4 mm kernel smoothed images for higher localization accuracy. In this analysis, the cluster was in adjacent location, but more superior in the anterior CA3/dentate gyrus area (Fig. 3). In a distinct cluster in left middle hippocampus (8 mm kernel only), high BOLD signal related to high threat magnitude. Exploratory analysis of the remaining brain coverage (Table 3) revealed that approach behavior related to BOLD signal in two large clusters encompassing bilateral cerebellum and extending from bilateral thalamus to striatum and midbrain structures. Furthermore, approach behavior related to activation in left substantia nigra, bilateral anterior cingulate cortex (ACC) and dorsomedial PFC, anterior short gyrus of right insula, opercular part of right inferior frontal gyrus, and precentral cortex. These clusters showed partial overlap with impact of low threat magnitude as shown in P1, in bilateral cerebellum, left putamen, and anterior insula, as well as with impact of low threat probability in right cerebellum and left putamen. After controlling for behavior in P2, no linear effects of threat probability were observed. A quadratic modulating effect of threat probability emerged in the left lateral amygdaloid nucleus (8 mm kernel only; i.e., high activation for low and for high, but not for medium threat probability) (Fig. 3). High threat magnitude was related to high BOLD signal in right anterior insula (anterior and middle short gyrus) and frontal operculum. A second adjacent cluster in the right frontal operculum showed a linear relation of BOLD signal with threat magnitude specifically in combination with approach (interaction threat magnitude × behavior). This effect, however, was only estimable in 16 subjects.

In P3, there were no hippocampus or amygdala clusters at whole-brain or SVC. Exploratory analysis of the remaining brain coverage (Table 4) revealed that, for approach trials, high threat magnitude was associated with high BOLD signal in bilateral anterior short gyrus of insula, opercular part of inferior frontal gyrus, and bilateral anterior cingulate. Overlap with activation related to approach behavior in P2 was seen primarily in anterior cingulate, whereas overlap with activation related to high threat magnitude in P2 was seen in right anterior insula (replicating the previous finding). Linear interaction of high threat probability and magnitude in approach trials (estimable in 17 subjects) furthermore related to BOLD signal in right superior colliculus (partial overlap with approach-related activation from P2) and a cluster extending from left brachium of inferior colliculus into the medial geniculate nucleus. Specifically, BOLD response increased with threat magnitude for medium and high threat probabilities, but not for low probability.

Effects in avoidance trials were only partially estimable because of unequal distribution (i.e., relative scarcity of avoidance trials across participants) and yielded no significant results.

Finally, we found a positive relation of approach latency with BOLD activation in left ACC and right anterior insula over all approach trials. When controlling for withdrawal latency in the subset of trials where the player was not caught, neither this nor any other relation was seen. There were no significant clusters in relation to withdrawal latency independent of approach latency or threat features.

ROI analysis results

A priori ROI analysis was conducted across both anterior hippocampi, and across both amygdalae. Results were correctedfor multiple comparisons across the two ROIs using Holm-Bonferroni adjusted significance level (Table 5).

View this table:
  • View inline
  • View popup
Table 5.

ROI analyses in aHC and amygdala

In the aHC ROI, we observed a linear main effect of behavioral response and quadratic main effects for threat probability and magnitude. BOLD signal was higher for avoidance than for approach trials. Similar to left lateral amygdala in parametric analysis P2, aHC also responded to low and high threat probability (positive quadratic effect). Strikingly, this effect seemed to be behavior-dependent and lateralized as left aHC responded to high threat probability and right hippocampus activation related to low threat probability, both during avoidance only (quadratic × linear interaction of threat probability and behavior, and linear interaction of threat probability × behavior × hemisphere) (Fig. 4). Moreover, for zero threat magnitude, hippocampus BOLD signal was low, while increasing to peak levels for low to intermediate levels and falling again with higher magnitude, resulting in a significant negative quadratic pattern. Finally, aHC exhibited a complex linear interaction of threat features and hemisphere: BOLD response showed a negative linear relation with threat magnitude for high threat probability and for left hemisphere only.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

ROI analyses for aHC (A), amygdala versus aHC (B), and anterior CA1 versus anterior CA2/3 (B). A, Interaction effect of threat probability, approach, and hemisphere: that is, estimated condition × condition BOLD response amplitudes ± SEM defined as SD of mixed-effects model residuals divided by square root of number of data points. B, Interaction effect of behavior × ROI for amygdala versus aHC and anterior CA1 versus anterior CA2/3 (condition × condition BOLD response ± SEM).

The response to low and high threat probability seen in lateral amygdala after control for behavior was replicated in the amygdala ROI analysis (positive quadratic main effect), while interactions between threat probability and behavior were not detected. Left hemisphere showed overall higher BOLD responses in the amygdala.

A combined analysis of amygdala and aHC revealed distinct activation patterns in relation to behavior (Fig. 4; Table 5). While aHC was clearly more active during avoidant behavior, amygdala exhibited a slightly higher BOLD response during approach (behavior × ROI interaction). The quadratic response to threat magnitude appeared to be specific to aHC; moreover, aHC was different from amygdala in its lateralized response to threat probability during avoidance (linear threat probability × magnitude × behavior × ROI interaction).

In addition to planned ROI analysis, exploratory follow-up analyses were conducted in bilateral aHC subfields CA1 and combined CA2/3 as well as basolateral and centrocortical amygdala ROIs. Results were corrected for multiple comparisons across four ROIs using Holm-Bonferroni method (Table 5).

Subfield analysis in anterior CA1 revealed complex and interacting effects of threat dimensions with distinct activation patterns for approach and avoidant behavior and depending on hemisphere. As in entire aHC, a relation to low and high threat probabilities was seen during avoidance only (quadratic × linear interaction). A complex linear interaction effect of threat features and hemisphere also similar to entire aHC was observed.

Activation in combined hippocampal subfield CA2/3 was higher for avoidance than approach behavior, reflecting the main effect found in the combined aHC ROI. The difference (or increase) in BOLD response for avoidance compared with approach conditions was furthermore higher for anterior CA2/3 (mean = 1.31, SD = 1.28) than for subfield CA1 (mean = 0.22, SD = 1.17) in a post hoc paired sample t test (t(17) = −3.31, p = 0.004), underlining the difference between the two subfields. A combined model for the subfields confirmed this distinction with a significant behavior × ROI interaction effect (Fig. 4; Table 5). At the suggestion of a reviewer, we analyzed BOLD responses in anterior dentate gyrus, based on findings that this area may have a role similar to that of CA3. However, we did not find a significant relation with avoidant behavior here.

Further exploratory analysis in amygdala subnuclei using a probabilistic amygdala mask from a previous study (Abivardi and Bach, 2017) revealed activation of basolateral amygdala with increasing threat magnitude (linear main effect). This effect was not seen for entire amygdala. Left basolateral and centrocortical amygdala were more active than amygdala of the right hemisphere as also seen for entire amygdala. Also, centrocortical amygdala exhibited heightened BOLD response to intermediate threat magnitudes, especially during avoidance (quadratic main effect + quadratic × linear interaction). We note that the centrocortical amygdala parcellation was defined by structural connectivity with lateral orbitofrontal cortex (Bach et al., 2011) based on preferred projections to central, medial, and cortical amygdala in rodents and nonhuman primates (Carmichael and Price, 1995; McDonald, 1998; Pitkänen, 2000). Morphologically, this group parcellation (resulting from a sample with size of n = 50) (Abivardi and Bach, 2017) probably includes central, medial, cortical, as well as basomedial nuclei.

Discussion

Rodent and human ventral or aHC is crucial to cautious behavior in AAC tests (Ito and Lee, 2016). However, how distinct threat features are represented and integrated has only recently received attention (Korn and Bach, 2019). Harnessing a human operant AAC computer game during high-resolution fMRI, we investigated representation of threat probability and threat magnitude, and of approach or avoidance behavior, in aHC and amygdala. Two key findings emerged. First, aHC BOLD activity was related to behavioral avoidance, particularly for CA2/3 but not for CA1. Second, there was no evidence that aHC unambiguously represents elementary threat features in a linear manner. Similarly, exploratory analyses of further brain areas within our limited coverage did not reveal a coherent linear representation of threat probability or magnitude.

In mass-univariate analysis, we observed that BOLD signal in left aHC/entorhinal cortex, specifically the subiculum-entorhinal area, was related to the combination of high probability and magnitude of threat (analysis P1), both of which result in more avoidant behavior. After controlling for behavior (P2), no such relation was found. Instead, neural activity in a slightly more posterior cluster was related to avoidant behavior. Using a smaller smoothing kernel to fully harness high spatial resolution, we localized this second cluster to the anterior CA3/dentate gyrus area. A priori ROI analysis confirmed these findings:averaged aHC BOLD signal was increased during avoidance. Follow-up analysis of anterior subfields revealed that this avoidance-related increase occurred in CA2/3 but not CA1. This finding resonates with a rat experiment by Schumacher et al. (2018) who demonstrated that selective pharmacological inactivation of ventral CA3 increased approach behavior. A role paralleling CA3 has been recently described for rodent ventral dentate gyrus (Yeates et al., 2020). We note that it remains possible that our CA2/3 parcellations contain individual voxels belonging to bordering dentate gyrus. Nonetheless, exploratory ROI analysis in dentate gyrus did not detect a similar effect here. On the other hand, CA1 activity in our study showed no simple relationship with threat features or behavior, whereas selective pharmacological ventral CA1 inactivation increased avoidance in a previous rat experiment (Schumacher et al., 2018).

Our finding of aHC activity relating to avoidance is in keeping with a previous human fMRI study involving abstract AAC decisions, which reported inferior aHC BOLD activity during avoidance (Loh et al., 2017), in proximity to the left aHC cluster relating to avoidance here. We note that, in this previous study, most voxels in this cluster were labeled as belonging to CA1; however, the authors noted that anatomic specificity might have been limited because of lower spatial resolution (3 mm), as opposed to the present approach.

In a lesion study with the same paradigm as used here, we found that hippocampus lesions impaired approach–avoidance decisions, whereas impact of threat on other behaviors remained intact (Bach et al., 2019), further suggesting a specific role of aHC in generating avoidance behavior. Selective amygdala and hippocampus lesions were moreover associated with shorter approach latency, but not with a different relationship between threat and approach latency. This may suggest that these regions do not contribute to parametric variation in approach latency. In keeping with this, we presently found that variation in approach latency did not relate to hippocampus or amygdala signal. A previous magnetoencephalography study reported a relation between approach latency and posterior hippocampus activity (Khemka et al., 2017), not observed here.

Regarding threat feature representation, ROI analysis revealed a more complicated picture than previously assumed. Although we observed significant responses of aHC to low and intermediate threat magnitude levels, forming a quadratic pattern, BOLD signal also depended on interactions between threat features and behavior, with some effects strikingly different between hemispheres. Specifically, left aHC responded to high threat probability, whereas right aHC related to low probability during avoidance. In humans, left hippocampus has been implicatedin contextual and spatial memory encoding, whereas right hippocampus has been linked to navigation accuracy (Maguire et al., 1998; Spiers et al., 2001). Hemisphere-specific connectivity profiles in human aHC (Robinson et al., 2016) and task-related activity in rat ventral HC (Sakaguchi and Sakurai, 2017) have been reported. However, we note the historical and ongoing debate on lateralization of emotional functions, which is based on partly contradicting observations (Gainotti, 2019). It would therefore appear useful to replicate our findings in an independent sample.

In contrast, a previous fMRI study (Bach et al., 2014) using a more ethological paradigm reported linearly increasing activity in left aHC with higher threat probability. Accounting for the influence of behavior in this temporally extended paradigm was, however, difficult. Furthermore, previous threat probabilities were higher (0.2/0.5/0.8) than the current ones (0.1/0.2/0.3). Also, this previous study did not explicitly control threat magnitude, which we achieved here. Another fMRI study involving more abstract foraging decisions under predation (Korn and Bach, 2019) found a cluster in which aHC signal increased with threat probability (0.1–0.4) but a partly overlapping cluster in which aHC signal decreased from 0.1 to 0.3 and increased from 0.3 to 0.4, yielding an overall quadratic pattern. To reconcile these findings, it appears necessary to cover a larger probability range.

As a further finding, BOLD signal in left lateral amygdala related to low and high, but not intermediate, threat probability independent of behavior (P2). ROI analysis in entire amygdala replicated this behavior-independent activation pattern. The role of amygdala in AAC is reported more controversially than for aHC (Kirlic et al., 2017); nevertheless, a recent human lesion study suggested specific involvement in controlling vigor of return to safety (Bach et al., 2019).

Results from exploratory focused brain analysis revealed several clusters with complex and differential relation with threat features and behavior. Left dlPFC response was related to low threat probability (P1), resonating with reports that anxiety is inversely correlated with dlPFC activity (Balderston et al., 2017). Right vlPFC has been implicated in motor inhibition and characterized as a “brake,” which however has been debated (Aron et al., 2004, 2014; Swick and Chatham, 2014). Here, right vlPFC activity related to low threat magnitude (P1), approach behavior (P2), and high magnitude during approach (P3). While we did observe behavioral inhibition during approach trials relating to threat, the relation to approach behavior seems at odds with pure motor inhibition. Swick and Chatham (2014) propose that vlPFC monitors action-relevant situational changes, compatible with response to threat magnitude here.

In a recent optogenetic study, ACC activation decreased rodent freezing behavior via input to basolateral amygdala (Jhang et al., 2018). Anterior cingulate also appears to signal value predictions of rewards and punishments (Monosov, 2017). Conceptually, dorsal anterior cingulate has been theorized to monitor conflict (or expected value of top-down control) (Botvinick et al., 1999; Shenhav et al., 2016) or to adaptively track context-relevant and action-guiding variables (Heilbronner and Hayden, 2016). Here, dorsal anterior cingulate related to approach behavior (P2) while also relating to rises in threat magnitude during approach trials (P3). The former finding matches anterior cingulate role in freezing in mice and supports a more active role arbitrating behavior. The latter finding may equally well constitute measurement of conflict, context-relevant variable tracking or punishment-related value predictions.

Anterior insular cortex activity was related to approach decisions and both threat features. Left insula related to low threat magnitude and probability before accounting for behavior (P1), right anterior insula activation was related to approach (P2) and bilateral insula to high threat magnitude in separated approach trials (P3). This contrasts reports from a study implicating anterior insula activation in avoidance decisions (Aupperle et al., 2015). Overall, insula showed similar responses to anterior cingulate, adding to evidence of their close functional link (Medford and Critchley, 2010).

Limitations of our study include the use of a limited FOV as a necessary compromise for higher-resolution imaging of ROIs. Furthermore, a relative scarcity of avoidance decisions across participants, compared with previous studies with cumulative token collection (Bach, 2015, 2017; Bach et al., 2019) hindered analysis of threat representation during avoidance and reduced power to detect brain areas involved in avoidant decision-making. A focus on single-stage decisions in the present study precludes analyzing to what extent assumptions about future foraging attempts may prompt avoidance on the current one (Korn and Bach, 2019; Zorowitz et al., 2020). Last, orthogonalization in SPM12 penalizes parametric modulators in a serial manner along the design matrix, which demands careful interpretation of results (Mumford et al., 2015).

In conclusion, in this study, we disambiguated a relation of neural tissue activity with behavior and situational threat features. AHC BOLD signal, in particular in CA2/3, increased when participants avoided threat. Representation of threat features showed a complicated pattern, and for threat probability depended on behavior. This is in line with a notion that hippocampus does not linearly represent threat features but retrieves them, possibly in a manner that changes over time, to compute decisions. It would be useful to increase the range of these threat features, as well as improve both spatial and temporal precision of recording, for example using electrophysiology, to understand how these computations emerge over time in different hippocampal subfields.

Footnotes

  • A.A. was supported by Wilhelm-Hurka Foundation 3389322 and EMDO Foundation 1045. D.R.B. was supported by European Research Council, European Union's Horizon 2020 Research and Innovation Program Grant Agreement ERC-2018 CoG-816564 ActionContraThreat. The Wellcome Trust Center for Human Neuroimaging was supported by core funding from the Wellcome Trust (203147/Z/16/Z).

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Aslan Abivardi at aslan.abivardi{at}puk.zh.ch

This is an open-access article distributed under the terms of the Creative Commons Attribution License Creative Commons Attribution 4.0 International, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    1. Abivardi A,
    2. Bach DR
    (2017) Deconstructing white matter connectivity of human amygdala nuclei with thalamus and cortex subdivisions in vivo. Hum Brain Mapp 38:3927–3940. doi:10.1002/hbm.23639 pmid:28512761
    OpenUrlCrossRefPubMed
  2. ↵
    1. Abivardi A,
    2. Khemka S,
    3. Bach DR
    (2020) Hippocampal representation of threat features and behavior in a human approach–avoidance conflict anxiety task: unthresholded SPM activation maps [dataset]. https://doi.org/10.5281/zenodo.3893442.
  3. ↵
    1. Adhikari A,
    2. Topiwala MA,
    3. Gordon JA
    (2010) Synchronized activity between the ventral hippocampus and the medial prefrontal cortex during anxiety. Neuron 65:257–269. doi:10.1016/j.neuron.2009.12.002 pmid:20152131
    OpenUrlCrossRefPubMed
  4. ↵
    1. Andersson JL,
    2. Hutton C,
    3. Ashburner J,
    4. Turner R,
    5. Friston K
    (2001) Modeling geometric deformations in EPI time series. Neuroimage 13:903–919. doi:10.1006/nimg.2001.0746 pmid:11304086
    OpenUrlCrossRefPubMed
  5. ↵
    1. Aron AR,
    2. Robbins TW,
    3. Poldrack RA
    (2004) Inhibition and the right inferior frontal cortex. Trends Cogn Sci 8:170–177. doi:10.1016/j.tics.2004.02.010
    OpenUrlCrossRefPubMed
  6. ↵
    1. Aron AR,
    2. Robbins TW,
    3. Poldrack RA
    (2014) Inhibition and the right inferior frontal cortex: one decade on. Trends Cogn Sci 18:177–185. doi:10.1016/j.tics.2013.12.003 pmid:24440116
    OpenUrlCrossRefPubMed
  7. ↵
    1. Aupperle RL,
    2. Paulus MP
    (2010) Neural systems underlying approach and avoidance in anxiety disorders. Dialogues Clin Neurosci 12:517–531. pmid:21319496
    OpenUrlPubMed
  8. ↵
    1. Aupperle RL,
    2. Melrose AJ,
    3. Francisco A,
    4. Paulus MP,
    5. Stein MB
    (2015) Neural substrates of approach–avoidance conflict decision-making. Hum Brain Mapp 36:449–462. doi:10.1002/hbm.22639 pmid:25224633
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bach DR
    (2015) Anxiety-like behavioural inhibition is normative under environmental threat-reward correlations. PLoS Comput Biol 11:e1004646. doi:10.1371/journal.pcbi.1004646 pmid:26650585
    OpenUrlCrossRefPubMed
  10. ↵
    1. Bach DR
    (2017) The cognitive architecture of anxiety-like behavioral inhibition. J Exp Psychol Hum Percept Perform 43:18–29. doi:10.1037/xhp0000282 pmid:27797550
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bach DR,
    2. Behrens TE,
    3. Garrido L,
    4. Weiskopf N,
    5. Dolan RJ
    (2011) Deep and superficial amygdala nuclei projections revealed in vivo by probabilistic tractography. J Neurosci 31:618–623. doi:10.1523/JNEUROSCI.2744-10.2011 pmid:21228170
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Bach DR,
    2. Guitart-Masip M,
    3. Packard PA,
    4. Miro J,
    5. Falip M,
    6. Fuentemilla L,
    7. Dolan RJ
    (2014) Human hippocampus arbitrates approach–avoidance conflict. Curr Biol 24:541–547. doi:10.1016/j.cub.2014.01.046 pmid:24560572
    OpenUrlCrossRefPubMed
  13. ↵
    1. Bach DR,
    2. Korn CW,
    3. Vunder J,
    4. Bantel A
    (2018) Effect of valproate and pregabalin on human anxiety-like behaviour in a randomised controlled trial. Transl Psychiatry 8:157. doi:10.1038/s41398-018-0206-7 pmid:30115911
    OpenUrlCrossRefPubMed
  14. ↵
    1. Bach DR,
    2. Hoffmann M,
    3. Finke C,
    4. Hurlemann R,
    5. Ploner CJ
    (2019) Disentangling hippocampal and amygdala contribution to human anxiety-like behaviour. J Neurosci 39:8517–0419. doi:10.1523/JNEUROSCI.0412-19.2019
    OpenUrlCrossRef
  15. ↵
    1. Balderston NL,
    2. Liu J,
    3. Roberson-Nay R,
    4. Ernst M,
    5. Grillon C
    (2017) The relationship between dlPFC activity during unpredictable threat and CO2-induced panic symptoms. Transl Psychiatry 7:1266. doi:10.1038/s41398-017-0006-5 pmid:29213110
    OpenUrlCrossRefPubMed
  16. ↵
    1. Bannerman DM,
    2. Grubb M,
    3. Deacon RM,
    4. Yee BK,
    5. Feldon J,
    6. Rawlins JN
    (2003) Ventral hippocampal lesions affect anxiety but not spatial learning. Behav Brain Res 139:197–213. doi:10.1016/S0166-4328(02)00268-1
    OpenUrlCrossRefPubMed
  17. ↵
    1. Botvinick M,
    2. Nystrom LE,
    3. Fissell K,
    4. Carter CS,
    5. Cohen JD
    (1999) Conflict monitoring versus selection-for-action in anterior cingulate cortex. Nature 402:179–181. doi:10.1038/46035 pmid:10647008
    OpenUrlCrossRefPubMed
  18. ↵
    1. Calhoon GG,
    2. Tye KM
    (2015) Resolving the neural circuits of anxiety. Nat Neurosci 18:1394–1404. doi:10.1038/nn.4101 pmid:26404714
    OpenUrlCrossRefPubMed
  19. ↵
    1. Carmichael ST,
    2. Price JL
    (1995) Limbic connections of the orbital and medial prefrontal cortex in macaque monkeys. J Comp Neurol 363:615–641. doi:10.1002/cne.903630408 pmid:8847421
    OpenUrlCrossRefPubMed
  20. ↵
    1. Chudasama Y,
    2. Wright KS,
    3. Murray EA
    (2008) Hippocampal lesions in rhesus monkeys disrupt emotional responses but not reinforcer devaluation effects. Biol Psychiatry 63:1084–1091. doi:10.1016/j.biopsych.2007.11.012 pmid:18191111
    OpenUrlCrossRefPubMed
  21. ↵
    1. Dale AM,
    2. Fischl B,
    3. Sereno MI
    (1999) Cortical surface-based analysis: I. Segmentation and surface reconstruction. Neuroimage 9:179–194. doi:10.1006/nimg.1998.0395 pmid:9931268
    OpenUrlCrossRefPubMed
  22. ↵
    1. Desikan RS,
    2. Segonne F,
    3. Fischl B,
    4. Quinn BT,
    5. Dickerson BC,
    6. Blacker D,
    7. Buckner RL,
    8. Dale AM,
    9. Maguire RP,
    10. Hyman BT,
    11. Albert MS,
    12. Killiany RJ
    (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31:968–980. doi:10.1016/j.neuroimage.2006.01.021 pmid:16530430
    OpenUrlCrossRefPubMed
  23. ↵
    1. Evans DA,
    2. Stempel AV,
    3. Vale R,
    4. Branco T
    (2019) Cognitive control of escape behaviour. Trends Cogn Sci 23:334–348. doi:10.1016/j.tics.2019.01.012 pmid:30852123
    OpenUrlCrossRefPubMed
  24. ↵
    1. Evenden J,
    2. Ross L,
    3. Jonak G,
    4. Zhou J
    (2009) A novel operant conflict procedure using incrementing shock intensities to assess the anxiolytic and anxiogenic effects of drugs. Behav Pharmacol 20:226–236. doi:10.1097/fbp.0b013e32832a8110 pmid:19455771
    OpenUrlCrossRefPubMed
  25. ↵
    1. Fischl B,
    2. Sereno MI,
    3. Dale AM
    (1999a) Cortical surface-based analysis: II. Inflation, flattening, and a surface-based coordinate system. Neuroimage 9:195–207. doi:10.1006/nimg.1998.0396 pmid:9931269
    OpenUrlCrossRefPubMed
  26. ↵
    1. Fischl B,
    2. Sereno MI,
    3. Tootell RB,
    4. Dale AM
    (1999b) High-resolution intersubject averaging and a coordinate system for the cortical surface. Hum Brain Mapp 8:272–284. doi:10.1002/(SICI)1097-0193(1999)8:4<272::AID-HBM10>3.0.CO;2-4
    OpenUrlCrossRefPubMed
  27. ↵
    1. Fischl B,
    2. Salat DH,
    3. Busa E,
    4. Albert M,
    5. Dieterich M,
    6. Haselgrove C,
    7. van der Kouwe A,
    8. Killiany R,
    9. Kennedy D,
    10. Klaveness S,
    11. Montillo A,
    12. Makris N,
    13. Rosen B,
    14. Dale AM
    (2002) Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron 33:341–355. doi:10.1016/s0896-6273(02)00569-x pmid:11832223
    OpenUrlCrossRefPubMed
  28. ↵
    1. Fischl B,
    2. Rajendran N,
    3. Busa E,
    4. Augustinack J,
    5. Hinds O,
    6. Yeo BT,
    7. Mohlberg H,
    8. Amunts K,
    9. Zilles K
    (2008) Cortical folding patterns and predicting cytoarchitecture. Cereb Cortex 18:1973–1980. doi:10.1093/cercor/bhm225 pmid:18079129
    OpenUrlCrossRefPubMed
  29. ↵
    1. Gainotti G
    (2019) A historical review of investigations on laterality of emotions in the human brain. J Hist Neurosci 28:23–41. doi:10.1080/0964704X.2018.1524683 pmid:30475661
    OpenUrlCrossRefPubMed
  30. ↵
    1. Gordon JA,
    2. Lacefield CO,
    3. Kentros CG,
    4. Hen R
    (2005) State-dependent alterations in hippocampal oscillations in serotonin 1A receptor-deficient mice. J Neurosci 25:6509–6519. doi:10.1523/JNEUROSCI.1211-05.2005 pmid:16014712
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Gray J,
    2. McNaughton N
    (2000) The neuropsychology of anxiety, pp 72–82. Oxford: Oxford UP.
  32. ↵
    1. Heilbronner SR,
    2. Hayden BY
    (2016) Dorsal anterior cingulate cortex: a bottom-up view. Annu Rev Neurosci 39:149–170. doi:10.1146/annurev-neuro-070815-013952 pmid:27090954
    OpenUrlCrossRefPubMed
  33. ↵
    1. Holm S
    (1979) A simple sequentially rejective multiple test procedure. Scand J Stat 6:65–70.
    OpenUrl
  34. ↵
    1. Hutton C,
    2. Bork A,
    3. Josephs O,
    4. Deichmann R,
    5. Ashburner J,
    6. Turner R
    (2002) Image distortion correction in fMRI: a quantitative evaluation. Neuroimage 16:217–240. doi:10.1006/nimg.2001.1054 pmid:11969330
    OpenUrlCrossRefPubMed
  35. ↵
    1. Iglesias JE,
    2. Augustinack JC,
    3. Nguyen K,
    4. Player CM,
    5. Player A,
    6. Wright M,
    7. Roy N,
    8. Frosch MP,
    9. McKee AC,
    10. Wald LL,
    11. Fischl B,
    12. Van Leemput K
    , Alzheimer's Disease Neuroimaging Initiative. (2015) A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: application to adaptive segmentation of in vivo MRI. Neuroimage 115:117–137. doi:10.1016/j.neuroimage.2015.04.042 pmid:25936807
    OpenUrlCrossRefPubMed
  36. ↵
    1. Ito R,
    2. Lee AC
    (2016) The role of the hippocampus in approach–avoidance conflict decision-making: evidence from rodent and human studies. Behav Brain Res 313:345–357. doi:10.1016/j.bbr.2016.07.039 pmid:27457133
    OpenUrlCrossRefPubMed
  37. ↵
    1. Jhang J,
    2. Lee H,
    3. Kang MS,
    4. Lee HS,
    5. Park H,
    6. Han JH
    (2018) Anterior cingulate cortex and its input to the basolateral amygdala control innate fear response. Nat Commun 9:2744. doi:10.1038/s41467-018-05090-y pmid:30013065
    OpenUrlCrossRefPubMed
  38. ↵
    1. Khemka S,
    2. Barnes G,
    3. Dolan RJ,
    4. Bach DR
    (2017) Dissecting the function of hippocampal oscillations in a human anxiety model. J Neurosci 37:6869–6876. doi:10.1523/JNEUROSCI.1834-16.2017 pmid:28626018
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Kirlic N,
    2. Young J,
    3. Aupperle RL
    (2017) Animal to human translational paradigms relevant for approach avoidance conflict decision making. Behav Res Ther 96:14–29. doi:10.1016/j.brat.2017.04.010 pmid:28495358
    OpenUrlCrossRefPubMed
  40. ↵
    1. Kjelstrup KG,
    2. Tuvnes FA,
    3. Steffenach HA,
    4. Murison R,
    5. Moser EI,
    6. Moser MB
    (2002) Reduced fear expression after lesions of the ventral hippocampus. Proc Natl Acad Sci USA 99:10825–10830. doi:10.1073/pnas.152112399 pmid:12149439
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Korn CW,
    2. Bach DR
    (2015) Maintaining homeostasis by decision-making. PLOS Comput Biol 11:e1004301. doi:10.1371/journal.pcbi.1004301 pmid:26024504
    OpenUrlCrossRefPubMed
  42. ↵
    1. Korn CW,
    2. Bach DR
    (2018) Heuristic and optimal policy computations in the human brain during sequential decision-making. Nat Commun 9:325. doi:10.1038/s41467-017-02750-3 pmid:29362449
    OpenUrlCrossRefPubMed
  43. ↵
    1. Korn CW,
    2. Bach DR
    (2019) Minimizing threat via heuristic and optimal policies recruits hippocampus and medial prefrontal cortex. Nat Hum Behav 3:733–745. doi:10.1038/s41562-019-0603-9 pmid:31110338
    OpenUrlCrossRefPubMed
  44. ↵
    1. Korn CW,
    2. Vunder J,
    3. Miró J,
    4. Fuentemilla L,
    5. Hurlemann R,
    6. Bach DR
    (2017) Amygdala lesions reduce anxiety-like behavior in a human benzodiazepine-sensitive approach–avoidance conflict test. Biol Psychiatry 82:522–531. doi:10.1016/j.biopsych.2017.01.018
    OpenUrlCrossRefPubMed
  45. ↵
    1. Loh E,
    2. Kurth-Nelson Z,
    3. Berron D,
    4. Dayan P,
    5. Duzel E,
    6. Dolan R,
    7. Guitart-Masip M
    (2017) Parsing the role of the hippocampus in approach–avoidance conflict. Cereb Cortex 27:201–215. doi:10.1093/cercor/bhw378 pmid:27993819
    OpenUrlCrossRefPubMed
  46. ↵
    1. Luke SG
    (2017) Evaluating significance in linear mixed-effects models in R. Behav Res Methods 49:1494–1502. doi:10.3758/s13428-016-0809-y pmid:27620283
    OpenUrlCrossRefPubMed
  47. ↵
    1. Machado CJ,
    2. Kazama AM,
    3. Bachevalier J
    (2009) Impact of amygdala, orbital frontal, or hippocampal lesions on threat avoidance and emotional reactivity in nonhuman primates. Emotion 9:147–163. doi:10.1037/a0014539 pmid:19348528
    OpenUrlCrossRefPubMed
  48. ↵
    1. Maguire EA,
    2. Burgess N,
    3. Donnett JG,
    4. Frackowiak RS,
    5. Frith CD,
    6. Keefe J
    (1998) Knowing where and getting there: a human navigation network. Science 280:921–924. doi:10.1126/science.280.5365.921 pmid:9572740
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Mai JK,
    2. Majtanik M,
    3. Paxinos G
    (2016) Atlas of the human brain, Ed 4. San Diego: Academic.
  50. ↵
    1. McDonald AJ
    (1998) Cortical pathways to the mammalian amygdala. Prog Neurobiol 55:257–332. doi:10.1016/s0301-0082(98)00003-3 pmid:9643556
    OpenUrlCrossRefPubMed
  51. ↵
    1. McHugh SB,
    2. Deacon RM,
    3. Rawlins JN,
    4. Bannerman DM
    (2004) Amygdala and ventral hippocampus contribute differentially to mechanisms of fear and anxiety. Behav Neurosci 118:63–78. doi:10.1037/0735-7044.118.1.63 pmid:14979783
    OpenUrlCrossRefPubMed
  52. ↵
    1. Medford N,
    2. Critchley HD
    (2010) Conjoint activity of anterior insular and anterior cingulate cortex: awareness and response. Brain Struct Funct 214:535–549. doi:10.1007/s00429-010-0265-x pmid:20512367
    OpenUrlCrossRefPubMed
  53. ↵
    1. Mikl M,
    2. Mareček R,
    3. Hluštík P,
    4. Pavlicová M,
    5. Drastich A,
    6. Chlebus P,
    7. Brázdil M,
    8. Krupa P
    (2008) Effects of spatial smoothing on fMRI group inferences. Magn Reson Imaging 26:490–503. doi:10.1016/j.mri.2007.08.006 pmid:18060720
    OpenUrlCrossRefPubMed
  54. ↵
    1. Monosov IE
    (2017) Anterior cingulate is a source of valence-specific information about value and uncertainty. Nat Commun 8:134. doi:10.1038/s41467-017-00072-y pmid:28747623
    OpenUrlCrossRefPubMed
  55. ↵
    1. Mumford JA,
    2. Poline JB,
    3. Poldrack RA
    (2015) Orthogonalization of regressors in fMRI models. PLoS One 10:e0126255. doi:10.1371/journal.pone.0126255 pmid:25919488
    OpenUrlCrossRefPubMed
  56. ↵
    1. O'Neil EB,
    2. Newsome RN,
    3. Li IH,
    4. Thavabalasingam S,
    5. Ito R,
    6. Lee AC
    (2015) Examining the role of the human hippocampus in approach–avoidance decision making using a novel conflict paradigm and multivariate functional magnetic resonance imaging. J Neurosci 35:15039–15049. doi:10.1523/JNEUROSCI.1915-15.2015 pmid:26558775
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Oberrauch S,
    2. Sigrist H,
    3. Sautter E,
    4. Gerster S,
    5. Bach DR,
    6. Pryce CR
    (2019) Establishing operant conflict tests for the translational study of anxiety in mice. Psychopharmacology (Berl) 236:2527–2541. doi:10.1007/s00213-019-05315-y
    OpenUrlCrossRefPubMed
  58. ↵
    1. Padilla-Coreano N,
    2. Bolkan SS,
    3. Pierce GM,
    4. Blackman DR,
    5. Hardin WD,
    6. Garcia-Garcia AL,
    7. Spellman TJ,
    8. Gordon JA
    (2016) Direct ventral hippocampal-prefrontal input is required for anxiety-related neural activity and behavior. Neuron 89:857–866. doi:10.1016/j.neuron.2016.01.011 pmid:26853301
    OpenUrlCrossRefPubMed
  59. ↵
    1. Pitkänen A
    (2000) Connectivity of the rat amygdaloid complex. In: The amygdala: a functional analysis (JP A, ed), pp 31–116. Oxford: Oxford UP.
  60. ↵
    1. Robinson JL,
    2. Salibi N,
    3. Deshpande G
    (2016) Functional connectivity of the left and right hippocampi: evidence for functional lateralization along the long-axis using meta-analytic approaches and ultra-high field functional neuroimaging. Neuroimage 135:64–78. doi:10.1016/j.neuroimage.2016.04.022 pmid:27132046
    OpenUrlCrossRefPubMed
  61. ↵
    1. Sakaguchi Y,
    2. Sakurai Y
    (2017) Left–right functional asymmetry of ventral hippocampus depends on aversiveness of situations. Behav Brain Res 325:25–33. doi:10.1016/j.bbr.2017.02.028 pmid:28235588
    OpenUrlCrossRefPubMed
  62. ↵
    1. Schumacher A,
    2. Villaruel FR,
    3. Ussling A,
    4. Riaz S,
    5. Lee AC,
    6. Ito R
    (2018) Ventral hippocampal CA1 and CA3 differentially mediate learned approach–avoidance conflict processing. Curr Biol 28:1318–1324.e1314. doi:10.1016/j.cub.2018.03.012 pmid:29606418
    OpenUrlCrossRefPubMed
  63. ↵
    1. Segonne F,
    2. Grimson E,
    3. Fischl B
    (2005) A genetic algorithm for the topology correction of cortical surfaces. Inf Process Med Imaging 19:393–405.
    OpenUrlPubMed
  64. ↵
    1. Shenhav A,
    2. Cohen JD,
    3. Botvinick MM
    (2016) Dorsal anterior cingulate cortex and the value of control. Nat Neurosci 19:1286–1291. doi:10.1038/nn.4384 pmid:27669989
    OpenUrlCrossRefPubMed
  65. ↵
    1. Sladky R,
    2. Friston KJ,
    3. Trostl J,
    4. Cunnington R,
    5. Moser E,
    6. Windischberger C
    (2011) Slice-timing effects and their correction in functional MRI. Neuroimage 58:588–594. doi:10.1016/j.neuroimage.2011.06.078 pmid:21757015
    OpenUrlCrossRefPubMed
  66. ↵
    1. Spiers HJ,
    2. Burgess N,
    3. Maguire EA,
    4. Baxendale SA,
    5. Hartley T,
    6. Thompson PJ,
    7. O'Keefe J
    (2001) Unilateral temporal lobectomy patients show lateralized topographical and episodic memory deficits in a virtual town. Brain 124:2476–2489. doi:10.1093/brain/124.12.2476 pmid:11701601
    OpenUrlCrossRefPubMed
  67. ↵
    1. Strange BA,
    2. Witter MP,
    3. Lein ES,
    4. Moser EI
    (2014) Functional organization of the hippocampal longitudinal axis. Nat Rev Neurosci 15:655–669. doi:10.1038/nrn3785 pmid:25234264
    OpenUrlCrossRefPubMed
  68. ↵
    1. Swick D,
    2. Chatham CH
    (2014) Ten years of inhibition revisited. Front Hum Neurosci 8:329. doi:10.3389/fnhum.2014.00329 pmid:24904369
    OpenUrlCrossRefPubMed
  69. ↵
    1. Tzourio-Mazoyer N,
    2. Landeau B,
    3. Papathanassiou D,
    4. Crivello F,
    5. Etard O,
    6. Delcroix N,
    7. Mazoyer B,
    8. Joliot M
    (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15:273–289. doi:10.1006/nimg.2001.0978 pmid:11771995
    OpenUrlCrossRefPubMed
  70. ↵
    1. Wallis CU,
    2. Cockcroft GJ,
    3. Cardinal RN,
    4. Roberts AC,
    5. Clarke HF
    (2019) Hippocampal interaction with area 25, but not area 32, regulates marmoset approach–avoidance behavior. Cereb Cortex 29:4818–4830. doi:10.1093/cercor/bhz015
    OpenUrlCrossRefPubMed
  71. ↵
    1. Worsley KJ,
    2. Evans AC,
    3. Marrett S,
    4. Neelin P
    (1992) A three-dimensional statistical analysis for CBF activation studies in human brain. J Cereb Blood Flow Metab 12:900–918. doi:10.1038/jcbfm.1992.127 pmid:1400644
    OpenUrlCrossRefPubMed
  72. ↵
    1. Yeates DC,
    2. Ussling A,
    3. Lee AC,
    4. Ito R
    (2020) Double dissociation of learned approach–avoidance conflict processing and spatial pattern separation along the dorsoventral axis of the dentate gyrus. Hippocampus 30:596–609. doi:10.1002/hipo.23182
    OpenUrlCrossRef
  73. ↵
    1. Zorowitz S,
    2. Momennejad I,
    3. Daw ND
    (2020) Anxiety, avoidance, and sequential evaluation. Comput Psychiatry 4:1–17. doi:10.1162/cpsy_a_00026
    OpenUrlCrossRef
Back to top

In this issue

The Journal of Neuroscience: 40 (35)
Journal of Neuroscience
Vol. 40, Issue 35
26 Aug 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Hippocampal Representation of Threat Features and Behavior in a Human Approach–Avoidance Conflict Anxiety Task
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Hippocampal Representation of Threat Features and Behavior in a Human Approach–Avoidance Conflict Anxiety Task
Aslan Abivardi, Saurabh Khemka, Dominik R. Bach
Journal of Neuroscience 26 August 2020, 40 (35) 6748-6758; DOI: 10.1523/JNEUROSCI.2732-19.2020

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Hippocampal Representation of Threat Features and Behavior in a Human Approach–Avoidance Conflict Anxiety Task
Aslan Abivardi, Saurabh Khemka, Dominik R. Bach
Journal of Neuroscience 26 August 2020, 40 (35) 6748-6758; DOI: 10.1523/JNEUROSCI.2732-19.2020
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • anterior hippocampus
  • decision-making under predation
  • high-resolution fMRI
  • hippocampal subfields
  • lateral amygdala
  • operant conflict test

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Oxidative stress-induced damage to the developing hippocampus is mediated by GSK3beta
  • Disruption of endosomal sorting in Schwann cells leads to defective myelination and endosomal abnormalities observed in Charcot-Marie-Tooth disease
  • Depolarizing NaV and hyperpolarizing KV channels are co-trafficked in sensory neurons
Show more Research Articles

Behavioral/Cognitive

  • Spontaneous Alpha-Band Oscillations Bias Subjective Contrast Perception
  • The role of visual experience in individual differences of brain connectivity
  • A Neurodevelopmental Shift in Reward Circuitry from Mother's to Nonfamilial Voices in Adolescence
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2022 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.