Document Detail

Brain activation during fear conditioning in humans depends on genetic variations related to functioning of the hypothalamic-pituitary-adrenal axis: first evidence from two independent subsamples.
Jump to Full Text
MedLine Citation:
PMID:  22410078     Owner:  NLM     Status:  Publisher    
BACKGROUND: Enhanced acquisition and delayed extinction of fear conditioning are viewed as major determinants of anxiety disorders, which are often characterized by a dysfunctional hypothalamic-pituitary-adrenal (HPA) axis.MethodIn this study we employed cued fear conditioning in two independent samples of healthy subjects (sample 1: n=60, sample 2: n=52). Two graphical shapes served as conditioned stimuli and painful electrical stimulation as the unconditioned stimulus. In addition, guided by findings from published animal studies on HPA axis-related genes in fear conditioning, we examined variants of the glucocorticoid receptor and corticotropin-releasing hormone receptor 1 genes. RESULTS: Variation in these genes showed enhanced amygdala activation during the acquisition and reduced prefrontal activation during the extinction of fear as well as altered amygdala-prefrontal connectivity. CONCLUSIONS: This is the first demonstration of the involvement of genes related to the HPA axis in human fear conditioning.
S Ridder; J Treutlein; F Nees; S Lang; S Diener; M Wessa; A Kroll; S Pohlack; R Cacciaglia; P Gass; G Schütz; G Schumann; H Flor
Related Documents :
22479568 - The flagellar regulator flit represses salmonella pathogenicity island 1 through flhdc ...
22582958 - Redox-responsive regulation of denitrification genes in brucella.
22479468 - Functional identification and characterization of the brassica napus transcription fact...
22390678 - New prognostic markers, determined using gene expression analyses, reveal two distinct ...
2458298 - Dna methylation: a secondary event in globin gene switching?
11282888 - Gene expression patterns in the lungs of patients with primary pulmonary hypertension: ...
Publication Detail:
Type:  JOURNAL ARTICLE     Date:  2012-3-12
Journal Detail:
Title:  Psychological medicine     Volume:  -     ISSN:  1469-8978     ISO Abbreviation:  -     Publication Date:  2012 Mar 
Date Detail:
Created Date:  2012-3-13     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  1254142     Medline TA:  Psychol Med     Country:  -    
Other Details:
Languages:  ENG     Pagination:  1-11     Citation Subset:  -    
Department of Cognitive and Clinical Neuroscience, Central Institute of Mental Health, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Psychol Med
Journal ID (iso-abbrev): Psychol Med
Journal ID (publisher-id): PSM
ISSN: 0033-2917
ISSN: 1469-8978
Publisher: Cambridge University Press, Cambridge, UK
Article Information
Download PDF
Copyright © Cambridge University Press 2012 The online version of this article is published within an Open Access environment subject to the conditions of the Creative Commons Attribution-NonCommercial-ShareAlike licence <>. The written permission of Cambridge University Press must be obtained for commercial re-use.
Received Day: 27 Month: 7 Year: 2011
Revision Received Day: 04 Month: 2 Year: 2012
Accepted Day: 06 Month: 2 Year: 2012
Print publication date: Month: 11 Year: 2012
Electronic publication date: Day: 12 Month: 3 Year: 2012
Volume: 42 Issue: 11
First Page: 2325 Last Page: 2335
PubMed Id: 22410078
ID: 3466052
DOI: 10.1017/S0033291712000359
Publisher Item Identifier: S0033291712000359
Publisher Id: 00035

Brain activation during fear conditioning in humans depends on genetic variations related to functioning of the hypothalamic–pituitary–adrenal axis: first evidence from two independent subsamples Alternate Title:S. Ridder et al. Alternate Title:Brain activation during fear conditioning
S. Ridder1
J. Treutlein2
F. Nees1
S. Lang3
S. Diener1
M. Wessa1
A. Kroll1
S. Pohlack1
R. Cacciaglia1
P. Gass4
G. Schütz5
G. Schumann6
H. Flor1*
1Department of Cognitive and Clinical Neuroscience, Central Institute of Mental Health, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
2Department of Genetic Epidemiology in Psychiatry, Central Institute of Mental Health, Mannheim, Germany
3Institute of Psychology, Department of Clinical Psychology and Psychotherapy, Heidelberg University, Heidelberg, Germany
4Research Group Behavioral Biology, Central Institute of Mental Health, Mannheim, Germany
5Division of Molecular Biology of the Cell I, German Cancer Research Center, Heidelberg, Germany
6Section of Addiction Biology, Division of Psychological Medicine and Psychiatry, Institute of Psychiatry, King's College London, London, UK
*Address for correspondence: H. Flor, Ph.D., Department of Cognitive and Clinical Neuroscience, Central Institute of Mental Health, Medical Faculty Mannheim, Heidelberg University, Square J5, 68159 Mannheim, Germany. (Email:


One of the best-understood and most studied learning paradigms is classical cued fear conditioning where a neutral stimulus is repeatedly paired with an aversive unconditioned stimulus (US) such as an electric shock. Over time the neutral stimulus alone elicits a learnt or conditioned fear response that resembles innate or unconditioned responses. Repeated presentation of the unreinforced conditioned stimulus (CS) leads to a decrease of the acquired fear response or extinction. A key brain region for the acquisition of conditioned fear is the amygdala (LeDoux, 2000) whereas extinction has been related to the inhibition of amygdalar activation by the prefrontal cortex (Herry et al.2010). In post-traumatic stress disorder (PTSD) and other anxiety disorders the enhanced acquisition and delayed extinction of fear responses are viewed as a major aetiological mechanism (Bremner et al.2005; Delgado et al.2006; Rauch et al.2006; Blechert et al.2007). Recent evidence further suggests that insufficient prefrontal activation during extinction might be important for the maintenance of fear symptoms and stimulation of these brain areas could have therapeutic effects (Amano et al.2010). Anxiety disorders such as PTSD are also characterized by alterations in the function of the hypothalamic–pituitary–adrenal (HPA) axis to a hypoactive mode with reduced basal and stress-related cortisol levels as well as increased suppression after the dexamethasone-suppression test (Yehuda, 2006; Wessa & Rohleder, 2007).

The heritability of the components of fear conditioning has been estimated at 35–45% (Hettema et al.2003) and several genes involved in the regulation of plasticity and emotional reactivity have been found associated with the conditionability of fear responses (Lonsdorf et al.2009). Besides altered fear conditioning, anxiety disorders such as PTSD are also characterized by a dysfunction of the HPA axis, which is regulated by corticotropin-releasing hormone (CRH) and glucocorticoids. They have been implicated in fear conditioning in animal models as well as in PTSD in humans (de Kloet et al.2008; Rodrigues et al.2009). For example, CRH injection into the central nucleus of the amygdala and CRH injection into the medial prefrontal cortex both result in increased anxiety-related behaviour (Timpl et al.1998), while mice lacking the CRH receptor (CRHR1) gene display reduced anxiety-related behaviour (Müller et al.2003). In PTSD, elevation of CRH was found in the cerebrospinal fluid, indicating a possible role in altered brain activity (Baker et al.1999; Jaferi & Bhatnagar, 2007; de Kloet et al.2008) in fear learning. These observations strengthen the assumption of a tight connection between the amygdaloid CRH system and HPA functioning in fear conditioning as well as in stress-related disorders, such as PTSD, in humans.

The glucocorticoid receptor (GR) is a main regulator of the HPA axis and has been shown to influence endocrine and behavioural measures of fear in various animal studies (Bremner et al.1997; de Kloet et al.1998). Pharmacological GR (NR3C1) manipulations, for example, of amygdaloid GRs were shown to modulate processes of fear conditioning (Yang et al.2006) and to be related to the establishment of fear memories (e.g. Oitzl et al.1998), and to facilitated extinction of conditioned fear (Yang et al.2006). GR receptors in the prefrontal cortex might be related to a failure to extinguish established fear responses (Tronel & Alberini, 2007) and to an enhancement of emotional memory consolidation, suggesting that these effects reflect an interaction of the medial prefrontal cortex and the basolateral amygdala (Roozendaal et al.2009).

However, it is not known if the reported changes in brain activity during fear conditioning and the altered HPA axis activity seen in PTSD patients reflect consequences of the traumatic event or might be interpreted as predisposing vulnerability factors increasing the susceptibility for the disorder. The latter has been proposed in a study that examined firefighters before the experience of traumatic events. Reduced extinction of an aversively conditioned corrugator electromyographic response predicted more than 30% of the variance of PTSD symptoms following trauma (Guthrie & Bryant, 2006). A twin study conducted with combat veterans and their non-combat-exposed sibling revealed that structural changes in the brain are pre-existing familial vulnerability factors (Gilbertson et al.2002; Pitman et al.2006). As a formal genetic study attributes moderate heritability to all components of the fear-conditioning process (Hettema et al.2003), it is of interest to identify the genetic variations responsible at the molecular level. The identified genes could be regarded as vulnerability factors making subjects more susceptible for stress-related disorders.

We used an imaging genetics approach in two independent samples (sample 1: n=60, sample 2: n=52) of healthy individuals to investigate the influence of genetic variation of HPA axis-related genes [CRH receptor 1 (CRHR1), GR (NR3C1)] on the acquisition and extinction of cued fear responses. The products of both genes are among the most important molecules of the HPA axis and were repeatedly reported to be involved in processes of fear conditioning in animal studies (CRHR1, e.g. Radulovic et al.1999, Kikusui et al.2000, Otagiri et al.2000; NR3C1, e.g. Cordero et al.2002, Yang et al.2006, Kohda et al.2007, Kolber et al.2008), which was the reason to include these two in the present study. In addition, the first genetic findings for PTSD exist for the Bcl1 variant (rs41423247) of the GR gene. The GG genotype of this single nucleotide polymorphism (SNP) was associated with low basal cortisol levels in PTSD (Bachmann et al.2005) and with more long-term traumatic memories and higher PTSD symptom scores (Hauer et al.2011). We employed functional magnetic resonance imaging (fMRI) and used a linear regression approach (according to Ressler et al.2011), which permitted the comparison of groups with respect to degree of genetic variation with and without minor alleles. We examined amygdala activation as an indicator for acquisition and prefrontal activation as an indicator for extinction and also included connectivity analyses. To control for successful conditioning, we assessed skin conductance responses (SCRs) and self-report measures.


A total of sixty persons [38 male, 22 female, mean age 21.25 (S.D.=3.02, range 19–37) years, all right-handed] recruited in schools for ambulance rescue workers as part of a longitudinal study were examined in study 1 (sample 1) and 52 persons [32 male, 20 female, mean age 22.27 (S.D.=3.67, range 18–37) years, all right-handed] from the same population were recruited for study 2 (sample 2). Participants with past traumatic events as assessed by the German version of the Posttraumatic Diagnostic Scale (Griesel et al.2006) were excluded. All participants were medication free and had no physical or mental disorders. The study was approved by the Ethics Committee of the Medical Faculty Mannheim, Heidelberg University, and written informed consent was obtained from all participants, who were paid for participation.

Psychological assessment

Standardized clinical assessment with the Structured Clinical Interview for DSM-IV Axis I (First et al.1997b) and Axis II (First et al.1997a) was performed to exclude persons with mental disorders. Psychological assessment was identical for all subjects and included the Center for Epidemiological Studies Depression Scale (Radloff, 1977), the trait section of the State-Trait Anxiety Inventory (Spielberger et al.1970), the Perceived Stress Scale (Cohen et al.1983) and the Childhood Trauma Questionnaire (Bernstein & Fink, 1998) (see Table 1). Since we did not find any differences in scores on the psychological measures depending on the genotype, we did not focus on these data in the Results section.

SCR and self-report data

The SCRs were recorded from two electrodes placed on the thenar and hypothenar eminence of the participants' right hand using a sampling rate of 16 Hz and a VarioPort recording system (BECKER MEDITEC, Germany). Data analysis was performed using EDA-PARA software (F. Schäfer, Germany) and followed the guidelines of Fowles et al. (1981) . Trials were visually inspected for artifacts and SCR amplitudes were quantified as the maximum response in the time window of 1–4 s (first interval response) and 5–9 s (second interval response; Prokasy & Ebel, 1967) after stimulus onset and were measured in microSiemens (μS). SCR amplitudes below 0.05 μS were classified as zero responses. SCR data were normalized using a logarithmic [log (1+SCR)] transformation.

After each conditioning phase, participants verbally rated the emotional valence and arousal of the CSs (1=very calm to 9=very arousing, 1=very pleasant to 9=very unpleasant) as well as the CS–US contingency (1=no CS–US contingency to 9=perfect CS–US contingency). All auditory or visual instructions for the experimental procedure were standardized. Communication was realized via headphones with attached microphones.

SCRs and self-report data were analysed separately using Predictive Analytic Software (PASW) for Windows, version 18.0.1 (SPSS Inc., USA). Both SCRs and self-reports showed successful conditioning and extinction in samples 1 and 2. Since differences in the genotype groups could not be observed for either measure, we present only the fMRI analyses in the Results section.

DNA extraction, selection of SNPs and genotyping

Venous blood samples were obtained from all participants. Genomic DNA was isolated with the QIAamp DNA extraction kit ( For genetic characterization of the NR3C1 and CRHR1 genes, we selected SNPs with potential functionality from the literature as well as tagging SNPs from the HapMap database and literature. For the NR3C1 gene, we chose the potentially functional variants N363S (rs6195) (e.g. Jewel & Cidlowski, 2007), BclI (rs41423247) (Stevens et al.2004) and Tth111I (rs10052957) (Rosmond et al.2000). Tagging SNPs for the NR3C1 gene were selected by a blockwise strategy from HapMap data, using haplotypes above 5% frequency in HaploView (e.g. Barrett et al.2005). NR3C1 transcript NM_000176, which covers 123.8 kbp on chr5, contained only one large haplotype block, which is tagged by four haplotype tagging SNPs, i.e. rs33389, rs4986593, rs10482672 and rs190488 (HapMap Rel 16c, NCBI B34 assembly, dbSNP b124). Tagging SNPs for the CRHR1 gene, rs1876831 and rs242938, were selected from the literature, based on detailed linkage disequilibrium information of CRHR1 SNPs from several publications (e.g. Treutlein et al.2006; Wassermann et al.2009; Grabe et al.2010; Nelson et al.2010). The tagging markers were shown to be sufficient to provide information on the brain activation during fear conditioning and captured 66% of the markers of the NR3C1 and 45% of SNPs of the CRHR1 genes according to HapMap release 24 (threshold r2 ⩾ 0.8, allele frequency ⩾5%). Genotyping was performed using an Applied Biosystems Prism 7900HT RealTime PCR system. Table S1 (see Supplementary material) lists the IDs of Applied Biosystems Assays-on-DemandTM and Primer/Probe sequences of Assays-by-Design. Allele frequencies and genotype counts for the polymorphisms that contributed significantly to the association signals in the two samples are shown in Table 2.

Design of the fear-conditioning experiment

Perception and pain threshold as well as pain tolerance (three repetitions in an ascending series) were assessed for the calculation of the intensity of the US, which was painful electrical stimulation (Digitimer, UK) at the thumb of the right hand set to a level of 80% of pain tolerance. Two different geometric shapes (a square and a rhombus) in different colours (blue or yellow) were presented as cue stimuli with colours and shapes counterbalanced across subjects. During habituation 10 CS+ (neutral stimulus that would later be followed by the aversive US), 10 CS− (neutral stimulus that would later be followed by the absence of the US) and four US were presented in random order. The acquisition was divided into two phases, each consisting of nine CS+ paired with the US, nine CS+ not paired with the US and 18 CS− (never paired with the US) in random order. The extinction phase consisted of 18 CS+ and 18 CS− trials. The CS+ was reinforced in 50% of the trials with a shock duration of 2.7 s before the end of the CS+ projection (see Fig. 1). After each phase the participants rated valence and arousal on a self-assessment manikin. CS–US contingency was rated on a scale from 0 (US will not follow the CS) to 100 (US will definitely follow the CS).

Statistical analysis

Effects of genetic markers on amygdala activation during fear conditioning and medial prefrontal activation during fear extinction were assessed using a linear regression approach (Ressler et al.2011). SNPs were coded using an additive model, i.e. using the number of minor alleles of the respective markers as predictors. Age and gender were also tested but revealed no significant effects and were thus not included in the further analyses. For the first sample a model was built including all of the above-mentioned SNPs as predictors. For the second sample the same model was tested, but this time omitting the markers that had not obtained a nominally significant p value in the first sample. In these models the minor alleles of the nominally significant markers uniformly showed up as associated alleles. In order to build a summary score of all SNPs, genotype was coded by the total number of minor alleles across all markers. Scores were built across both genes and for the NR3C1 gene separately. This permitted the comparison of groups with respect to the degree of genetic variation with and without minor alleles for brain activation (see the legend of Fig. 2 for details). All significance levels were set to p<0.05.


Neuroimaging was performed during classical aversive delay cued conditioning in a 1.5 T Magnetom Vision scanner (Siemens Medical Solutions, Germany). Contiguous transversal T2*-weighted echo-planar images (EPI) with blood oxygenation level-dependent (BOLD) contrast were used (echo time 45 ms, flip angle 90°) that covered the whole brain (35 slices, slice thickness 3 mm, 1 mm gap, field of view=220×220 mm2, 64×64 matrix). The effective repetition time was 3.77 s per volume. A total of 560 volumes were recorded. Statistical parametric mapping (SPM) software was used for image processing and analysis (SPM2, The images were slice-time corrected for phase shift during volume acquisition, realigned to the first image for motion artefacts, spatially normalized to a standard EPI template, spatially smoothed with an isotropic Gaussian kernel with a full width at half-maximum of 10 mm, and temporal high-pass filtered (cut-off 128 s). Specific effects were tested by applying linear contrasts to the parameter estimated for each event. Contrast images of interest were calculated for each subject, and the resulting contrast images were entered into a second-level random-effects analysis to produce group results (one-sample t test). For each phase (early acquisition, late acquisition, and extinction) contrasts between CS+ and CS−, i.e. CS+>CS−, were calculated. We used a threshold of p<0.001 for the entire brain (uncorrected, extent threshold k=5 voxels). Additionally, according to our a priori hypothesis that alterations of amygdala activation during acquisition and prefrontal cortex activation during extinction should be associated with HPA axis-related genes (e.g. Yang et al.2006; Tronel & Alberini, 2007), we adopted a region-of-interest (ROI) approach with small volume correction [p<0.05, family-wise error (FWE) corrected], also corrected for the number of ROIs. Regions were defined using the MARINA software package ( To detect the association between genotype and fMRI activation in the amygdala and in the prefrontal cortex on a voxel-by-voxel basis, the contrast images of all subjects (percentage signal change of CS+ versus CS−) were included in a regression analysis with SPM.

Functional coupling analysis

Connectivity analyses were performed to determine genotype-dependent changes in functional coupling between the amygdala, the hippocampus and the prefrontal cortex that are related to the acquisition and extinction of the learned response. Functional coupling between seed regions (spheres with radius=6 mm, coordinates based on the current sample for acquisition phase: x, y, z: −24, −6, −18; for extinction phase: x, y, z: −36, 60, 0) and target regions was determined for the early acquisition phase and the extinction phase separately using standard SPM methods. For the early acquisition phase the left amygdala served as the seed region and the prefrontal cortex as the target region. For the extinction phase, the left prefrontal cortex [Brodmann area (BA) 10] (seed region) and left amygdala (target region) were used for the functional coupling analysis.

From both seed regions fMRI time series were extracted and used as regressors in a subsequent single-subject analysis, where the movement-related covariates were additionally included. These contrasts were used to carry out a random-effects analysis to determine functional coupling between groups assigned to different genotypes using two-sample t tests.

We report T values small-volume corrected using the following ROIs: prefrontal cortex (coordinates x=0, y=52, z=−3; sphere with radius=9 mm, see Heinz et al.2007) and left amygdala (created using the MARINA software package). Peak activations were correlated with subjective and endocrinological variables using Pearson correlations.


Acquisition- and extinction-related significant BOLD signal differences between CS+ and CS− were found in a number of brain regions (see Supplementary Table S2 and Fig. 2).

A linear regression (Ressler et al.2011) including all tested variants yielded highly significant associations of the minor alleles of several SNPs with fear conditioning. These SNPs were then included in a score-based analysis (see Fig. 2 and Table 3). We observed in the left, but not the right, amygdala a positive correlation between the peak BOLD signal change in the amygdala elicited by CS+ (danger signal) versus CS− (safety signal) during early acquisition and the number of minor alleles of NR3C1 rs33389, NR3C1 Bcl1 (rs41423247) and NR3C1 rs4986593 [r=0.48, F(3, 55)=6.01, p=0.001, see also Fig. 2]. For extinction, lack of differential activation in the prefrontal cortex was associated with an increasing number of minor alleles of CRHR1 rs242938, CRHR1 rs1876831, NR3C1 rs6195 and NR3C1 Tth111I (rs10052957), suggesting impaired extinction with increased genetic minor allele variants [r=0.60, F(4, 54)=6.24, p<0.001].

During early acquisition, the subjects with several minor alleles showed significantly more functional coupling of the BOLD signal between the left amygdala and prefrontal cortex compared with those with only one or no minor allele on the NR3C1 gene. For extinction, functional coupling strength between the left prefrontal cortex and left amygdala was most pronounced for the subjects with at least two minor alleles on the NR3C1 and CRHR1 genes. No significant association with single genetic markers was seen.

We sought to replicate the association of these polymorphisms with amygdala and prefrontal cortex activation during acquisition and extinction in a second independent sample of 52 individuals from the same population (for sample characteristics, see Tables 1 and 2). Of the three NR3C1 polymorphisms, two were again significant predictors of differential amygdala activation in the early acquisition phase [NR3C1 rs33389, NR3C1 rs4986593: r=0.45, F(2, 42)=5.45, p=0.008] when all variables were entered, but NR3C1 Bcl1 (rs41423247) was also significant by itself when a stepwise solution was employed [r=0.41, F(1, 44)=8.82, p=0.005]. The same four SNPs [CRHR1 rs242938, CRHR1 rs1876831, NR3C1 N363S (rs6195), NR3C1 Tth111I (rs10052957)] were again significant for differential activation (i.e. delayed extinction) in the prefrontal cortex (BA10) during extinction [r=0.47, F(4, 43)=2.79, p=0.04].

These findings could also be found when combining both samples. Genetic variants were unrelated to skin conductance or self-report data of fear conditioning.


The present study is the first to show a significant contribution of HPA axis-related genes to brain activation during fear conditioning and extinction in humans. Thus, more minor alleles in the NR3C1 gene were related to enhanced amygdala activation during conditioning and more minor alleles in the NR3C1 and CRHR1 genes were associated with reduced prefrontal activation during extinction, suggesting an accumulation of the effects. In addition, enhanced amygdala–prefrontal connectivity during acquisition, suggesting better fear memory consolidation, and higher amygdala–prefrontal interaction during extinction, suggesting sustained fear, were genotype related.

The development of PTSD after trauma exposure might depend on altered processes of fear learning. An inability to habituate to aversive stimuli and a reduction in inhibition of fear memories have often been reported (Rauch et al.2006; Yehuda & LeDoux, 2007; Jovanovic et al. 2010; Jovanovic & Ressler, 2010; Shin & Liberzon, 2010). The amygdala was previously shown to be involved in the elaboration of conditioned fear responses (Davis, 1992), while the prefrontal cortex was shown to regulate the activity of the amygdala in a top-down process and to inhibit the extinction of conditioned fear responses (Milad & Quirk, 2002; Peters et al.2009). The present data are in accordance with animal data that reported a close association of CRH1 receptors in the amygdala and fear conditioning (Yang et al.2006) and a close association of GRs and fear acquisition (Kolber et al.2008). The finding that significant amygdala activation was only present during early acquisition also supports previous reports on time-dependent activation of the amygdala during conditioning (e.g. Büchel et al.1998). Thus, the amygdala seems to play a larger role in the initiation than the maintenance of the response. In addition, we did not find similar genotype-related differences for subjective indicators of conditioning or skin conductance as we found for brain activation. This is not unusual and related to the fact that brain activation may be more closely related to genetic factors than subjective ratings or peripheral responses and may thus be more suitable as an intermediate phenotype (e.g. Meyer-Lindenberg, 2010). Our data strengthen the hypothesis that HPA axis functioning and stress play a role in the development of anxiety disorders (McFarlane et al.2011) and are in accordance with the findings described above that genes involved in glucocorticoid signalling are differentially expressed in PTSD (Yehuda et al.2009). The observed pattern corresponds to the deficits in patients with PTSD who retain responses to trauma cues, fail to extinguish cue-relevant associations and show altered genotypes of HPA axis-related genes (Rauch et al.2006; Yehuda et al.2010).

Although the minor allele cannot be regarded as relevant to the analysed conditioning and clinical characteristics a priori, the distribution of alleles in the model applied to the first sample showed that the minor alleles accumulated for the observed neural activation patterns during conditioning. The reason for our observation that the associated alleles are the minor alleles (allele frequency of ⩽40.4%) is unclear. From an evolutionary point of view, both minor and major alleles can confer susceptibility to a trait, and based on selection for a trait or drift, the allele frequency of a marker is subject to change over time (Jobling et al.2003; Keller & Miller, 2006). However, for two of the analysed markers [NR3C1 Bcl1 (van Rossum et al.2003), NR3C1 N363S (Huizenga et al.1998)] the minor alleles have been reported to increase glucocorticoid sensitivity (which was suggested to occur in PTSD, e.g. Rohleder et al.2004; Yehuda et al.2004), but this seems not to apply to all tissues (Koper et al.1997; Kumsta et al.2008). This strengthens our observation of the association of the minor alleles with trait-like factors of increased fear memory consolidation and failed extinction reflected by amygdala–prefrontal activation patterns.


Support for this study was provided by grant no. SFB636/C1 from the Deutsche Forschungsgemei-nschaft to H.F.

Supplementary material

For supplementary material accompanying this paper, visit

Declaration of Interest


Amano T,Unal CT,Paré D, (. Year: 2010Synaptic correlates of fear extinction in the amygdalaNature Neuroscience13489494
Bachmann AW,Sedgley TL,Jackson RV,Gibson JN,Young RM,Torpy DJ, (. Year: 2005Glucocorticoid receptor polymorphisms and post-traumatic stress disorderPsychoneuroendocrinology3029730615511603
Baker DG,West SA,Nicholson WE,Ekhator NN,Kasckow JW,Hill KK,Bruce AB,Orth DN,Geracioti TD Jr., (. Year: 1999Serial CSF corticotropin-releasing hormone levels and adrenocortical activity in combat veterans with posttraumatic stress disorderAmerican Journal of Psychiatry15658558810200738
Barrett JC,Fry B,Maller J,Daly MJ, (. Year: 2005Haploview: analysis and visualization of LD and haplotype mapsBioinformatics2126326515297300
Bernstein DP,Fink L, (. Year: 1998Childhood Trauma Questionnaire: A Retrospective Self-Report Questionnaire and ManualThe Psychological CorporationSan Antonio, TX
Blechert J,Michael T,Vriends N,Margraf J,Wilhelm FH, (. Year: 2007Fear conditioning in posttraumatic stress disorder: evidence for delayed extinction of autonomic, experiential, and behavioural responsesBehavioral Research Therapy4520192033
Bremner JD,Licinio J,Darnell A,Krystal JH,Owens MJ,Southwick SM,Nemeroff CB,Charney DS, (. Year: 1997Elevated CSF corticotropin-releasing factor concentrations in posttraumatic stress disorderAmerican Journal of Psychiatry1546246299137116
Bremner JD,Vermetten E,Schmahl C,Vaccarino V,Vythilingam M,Afzal N,Grillon C,Charney DS, (. Year: 2005Positron emission tomographic imaging of neural correlates of a fear acquisition and extinction paradigm in women with childhood sexual-abuse-related post-traumatic stress disorderPsychological Medicine3579180615997600
Büchel C,Morris J,Dolan RJ,Friston KJ, (. Year: 1998Brain systems mediating aversive conditioning: an event-related fMRI studyNeuron209479579620699
Cordero MI,Kruyt ND,Merino JJ,Sandi C, (. Year: 2002Glucocorticoid involvement in memory formation in a rat model for traumatic memoryStress5737912171770
Cohen S,Kamarck T,Mermelstein R, (. Year: 1983A global measure of perceived stressJournal of Health and Social Behavior243853966668417
Davis M, (. Year: 1992The role of the amygdala in fear and anxietyAnnual Review of Neuroscience15353375
de Kloet CS,Vermetten E,Geuze E,Lentjes EG,Heijnen CJ,Stalla GK,Westenberg HG, (. Year: 2008Elevated plasma corticotrophin-releasing hormone levels in veterans with posttraumatic stress disorderProgress in Brain Research16728729118037027
de Kloet ER,Vreugdenhil E,Oitzl MS,Joëls M, (. Year: 1998Brain corticosteroid receptor balance in health and diseaseEndocrinological Review19269301
Delgado MR,Olsson A,Phelps EA, (. Year: 2006Extending animal models of fear conditioning to humansBiological Psychology73394816472906
First MB,Gibbon M,Spitzer RL,Williams JBW,Benjamin LS, (. Year: 1997aUsers Guide for the Structured Clinical Interview for DSM-IV Axis II Personality Disorders (SCID-II)American Psychiatric PressWashington, DC
First MB,Spitzer RL,Gibbon M,Williams JBW, (. Year: 1997bUsers Guide for the Structured Clinical Interview for DSM-IV Axis I Disorders (SCID-I)-Clinical VersionAmerican Psychiatric PressWashington, DC
Fowles DC,Christie MJ,Edelberg R,Grings WW,Lykken DT,Venables PH, (. Year: 1981Committee report. Publication recommendations for electrodermal measurementsPsychophysiology182322397291438
Gilbertson MW,Shenton ME,Ciszewski A,Kasai K,Lasko NB,Orr SP,Pitman RK, (. Year: 2002Smaller hippocampal volume predicts pathologic vulnerability to psychological traumaNature Neuroscience512421247
Grabe HJ,Schwahn C,Appel K,Mahler J,Schulz A,Spitzer C,Fenske K,Barnow S,Lucht M,Freyberger HJ,John U,Teumer A,Wallaschofski H,Nauck M,Völzke H, (. Year: 2010Childhood maltreatment, the corticotropin-releasing hormone receptor gene and adult depression in the general populationAmerican Journal of Medical Genetics Part B Neuropsychiatric Genetics153B14831493
Griesel D,Wessa M,Flor H, (. Year: 2006Psychometric qualities of the German version of the Posttraumatic Diagnostic Scale (PTDS)Psychological Assessment1826226816953729
Guthrie RM,Bryant RA, (. Year: 2006Extinction learning before trauma and subsequent posttraumatic stressPsychosomatic Medicine6830731116554398
Hauer D,Weis F,Papassotiropoulos A,Schmoeckel M,Beiras-Fernandez A,Lieke J,Kaufmann I,Kirchhoff F,Vogeser M,Roozendaal B,Briegel J,de Quervain D,Schelling G, (. Year: 2011Relationship of a common polymorphism of the glucocorticoid receptor gene to traumatic memories and posttraumatic stress disorder in patients after intensive care therapyCritical Care Medicine3964365021169818
Heinz A,Wrase J,Kahnt T,Beck A,Bromand Z,Grüsser SM,Kienast T,Smolka MN,Flor H,Mann K, (. Year: 2007Brain activation elicited by affectively positive stimuli is associated with a lower risk of relapse in detoxified alcoholic subjectsAlcoholism: Clinical and Experimental Research3111381147
Herry C,Ferraguti F,Singewald N,Letzkus JJ,Ehrlich I,Lüthi A, (. Year: 2010Neuronal circuits of fear extinctionEuropean Journal of Neuroscience3159961220384807
Hettema JM,Annas P,Neale MC,Kendler KS,Fredrikson M, (. Year: 2003A twin study of the genetics of fear conditioningArchives of Genetic Psychiatry60702708
Huizenga NA,Koper JW,De Lange P,Pols HA,Stolk RP,Burger H,Grobbee DE,Brinkmann AO,De Jong FH,Lamberts SW, (. Year: 1998A polymorphism in the glucocorticoid receptor gene may be associated with an increased sensitivity to glucocorticoids in vivoJournal of Clinical Endocrinology and Metabolism831441519435432
Jaferi A,Bhatnagar S, (. Year: 2007Corticotropin-releasing hormone receptors in the medial prefrontal cortex regulate hypothalamic–pituitary–adrenal activity and anxiety-related behavior regardless of prior stress experienceBrain Research118621222318001698
Jewell CM,Cidlowski JA, (. Year: 2007Molecular evidence for a link between the N363S glucocorticoid receptor polymorphism and altered gene expressionJournal of Clinical Endocrinology and Metabolism923268327717535992
Jobling MA,Hurles M,Tyler-Smith C, (. Year: 2003Human Evolutionary Genetics: Origins, Peoples and DiseaseGarland ScienceNew York
Jovanovic T,Ressler KJ, (. Year: 2010How the neurocircuitry and genetics of fear inhibition may inform our understanding of PTSDAmerican Journal of Psychiatry664866220231322
Jovanovic T,Norrholm SD,Blanding NQ,Davis M,Duncan E,Bradley B,Ressler KJ, (. Year: 2010Impaired fear inhibition is a biomarker of PTSD but not depressionDepression & Anxiety2724425120143428
Keller MC,Miller G, (. Year: 2006Resolving the paradox of common, harmful, heritable mental disorders: which evolutionary genetic models work best?Behavioral and Brain Sciences2938540417094843
Kikusui T,Takeuchi Y,Mori Y, (. Year: 2000Involvement of corticotropin-releasing factor in the retrieval process of fear-conditioned ultrasonic vocalization in ratsPhysiology and Behavior7132332811150564
Kohda K,Harada K,Kato K,Hoshino A,Motohashi J,Yamaji T,Morinobu S,Matsuoka N,Kato N, (. Year: 2007Glucocorticoid receptor activation is involved in producing abnormal phenotypes of single-prolonged stress rats: a putative post-traumatic stress disorder modelNeuroscience148223317644267
Kolber BJ,Roberts MS,Howell MP,Wozniak DF,Sands MS,Muglia LJ, (. Year: 2008Central amygdala glucocorticoid receptor action promotes fear-associated CRH activation and conditioningProceedings of the National Academy of Sciences USA1051200412009
Koper JW,Stolk RP,de Lange P,Huizenga NA,Molijn GJ,Pols HA,Grobbee DE,Karl M,de Jong FH,Brinkmann AO,Lamberts SW, (. Year: 1997Lack of association between five polymorphisms in the human glucocorticoid receptor gene and glucocorticoid resistanceHuman Genetics996636689150737
Kumsta R,Entringer S,Koper JW,van Rossum EF,Hellhammer DH,Wüst S, (. Year: 2008Glucocorticoid receptor gene polymorphisms and glucocorticoid sensitivity of subdermal blood vessels and leukocytesBiological Psychology7917918418502562
LeDoux JE, (. Year: 2000Emotion circuits in the brainAnnual Review of Neuroscience23155184
Lonsdorf TB,Weike AI,Nikamo P,Schalling M,Hamm AO,Ohman A, (. Year: 2009Genetic gating of human fear learning and extinction: possible implications for gene–environment interaction in anxiety disorderPsychological Sciences20198206
McFarlane AC,Barton CA,Yehuda R,Wittert G, (. Year: 2011Cortisol response to acute trauma and risk of posttraumatic stress disorderPsychoneuroendocrinology3672072721093988
Meyer-Lindenberg A, (. Year: 2010Intermediate or brainless phenotypes for psychiatric research?Psychological Medicine401057106220540175
Milad MR,Quirk GJ, (. Year: 2002Neurons in medial prefrontal cortex signal memory for fear extinctionNature420707412422216
Müller MB,Zimmermann S,Sillaber I,Hagemeyer TP,Deussing JM,Timpl P,Kormann MS,Droste SK,Kühn R,Reul JM,Holsboer F,Wurst W, (. Year: 2003Limbic corticotropin-releasing hormone receptor 1 mediates anxiety-related behavior and hormonal adaptation to stressNature Neuroscience611001107
Nelson EC,Agrawal A,Pergadia ML,Wang JC,Whitfield JB,Saccone FS,Kern J,Grant JD,Schrage AJ,Rice JP,Montgomery GW,Heath AC,Goate AM,Martin NG,Madden PA, (. Year: 2010H2 haplotype at chromosome 17q21.31 protects against childhood sexual abuse-associated risk for alcohol consumption and dependenceAddiction Biology1511119878140
Oitzl MS,Fluttert M,Sutanto W,de Kloet ER, (. Year: 1998Continuous blockade of brain glucocorticoid receptors facilitates spatial learning and memory in ratsEuropean Journal of Neuroscience10375937669875354
Otagiri A,Wakabayashi I,Shibasaki T, (. Year: 2000Selective corticotropin-releasing factor type 1 receptor antagonist blocks conditioned fear-induced release of noradrenaline in the hypothalamic paraventricular nucleus of ratsJournal of Neuroendocrinology121022102611012844
Peters J,Kalivas PW,Quirk GJ, (. Year: 2009Extinction circuits for fear and addiction overlap in prefrontal cortexLearning and Memory1627928819380710
Pitman RK,Gilbertson MW,Gurvits TV,May FS,Lasko NB,Metzger LJ,Shenton ME,Yehuda R,Orr SP,Harvard/VA PTSD Twin Study Investigators (. Year: 2006Clarifying the origin of biological abnormalities in PTSD through the study of identical twins discordant for combat exposureAnnals of the New York Academy of Sciences107124225416891575
Prokasy WF,Ebel HC, (. Year: 1967GSR conditioning and sensitization as a function of intertrial intervalJournal of Experimental Psychology73247256
Radloff LS, (. Year: 1977The CES-D Scale: a self report depression scale for research in the general populationApplied Psychological Measurement1385401
Radulovic J,Rühmann A,Liepold T,Spiess J, (. Year: 1999Modulation of learning and anxiety by corticotropin-releasing factor (CRF) and stress: differential roles of CRF receptors 1 and 2Journal of Neuroscience195016502510366634
Rauch SL,Shin LM,Phelps EA, (. Year: 2006Neurocircuitry models of posttraumatic stress disorder and extinction: human neuroimaging research – past, present, and futureBiological Psychiatry1537638216919525
Ressler KJ,Mercer KB,Bradley B,Jovanovic T,Mahan A,Kerley K,Norrholm SD,Kilaru V,Smith AK,Myers AJ,Ramirez M,Engel A,Hammack SE,Toufexis D,Braas KM,Binder EB,May V, (. Year: 2011Post-traumatic stress disorder is associated with PACAP and the PAC1 receptorNature47049249721350482
Rodrigues SM,LeDoux JE,Sapolsky RM, (. Year: 2009The influence of stress hormones on fear circuitryAnnual Review of Neuroscience32289313
Rohleder N,Joksimovic L,Wolf JM,Kirschbaum C, (. Year: 2004Hypocortisolism and increased glucocorticoid sensitivity of pro-inflammatory cytokine production in Bosnian war refugees with posttraumatic stress disorderBiological Psychiatry5574575115039004
Roozendaal B,McReynolds JR,Van der Zee EA,Lee S,McGaugh JL,McIntyre CK, (. Year: 2009Glucocorticoid effects on memory consolidation depend on functional interactions between the medial prefrontal cortex and basolateral amygdalaJournal of Neuroscience29142991430819906977
Rosmond R,Chagnon YC,Chagnon M,Pérusse L,Bouchard C,Björntorp P, (. Year: 2000A polymorphism of the 5'-flanking region of the glucocorticoid receptor gene locus is associated with basal cortisol secretion in menMetabolism491197119911016903
Shin LM,Liberzon I, (. Year: 2010The neurocircuitry of fear, stress, and anxiety disordersNeuropsychopharmacology3516919119625997
Spielberger CD,Gorsuch RL,Lushene RE, (. Year: 1970Manual for the State-Trait Anxiety InventoryConsulting Psychologists PressPalo Alto
Stevens A,Ray DW,Zeggini E,John S,Richards HL,Griffiths CE,Donn R, (. Year: 2004Glucocorticoid sensitivity is determined by a specific glucocorticoid receptor haplotypeJournal of Clinical Endocrinology and Metabolism8989289714764810
Timpl P,Spanagel R,Sillaber I,Kresse A,Reul JM,Stalla GK,Blanquet V,Steckler T,Holsboer F,Wurst W, (. Year: 1998Impaired stress response and reduced anxiety in mice lacking a functional corticotropin-releasing hormone receptor 1Nature Genetics191621669620773
Treutlein J,Kissling C,Frank J,Wiemann S,Dong L,Depner M,Saam C,Lascorz J,Soyka M,Preuss UW,Rujescu D,Skowronek MH,Rietschel M,Spanagel R,Heinz A,Laucht M,Mann K,Schumann G, (. Year: 2006Genetic association of the human corticotropin releasing hormone receptor 1 (CRHR1) with binge drinking and alcohol intake patterns in two independent samplesMolecular Psychiatry1159460216550213
Tronel S,Alberini CM, (. Year: 2007Persistent disruption of a traumatic memory by postretrieval inactivation of glucocorticoid receptors in the amygdalaBiological Psychiatry62333917207472
van Rossum EF,Koper JW,van den Beld AW,Uitterlinden AG,Arp P,Ester W,Janssen JA,Brinkmann AO,de Jong FH,Grobbee DE,Pols HA,Lamberts SW, (. Year: 2003Identification of the BclI polymorphism in the glucocorticoid receptor gene: association with sensitivity to glucocorticoids in vivo and body mass indexClinical Endocrinology (Oxford)5958559214616881
Wasserman D,Wasserman J,Rozanov V,Sokolowski M, (. Year: 2009Depression in suicidal males: genetic risk variants in the CRHR1 geneGenes Brain and Behavior87279
Wessa M,Rohleder N, (. Year: 2007Endocrine and inflammatory alterations in posttraumatic stress disorderExpert Review of Endocrinology and Metabolism291
Wigginton JE,Cutler DJ,Abecasis GR, (. Year: 2005A note on exact tests of Hardy–Weinberg equilibriumAmerican Journal of Human Genetics7688789315789306
Yang YL,Chao PK,Lu KT, (. Year: 2006Systemic and intra-amygdala administration of glucocorticoid agonist and antagonist modulate extinction of conditioned fearNeuropsychopharmacology3191292416205786
Yehuda R, (. Year: 2006Advances in understanding neuroendocrine alterations in PTSD and their therapeutic implicationsAnnals of the New York Academy of Sciences107113716616891568
Yehuda R,Cai G,Golier JA,Sarapas C,Galea S,Ising M,Rein T,Schmeidler J,Müller-Myhsok B,Holsboer F,Buxbaum JD, (. Year: 2009Gene expression patterns associated with posttraumatic stress disorder following exposure to the World Trade Center attacksBiological Psychiatry6670871119393990
Yehuda R,Flory JD,Pratchett LC,Buxbaum J,Ising M,Holsboer F, (. Year: 2010Putative biological mechanisms for the association between early life adversity and the subsequent development of PTSDPsychopharmacology21240541720706708
Yehuda R,Golier JA,Yang RK,Tischler L, (. Year: 2004Enhanced sensitivity to glucocorticoids in peripheral mononuclear leukocytes in posttraumatic stress disorderBiological Psychiatry551110111615158431
Yehuda R,LeDoux J, (. Year: 2007Response variation following trauma: a translational neuroscience approach to understanding PTSDNeuron56193217920012


[Figure ID: fig001]
Fig. 1 

Fear-conditioning paradigm. CS+, Neutral conditioned stimulus (yellow rhombus) later followed by an aversive unconditioned stimulus (US), an electric shock (red flash); ITI, inter-trial interval; CS−, neutral stimulus (blue square) later followed by the absence of the US.

[Figure ID: fig002]
Fig. 2 

(a) Increased activity in the left amygdala elicited by CS+ versus CS- in the first half of the acquisition as a function of NR3C1 genotype, coded 0 for no minor allele ( ), 1 for one or two minor alleles (□), 2 for more than two minor alleles ( ) (for sample 1, group 0: n=26, group 1: n=26, group 2: n=8; for sample 2, group 0: n=16, group 1: n=30, group 2: n=6). CS+, neutral conditioned stimulus later followed by the aversive unconditioned stimulus (US); CS-, neutral stimulus later followed by the absence of the US. Values are means, with 95% confidence intervals (CIs) represented by vertical bars. * p<0.05. (b) Genotype-dependent differential activation of the prefrontal cortex during extinction involving CRHR1 and NR3C1 genotypes (coded 0 for no minor allele ( ), 1 for one minor allele (□), 2 for more than one minor allele ( ) (for sample 1, group 0: n=12, group 1: n=24, group 2: n=23; for sample 2, group 0: n=13, group 1: n=20, group 2: n=19). BA, Brodmann area. Values are means, with 95% CIs represented by vertical bars. * p<0.05. (c) T maps revealing increases in functional coupling for the contrasts between genotype group 2 versus groups 0 and 1 during the early acquisition phase (left panel) and genotype-dependent functional coupling during early acquisition between the left amygdala and prefrontal cortex (right panel). Group 0, no minor allele ( ); group 1, one minor allele (□); group 2, more than one minor allele ( ); AU, arbitrary units at the target-region peak voxels. Values are means, with 95% CIs represented by vertical bars. (d) t Maps revealing increases in functional coupling for the contrasts between the genotype groups during the extinction phase (left panel) and coupling strength for the extinction phase between the left prefrontal cortex and left amygdala (right panel). Group 0, no minor allele ( ); group 1, one minor allele (□); group 2, more than one minor allele ( ). Values are means, with 95% CIs represented by vertical bars.

[TableWrap ID: tab001] Table 1 

Clinical characteristics of the participants

tfn001_001Data are given as mean (standard deviation).

[TableWrap ID: tab002] Table 2 

Allele frequencies and genotype counts for the polymorphisms that contributed significantly to the association signals in the two samples

tfn002_000HWE, Hardy–Weinberg equilibrium.

aTest for deviation from HWE (Wigginton et al.2005).

[TableWrap ID: tab003] Table 3 

Peak voxel values of functional coupling analysis in sample 1 showing increased coupling for the genotype group coded 2 as compared with genotype groups 0 and 1a

tfn003_000BA, Brodmann area; MNI, Montreal Neurological Institute.

aCoordinates are stated in MNI space.

Article Categories:
  • Original Articles

Keywords: Fear conditioning, functional magnetic resonance imaging, genes, hypothalamic–pituitary–adrenal axis, PTSD.

Previous Document:  The inferior mesenteric vein to the left gonadal vein shunt for gastroesophageal varices and extrahe...
Next Document:  Infant mortality in Brazil, 1980-2000: a spatial panel data analysis.