Trait neuroticism is a relatively stable personality domain reflecting individual differences in threat and punishment sensitivity as well as the tendency toward negative affect, including emotions such as fear, anger, worry, frustration, sensitivity to criticism, hostility, vulnerability, self-consciousness, and frustration (DeYoung, Reference DeYoung2015; Widiger, Reference Widiger, Leary and Hoyle2009). Individuals high in neuroticism tend to interpret events as more threatening, react more negatively to events, and utilize avoidant and defensive coping strategies, such as anger, irritability, depression, panic, and anxiety, in the face of negative events (McCrae & Costa, Reference McCrae and Costa2003; Tackett & Lahey, Reference Tackett and Lahey2016). There is continuity between personality/temperament constructs assessed from childhood into adulthood (e.g., Caspi et al., Reference Caspi, Harrington, Milne, Amell, Theodore and Moffitt2003), and trait neuroticism in adulthood has its origins in related constructs that have been well documented in childhood (e.g., behavioral inhibition/trait anxiety; Caspi et al., Reference Caspi, Harrington, Milne, Amell, Theodore and Moffitt2003; Kagan, Reznick, Clarke, Snidman, & Garcia-Coll, Reference Kagan, Reznick, Clarke, Snidman and Garcia-Coll1984).
Trait neuroticism has been associated with increased rates of transdiagnostic psychopathology in childhood, adolescence, and adulthood (e.g., Ormel, Jeronimus, et al., Reference Ormel, Jeronimus, Kotov, Riese, Bos, Hankin and Oldehinkel2013; Ormel, Rosmalen, & Farmer, Reference Ormel, Rosmalen and Farmer2004; Tackett, Reference Tackett2006). In addition, this proneness toward negative emotionality and maladaptive behavior that is associated with neuroticism predicts a multitude of adverse functional outcomes; neuroticism is inversely correlated with marital satisfaction, occupational success, and quality of life (Lahey, Reference Lahey2009). As a result, trait neuroticism has an enormous economic cost on society, a cost that exceeds even that of common mental disorders (Cuijpers et al., Reference Cuijpers, Smit, Penninx, de Graaf, ten Have and Beekman2010). Given its important implications for both individuals and society, there has been considerable interest in investigations of how neuroticism relates to both psychological processes and biological systems (DeYoung et al., Reference DeYoung, Hirsh, Shane, Papademetris, Rajeevan and Gray2010; Ormel, Bastiaansen, et al., Reference Ormel, Bastiaansen, Riese, Bos, Servaas, Ellenbogen and Aleman2013; Patrick, Curtin, & Tellegen, Reference Patrick, Curtin and Tellegen2002). Substantial evidence supports the heritability of personality traits, including neuroticism (Lahey, Reference Lahey2009; Van Den Berg et al., Reference Van Den Berg, De Moor, McGue, Pettersson, Terracciano, Verweij and Boomsma2014; Widiger, Reference Widiger, Leary and Hoyle2009), suggesting the likely existence of underlying biological processes that explain individual differences in trait neuroticism.
Early Biological Studies on Trait Neuroticism
Early studies into the biological underpinnings of trait neuroticism were largely based on Eysenck's arousal theory of personality and Gray's reinforcement sensitivity theory of personality (Allen & Deyoung, Reference Allen and Deyoung2016; De Fruyt et al., Reference De Fruyt, De Clercq, De Bolle, Wille, Markon and Krueger2013; Eysenck, Reference Eysenck1967; Gray Reference Gray1982, Reference Gray and Madden1991). These theories suggest that neuroticism stems from the hyperarousal of certain key brain systems (e.g., limbic structures) with downstream effects on autonomic arousal and the hypothalamic–pituitary–adrenal axis. This system is responsible for regulating the body's stress reaction in response to threats or danger. Both Eysenck and Gray's early theories, as well as more recent conceptualizations of trait neuroticism (e.g., DeYoung & Gray, Reference DeYoung and Gray2009), posit a key role for altered functioning in the amygdala. The amygdala, an almond shaped mass of gray matter, is found in the medial temporal lobe of the brain, and is the brain region most closely associated with fear processing (Aggleton, Reference Aggleton1992). A wealth of research from the animal literature has provided the basis for understanding the association between the amygdala and its projections to other brain areas in response to threatening stimuli (Davis, Reference Davis1992). The amygdala connects to downstream brain regions involved in coordinating behavioral, neuroendocrine, and autonomic responses to emotional stimuli, including the hypothalamus and the brain stem (Depue, Reference Depue2009; Heimer, Reference Heimer2003; LeDoux, Reference LeDoux1998; Ormel, Bastiaansen, et al., Reference Ormel, Bastiaansen, Riese, Bos, Servaas, Ellenbogen and Aleman2013). The amygdala also has connections to higher cortical brain regions including the anterior cingulate cortex (ACC) and regions in the prefrontal cortex (PFC) involved in self-referential processing and the cognitive control of emotions (Kim et al., Reference Kim, Loucks, Palmer, Brown, Solomon, Marchante and Whalen2011; Ochsner & Gross, Reference Ochsner and Gross2005).
Early biological studies testing Eysenck and Gray's theories primarily utilized electrophysiological methods, such as skin conductance and salivary and urinary cortisol, which serve as global measures of central nervous system arousal and were understood to be proxies for amygdala activation (Cheng, Richards, & Helmstetter, Reference Cheng, Richards and Helmstetter2007; Urry et al., Reference Urry, van Reekum, Johnstone, Kalin, Thurow, Schaefer and Davidson2006). These studies found that individuals higher in neuroticism exhibited greater reactivity in a variety of psychophysiological markers, such as skin conductance (Norris, Larsen, & Cacioppo, Reference Norris, Larsen and Cacioppo2007), greater startle reactions to fearful stimuli (Wilson, Kumari, Gray, & Corr, Reference Wilson, Kumari, Gray and Corr2000), and greater salivary cortisol upon waking (Portella, Harmer, Flint, Cowen, & Goodwin, Reference Portella, Harmer, Flint, Cowen and Goodwin2005). Taken together, these studies suggest that there is a higher level of arousability and reactivity among individuals higher in neuroticism and suggest that neuroticism might be associated with hyperarousal in the upstream limbic system.
Magnitude of Brain Activation and Neuroticism
Basic emotion research with humans has provided evidence for the important role of the amygdala in emotional learning and memory, negative emotion processing, and threat appraisal (Britton, Lissek, Grillon, Norcross, & Pine, Reference Britton, Lissek, Grillon, Norcross and Pine2011; Davis & Whalen, Reference Davis and Whalen2001; Morris et al., Reference Morris, Friston, Büchel, Frith, Young, Calder and Dolan1998). In turn, a growing body of research has employed functional magnetic resonance imaging (fMRI) to identify functional neural correlates of trait neuroticism, with a focus on the amygdala, particularly during negative emotion processing (Ormel, Bastiaansen, et al., Reference Ormel, Bastiaansen, Riese, Bos, Servaas, Ellenbogen and Aleman2013; Servaas et al., 2013). A number of fMRI studies have found evidence of a positive correlation between trait neuroticism and amygdala activation, measured by quantifying the magnitude in task-related change in the blood oxygen level dependent (BOLD) signal. These studies employ diverse paradigms involving affective content, such as emotional scenes and faces, the emotional Stroop task, and an emotional prosody task (Brück, Kreifelts, Kaza, Lotze, & Wildgruber, Reference Brück, Kreifelts, Kaza, Lotze and Wildgruber2011; Chan, Norbury, Goodwin, & Harmer, Reference Chan, Norbury, Goodwin and Harmer2009; Cunningham, Arbuckle, Jahn, Mowrer, & Abduljalil, Reference Cunningham, Arbuckle, Jahn, Mowrer and Abduljalil2011; Haas, Omura, Constable, & Canli, Reference Haas, Omura, Constable and Canli2007; Harenski, Kim, & Hamann, Reference Harenski, Kim and Hamann2009). However, other studies using similar fMRI paradigms have failed to find such an association (Cremers et al., Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010; Drabant, McRae, Manuck, Hariri, & Gross, Reference Drabant, McRae, Manuck, Hariri and Gross2009; Haas, Constable, & Canli, Reference Haas, Constable and Canli2008; Hyde, Gorka, Manuck, & Hariri, Reference Hyde, Gorka, Manuck and Hariri2011; Thomas et al., Reference Thomas, Elliott, McKie, Arnone, Downey, Juhasz and Anderson2011). Most of these studies have been conducted in relatively small samples, which may account for these inconsistent findings (Yarkoni, Reference Yarkoni2009). However, a recent quantitative, parametric coordinate-based meta-analysis of fMRI and positron emission tomography studies investigating neural activity associated with neuroticism also failed to find a significant positive association between likelihood of amygdala activation and neuroticism across 15 studies with a total sample size of 485 participants (Servaas et al., Reference Servaas, van der Velde, Costafreda, Horton, Ormel, Riese and Aleman2013). These inconsistencies in the literature compelled us to ask whether there might be an association between neuroticism and other amygdala-based indicators, such as amygdala habituation or amygdala functional connectivity, which might support the hypothesized relationship between neuroticism and amygdala function.
Time Course of Activation and Neuroticism
Some recent research has shifted focus from examining the magnitude of brain activation in regions of interest, such as the amygdala, to studying the temporal dynamics of activation in brain regions as a possible neural marker of neuroticism. For example, during sustained processing of negative information, Haas et al. (Reference Haas, Constable and Canli2008) found neuroticism to be associated with ongoing activation in the medial PFC while responding to sad emotional faces, but not fearful or happy faces. The medial PFC has been associated with higher order cognitive control of emotional reactions. One study by Schuyler et al. (Reference Schuyler, Kral, Jacquart, Burghy, Weng, Perlman and Davidson2014), in a relatively large sample of 120 individuals, found that while initial amygdala activation magnitude after seeing negative images was not predictive of trait neuroticism, a slower recovery time for the amygdala to return to baseline was predictive of neuroticism. Time to recovery, or habituation, has been defined as differential response amplitude to repeated stimuli over time. Neural habituation enables an individual to ignore known information and to focus on novel information. Failure to habituate indicates that the individual may have difficulty learning that an environment is familiar or predictable. Likely, a longer time course for habituation of amygdala activation may suggest greater difficulty recovering from emotionally evocative stimuli, which is in line with the phenotypic presentation of trait neuroticism.
There is some initial evidence that habituation may be a sensitive individual difference marker. One study found that differences in magnitude between left and right amygdala activation could actually be explained by differences in habituation rates (Phillips et al., Reference Phillips, Medford, Young, Williams, Williams, Bullmore and Brammer2001). Differences in habituation may, partially, explain genetic differences seen in individuals with certain polymorphisms (5-HTTLPR genotype group) associated with amygdala reactivity (Lonsdorf et al., Reference Lonsdorf, Golkar, Lindstöm, Fransson, Schalling, Öhman and Ingvar2011). One recent study found that amygdala habituation (during negative emotion processing) was a more reliable neural marker, exhibiting higher within-subject reliability in test–retest, than magnitude of amygdala activation. This finding suggests that amygdala habituation might be more well suited than magnitude of amygdala activation to individual difference research on dimensional constructs (Plichta et al., Reference Plichta, Grimm, Morgen, Mier, Sauer, Haddad and Meyer-Lindenberg2014).
Brain Functional Connectivity and Neuroticism
Multiple brain regions are often involved in cognitive and affective processes. As such, interactions between brain regions during task-based fMRI may reflect reliable patterns of neural activity associated with psychological constructs (Buckholtz & Meyer-Lindenberg, Reference Buckholtz and Meyer-Lindenberg2012; Friston, Reference Friston2005; Sporns, Reference Sporns2011). Network connectivity models, instead of studying the specialized processing occurring in specific brain regions, examine the coactivation of distributed brain systems (Mesulam, Reference Mesulam1998). Evidence using both structural and functional data suggests that dysconnectivity between limbic and prefrontal regions may be a useful neural correlate for examining individual differences in neuroticism (Bjørnebekk et al., Reference Bjørnebekk, Fjell, Walhovd, Grydeland, Torgersen and Westlye2013; Cremers et al., Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010; Servaas et al., Reference Servaas, Geerligs, Renken, Marsman, Ormel, Riese and Aleman2015; Servaas et al., Reference Servaas, van der Velde, Costafreda, Horton, Ormel, Riese and Aleman2013; Xu & Potenza, Reference Xu and Potenza2012). These studies suggest that inefficiencies in the prefrontal regions associated with top-down control of emotion generating regions (such as the amygdala) might account for individual differences in levels of trait neuroticism.
There is strong support for an association between amygdala–prefrontal connectivity, as it relates to emotion dysregulation in the developmental psychopathology literature, suggesting that this might be a transdiagnostic vulnerability factor for the development in childhood and adolescence and maintenance into adulthood of both internalizing and externalizing psychopathology (Beauchaine & Zisner, Reference Beauchaine and Zisner2017). Beauchaine and Gatke-Kopp (Reference Beauchaine and Gatzke-Kopp2012) define emotion dysregulation as “a pattern of emotional experience and/or expression that interferes with appropriate goal-directed behavior.” In this framework, studies have found that altered connectivity between the amygdala and regions of the frontal cortex, such as the medial PFC or the orbitofrontal cortex have been implicated in emotional lability and failures in self-regulation (Churchwell, Morris, Heurtelou, & Kesner, Reference Churchwell, Morris, Heurtelou and Kesner2009; Hilt, Hanson, & Pollak, Reference Hilt, Hanson, Pollak, Brown and Prinstein2011). These altered patterns of connectivity have also been associated with adolescent psychopathology, such as higher anxiety and generalized anxiety disorder (Kujawa et al., Reference Kujawa, Wu, Klumpp, Pine, Swain, Fitzgerald and Phan2016; Monk et al., Reference Monk, Telzer, Mogg, Bradley, Mai, Louro and Pine2008).
In the adult literature, studies have found evidence for an association between the patterns of connectivity in the amygdala and certain prefrontal regions as they relate to neuroticism and other related constructs. Using connectivity analyses (psychophysiological interactions [PPI]) during an event-related negative emotion processing task, Cremers et al. (Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010) found that amygdala–ACC connectivity was inversely correlated with trait neuroticism in a community sample of 60 individuals. Another study, examining trait anxiety, found a similar association between amygdala–ACC dysconnectivity in a sample of 13 men (Kienast et al., Reference Kienast, Hariri, Schlagenhauf, Wrase, Sterzer, Buchholz and Heinz2008). Given the important role that these prefrontal brain regions play in cognitive control of emotions, altered connectivity between the amygdala and areas of the PFC could provide a neural basis for the negative emotionality that is characteristic of neuroticism (Ochsner & Gross, Reference Ochsner and Gross2005). Such findings suggest that dysconnectivity between brain regions involved in emotion generation and brain regions involved in emotion regulation might underlie the experience of negative emotions characterized by neuroticism. These functional dysconnectivity models are supported by structural connectivity studies, typically measured using diffusion tensor imaging, which show that there is a similar pattern of structural dysconnectivity associated with neuroticism (Xu & Potenza, Reference Xu and Potenza2012).
The Current Study
Activation, habituation, and dysconnectivity of the amygdala are all hypothesized neural correlates of trait neuroticism, but the evidence, to date, on the associations between these neural phenomena and trait neuroticism remains, at times, contradictory and inconclusive. The goal of the current study was to clarify the relationship between neuroticism and these three hypothesized neural markers, each of which have been found to play an important role in the experience and regulation of emotionality, using by far the largest sample (N = 663 twins) ever examined for this purpose. Twin participants underwent a neuroimaging assessment and completed a battery of self-report psychological measures, including several measures of neuroticism, allowing us to produce a robust index of trait neuroticism. Participants completed an emotional face-matching fMRI task known to activate the amygdala. We separated twin pairs for the group-level analyses, resulting in two large twin groups (both with n > 270 individuals) and enabling a within-study replication attempt of the neuroticism–neural marker associations in both twin groups.
Based on previous research, we expected to find an association between one or more of the investigated neural markers and trait neuroticism. If we confirmed an association between neuroticism and altered patterns of amygdala activation, this would suggest that the trait might be characterized by differences in perception and processing of emotional information, whereas if we found an association between neuroticism and alterations in habituation or in patterns of brain connectivity to other brain regions, this would suggest deficiencies in the emotion regulation neurocircuitry across the neuroticism continuum. Better characterization of the neural correlates of neuroticism could help improve the current understanding of this key transdiagnostic trait by incorporating meaningful neurobiological data into a construct typically measured using self-report.
Method
Participants and procedures
The current study included 663 same-sex male and female twin participants from the Minnesota Center for Twin and Family Research, which comprises several ongoing, population-based longitudinal twin and family studies. Participants were drawn from two twin cohorts in the Minnesota Twin Family Study (Iacono, Carlson, Taylor, Elkins, & McGue, Reference Iacono, Carlson, Taylor, Elkins and McGue1999; Iacono & McGue, Reference Iacono and McGue2002; Iacono, McGue, & Krueger, Reference Iacono, McGue and Krueger2006; Keyes et al., Reference Keyes, Malone, Elkins, Legrand, McGue and Iacono2009). Both cohorts were first assessed at age 11, with follow-up assessments every 3 to 7 years. The average age across the sample in the current study was 30.40 (range = 23.36–37.48; SD = 5.19), and the average age of the older cohort was 34.79 (range = 33.01–37.48; SD = 1.18) and of the younger cohort was 24.49 (range = 23.36–26.27; SD = 0.66). The sample is 54% female. Both cohorts underwent a comparable comprehensive, multimodal assessment that included interview, questionnaire, laboratory, and neuroimaging components. The neuroimaging assessment included a structural scan and an fMRI scan, during which participants completed an emotion processing task. The study was approved by the University of Minnesota's institutional review board, and all participants were given monetary compensation for their participation.
Trait neuroticism
Participants completed three self-report personality measures that assess indicators of trait neuroticism: the Inventory for Depression and Anxiety Symptoms-II (IDAS-II; Watson et al., Reference Watson, O'Hara, Naragon-Gainey, Koffel, Chmielewski, Kotov and Ruggero2012), the Multidimensional Personality Questionnaire (MPQ; Tellegen & Waller, Reference Tellegen and Waller2008), and the Personality Inventory for DSM-5 (PID-5; Krueger, Derringer, Markon, Watson, & Skodol, Reference Krueger, Derringer, Markon, Watson and Skodol2012).
IDAS-II
We used a shortened version of the IDAS-II (Watson et al., Reference Watson, O'Hara, Naragon-Gainey, Koffel, Chmielewski, Kotov and Ruggero2012) that includes 38 items (rated from 1 = not at all to 5 = extremely) that yielded scores on dysphoria, panic, suicidality, social anxiety, and traumatic avoidance and traumatic intrusions associated with posttraumatic stress disorder. We selected the dysphoria scale (which included questions such as “I found myself worrying all the time”) as an indicator of trait neuroticism; the dysphoria scale reflects general distress and negative affect and correlates highly with Big Five measures of neuroticism (Simms, Grös, Watson, & O'Hara, Reference Simms, Grös, Watson and O'Hara2008; Watson, Gamez, & Simms, Reference Watson, Gamez and Simms2005; Watson & Naragon-Gainey, Reference Watson and Naragon-Gainey2014). Internal consistency for the dysphoria scale in the present sample was excellent (Cronbach's α = 0.92).
MPQ
We used a brief version of the MPQ (Tellegen & Waller, Reference Tellegen and Waller2008) that includes 138 items (rated from 1 = definitely true to 4 = definitely false) that yielded scores on stress reactivity, alienation, aggression, control, harm avoidance, and well-being. We selected the stress reactivity scale (which included questions such as “I easily get upset”) as an indicator of trait neuroticism. The stress reactivity scale reflects a tendency toward being easily upset, having unaccountable mood changes, being nervous/tense, being prone to feeling guilty, being sensitive/vulnerable, and being worry-prone/anxious and correlates highly with Big Five measures of neuroticism (Hankin, Lakdawalla, Carter, Abela, & Adams, Reference Hankin, Lakdawalla, Carter, Abela and Adams2007; Tellegen & Waller, Reference Tellegen and Waller2008). Internal consistency for the stress reactivity scale in the present sample was excellent (Cronbach's α = 0.91).
PID-5
The PID-5 (Krueger et al., Reference Krueger, Derringer, Markon, Watson and Skodol2012) includes 220 items (rated from 0 = very false or often false to 3 = very true or often true) that yielded scores on the domains negative affect, detachment, antagonism, disinhibition, and psychoticism. We selected the negative affect scale (which included questions such as “I worry a lot about terrible things that might happen”) as an indicator of trait neuroticism; the negative affect scale is composed of the facets of anxiousness, emotional lability, hostility, perseveration, (lack of) restricted affectivity, separation insecurity, depressivity, suspiciousness, and submissiveness and the negative affect scale correlates highly with Big Five measures of neuroticism (Krueger et al., Reference Krueger, Derringer, Markon, Watson and Skodol2012). Internal consistency for the negative affect scale in the present sample was excellent (Cronbach's α = 0.94).
Neuroticism composite
We computed a neuroticism composite score comprising scores on the IDAS dysphoria, MPQ stress reactivity, and PID-5 negative affect scales. Based on the instructions delineated in the PID-5, facet and domain-level scores were only calculated if at least 75% of items for a given facet or domain were completed (American Psychiatric Association, 2013). This same rule was also applied to the IDAS-II scales and the MPQ facets, though not explicitly specified in the scoring instructions for those measures, to ensure continuity across measures. Participants were excluded if we could not derive factor scores (i.e., due to one or more missing domain scores) through a factor analysis without data imputation (n = 50). The version of the PID-5 used in the Minnesota Center for Twin and Family Research protocol includes two validity items that are unlikely to be endorsed to ensure careful responding (e.g., “I was born on the moon” and “Two plus two equals five”). Because there were no validity items in the other two measures (IDAS-II or MPQ), but all three measures were typically completed during the same assessment session, if a participant failed to respond correctly to one or both of the PID-5 validity items, his or her data for all three personality measures was removed from the final data set (n = 25).
Not surprisingly, because the maladaptive traits in the PID-5 and symptoms of anxiety and depression from the IDAS-II are unlikely to be endorsed in a normative population, IDAS-II dysphoria and PID-5 negative affectivity were significantly, positively skewed as measured by the Kolmogorov–Smirnov test. In addition, even though it was developed for use with normative populations, MPQ stress reactivity was also positively skewed, using the Kolmogorov–Smirnov test. To correct for robust positive skew, the Blom transformation, a rank-based transformation, was applied to IDAS-II dysphoria, MPQ stress reactivity, and PID-5 negative affectivity scores (Blom, Reference Blom1958). Before and after skew values were IDAS-II dysphoria (pre =.699, post = .418), MPQ stress reactivity (pre = .466, post = .379), and PID-5 negative affectivity (pre = .907, post = .085).
Zero-order correlations among the three neuroticism indicators were moderate to strong (IDAS dysphoria and PID-5 negative affect, r = .48, p < .001; MPQ stress reactivity and PID-5 negative affect, r = .79, p < .001; IDAS dysphoria and MPQ stress reactivity, r = .54, p < .001). A neuroticism factor score was derived from the three neuroticism indicators using exploratory factor analyses conducted in R (R Core Team, 2015). IDAS-II dysphoria had a factor loading of 0.57, MPQ stress reactivity had a factor loading of 0.94, and PID-5 negative affectivity had a factor loading of 0.83 on the general neuroticism factor. Regression analyses were conducted to control for the linear effects of age and sex on the neuroticism factor score and subsequent analyses used this age/sex adjusted neuroticism score. This decision was based on previous studies that have shown age- and gender-based differences in levels of neuroticism (Soto, John, Gosling, & Potter, Reference Soto, John, Gosling and Potter2011; Weisberg, DeYoung, & Hirsh, Reference Weisberg, De Young and Hirsh2011). Individuals’ factor scores from this age- and sex-corrected neuroticism factor were extracted and included in subsequent fMRI analyses, to determine whether robustly measured trait neuroticism related to the hypothesized neural markers of interest.
fMRI acquisition and preprocessing
Imaging was performed using 3T Siemens Trio (n = 102) and Prisma (n = 561) MRI scanners at the Center for Magnetic Resonance Research at the University of Minnesota. Images were collected with a 32-channel head coil, with foam placed between the participant and the coil to reduce head motion. High-resolution structural scans were collected and were used in the current study for localization of function. The images were acquired using the following sequence: magnetization-prepared rapid acquisition with gradient echo, echo time = 3.65 ms, repetition time = 2530 ms, flip angle = 7 degrees, field of view = 256 mm, matrix = 256 × 256, in-plane resolution = 1.0 mm × 1.0 mm, slice thickness = 1 mm, 240 slices, acceleration factor of 2 (GRAPPA). During the structural scan, participants watched a movie of their choice. Each participant's structural scan was reviewed for radiological abnormalities, and participants with clinically significant findings, as determined by radiological review, were removed from the final data set, as were participants with brain anatomical deviations deemed to be significant enough to potentially alter brain functionality (e.g., sizable visible cysts; n = 32). In addition, n = 8 participants were removed from the final data set because of problems associated with the fMRI task (e.g., the task did not start properly or the participant was unable to see the task due to poor eyesight).
A mirror on the head coil enabled participants to view the behavioral task projected onto a rear-projection screen at the head of the scanner bore. Functional scans were collected on the Trio and Prisma scanners using a T2*-sensitive echo planar sequence (EPI, echo time = 27 ms, repetition time = 2.5 s, flip angle = 80 degrees, field of view = 200 mm, matrix = 64 × 64, slice thickness = 3.1mm with a 20% gap, in-plane resolution = 3.1 mm × 3.1 mm, 43 transversal slices, interleaved slice acquisition, 120 volumes). The phase encoding direction was posterior to anterior. Immediately prior to acquiring the task data, a short, 10-volume echo-planar scan was collected using the same parameters and positioning as the task data, but with opposite phase encoding (anterior to posterior). This opposite phase encoded scan was later used to correct geometric distortions in the task data.
Functional data were analyzed using FMRIB's Software Library (FSL 5.0.9; www.fmrib.ox.ac.uk/fsl/). Each functional data set was registered to the relevant structural data set using a rigid-body (6 degrees of freedom) linear transformation, and each participant's structural data were registered to a standard coordinate space (Montreal Neurological Institute's MNI 152 2-mm volume) using a full affine transformation to allow cross-participant comparisons in a common space.
Functional data were preprocessed using the following steps: (a) geometric distortion correction using FSL's topup and applytopup (for this process, the first 10 images from the task data were extracted from the data set, and they were combined with the 10 images from the AP data set and submitted to the topup algorithm to create an unwarping field; that unwarping field was then applied to the full task data set to produce an unwarped data set); (b) motion correction using the first volume in the functional series as the reference volume; (c) slice timing correction; (d) skull stripping; (e) spatial smoothing using a 6-mm full width at half maximum Gaussian filter; (f) grand-mean scaling; and (g) high-pass temporal filtering with a 100s cutoff. After preprocessing, each participant's data run was evaluated for excessive motion using a tool in FSL, fsl_motion_outliers, which can detect time points in an fMRI data set that have been corrupted by large motion (any volume exceeding the 75th percentile + 1.5x the interquartile range). The default metric was utilized, which involves examining the root mean squared intensity difference of volume N relative to the reference volume. No participants showed motion outliers on greater than 25% of the volumes across all task conditions (the predetermined exclusion cutoff); as such, no participants were excluded for excessive motion. In addition, a covariate of noninterest file was produced for each participant that included the three linear and three rotational motion estimates produced by the motion correction step.
Brain masks were prepared for use in subsequent data extraction and small volume correction analyses. A mask was created for the bilateral amygdala and was defined structurally by the Harvard–Oxford Subcortical Atlas (at 100% probability). Based on previous studies, which have found altered neuroticism to be associated with altered connectivity between the amygdala and the medial PFC, a medial PFC region of interest was selected using coordinates from a study of the regulation of negative emotions (Diekhof, Geier, Falkai, & Gruber, Reference Diekhof, Geier, Falkai and Gruber2011). This region in the ventral-medial PFC (vmPFC) is the result of a coordinate-based quantitative meta-analysis on 49 studies examining neural correlates of emotion regulation in response to negative emotions. In the meta-analysis, the authors found that downregulation of emotion in this region was associated with reduced activation in the amygdala, suggesting that it might be a central hub for amygdala–cortical emotion regulation circuitry. Though studies looking at emotion regulation circuitry in the developmental literature have identified alternative regions of interest such as the ACC or the ventrolateral PFC, the region identified by Diekhof et al. represents a consensus in healthy adults of a primary region associated with the downregulation of emotion from the PFC to the amygdala across various emotion regulation strategies (Kujawa et al., Reference Kujawa, Wu, Klumpp, Pine, Swain, Fitzgerald and Phan2016; Monk et al., Reference Monk, Telzer, Mogg, Bradley, Mai, Louro and Pine2008). Using FSL, a spherical mask of the region identified by Diekhof et al. was constructed, with an 8-mm radius around the coordinate (x = 0, y = 40, z = –18 in MNI space).
fMRI task
During the fMRI session, participants completed a version of an emotion processing task adapted from Hariri et al. (2002). The task has been shown to reliably activate the amygdala in response to negative emotion viewing (Sauder, Hajcak, Angstadt, & Phan, Reference Sauder, Hajcak, Angstadt and Phan2013). During the task, participants completed shape-matching and emotional face-matching trials in 30-s blocks. During shape trials, participants saw a shape at the top of the screen and were instructed to choose which of two shapes at the bottom of the screen matched the shape at the top of the screen. Stimuli were solid black shapes (circle, horizontal ellipse, and vertical ellipse) on a backing white rectangle. Participants used their right hand to respond, pressing one of two buttons on a Current Designs button box to indicate whether the left or right stimulus matched the top stimulus. Emotional face-matching blocks were identical to shape-matching blocks other than the stimuli presented. Face stimuli included black-and-white photos of White male and female actors posing facial expressions of anger or fear (Ekman & Friesen, Reference Ekman and Friesen1977). Participants could not use identity to match the faces as three separate actors (same sex) were presented on each trial. Instead, participants were instructed to choose the actor on the bottom of the screen “who feels the same way” as the target actor at the top of the screen. Only fear and anger expressions were included. A total of four female actors and four male actors were used, with both anger and fear expressions for each actor. The selection and order of stimuli was random across participants, with the constraint that an equal number of fear and anger targets and an equal number of male and female actors were included in each block. Each block included six trials of the same condition (shape matching or emotional face matching). On each trial, the three stimulus images were presented on a solid gray background for 4500 ms. Each trial ended with a 500-ms interstimulus interval that included only the gray background. Responses were allowed for 4850 ms from the onset of the stimulus. In addition to shape and emotional face trials, three 20-s blocks of rest trials were included. Rest blocks included four trials. Stimuli included three white rectangles, identical in size and location to the shape and face stimuli. Rest stimuli were presented for 4500 ms, followed by a 500-ms interstimulus interval presenting a solid gray screen. No response was required from participants during rest trials. Block types alternated throughout the task, following a fixed order of rest (R), shape (S), and emotional face (E) blocks across all participants (R-S-E-S-E-R-S-E-S-E-R). The task lasted 5 min in total. An example image of the task is displayed in Figure 1.
Figure 1. Experimental paradigm. Participants performed a matching task, with two trial types, emotion trials and shape trials, while undergoing functional magnetic resonance imaging. In order to identify amygdala responses to facial expressions and to compare these patterns to amygdala responses during a neutral shape matching condition, participants were asked to match (a) during the emotion blocks, the face on the bottom of the screen that matches the emotion of the face at the top of the screen (e.g., the face at the bottom left of the image matches the face at the top of the A side of the image), and (b) during the shape blocks, the shape on the bottom of the screen that matches the shape on the top of the screen (e.g., the oval on the bottom left of the image that matches the oval at the top of the B side of the image).
Subject- and group-level fMRI task analyses
The final sample for subject-level analyses was n = 548 individuals after excluding participants for brain anatomical abnormalities, missing personality data, invalid personality data, and rare problems with the fMRI task or data collection (see above for more information about exclusion criteria). There were commonalities across the processing stream used for all subject-level data, regardless of which neural marker (response magnitude, habituation, or functional dysconnectivity) we were analyzing. General linear model (GLM) analyses were used to model the task, and while the modeling of the task differed across the three neural markers of interest, rest trials were not modeled in any of the analyses so that they could serve as a baseline. Task predictors included in the different models were convolved with a prototypical gamma-function approximation of the hemodynamic response. The temporal derivative for each task predictor was also added to the models. Predictors of no interest included three linear translation and three rotation motion predictors, as well as participant-specific nuisance predictors for each volume for that participant that exceeded the motion criterion (mean number of volumes that exceeded the motion criteria across participants = 5.2, SD = 3.2). Task predictors were coded as the start and duration of all of the emotion and shape portions of the task. All task predictors were coded in seconds.
The twin data enabled us to conduct a within-study replication attempt of the analyses, in two large subsamples. That is, once analyses were completed at the subject level, they were submitted to separate group level GLM analyses for the three neural markers of interest in the first-born twin group (n = 275) and second-born twin group (n = 273). While not independent samples, confirming the findings in both of these twin groups strengthens the results.
This analytic approach resulted in six total group-level GLMs, with two group-level analyses for each neural marker of interest (magnitude of amygdala activation across the task, habituation of amygdala activation across the task, and amygdala–prefrontal dysconnectivity) in each twin group. Three regressors were included in each of these analyses: a constant regressor of all 1s (to assess the group mean), as well as one predictor of interest, each individual's neuroticism factor score, and one regressor of noninterest for scanner type (Prisma or Trio), to control for the fact that two different scanner types were used over the course of the study (Han et al., Reference Han, Jovicich, Salat, van der Kouwe, Quinn, Czanner and Fischl2006; Jovicich et al., Reference Jovicich, Czanner, Han, Salat, Kouwe, Van Der, Quinn and Fischl2009). Scanner type was identified by a 1 or a 0. The values for neuroticism and scanner type were grand mean-centered. Both the activation and habituation models were performed using a small-volume correction with the mask of the bilateral amygdala (defined structurally using the Harvard–Oxford Subcortical Atlas with a probabilistic threshold of 1.0) and the group-level statistical maps were then cluster-thresholded at z = 2.3, p < .05. The amygdala–prefrontal dysconnectivity analysis was performed using a small-volume correction with the mask of the vmPFC (defined based on the coordinates from Diekhof et al., Reference Diekhof, Geier, Falkai and Gruber2011) and the group-level statistical map was cluster-thresholded at z = 2.3, p < .05.
Magnitude of amygdala activation
The first subject-level analysis modeled magnitude of amygdala activation across the task. In this model, two task predictors of interest were included: shape and emotion. The contrast of interest compared activation during emotion blocks to activation during shape blocks (emotion > shape). The within-subject output from the first subject-level GLM was passed to two separate group-level GLMs, one composed of first-born twins and one composed of second-born twins. Next, these analyses were run with trait neuroticism as a correlate and controlling for scanner type.
Habituation of amygdala activation
The second subject-level analysis modeled habituation of amygdala activation over the course of the task. Habituation is defined as the reduction in response over time. In this case, it refers to a decrease in BOLD signal in the amygdala over the course of the emotion blocks during the task. In order to model this, the task session was divided into four emotion blocks and the BOLD response for each of the four blocks was averaged within each block. The four shape blocks during the session were averaged into a single value, based on the goal of reducing the number of predictors in the model and the premise that any neural habituation to the shape stimuli over the course of the task was not of interest in the analysis. In this model, five predictors of interest were included: one shape predictor (accounting for all four shape blocks) and four weighted predictors modeling a different level of emotional arousal (1.50, 0.50, –0.50, –1.50) separately for each emotion block. Task predictors were coded as the start and duration of each of the four emotion blocks and all of the shape blocks. The contrast of interest examined the negative linear trend of the data (habituation), modeling the response to emotion decreasing, relative to fixation, across blocks. Figure 3 depicts a graphical representation of mean activation across the four emotion blocks in the bilateral amygdala in all participants included in fMRI analyses. As hypothesized, it indicates a generally linear trajectory of BOLD activation decreasing across emotion blocks over the course of the task.
Next, the within-subject output was passed to two separate group-level GLMs, one composed of first-born twins and one composed of second-born twins. Next, these analyses were run with trait neuroticism as a correlate and controlling for scanner type.
Amygdala–vmPFC PPI analysis
The third subject-level analysis modeled amygdala–vmPFC dysconnectivity during the emotion block relative to shape blocks using a PPI model. PPIs measure the relationship between activity in a seed region and activity throughout the entire brain related to a specific task component (O'Reilly, Woolrich, Behrens, Smith, & Johansen-Berg, Reference O'Reilly, Woolrich, Behrens, Smith and Johansen-Berg2012). It is one of the primary ways to measure task-relevant functional connectivity. Based on the previous literature, the amygdala was selected as the seed region. In order to conduct the PPI, the bilateral structural amygdala mask was transformed into each participant's functional space using the FLIRT tool found in FSL. The mean time series within the amygdala mask was then extracted for each individual (using the fslmeants tool found in FSL). Subject-level connectivity analyses were conducted using the FEAT tool in FSL. Using the same preprocessing as above, three task predictors were included in the model: (a) a single task predictor for the start and duration times for all emotion and shape blocks (where emotion blocks were coded with a 1 and shape blocks were coded with a –1), (b) a predictor with each individual's bilateral amygdala mean time series, and (c) a predictor quantifying the interaction between the amygdala time course and the task predictor. The predictor of interest results from the interaction between the amygdala time course and the task predictor. This predictor identifies regions that display stronger functional connectivity (or task-related coactivation) with the amygdala for negative emotional faces compared to shapes.
Next, the within-subject output was passed to two separate group-level GLMs, one composed of first-born twins and one composed of second-born twins. These analyses produced mean statistical maps of brain regions showing greater connectivity with the bilateral amygdala during the emotion blocks relative to the shape blocks. Next, these analyses were run with trait neuroticism as a correlate and controlling for scanner type, with the small-volume correction of the vmPFC.
Results
Brainwise maps were produced for both twin groups for magnitude of activation, habituation, and amygdala–whole brain dysconnectivity during the emotion > shape trials of the task. The results of these are presented in Figure 3a, 3b, and 3c. In addition, the regions that displayed significant activation, habituation, and amygdala–whole brain dysconnectivity, along with the maximum Z-value, and the coordinates of peak significance are identified in Table 1. These maps identify activation, habituation, and amygdala–whole brain connectivity in regions that would be expected across large groups of individuals during the task, in the absence of hypothesized predictors (such as neuroticism) or constrained, small volume-correction analyses. These findings suggest that the task was performing as to be expected for all three task analysis methods in both twin groups.
Table 1. Brainwise patterns of activation, habituation, and connectivity in both twin groups
Note: Coordinates listed are in Montreal Neurological Institute (MNI) space.
Amygdala activation
Using a whole-brain analysis and a cluster correction of z = 2.3, p = .050 for the emotion > shape contrast, we observed significant activation in the amygdala, along with multiple other brain regions of the task (suggesting that the task was activating relevant regions, as it has in previous studies). These analyses produced mean statistical maps and clusters of significant regions for each twin group (Figure 3a and Table 1). However, using a small-volume correction (from the bilateral amygdala mask based on the Harvard–Oxford Subcortical Atlas) with a cluster correction of z = 2.3, p = .050, and with neuroticism included in the model as a correlate (and controlling for scanner type), no voxels were shown to be activated at a significant cluster threshold in either twin group. Therefore, we determined no significant relationship between trait neuroticism and activation in the amygdala to the emotion > shape contrast.
Amygdala habituation
Using a whole-brain analysis and a cluster correction of z = 2.3, p = .050 for the emotion > shape contrast, we observed a number of brain regions, including the amygdala, involved in habituation, and these analyses produced mean statistical maps and clusters of significant regions for each twin group (Figure 3b and Table 1). In addition, a negative linear trend of amygdala activation was observed, by examining change in mean amygdala activation over the four emotion blocks of the task (suggesting that, as hypothesized, amygdala habituation was occurring to the emotional faces across participants over the course of the task; Figure 2). However, using a small-volume correction (from the bilateral amygdala mask based on the Harvard–Oxford Subcortical Atlas) with a cluster correction of z = 2.3, p = .050, and with neuroticism included in the model as a correlate (and controlling for scanner type), no cluster-corrected voxels showed a significant association between trait neuroticism and amygdala habituation in either twin group. Therefore, we determined no significant relationship between trait neuroticism and habituation of the amygdala over the course of the task during emotion > shape trials.
Figure 2. (a) Functional activation in emotion > shape contrast in first-born twins (top) and second-born twins (bottom). (b) Model of habituation over four emotion > shape blocks in first-born twins (top) and second-born twins (bottom). (c) Functional connectivity seeding from the amygdala using a psychophysiological interaction in emotion > shape contrast in first-born twins (top) and second-born twins (bottom).
Figure 3. Hariri emotion task habituation across all participants. Graphical depiction of mean activation across all participants from the four emotion trials in the Hariri task. As predicted, we observed a generally linear decrease in blood oxygen level dependent activation across the emotion blocks over the course of the task.
Amygdala–vmPFC functional dysconnectivity
Using a whole-brain analysis and a cluster correction of z = 2.3, p < .050, we observed significant connectivity between the time course of the amygdala and the bilateral occipital cortex during the emotion component of the task. These analyses produced mean statistical maps and clusters of significant connectivity for each twin group (Figure 3c and Table 1). Using a small-volume correction (defined functionally from the coordinates from Diekhof et al., Reference Diekhof, Geier, Falkai and Gruber2011) with a cluster correction of z = 2.3, p = .050 and with neuroticism included in the model as a correlate (and controlling for scanner type), we observed increased task-related functional connectivity during the emotion component of the task between the amygdala and the vmPFC, which was significantly correlated with increasing levels of trait neuroticism. This finding replicated across both twin groups. This finding can be seen visually in Figure 4, which displays an axial brain slice. In this image, the vmPFC region that was constructed based on the region found to be responsible for negative emotion regulation across studies from a recent meta-analysis (Diekhof et al., Reference Diekhof, Geier, Falkai and Gruber2011) is outlined in red. Overlaid on this region, the green area represents the subregion within this vmPFC region in which neuroticism correlated with significant functional connectivity with the amygdala in the first-born twin group, during the negative emotional task component (x = 6, y = 40, z = –18). In addition, within the vmPFC region, the blue subregion represents the area of the vmPFC in which neuroticism correlated with significant functional connectivity with the amygdala in the second-born twin group, during the emotional task component (x = 0, y = 32, z = –18). Thus, we determined a significant association between trait neuroticism and amygdala–prefrontal dysconnectivity over the course of the task during emotion > shape contrasts.
Figure 4. Amygdala–vmPFC dysconnectivity. In this figure, the red indicates the vmPFC region of interest based on the meta-analysis of emotion regulation (Diekhof et al., Reference Diekhof, Geier, Falkai and Gruber2011). Overlaid are the regions of significant task-related functional connectivity between the amygdala and the vmPFC during the emotion component of the task that correlate with trait neuroticism in each of the twin groups.
Discussion
The present study examined whether neural markers (amygdala activation, amygdala habituation, and amygdala–vmPFC dysconnectivity) derived from task-based fMRI during negative emotion processing were related to trait neuroticism, as defined using multiple self-report measures. These neural markers were selected based on previous literature, which has been limited by small sample sizes and inconsistent results. In a much larger sample than has typically been used in analyses of these kind, we failed to find evidence of an association between trait neuroticism and activation magnitude or habituation of amygdala activation during negative emotion processing. However, we did find evidence of increased task-related functional connectivity between the amygdala and the vmPFC during the emotional component of the face-matching task, which correlated with increasing trait neuroticism. Taken together, results of the present study suggest that increasing levels of neuroticism represent alterations in top-down control and regulation of emotions, as evidenced by greater amygdala–vmPFC dysconnectivity, rather than from overactive emotion reaction processes per se.
Magnitude of amygdala activation
Although several studies have reported significant associations between magnitude of amygdala activation and trait neuroticism (Brück et al., Reference Brück, Kreifelts, Kaza, Lotze and Wildgruber2011; Chan et al., Reference Chan, Norbury, Goodwin and Harmer2009; Cunningham et al., Reference Cunningham, Arbuckle, Jahn, Mowrer and Abduljalil2011; Haas et al., Reference Haas, Omura, Constable and Canli2007; Harenski et al., Reference Harenski, Kim and Hamann2009), a number have also reported nonsignificant results (Cremers et al., Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010; Drabant et al., Reference Drabant, McRae, Manuck, Hariri and Gross2009; Haas et al., Reference Haas, Constable and Canli2008; Hyde et al., Reference Hyde, Gorka, Manuck and Hariri2011; Thomas et al., Reference Thomas, Elliott, McKie, Arnone, Downey, Juhasz and Anderson2011). Much of the existing research has been conducted in relatively small samples, which increases the likelihood of both Type I and Type II errors, and may account for these inconsistent findings. The nonsignificant association between magnitude of amygdala activation and trait neuroticism found in the current study, which included a much larger sample than previous studies, is consistent with results of a recent meta-analysis that also failed to find an association across 18 studies examining neuroticism and brain activation during emotion processing tasks (Servaas et al., Reference Servaas, van der Velde, Costafreda, Horton, Ormel, Riese and Aleman2013). Amygdala activation indexes the change in BOLD signal in the amygdala in response to the emotional face stimuli, and higher activation is understood to reflect heightened emotional reactivity in response to emotional stimuli. Servaas et al. suggest that while the amygdala plays an important role in threat detection in response to salient stimuli in the environment, the regions that show alterations across the studies in the meta-analysis are those involved in fear learning (e.g., hippocampus and parahippocampus), anticipation of aversive stimuli (e.g., anterior cingulate cortex and posterior cingulate cortex), and emotion processing (e.g., middle cingulate gyrus and dorsal-medial PFC).
Amygdala habituation
In the current study, we found no association between amygdala habituation and trait neuroticism. Only one previous study found an association between neuroticism and habituation (Schuyler et al., Reference Schuyler, Kral, Jacquart, Burghy, Weng, Perlman and Davidson2014). Habituation refers to the reduction in neural response to a repeatedly presented stimulus over time, and it can be understood as a type of neural regulatory process, as activation goes down as the brain grows accustomed to a certain stimulus (Thompson & Spencer, Reference Thompson and Spencer1966). Though previous research has shown that amygdala reactivity is generally consistent over multiple scanning sessions to the same stimulus (Johnstone et al., Reference Johnstone, Somerville, Alexander, Oakes, Davidson, Kalin and Whalen2005), there is also evidence that within a scan session, amygdala reactivity reduces with time to emotional stimuli (Breiter et al., Reference Breiter, Etcoff, Whalen, Kennedy, Rauch, Buckner and Rosen1996; Strauss et al., Reference Strauss, Makris, Aharon, Vangel, Goodman, Kennedy and Breiter2005). Breiter et al. found that the amygdala habituates rapidly (i.e., within 1 min) to happy and fearful faces and that these habituation effects were maintained after an interstimulus delay of 4 min. Thus, it is possible that the length of the task in the current study obscures the habituation process, which may be occurring more rapidly. However, given that other studies examining the relationship between habituation and psychopathology have used a longer time window and have found positive associations (e.g., Koeningsberg et al., Reference Koeningsberg, Denny, Fan, Liu, Guerreri, Mayson and Siever2014), the current findings, in this large sample, suggest that amygdala habituation does not significantly relate to trait neuroticism.
Amygdala–vmPFC dysconnectivity
The results of the PPI analysis examining the coactivation of the amygdala and the vmPFC were significant, suggesting that the time courses of activation in these two brain regions are correlated and that this correlation is positively associated with individual differences in trait neuroticism. As the vmPFC–amygdala pathway is understood to be central to the neural mechanisms associated with emotion regulation, this significant finding suggests that altered connectivity between brain regions with emotion perception and brain regions associated with emotion regulation may be key to the neural underpinnings of trait neuroticism (Delgado, Nearing, LeDoux, & Phelps, Reference Delgado, Nearing, LeDoux and Phelps2008; Phelps & LeDoux, Reference Phelps and LeDoux2005; Stein et al., Reference Stein, Wiedholz, Bassett, Weinberger, Zink, Mattay and Meyer-Lindenberg2007). This finding is also consistent with the results of the Servaas et al. (Reference Servaas, van der Velde, Costafreda, Horton, Ormel, Riese and Aleman2013) meta-analysis. Servaas et al. offer a model of the relationship between neuroticism and neural activation, based on their meta-analysis, in which individuals high in neuroticism have an overactive fear learning system coupled with difficulties anticipating or predicting negative outcomes. This combination of neural patterns of activation results in uncertainty and higher levels of stress, or a “neurotic cascade” (Suls & Martin, Reference Suls and Martin2005), which is characterized by increased daily problems, higher emotional reactivity to these problems, more mood “spillover” from previous problems, and stronger emotional reactions to recurring problems. As a result of this increased emotional reactivity, there is a need for greater regulatory control in individuals higher in neuroticism. In such a model, the neural activation patterns associated with neuroticism might be better conceptualized as relating to the brain systems involved in downregulating emotions, as opposed to primarily emotion perception and lower order systems, such as threat detection, centered in the amygdala (Haas et al., Reference Haas, Constable and Canli2008; Lemogne et al., Reference Lemogne, Gorwood, Bergouignan, Pélissolo, Lehéricy and Fossati2011; Williams et al., Reference Williams, Brown, Palmer, Liddell, Kemp, Olivieri and Gordon2006).
Previous studies using task-based fMRI provide support for this pattern of neural connectivity associated with individual differences in trait neuroticism. Using fMRI during negative emotion processing in a healthy sample, Heinz et al. (Reference Heinz, Braus, Smolka, Wrase, Puls, Hermann and Büchel2005) found that individuals with the short (s) allele of the human serotonin transporter gene (SLC6A4), a specific polymorphism associated with major depression, showed increased coupling of the amygdala and the vmPFC. In a sample of healthy adolescents, increased coupling of the amygdala and the vmPFC (as well as the amygdala and the ACC and the dorsal-lateral PFC) was associated with trait neuroticism during fear learning (Tzschoppe et al., Reference Tzschoppe, Nees, Banaschewski, Barker, Büchel, Conrod and Flor2014). Cremers et al. (Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010) found that trait neuroticism was positively associated with amygdala and dorsal-medial PFC connectivity during angry and fearful face processing, though they found an inverse correlation between amygdala–ACC connectivity and neuroticism during sad face processing. In individuals with borderline personality disorder, a disorder that has been characterized as an extreme maladaptive version of neuroticism (Samuel, Carroll, Rounsaville, & Ball, Reference Samuel, Carroll, Rounsaville and Ball2013), increased vmPFC–amygdala connectivity was found during fear processing, suggesting one hypothesized mechanism of exaggerated amygdala response, as the vmPFC may not be functioning in its inhibitory role over amygdala activity (Kamphausen et al., Reference Kamphausen, Schröder, Maier, Bader, Feige, Kaller and Tüscher2013).
Another possible explanation for this positive association between the amgydala–vmPFC and trait neuroticism and associated psychopathology is that the process of emotion regulation involves more effortful cognitive control in individuals higher in neuroticism and that the regulatory function may not be as effective at tempering amygdala reactivity, relative to individuals lower in neuroticism (Johnstone, van Reekum, Urry, Kalin, & Davidson, Reference Johnstone, van Reekum, Urry, Kalin and Davidson2007; Ochsner, Silvers, & Buhle, Reference Ochsner, Silvers and Buhle2012; Urry et al., Reference Urry, van Reekum, Johnstone, Kalin, Thurow, Schaefer and Davidson2006). It may also be that this coupling reveals a more intertwined set of functional processes between the limbic emotion centers of the brain and the vmPFC, suggestive of more overthinking, self-referential processing, and ruminating about emotional stimuli (Blair et al., Reference Blair, Geraci, Devido, McCaffrey, Chen, Vythilingam and Pine2008; Cremers et al., Reference Cremers, Demenescu, Aleman, Renken, van Tol, van der Wee and Roelofs2010; Disner, Beevers, Haigh, & Beck, Reference Disner, Beevers, Haigh and Beck2011). These failures in neural processes of emotion regulation associated with neuroticism align with a relatively recent and growing understanding of transdiagnostic psychopathology in some part stemming from emotion dysregulation (Campbell-Sills & Barlow, Reference Campbell-Sills, Barlow and Gross2007). The trend toward this thinking has been underscored by the growing use of emotion regulation-based treatments for various forms of psychopathology, such as dialectical behavioral therapy, emotion regulation skills, and emotion-focused therapy (Berking et al., Reference Berking, Wupperman, Reichardt, Pejic, Dippel and Znoj2008; Gratz & Tull, Reference Gratz, Tull and Baer2010; Greenberg, Reference Greenberg2017; Lynch, Trost, Salsman, & Linehan, Reference Lynch, Trost, Salsman and Linehan2007). As such, the relationship among neuroticism, transdiagnostic psychopathology, and emotion dysregulation could be explained, in part, by alterations in this neural circuitry associated with emotion regulation.
Altered coupling between the amygdala and the PFC is increasingly being understood as part of the neural pathway associated with emotion dysregulation and as a transdiagnostic marker for psychopathology both in adulthood and during development (Beauchaine & Zisner, Reference Beauchaine and Zisner2017; Bruhl, Delsignore, Komossa, & Weidt, Reference Bruhl, Delsignore, Komossa and Weidt2014; Hardee et al., Reference Hardee, Benson, Bar-Haim, Mogg, Bradley, Chen and Perez-Edgar2013; Kim et al., Reference Kim, Loucks, Palmer, Brown, Solomon, Marchante and Whalen2011; Monk et al., Reference Monk, Telzer, Mogg, Bradley, Mai, Louro and Pine2008). Still, the exact prefrontal regions implicated and the directionality of the relationship have been debated, and there is some research that suggests that there may also be a developmental influence on the directionality of the relationship. Previous studies have found hypoconnectivity between the amygdala and the vmPFC in adolescents with anxiety disorders (Hamm et al., Reference Hamm, Jacobs, Johnson, Fitzgerald, Fitzgerald and Phan2014) during a resting-state fMRI study. In addition, one previous study found a similar negative coupling between the amygdala and the vmPFC in anxious children (ages 11–19) during a task that required paying attention to threat versus nonthreat features, whereas in anxious adults (ages 24–48), they found positive coupling between the amygdala and the vmPFC (Gold et al., Reference Gold, Shechner, Farber, Spiro, Leibenluft, Pine and Britton2016). This finding contrasts with the current study and other studies that have found a negative coupling between the amygdala and the vmPFC in adults higher in anxiety (Kim et al., Reference Kim, Loucks, Palmer, Brown, Solomon, Marchante and Whalen2011). This difference in directionality between the findings in the current study and other resting-state studies, which have found negative coupling in adults to be associated with higher levels of anxiety or neuroticism, as opposed to the positive coupling between the amygdala and the vmPFC found by Gold et al. (2016), might reflect different altered pathways associated with passive observation of emotional information (as in the current study) as opposed to a more active condition of threat detection (as in the study by Gold et al.).
Strengths and limitations
The present study had a number of strengths, including its large, population-based sample. Because of both the size and the composition of the sample, the data are well suited for answering questions about neuroticism, which has been shown to be normally distributed in the general population. In addition, this study uses multiple self-report measures to assess trait neuroticism, in conjunction with multiple fMRI indicators. In this way, it is strengthened by its use of multiple measures across two methodological domains and serves to improve our understanding of the underlying construct of neuroticism. Particularly for the research area of neuroticism, this approach is novel, as the lion's share of extant research on trait neuroticism relies on self-report alone. The current approach is in line with methods advocated by the research domain criteria for the National Institute of Mental Health funding and lays the ground work for additional research integrating behavioral constructs and neurobiological approaches (Patrick et al., Reference Patrick, Venables, Yancey, Hicks, Nelson and Kramer2013).
The study design allowed for a within-study replication attempt, in that we were able to attempt to reproduce results found in first-born twins in the second-born twins from the twin pair. Though not based on independent samples, this method strengthens confidence in the conclusions that can be drawn from this project, as findings are generally consistent across the first-born and second-born twins. Of note, even when separated into first- and second-born twin groups, the present sample is much larger than the extant studies on amygdala activation, habituation, and dysconnectivity associated with neuroticism. Still, conducting analyses separately in the first-born twin group and the second-born twin group does not represent a true replication attempt, as the twins are related and therefore nonindependent data points. Of note, there are currently limited tools for fitting biometrical models to interdependent twin data, and there is no way to correct for the nonindependent nature of twin data when all twins are included in a whole-group neuroimaging analysis.
It is important to acknowledge additional limitations of this study method. While the sample was much larger than typical fMRI samples, it was a predominantly White sample, mapping onto the demographics of Minnesota at the time the twins were born. Therefore, the results might not be generalizable to non-White populations. In addition, the sample is population based, suggesting that it might be generalizable to most individuals, but might not necessarily reflect the extreme poles of trait neuroticism in the population. This is an important caveat for trait neuroticism, which might have different neural features at the higher ranges, as some studies have shown the opposite relationship between amygdala–vmPFC in individuals with diagnosed psychopathology associated with high levels of neuroticism (Kim & Whalen, Reference Kim and Whalen2009).
In addition, in the current study, the personality measures employed to derive the neuroticism factor scores were not standard Big Five measures. Rather, the PID-5 negative affect scale was designed to assess maladaptive personality traits and the IDAS-II dysphoria scale was developed to assess dysphoria associated with DSM-IV diagnoses of depression and anxiety. While evidence broadly finds that these measures are associated with normal-range personality measures (Gore & Widiger, Reference Gore and Widiger2013; Watson, Stasik, Ro, & Clark, Reference Watson, Stasik, Ro and Clark2013), due to their clinical focus, these measures might fail to adequately characterize the range of neuroticism-like features in the community sample used in the current study.
Conclusion
Using multiple methods across two measurement domains, we found evidence for a relationship between trait neuroticism and vmPFC–amygdala coactivation during negative emotion processing in a large sample, while failing to find evidence of a relationship between neuroticism and amygdala activation or habituation. These findings suggest that a key element in the neural underpinnings of neuroticism is related to the interplay between the vmPFC, a central brain region involved in emotion regulation, and the amygdala, a brain region involved in emotion perception and processing. This neural interaction indicates that a crucial mechanism underlying neuroticism may be a failure in emotion regulation, rather than a problem with gating emotional information. Failures in emotion regulation have been shown to confer vulnerability for psychopathology, especially in the context of environmental risk. The present study may help explain the utility and increased use of treatments focusing on emotion regulation for various transdiagnostic mental illnesses that are associated with elevated levels of neuroticism. This finding further suggests that focusing on mechanisms of emotion regulation in relationship to trait neuroticism may prove fruitful for identifying additional neuroscientifically derived biomarkers.
Financial support
This research was supported by resources from the University of Minnesota's Center for Magnetic Resonance Research (Grants P41-RR008079, P41-EB015894, and P30-NS076408) and the Minnesota Supercomputing Institute. Research reported in this paper was supported by the National Science Foundation Grant NSF00039202 (to M.H.S.) and the National Institute on Drug Abuse of the National Institutes of Health Grants K01DA037280 (to S.W.), R01DA036216 (to W.G.I.), and R37DA005147 (to W.G.I.).
