ORIGINAL RESEARCH article

Front. Neurosci., 25 October 2019

Sec. Brain Imaging Methods

Volume 13 - 2019 | https://doi.org/10.3389/fnins.2019.01069

Negatively Linking Connector Networks in Cognitive Control of Affective Pictures

  • 1. College of Education Science, The Big Data Centre for Cognitive Neuroscience and AI, Hengyang Normal University, Hengyang, China

  • 2. Research Center of Brain and Cognitive Neuroscience, Liaoning Normal University, Dalian, China

  • 3. Institute for Brain and Psychological Sciences, Sichuan Normal University, Chengdu, China

  • 4. School of Education, Anyang Normal University, Anyang, China

  • 5. College of Psychology and Sociology, Shenzhen University, Shenzhen, China

Abstract

Cognitive control of emotions depends on intermodular long-distance communications. However, negative connections between connector hubs are removed by traditional hard-thresholding approach in graph-theoretical research. Using soft-thresholding approach to reserve negative links, we explore time-varying features of connector hubs in intermodular communications during cognitive control of affective pictures. We develop a novel approach to sparse functional networks and construct negatively linking connector networks for positive, negative, and neutral pictures. We find that consisting of flexible hubs, the frontoparietal system dynamically top–down inhibits neural activities through negative connections from the salience subnetwork and visual processing area. Moreover, the shared connectors form functional backbones that dynamically reconfigure according to differently-valenced pictures in order to coordinate both stability and flexibility of cognitive connector networks. These results reveal the necessity of conserving negative links between intermodular communications in chronnectome research and deepen the understanding of how connector networks dynamically evolute during cognitive control of affective processing.

Introduction

The ability to cognitive control of emotions is critical to behavioral flexibility and well-being (Ochsner et al., 2012; Marek et al., 2015). Although patterns of activation during cognitive control of emotions have been well-characterized, graph theory-based connectivity research has shown that the brain network composes of functionally separate subnetworks or modules (Buckner et al., 2009; Barret and Satpute, 2013; Sripada et al., 2014). Furthermore, functional connectome research has revealed that affective processing at least contains frontoparietal subnetwork (e.g., dorsolateral prefrontal cortex, inferior and superior parietal lobule), salience subnetwork (e.g., insula, middle, and posterior cingulate cortex), visual subnetwork (within the occipital lobe), and cerebellum (Wolf et al., 2009; Lindquist and Barrett, 2012; Oosterwijk et al., 2012). However, with the emergence of chronnectome (Calhoun et al., 2014), little is known about how affective information is dynamically represented through distributed large-scale brain networks.

During affective processing, cognitive control relies on not only short-distance intramodular processing but also long-distance intermodular communications (Mesulam, 1990; Wager et al., 2008; Barret and Satpute, 2013). As a matter of fact, whole-brain networks contain some important long-distance communications (e.g., long intra/inter-hemispheric connections) (Sepulcre et al., 2010; Hermundstad et al., 2013). In graph-theoretic research, hubs integrate and distribute information in powerful ways and are further classifies into provincial and connector hubs (Power et al., 2013; Van den Heuvel and Sporns, 2013). Different from provincial hubs that connect nodes within a module, connector hubs with high participant coefficients connect to multiple communities and have access to a variety of types of different information (Power et al., 2011). Particularly, the focal damage to a connector hub (e.g., a frontoparietal hub) caused large-scale long-distance functional deficits (Gratton et al., 2012). Therefore, focusing on the role of connector hubs is crucial for understanding of intermodular communication during cognitive control of affective pictures. However, static connectomics did not capture temporal structure of brain activity (Calhoun et al., 2014). Recently, more and more research has become interested in understanding network dynamics that support emotional or cognitive states (Fornito et al., 2012; Shine et al., 2016; Betzel et al., 2017), learning (Mantzaris et al., 2013), development (Gu et al., 2015), aging (Betzel et al., 2014), and disease progression (Raj et al., 2015). Time-vary features of functional network (e.g., connector hubs) are import for understanding how functional network dynamically evolute during information communication (Sizemore and Bassett, 2018), which further deepens our understanding the nature of temporal dynamics of brain activity (Vidaurre et al., 2018).

To obtain a sparse adjacent matrix in graph theoretic analysis, the traditional hard-thresholding approach eliminates weak and negative couplings between nodes because these negative links are often conceived as noise or minor (Rubinov and Sporns, 2010, 2011). However, not all negative connections are biologically meaningless (Chang and Glover, 2009; Fornito et al., 2013). For instance, during cognitive reappraisal of affective pictures, there exist strong anticorrelations between prefrontal cortex (PFC) (e.g., dorsomedial PFC or dorsolateral PFC) and salience/limbic areas (e.g., insula or amygdala) (Ochsner et al., 2002; Urry et al., 2006; Drabant et al., 2009). Undeniably, negative connections between modules are quite pervasive in functional network architectures and play an important role in behavior (Fox et al., 2005; De pisapia et al., 2012; Fornito et al., 2013), particularly in the task-evoked networks, although they were rarely mentioned. Simply removing these negative links lead to inaccurate representation of affective networks. The reservation of internetwork anti-correlations are critical for the precise representation of connector hubs in graph theoretical analysis because connector hubs are responsible for intermodular communication (Buckner et al., 2009; Van den Heuvel and Sporns, 2011). In contrast, the soft-thresholding approaches to mapping [−1,1] to [0, 1] not only avoid issues related to network fragmentation but also simply suppress rather than remove weaker connections (Schwarz and McGonigle, 2011). Here we aimed to use soft-thresholding approach to explore time-varying features of connector hubs in intermodular communication during affective information processing.

This study aimed to explore how connector networks dynamically evolute during cognitive control of affective processing. We collected fMRI data from 26 right-handed healthy undergraduates while they performed a cognitive task (silently counting backward) during the viewing of affective (positive, negative, and neutral) pictures. After preprocessing, time-series of 278 whole-brain parcels were extracted and subjected to a wavelet decomposition (Shen et al., 2013). Then, we estimated the dynamic correlation coefficients (DCCs) of all possible pairs of 278 parcels (Lindquist et al., 2014) and obtained 27 group-averaged DCC matrices. These DCC matrices were subsequently Fisher z-transformed and subjected to soft-thresholding. Next, we used the BCT toolbox1 to calculate graph theoretic features of these full-connected matrices. Finally, we constructed and characterized connector networks for affective pictures across 9 time points. Based on the previous study of dynamic networks (Fornito et al., 2012; Shine et al., 2016; Betzel et al., 2017; Sizemore and Bassett, 2018; Vidaurre et al., 2018), we expected that connector networks would exhibit valence-dependent dynamic connectivity patterns over times; in order to coordinate the stability and the flexibility of connector networks during dynamic affective processing, we expected that the shared connector would keep robust connections and form the functional backbone correspondent with structural backbone (Hagmann et al., 2008; Liang et al., 2013; Zhang et al., 2015).

Materials and Methods

Participants

The participants consisted of 26 right-hand undergraduates (14 females, 12 males; mean age 21.25 ± 2.34 years). The participants self-reported having normal or corrected-to-normal visual acuity and had no prior history of psychiatric or neurological disorders. Data from one male and one female were ruled out because of poor-quality data collection. All the participants provided the written informed consent before the experiment and were paid for their participation after the experiment. All experiments were performed in accordance with relevant guidelines and regulations. This study was approved by the Ethical Committee of School of Psychology at East China Normal University.

Materials

Stimuli consisted of 81 pictures from the Chinese Affective Picture System (Bai et al., 2005), which is a collection of standardized photographic materials from the International Affective Picture System. The images were grouped into three conditions: 27 positive pictures (e.g., a smiling face), 27 negative pictures (e.g., a wreckage), and 27 neutral pictures (e.g., a houshold object). There exist significant valence differences among the three picture categories [F(2,78) = 90.25, p < 0.001; M ± SD: Positive = 7.35 ± 0.14; Neutral = 4.79 ± 0.07; Negative = 2.41 ± 0.12]. The positive and negative pictures significantly differ from the neutral ones [F(1,75) = 54.62, 59.21, p < 0.001] but do not significantly differ from each other [F(1,75) = 1.28, p > 0.05; M ± SD: Positive = 5.69 ± 0.39; Neutral = 4.45 ± 0.41; Negative = 5.58 ± 0.36]. We used E-prime1.1 software (Psychology Software Tools, Pittsburgh, PA, United States) to present these images and control their timing.

Task

On full trials (see Figure 1), a cue of the Chinese word “counting silently” under a fixation cross was first presented for 2 s. Then, three images with the same valence were continuously presented for 18 s, 6 s per picture. Once the cue emerged on the screen, the participants were required to count backward from 100 by three’s, i.e., 100, 97, 94, …, until all the pictures vanished. During counting silently, the participants always viewed the pictures. On the next screen the Chinese word “Assessment” appeared for 2 s, which prompted the participants to rate the arousal level of the pictures on a nine-point scale (1 = extremely weak to 9 = extremely strong) by pressing a button. Following the rating, a 4-, 6-, or 8-s jittered blank screen concluded the trial. The experiment contained 27 trials, 9 trials per picture type. The trial order was randomized and counterbalanced across subjects. After the fMRI experiment completed, the participants were asked to offline rate the arousal level using the similar paradigm with simply viewing 81 pictures outside the scanner. Before entering the scanner, the participants were familiar with the experimental procedure.

FIGURE 1

fMRI Acquisition

Structural and functional MRI data were obtained with a 3-T Trio Tim Magnetic Resonance Imaging scanner (Siemens company) at East China Normal University. Foam padding was used to minimize head movement. Whole-brain functional data were acquired using a BOLD contrast-sensitive echo-planar T2-weighted imaging (EPI) sequence (32 axial slices, slice thickness = 5 mm; TR = 2000 ms, TE = 30 ms, flip angle = 90°, field of view = 240 mm × 240 mm, matrix size = 64 × 64, voxel size = 3.75 mm × 3.75 mm × 5 mm). The first two functional volumes were discarded to allow for equilibration effects. A T1-weighted anatomical image was obtained for each participants with 1 mm × 1 mm × 1 mm voxels (TR = 1900 ms, TE = 3.43 ms, flip angle = 7°).

fMRI Analysis

The preprocessing of functional data was carried out using the SPM8 software package2. After functional data were corrected for slice timing and head motion, the T1-weighted images were coregistered to the mean EPI image. The coregistered data were spatially normalized into the standard MNI space and resampled to 3 mm isotropic voxels. Images were high-filtered with a cut-off of 128 s to remove low-frequency signal drift.

Network Construction

The whole brain was parcellated into a set of subunits with 2 mm resolution that correspond to the 278 cortical and subcortical structures (Shen et al., 2013). These subunits were obtained using resting-state fMRI and have been shown to be highly functionally homogeneous and spatially coherent. For each individual fMRI dataset, regional mean BOLD time-series of each subunit were extracted using the marsbar toolbox3 separately for the positive, negative, and neutral picture conditions. These regional time-series were then subjected to a wavelet decomposition to reconstruct wavelet coefficients in the 0.03–0.07 Hz range using the wmtsa toolbox4. Unlike the traditional assumption that the functional connectivity between time series from distinct brain regions is constant across time, the DCC provides a model-based approach toward estimating dynamic correlations and achieves the best overall balance between sensitivity and specificity in detecting dynamic changes in correlations (Lindquist et al., 2014). Therefore, we used the DCC toolbox (Lindquist et al., 2014) to estimate DCCs of these wavelet coefficients time-series between all possible pairs of 278 brain subunits and averaged the 278 × 278 functional connectivity matrices over all trials and all participants. The resulting 27 group-averaged DCC matrices (9 TRs × 3 picture conditions) were Fisher z-transformed.

Graph-Theoretic Analysis

To conserve negative and weak edges, we used the soft-thresholding approach to map a correlation coefficients rij[−1, 1] to a weight wij[0, 1]. These mappings have recently been evaluated in a study of the modular organization of the brain from functional connectivity data (Mumford et al., 2010). This approach can replace the hard thresholding with a continuous mapping of all correlation values to edge weights, suppressing rather than removing weaker connections and avoiding issues related to network fragmentation. Here, a soft-thresholding adjacency function was defined by a power function (Schwarz and McGonigle, 2011), wij = [(rij + 1)/2]12. We used the BCT toolbox to compute the following four graph-theoretic characteristics of parcel-wise whole-brain networks. The first analysis used Louvain algorithm to identify community structures of full-connection networks of 278 nodes at each time point. Then, we computed the participation coefficients to examine the role of hubs in interconnecting distinct modules. We also calculated nodal betweenness centrality and edge betweenness centrality. In this study, the hubs with participation coefficient > 0.5 and both node betweenness centrality and edge betweenness centrality ≠ 0 at each time point were classified into connector hubs. These connector hubs describe the level of intermodule connectivity and facilitate the integration of multiple types of information during intermodular communication. To sparse the full-connected DCC matrices, we binarized the matrices of edge betweenness centrality and then dot-multiplied them by the correspondent DCC matrices. That is, we retained the connectivity values between two connectors in the DDC matrices only if their edge betweenness centrality value was not equal to zero. Thus, we obtained the connector networks.

Results

Behavioral Data

In order to verify if silently counting validly reduce the emotional responses evoked by affective picture, a two-way repeated measure ANOVA with valence (positive, neutral, and negative) × task (silently counting and simple viewing) revealed a significant main effect of valence on arousal rating [F(2,46) = 10.25, = 0.39, p < 0.001] and a significant interaction between valence and task [F(2,46) = 9.43, = 0.32, p < 0.001]. Post hoc tests indicated that the counting task produced a lower arousal rating for positive and negative pictures relative to simple viewing (p < 0.01), but not for neutral pictures (p > 0.05) (see Figure 2). These results suggest that participants were involved in the cognitive control of affective picture processing within the scanner.

FIGURE 2

Identification of Connector Hubs

We calculated several graph-theoretic features to examine the role of hubs in interconnecting distinct modules. Participant coefficient is a measure of diversity of intermodular connections of individual nodes (Van den Heuvel and Sporns, 2011). Nodes with high participant coefficient display high correlation with multiple modules, which facilitates integration of multiple types of information (Warren et al., 2014). Nodal betweenness centrality measures the influence of a node on information flow between additional nodes in the network (Hagmann et al., 2008). Similarly, edge betweenness centrality measures the influence of an edge on information flow (Brandes, 2001). At all 9 time points, nodes with participation coefficient > 0.5 and both nodal betweenness centrality and edge betweenness centrality ≠ 0 were classified as connector hubs. We retained the connector hubs that have at least a connection with other connector hubs in different communities. Finally, we binarized edge betweenness centrality matrices and dot-multiplied them by the correspondent DCC matrices. Thus, we sparsed the full-connected DCC matrices into connector networks (Supplementary Figure S1).

Figure 3 shows 65 connector hubs for positive, negative, and neutral pictures for their coordinates. The connector networks contained 31 nodes for positive pictures (see Table 1), mainly including the frontoparietal nodes (e.g., inferior and middle frontal gyrus, inferior and superior parietal lobule), the salience nodes (e.g., insula, putamen, caudate, amygdala), and the visual nodes (e.g., cuneus, middle and superior occipital gyrus). For negative pictures, the connector network contained 42 nodes (see Table 2) and their mayjority were located in the frontoparietal area (e.g., inferior and superior frontal gyrus, inferior and superior parietal lobe, precentral gyrus, angular gyrus), the salience system (e.g., insula, putamen, cingulate cortex, hippocampus), and the cerebellum. For neutral pictures, the connector network included 23 nodes (see Table 3) and were mainly found in the frontoparietal area (inferior frontal gyrus, superior parietal lobe, post-central gyrus), the salience system (insula, putamen, posterior cingulate cortex), and the visual area (cuneus, superior occipital gyrus). Moreover, all the connector networks were always composed of seven modules (see Supplementary Figure S2) and each module possessed more than one connectors.

FIGURE 3

TABLE 1

NodeSubunitBAAALTalairach

XYZ
1115Pons_R7−23−26
2133R.BA53.1Amygdala_R25−2−16
3149L.BA49.1Putamen_L−2342
4154L.BA4.1Pretcentral_L−57−928
5159Culmen_L−44−53−36
6165L.BA18.8Cuneus_L−4−9015
7169L.BA45.2Insula_L−32224
8174Cerebellum_8_L−20−65−40
9177L.BA6.4Inferior Frontal Gyrus−41931
10190L.BA7.6Occipital_Sup_L−14−6829
11195L.BA40.1Inferior Parietal Lobule−57−2833
12201L.BA37.1Temporal_Mid_L−57−48−9
13203L.BA6.6Middle Frontal Gyrus−241151
14211Cerebellar Tonsil−27−49−41
15216L.BA13.2Insula_L−367−1
16223L.BA19.9Cuneus_L−10−8332
17240L.BA7.4Parietal_Sup_L−21−5657
18244L.BA38.1Temporal_Sup_L−460−12
19245L.BA40.3Inferior Parietal Lobule−55−3823
20249L.BA7.9Precuneus_L−3−6738
21255L.BA7.7Precuneus_L−12−3964
22260L.BA48.1Caudate_L−14−1520
23269L.BA44.2Insula_L−35511
24277L.BA6.7Supp_Motor_Area_L−14−263
2529R.BA7.6Superior Parietal Lobule24−4661
2639R.BA19.5Occipital_Mid_R43−7327
2745R.BA7.3Paracentral_Lobule_R7−4062
2860R.BA7.5Superior Parietal Lobule23−6649
2966R.BA20.3Cerebelum_Crus1_R46−53−33
3080R.BA7.4Parietal_Sup_R21−5655
3196R.BA18.5Occipital_Inf_R23−90−3

Connector hubs for positive pictures.

The first column denotes the nodal number in the Pajek graphs. The second column represents the subunit codes of whole-brain parcellation by Shen et al. (2013). BA, Broadmann’s area; AAL, automated anatomical labeling.

TABLE 2

NodeSubunitBAAALTalairach

XYZ
1111R.BA44.2Frontal_Inf_Oper_R521312
2117R.BA19.1Lingual_R14−580
3123R.BA38.2Temporal_Pole_Sup_R3712−24
4134R.BA47.2Frontal_Inf_Orb_R2633−13
5136R.BA36.2Precuneus_R15−433
6149L.BA49.1Putamen_L−2342
7154L.BA4.1Pretcentral_L−57−928
8158Pons−4−21−27
9159Culmen_L−44−53−36
‘10162L.BA37.3Cerebellum_Crus2_L−46−63−23
11165L.BA18.8Cuneus_L−4−9015
12169L.BA45.2Insula_L−32224
13171L.BA23.1Posterior Cingulate−8−5621
14174Cerebelum_8_L−20−65−40
15177L.BA6.4Inferior Frontal Gyrus−41931
16190L.BA7.6Occipital_Sup_L−14−6829
17194L.BA39.5Angular_L−39−6144
18195L.BA40.1Inferior Parietal Lobule−57−2833
19200L.BA10.4Frontal_Sup_L (aal)−16603
20202L.BA19.6Cerebellum_6_L−26−67−20
21206L.BA6.1Inferior Frontal Gyrus−53123
22207L.BA13.3Insula_L−35−120
23216L.BA13.2Insula_L−367−1
24223L.BA19.9Cuneus_L−10−8332
25225L.BA10.7Cingulum_Ant_L−7470
26238L.BA21.3Temporal_Mid_L−56−29−13
27240L.BA7.4Parietal_Sup_L−21−5657
28244L.BA38.1Temporal_Sup_L−460−12
2924R.BA39.2Temporal_Sup_R48−5934
30251L.BA8.2Frontal_Sup_L−212540
31255L.BA7.7Precuneus_L−12−3964
32267L.BA10.6Frontal_Sup_Medial_L−105720
33269L.BA44.2Insula_L−35511
34273L.BA48.2Caudate_L−12713
3535R.BA37.8Temporal_Mid_R54−560
3642R.BA54.2Hippocampus_R30−360
3755R.BA31.2Cingulum_Mid_R9−2845
3875R.BA39.4Inferior Parietal Lobule58−4724
3981R.BA10.1Cingulum_Ant_R8422
4086R.BA37.7Cerebellum_Crus1_R47−72−24
4193R.BA13.3Insula_R40−614
4294R.BA40.3Parietal_Inf_R51−4640

Connector hubs for negative pictures.

TABLE 3

NodeSubunitBAAALTalairach

XYZ
1136R.BA36.2Precuneus_R15−433
2149L.BA49.1Putamen_L−2342
3154L.BA4.1Pretcentral_L−57−928
4165L.BA18.8Cuneus_L−4−9015
5169L.BA45.2Insula_L−32224
6171L.BA23.1Posterior Cingulate−8−5621
7174Cerebellum_8_L−20−65−40
8190L.BA7.6Occipital_Sup_L−14−6829
9198L.BA20.1Temporal_Inf_L−49−11−26
10206L.BA6.1Inferior Frontal Gyrus−53123
11216L.BA13.2Insula_L−367−1
12223L.BA19.9Cuneus_L−10−8332
13224L.BA31.1Precuneus_L−8−4047
14240L.BA7.4Parietal_Sup_L−21−5657
15244L.BA38.1Temporal_Sup_L−460−12
16252L.BA37.6Fusiform_L−32−46−24
17263L.BA48.3Caudate_L−11182
18266L.BA54.1−28−20−9
19269L.BA44.2Insula_L−35511
2045R.BA7.3Paracentral_Lobule_R7−4062
215R.BA40.4Postcentral_R57−2031
2280R.BA7.4Parietal_Sup_R21−5655
2390R.BA37.9Temporal_Inf_R54−48−14

Connector hubs for neutral pictures.

Time-Varying Connector Networks in Affective Processing

Across 9 time points (from TR = 1 to 9), the connector networks displayed different dynamic connectivity patterns when participants performed cognitive processing of positive, negative, and neutral pictures (see Figures 46, respectively). The DCC matrices of connector hubs exhibited time-varying patterns (see Supplementary Video S1 for positive, Supplementary Video S2 for negative, and Supplementary Video S3 for neutral online). However, regardless of picture valences, the connection values between connectors were always negative when the connectors belongs to different communities. Their affiliative communities also reconfigured dynamically over times (Figures 7A,E,I).

FIGURE 4

FIGURE 5

FIGURE 6

FIGURE 7

Next, we used the Gretna software5 to calculate the small-worldness value of connector networks. Small-worldness is a property of a network with high clustering but low characteristic path length (Yan et al., 2011). A real network is considered small-world if the ratio of cluster coefficient to characteristic path length is greater than one relative to random networks. Previous research has demonstrated that the small-world feature are quite universal in structural and functional networks (Varoquaux and Craddock, 2013). However, small-worldness in this study was less than one across nine time points (all ps < 0.01) because standardized cluster coefficients (gamma value) were significantly less than one (all ps < 0.01) but standardized characteristic path length (lambda value) were almost close to one. Thus, small-worldness did not emerge in connector networks (Figure 8).

FIGURE 8

Positive Pictures

In the frontoparietal system, the left superior parietal lobule (Node 17, BA7.4) had larger betweenness centrality value than other 12 connector hubs (p < 0.001, Figure 7B). In the salience system, betweenness centrality value of the left insula (Node 23, BA44.2) was larger than other five salience hubs (p < 0.01, Figure 7C). In the visual system, the left cuneus (Node 6, BA18.8) exhibited larger betweenness centrality value than other four hubs (p < 0.01, Figure 7D).

Negative Pictures

In the frontoparietal system, the superior medial frontal gyrus (Node 32, BA10.6) in the left hemisphere had larger betweenness centrality value than other 13 connector hubs (p < 0.001, Figure 7F). Of 13 hubs in the salience system, the left insula (Node 33, BA44.2) had the largest betweenness centrality value (p < 0.001, Figure 7G). The cerebellum had five connector hubs that involve negative processing, of which Node 20 exhibited the largest betweenness centrality value (p < 0.01, Figure 7H).

Neutral Pictures

In the frontoparietal system, the betweenness centrality value of the left superior parietal lobe (Node 14, BA7.4) was larger than other seven connector hubs (p < 0.001, Figure 7J). The left insula (Node 19, BA44.2) exhibited the largest betweenness centrality value in the salience sytem (p < 0.01, Figure 7K). In the visual system, the left cuneus had larger betweenness centrality value than other two hubs (p < 0.001, Figure 7L).

Shared Hubs in the Functional Backbone

The connector networks for positive, negative, and neutral pictures always shared 11 hubs in the left hemisphere (Figure 9). They were involved in affective processing (e.g., insula, putamen, and cerebellum), visual processing (e.g., cuneus and superior occipital gyrus), and cognitive processing (e.g., superior parietal lobule, precentral gyrus, and superior temporal gyrus). After binarizing edge betweenness centrality matrices composed of 11 nodes, these matrices always kept invariant functional connectivity patterns on the same picture condition at all 9 time points (solid lines in Figures 46). Moreover, the shared 11 connectors always kept relatively robust during full trials. These shared connectors kept largely correspondent to structural cores with high metabolic activity in the structural backbone. Somewhat similar to a structural backbone in structural connectivity networks (Hagmann et al., 2008), the shared connector networks were tentatively called functional backbone (see secion “Discussion”), which is necessary for coordinating both the stability and the flexibility of cognitive control networks during dynamic affective processing.

FIGURE 9

Discussion

The current study used graph theory approach to explore how connector hubs dynamically evolute during cognitive control of affective pictures. We used soft-thresholding to reserve negative connections (Schwarz and McGonigle, 2011). Results revealed that connector networks consisted of cross-modular negative connectivities, which is impossible to be unveiled by hard-thresholding. This might be because we were the first time to use edge betweenness centrality to sparse connector networks, different from traditional hard-thresholding that just removes weak links (Rubinov and Sporns, 2010, 2011). Edge betweenness centrality is the fraction of all shortest paths in the network that contain a given edge and reflects an edge’ centrality of in a network but not link strength (Brandes, 2001). Thus, negative links between connector hubs is not noise or spurious connections but take an important role in intermodular long-distance communications. Overall, this study not only provides a novel method to sparse functional networks but also reveals that there exist negatively linking connector networks during intermodular integration of affective information.

Connector networks exhibited valence-related dynamic connectivity patterns over times, e.g., dynamic community affiliations and time-varying connections with other parts of whole-brain networks. These dynamic changes facilitate intermodular communications of affective information. However, connector networks did not exhibit small-worldness because of low cluster coefficient but high characteristic path length, regardless of picture valences. This implies that these connector networks contain less direct links but more indirect links. For example, the frontoparietal subnetwork does not receive direct afferents from the visual system, but executes top–down control indirectly via the ventral attention/salience network, which contains the anterior insula (Seeley et al., 2007; Vossel et al., 2012; Palaniyappan et al., 2013). This study requested participants to perform the silent-counting task while they viewed pictures on the screen. Participants needed to inhibit visually affective information when pictures appeared on the screen. Behavioral results indicated that the counting task significantly decreased arousal rating for positive and negative pictures but not neutral pictures compared with offline viewing. Consistent with behavioral findings, the frontoparietal system (e.g., superior frontal gyrus and superior parietal lobe) through negative connections dynamically top–down inhibited the neural activities from salience subnetwork (e.g., insula) and visual area (e.g., cuneus).

Affective picture processing also exhibited different nodal betweenness patterns in connector networks. The present study found that the left superior parietal lobe (BA7.4) in the frontoparietal system exhibited larger nodal betweenness for positive and neutral pictures, although connector networks contained less hubs in the neutral picture condition than in the positive picture condition. In contrast, connector networks of negative pictures included larger nodal betweenness in the superior medial frontal gyrus (BA10.6). These results were consistent with the frontoparietal control system consisting of flexible hubs that integrate information from the external environment (e.g., visual stimuli) with stored internal representations (e.g., silently counting) and regulate distributed systems (e.g., visual, salience) according to affective valences (Vincent et al., 2008; Woolgar et al., 2015). Moreover, the left insula (BA44.2) had larger nodal betweenness centrality in the salience system of all connector networks. This might be attributed to that the insula regulates the interaction between the salience of the selective attention created to achieve a task and the salience of emotional arousal and has an important role in the experience of emotions (Taylor et al., 2009; Quarto et al., 2016).

In particular, five cerebellar hubs participated in network processing of negative pictures while fewer cerebellar connector hubs for positive and neutral pictures. As being famous for sensorimotor control (Bastian, 2011), the human cerebellum has connections with other brain areas involved in affective regulation, including the insula, amygdala, hypothalamus, as well as the neocortex (Baumann and Mattingley, 2012). Patients with cerebellar diseases had significant higher depressive scores relative to healthy comparison groups (Wolf et al., 2009). Neuroimaging studies revealed that the cerebellum was involved in control processing of aversive stimuli and other basic emotions (Moulton et al., 2011). However, cerebellar hubs were ignored by several cortical parcellations (e.g., Power et al., 2011; Yeo et al., 2011; Gordon et al., 2014). Future studies should follow in interest the centrality of cerebellar hubs in whole-brain network processing of affective information.

The shared 11 connectors kept robust connections and formed the functional backbone. The functional backbone mainly consisted of connector hubs with high metabolic activity and might be constrained by cortical structural cores (Hagmann et al., 2008). Different from structural connections, functional connectivities reflect statistical dependence between nodes and thus functional hubs have unnecessarily direct structural connection with other hubs (Sporns, 2014). Under the disturbing of active tasks, metabolic energy is necessarily redistributed more to task-specific hubs to complete tasks (Liang et al., 2013; Zhang et al., 2015). The present functional backbone exhibited different connectivity patterns among positive, negative, and neutral pictures, which reconfigured during cognitive control of differently-valenced affective processing.

From a perspective of temporal dynamics, the existence of the functional backbone helps not only maintain network integrity (or task set) across full trials but also flexibly integrate different types of information from different modules during cognitive processing. Therefore, the functional backbone plays a critical role in coordinating both the stability and the flexibility of cognitive control networks during dynamic processing of affective information. In particular, the connector lesions influence processing in multiple systems and produces impairment across several cognitive domains (Gratton et al., 2012). Therefore, the functional backbone concept has important implications for the understanding the relationship between disruption of whole-brain network architecture and functional deficits in neurological or psychiatric disorders. Future work will need to further examine the importance of functional backbone of connector networks in clinical research with other experimental paradigms.

The limitations in this study must be mentioned as follows. First, the sample size is small and future studies should use larger sample to confirm the results. Second, dynamic networks are a very active and fast field. Recently, some novel approaches on temporal dynamic of brain networks have been put forward (Sizemore and Bassett, 2018; Vidaurre et al., 2018). Future studies should use these new approaches to further investigate how connector networks dynamically evolute and verify the importance of the functional backbone. Third, since the hard- and soft-thresholding have different ideas about how to deal with weak and negative connections, future studies might use the simulating methods to completely compare the soft-thresholding approach with the traditional approach so as to further clarify their appropriateness.

In the summary, the present findings lay stress on the necessity of conserving negative connections between multiple modules in chronnectome research, and provide evidence for the existence of functional backbone during dynamic evolution of the connector networks. These results deepen the understanding of the relationship between structural and functional connectivity during dynamic processing of affective information.

Statements

Data availability statement

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics statement

All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Ethical Committee of School of Psychology at East China Normal University.

Author contributions

QL, GR, and WZ designed and performed the study and analyzed the data. WZ wrote the manuscript. FT, HL, GR, and JC reviewed and modified the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (31470997 and 81171289).

Conflict of interest

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.01069/full#supplementary-material

References

  • 1

    BaiL.MaH.HuangY. X. (2005). The development of native Chinese affective picture system—a pretest in 46 college students.Chin. Ment. Heal. J.19719722.

  • 2

    BarretL. F.SatputeA. B. (2013). Large-scale brain networks in affective and social neuroscience: towards an integrative functional architecture of the brain.Curr. Opin. Neurobiol.23361372. 10.1016/j.conb.2012.12.012

  • 3

    BastianA. J. (2011). Moving, sensing and learning with cerebellar damage.Curr. Opin. Neurobiol.21596601. 10.1016/j.conb.2011.06.007

  • 4

    BaumannO.MattingleyJ. B. (2012). Functional topography of primary emotion processing in the human cerebellum.Neuroimage.61805811. 10.1016/j.neuroimage.2012.03.044

  • 5

    BetzelR. F.ByrgeL.HeY.GoniJ.ZuoX.SpornsO. (2014). Changes in structural and functional connectivity among resting-state networks across the human lifespan.Neuroimage102345357. 10.1016/j.neuroimage.2014.07.067

  • 6

    BetzelR. F.SatterthwaiteT. D.GoldJ. I.BassettD. S. (2017). Positive affect, surprise, and fatigue are correlates of network flexibility.Sci. Rep.7:520.

  • 7

    BrandesU. (2001). A faster algorithm for betweenness centrality.J. Math. Sociol.25163177. 10.1080/0022250x.2001.9990249

  • 8

    BucknerR. L.SepulcreJ.TalukdaT.KrienenF. M.LiuH.HeddenT.et al (2009). Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer’s disease.J. Neurosci.2918601873. 10.1523/jneurosci.5062-08.2009

  • 9

    CalhounV. D.MillerR.PearlsonG.AdaliT. (2014). The chronnectome: time-varying connectivity networks as the next frontier in f MRI data discovery.Neuron84262274. 10.1016/j.neuron.2014.10.015

  • 10

    ChangC.GloverG. H. (2009). Effects of model-based physiological noise correction on default mode network anti-correlations and correlations.Neuroimage4714481459. 10.1016/j.neuroimage.2009.05.012

  • 11

    De pisapiaN.TurattoM.LinP.JovicichJ.CaramazzaA. (2012). Unconscious priming instructions modulate activity in default and executive networks of the human brain.Cereb. Cortex22639649. 10.1093/cercor/bhr146

  • 12

    DrabantE. M.McRaeK.ManuckS. B.HaririA. R.GrossJ. J. (2009). Individual differences in typical reappraisal use predict amygdala and prefrontal responses.Biol. Psychiatry.65367373. 10.1016/j.biopsych.2008.09.007

  • 13

    FornitoA.HarrisonB. J.ZaleskyA.SimonsJ. S. (2012). Competitive and cooperative dynamics of large-scale brain functional networks supporting recollection.Proc. Natl. Acad. Sci. U.S.A.1091278812793. 10.1073/pnas.1204185109

  • 14

    FornitoA.ZaleskyA.BreakspearM. (2013). Graph analysis of the human connectome: promise, progress, andpitfalls.Neuroimage.80426444. 10.1016/j.neuroimage.2013.04.087

  • 15

    FoxM. D.SnyderA. Z.VincentJ. L.CorbettaM.Van EssenD. C.RaichleM. E. (2005). The human brain is intrinsically organized into dynamics, anticorrelated functional networks.Proc. Natl. Acad. Sci. U.S.A.10296739678. 10.1073/pnas.0504136102

  • 16

    GordonE. M.LaumannT. O.AdeyemoB.HuckinsJ. F.KelleyW. M.PetersenS. E. (2014). Generation and evaluation of a cortical area parcellation from resting-state correlations.Cereb. Cortex.26288303. 10.1093/cercor/bhu239

  • 17

    GrattonC.NomuraE. M.PérezF.D’EspositM. (2012). Focal brain lesions to critical locations cause widespread disruption of the modular organization of the brain.J. Cogn. Neurosci.2412751285. 10.1162/jocn_a_00222

  • 18

    GuS.SatterthwaiteT. D.MedagliaJ. D.YangM.GurR. E.GurR. C.et al (2015). Emergence of system roles in normative neurodevelopment.Proc. Natl. Acad. Sci. U.S.A.1121368113686. 10.1073/pnas.1502829112

  • 19

    HagmannP.CammounL.GigandetX.MeuliR.HoneyC. J.WedeenV. J.et al (2008). Mapping the structural core of human cerebral cortex.PLoS. Biol.6:e159. 10.1371/journal.pbio.0060159

  • 20

    HermundstadA. M.BassettD. S.BrownK. S.AminoffE. M.ClewettD.FreemanS.et al (2013). Structural foundations of resting-state and task-based functional connectivity in the human brain.Proc. Natl. Acad. Sci. U.S.A.11061696174. 10.1073/pnas.1219562110

  • 21

    LiangX.ZouQ.HeY.YangY. (2013). Coupling of functional connectivity and regional cerebral blood flow reveals a physiological basis for network hubs of the human brain.Proc. Natl. Acad. Sci. U.S.A.11019291934. 10.1073/pnas.1214900110

  • 22

    LindquistK. A.BarrettL. F. (2012). A functional architecture of the human brain: emerging insights from the science of emotion.Trends. Cogn. Sci.26533540. 10.1016/j.tics.2012.09.005

  • 23

    LindquistM. A.XuY.NebelM. B.CaffoB. S. (2014). Evaluating dynamic bivariate correlations in resting-state f MRI: a comparison study and a new approach.Neuroimage101531546. 10.1016/j.neuroimage.2014.06.052

  • 24

    MantzarisA. V.BassettD. S.WymbsN. F.EstradaE.PorterM. A.MuchaP. J.et al (2013). Dynamic network centrality summarizes learning in the human brain.J. Complex Netw.18392. 10.1093/comnet/cnt001

  • 25

    MarekS.HwangK.ForanW.HallquistM. N.LunaB. (2015). The contribution of network orgnization and integration to the development of cognitive control.PloS Biol.13:e1002326. 10.1371/journal.pbio.1002328

  • 26

    MesulamM. M. (1990). Large-scale neurocognitive networks and distributed processing for attention, language, and memory.Ann. Neurolog.28597613. 10.1002/ana.410280502

  • 27

    MoultonE. A.ElmanI.PendseG.SchmahmannJ.BecerraL.BorsookD. (2011). Aversion-related circuitry in the cerebellum: responses to noxious heat and unpleasant images.J. Neurosci.3137953804. 10.1523/JNEUROSCI.6709-10.2011

  • 28

    MumfordJ.HorvathS.OldhamM.LangfelderP.GeschwindD.PoldrackR. (2010). Detecting network modules in fMRI time series: a weighted network analysis approach.Neuroimage.5214561476. 10.1016/j.neuroimage.2010.05.047

  • 29

    OchsnerK.SilversJ.BuhleJ. (2012). Functional imaging studies of emotion regulation: a synthetic review and evolving model of the cognitive control of emotion.Ann. N. Y. Acad. Sci.1251E1E24. 10.1111/j.1749-6632.2012.06751.x

  • 30

    OchsnerK. N.BungeS. A.GrossJ. J.GabrieliJ. D. (2002). Rethinking feelings: an fMRI study of the cognitive regulation of emotion.J. Cogn. Neurosci.1412151229. 10.1162/089892902760807212

  • 31

    OosterwijkS.LindquistK. A.AndsonE.DautoffR.MoriguchiY.BarrettL. F. (2012). States of mind: emotions, body feelings and thoughts share distributed neural networks.Neuroimage6221102128. 10.1016/j.neuroimage.2012.05.079

  • 32

    PalaniyappanL.SimmoniteM.WhiteT. P.LiddleE. B.LiddleP. F. (2013). Neural primacy of the salience processing system in schizophrenia.Neuron79814828. 10.1016/j.neuron.2013.06.027

  • 33

    PowerJ. D.CohenA. L.NelsonS. M.WigG. S.BarnesK. A.ChurchJ. A.et al (2011). Functional network organization of the human brain.Neuron72665678.

  • 34

    PowerJ. D.SchlaggarB. L.Lessov-SchlaggarC. N.PetersenS. E. (2013). Evidence for hubs in human functional brain networks.Neuron79798813. 10.1016/j.neuron.2013.07.035

  • 35

    QuartoT.BlasiG.MaddalenaC.ViscantiG.LancianoT.SoletiE.et al (2016). Association between ability emotional intelligence and left insula during social judgment of facial emotions.PLoS One11:e0148621. 10.1371/journal.pone.0148621

  • 36

    RajA.LoCastroE.KuceyeskiA.TosunD.RelkinN.WeinerM.et al (2015). Network diffusion model of progression predicts longitudinal patterns of atrophy and metabolism in Alzheimer’sdisease.Cell Rep.10359369. 10.1016/j.celrep.2014.12.034

  • 37

    RubinovM.SpornsO. (2010). Complex network measures of brain connectivity, Uses and interpretations.Neuroimage5210591069. 10.1016/j.neuroimage.2009.10.003

  • 38

    RubinovM.SpornsO. (2011). Weight-conserving characterization of complex functional brain networks.Neuroiamge5620682079. 10.1016/j.neuroimage.2011.03.069

  • 39

    SchwarzA. J.McGonigleJ. (2011). Negative edges and soft thresholding in complex network analysis of resting state functional connectivity data.Neuroimage5511321146. 10.1016/j.neuroimage.2010.12.047

  • 40

    SeeleyW. W.MenonV.SchatzbergA. F.KellerJ.GloverG. H.KennaH.et al (2007). Dissociable intrinsic connectivity networks for salience processing and executive control.J. Neurosci.2723492356. 10.1523/jneurosci.5587-06.2007

  • 41

    SepulcreJ.LiuH.TalukdarT.MartincorenaI.YeoB. T. T.BucknerR. L. (2010). The organization of local and distant functional connectivity in the human brain.PLoS. Comput. Biol.6:e1000808. 10.1371/journal.pcbi.1000808

  • 42

    ShenX.TokoqluF.PapademetrisX.ConstableR. T. (2013). Groupwise whole-brain parcellation from resting-state fMRI data for network node identification.Neuroimage82403415. 10.1016/j.neuroimage.2013.05.081

  • 43

    ShineJ. M.KoyejoO.PoldrackR. A. (2016). Temporal metastates are associated withdifferential patterns of time-resolved connectivity, network topology, and attention.Proc. Natl. Acad. Sci. U. S. A.11398889891. 10.1073/pnas.1604898113

  • 44

    SizemoreA. E.BassettD. S. (2018). Dynamic graph metrics: tutorial, toolbox, and tale.Neuroimage180417427. 10.1016/j.neuroimage.2017.06.081

  • 45

    SpornsO. (2014). Contributions and challenges for network models in cognitive neuroscience.Nat. Neurosci.17652660. 10.1038/nn.3690

  • 46

    SripadaC.AngstadtM.KesslerD.PhanK. L.LiberzonI.EvansG. W.et al (2014). Volitional regulation of emotions produces distributed alterations in connectivity between visual, attention control, and default networks.Neuroimage89110121. 10.1016/j.neuroimage.2013.11.006

  • 47

    TaylorK. S.SeminowiczD. A.DavisK. D. (2009). Two systems of resting state connectivity between the insula and cingulate cortex.Hum. Brain. Mapp.3027312745. 10.1002/hbm.20705

  • 48

    UrryH. L.Van-ReekumC. M.JohnstoneT.KalinN. H.ThurowM. E.SchaeferH. S.et al (2006). Amygdala and ventromedial prefrontal cortex are inversely coupled during regulation of negative affect and predict the diurnal pattern of cortisol Secretion among older adults.J. Neurosci.2644154425. 10.1523/jneurosci.3215-05.2006

  • 49

    Van den HeuvelM. P.SpornsO. (2011). Rich-club organization of the human connectome.J. Neurosci.111577515786. 10.1523/JNEUROSCI.3539-11.2011

  • 50

    Van den HeuvelM. P.SpornsO. (2013). Network hubs in the human brain.Trends Cogn. Sci.17683696. 10.1016/j.tics.2013.09.012

  • 51

    VaroquauxG.CraddockR. C. (2013). Learning and comparing functional connectomes across subjects.Neuroimage80405415. 10.1016/j.neuroimage.2013.04.007

  • 52

    VidaurreD.AbeysuriyaR.BeckerR.QuinnA. J.AlmagroF. A.SmithS. M.et al (2018). Discovering dynamic brain networks from big data in rest and task.Neuroimage180646656. 10.1016/j.neuroimage.2017.06.077

  • 53

    VincentJ. L.KahnI.SnyderA. Z.RaichleM. E.BucknerR. L. (2008). Evidence for a frontoparietal control system revealed by intrinsic functional connectivity.J. Neurophysiol.10033283342. 10.1152/jn.90355.2008

  • 54

    VosselS.WeidnerR.DriverJ.FristonK. J.FinkG. R. (2012). Deconstructing the architecture of dorsal and ventral attention systems with dynamic causal modeling.J. Neurosci.321063710648. 10.1523/JNEUROSCI.0414-12.2012

  • 55

    WagerT. D.DavidsonM. L.HughesB. L.LindquistM. A.OchsnerK. N. (2008). Prefrontal-subcortical pathways mediating successful emotion regulation.Neuron5910371050. 10.1016/j.neuron.2008.09.006

  • 56

    WarrenD. E.PowerJ. D.BrussJ.DenburgN. L.WaldronE. J.SunH.et al (2014). Network measures of predict neuropsychological outcome after brain injury.Proc. Natl. Acad. Sci. U.S.A.1111424714252. 10.1073/pnas.1322173111

  • 57

    WolfU.RapoportM. J.SchweizerT. A. (2009). Evaluating the affective component of the cerebellar cognitive affective syndrome.J. Neuropsychiatr. Clin. Neurosci.21245253. 10.1176/appi.neuropsych.21.3.245

  • 58

    WoolgarA.AfsharS.WilliamsM.WilliamsM. A.RichA. N. (2015). Flexible coding of task rules in frontoparietal cortex: an adaptive system for flexible cognitive control.J. Cogn. Neurosci.2718951911. 10.1162/jocn_a_00827

  • 59

    YanC.GongG.WangJ.WangD.LiuD.ZhuC.et al (2011). Sex- and brain size-related small-world structural cortical networks in young adults: a DTI tractography study.Cereb. Cortex21449458. 10.1093/cercor/bhq111

  • 60

    YeoB. T. T.KrienenF. M.SepulcreJ.SabuncuM. R.LashkariD.HollinsheadM.et al (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity.J. Neurophysiol.10611251165. 10.1152/jn.00338.2011

  • 61

    ZhangW.LiH.PanX. (2015). Positive and negative affective processing exhibit dissociable functional hubs during the viewing of affective pictures.Hum. Brain. Mapp.36415426. 10.1002/hbm.22636

Summary

Keywords

affective pictures, cognitive control, connector hubs, graph theory, negative connections

Citation

Zhang W, Chen J, Ren G, Tang F, Liu Q and Li H (2019) Negatively Linking Connector Networks in Cognitive Control of Affective Pictures. Front. Neurosci. 13:1069. doi: 10.3389/fnins.2019.01069

Received

13 May 2019

Accepted

24 September 2019

Published

25 October 2019

Volume

13 - 2019

Edited by

Feng Liu, Tianjin Medical University General Hospital, China

Reviewed by

Guangheng Dong, Hangzhou Normal University, China; Xiangyi Zhang, Hunan Normal University, China

Updates

Copyright

*Correspondence: Wenhai Zhang, Hong Li,

Co-first authors

This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics