Original Research ARTICLE
Emotion Regulation and Complex Brain Networks: Association Between Expressive Suppression and Efficiency in the Fronto-Parietal Network and Default-Mode Network
- Department of Psychology, Sun Yat-sen University, Guangzhou, China
Emotion regulation (ER) refers to the “implementation of a conscious or non-conscious goal to start, stop or otherwise modulate the trajectory of an emotion” (Etkin et al., 2015). Whereas multiple brain areas have been found to be involved in ER, relatively little is known about whether and how ER is associated with the global functioning of brain networks. Recent advances in brain connectivity research using graph-theory based analysis have shown that the brain can be organized into complex networks composed of functionally or structurally connected brain areas. Global efficiency is one graphic metric indicating the efficiency of information exchange among brain areas and is utilized to measure global functioning of brain networks. The present study examined the relationship between trait measures of ER (expressive suppression (ES) and cognitive reappraisal (CR)) and global efficiency in resting-state functional brain networks (the whole brain network and ten predefined networks) using structural equation modeling (SEM). The results showed that ES was reliably associated with efficiency in the fronto-parietal network and default-mode network. The finding advances the understanding of neural substrates of ER, revealing the relationship between ES and efficient organization of brain networks.
Emotion regulation (ER) refers to “the processes that influence which emotions we have, when we have them, and how we experience and express them” (Gross, 1998). People adopt a wide variety of ER strategies, such as situation selection, situation modification, attentional deployment, cognitive change and response modulation (Gross, 1998). In face of the complexity of ER, Gross and John (2003) suggested to focus on a smaller number of well-defined strategies: “Our focus on two specific, well-defined processes is predicated on the belief that our understanding of complex emotion regulatory processes is best advanced if we focus intensively on one or two processes at a time”. The two major ER strategies or forms as suggested by Gross and John (2003) are expressive suppression (ES) and cognitive reappraisal (CR). ES refers to the alteration of one’s response to an emotional incident, while CR refers to the change of how one views the emotional incident in order to alter one’s feelings. For example, when feeling sad after a heart-broken breakup, one could withhold any sad expression in order to control the intense feelings (i.e., ES); or he/she could persuade himself into thinking that the finished relationship is perhaps actually good for both (i.e., CR). Using ER tasks that typically instruct participants to use the ES or CR strategy when presented with emotionally negative stimuli or adopting the ER questionnaire (ERQ) that was specifically developed to assess the habitual ES and CR usage (Gross and John, 2003), CR has been shown to produce affective, cognitive and social consequences that are more beneficial, whereas ES has been consistently linked to more detrimental consequences such as depressive symptoms (Gross, 2002; Gross and John, 2003; Goldin et al., 2008).
Neural mechanisms of ER have been under investigation using brain-imaging approaches. Task-based functional magnetic resonance imaging (fMRI) studies have shown the involvement of the ventrolateral prefrontal cortex (vlPFC), inferior frontal gyrus (IFG), insula and amygdala in ES (Goldin et al., 2008; Lee et al., 2008) and the engagement of the dorsomedial PFC (dmPFC), dorsolateral PFC (dlPFC), vlPFC, insula, temporal cortex, parietal cortex and amygdala in CR (Kalisch, 2009; Diekhof et al., 2011; Buhle et al., 2014; Gross, 2015). Using the ERQ, studies analyzing gray matter volume and surface thickness have shown that brain structural variations in the ventromedial PFC (vmPFC), dorsal anterior cingulate cortex (dACC), dmPFC, superior frontal gyrus (SFG) and insula are associated with ES (Welborn et al., 2009; Giuliani et al., 2011a; Hermann et al., 2014; Wang et al., 2017), and CR is related to structural variations in the superior frontal cortex (SFC), vmPFC, dACC and amygdala (Welborn et al., 2009; Giuliani et al., 2011b; Hermann et al., 2014; Moore et al., 2016). Moreover, using the ERQ and resting-based fMRI, Wang et al. (2017) found that functional connections between the SFG and regions including the medial PFC (mPFC), precuneus and parahippocampal gyrus were related to ES gender difference.
Previous brain-imaging findings thus suggest the involvement of multiple brain areas in ER and indicate that ER may be related to the functioning of brain networks, which are composed of structurally or functionally connected brain regions (van den Heuvel and Hulshoff Pol, 2010; Smith et al., 2013). Whereas prior studies have investigated functions of individual brain regions in ER, it remains unknown whether and how ER is supported by global functioning of brain networks. To this end, we used global efficiency, one of the graph theory’s global metrics, to assess the role of global functioning of brain networks in ER. There has been a growing interest to investigate the structural and functional organization of the brain by graph-theory based analyses, which consider the brain as organized into complex networks consisting of nodes (brain regions) and edges (structural or functional connectivity between regions; Rubinov and Sporns, 2010; De Vico Fallani et al., 2014; Mears and Pollard, 2016). The topological properties of complex networks can be assessed by a variety of measures (for reviews, see Bullmore and Sporns, 2009; Rubinov and Sporns, 2010; Wang et al., 2010; De Vico Fallani et al., 2014), among which network efficiency is considered to be a “more biologically relevant metric” describing networks in terms of information exchange among brain regions (Wang et al., 2010). Network efficiency can be measured by global efficiency and local efficiency at the global and local level, respectively (Wang et al., 2010; De Vico Fallani et al., 2014; Stanley et al., 2015). Note that global efficiency is inversely related to path length and is suggested to be “easier to estimate than path length when studying sparse networks” (Bullmore and Sporns, 2009); and local efficiency is positively associated with clustering coefficient. Given the purpose of investigating global functioning of brain networks, the present study focused on the metric of global efficiency, which is defined as the average inverse shortest path length in the network (Beaty et al., 2016). Global efficiency has been found to be related to cognitive (e.g., intelligence (Li et al., 2009; van den Heuvel et al., 2009), working memory (Alavash et al., 2015; Stanley et al., 2015)) and social (e.g., personality trait (Beaty et al., 2016) functions, as well as mental disorders (e.g., major depressive disorder (Meng et al., 2014), bipolar disorder (Collin et al., 2016), attention deficit/hyperactivity disorder (Wang et al., 2009)).
Global efficiency can be analyzed for the whole brain network (Li et al., 2009; Wang et al., 2009; Meng et al., 2014; Stanley et al., 2015) or subnetworks (Sheffield et al., 2015; Beaty et al., 2016). For subnetworks, research of resting-state functional connectivity has shown that the brain can be organized into networks composed of functionally connected brain regions (van den Heuvel and Hulshoff Pol, 2010; Power et al., 2011; Raichle, 2011; Cole et al., 2013; Smith et al., 2013). In brain connectivity research the brain is usually divided into a set of non-overlapping regions and the definition of brain networks is related to the brain parcellation (note that a validated parcellation strategy is still lacking; Power et al., 2011; Wig et al., 2011). Power et al. (2011) divided the brain into 264 putative regions and showed that 13 networks based on the 264 regions are in good agreement with major functional systems of the brain. Cole et al. (2013) further reduced the 13 networks into 10 networks by excluding two networks with less clear functionality and combing the “hand” and “face” networks based on consensus of a unified primary motor system. The 10 networks based on the 264-region parcellation were the somato-motor network (SMN), cingulo-opercular network (CON), auditory network (Aud.), default mode network (DMN), visual network (Vis.), fronto-parietal network (FPN), salience network (SAN), subcortical network (Sub.), ventral attention network (VAN) and dorsal attention network (DAN); and have been widely adopted in studies investigating relationships between network properties and cognitive functions (Cole et al., 2014b; Thompson and Fransson, 2015; Uddin, 2015; Vatansever et al., 2015; Long et al., 2016). For the major brain areas that have been found to be involved in ER, the lateral PFC is a core region of the FPN (Cole et al., 2013), the media PFC is a key region of the DMN (Buckner et al., 2008; Broyd et al., 2009), and the ACC is an important component of the CON (Power et al., 2011; Sadaghiani and D’Esposito, 2015). Based on the relationships between individual brain areas and ER and the relationships between individual brain areas and networks, it could be indicated that global efficiency in the FPN, DMN and CON is more likely to be associated with ER. This inference, however, remains to be verified by experimental data and would be difficult to be considered as a strong prior hypothesis. For example, a study examining the Big Five personality traits showed that extraversion and agreeableness were correlated with activities in the midline regions of the DMN and neuroticism, openness and conscientiousness were correlated with activities in the parietal regions of the DMN (Sampaio et al., 2014). However, Beaty et al. (2016) found that only openness was associated with global efficiency in the DMN and suggested that the involvement of a brain region in a cognitive process does not necessarily imply the engagement of global functioning of the network to which that brain region belongs, because individual regions “are related to a wide range of cognitive, behavioral and emotional variables”.
Taken together, the present study aimed to examine the relationship between ER and global functioning of brain networks. The ERQ was used to assess trait measures of the two ER strategies, ES and CR (Gross and John, 2003; Wang et al., 2007). Graph theory-based methods were applied to resting-state fMRI data (Bullmore and Sporns, 2009), and global functioning of the whole brain network and the 10 predefined networks (Cole et al., 2013) was examined by analyzing global efficiency. The association between the ER strategies and global efficiency in the brain networks was evaluated using structural equation modeling (SEM) that models error variance separately from true measurement variance (Beaty et al., 2016). We hypothesized that ER would be associated with global functioning of brain networks as indexed by global efficiency. More specifically, we inferred that global efficiency in the FPN, DMN and CON is likely to be associated with ER. In addition to the FPN, DMN and CON, the whole brain network and other networks in the 10 predefined networks were also analyzed, with the purpose of a more comprehensive understanding of the relationship between global efficiency in networks and ER.
Materials and Methods
The present study recruited 54 participants from Sun Yat-sen University. Two participants were excluded for excessive motion in the MRI scan and four participants were excluded for missing behavioral data, and finally 48 participants (all right-handed, 13 male, mean age ± SD 22.77 ± 1.59) were included in the subsequent analyses. All participants were healthy and reported no history of neurological or psychiatric disorders, or cognitive or affective impairments. This study was carried out in accordance with the recommendations of research protocol approved by the Institutional Review Board of Psychology Department of Sun Yat-sen University with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Institutional Review Board of Psychology Department of Sun Yat-sen University.
MRI Data Acquisition and Preprocessing
The participants were scanned using a Siemens 3.0 Tesla MRI scanner (Siemens, Erlangen, Germany) located at South China Normal University (Guangzhou, China). Whole brain T2*-weighted resting-state functional images were acquired for 8 min as participants relaxed (stayed awake without thinking anything) with eyes closed, using an echo-planar imaging (EPI) sequence: repetition time (TR) = 2000 ms, echo time (TE) = 30 ms, flip angle (FA) = 90°, field of view (FOV) = 224 × 224 mm2, slices = 32, matrix = 64 × 64, slice thickness = 3.5 mm, voxel size = 3.5 × 3.5 × 3.5 mm3, 240 volumes, and interleaved slice ordering. Whole brain T1-weighted structural images were obtained in a sagittal orientation using the magnetization-prepared rapid gradient-echo (MPRAGE) sequence: TR = 2300 ms, TE = 3.24 ms, FA = 9°, FOV = 256 × 256 mm2, inversion time = 900 ms, matrix = 256 × 256, slices = 176, slice thickness = 1 mm, and voxel size = 1 × 1 × 1 mm3.
The data were preprocessed using Statistical Parametric Mapping (SPM121) and Data Processing Assistant for Resting-State fMRI (DPARSF; Yang and Zang, 2010). Preprocessing consisted of standard resting-state functional connectivity preprocessing procedures as implemented in DPARSF, including removing the first 10 volumes of functional images, slice timing correction, motion correction (data of two participants were excluded under the criterion of 2 mm displacement or 2° rotation in any direction), coregistration of structure images to functional images, segmentation with the DARTEL method (Ashburner, 2007), normalization to the standard MNI space with the DARTEL method and resampling functional images at a voxel size of 3 × 3 × 3 mm3, removing linear trends, regressing out nuisance variables (head motion parameters, white matter signals, and cerebrospinal fluid signals), filtering (0.01–0.1 Hz), and spatial smoothing (4-mm FHWM). Following previous work analyzing functional connectivity using graphic approaches (Zhao et al., 2012; Cole et al., 2013; Arnold et al., 2014; Santarnecchi et al., 2014; Beaty et al., 2016), in the present study whole brain signal was not regressed out as a nuisance variable because of the current controversy over global signal regression (Murphy et al., 2009) and the potential impact of global signal removal on topological properties of brain networks (Santarnecchi et al., 2014).
Network Construction and Graphic Analyses
Graph analyses were conducted using graph theoretical network analysis (GRETNA; Wang et al., 2015). The present study adopted the 10 predefined networks based on the parcellation of the brain with 264 cortical and subcortical 10-mm diameter spherical regions (for detailed information, see Power et al., 2011; Cole et al., 2013). For each of the 10 networks in each participant’s dataset, regional time series were calculated by averaging voxel time series in each of the N regions in that network, and a N × N functional connectivity matrix was constructed by computing the Pearson correlation (with Fisher transformation) between each pair of regional time series of the N regions. The brain regions represented network nodes and the correlations between nodes represented network edges. The adjacency matrix was then computed from the correlation matrix by applying a threshold, which resulted in the binary undirected graph. Note that different thresholds would generate graphs of different connection density or sparsity, and network properties are suggested to be examined over a wide range of thresholds since currently there is no “good” threshold for graphic analysis (Bullmore and Sporns, 2009; Langer et al., 2012; Power et al., 2013). Therefore network thresholding in the present study used Fisher transformed r values ranging from 0.15 to 0.85 (Achard et al., 2006; He et al., 2007; Langer et al., 2012; Arnold et al., 2014; Santarnecchi et al., 2014), with the aim of a more comprehensive analysis. Moreover, whether a network would have the property of the small world was evaluated for each threshold. The small-world topology is characterized through the characteristic path length L and clustering coefficient C (Watts and Strogatz, 1998). The characteristic path length is defined as the average of the shortest path lengths between any pair of nodes in the network. The clustering coefficient is defined as the average of the clustering coefficients over all nodes, where the clustering coefficient of a node is defined as the proportion of possible connections that actually exist between the nearest neighbors of a node. Compared with random networks that have the same number of nodes, same mean degree and same degree distribution as the real network (e.g., brain network) of interest, if the real network has higher clustering coefficient and similar characteristic path length (i.e., γ = Cnet/Cran > 1, λ = Lnet /Lran~ 1, with Cnet and Lnet indicating the clustering coefficient and characteristic path length of the real network, respectively, and Cran and Lran indicating the clustering coefficient and characteristic path length of the random network, respectively), the real network would be considered as a small world (Watts and Strogatz, 1998). The ratio between the two parameters γ and λ (γ/λ) can be defined as a parameter σ to measure the small-word-ness, with σ > 1 indicating that a network exhibits small-world property (Humphries et al., 2006). Moreover, before calculating the small-world-ness, it is recommended to first evaluate whether small word property of a network is estimable, which is indicated when the mean degree (k) of the network is larger than the log of the number (n) of the network’s nodes (i.e., k > log(n); Watts and Strogatz, 1998). Note that small-world property tends to be not estimable for higher thresholds. This is because as the threshold increases, more weak connections will be removed from the network and the network will become sparser. As a result, the mean degree will decrease and is more likely to be less than the log of the number of the network’s nodes (Watts and Strogatz, 1998; Achard et al., 2006; He et al., 2007). In the present study investigating brain networks, for each of the whole brain network and the 10 predefined networks, following the procedure described by Jäncke and Langer (2011) and Langer et al. (2012), the correlation matrices of all participants were first averaged to obtain an average correlation matrix and the corresponding adjacency matrix was calculated. Then for each threshold, whether small word property of a network is estimable was evaluated, and the small-word-ness (σ) was calculated. Note that the following analyses only included the thresholds at which small word property of a network was estimable (i.e., k > log(n)) and existed (σ > 1). After that, global efficiency was calculated for each network at each threshold for each participant, which was mathematically expressed as the inverse of the average shortest path length in the network (Wang et al., 2010; Beaty et al., 2016).
Participants completed a Chinese version of the ERQ (Gross and John, 2003; Wang et al., 2007). The Chinese version of ERQ has satisfactory internal consistency and test-retest reliability (Wang et al., 2007). The ERQ had 10 items, with four items measuring ES and the other six items measuring CR. An example item for ES was: “I control my emotions by not expressing them”; and an example item for CR was “When I want to feel more positive emotion (such as joy or amusement), I change what I’m thinking about”. Participants answered to these items on a 1-to-7 Linkert scale (“1” = “strongly disagree”; “7” = “strongly agree”), to indicate their habitual use of these two ER strategies. Participants further completed a Chinese version of negative affect (NA) subscale of the Positive and NA Schedule (PANAS; Watson et al., 1988; Qiu et al., 2008). The Chinese version of NA subscale had nine items, and participants indicated their general frequency of negative mood states on a 1-to-5 Likert scale (“1” = “not at all”, “5” = “extremely”).
Structural Equation Modeling (SEM)
The relationship between trait measures of ER (ES/CR) and global efficiency of a network was examined using SEM, which used latent variables to model error variance separately from true measurement variance (Skrondal and Rabe-Hesketh, 2004; Beaty et al., 2016; Figure 1A). ES and CR were formed as latent variables by specifying their corresponding questionnaire items (four items for ES and six items for CR) as indicators. Age and gender were modeled as observed variables. Note that in order to control for NA, the NA score was also modeled as an observed variable (Giuliani et al., 2011a,b). The analyses without controlling for the NA score are presented in the Supplementary Materials (see “Analyses without modeling the NA score as an observed variable”, Supplementary Tables S1, S2). Standardized regression coefficients were reported. SEM analyses were conducted with Mplus 8 using maximum likelihood robust estimation.
The mean scores of the two subscales of ERQ were: 2.92 (SD = 1.02) for ES and 4.99 (SD = 0.89) for CR. Consistent with previous results, the ES and CR scores did not statistically correlate with each other (r = 0.037, p > 0.05), and the internal consistency was acceptable (Cronbach’s α was 0.73 and 0.86 for ES and CR, respectively; Gross and John, 2003; Kühn et al., 2011; Picó-Pérez et al., 2017). The mean score of the NA subscale was 19.88 (SD = 5.28). Consistent with previous reports (Giuliani et al., 2011a,b), the NA score did not correlate with the ES score (r = 0.14, p = 0.32) or the CR score (r = 0.018, p = 0.90), and no NA outlier (>3 SD) was found.
The relationship between ES/CR and global efficiency of each of the whole brain network and the 10 networks was examined using SEM over the 15 network thresholds. In sum, the data fit the model well regarding the model fit indices, e.g., comparative fit index (CFI), Tucker–Lewis index (TLI), root mean square error of approximation (RMSEA) and standardized root mean square residual (SRMR; Hu and Bentler, 1998; Brown, 2006; see “Supplementary description of SEM” in the Supplementary Materials for further details of the SEM method), providing support for the proposed theoretical models. Whether a network would have small world property was first evaluated for each threshold, and the following analyses only included the thresholds at which small word property of the network was estimable (i.e., k > log(n)) and was presented (σ > 1; see the “Materials and Methods” section for details). In general, consistent with previous findings (Watts and Strogatz, 1998; Achard et al., 2006; He et al., 2007), networks lacked small word property for higher thresholds (Figure 1). Connection density was assessed and listed in Table 1. Global efficiency distributions under different thresholds are listed for all networks in Supplementary Figure S1.
Figure 1. Illustrating the small-world-ness of networks under different thresholds. Whole brain network (WBN), somato-motor network (SMN), cingulo-opercular network (CON), auditory network (Aud.), default mode network (DMN), visual network (Vis.), frontoparietal network (FPN), salience network (SAN), subcortical network (Sub.), ventral attention network (VAN) and dorsal attention network (DAN) represent the WBN, SMN, CON, Aud., DMN, Vis., FPN, SAN, Sub., VAN and DAN, respectively. For each threshold (Fisher transformed r ranging from 0.15 to 0.85 in 0.05 steps; totally 15 thresholds), it was first evaluated whether small word property of a network was estimable. The thresholds at which small word property of the network was not estimable were not included in following analyses and are not presented in the figure. The small-word-ness (σ) was then calculated with σ > 1 indicating that a network exhibits small-world property. The thresholds at which σ > 1 are indicated by black circles, and the thresholds at which σ ≤ 1 are indicated by gray circles and were not included in following analyses.
The results of ES are listed in Table 2. Note that it has been suggested to examine network properties over a wide range of thresholds since currently there is no “good” threshold for graphic analysis, and findings observed across a range of thresholds (in contrast to those occasionally observed at only a single or a few thresholds) would be considered more reliable (Stam et al., 2006; Bullmore and Sporns, 2009; van den Heuvel et al., 2009; Langer et al., 2012; Power et al., 2013; Xu et al., 2015). In the present results the most reliable association was found between ES and FPN, which was statistically significant (pcorrected < 0.05, corrected for 10 networks using Bonferroni correction) for thresholds from 0.15 to 0.4 (although marginally significant at threshold 0.2; Table 2 and Figure 2). The association between ES and DMN would also be considered reliable, which was statistically significant (pcorrected < 0.05) for thresholds 0.15 and 0.2 and was marginally significant for thresholds 0.25 and 0.3 (Table 2 and Figure 3). Moreover, ES was statistically (pcorrected < 0.05) associated with Aud. at threshold 0.15 (the association was marginally significant for thresholds 0.2 and 0.3) and was marginally associated with the whole brain network at thresholds 0.15, 0.2 and 0.25; which may not be as reliable as the associations between ES and FPN and between ES and DMN.
Figure 2. Illustration of association between expressive suppression (ES) and global efficiency in the FPN. Shown are the results for thresholds 0.15 (top row), 0.25 (middle row) and 0.35 (bottom row). For each threshold, structural equation modeling (SEM) is presented on the left, which shows effects of the latent emotion regulation (ER) variables on FPN efficiency. es1–4 refer to the four questionnaire items for ES and cr1–6 refer to the six questionnaire items for cognitive reappraisal (CR). The paths are standardized coefficients. NA indicates negative affect. Scatter plot of the association between latent ES and FPN efficiency is presented in the middle. Note that all variables were standardized by Z-transformation. Nodes and edges that were used to define FPN are presented on the right for a representative participant (the size of nodes represents the degree). Other conventions are as in Table 2.
Figure 3. Illustration of association between ES and global efficiency in the DMN. Results for thresholds 0.15 (top row) and 0.2 (bottom row) are shown. Other conventions are as in Figure 2.
The results of CR are listed in Table 3. In contrast to ES results in which reliable associations across thresholds were found for the FPN and DMN (and less reliable associations for the Aud. and the whole brain network), CR was not statistically associated with efficiency in any network at any threshold.
The present study examined the relationship between trait measures of ER (ES and CR) and global efficiency in resting-state brain networks (the whole brain network and 10 predefined networks) using SEM over 15 network thresholds (although in general networks lacked small word property for higher thresholds). The results showed that ES was reliably associated with efficiency in the fronto-parietal network (FPN) and default-mode network (DMN).
The FPN includes primarily portions of the frontal cortex and parietal cortex, in which the lateral PFC is a core region (Vincent et al., 2008; Dodds et al., 2011; Cole et al., 2013; Smith et al., 2013). The lateral PFC is an essential brain area for flexible control of thoughts and actions (MacDonald et al., 2000; Kerns et al., 2004; Egner and Hirsch, 2005) and is involved in various cognitive control tasks (for reviews, see Matsumoto and Tanaka, 2004; Badre and D’Esposito, 2009). The FPN has been shown to be related to task control processes (Dodds et al., 2011; Power et al., 2011; Cole et al., 2013, 2014b). The involvement of the lateral PFC in ES as shown in previous studies (Goldin et al., 2008) would be in line with the current finding of the association between ES and efficiency in the FPN. The DMN comprises majorly the media PFC, posterior cingulate cortex (PCC) and inferior parietal cortex (Buckner et al., 2008; Broyd et al., 2009; Power et al., 2011). The DMN was initially observed when participants are not focused on an specific task (e.g., in a resting state) and following studies have suggested that the DMN is related to a wide range of spontaneous and self-generated processes, such as episodic future thinking, autobiographical memory processing, mind wandering and thinking about others (Anticevic et al., 2012; Andrews-Hanna et al., 2014; Fox et al., 2015; Hamilton et al., 2015; Stawarczyk and D’Argembeau, 2015). Note that it has been suggested that the DMN is associated with reflecting about one’s own emotional state, and understanding others’ emotions (Andrews-Hanna, 2012; Andrews-Hanna et al., 2014). As a key region of the DMN, the media PFC has been found to be engaged in ES (Goldin et al., 2008; Welborn et al., 2009), which would be in agreement with the present finding of the association between ES and efficiency in the DMN. Combining these findings, it could be possible that for ES, efficiency in the DMN is related to emotional processing of an incident and efficiency in the FPN is associated with control of emotional responses. Furthermore, it was worth noting that an association between ES and efficiency in the CON was not observed, given dACC volume has been found to be related to ES (Dosenbach et al., 2007; Power et al., 2011; Hermann et al., 2014). This, in line with previous results (Sampaio et al., 2014; Beaty et al., 2016), could indicate that global functioning of a brain network may not be engaged in a cognitive process even when a brain region within the network is related to the process. Moreover, the present study systematically explored the whole brain network and 10 predefined networks and ES showed less reliable association with the Aud., which would require further investigation.
The ERQ assesses the habitual ES and CR usage: the more frequent one utilizes ES/CR to regulate his/her emotions, the higher this person will score on the ES/CR scale; and vice versa (Gross and John, 2003). In previous brain structural studies using the ERQ, positive relationships between ES and brain region volumes were often reported (Giuliani et al., 2011a; Kühn et al., 2011; Hermann et al., 2014; Wang et al., 2017); whereas Welborn et al. (2009) did not observe such positive associations but rather found a negative relationship between the vmPFC volume and ES. It has been suggested that the discrepancies among results of different studies could be due to different methodological factors such as whole-brain vs. ROI analysis (Hermann et al., 2014) and voxel- vs. surface-based analysis (Moore et al., 2016). For the current analyses of global functioning of brain networks, ES was negatively associated with efficiency in the FPN and DMN (the association with the Aud. was less reliable). Global efficiency in the whole brain networks has been found negatively associated with some mental disorders (Wang et al., 2009; Meng et al., 2014; Collin et al., 2016), such as major depressive disorder (Meng et al., 2014); and notably, ES has been consistently connected to negative affective and social consequences, including depressive symptoms (Gross, 2002). This could be one plausible clue for the currently observed negative association between ES and network efficiency.
Moreover, association between CR and efficiency in a network was not observed. While the lack of CR association was unexpected, it appeared to be in agreement with a report showing that CR usage as measured with the ERQ was associated with fewer brain regions as compared to ES usage; while ES was found to be related to multiple brain regions including the vmPFC, dmPFC and dACC, CR was only associated with the amygdala (Hermann et al., 2014). The discrepancies in structural (Hermann et al., 2014) and functional network (as in the present study) properties between ES and CR usages require to be further clarified.
The present study would be considered as a preliminary attempt to investigate the relationship between ER and complex brain networks and has many limitations. Besides the limitations as discussed above, the topological properties of complex networks can be assessed by a wide variety of measures (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010; Wang et al., 2010; De Vico Fallani et al., 2014). While the present study focused on the global metric of global efficiency, the topological properties of complex networks would be investigated comprehensively in future ER studies. In particular, roles of individual nodes (regions) in ER remain to be addressed by using metrics such as centrality. Accordingly, the currently adopted 10 predefined brain networks were defined based on 264 putative brain regions. While the 264-region parcellation and the 10 networks are helpful for investigation of network properties in general, investigation of emotion-related functions would require further work. Particularly, an emotion-related network is not defined and amygdala -a key area in emotion processing- is not specifically defined (Power et al., 2011; Cole et al., 2013). Future research could define the region of interest (ROI) of the amygdala structurally or functionally (Poldrack, 2007) and define an emotion-related network by calculating correlation between the amygdala ROI with all other voxels in the brain (van den Heuvel and Hulshoff Pol, 2010). Then topological properties of the emotion-related network could be calculated and their relationships with ER could be examined.
In sum, ER is essential for human adaptive functioning (Gross, 1998; Ochsner and Gross, 2005). Whereas previous brain-imaging studies have investigated functions of individual brain areas in ER, the role of global functioning of brain networks remains unknown. The present finding of the association of ES with global efficiency in the FPN and DMN suggests that efficient organization of specific brain networks is a fundamental neural mechanism for ER. Meanwhile, whereas the results of current graphic analyses of resting-state functional networks showed consistency with the regional results from previous task-based fMRI and structural MRI studies (e.g., association between ES and areas in the FPN and DMN), discrepancies in the results were also revealed and would be addressed in future research (e.g., lack of association between ES and efficiency in the CON, lack of association between CR and network efficiency). Therefore, while the present findings suggest the important role of global functioning of brain networks in ER, it is also indicated that combining different research approaches is required for a better understanding of neural mechanisms underlying ER.
The data generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
LZ, CH, MH and XW designed the research. SZ, LG, JY and YH collected MRI data. LZ, CW, JY and QW analyzed MRI data. CH, XX, QC and HZ collected questionnaire data. CH analyzed questionnaire data. JP and LZ performed SEM analysis. JP, LZ, CH, MH and XW wrote the manuscript. All authors commented on the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by National Natural Science Foundation of China (31371129) and Research Project of Sun Yat-sen University (26000-31620003). We thank Zhengjia Dai for the discussion of MRI data acquisition and analysis.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2018.00070/full#supplementary-material
Achard, S., Salvador, R., Whitcher, B., Suckling, J., and Bullmore, E. (2006). A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J. Neurosci. 26, 63–72. doi: 10.1523/JNEUROSCI.3874-05.2006
Alavash, M., Doebler, P., Holling, H., Thiel, C. M., and Gießing, C. (2015). Is functional integration of resting state brain networks an unspecific biomarker for working memory performance? Neuroimage 108, 182–193. doi: 10.1016/j.neuroimage.2014.12.046
Andrews-Hanna, J. R., Smallwood, J., and Spreng, R. N. (2014). The default network and self-generated thought: component processes, dynamic control, and clinical relevance: the brain’s default network. Ann. N Y Acad. Sci. 1316, 29–52. doi: 10.1111/nyas.12360
Anticevic, A., Cole, M. W., Murray, J. D., Corlett, P. R., Wang, X.-J., and Krystal, J. H. (2012). The role of default network deactivation in cognition and disease. Trends Cogn. Sci. 16, 584–592. doi: 10.1016/j.tics.2012.10.008
Arnold, A. E. G. F., Protzner, A. B., Bray, S., Levy, R. M., and Iaria, G. (2014). Neural network configuration and efficiency underlies individual differences in spatial orientation ability. J. Cogn. Neurosci. 26, 380–394. doi: 10.1162/jocn_a_00491
Beaty, R. E., Kaufman, S. B., Benedek, M., Jung, R. E., Kenett, Y. N., Jauk, E., et al. (2016). Personality and complex brain networks: the role of openness to experience in default network efficiency: openness and the default network. Hum. Brain Mapp. 37, 773–779. doi: 10.1002/hbm.23065
Broyd, S. J., Demanuele, C., Debener, S., Helps, S. K., James, C. J., and Sonuga-Barke, E. J. S. (2009). Default-mode brain dysfunction in mental disorders: a systematic review. Neurosci. Biobehav. Rev. 33, 279–296. doi: 10.1016/j.neubiorev.2008.09.002
Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. N Y Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011
Buhle, J. T., Silvers, J. A., Wager, T. D., Lopez, R., Onyemekwu, C., Kober, H., et al. (2014). Cognitive reappraisal of emotion: a meta-analysis of human neuroimaging studies. Cereb. Cortex 24, 2981–2990. doi: 10.1093/cercor/bht154
Cole, M. W., Bassett, D. S., Power, J. D., Braver, T. S., and Petersen, S. E. (2014a). Intrinsic and task-evoked network architectures of the human brain. Neuron 83, 238–251. doi: 10.1016/j.neuron.2014.05.014
Cole, M. W., Reynolds, J. R., Power, J. D., Repovs, G., Anticevic, A., and Braver, T. S. (2013). Multi-task connectivity reveals flexible hubs for adaptive task control. Nat. Neurosci. 16, 1348–1355. doi: 10.1038/nn.3470
Collin, G., van den Heuvel, M. P., Abramovic, L., Vreeker, A., de Reus, M. A., van Haren, N. E. M., et al. (2016). Brain network analysis reveals affected connectome structure in bipolar I disorder. Hum. Brain Mapp. 37, 122–134. doi: 10.1002/hbm.23017
De Vico Fallani, F., Richiardi, J., Chavez, M., and Achard, S. (2014). Graph analysis of functional brain networks: practical issues in translational neuroscience. Philos. Trans. R. Soc. Lond. B Biol. Sci. 369:20130521. doi: 10.1098/rstb.2013.0521
Dodds, C. M., Morein-Zamir, S., and Robbins, T. W. (2011). Dissociating inhibition, attention, and response control in the frontoparietal network using functional magnetic resonance imaging. Cereb. Cortex 21, 1155–1165. doi: 10.1093/cercor/bhq187
Dosenbach, N. U. F., Fair, D. A., Miezin, F. M., Cohen, A. L., Wenger, K. K., Dosenbach, R. A. T., et al. (2007). Distinct brain networks for adaptive and stable task control in humans. Proc. Natl. Acad. Sci. U S A 104, 11073–11078. doi: 10.1073/pnas.0704320104
Fox, K. C. R., Spreng, R. N., Ellamil, M., Andrews-Hanna, J. R., and Christoff, K. (2015). The wandering brain: meta-analysis of functional neuroimaging studies of mind-wandering and related spontaneous thought processes. Neuroimage 111, 611–621. doi: 10.1016/j.neuroimage.2015.02.039
Giuliani, N. R., Drabant, E. M., Bhatnagar, R., and Gross, J. J. (2011a). Emotion regulation and brain plasticity: expressive suppression use predicts anterior insula volume. Neuroimage 58, 10–15. doi: 10.1016/j.neuroimage.2011.06.028
Goldin, P. R., McRae, K., Ramel, W., and Gross, J. J. (2008). The neural bases of emotion regulation: reappraisal and suppression of negative emotion. Biol. Psychiatry 63, 577–586. doi: 10.1016/j.biopsych.2007.05.031
Gross, J. J., and John, O. P. (2003). Individual differences in two emotion regulation processes: implications for affect, relationships, and well-being. J. Pers. Soc. Psychol. 85, 348–362. doi: 10.1037/0022-35184.108.40.2068
Hamilton, J. P., Farmer, M., Fogelman, P., and Gotlib, I. H. (2015). Depressive rumination, the default-mode network and the dark matter of clinical neuroscience. Biol. Psychiatry 78, 224–230. doi: 10.1016/j.biopsych.2015.02.020
Hermann, A., Bieber, A., Keck, T., Vaitl, D., and Stark, R. (2014). Brain structural basis of cognitive reappraisal and expressive suppression. Soc. Cogn. Affect. Neurosci. 9, 1435–1442. doi: 10.1093/scan/nst130
Hu, L., and Bentler, P. M. (1998). Fit indices in covariance structure modeling: sensitivity to under parameterized model misspecification. Psychol. Methods 3, 424–453. doi: 10.1037//1082-989x.3.4.424
Jäncke, L., and Langer, N. (2011). A strong parietal hub in the small-world network of coloured-hearing synaesthetes during resting state EEG: small-world network of coloured-hearing synaesthetes. J. Neuropsychol. 5, 178–202. doi: 10.1111/j.1748-6653.2011.02004.x
Kerns, J. G., Cohen, J. D., MacDonald, A. W. III., Cho, R. Y., Stenger, V. A., and Carter, C. S. (2004). Anterior cingulate conflict monitoring and adjustments in control. Science 303, 1023–1026. doi: 10.1126/science.1089910
Langer, N., Pedroni, A., Gianotti, L. R. R., Hänggi, J., Knoch, D., and Jäncke, L. (2012). Functional brain network efficiency predicts intelligence. Hum. Brain Mapp. 33, 1393–1406. doi: 10.1002/hbm.21297
Lee, T.-W., Dolan, R. J., and Critchley, H. D. (2008). Controlling emotional expression: behavioral and neural correlates of nonimitative emotional responses. Cereb. Cortex 18, 104–113. doi: 10.1093/cercor/bhm035
Long, J., Xie, Q., Ma, Q., Urbin, M. A., Liu, L., Weng, L., et al. (2016). Distinct interactions between fronto-parietal and default mode networks in impaired consciousness. Sci. Rep. 6:38866. doi: 10.1038/srep38866
MacDonald, A. W. III., Cohen, J. D., Stenger, V. A., and Carter, C. S. (2000). Dissociating the role of the dorsolateral prefrontal and anterior cingulate cortex in cognitive control. Science 288, 1835–1838. doi: 10.1126/science.288.5472.1835
Mears, D., and Pollard, H. B. (2016). Network science and the human brain: using graph theory to understand the brain and one of its hubs, the amygdala, in health and disease. J. Neurosci. Res. 94, 590–605. doi: 10.1002/jnr.23705
Meng, C., Brandl, F., Tahmasian, M., Shao, J., Manoliu, A., Scherr, M., et al. (2014). Aberrant topology of striatum’s connectivity is associated with the number of episodes in depression. Brain 137, 598–609. doi: 10.1093/brain/awt290
Moore, M., Iordan, A. D., Hu, Y., Kragel, J. E., Dolcos, S., and Dolcos, F. (2016). Localized or diffuse: the link between prefrontal cortex volume and cognitive reappraisal. Soc. Cogn. Affect. Neurosci. 11, 1317–1325. doi: 10.1093/scan/nsw043
Murphy, K., Birn, R. M., Handwerker, D. A., Jones, T. B., and Bandettini, P. A. (2009). The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage 44, 893–905. doi: 10.1016/j.neuroimage.2008.09.036
Picó-Pérez, M., Alonso, P., Contreras-Rodríguez, O., Martínez-Zalacaín, I., López-Solà, C., Jiménez-Murcia, S., et al. (2017). Dispositional use of emotion regulation strategies and resting-state cortico-limbic functional connectivity. Brain Imaging Behav. doi: 10.1007/s11682-017-9762-3 [Epub ahead of print].
Power, J. D., Cohen, A. L., Nelson, S. M., Wig, G. S., Barnes, K. A., Church, J. A., et al. (2011). Functional network organization of the human brain. Neuron 72, 665–678. doi: 10.1016/j.neuron.2011.09.006
Santarnecchi, E., Galli, G., Polizzotto, N. R., Rossi, A., and Rossi, S. (2014). Efficiency of weak brain connections support general cognitive functioning. Hum. Brain Mapp. 35, 4566–4582. doi: 10.1002/hbm.22495
Sheffield, J. M., Repovs, G., Harms, M. P., Carter, C. S., Gold, J. M., MacDonald III, A. W., et al. (2015). Fronto-parietal and cingulo-opercular network integrity and cognition in health and schizophrenia. Neuropsychologia 73, 82–93. doi: 10.1016/j.neuropsychologia.2015.05.006
Smith, S. M., Vidaurre, D., Beckmann, C. F., Glasser, M. F., Jenkinson, M., Miller, K. L., et al. (2013). Functional connectomics from resting-state fMRI. Trends Cogn. Sci. 17, 666–682. doi: 10.1016/j.tics.2013.09.016
Stam, C. J., Jones, B. F., Nolte, G., Breakspear, M., and Scheltens, P. (2006). Small-world networks and functional connectivity in Alzheimer’s disease. Cereb. Cortex 17, 92–99. doi: 10.1093/cercor/bhj127
Stanley, M. L., Simpson, S. L., Dagenbach, D., Lyday, R. G., Burdette, J. H., and Laurienti, P. J. (2015). Changes in brain network efficiency and working memory performance in aging. PLoS One 10:e0123950. doi: 10.1371/journal.pone.0123950
Stawarczyk, D., and D’Argembeau, A. (2015). Neural correlates of personal goal processing during episodic future thinking and mind-wandering: an ALE meta-analysis. Hum. Brain Mapp. 36, 2928–2947. doi: 10.1002/hbm.22818
Thompson, W. H., and Fransson, P. (2015). The frequency dimension of fMRI dynamic connectivity: network connectivity, functional hubs and integration in the resting brain. Neuroimage 121, 227–242. doi: 10.1016/j.neuroimage.2015.07.022
van den Heuvel, M. P., and Hulshoff Pol, H. E. (2010). Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 8, 519–534. doi: 10.1016/j.euroneuro.2010.03.008
van den Heuvel, M. P., Stam, C. J., Kahn, R. S., and Hulshoff Pol, H. E. (2009). Efficiency of functional brain networks and intellectual performance. J. Neurosci. 29, 7619–7624. doi: 10.1523/jneurosci.1443-09.2009
Vatansever, D., Menon, D. K., Manktelow, A. E., Sahakian, B. J., and Stamatakis, E. A. (2015). Default mode dynamics for global functional integration. J. Neurosci. 35, 15254–15262. doi: 10.1523/jneurosci.2135-15.2015
Vincent, J. L., Kahn, I., Snyder, A. Z., Raichle, M. E., and Buckner, R. L. (2008). Evidence for a frontoparietal control system revealed by intrinsic functional connectivity. J. Neurophysiol. 100, 3328–3342. doi: 10.1152/jn.90355.2008
Wang, K., Huang, H., Chen, L., Hou, X., Zhang, Y., Yang, J., et al. (2017). MRI correlates of interaction between gender and expressive suppression among the Chinese population. Neuroscience 347, 76–84. doi: 10.1016/j.neuroscience.2017.01.042
Wang, J., Wang, X., Xia, M., Liao, X., Evans, A., and He, Y. (2015). Corrigendum: GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front. Hum. Neurosci. 9:458. doi: 10.3389/fnhum.2015.00458
Wang, L., Zhu, C., He, Y., Zang, Y., Cao, Q., Zhang, H., et al. (2009). Altered small-world brain functional networks in children with attention-deficit/hyperactivity disorder. Hum. Brain Mapp. 30, 638–649. doi: 10.1002/hbm.20530
Watson, D., Clark, L., and Tellegen, A. (1988). Development and validation of brief measures of positive and negative affect: the PANAS scales. J. Pers. Soc. Psychol. 54, 1063—1070. doi: 10.1037//0022-35220.127.116.113
Welborn, B. L., Papademetris, X., Reis, D. L., Rajeevan, N., Bloise, S. M., and Gray, J. R. (2009). Variation in orbitofrontal cortex volume: relation to sex, emotion regulation and affect. Soc. Cogn. Affect. Neurosci. 4, 328–339. doi: 10.1093/scan/nsp028
Xu, J., Yin, X., Ge, H., Han, Y., Pang, Z., Tang, Y., et al. (2015). Attentional performance is correlated with the local regional efficiency of intrinsic brain networks. Front. Behav. Neurosci. 9:200. doi: 10.3389/fnbeh.2015.00200
Keywords: DMN, emotion regulation, FPN, graph theory, resting-state fMRI
Citation: Pan J, Zhan L, Hu C, Yang J, Wang C, Gu L, Zhong S, Huang Y, Wu Q, Xie X, Chen Q, Zhou H, Huang M and Wu X (2018) Emotion Regulation and Complex Brain Networks: Association Between Expressive Suppression and Efficiency in the Fronto-Parietal Network and Default-Mode Network. Front. Hum. Neurosci. 12:70. doi: 10.3389/fnhum.2018.00070
Received: 05 December 2017; Accepted: 07 February 2018;
Published: 16 March 2018.
Edited by:Delin Sun, Duke University, United States
Copyright © 2018 Pan, Zhan, Hu, Yang, Wang, Gu, Zhong, Huang, Wu, Xie, Chen, Zhou, Huang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
† These authors have contributed equally to this work.