INTRODUCTION

Emotional dysregulation is a cardinal feature of bipolar disorder (BD). The ventral anterior cingulate cortex (vACC) and the amygdala are highly interconnected structures that serve as essential nodes in emotional processing (Bush et al, 2000; Devinsky et al, 1995; Etkin et al, 2006; Maren and Quirk, 2004; Paus, 2001; Pezawas et al, 2005; Rosenkranz et al, 2003; Whalen et al, 1998a; Vogt et al, 1992). Neuroimaging studies of BD implicate aberrancies in vACC structure (Drevets et al, 1997; Lyoo et al, 2004; Sassi et al, 2004) and function (Blumberg et al, 2005; Drevets et al, 1997; Gruber et al, 2004; Kruger et al, 2002; Lawrence et al, 2004) including abnormalities in response to faces depicting positive, negative, or neutral emotional expressions (Blumberg et al, 2005; Lawrence et al, 2004). Heightened amygdala responses to faces have also been reported in BD (Altshuler et al, 2005; Blumberg et al, 2005; Chen et al, 2006; Ketter et al, 2001; Lawrence et al, 2004; Pavuluri et al, 2007; Rich et al, 2006; Yurgelun-Todd et al, 2000).

BD is a highly heritable disorder (Bertelsen et al, 1977; McGuffin et al, 2003; Smoller and Finn, 2003), and although several genetic loci have been linked to or associated with BD, many of these findings have not been consistently replicated (DePaulo, 2004). This suggests that BD may best be appreciated as a spectrum of illnesses consisting of several overlapping but possibly differentiable entities, with different sets of genes predisposing to particular phenotypes within this spectrum (Kelsoe, 2003). To understand this phenotypic heterogeneity in the context of a genetically complex background, researchers have begun to focus on endophenotypes within BD and on elucidating their underlying genes (Hasler et al, 2006).

The short‘s’ allele, as opposed to the long ‘l’ allele, of a functional polymorphism (5-HTTLPR) within the promoter region of the serotonin transporter protein gene (locus SLC6A4) has been associated with younger age of BD onset (Bellivier et al, 2002), differences in mood-state cycling (Cusin et al, 2001; Rousseva et al, 2003), psychotic symptoms (Ho et al, 2000; Ospina-Duque et al, 2000), and violent suicidal behavior (Bellivier et al, 2000). These findings in BD suggest that variation at this locus may be associated with phenotypic features of the disorder. In addition, Pezawas et al (2005), using structural and functional magnetic resonance imaging (fMRI), found that healthy individuals carrying the s allele had reduced gray matter volume in the ACC as well as impaired structural and functional connectivity from the vACC to the amygdala. Given that the vACC contains the richest concentration of serotonergic neurons of any cortical structure (Varnas et al, 2004), it is not surprising that the vACC is particularly influenced by 5-HTTLPR. In BD, abnormalities in 5-HTT binding in the vACC were reported in a positron emission tomography study (Cannon et al, 2006). Also, Hariri et al (2002) found that healthy s carriers had greater amygdala activation to negatively valenced stimuli compared to ll homozygotes, a finding replicated in other healthy samples (Bertolino et al, 2005; Canli et al, 2005; Heinz et al, 2005) and among individuals with mood and anxiety disorders (Dannlowski et al, 2007; Domschke et al, 2006; Furmark et al, 2004). Taken together, these findings suggest that 5-HTTLPR variation might contribute to abnormalities in vACC and amygdala response to emotional stimuli in BD.

The aim of this study was to investigate whether 5-HTTLPR variation influences vACC and amygdala dysfunction in BD. This was measured by blood oxygen level-dependent (BOLD) response to emotional face stimuli that were found previously to differentially activate this neural system in a BD group compared to a healthy comparison (HC) group (Blumberg et al, 2005). Euthymic BD individuals were studied to assess effects associated with the BD trait. It was hypothesized that in the vACC, BOLD response to face stimuli would be decreased in BD participants compared to HC individuals, and decreased amongst s allele carriers compared to ll homozygotes within each diagnostic group, such that the BD s carriers would have the greatest magnitude of dysfunction. It was similarly hypothesized that in the amygdala, BOLD response to faces would be increased in BD participants compared to HC individuals, and increased amongst s allele carriers compared to ll homozygotes within each diagnostic group. As differing effects within this vACC–amygdala neural system in association with specific facial emotions have previously been described (Blumberg et al, 2005; Lawrence et al, 2004), these effects were also explored in this study.

METHODS

Participants

The HC group included 48 participants, 58% of whom were women, with mean age and standard deviation of 28.2±10.0 years, whereas the BD group included 30 participants, 67% women, age of 31.9±13.5 years. A total of 65% of the HC group were European American (EA), 19% were African American, 6% were Asian American, and 10% were of other ethnicity, whereas 70% of the BD group were EA, 27% were African American, and one individual (3%) was of other ethnicity (see Table 1 for further details). BD participants were recruited from the Veterans Affairs Connecticut Healthcare System in West Haven, CT, other medical centers affiliated with the Yale University School of Medicine, and from surrounding communities. HC participants were recruited from surrounding communities. Exclusion criteria for both groups included: history of neurological disease, loss of consciousness for 5 min, any medical condition that may potentially affect neurovascular function including hypertension, or current substance abuse or dependence. Presence or absence of psychiatric illness was assessed using the Structured Clinical Interview for DSM-IV Axis I Disorders (SCID) Version 2.0 (First et al, 1995) which also was used to determine if BD participants were euthymic at scanning. In addition, scores on the Young Mania Rating Scale (YMRS (Young et al, 1978)) and the 29-item Hamilton Depression Rating Scale (HDRS (Williams et al, 1994)) were obtained from each participant except for one who did not complete these assessments. Exclusion criteria for the HC group included personal history of DSM-IV Axis I diagnosis or first-degree family member history of such illnesses as assessed by the Family History Screen for Epidemiologic Studies (FHE (Lish et al, 1995)). Written informed consent was obtained before participation in accordance and with approval of the human investigation committees of the Yale School of Medicine and the Veterans Affairs Connecticut Healthcare System.

Table 1 Features of the Bipolar Disorder and Healthy Comparison Groups

A total of 14 (47%) of the BD participants met criteria for rapid cycling and 9 (30%) had a history of psychosis, but none were psychotic at time of scan. Four (13%) were off medications at time of scan, all for greater than 1 year before scan, with the remaining participants taking psychotropic medications at time of scan, including lithium (33%, 10 of 30), anticonvulsants (60%, 18), atypical antipsychotics (43%, 13), antidepressants (50%, 15), and benzodiazepines (10%, 3). Ten (33%) BD participants had a prior history of substance related disorders. This included nine (30%) BD participants with a history of alcohol abuse or dependence, and eight (27%) with history of substance abuse or dependence; in all instances, individuals had been in remission for greater than 6 months at time of scan. The majority of BD subjects (80%) reported experiencing significant mood symptoms in childhood or adolescence, and all were symptomatic by young adulthood. A total of 22 (73%) of the BD participants had been hospitalized for psychiatric treatment at least once, and 6 (20%) had a history of suicide attempt.

Genotyping

Blood samples (10 ml) were drawn from participants for DNA extraction, which was accomplished by standard methods. PCR amplification was completed using an MJR tetrad cycler using methods and PCR primers previously described and designed by Gelernter et al (1997). As the s allele has been reported to be dominant in transcriptional activity vs the ll genotype (Heils et al, 1996), genotyping led to the creation of two genotype comparison groups, ll individuals and s carriers (ls or ss).

MRI Data Acquisition and Emotional Face Paradigm

Participants were scanned using a 3-Tesla Siemens Trio MR scanner (Siemens, Erlangen, Germany) at the Yale Magnetic Resonance Research Center in New Haven, CT. fMRI data were acquired with a single-shot echo planar imaging (EPI) sequence in alignment with the anterior commissure–posterior commissure (AC-PC) plane to create thirty-two 3-mm thick slices without gap with the following parameters: TR=2000 ms, TE=25 ms, matrix=64 × 64, FOV=240 × 240 mm2, and flip angle=80°.

During the functional runs, an event-related emotional face task was completed by each participant. Faces from the Ekman series (Ekman and Friesen, 1979) depicting expressions of fear, happiness, or neutrality were shown to the participants via the PsyScope software package (Cohen et al, 1993) on a computer attached to a projector. Participants were asked to make a male–female discrimination via a two-response button box that they were oriented to before scanning. This task was chosen to induce implicit rather than explicit processing of the emotional face as the former has been found to activate the emotional processing neural circuit more readily (Whalen et al, 1998b). Each face was presented for 2 s and separated by 4, 8, or 12-s intervals during which time a cross-hair fixation point was displayed. Face stimuli included images of 10 actors (five of each sex), with each individual exhibiting all three of the expressions for a total of 30 faces per run. Ordering of face stimuli was varied systematically to control for sequential dependencies and counterbalanced for facial expression, sex, the identity of the face, and the length of the interval between stimuli. Each run lasted 4 min and 50 s and data were compiled and averaged over four runs.

fMRI Data Processing

Raw data preprocessing was completed using SPM99 software (Statistical Parametric Mapping; http://www.fil.ion.ucl.ac.uk/spm). Two images at the beginning of each fMRI run were discarded to account for the approach of the hemodynamic response to steady state. The functional scans were realigned to the first volume to correct for interscan motion. The functional data were then spatially normalized to a standard EPI template from the Montreal Neurological Institute (MNI) with resampling to 4 mm3 voxels, and spatially smoothed with a 12 mm full width at half maximum (FWHM) Gaussian kernel.

SPM99 was also used for the model specification and estimation. At the individual subject level, event-related response amplitudes were estimated using the general linear model (Friston et al, 1995) for each of the three event types: fearful, happy, and neutral expression. This created statistical images of the BOLD signal change associated with each of the three emotional face types vs the baseline fixation cross-hair control for each individual subject.

Region of Interest Definition

The vACC region of interest (ROI) was defined by the WFU toolbox (Wake Forest University PickAtlas Region of Interest Toolbox; http://www.fmri.wfubmc.edu/download.htm) predefined anatomic label excluding the dorsal component (dorsal to the plane at z=0), creating a region that was 5496 mm3 (86 voxels) in volume with center at MNI coordinates of x=0 mm, y=36 mm, z=−4 mm. This division approximates the region of the ACC most associated with emotional processing (Bush et al, 2000; Devinsky et al, 1995; Vogt et al, 1992). Bilateral amygdala ROIs were similarly defined by the WFU toolbox, each 1920 mm3 in volume (30 voxels). The mean percent signal change was extracted from the ROI using the MarsBar ROI toolbox (http://marsbar.sourceforge.net) from the BOLD contrast images for each of the three emotional stimulus types for each participant.

Statistical Analysis

Statistical procedures were performed in SAS software, version 9.1 (SAS Institute Inc., Cary, NC) using a linear mixed-model analysis, in which the mean signal change in the vACC ROI for each subject was the dependent measure. Diagnosis and genotype were included as between-subjects factors, whereas the three emotion types were included as within-subject factors. A second mixed-model analysis was performed using the mean signal change in the amygdala ROIs with an additional left vs right hemisphere within-subject factor. Age and sex were entered as covariates in both initial models and were removed if they were not found to have significant effect. All two- and three-way interactions were tested. The correlation structure of the data was modeled by an unstructured variance–covariance matrix for observations on the three emotions. The latter variance–covariance structure was the best fitting according to the Akaike Information Criterion. Least square means (ls means) and standard errors (SE) were calculated in the mixed-model to allow interpretation of significant interactions. Planned post hoc analyses were performed within the main effect analysis for genotype and diagnosis to allow for evaluation of the potential differential effects of the varying emotional face types. SPM99 was used to afford visual localization of significant difference clusters from the above analyses.

Post hoc exploratory analyses were performed for potential main effects of clinical variables among BD participants on mean signal change values in the vACC and amygdala ROIs. Clinical factors examined included presence or absence of rapid cycling, history of psychosis, history of alcohol or substance abuse or dependence, history of psychiatric hospitalization, history of suicide attempt, and presence or absence of medication at time of scanning for specific medication subclasses.

A pertinent consideration in this study was the relative ethnic heterogeneity within the overall participant sample as there is evidence of variability in the frequency of the l and s alleles across ethnic groups (Gelernter et al, 1997), including EA and African-American populations. Thus, it was important to attempt to disentangle, as much as possible, the effect of a specific allele vs a specific ethnicity. To address this, t-tests were performed comparing activation in the vACC and amygdala between the EA and African-American participants within each gene-by-diagnosis subgroup, and secondary within-population EA only subgroup analyses were performed to verify that similar effects could be observed in a more homogeneous population (there were too few participants of the other population groups to allow for statistically meaningful analyses).

Whole Brain Analyses

To explore group differences not hypothesized a priori, SPM99 was used to perform whole brain voxel-based t-test comparisons between-genotype and between-diagnosis groups for each emotional face type. Findings were considered significant at a threshold of p<0.05, corrected for false discovery rate (FDR (Genovese et al, 2002)).

RESULTS

Participants

The HC and BD groups did not differ statistically in age (p=0.18) or sex distribution (p=0.39). Genotyping at the SLC6A4 locus led to four diagnosis–genotype groups: 14 HC ll participants, 34 HC s carriers (29 ls and 5 ss), 10 BD ll, and 20 BD s carriers (17 ls and 3 ss). The frequency of the l allele in this sample was 60.3%, and the frequency of the s allele was 39.7%. Clinical characteristics for these groups are summarized in Table 1. There was no significant difference in the distribution of EA vs non-EA participants across the groups (χ2=4.84, p=0.18). The BD group had significantly higher YMRS and HDRS scores than the HC group (p<0.001 for each). There was no significant difference in scores between the HC ll and s carriers or between the BD ll and s carriers.

The magnitude of movement during scanning was well within standards in the literature (Caparelli et al, 2003) and the groups did not have significantly different physical head motion in terms of translational (HC ll mean translation 0.62±standard deviation 0.27 mm, HC s car 0.56±0.25 mm, BD ll 0.61±0.20 mm, and BD s car 0.71±0.28 mm; F=1.57, p=0.20) or rotational movement (HC ll having 0.51±0.34°, HC s car 0.45±0.27°, BD ll 0.50±0.23°, and BD s car 0.55±0.44°; F=0.45, p=0.72).

Accuracy and button box reaction time averages were calculated for each participant in response to the different face types and these data are also shown in Table 1. There were no significant differences in mean reaction time or accuracy between the four diagnosis-by-genotype groups.

Ventral Anterior Cingulate Cortex

Age and sex were initially entered as covariates but were removed from the model as they did not have a significant effect on the results. There was a significant effect for genotype (F(1,74)=12.78, p=0.0006). The interaction between genotype and emotional face type was not significant (F(2,74)=1.89, p=0.16) and our planned post hoc analyses supported a genotype main effect that was present amongst ll individuals compared to s allele carriers across emotional face types: response to fearful (F(1,74)=15.2, p=0.0002), happy (F(1,74)=5.25, p=0.025), and neutral faces (F(1,74)=6.55, p=0.013). There was also a significant main effect of diagnosis (F(1,74)=6.58, p=0.012). Here results of our post hoc analyses for emotional face type showed that, although the interaction between diagnosis and emotional face type was not significant (F(2,74)=1.89, p=0.47), activation increases in the HC group compared to the BD group were significant in response to fearful (F(1,74)=7.85, p=0.007) and happy faces (F(1,74)=4.8, p=0.032) but not in response to neutral faces (F(1,74)=1.75, p=0.19). Figure 1 displays a graph depicting the ls mean signal change in each diagnosis-by-genotype group in response to the face stimuli in the vACC and Figure 2 shows BOLD contrast maps from SPM99 demonstrating the localization of clusters of difference within the vACC ROI for the significant comparisons from the mixed model. There were no significant interactions between genotype, diagnosis, or emotion type within the vACC mixed model.

Figure 1
figure 1

Blood oxygen level-dependent (BOLD) response in ventral anterior cingulate cortex (vACC) to faces. The graph displays the BOLD signal change to faces depicting fearful, happy, and neutral expressions (ls mean±standard error) by group (healthy comparison (HC) and bipolar disorder (BD)) and genotype (homozygous for the 5-HTTLPR ‘l’ allele (ll) and s carriers heterozygous and homozygous for the ‘s’ allele (s car)) for the vACC demonstrating the lowest signal in BD s carriers.

Figure 2
figure 2

Localization of blood oxygen level-dependent (BOLD) response differences in ventral anterior cingulate cortex (vACC) to faces depicting fearful expressions. The midsagittal images (x=0 mm) show the vACC region where activation to fear faces was decreased in (a) the bipolar disorder (BD) group compared to the healthy comparison (HC) group, (b) the HC s carrier group compared to the HC ll group, and (c) the BD s carriers compared to the BD ll group (all at p<0.01, uncorrected) displayed on an average T1 structural magnetic resonance imaging (MRI) template.

Post hoc analyses revealed that, among the BD participants, no significant effects on vACC response to any of the face types were detected in relation to clinical factors including presence of rapid cycling, history of psychosis, history of alcohol abuse/dependence, history of substance abuse/dependence, history of psychiatric hospitalization, history of suicide attempt, or medication subclass including lithium, anticonvulsants, antipsychotics, antidepressants, and benzodiazepines.

Amygdala

Age and sex were removed from the model as they did not have a significant effect on results. There was a significant main effect of hemisphere as left amygdala activation was higher compared to right amygdala activation to all face types (F(1,370)=10.2, p=0.0015). There was also a significant main effect of emotion as activation to fearful faces was greater than activation to both happy and neutral faces (F(2,370)=9.25, p=0.0001). Graphs depicting ls mean signal change in the left and right amygdala in response to face stimuli are shown in Figure 3. There were no significant main or interaction effects of genotype, diagnosis, emotion type, or hemispheric laterality within the amygdala mixed model.

Figure 3
figure 3

Blood oxygen level-dependent (BOLD) response in amygdala to faces. The graphs display the BOLD signal change to faces depicting fearful, happy, and neutral expressions (ls mean±standard error) by group (healthy comparison (HC) and bipolar disorder (BD)) and genotype (homozygous for the 5-HTTLPR ‘l’ allele (ll) and s carriers heterozygous and homozygous for the ‘s’ allele (s car)) for the (a) left amygdala and (b) right amygdala.

As in the vACC, amongst the BD participants, none of the clinical factors discussed above were found to have a significant effect on response in the amygdala to any of the emotional face types on either the left or right side.

Ethnicity Analyses

There were no statistically significant differences between EA and African-American participants in response in the vACC or amygdala ROIs within any of the diagnosis-by-genotype subgroups (all p-values >0.05). Consistent with the findings in the overall sample, in the EA only subgroup analyses for the vACC ROI, there was a significant main effect for genotype (F(1,48)=8.03, p=0.007), derived from an overall greater activation amongst ll individuals compared to s allele carriers in response to fear (F(1,48)=11.61, p=0.001) and neutral faces (F(1,48)=6.28, p=0.016). There was no significant main effect for diagnosis, but the directionality of difference was the same as above.

Whole Brain Analyses

In the exploratory whole brain diagnosis comparisons, in response to fearful faces, an additional region of higher activation in the HC group than the BD was detected in the right parahippocampal gyrus (BA 37). There were no other significant regions for which activation was significantly greater between diagnosis groups or genotype within each diagnostic group in response to fearful, happy, or neutral faces.

DISCUSSION

In this study, fMRI analyses demonstrated decreased response in the vACC among BD participants compared to the HC group in response to fearful faces and happy faces. This reduction in vACC response during emotion processing in BD is consistent with previous reports (Blumberg et al, 2005; Kruger et al, 2002). As this sample included euthymic-only BD participants, the findings suggest that deficits in vACC response to emotionally valenced face stimuli may be trait features of the disorder.

In the same vACC region, response was also decreased amongst carriers of the s allele of the 5-HTTLPR polymorphism in both the HC and BD groups compared to the ll individuals within the respective diagnostic groups, with the greatest magnitude of dysfunction seen in the BD s carriers. While previous studies in healthy individuals reported an association of the s allele with disruption of functional connectivity between the vACC and the amygdala (Pezawas et al, 2005), this study suggests it influences functioning within the vACC. This may be the result of a direct effect within the vACC. It is also possible that it is the result of modulation by other structures, such as the parahippocampal gyrus detected in exploratory analyses. In addition, whereas function in the vACC was found to be impaired in BD participants as a group overall, those carrying the s alleles exhibited the greatest magnitude of impairment in this region. Although the interaction of diagnosis and genotype was not significant, and the results of greater vACC dysfunction in BD participants who were also s carriers could represent independent and addictive effects, the subgroup of BD s carriers may represent a biologically salient subtype of BD that might most benefit from treatments targeting serotonergic transmission or treatments that target the vACC.

Although the vACC differences were consistent with our initial hypotheses, the lack of amygdala differences was not expected given reports of increased amygdala activation to emotional stimuli in BD participants (Altshuler et al, 2005; Blumberg et al, 2005; Chen et al, 2006; Drevets et al, 2002; Ketter et al, 2001; Lawrence et al, 2004; Pavuluri et al, 2007; Rich et al, 2006; Yurgelun-Todd et al, 2000), as well as in HC s carriers compared to HC ll individuals (Bertolino et al, 2005; Canli et al, 2005; Hariri et al, 2002; Heinz et al, 2005). A recent meta-analysis did note that there are unpublished data of samples in which those carrying the s allele did not have greater amygdala activation compared to ll individuals (Munafó et al, 2008). Many of the published reports in this meta-analysis involved subtraction of the response to a neutral condition from an emotional stimulus experimental condition (Munafó et al, 2008); in the current study the fixation crosshair was used as the baseline as amygdala response to neutral faces has been suggested to be heightened in BD and other disorders of emotional regulation (Donegan et al, 2003); significant amygdala activation to neutral face stimuli was also detected in the data presented herein. The level of amygdala signal change suggests that the emotional face task was effective in eliciting amygdala response across the entire sample, but perhaps so robustly as to disallow detection of significant group differences that might be detected with paradigms including faces depicting milder emotional expressions (Lawrence et al, 2004). Given the vACC differences, and the known functional and structural connectivity between the vACC and amygdala, the dissociation between effects in vACC and amygdala also raises the possibility that there is a defect, perhaps, in vACC–amygdala communication that is most prominent in the BD s carriers that will be of interest to pursue in future work.

The lack of a between-diagnosis difference in amygdala response may also have been due to the fact that the majority (86.7%) of the BD participants were on medications, which has been suggested to blunt amygdala activation in BD individuals compared to participants not taking medications (Blumberg et al, 2005). Although there was no significant difference in amygdala activation overall between the unmedicated and medicated BD participants perhaps due to the small number of unmedicated participants in this sample, the direction of the difference is consistent with this blunting, with greater activation in the unmedicated BD group compared to the medicated group (percent signal change across all face emotion types 0.41 vs 0.25, respectively), and three of these four participants having mean BOLD response in the upper third of the total distribution. Also, previous reports of increased amygdala BOLD activation in BD adults included individuals in noneuthymic mood states (Altshuler et al, 2005; Chen et al, 2006; Lawrence et al, 2004; Yurgelun-Todd et al, 2000), whereas the current study included only individuals who were euthymic at time of scan, raising the possibility that dysfunction in the amygdala is more pronounced during acute mood episodes of BD.

Analyses were performed to attempt to disentangle the potential effect of ethnicity given that there is known variability in allelic frequency amongst different ethnic groups (Gelernter et al, 1997; Lotrich et al, 2003). There was no significant difference in ethnic distribution between the groups and there were no significant differences found in vACC or amygdala response to the emotional face stimuli between the EA and African-American participants within each of the four diagnosis-by-genotype subgroup. In addition, in analyses for the EA only group, the direction of the effects of genotype and group in the vACC were the same as in the overall sample; however, the effects only reached significance for the effect of genotype for the fear faces, likely owing to the decreased power caused by the smaller number of EA only subjects. Sample sizes for the other ethnic groups did not permit meaningful analyses with those subjects only, and thus replication in larger, more homogeneous samples will be important.

We did not detect significant influences on neural circuitry functioning of clinical factors such as presence of rapid cycling, history of psychosis, history of substance or alcohol abuse/dependence, onset of mood symptoms in adolescence/young adulthood, history of psychiatric hospitalization, history of suicide attempt, or medication subclasses. However, it was unlikely that we had sufficient power to detect such effects. Continued study of larger samples will hopefully allow characterization of a ‘gene expression’ to ‘neural circuit dysfunction’ to ‘clinical phenotype’ mechanistic bridge, and further characterize the BD s carrier subgroup. As one potential example, we did note that the majority (88.9%) of the BD participants with a history of psychosis were s carriers, consistent with prior epidemiological studies associating the s allele with psychosis in BD (Ho et al, 2000; Ospina-Duque et al, 2000), suggesting that future studies of 5-HTTLPR variation, vACC neural system functioning and psychosis may further clarify this potential relationship. In addition, reports of interactions between stress and other genetic variants such as those in brain-derived neurotrophic growth factor with 5-HTTLPR variation (Kaufman et al, 2006; Pezawas et al, 2008) suggest that future work investigating these interactions may help to elucidate how variations in several genes and environmental factors may converge to lead to the development of a BD phenotype.

Whole brain analyses were conducted to explore for regions in the brain besides the vACC and amygdala for which BOLD response to the emotional face task differed across diagnosis and genotype groups. The only significant finding was a decrease in activation in the BD group compared to the HC group in response to fearful faces in the right parahippocampal gyrus. This area is considered part of the limbic association cortex, which also includes the vACC (Ebert and Ebmeier, 1996), and has been implicated in the processing of social and emotional stimuli (Rich et al, 2008). In addition, Lawrence et al (2004) found that parahippocampal gyrus response was decreased in adults with BD in response to emotional face stimuli, and Rich et al (2008) have reported similar functional deficits in pediatric BD participants.

In the future, studies of children and adolescents with BD, and family members at high risk for BD, may allow examination of groups less likely to be affected by medications and may reveal how the 5-HTTLPR influences neurodevelopmental changes over time. Positron emission tomography studies of healthy adults (Parsey et al, 2006) and postmortem histological studies of suicide victims have not shown differences in serotonin transporter distribution or concentration in ventral frontal regions amongst individuals carrying the s allele compared to ll individuals (Mann et al, 2000). These findings suggest that variation in 5-HTTLPR may result in vACC dysfunction in adulthood due to neurodevelopmental effects rather than acute effects on 5-HTT distribution. This is particularly relevant to BD, as ACC structural differences have not been consistently found in pediatric samples (Blumberg et al, 2006; Chang et al, 2005; Kaur et al, 2005; Sanches et al, 2005; Wilke et al, 2004), again suggesting that this dysfunction may progress over the course of adolescence as the vACC continues to mature (Blumberg, 2007).

In conclusion, the present work provides evidence for an association of the s allele of the 5-HTTLPR with deficits in the response of the vACC to emotional stimuli, a region that also demonstrated dysfunction in a sample of euthymic individuals with BD, with dysfunction greatest in those in the BD group carrying the s allele. In combination, these findings suggest that pronounced vACC dysfunction may be a heritable neurobiological trait, which underlies a distinct subtype within the more heterogeneous BD clinical phenotype. Evidence is emerging that this polymorphism may influence response to psychotropic medications, with s allele homozygotes having less favorable response to selective serotonin reuptake inhibitors (Smits et al, 2004) but a greater benefit with lithium compared to ll homozygotes (Stamm et al, 2007). Together, with continued work on the effects of neural circuitry, further investigation of treatment response related to 5-HTTLPR variation may in the future direct treatments to those most likely to benefit from serotonergic treatments and those that target the vACC (Mayberg et al, 2005).