Abstract
Cognitively demanding goal-directed tasks in the human brain are thought to involve the dynamic interplay of several large-scale neural networks, including the default-mode network (DMN), salience network (SN), and central-executive network (CEN). Resting-state functional magnetic resonance imaging (rsfMRI) studies have consistently shown that the CEN and SN negatively regulate activity in the DMN, and this switching is argued to be controlled by the right anterior insula (rAI) of the SN. However, what remains to be investigated is the pattern of directed network interactions during difficult perceptual decision-making tasks. We recorded fMRI data while participants categorized the left–right motion of moving dots. We defined regions of interest, extracted fMRI time series, and performed directed connectivity analysis using Granger causality techniques. Our analyses demonstrated that the slow oscillation (0.07–0.19 Hz) mediated the interactions within and between the DMN, SN, and CEN nodes both for easier and harder decision-making tasks. We found that the rAI, a key node of the SN, played a causal control over the DMN and CEN for easier decision-making tasks. The combined effort of the rAI and dorsal anterior cingulate cortex of the SN had the causal control over the DMN and CEN for a harder task. These findings provide important insights into how a sensory signal organizes among the DMN, SN, and CEN during sensory information-guided, goal-directed tasks.
Introduction
Previous research has described a set of large-scale, intrinsically organized brain networks underlying a broad range of functions, from basic sensory and motor capacities to cognition and higher level functions such as selecting, switching, and attending to salience events in the surrounding (Becerra et al., 2011; Deco et al., 2011; Power et al., 2011; Seeley et al., 2007). Studies have identified the existence of at least three canonical networks: (1) the default-mode network (DMN) that includes the posterior cingulate cortex (PCC) and ventromedial prefrontal cortex (VMPFC); (2) the salience network (SN) that has the right anterior insula (rAI) and dorsal anterior cingulate cortex (DACC); and (3) the central-executive network (CEN), whose key regions include the posterior parietal cortex (PPC) and dorsolateral prefrontal cortex (DLPFC) (Chen et al., 2013; Uddin, 2015).
Resting-state functional magnetic resonance imaging (rsfMRI) studies have consistently shown that the CEN and/or SN negatively regulate activity in the DMN (Chen et al., 2013; Fox et al., 2005; Fransson, 2005). However, what remains to be understood is how the activity of these networks, identified in the rsfMRI, relates to their function in cognitively demanding goal-directed tasks such as perceptual decision-making (PDM).
There are two competing lines of evidence in trying to explain how sensory signal enters SN (Ham et al., 2013) and SN might mediate interactions between DMN and CEN. One line of evidence suggests that the DACC monitors performance and its role is important for behavioral adaptation (Ridderinkhof et al., 2004). Activity in the DACC indicates an enhanced cognitive control, and its interactions with the lateral prefrontal structures implement subsequent behavioral changes (Egner, 2009; Ridderinkhof et al., 2004). Another line of evidence suggests that the rAI is a cortical outflow hub of the SN and it coordinates a change in activity across multiple brain networks, including the DMN and CEN (Bonnelle et al., 2012; Menon and Uddin, 2010; Sridharan et al., 2008). In this study, we seek to examine the interactions among the DMN, SN, and CEN, with particular interest in the roles of the rAI and/or DACC in regulating these key networks and their modulations with difficulty of PDM.
We used the traditional left–right motion categorization of moving dots similar to the previous study (Heekeren et al., 2006), conducted an fMRI experiment, and applied the spectral Granger causality (GC) analysis method (Dhamala et al., 2008a, 2008b) to look at the interactions among the DMN, SN, and CEN. We hypothesized that the DMN and CEN would be under causal control from the SN, and if so, a node of the SN would regulate this control. We further hypothesized whether a node that achieved the causal control for easier task would remain the same or change for harder decision-making tasks.
Materials and Methods
Participants
Thirty-two human volunteers (16 males, 16 females; mean age: 27.6 years; standard deviation: 4.7 years) participated in the experiment. All participants reported neurologically healthy and had normal or corrected to normal vision. Written informed consent was collected from each participant before data collection. The experimental protocol was approved from the Institutional Review Board of Center for Advanced Brain Imaging (CABI), Georgia State University and Georgia Institute of Technology. Two participants were excluded from the final analyses because of missing behavior performance records.
Task paradigm
We created random green and red dots. They were presented on a black background and were drawn in a circular aperture. Dots were redrawn after some delay at either a random location or a neighboring spatial location to induce apparent motion to create the resultant motion effect of random dots. The coherence level was determined by the fraction of dots displaced in apparent motion. We presented dots with 12%, 20%, and 50% coherence level (coherence dots were red and green). We also presented 12%, 20%, and 50% coherence dots that had only green color, and they are named color coherence henceforth. The stimuli were presented for 24-sec block, followed by a question mark (‘?’), after which time participants responded with a keyboard button press, and then a rest period of 10 sec before starting the next stimulus. We had rest periods of 30 sec at the beginning and at the end of each run. Participants performed two functional runs, each of 810 sec long. The total number of trials was 36 (six trials for each condition). In both sessions, participants categorized random dots as either moving left or right and indicated their decisions by button responses. The presentation software (
Data acquisition
Imaging was done on a 3-T Siemens MRI scanner, available at CABI of Georgia State University and Georgia Institute of Technology. High-resolution T1-weighted anatomical images were acquired first and then followed by two whole-brain functional runs with 405 scans per run with the following parameters: echo-planar imaging, gradient recalled echo sequence; repetition time (TR) = 2000 msec; echo time (TE) = 30 msec; flip angle = 90°; 68 × 68 matrix, voxel dimensions = 3 × 3 × 3 mm3, 37 axial slices each of 3 mm thickness acquired in an interleaved manner with orientation parallel to the anterior–posterior commissural line for the measurement of the T2*-weighted blood oxygenation level-dependent effect. High-resolution anatomical T1-weighted images were acquired for anatomical references using a magnetization-prepared rapid acquisition gradient echo sequence with TR = 2250 msec; TE = 4.18 msec; inversion time = 900 msec; flip angle = 9°; and voxel dimensions = 1 × 1 × 1 mm3.
Data preprocessing
We preprocessed fMRI data using SPM8 (Wellcome Trust Centre for Neuroimaging, London, United Kingdom;
Selection of regions of interest and extraction of time series
We selected regions of interest (ROIs) using the MarsBaR software package (
Computation of spectral power and GC
Power spectra were computed using parametric and nonparametric approaches (Chand and Dhamala, 2014; Dhamala et al., 2008a, 2008b). GC spectral analyses are primarily used to examine the strengths, directions, and frequencies of interactions between dynamic processes. The GC spectrum from the second time series X2 to the first time series X1 (i.e., 2→1) at a frequency (f) is defined as (Dhamala et al., 2008a, 2008b)
where S is spectral power, H is transfer function, and Σ is noise covariance.
For N nodes, the frequency-specific causal outflow (F) at a node i can be defined as follows:
We used time series extracted from each ROI to compute GC among the nodes of DMN, SN, and CEN. To compute group level GC, time series from all subjects were considered. To obtain the net causal outflow, we first computed GC among all ROIs for each subject, used Equation (2) for each ROI for each subject, and then calculated the mean and standard error of the mean (SEM) over subjects.
The threshold value of GC, for statistical significance, was computed from surrogate data methods by using data permutation calculating GC values and a gamma-function to a distribution of maximum GC values from each permutation (Adhikari et al., 2014; Blair and Karniski, 1993). This threshold was designed to reject a null hypothesis of no interdependence at a significance level of p < 10−3.
Results
Behavioral results
The performance (%), which is a ratio of the number of correctly responded trials to the total number of presented trials multiplied by hundred, was calculated for coherence levels of moving dots. Based on t-tests, the performance (%) was significantly higher for 50% coherence level (represented as 50% henceforth) compared to the 12% coherence level (represented as 12%). The 20% coherence level did not show statistically significant difference with 50% and 12% coherence levels and hence discarded from subsequent analyses. Only stimuli with 12% color coherence (represented as 12% color coherence) showed statistically significant difference (paired t-test) with 12% coherence, and therefore, other color coherence levels (20% and 50% color coherences) were not considered for the subsequent analyses (Fig. 1).

Behavior responses.
fMRI results
Individual DMN, SN, and CEN for 50% coherence
Power spectra computed for the nodes of the DMN (the PCC and VMPFC), SN (the rAI and DACC), and CEN (the PPC and DLPFC), respectively, showed peak activities in slow frequency band (0.07–0.19 Hz), often referred as slow-3 band (Bajaj et al., 2013; Zuo et al., 2010), when the participants perceived 50% coherence dots (Fig. 2 and Supplementary Fig. S1; Supplementary Data are available online at

Selection of the ROIs:
Figure 3 presents GC spectra as a function of frequency, where horizontal lines (Fig. 3A–C) represent a statistically significant threshold value. In the DMN, only the causal flow from the PCC to VMPFC is higher than the threshold value (Fig. 3A). The causal flow is bidirectional between the rAI and DACC in the SN; however, a dominant flow is from the rAI to DACC (Fig. 3B). In the CEN, the causal flow from the PPC to DLPFC is significant but not the other way around. In addition, we also calculated a time-reversed GC (Haufe et al., 2013; Vinck et al., 2015) to cross-validate our GC results. Time-reversed GC calculations demonstrated an altered dominant direction within the DMN, SN, and CEN as expected (Supplementary Fig. S4).

GC spectra within the
Individual DMN, SN, and CEN for 50% versus 12% coherence
Power spectra computed at the nodes of the DMN (the PCC and VMPFC), SN (the rAI and DACC), and CEN (the PPC and DLPFC), respectively, also showed peak activity in slow-3 band (0.07–0.19 Hz) when the participants perceived 12% coherence dots (Supplementary Fig. S2).
For 12% coherence, dominant causal flow is still from the PCC to VMPFC in the DMN; however, the value is lower than that for 50% coherence (Fig. 4). In the SN, dominant causal flow is still from the rAI to DACC, but with suppressed GC compared to that of 50% coherence. The dominant causal flow is still from the PPC to DLPFC in the CEN.

GC spectra within the
Individual DMN, SN, and CEN for 12% coherence versus 12% color coherence
Power spectra computed at the nodes of the DMN (the PCC and VMPFC), SN (the rAI and DACC), and CEN (the PPC and DLPFC), respectively, also showed peak activity in slow-3 band (0.07–0.19 Hz) when the participants perceived dots with 12% color coherence (Supplementary Fig. S3).
For 12% color coherence, dominant causal flow is still from the PCC to VMPFC in the DMN; however, its strength is enhanced compared to that for 12% coherence (Fig. 5). In the SN, dominant causal flow is still from the rAI to DACC, but the connectivity strength is enhanced compared to that for 12% coherence. The dominant causal flow is from the PPC to DLPFC in the CEN in 12% color coherence.

GC spectra within the
Causal interaction among the DMN, SN, and CEN
In 50% coherence dots, we computed net causal outflow (out-in degree) at the key nodes of the SN, CEN, and DMN. The comparison of the net causal outflows revealed that the rAI had a significantly higher causal control over the DACC and the nodes of the DMN and CEN (Fig. 6A).

Net causal outflow of the key nodes of the SN, CEN, and DMN.
We further calculated net causal outflows for 12% coherence dots, which is a harder task compared to 50% coherence as reflected from behavioral measures (Fig. 1A). In 12% coherence dots (Fig. 6B), comparison of the net causal outflows among the key nodes of the SN, CEN, and DMN uncovered that, in addition to the rAI, the DACC had a significantly higher causal control over the CEN and DMN regions (Supplementary Fig. S5 for 20% coherence dots). The net causal outflow comparison between the rAI and DACC is not significantly different. We further calculated net causal outflow for 12% color coherence, which is an easier task compared to 12% coherence as reflected from behavioral measures (Fig. 1B). For 12% color coherence dots, comparison of the net causal outflows among the key nodes of the SN, CEN, and DMN revealed that the rAI had a significantly higher net causal outflow than the CEN and DMN regions (Fig. 6C). The net causal outflows of the rAI and DACC are not significantly different, but showed an enhanced causal outflow of the rAI and a suppressed causal outflow of the DACC compared to the 12% coherence (Fig. 6B, C).
Discussion
We demonstrated that an individual DMN, SN, and CEN had dominantly feed-forward network activity mediated by slow-frequency band (0.07–0.19 Hz) often referred as slow-3 band (Bajaj et al., 2013; Zuo et al., 2010). We found two important results that showed organization of sensory information among the DMN, SN, and CEN. First, the rAI of the SN played causal control over the DMN and CEN for easier decision, which is consistent with previous investigation for different sensory information-driven, goal-directed tasks (Sridharan et al., 2008). Second, the rAI played a relatively lower causal control and the rAI and DACC exerted a causal control over the DMN and CEN (Fig. 6B). Moreover, when color dots were present (12% color coherence), task became behaviorally easier and an increased role of rAI (compared to 12% coherence) was observed (Fig. 6C).
Recent studies have consistently demonstrated that the rAI of the SN drives the DMN and CEN for different tasks (Goulden et al., 2014; Sridharan et al., 2008). Our results for easier decision-making tasks also replicated those existing studies. In addition, our result for harder decision-making tasks demonstrated further that both the rAI and DACC of the SN causally control the DMN and CEN. Anatomically, the AI and anterior cingulate cortex (ACC) are a part of a network (Critchley et al., 2004) and also share a unique feature at the neuronal level. The von Economo neurons (VENs)—a special type of neuron exclusively localized to the AI and ACC—relay information processed within the AI and ACC to other parts of the brain (Allman et al., 2005; Watson et al., 2006). The control signals generated by the rAI (for easier decisions) and by both the DACC and rAI (for harder decisions) might be supported by the neuronal basis of control signals of the VENs. Previous studies of attention and cognitive control have reported coactivation of the rAI and DACC (Crottaz-Herbette and Menon, 2006; Dosenbach et al., 2006; Ham et al., 2013), but their differential functional roles are still in debate. The rAI is functionally connected to the networks responsible for adaptive behavior, including the SN (Seeley et al., 2007), as well as other parts of the frontoparietal control network (Vincent et al., 2008). Diffusion tensor imaging demonstrated that this cortical area has direct white matter connections to other key regions, including the DACC (van den Heuvel et al., 2009), the inferior parietal lobe (Uddin et al., 2010), and the temporoparietal junction (Kucyi et al., 2012), making the insula (rAI) well placed to perform its putative role of evaluating (Uddin et al., 2010), reorienting attention (Ullsperger et al., 2010), and switching between cognitive resources in response to salient events (Uddin and Menon, 2009), similar to the line that we found for easier decisions. In contrast, the activity in the DACC is known to signal the need for enhanced cognitive control (Ridderinkhof et al., 2004) and to implement behavioral changes (Egner, 2009; Ridderinkhof et al., 2004). Our results for harder decisions might implicate the role of the DACC for an enhanced cognitive control (Egner, 2009; Ridderinkhof et al., 2004) especially when sensory information is limited (harder decisions). Our findings for both easier and harder decisions therefore provide strong evidence for the neural basis of rAI–DACC network to initiate control signals in sensory information-driven, goal-directed tasks (Allman et al., 2005; Watson et al., 2006).
In conclusion, we evaluated the patterns of interactions among DMN, SN, and CEN in perceptual decisions of moving dots. We found that the slow-3 oscillation (0.07–0.19 Hz) mediated interaction within and between these networks, the rAI of the SN played causal control over the DMN and CEN for easier decisions, and finally, both the rAI and DACC of the SN sent causal control in harder decisions. These findings advance our current understanding of how sensory signal organizes among the DMN, SN, and CEN for sensory-information, goal-directed tasks, including easier and harder perceptual decisions.
Footnotes
Acknowledgments
The authors thank Crystal Smith for writing computer programs to control and present the stimuli, Drs. Bidhan Lamichhane and Bhim Adhikari for the data collection, and Timothy Jordan for proofreading the manuscript. A National Science Foundation CAREER award (BCS 0955037) was given to the author MD financially supporting this work.
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
