Abstract
Catechol-O-methyltransferase (COMT) degrades the catecholamine neurotransmitters dopamine, epinephrine, and norepinephrine. A functional polymorphism in the COMT gene (val158met) accounts for a fourfold variation in enzyme activity. The low-activity met158 allele has been associated with improved working memory but with higher risk for anxiety-related behaviors. Using functional magnetic resonance imaging, we assessed the effects of COMT genotype on brain activation by standardized affective visual stimuli (unpleasant, pleasant, and neutral) in 35 healthy subjects. The analysis of genotype effects was restricted to brain areas with robust activation by the task. To determine genedose effects, the number of met158 alleles (0, 1, or 2) was correlated with the blood oxygen level-dependent (BOLD) response elicited by pleasant or unpleasant stimuli compared with neutral stimuli. COMT genotype had no significant impact on brain activation by pleasant stimuli but was related to the neural response to unpleasant stimuli: reactivity to unpleasant stimuli was significantly positively correlated with the number of met158 alleles in the limbic system (left hippocampus, right amygdala, right thalamus), connected prefrontal areas (bilateral ventrolateral prefrontal cortex, right dorsolateral prefrontal cortex), and the visuospatial attention system (bilateral fusiform gyrus, left inferior parietal lobule). Genotype explained up to 38% of interindividual variance in BOLD response elicited by unpleasant stimuli. We conclude that (1) genetic variations can account for a substantial part of interindividual variance in task-related brain activation and that (2) increased limbic and prefrontal activation elicited by unpleasant stimuli in subjects with more met158 alleles might contribute to the observed lower emotional resilience against negative mood states.
Introduction
Catechol-O-methyltransferase (COMT) degrades the catecholamine neurotransmitters dopamine (DA), epinephrine (EPI), and norepinephrine (NE). The COMT gene contains a G to A missense variant (Lachman et al., 1996) that translates into a substitution of methionine for valine at codon 158 (val158met). The enzyme containing met158 is unstable at 37°C and has one-third to one-fourth of the activity of the val158 enzyme (Spielman and Weinshilboum, 1981; Lotta et al., 1995). The alleles are codominant, because heterozygous individuals have enzyme activity that is midway between homozygous individuals (Weinshilboum et al., 1999). Despite this rather strong effect of COMT genotype on neurotransmitter metabolism, effects on behavior are moderate: the val158 allele was associated with impaired working memory (Goldberg et al., 2003) and increased perseverative errors in an executive function task (Egan et al., 2001), but only 4% of variance was attributable to COMT genotype. Intermediate phenotypes such as task-dependent brain activation may more sensitively measure gene effects on the brain, and, indeed, functional magnetic resonance imaging (fMRI) revealed a greater efficiency in the prefrontal cortex (PFC) of carriers of the met158 allele (Egan et al., 2001). Egan et al. (2001) suggested that the COMT val158 allele is associated with faster prefrontal DA metabolism, lower DA concentrations, and a reduced PFC neuronal signal/noise ratio.
Given its disadvantages for PFC function, it is curious that the val158 allele is highly prevalent (∼50%) in various human populations (Palmatier et al., 1999). Recent findings indicate that the compensating advantage of the val158 allele may be an increase in emotional resilience against anxiety and dysphoric mood (Enoch et al., 2003a). Zubieta et al. (2003) demonstrated that individuals homozygous for the met158 allele of COMT displayed higher sensory and affective ratings of pain and a more negative affective response to sustained pain. The homozygous met158 genotype was associated with higher levels of dimensionally measured anxiety among women in two populations (Enoch et al., 2003b). Moreover, the met158 allele has been associated with several mental disorders characterized by negative mood states such as panic disorder (Woo et al., 2004), alcoholism (Tiihonen et al., 1999; Wang et al., 2001) as well as increased alcohol intake in social drinkers (Kauhanen et al., 2000), major depression (Ohara et al., 1998a), bipolar affective disorder (Mynett-Johnson et al., 1998; Papolos et al., 1998), and obsessive-compulsive disorder in males (Karayiorgou et al., 1999). However, Ohara et al. (1998b) found no significant association of COMT genotype and anxiety disorder, and Domschke et al. (2004) reported an association of panic disorder with the val158 allele in women.
If genotype effects on brain activation are stronger than on overt behavior, functional brain imaging should be an excellent tool to assess genotype effects on central processing of emotional stimuli (Hariri and Weinberger, 2003). Therefore, we examined the effects of the COMT val158met genotype on the brain response to standardized affective visual stimuli (Lang et al., 1988) in healthy volunteers using fMRI. We hypothesized that COMT activity determined by the val158met genotype would be correlated with brain activity elicited by emotional stimuli.
Materials and Methods
Subjects. Thirty-five right-handed healthy volunteers [9 women and 26 men; age, 40.6 ± 7.6 (mean ± SD) years] participated in the study after providing informed, written consent according to the Declaration of Helsinki. The Ethics Committee of the University of Heidelberg approved the study. All subjects were Caucasians of central European descent. Standardized clinical assessment with the Structured Clinical Interview I (First et al., 2001) was performed to exclude axes I and II psychiatric disorders according to Diagnostic and Statistical Manual of Mental Disorders, Revision IV and International Classification of Diseases of the World Health Organization, Revision 10. Pregnancy and present drug abuse were excluded with urine tests. Only subjects free of any medication were included. Additionally, levels of anxiety [Symptom Checklist 90 -Revised (SCL-90-R) (Derogatis, 1983); State-Trait-Anxiety Inventory (STAI) (Spielberger et al., 1970)] and depression [Center for Epidemiological Studies Depression Scale (CES-D) (Radloff, 1997); Hamilton Rating Scale for Depression (HAMD) (Hamilton, 1986)] as well as personality [Temperament and Character Inventory (TCI) (Cloninger et al., 1994)] and education were assessed (Table 1).
Imaging study. For emotion induction, we used affectively unpleasant, pleasant, and neutral pictures. Each category was represented by 18 pictures. Pleasant and unpleasant cues were taken from the International Affective Picture System (IAPS) (Lang et al., 1988), in which images are standardized across the dimensions of emotion, arousal, and valence (Lang, 1995). Participants were instructed to passively view the stimuli, because even simple rating tasks can alter brain activation pattern (Taylor et al., 2003).
The stimuli were presented for 750 ms using an event-related design and were arranged in an individually randomized order for each subject. To reconstruct the blood oxygen level-dependent (BOLD) event-related time course, it is necessary to sample data points at different peristimulus time points. This was achieved by a random jitter between intertrial interval and acquisition time, resulting in an equal distribution of data points after each single stimulus. The intertrial interval was randomized between three and six acquisition times (i.e., 9.9-19.8 s). During the intertrial interval, a fixation cross was presented.
Scanning was performed with a 1.5 T whole-body tomograph (Magnetom VISION; Siemens, Erlangen, Germany) equipped with a standard quadrature head coil. For fMRI, 24 slices were acquired every 3.3 s (4 mm thickness, 1 mm gap) using an EPI-Sequence [repetition time (TR), 1.8 ms; echo time (TE), 66 ms; α = 90°] with in-plane resolution of 64 × 64 pixels [field of view (FOV), 220 mm]. fMRI slices were oriented axially parallel to the anterior commissure-posterior commisure. A morphological three-dimensional T1-weighted magnetization prepared rapid gradient echo image data set (1 × 1 × 1 mm3 voxel size; FOV, 256 mm; 162 slices; TR, 11.4 ms; TE, 4.4 ms; α = 12°) covering the whole head was acquired for anatomical reference.
After the MR scan, a subset of the stimuli (8 of 18 per category) was again presented on a computer monitor for 6 s and assessed for arousal and valence according to the standardized procedure described by Lang et al. (1988).
Data analysis. Data were analyzed with Statistical Parametric Mapping (SPM2) (Welcome Department of Neurology, University College London, London, UK). The structural three-dimensional data set was coregistered to the first T2* image. The structural image was spatially normalized to a standard template using a 12-parameter affine transformation with additional nonlinear components. A nonlinear transformation was subsequently applied to the T2* data, and voxels were resampled at a resolution of 3 × 3 × 3 mm3. The functional data were smoothed using an isotropic Gaussian kernel for group analysis (12 mm full-width at half-maximum).
Statistical analysis was performed by modeling the different conditions (pleasant, unpleasant, and neutral pictures; δ functions convolved with a synthetic hemodynamic response function and its time derivative) as explanatory variables within the context of the general linear model on a voxel-by-voxel basis with SPM2. To detect association between COMT genotype and fMRI activation on a voxel-by-voxel basis, the contrast images of all subjects (signal change of pleasant vs neutral and unpleasant vs neutral pictures) were included in a second level regression analysis with SPM2. To model the assumed genedose effect, COMT genotype was coded as a covariate by the number of met158 alleles (0, 1, or 2). To control for potential confounding of genotype effects by gender, we performed an additional SPM analysis using a multiple regression with genotype and sex as independent variables.
The analysis of the genotype × task interaction was restricted to areas showing robust activation related to the task (emotional vs neutral stimuli) as defined by an F-contrast of p < 0.01 and a cluster size of at least 10 voxels.
A statistical threshold of p < 0.05 corrected for the volume of the entire mask was used for the regression analyses. Correction for multiple testing was performed with the False Discovery Rate (FDR) procedure. The FDR is the proportion of suprathreshold voxels that are false positives. In comparison, a traditional multiple comparisons procedure (e.g., Bonferroni or random field correction) controls Family-Wise Error Rate (FWER) at or below α. FWER is the chance of one or more false positives anywhere (not just among suprathreshold voxels). If there is truly no signal in the image anywhere, then an FDR procedure controls FWER, just as Bonferroni and random field methods do. If there is some signal in the image, an FDR method will be more powerful than the traditional methods (Genovese et al., 2002).
Genotyping. The genotype information and flanking sequences of the COMT polymorphism were collected using public single-nucleotide polymorphism (SNP) databases and analyzed using the following procedure: DNA was isolated from venous blood with the QIAamp DNA Mini kit (Qiagen, Hilden, Germany). PCR was performed with HotStarTaq-DNA Polymerase (Qiagen), 4 ng of template DNA in a total volume of 25 μl of PCR, and the primers COMTfor (5′-CACCTGTGCTCACCTCTCCT-3′) and COMTrev (5′-GGGTTTTCAGTGAACGTGGT-3′). PCR conditions were as follows: 95°C for 15 min initial heating, then 35 cycles of 94°C for 60 s, 58°C for 60 s, 72°C for 60 s, and final extension for 10 min at 72°C. After purification of the amplified products (QIAquick PCR Purification kit; Qiagen), the oligonucleotide primer 5′-CGGATGGTGGATTTCGCTGGC-3′ was extended according to the manufacturer's instructions (AcycloPrime-FP SNP detection kit; Perkin-Elmer Life Sciences, Emeryville, CA), and fluorescence polarization was measured using a Wallac Victor2 1420 Multilabel Counter (Perkin-Elmer). We assessed the genotype of the promoter of the serotonin transporter gene (5-HTT) with PCR using oligonucleotide primers (stpr5, 5′-GGCGTTGCCGCTCTGAATGC; int1, 5′-CAGGGGAGATCCTGGGAGGA). PCR amplification was performed in a final volume of 30 μl consisting of 50 ng of genomic DNA, 2.5 mm deoxyribonucleotides (dGTP/7-deaza-2′-dGTP; 1:1), 0.1 μg of sense and antisense primers, 10 mm tris-HCl, pH 8.3, 50 mm KCl, 1.5 mm MgCl2, and 1 U of TaqDNA polymerase. Annealing was performed at 61°C for 30 s, extension was performed at 72°C for 1 min, and denaturation was performed at 95°C for 30 s for 35 cycles. The resulting PCR product was electrophoretically separated on a 1.5% agarose gel.
Results
Genotyping
Of the 35 healthy subjects, 9 (3 females and 6 males) were homozygous for the met158 allele and 10 (2 females and 8 males) were homozygous for the val158 allele of the COMT gene; 16 subjects (4 females and 12 males) were heterozygous. Genotype and sex were not associated (χ2 = 0.499; df = 2; p = 0.799). We also assessed the promoter region of the serotonin transporter gene (5-HTTLPR), which has recently been associated with amygdala activation elicited by angry or fearful faces (Hariri et al., 2002a). No association between the two genotypes analyzed in this study was found (χ2 = 2.3; df = 4; p = 0.68) (for details, see supplemental Table S1, available at www.jneurosci.org as supplemental material).
Behavioral data
Data concerning anxiety, depression, and education are shown in Table 1. All scores of the rating scales were in the normal range and not significantly associated with COMT genotype. Subjective ratings of valence and arousal induced by the visual stimuli (unpleasant, neutral, and pleasant) showed that subjects experienced the positive pictures as arousing and affectively pleasant, whereas they rated the negative pictures as equally arousing (p = 0.19) but unpleasant (Table 2). Valence and arousal were not significantly correlated with COMT genotype (p > 0.18).
Main effect of the task regardless of genotype
Significant differences in brain activity elicited by emotional stimuli (pleasant and unpleasant) compared with neutral stimuli were found in a total of 4247 voxels (3 × 3 × 3 mm). The main effect of emotional stimuli comprised a distributed neuronal network of secondary and tertiary visual regions in the occipital and parietal lobe, bilateral amygdala, hippocampus and parahippocampus, dorsal cingulate, and regions in the left and right medial, ventrolateral, and dorsolateral PFC (compare Fig. 1). The following analyses of the genotype × task interactions were restricted to these brain regions [volume of interest (VOI)].
COMT val158met genotype and processing of pleasant visual stimuli
No significant positive (p > 0.2; corrected) or negative (p = 1; corrected) correlations between COMT genotype and brain activity elicited by pleasant stimuli were found inside the predefined VOI. Also, when gender was controlled for, no significant association between COMT genotype and brain activity was found.
COMT val158met genotype and processing of unpleasant visual stimuli
For unpleasant visual stimuli, regression analyses revealed a significant (p < 0.05 corrected for entire VOI) positive correlation between the number of met158 alleles and the BOLD fMRI response in limbic and prefrontal brain areas as well as in brain areas related to visuospatial attention (Table 3, Fig. 2). Brain activation in 573 of 4247 voxels comprised by the VOI (13%) was significantly affected by genotype. The highest correlations were found in the ventrolateral PFC [Brodmann area (BA) 47]: 38% of the interindividual variance of the BOLD response in the left ventrolateral PFC and 31% in the right ventrolateral PFC were explained by genotype. Besides being strong, the effect of genotype in the ventrolateral PFC was rather widespread. The clusters comprised 258/141 (left/right) voxels. Moreover, a significant correlation was found in the right dorsolateral PFC. The number of met158 alleles was also positively correlated with brain activation elicited by unpleasant cues in the left hippocampus, the right amygdala, and in the right thalamus. In these areas, 25-33% of interindividual variance was explained by genotype. Above all, brain activation elicited by unpleasant stimuli in the left and right fusiform gyrus and the inferior parietal lobule, which are secondary and tertiary visual areas, was also positively correlated with the dosage of met158 alleles. When gender was controlled for, the relationship of COMT genotype and the processing of unpleasant stimuli remained statistically significant in all brain regions described above (for details, see supplemental Table S2, available at www.jneurosci.org as supplemental material), indicating that the observed effects of COMT on brain activation are not considerably modulated by gender.
Notably, in none of the investigated brain areas, a significant negative correlation between the number of met158 alleles and BOLD response was found. This was also true when gender was controlled for (p = 1; corrected).
Discussion
Val158met, a functional polymorphism in the COMT gene, interacted substantially with the neural processing of unpleasant stimuli in the limbic system, connected PFC, and visuospatial attention system: up to 38% of the interindividual variance of the fMRI BOLD response to unpleasant emotional stimuli was explained by a gene-dose effect of the COMT val158met genotype. The magnitude of the BOLD response to unpleasant stimuli always increased with the number of met158 alleles; in none of these brain areas, a negative correlation was detected.
The COMT variant containing met158 has 25-70% of the activity of the val158 enzyme (Spielman and Weinshilboum, 1981; Chen et al., 2004). COMT is expressed in all regions of the human brain, including the limbic and prefrontal areas (Hong et al., 1998). Therefore, the positive correlation between the number of met158 alleles and the BOLD response suggests that the amplified activation of limbic and prefrontal areas elicited by unpleasant stimuli may be related to increased levels of DA and/or NE in these brain regions.
The impact of COMT genotype on the processing of unpleasant stimuli was most pronounced in the ventrolateral PFC. Studies in COMT-deficient mice and studies in rats with the COMT inhibitor tolcapone showed increased DA but unaltered NE levels in the PFC (Gogos et al., 1998; Tunbridge et al., 2004). Therefore, the modulating effect of COMT genotype on PFC function might primarily be attributable to differences in DA levels. Activation of the ventrolateral PFC has been reported during induction of sadness (Levesque et al., 2003a,b) and recall of sad autobiographical memories (Markowitsch et al., 2003; Pelletier et al., 2003). Increased activation in the ventrolateral PFC was also found in patients suffering from major depression (Drevets, 2000; Brody et al., 2001). In patients with obsessive-compulsive disorder, response to treatment was associated with metabolic decreases in the ventrolateral PFC (Saxena et al., 2002), and provocation of obsessive and compulsive symptoms in healthy controls activated the ventrolateral PFC (Mataix-Cols et al., 2003). These findings indicate that the ventrolateral PFC participates in the processing of negative emotions. This prefrontal region receives extensive sensory inputs as well as limbic afferences from the amygdala. Therefore, the observed activation of the ventrolateral PFC may be associated with the integration of sensory inputs and information regarding changes in the subjects' emotional state.
In this study, activation of the right amygdala by unpleasant stimuli also increased with the number of met158 alleles and may indicate identification of stimuli with emotional significance (Phan et al., 2002; Phillips et al., 2003). In rats, DA potentiates the response of the amygdala by attenuating the effect of inhibitory input from the PFC and augmenting the effect of excitatory input from sensory cortices (Rosenkranz and Grace, 1999). Consistent with this, two fMRI studies in humans showed that DA modulates the responsivity of the amygdala to emotional stimuli (Hariri et al., 2002b; Tessitore et al., 2002). The potentiated amygdala response we report may thus reflect DA gating of amygdala inputs. Moreover, COMT genotype interacted with activation in the right thalamus and left dorsal hippocampus elicited by unpleasant stimuli. The thalamus is closely linked to limbic areas and is involved in controlling attentional processes. The level of its activation is modified by NE (Coull, 1998). The hippocampus has been implicated in the inhibition of stress responses (Lopez et al., 1999) and is crucial in the regulation of arousal and affective states (McNaughton and Gray, 2000). COMT genotype also seems to modulate thalamic endorphin release and associated pain perception (Zubieta et al., 2003). Together with our findings, these data indicate that in carriers of the met158 allele, limbic networks may be more sensitive to unpleasant stimuli and signaling of emotional salience might be enhanced.
Activation in the fusiform gyrus and parietal lobe elicited by unpleasant stimuli was also linked to a higher dosage of the met158 allele. Activation of these secondary and tertiary visual association areas may reflect increased attentional processing of emotionally salient visual stimuli and has been observed in a series of brain imaging studies using the IAPS (Lang et al., 1998). Positive feedback from the amygdala to the fusiform gyrus might contribute to the observed effect in this brain area (Morris et al., 1998; Vuilleumier et al., 2002).
Interestingly, in our study, the COMT val158met genotype was only associated with processing of unpleasant but not affectively pleasant stimuli. Because stimuli of both categories were rated as equally arousing, this specific effect cannot be explained on the basis of arousal differences. Another explanation might be that only stressful but not pleasant events activate the locus ceruleus, which results in increased NE release in its projection sites, including the amygdala, hippocampus, and the PFC (Tanaka et al., 2000; Charney, 2004). Therefore, in our study, presentation of aversive stimuli may have induced NE efflux in limbic structures and the connected PFC, and the resulting extracellular NE concentrations may then be affected by differences in NE metabolism attributable to COMT genotype.
Because COMT val158met genotype interacted with brain activity elicited by unpleasant but not pleasant stimuli, it seems to be primarily relevant for the processing and regulation of negative affective states. It can be speculated that in carriers of the met158 allele, heightened signaling of emotional salience and increased allocation of attentional resources leads to an increased demand on prefrontal capacities for integration and regulation of emotional states. Interestingly, ratings of emotional valence and arousal did not differ between genotypes. Thus, genetically determined differences in emotional processing detected with fMRI may be compensated at the behavioral level. Viewing unpleasant, aversive, or disgusting pictures inside the scanner may only be a minor challenge for emotional brain circuits, and highly demanding situations in real life should have a much stronger impact. In subjects with one met158 allele, and progressively more in carriers of two copies, the limits for integration and regulation of emotional states might be reached earlier, resulting in a reduced resilience against negative mood states (Enoch et al., 2003b; Zubieta et al., 2003).
In the study by Egan et al. (2001), the met158 allele was associated with a more focused dorsolateral prefrontal response during performance of a test of executive function. In our study, the met158 allele was associated with increased activation of the ventrolateral PFC during passive viewing of unpleasant stimuli, indicating that this allele is not generally associated with a low BOLD response in the PFC and that the associated changes in DA or NE concentrations may have differential effects on the performance of executive functions and the passive perception of affective stimuli in different parts of the PFC.
Hariri et al. (2002a) found that carriers of the short allele of the serotonin transporter gene (5-HTTLPR) showed a significantly enhanced amygdala response to angry and fearful faces, which probably reflects hyperresponsiveness of their amygdala to relevant environmental stimuli. In carriers of the short allele, slightly increased anxiety-related personality traits have been found (Lesch et al., 1996), and stressful life events more often lead to mood disorder and suicidality (Caspi et al., 2003). The genetically driven increase of amygdala reactivity may thus contribute to higher vulnerability to stressful life events. fMRI paradigms like the ones used by Hariri et al. (2002a) and in our study may be useful tools to define intermediate phenotypes that are related to reduced resilience against negative mood states and therefore to vulnerability for anxiety or mood disorders.
Studies investigating associations of a genetic polymorphism with an (intermediate) phenotype are critically dependent on the assumption that the genotype groups are not different in factors other than genotype. To exclude confounding attributable to the effect of an occult substructure, we tested for differences in ethnicity, age, gender, education, anxiety or depression scores, 5-HTTLPR genotype, as well as arousal and valence ratings of the affective stimuli in the COMT genotype groups. None of these potentially confounding variables was associated with COMT genotype. Moreover, in our study, stratification attributable to population differences in clinical ascertainment is extremely unlikely because genotype is correlating with a difference in brain function, not with diagnosis.
The substantial influence of genetic polymorphisms on emotional as well as cognitive processing found by us and others indicates that effects of genetic variants on individual brain function might obscure other differences between subjects or groups of subjects. Therefore, controlling for genetic heterogeneity may reduce noise in data and facilitate detection of differences related to other variables of interest.
In summary, our data show that the processing of unpleasant visual stimuli increased with met158 allele dosage. These COMT-driven differences in emotional processing are possibly attributable to altered DA and NE levels. The increased activation in limbic circuits elicited by unpleasant stimuli was accompanied by increased activity in prefrontal regions that are involved in the regulation of emotions. When confronted with unpleasant cues, the load limit of systems concerned with emotional and cognitive behavior control may be reached earlier in met158 allele carriers, contributing to their lowered resilience against negative mood states.
Footnotes
This work was supported by the Deutsche Forschungsgemeinschaft (He 2597/4-2) and the German Ministry of Education and Research (01GS0117).
Correspondence should be addressed to Dr. Andreas Heinz, Department of Psychiatry and Psychotherapy, Charité Campus Mitte, Charité-University Medicine Berlin, Schumannstrasse 20-21, 10117 Berlin, Germany. E-mail: andreas.heinz{at}charite.de.
Copyright © 2005 Society for Neuroscience 0270-6474/05/250836-07$15.00/0