Introduction
============

Individuals differ widely in their ability to accurately perceive, understand, and reason about emotional information, and to effectively apply that knowledge to enhance cognitive and behavioral performance---a capacity known as emotional intelligence (EI) ([@nsx088-B73]). People who possess high EI are skilled at reasoning about emotional issues and are effective at using emotional information to solve problems and achieve goals ([@nsx088-B58]). Because EI capacities are often highly predictive of important aspects of social/interpersonal functioning and professional success ([@nsx088-B14]; [@nsx088-B15]; [@nsx088-B22]; [@nsx088-B87]; [@nsx088-B52]; [@nsx088-B65]), it is of critical importance to understand the neurobiological underpinnings of these capacities and the extent to which they can be improved through training. The goal of the present study is to address the existing knowledge gap regarding the neurocircuitry that underlies EI.

Since EI was first proposed in the early 1990s ([@nsx088-B57]; [@nsx088-B73]), two different approaches to conceptualizing and assessing this construct have emerged---the "trait" and the "ability" models ([@nsx088-B94]). The *Trait* model considers EI as a set of emotional self-perceptions and dispositions that are best assessed via introspection and self-report, similar to personality traits ([@nsx088-B69]; [@nsx088-B77]), while the *Ability* model considers EI as a set of cognitive capacities that permit an individual to accurately perceive, reason about, and use emotional information in adaptive ways, akin to other models of intelligence that focus on demonstrated knowledge and performance ability ([@nsx088-B61]; [@nsx088-B58]). Although both conceptual models lay claim to the term "emotional intelligence" and often appear to be addressing similar theoretical domains ([@nsx088-B68]), these two approaches are often only modestly correlated, if at all ([@nsx088-B16]; [@nsx088-B94]). The lack of statistical association between these two approaches has led to some confusion in the field and may be due to differences in measurement approaches or to core differences in the processes being measured ([@nsx088-B68]). Each approach accounts for unique aspects of emotional, personality and intellectual functioning ([@nsx088-B60]), with *Trait* indices accounting for relatively stable self-perceived competencies in handling emotional and social situations, similar to stable self-rated personality traits ([@nsx088-B7]; [@nsx088-B64]; [@nsx088-B66]), while the *Ability* approach emphasizes objective scoring and measured performance capacities with meaningful benchmarks or comparison metrics, and shares considerable variance with other indices of cognitive intelligence ([@nsx088-B94]). Interestingly, a recent meta-analysis showed that *Ability* EI was significantly associated with better performance on a variety of emotional (i.e. "hot" cognition) tasks whereas *Trait* EI measures were not associated with such outcomes ([@nsx088-B34]). On the other hand, *Trait* EI has been shown to predict career success and satisfaction over many years ([@nsx088-B3]). Because emotional processes are multifaceted and complex, and their adaptive use to achieve goals requires the integration of many capacities, including self-awareness, subjective perceptions, reasoning, and skilled behavioral responses, it is likely that any comprehensive approach to EI will require an understanding of both *Trait* and *Ability* aspects of EI. Therefore, in the present study, we include well-validated metrics of both the *Trait* and *Ability* models and assess their association with the functional connectivity of large-scale brain networks.

An influential neurocircuitry model of EI posited that these capacities emerge from the interaction of several key brain regions involved in emotional perception, experience, and decision-making ([@nsx088-B8]), including (but not limited to) the amygdala, insular cortex, ventromedial prefrontal cortex (vmPFC), and anterior cingulate cortex (ACC), which are key nodes of functional brain networks that are consistently linked with emotional experience, perception ([@nsx088-B49]; [@nsx088-B54]; [@nsx088-B33]), and affectively guided decision making ([@nsx088-B10]; [@nsx088-B23]). The first task-based neuroimaging study of EI found that higher scores on a standard measure of *Ability* EI were negatively associated with activation of nodes of this neurocircuitry, particularly the anterior medial prefrontal cortex (PFC), during a social reasoning task ([@nsx088-B72]), suggesting that EI may be directly related to measureable brain activation patterns. Our group reported similar findings for a measure of *Trait* EI in adolescents during a simple emotional face perception task ([@nsx088-B47]), suggesting a common neurocircuitry might be involved in both models of EI. On the other hand, subsequent task-based functional magnetic resonance imaging (fMRI) work from our lab found that only higher *Ability* (but not *Trait*) EI was associated with greater activation of the vmPFC, amygdala, and insula in response to facial cues reflecting changing levels of trustworthiness ([@nsx088-B45]), findings that were essentially supported in other recent work ([@nsx088-B71]). The role of these same brain regions in EI has also been supported by volumetric studies of brain structure, with *Ability* EI showing positive correlations with the volume of the insula and vmPFC ([@nsx088-B46]), while various components of *Trait* EI showed positive correlations with regions of the medial frontal cortex and anterior cingulate cortex ([@nsx088-B50]; [@nsx088-B86]; [@nsx088-B46]). Thus, task-related and structural volumetric studies have begun to map out a consistent set of brain regions that appear to be important for EI, though the specific regions appear to differ slightly depending on whether the *Trait* or *Ability* model is the focus.

As evident in the foregoing paragraphs, the *Trait* and *Ability* models both emphasize complex and highly inter-related capacities for emotion perception, comprehension and regulation, but each show different associations with behavior and brain function, suggesting that it is critical to study the neurocircuitry underlying both models. Because of the multifaceted nature of EI, it is unlikely to be easily reducible to simple patterns of activation in specific brain regions or discrete organization of brain volume and neurocircuitry. Rather, sophisticated capacities such as EI are more likely to be reflected in the *way* that emotion-related perceptual, interoceptive and experiential systems interact with higher-order integrative and problem solving systems of the brain. A first step in characterizing this complexity is to map the patterns of intrinsic connectivity within and among these systems. Even during spontaneous resting cognition, these brain regions are self-organized into intrinsic networks that interact with one another in meaningful ways that clarify the circuitry that underlies task-based functioning ([@nsx088-B80]).

[@nsx088-B85]) recently conducted the first investigation of the association between *Trait* EI and resting state functional connectivity and found that higher EI scores were associated with greater connectivity within the medial regions of the default mode network (DMN). Moreover, they showed that higher *Trait* EI was also associated with greater anti-correlation between certain regions of the DMN and a task positive network, or TPN (i.e. dorsolateral PFC). This is critical, as recent evidence suggests that normal brain development and the acquisition of superior cognitive performance is related to the magnitude of negatively correlated or "push--pull", activity among major cerebral networks. Notably, superior cognitive performance is associated with greater negative correlation (i.e. anti-correlation) between cortical TPNs that are typically involved in attention demanding cognition and the "task-negative" or DMN, which is typically suppressed during attentionally demanding or externally focused cognition ([@nsx088-B43]; [@nsx088-B35]; [@nsx088-B9]). Individuals with greater cognitive abilities tend to show stronger anti-correlations between these networks. In a parallel fashion, it has recently been shown that the strength of anti-correlation between TPNs and DMN increases over the developmental span from childhood into early adulthood ([@nsx088-B21]). On the other hand, lack of anti-correlation or even positive correlation between these normally anti-correlated networks has been associated with developmental immaturity ([@nsx088-B21]), affective disorders ([@nsx088-B56]), schizophrenia ([@nsx088-B97]), and risk for psychosis ([@nsx088-B100]). Thus, greater ability to suppress the DMN during important task-relevant cognitive processing while increasing activity in relevant TPNs appears to be related to better performance in several domains.

While the study of *Trait* EI by [@nsx088-B85] discussed above represents an important first step in addressing the role of these large-scale networks in EI, it did not examine *Ability* EI. That study also only examined functional connectivity between specific network regions of the DMN and TPN; it did not examine interactions within and between intrinsic networks themselves. Further, as it used a Japanese self-report instrument in a college student cohort, it remains unclear whether such anti-correlated networks would extend to the capacities measured by *Ability* EI assessments, and whether their findings would replicate for *Trait* EI among English speaking populations with a broader age range. To address these issues, we chose to examine intrinsic networks in relation to both *Trait* and *Ability* EI. We tested the hypothesis that a "push--pull" relationship between relevant resting state networks (RSNs) may also exist with regard to the networks involved in *Ability* EI. We further hypothesized that both higher *Trait* and *Ability* EI capacities would be directly related to the strength of the anti-correlation between the DMN and the primary networks involved in emotional processing (rather than traditional cognitive TPNs involved in attention and working memory), but that this pattern would be most strongly observed for the *Ability* model, given that it involves actual performance and is, therefore, most analogous to other cognitive ability studies. To provide an additional within-network test of this "push--pull" hypothesis, we also tested the prediction that, with higher *Trait* and *Ability* EI, these networks would also individually show stronger negative connectivity with a range of emotion-related brain regions---as might be expected in association with greater cognitive control of emotion. Based on prior findings from our structural and task-based functional studies of EI from this same sample ([@nsx088-B46]; [@nsx088-B45]; [@nsx088-B95]; [@nsx088-B2]; [@nsx088-B78]), we focused on networks comprised of emotion processing regions that are typically activated in response to emotional stimuli and which we have previously shown to be involved in EI (e.g. insular and orbitofrontal cortex, ventral striatum, amygdala); we also focused on the DMN, which is engaged during self-introspection and Theory of Mind ([@nsx088-B1]) processes that are critical to EI. To test these associations, we applied independent components analysis (ICA) with dual regression to functional magnetic resonance imaging (fMRI) data acquired while subjects were resting in the scanner. We then tested the relationships between network connectivity and total scores on well-established metrics of *Trait* and *Ability* EI.

Materials and methods
=====================

Participants
------------

Seventy right-handed adults (32 women, 38 men), ranging in age from 18--45 years (mean age = 30.9 years; SD = 8.4 years) from the Boston metropolitan area provided written informed consent to participate in a neuroimaging study. Participants underwent initial telephone screening to rule out any significant history of medical, neurological, psychiatric, or substance use disorders. Following consent, eight participants withdrew voluntarily or were dropped from the study due to failure to comply with study requirements. Data from an additional eight participants were excluded due to excess motion or other scan-related artifacts. The final sample consisted of 54 healthy right-handed adults (28 women, 26 men), with an average age of 30.1 (SD = 7.5) years. Participants were well educated, with an average of 15.2 years of formal education (SD = 2.2 years), and all were native English speakers. Some prior task-based fMRI and behavioral data from subsets of these same participants have been published elsewhere ([@nsx088-B46]; [@nsx088-B45]; [@nsx088-B95]; [@nsx088-B2]; [@nsx088-B78]), but the current EI-resting-state network correlational findings are novel and have not been presented previously. Participants were provided nominal financial compensation for their time. This research protocol was reviewed and approved by the McLean Hospital Institutional Review Board and the U.S. Army Human Research Protection Office.

### EI instruments

Prior to scanning, participants completed two computer-administered measures of EI. As a measure of *Trait* EI, participants completed the Bar-On Emotional Quotient Inventory (EQ-i) ([@nsx088-B6]), a 125-item self-report measure that provides a global measure of EI (*Total EQ*). Total EQ was the primary variable of interest from this scale, but five composite subscales were also available, measuring the ability to understand emotions in others (*Interpersonal*), emotional self-awareness and self-confidence (*Intrapersonal*), problem solving and emotional flexibility (*Adaptability*), coping capacity under stress (*Stress Management*), and general optimism and positive outlook (*General Mood*). Raw scores were transformed into standard scaled scores based on the general population norms available in the computerized scoring program ([@nsx088-B6]). Subscale scores were only included in ancillary analyses if the *Total* EQ was found to be significantly related to brain function (which it was not in this study). While other measures of *Trait* EI were available, we selected the EQ-i given that it is one of the most commonly used self-report measures of EI, it is commercially available, and its model has generated considerable research interest ([@nsx088-B25]). In addition, *Ability* EI was assessed with the Mayer--Salovey--Caruso Emotional Intelligence Test (MSCEIT) ([@nsx088-B59]), which measures emotional reasoning and problem solving. The primary outcome variable was the MSCEIT *Total* score. Additionally, four branch scores known as *Perceiving Emotions*, *Facilitating Emotions, Understanding Emotions*, and *Managing Emotions* were also calculated, but only included in ancillary analyses if the MSCEIT *Total* was found to be significantly related to brain function. Scaled scores were calculated and normalized based on the general consensus scoring method described in the test manual ([@nsx088-B59]). The MSCEIT was selected for this study, as it is the preeminent and most widely used *Ability* measure of EI ([@nsx088-B25]).

Functional neuroimaging
-----------------------

All MRI data were acquired on a Siemens Tim Trio 3 Tesla scanner (Erlangen, Germany) with a 12-channel head coil. High-resolution structural images were acquired using a magnetization-prepared rapid acquisition with gradient echo (MPRAGE) pulse sequence with the following parameters: TR/TE = 2.1 s/2.25 ms, slices = 128, matrix = 256 × 256, flip angle = 12°, resolution = 1.0 × 1.0 × 1.33 mm. Gradient echo echo-planar images sensitive to BOLD contrast were acquired using the following parameters: TR/TE = 2.0 s/30 ms, flip angle = 90°, slices = 34, voxel size = 3.5 mm isotropic. During the 6:04 minute resting state fMRI scan participants were asked to remain awake with their eyes open.

Data pre-processing
-------------------

Most data processing was done using the FMRIB Software Library ([@nsx088-B81]; [@nsx088-B99]; [@nsx088-B41]), FSL release 4.1 (FMRIB Analysis Group; Oxford University, UK; <http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/>) and MATLAB R2012a (The Mathworks, Inc, Natick, MA).

First, the raw data underwent quality assurance (QA) to assess for scanner artifacts and subject motion. A de-spiking algorithm (using an in-house program) was used to correct any datasets that had spikes (e.g. slice dropouts or slices/volumes with spikes in intensity). Motion correction was accomplished with MCFLIRT ([@nsx088-B40]) to do the initial QA. Data were excluded from further analysis if there was greater than one voxel motion in any direction or if the raw FMRI data showed variations in image intensity across the axial planes (e.g. intensity rolling due to interactions of motion with RF pulses), no matter how small the movement as these effects are not correctable with any widely available motion correction scheme. Based on the QA, data from eight participants were removed from further analysis, with one dataset showing excessive non-stationary ghosting, one large male subject with excess pulsations (e.g. physiological noise due to respiration), two with excessive scanner electronic noise and four with excessive motion. Pre-processing included: motion correction, slice timing correction, brain extraction using BET ([@nsx088-B79]), spatial smoothing with a Gaussian kernel of full-width half-maximum 7 mm and high-pass temporal filtering with Gaussian-weighted least-squares straight-line fitting with σ = 150 s. Functional MRI data were aligned to the high-resolution MPRAGE scan with 12 degrees of freedom. The MPRAGE scan was aligned with the MNI152 standard brain using non-linear registration. Transformation of the functional results into MNI space was done following concatenation of the two alignments into a single matrix. All spatially normalized fMRI data were re-sampled to 2 mm^3^ resolution.

Statistical modeling
====================

Independent component analysis (ICA)
------------------------------------

FSL MELODIC \[Multivariate Exploratory Linear Decomposition into Independent Components ([@nsx088-B12])\] was used to perform an independent component analysis of the fMRI data from the sample of subjects. This method decomposes the spatially normalized fMRI data from all subjects (e.g. concatenated into a single data matrix) into a set of "independent components" (ICs), with each IC being a distributed set of brain regions with a temporal trace that describes the evolution of that particular spatial pattern over time. The number of spatio-temporal patterns estimated by the group ICA was set to 30, run eight times, followed by a meta-level ICA fed by all of the spatial maps from the eight decompositions following the procedures used in [@nsx088-B80], which identified 26 stable independent components. In this case, each independent component is comprised of a spatial map and a corresponding timecourse and represents a spatio-temporal pattern of activity that is common to the entire set of participants.

Dual regression
---------------

To assess inter-subject variability in the independent components, the group ICA maps are "back-projected" into each subject's fMRI data to identify their unique spatial and temporal pattern for each ICA map. These "subject-specific spatial maps" (SSSMs) and "subject-specific time courses" (SSTCs) corresponding to each IC from the group ICA were estimated using the dual regression procedure with timecourse normalization ([@nsx088-B11]; [@nsx088-B28]). The SSSMs were used to conduct a "within-network" functional connectivity analysis to identify how the shape (e.g. spatial pattern) of a network might relate to EI covariates of interest. This analysis proceeded similarly to higher level group analyses of task-based fMRI, only instead of activation contrast maps for the subjects being passed to the group level analysis that would be done in a task fMRI study, in our resting state study, the SSSMs for the subjects were collected together (into a 4D data file, with subject being the fourth dimension) and submitted to a higher level general linear model of the group effects.

Additionally, the SSTCs were used to conduct a "between-network" functional connectivity analysis. In this case, for the RSN pairs of interest, the timecourses of the networks in each subject were used to calculate the within-subject partial correlation between the networks (e.g. taking all other components into account) to give one correlation value between each network pair in each individual. These between-network correlation values were collected together across all the subjects and tested to see if the strength of the correlation between a pair of networks was related to *Trait* and *Ability* EI (i.e. EQ-i and MSCEIT scores, respectively).

Estimation and inference
------------------------

To test for statistically significant linear relationships between within-network functional connectivity and EI, the 4D collection of SSSMs for each RSN of interest was submitted to a group general linear model (GLM) with the design matrix setup recommended on the FSL website (<http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/GLM#Single-Group_Average_with_Additional_Covariate>), using two-sided tests. The model included the two covariates as measures of EI: EQ-i *Total* score and MSCEIT *Total* score (which were not correlated in this sample, *r* = 0.05, *P* = 0.71). As the between-subjects effects may be non-Gaussian and/or non-stationary, which can greatly decrease the performance of standard parametric inference methods ([@nsx088-B36]; [@nsx088-B37]), we employed non-parametric permutation testing for estimation and inference ([@nsx088-B29]) using FSL Randomise. A novel permutation-based method to correct for the number of tests across both number of voxels and number of components was implemented ([@nsx088-B53]) to control the FWE across voxels and components at *P *\< 0.05 for the within-network functional connectivity analysis (which was assessed for four different networks). The procedure was implemented as follows: (i) for each component, the same seed was set in the call to Randomise to ensure that the permutations would be identical for all 5000 permutations for every component, with cluster-mass based thresholding (with cluster-forming threshold of 2.3, e.g. -C 2.3 options), and --N and −P options set to output the null distribution text file for each component (ii) the output null distribution text files for all of the components that were tested (e.g. a subset of the full set of independent components) were loaded into Matlab (The Mathworks, Natick, MA). Each text file contained a 5000-element long vector of the maximum cluster-mass statistic over all voxels obtained for each permutation (which are the values to construct the null distribution that is used to control the family-wise error over all voxels), (iii) a maximum of the maximums data matrix (a vector containing 5000 values) was created by recording for each permutation, the max (cluster-mass) value obtained across components, with (iv) the resulting max of the max values were used to construct the null distribution for the maximum cluster-mass statistic across voxels and components. This information then was used to control the FWE across voxels and components at *P* \< 0.05 by identifying the threshold corresponding to the 95th percentile.

The primary hypotheses were restricted to *Total* EI scores for the MSCEIT and EQ-i, not the individual sub-scales of each measure. However, following a hierarchical analysis strategy, only RSNs that showed significant relationships between within-network functional connectivity and *Total* EI scores were further probed to identify which subscales of each test (e.g. four MSCEIT branch subscales; five EQ-i facet subscales) were driving the results. These ancillary analyses were included only if the parent scale was significant and are provided merely for purposes of clarification, but no a priori hypotheses were posed for the subscales. As described below, since only the MSCEIT *Total* score was significant, only the four branch sub-scales of that test were subsequently analyzed. For each of these sub-scales, we show plots of the regression weights averaged over the regions showing a statistically significant linear relationship with Total EI (e.g. "Functional Connectivity") versus the sub-scale value.

Relationships between covariates of interest and between-network functional connectivity were assessed using the FSLNets toolbox ([@nsx088-B83]) (<http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSLNets>). Briefly, for each subject, the partial correlation coefficient (controlling for the contributions of other effects represented in other SSTCs) between the SSTCs for each pair of RSNs of interest was calculated. For each pair of networks, the relationship between covariates of interest (EQ-i and MSCEIT *Total* scores) and the between-network partial correlation coefficients was assessed using FSL Randomise (with a design matrix setup the same as for the spatial map group GLM), accounting for FWE by correcting for multiple tests (e.g. network pairs and contrasts using a Bonferroni correction) with the significance set at *P* \< 0.05.

Results
=======

Association between trait and ability measures
----------------------------------------------

The EQ-i (Trait EI) and MSCEIT (Ability EI) were first examined with a zero-order correlation to determine their association in this sample. As expected from prior work ([@nsx088-B16]; [@nsx088-B94]), the correlation between EI measures was non-significant (*r* = 0.052, *P* = 0.71).

ICA
---

RSNs of interest in the current study were those that contained ventral PFC and limbic/paralimbic regions, including basal ganglia, amygdala, anterior cingulate, hippocampus, anterior thalamus, orbitofrontal cortex (OFC) and nucleus accumbens (NAcc), which are involved in emotion processing ([@nsx088-B63]; [@nsx088-B46]; [@nsx088-B19]; [@nsx088-B45]; [@nsx088-B85]; [@nsx088-B95]; [@nsx088-B2]; [@nsx088-B78]), as well as two sub-networks of the DMN that have been reported to be involved with episodic memory retrieval, mentalizing, Theory of Mind, and autobiographical thought ([@nsx088-B5]). Prior to assessing for statistically significant relationships between network connectivity and EI, four networks, shown in [Figure 1](#nsx088-F1){ref-type="fig"}, were selected that contained these sub-regions:

![3D views and several axial slices are shown for the: (A) *Anterior default mode network (A-DMN)*: ventral and dorsal medial prefrontal cortex (PFC), anterior and posterior cingulate (A/PCC), precuneus, caudate and posterior parietal regions, (B) *Basal ganglia/limbic network*: limbic-paralimbic-striatal regions including inferior insula, basal ganglia (striatum, pallidum, NAcc, subthalamic nucleus, and substantia nigra), amygdala, and thalamus, (C) *Posterior default mode network (P-DMN):* precuneus/PCC, posterior parietal, and frontal cortex (including Brodmann area 8), (D) *Reward learning network:* ventromedial and orbital PFC, caudate, and nucleus accumbens. RSNs shown in red-yellow and overlaid onto the MNI 152 2 mm standard brain (*Z* \> 2.6).](nsx088f1){#nsx088-F1}

1.  *Anterior default mode network* ([@nsx088-B24]; [@nsx088-B90]): recent findings show the existence of several "default mode" sub-networks that can be parcellated from resting state data ([@nsx088-B17]; [@nsx088-B4]; [@nsx088-B82]). Consistent with prior observations ([@nsx088-B4]), here we find a sub-network comprising ventral and dorsal medial PFC, anterior and posterior cingulate (A/PCC), precuneus, caudate, and posterior parietal regions that is engaged during autobiographical thought ([@nsx088-B4]).

2.  *Basal ganglia/limbic network* ([@nsx088-B55]; [@nsx088-B51]; [@nsx088-B39]; [@nsx088-B53]; [@nsx088-B13]): a limbic-paralimbic-striatal network comprising basal ganglia (striatum, pallidum, NAcc, subthalamic nucleus and substantia nigra), inferior insula, amygdala, and thalamus ([@nsx088-B74]; [@nsx088-B89]) that has been associated with activation during reward and emotional tasks ([@nsx088-B51]).

3.  *Posterior default mode network:* We identified another subnetwork of the default mode comprising precuneus, PCC, posterior parietal, medial frontal cortex (including Brodmann area 8), retrosplenial cortex, and mesial temporal lobe, including amygdala. This network is a memory-related network associated with self-relevant recall ([@nsx088-B76]; [@nsx088-B5]) and evaluation of emotional valence and arousal ([@nsx088-B88])

4.  *Reward Learning* ([@nsx088-B51]): A network comprising ventromedial and orbitofrontal cortex, caudate, and nucleus accumbens that is also activated during reward and emotional tasks ([@nsx088-B51]; [@nsx088-B38]).

Relationships between functional connectivity, EQ-I, and MSCEIT
===============================================================

*Within-network functional connectivity*. There were no significant effects associated with EQ-i (reflecting *Trait* EI) for any of the primary networks of interest. On the other hand, the functional connectivity of each network was related to MSCEIT *Total* (reflecting *Ability* EI) and its branch sub-scales as follows: *Anterior default mode network (A-DMN)*. *Total* scores on the MSCEIT were unrelated to functional connectivity of the A-DMN.*Basal ganglia/limbic network (BGN)*. A negative linear relationship between within-network functional connectivity and MSCEIT *Total* was found for the BGN. [Figure 2A and B](#nsx088-F2){ref-type="fig"} shows brain regions in which lower MSCEIT *Total* scores were related to the strength of functional connectivity with the BGN, including hypothesized emotion processing regions in vmPFC, AI, amygdala, lateral OFC and other widespread regions including pallidum, thalamus, midbrain, subcallosal cortex, paracingulate, ventral ACC, and white matter. [Table 1](#nsx088-T1){ref-type="table"} shows information regarding the clusters showing a statistically significant negative linear relationship with MSCEIT *Total* in [Figure 2A and B](#nsx088-F2){ref-type="fig"}. The plots in [Figure 2C](#nsx088-F2){ref-type="fig"} show the average of the parameter estimate values from the dual regression for each subject versus MSCEIT *Total* scores. From the plots, it is evident that increased connectivity of the BGN with these brain regions in individuals with lower *Ability* EI was driven primarily by scores on the MSCEIT *Managing Emotions* and *Facilitating Thought* branches. These analyses were repeated with the two data points removed that appear to be outliers in the plots. The results held without these two subjects. Table 1.Clusters with cluster size 10 or larger found in regions whose functional connectivity with the basal ganglia network (BGN) were negatively correlated with MSCEIT total scoresCluster nameVoxelsMAXMAX X (mm)MAX Y (mm)MAX Z (mm)R. Brainstem/hippocampus27275.1312−10−12R. Inferior orbitofrontal gyrus9554.933828−16L. Precentral gyrus3114.09−40834R. Caudate1623.94222814L. Hippocampus1373.97−40−34−8R. Hippocampus653.2240−28−10L. Pallidum432.94−1602R. Heschl's gyrus403.3238−2418L. Postcentral gyrus373.26−46−834L. Medial orbiofrontal gyrus302.98−456−8L. Amygdala303.08−160−12L. Superior temporal gyrus303.32−52−2−12R. Inferior temporal gyrus/hippocampus223.4748−14−20L. Middle temporal gyrus163.4−52−16−2R. Middle temporal gyrus113.6748−20−10[^1]*Posterior default mode network (P-DMN)*. Within-network connectivity of the P-DMN was also related to MSCEIT *Total* scores in that the functional connectivity between P-DMN and a spatially distributed pattern of regions, including hypothesized emotion processing structures such as the amygdala, vmPFC, AI, PI, OFC, and other non-hypothesized regions such as the temporal pole, middle temporal gyrus, hippocampus, nucleus accumbens, and the ACC, precuneus/posterior cingulate cortex, and superior parietal regions, was negatively linearly related to MSCEIT *Total*. Interestingly, higher MSCEIT *Total* scores were associated with stronger *anti-correlations* between the P-DMN and these regions (e.g. with negative average regression weights), while lower MSCEIT *Total* was associated with *positive correlations* between the *P*-DMN and these various regions (see [Figure 3](#nsx088-F3){ref-type="fig"}). The correlations with MSCEIT *Total* were primarily driven by the *Managing Emotions* branch. All analyses were repeated with removal of the single data point that appears to be an outlier in the plots; all results held without this subject. [Table 2](#nsx088-T2){ref-type="table"} shows information regarding the clusters showing a statistically significant negative linear relationship with MSCEIT *Total* in [Figure 3A and B](#nsx088-F3){ref-type="fig"}. Table 2.Clusters with cluster size 10 or larger found in regions whose functional connectivity with the posterior default mode network (*P*-DMN) were negatively correlated with MSCEIT Total scoresCluster nameVoxelsMAXMAX X (mm)MAX Y (mm)MAX Z (mm)R. Hippocampus58716.1626−18−8L. Middle temporal gyrus13905.16−38−502L. Insula10334.84−34216R. Middle occipital gyrus8144.8530−6636L. Cuneus5433.99−20−5428L. Precentral gyrus4224.61−361042L. Middle cingulate gyrus2564.24−8−1842L. Inferior occipital gyrus2564.53−22−104−6L. Cerebellum crus 12033.76−36−84−18L. Supramarginal gyrus1874.75−44−3430R. Paracentral lobule1824.34−4272R. Parahippocampal gyrus1633.7324−40−4L. Middle occipital gyrus1604.01−38−7618L. Hemisphere white matter1433.96−24−232R. Cerebellum area 61214.6936−58−20L. Cerebellum area 61174.14−34−62−22R. Middle cingulate gyrus1154.2312−2840L. Middle occipital gyrus1143.55−30−7030L. Hemisphere white matter873.48−22−4042L. Superior occipital gyrus783.44−18−9820R. Gyrus rectus763.441020−12L. Hippocampus724.73−36−8−14R. Postcentral gyrus714.0328−3668R. Middle frontal gyrus643.91344412L. Cerebellum Area 10623.52−22−32−42R. Middle frontal gyrus584.2624462R. Posterior cingulate gyrus533.618−4014L. Thalamus524.39−16−64R. Superior frontal gyrus503.49185018L. Brainstem433.36−16−12−8R. Middle temporal gyrus273.1952−24−4L. Middle frontal gyrus253.12−244032R. Caudate253.312242R. Superior temporal gyrus243.6450−4−2L. Hemisphere white matter222.96−18−1634L. Rolandic operculum203.19−6266R. Precuneus173.456−5014L. Cerebellum area 6163.19−26−52−20L. Middle occipital gyrus153.74−28−7842L. Precentral gyrus153.2−56438R. Caudate153.186160L. Superior temporal gyrus143.26−54−66L. Middle frontal gyrus133.31−342058R. Precuneus123.1110−5670R. Calcarine cortex122.8930−6012R. Superior parietal cortex113.120−6470L. Calcarine cortex113.05−14−5012L. Precuneus113.07−4−4410L. Pallidum113.36−26−6−4R. Insula113.333814−16[^2]*Reward Learning (RL)*. Total scores on the MSCEIT were unrelated to functional connectivity with this network.

![Basal ganglia/limbic network. (A) Three orthogonal slices showing regions whose functional connectivity with the basal ganglia network (green) is negatively linearly associated with MSCEIT scores (i.e. reflecting *Ability* EI; red--yellow, *P* \< 0.05 corrected). The crosshairs show one region in the medial PFC. (B) 3D view to better illustrate the spatial extent/pattern. (C) Functional connectivity values (regression coefficients averaged over voxels in red--yellow) plotted against MSCEIT: Total and MSCEIT subscales for Facilitating, Managing, Perceiving, and Understanding. Partial correlation coefficients (PCC) between functional connectivity values and MSCEIT scores are shown for each subscale.](nsx088f2){#nsx088-F2}

![Posterior DMN network. Two orthogonal slices for two different x locations showing regions whose functional connectivity with the DMN is negatively linearly related to MSCEIT scores (i.e. reflecting *Ability* EI; red--yellow, *P* \< 0.05, corrected). Crosshairs in (A) bisect the amygdala and crosshairs in (B), the nucleus accumbens. (C) 3D view to better illustrate the spatial extent/pattern. (D) Functional connectivity values plotted against MSCEIT scores (PCCs are shown for each subscale).](nsx088f3){#nsx088-F3}

*Between-network functional connectivity*. For this sample of healthy adults, the A-DMN was strongly significantly correlated (*P* \< 0.05 corrected) with both the BGN (partial correlation coefficient or PCC = 0.93) and the Reward Learning network (PCC = 0.95). In addition, the P-DMN was significantly anti-correlated with the BGN network (PCC = −0.75). These between network relationships are illustrated in [Figure 4](#nsx088-F4){ref-type="fig"}. Furthermore, the correlations among these networks showed statistically significant relationships with MSCEIT subscales as follows: *A-DMN with BGN*: partial correlations between these two networks were negatively linearly related to MSCEIT *Total* (*P* = 0.05, corrected) and MSCEIT *Perceiving* (*P* = 0.004, corrected), such that the greater the positive correlation strength between the two networks, the lower the MSCEIT *Perceiving Emotions* score, whereas near-zero correlations between the two networks were associated with the highest MSCEIT scores.*BGN and P-DMN*: negative partial correlation (e.g. anti-correlation) between these two networks was linearly related to MSCEIT *Total* (*P* = 0.03, corrected) and MSCEIT *Managing Emotions*, (*P* = 0.03, corrected), such that the greater the strength of the anti-correlation, the higher the MSCEIT score. The lowest MSCEIT scores were associated with partial correlations near zero.

![(Left) Hierarchical clustering shows between-network correlations. The dendrogram is shown in the top of the figure, the RSNs in the middle, and the correlation matrix at the bottom (correlation coefficients shown in the bottom triangular portion and partial correlations shown in the top triangular portion). (Right) Scatterplots showing the significant association between Ability EI (MSCEIT) scales and the major network connectivity (represented as *z*-transformed partial correlation coefficients in the *y*-axis). A-DMN, anterior default mode network; BGN, basal ganglia/limbic network; P-DMN, posterior default mode network; PCC, partial correlation coefficient; RLN, reward learning network.](nsx088f4){#nsx088-F4}

Discussion
==========

Consistent with our "push--pull" hypothesis with respect to *Ability* EI, we found that higher MSCEIT *Total scores* were associated with more negative intrinsic connectivity values between specific "task negative" and "task positive" networks, and between those networks and many brain regions implicated in emotion processing. However, contrary to our expectations, the Bar-On Total EQ-i (*Trait* EI) was unrelated to resting state connectivity within or between the networks we investigated. Specifically, based on findings from our prior functional and structural studies in EI ([@nsx088-B46]; [@nsx088-B45]; [@nsx088-B95]; [@nsx088-B2]; [@nsx088-B78]), we selected four brain networks associated with key emotional and integrative processes from the results of a group ICA of resting state fMRI data that included: an (i) *Anterior Default Mode* network, (ii) *Basal Ganglia/Limbic* network, (iii) *Posterior Default Mode* network, and (iv) a *Reward Learning* network. The strength of resting state functional connectivity within each of these four networks (within-network connectivity) and with each other (between-network connectivity) was then statistically evaluated against *Trait* and *Ability* measures of EI using a dual regression procedure. Whereas *Trait* EI (EQ-i) was unrelated to brain connectivity in these four networks, *Ability* EI (MSCEIT) was significantly negatively associated with the functional connectivity of the *Basal Ganglia/Limbic* network and the *Posterior DMN* network with other brain regions involved in a wide range of emotion-related processes, including primary threat/salience assessment (e.g. amygdala), visceral sensations (e.g. insula), autobiographical memory retrieval (e.g. posterior cingulate), behavioral control (e.g. lateral orbitofrontal cortex), conflict monitoring (e.g. anterior cingulate cortex), reward processing (e.g. ventral striatum), and many others. In contrast, neither the *Anterior DMN* nor *Reward Learning* network connectivity patterns were associated with *Abilit*y EI. Given that the *Trait* and *Ability* models of EI were uncorrelated here, these findings suggest that the ability to reason about emotional information (*Ability* EI), but not necessarily perceived self-efficacy with regard to emotional competencies (*Trait* EI), is associated with greater negative functional connectivity between diverse brain regions and networks involved in the rapid assessment and interoceptive aspects of emotional experience, autobiographical memory, and regulatory control, but not significantly with intrinsic connectivity of areas involved in reward learning or self-referential processing.

Interestingly, in the P-DMN, a network often associated with autobiographical memory retrieval ([@nsx088-B96]; [@nsx088-B70]), we found that higher *Ability EI* (MSCEIT) was associated with synchronous, but *negative*, correlations between the DMN, and other key affective brain regions, while lower *Ability EI* appeared to involve *positively* correlated signals between this network and the same brain regions. In other words, individuals with high MSCEIT *Total* scores showed anti-correlations between the *P-DMN* and distributed brain regions involved in the perception, experience and regulation of emotion, while individuals with low MSCEIT *Total* scores showed positive correlations between the *P-DMN* and these regions, a finding generally consistent with another recent brain imaging study of *Trait* EI ([@nsx088-B85]). Although the interpretation of negative correlations in resting state data continues to be a subject of debate, when data are processed with the group ICA with dual regression, negative correlations in spontaneous fluctuations across spatially distributed regions may indeed reflect valid neurobiological phenomena ([@nsx088-B30]; [@nsx088-B20]) owing to the fact that this approach does not bring to bear the issue of global signal regression. As such, our findings suggest that, similar to observations of superior task performance in individuals with strongest anti-correlations between the *P-DMN* and Task Positive networks ([@nsx088-B21]), superior *Ability* EI is associated with strong anti-correlations between the *P-DMN* and key affective brain regions, while conversely, individuals with poor *Ability* EI show positive correlations between *P-DMN* and widespread cortico-limbic and striatal regions involved in emotional processing.

The present findings are consistent with numerous other studies that show that superior cognitive performance is often associated with greater suppression of the DMN, or a kind of "push--pull" association between DMN and TPN regions, while sub-optimal, developmentally immature, or clinically aberrant patterns of performance tend to be associated with reduced anti-correlation or even positive correlations between these networks ([@nsx088-B9]; [@nsx088-B21]). For example, among healthy adults, greater anti-correlation between DMN and TPN is associated with better attentional performance ([@nsx088-B43]). Similarly, healthy adults show greater anti-correlation between DMN and TPN than children or adolescents, a pattern that correlates with greater inhibitory control, but the strength of these anti-correlations declines again in older age ([@nsx088-B84]). By extension, even more severe versions of this pattern may be seen among individuals with various forms of psychopathology, including affective disorders ([@nsx088-B56]; [@nsx088-B62]; [@nsx088-B98]) or schizophrenia ([@nsx088-B97]). This is not to say that decreased anti-correlations are necessarily directly associated with psychopathology, but they tend to be associated with lower EI capacities, which in itself may serve as a potential risk factor ([@nsx088-B67]). This possibility remains an open question for further research.

Consistent with these prior studies, we propose that *Ability* EI may also involve a similar relationship to the patterns of connectivity between the *P-DMN* and the emotion-related systems described above. Specifically, our data suggest that a person with high *Ability* EI would be expected to show largely anticorrelated activation between the emotion-related and self-referential P-DMN regions. Such suppression of self-reflective systems when emotion-processing systems are active may facilitate clear boundaries between emotional responses and self-relevant cognitions. On the other hand, our data suggest that a person with low *Ability* EI would be expected to show positive coupling between these systems, perhaps contributing to difficulties separating emotional responses from self-reflective processing. We posit that this neurobiological obscuring of the boundaries between internal self-referential cognitions and emotional responses may be associated ultimately with a difficulty reasoning objectively about emotions and managing them in accord with external constraints (i.e. difficulty maintaining "emotional perspective"). Our findings are consistent with such a conceptualization, but replication in other domains of emotional behavior will be necessary to determine the extent to which these patterns of connectivity contribute to the ability to deal effectively with real-world social and emotional phenomena.

For the *BGN*, individuals with the lowest *Ability* EI (MSCEIT) had the strongest positive connectivity within the network itself and between the network and regions in vmPFC, anterior insula (AI), amygdala and lateral OFC. Prior research suggests that these regions are important in the regulation of emotion. Specifically, the vmPFC and contiguous regions of the rACC/sgACC have been proposed to modulate activation within key affective regions such as the amygdala ([@nsx088-B48]; [@nsx088-B91]; [@nsx088-B26]; [@nsx088-B27]; [@nsx088-B75]; [@nsx088-B93]; [@nsx088-B101]), whereas the lateral OFC appears to be involved in regulating emotion through reappraisal processes ([@nsx088-B32]; [@nsx088-B18]; [@nsx088-B92]). The less functionally connected these regions were with the *BGN*, the higher an individual's *Ability EI*. Furthermore, the *BGN* was more strongly functionally connected to the *A-DMN* in individuals with the lowest *Ability* EI. Thus, similar to the *P-DMN* findings above, strong positive coupling between the BGN and affective regions was associated with lower *Ability* EI.

The positive connectivity between emotional regulation regions and the *BGN* may suggest that low EI individuals may have difficulty regulating their emotional experiences and interoceptive sensations via prefrontal control systems. This is consistent with research showing that activation of the vmPFC is negatively related to activation in the amygdala in healthy individuals but is positively coupled in depressed patients ([@nsx088-B42]) as well as among individuals who have not obtained adequate sleep ([@nsx088-B44]). Notably, the association between *Total* EI and the connectivity among these regions was driven most strongly by the *Managing Emotions* branch of the MSCEIT, consistent with the hypothesis that poorer emotion regulation is associated with greater in-phase correlation between the aforementioned regulatory regions and the *BGN* (i.e. when prefrontal activity is high, activity within emotional and interoceptive regions also remains high, rather than being suppressed or effectively modulated). While causal direction cannot be inferred because these data are strictly correlational, they are nevertheless in accord with the growing consensus that the vmPFC plays a modulatory role over primary emotionally responsive regions within the limbic system (e.g. amygdala) ([@nsx088-B48]; [@nsx088-B91]; [@nsx088-B42]). Thus, we speculate that one of the neurobiological mechanisms underlying lower EI capacities may involve ineffective modulation of the *BGN* by prefrontal control regions.

It is also worthwhile to note that while we found significant associations between *Ability* EI (MSCEIT) and functional connectivity of the *BGN* and *P-DMN*, *Trait* EI (EQ-i) was not significantly related to connectivity patterns among the four the networks we studied. This stands in contrast to a recently published seed-based functional connectivity study that found significant correlations between a Japanese measure of *Trait* EI \[i.e. Emotional Intelligence Scale ([@nsx088-B31])\] and functional connectivity between several regions involved in social and emotional processing ([@nsx088-B85]). However, that study differed significantly from the present investigation in several important respects, including their use of a seed-based analysis to investigate the specific connectivity of the medial PFC and the left anterior insula as opposed to our use of ICA with dual regression to investigate the connectivity of the *networks* that contain these regions. In addition, they used a different *Trait* EI measure with a different factor structure, presented in a different language to a more homogeneous group of younger and more highly educated university students. Further research will be necessary to determine the role of these factors in the observed differences between studies. Nonetheless, the present results strongly suggest that higher *Ability* EI capacities appear to be associated with stronger anti-correlations between key networks involved in affective control versus emotional responses, and between emotionally responsive networks and those involved in self-reflective cognition.

Funding
=======

This study was supported by a grant from the U.S. Army Medical Research Acquisition Activity (USAMRAA) to WDSK (W81XWH-09-1-0730) and R01 DA037265 to LDN.

*Conflict of interest*. None declared.

[^1]: Cluster labeling was based on the nearest gray matter label to the cluster maximum according to the Automated Anatomical Labeling Atlas. R = right; L = left.

[^2]: Cluster labeling was based on the nearest gray matter label to the cluster maximum according to the Automated Anatomical Labeling Atlas. L, left; R, right.
