Cognitive Control in the Controllable Connectome
Cognition is supported by neurophysiological processes that occur both in local anatomical neighborhoods and in distributed large-scale circuits. Recent evidence from network control theory suggests that white matter pathways linking large-scale brain regions provide a critical substrate constraining the ability of single areas to affect control on those processes. Yet, no direct evidence exists for a relationship between brain network controllability and cognitive control performance. Here, we address this gap by constructing structural brain networks from diffusion tensor imaging data acquired in 125 healthy adult individuals. We define a simplified model of brain dynamics and simulate network control to quantify modal and boundary controllability, which together describe complementary features of a region’s theoretically predicted preference to drive the brain into different cognitive states. We observe that individual differences in these control features derived from structural connectivity are significantly correlated with individual differences in cognitive control performance, as measured by a continuous performance attention test, a color/shape switching task, the Stroop inhibition task, and a spatial n-back working memory task. Indeed, control hubs like anterior cingulate are distinguished from default mode and frontal association areas in terms of the relationship between their control properties and individual differences in cognitive function. These results provide the first empirical evidence that network control forms a fundamental mechanism of cognitive control.
Cognitive control refers to the ability to adaptively vary information processing and subsequent behavior to achieve specific goals. This ability enables humans to link information to solve problems, inhibit inappropriate behavioral responses, switch among tasks, and actively update, select, and maintain behaviorally relevant information. Failures in cognitive control are ubiquitously observed in neuropsychiatric conditions and predict long-term outcomes ranging from school performance [1, 2] to career achievement .
Concerted efforts in cognitive neuroscience have demonstrated that these functions have specific anatomical bases in the brain, including a distributed set of areas in frontal, medial, and parietal cortices [3, 4, 5]. Yet, how these regions work in concert with one another remains far from understood . One emerging perspective is that cognitive control is driven by dynamic interactions between large-scale neural circuits or networks [7, 8, 9]. These interactions can take the form of coordination or competition , and can change over time to enable behavioral outputs such as button presses and spoken words. Indeed, the network-level processes of cognitive control may facilitate the flexible management of large-scale brain activity in support of successful performance in a variety of effortful tasks [3, 10, 11].
While cognitive control is thought of as a network-level process, putative mechanisms of this process have not been identified. Recent theoretical work posits that network control theory - a nascent field of engineering - offers exactly such a fundamental mechanism . Network control theory is the study of how to design control strategies for networked systems , in which a set of nodes are connected by edges and in which the activity of nodes can be modeled using simplified linear dynamics. These tools predict that different nodes serve dissociable control roles defined by their location on the network . In the context of the brain, such a mechanism would suggest that brain regions are predisposed to drive or modulate neurophysiological dynamics in a manner consistent with their specific topological signatures of white matter connectivity [8, 14].
Here we test whether network control is a putative mechanism of cognitive control by asking whether the theoretically predicted control features of brain regions are related to cognitive performance on so-called cognitive control tasks. We restrict our attention to two distinct network control features known as modal controllability and boundary controllability. Intuitively, modal controllability describes the ability of a node (here, a brain region) to drive a network into difficult-to-reach states on an energy landscape, thus providing a means for switching between disparate tasks. Based on prior work locating modal controllers in fronto-parietal cognitive control systems , we hypothesize that individual differences in the modal controllability of these areas will be associated with individual differences in performance on demanding cognitive control tasks such as those with high working memory or inhibition demands, or switching costs. To complement these intuitions, we also study boundary controllability, which describes the ability of a node to steer the system into states where modules are either coupled or decoupled, thus providing a means for coordination and competition. Based on prior work locating boundary controllers in attention systems , we hypothesize that individual differences in the boundary controllability of these areas will be associated with individual differences in the performance of attention-demanding tasks.
To test these hypotheses, we construct structural brain networks from diffusion tensor imaging (DTI) data acquired in each of 125 healthy adult subjects (Fig. 1A). Each network contains 83 brain regions defined by the Lausanne anatomical parcellation (Fig. 1B), and each pair of regions is connected by an edge weighted by the number of streamlines linking those regions (Fig. 1C). We define a simplified model of brain dynamics and simulate network control to quantify modal and boundary controllability (Fig. 1D). We reveal intersubject variability in node controllability that is associated with different cognitive control functions , either being advantageous or disadvantageous for cognitive control performance (Fig. 1E). Performance is measured by the number of true positives on the continuous performance test , color/shape median switch cost (reaction time) , median Stroop effect (reaction time) , and a visuo-spatial n-back . The continuous performance test measures sustained attention or vigilance. The color/shape task median switch cost measures an individual’s ability to switch between task sets efficiently. The Stroop effect measures an individual’s ability to inhibit prepotent responses to semantic information. The visuo-spatial n-back measures the ability to hold and update a limited amount of information in mind. By examining the relationship between interindividual variability in controllability and distinct cognitive control functions, we establish a bridge between a neuroscientist’s notion of cognitive control and an engineer’s notion of network control, suggesting that the latter forms a mechanism for the former.
Individual Variability in Cognitive Measures
All cognitive measures were collected as part of a testing battery designed to assess cognitive control and decision making (see Methods for details; see Fig. 2). The four tests we study were selected based on (1) their statistical independence from one another, (2) their use in the measurement of distinct cognitive control processes, and (3) their sensitivity to disease diagnoses. We observed that these four tests contained little shared variance with one another in explaining individual differences in cognitive control performance (). These results indicate that the four measures represent unique information that can assist in distinguishing the cognitive control abilities of one person from those of another.
Cognitive Control Performance and Network Controllability
We next ask whether individual differences in cognitive control performance (as measured by these tasks) are related to individual differences in white matter architecture, the substrate of an engineer’s notion of structural network control. To address this question, we calculate the Spearman correlation coefficient between regional controllability values and performance. To determine statistical significance, we apply a spatial test to determine non-trivial structure in the anatomical distributions of these correlations (see Methods). Briefly, we determine the statistical similarity (Pearson correlation) between the inverse Euclidean distance between node centroid pairs and the mean absolute value of cognition-control correlations for the same node pair. A positive correlation demonstrates that the spatial proximity of nodes is related to the strength of their controllability relationship to cognition.
We report results for each of the four cognitive control performance measures, as well as for each of the two regional controllability statistics (modal and boundary). Across all eight analyses, we observed that each map demonstrated highly significant spatial organization in the relationships between individual differences in controllability and cognitive performance. For modal controllability, we observed the following relationships: CPT: , , ; Color/Shape: , , ; Stroop: , , ; N-Back: , , . For boundary controllability, we observed the following relationships: CPT: , , ; Color/Shape: , , ; Stroop: , , ; N-Back: , , . These results demonstrate that brain regions whose controllability values are associated with individual differences in cognitive control performance display non-trivial spatial clustering.
Importantly, these maps reveal differential relationships between controllability and cognitive control: increases in controllability in some regions are advantageous for cognitive control performance, whereas increases in others are disadvantageous. In the following sections, we examine the maps for each cognitive control performance measure separately, report significant regional predictors of performance, and interpret their anatomical locations.
Regional Controllers in the Continuous Performance Task
We begin by examining the anatomical distribution of correlations between controllability values estimated from structural brain networks and the number of true positives on the continuous performance test (Fig. 3). We interpret regions for which -to- transformed values of correlation surpass . We observe that modal controllability in the left inferior frontal gyrus, posterior cingulate, and right precuneus are most positively associated with performance, whereas controllability in the bilateral medial frontal gyri were most negatively associated with performance. Boundary controllability in the bilateral caudal anterior cingulate are advantageous for performance, whereas increasing strength in the right inferior parietal lobe, lateral temporal lobe and amygdala are relatively disadvantageous.
The anatomical loci of these control-cognition relationships are particularly interesting in light of the posited neuroanatomical basis of continuous performance task execution. For example, recent functional neuroimaging studies have identified activation in the cingulate cortex and precuneus during continuous performance tests . More recently, data has further supported an association between disrupted posterior cingulate function and arousal states, including states of internal or externally focused attention, and states modulating breadth of mental focus . The precuneus is typically suppressed during attention processing, and its deactivation is associated with stability of attention . Moreover, sustained attention is also associated with activation in anterior cingulate cortex, which is presumed to contribute to associated error monitoring . More generally, the inferior parietal lobe and medial frontal gyrus have been associated with attention shifts , and reduced activity in the lateral temporal lobe has been associated with reduced attention switching and inhibition performance . Furthermore, the ventral nucleus of the amygdala has been associated with attentional function in fear conditioning .
In light of this body of work in cognitive neuroscience, our results offer a fresh perspective on the potential structural mechanisms subserving these observations. Specifically, our findings suggest a tradeoff between advantageous and disadvantageous regional control roles: if regions that contribute to attention switching and emotional reactivity have an increased role in controlling the brain, maintaining sustained attention is more difficult. Conversely, stronger controllability situated in regions associated with sustaining attention leads to more stable attention in the continuous performance test. Sustained attention contributes to many domains of cognitive function, and is ubiquitously effected in neuropsychiatric disorders . Here, we find that increasing strength in the role of emotion processing and attention switching regions in driving brain wide dynamics may be one basis for these observed deficits.
Regional Controllers in the Color Shape Task
Next, we examine the anatomical distribution of correlations between controllability values estimated from structural brain networks and the median switch costs on the color shape task (Fig. 4). Again, we interpret regions for which -to- transformed values of correlation surpass . We observe that modal controllability in the right putamen and pallidum is positively associated with performance (reduced switch costs), whereas modal controllability in the left temporal pole, caudate, and insula is negatively associated with performance (enhanced switch costs). Boundary controllability in the bilateral nucleus accumbens is positively associated with performance, whereas boundary controllability in the right hippocampus, amygdala, and thalamus is negatively associated with performance.
The posited neuroanatomical basis of task switching paradigms traditionally place a special emphasis on the role of prefrontal regions [3, 4]. However, a broader view acknowledges that subcortical systems heavily interact with fronto-parietal areas supporting goal maintenance and switching [28, 29, 30]. Moreover, subcortical systems are critical for parsing and consolidating task representations during multitasking , and focal lesions in the basal ganglia result in errors when trial responses require the application of the correct task rule . Indeed, more generally, the basal ganglia are anatomically well-positioned to regulate information flow among multiple cortical regions, contributing highly regular computational functions in parallel-projecting segregated loops .
Our results complement these prior efforts by suggesting that strengthening in the role of the putamen and pallidum in driving the brain into difficult to reach states is associated with reduced switch costs, whereas strengthening in the role of the nucleus accumbens into integrated and segregated states is associated with reduced switch costs. In light of these findings, it is intuitively plausible that task-switching deficits in classically subcortical neurodegenerative diseases such as Parkinson’s disease may be attributable to disruption of the basal ganglia’s role in mediating trajectories that switch cognitive sets . Interestingly, in contrast to the supportive role of the putamen and pallidum, we find that high modal controllability in the caudate is disadvantageous for task switching, as is high boundary controllability in the hippocampus, amygdala, and thalamus. These findings indicate that as regions that are relatively specialized for long term memory storage and consolidation, reinforcement and reward learning, sensory processing and relays, and planned action in unpredictable environments become more dominant in the control hierarchy, fluid cognitive transitions may become more difficult. It is conceivable that relative increases in their role regulating intermodular communication competes with the role of the basal ganglia in regulating cortical dynamics during cognitive transitions. This may potentially explain observations of deficits during task switching in anxiety , which is associated with altered function and structure in the hippocampus, thalamus, and amygdala .
Regional Controllers in the Stroop Task
Next, we examine the anatomical distribution of correlations between controllability values estimated from structural brain networks and the median response time on the Stroop task (Fig. 5). We observe that modal controllability in the right inferior temporal and left lateral occipital lobe is positively associated with performance, whereas modal controllability in the left cuneus is negatively associated with performance. Boundary controllability in the bilateral precuneus and cuneus is positively associated with performance, whereas boundary controllability in the left hippocampus and right middle temporal gyrus is negatively associated with performance.
The anatomical loci of these control-cognition relationships are directly relevant to Stroop activations, as well as to regions known to process different facets of the Stroop stimuli. For example, the left lateral occipital cortex supports right hemifield processing, is selectively activated in processing letter strings relative to face and texture processing , and is involved in accessing orthographic whole-word representations . The right inferior temporal gyrus participates in the ventral stream pathway, integrating color and shape information during object recognition and recall . The current results suggest that the ease in suppressing the prepotent effects of lexical information in the Stroop task is partially attributable to the ability of these regions to propel the brain into difficult-to-reach states.
Prior work using fMRI to investigate Stroop task performance have demonstrated the recruitment of the precuneus in the presence of many or few distractors . The current results suggest that the dynamic role of the precuneus and its support of performance depends on its underlying structural profile: increased involvement in regulating intermodule processing, potentially facilitating information suppression early after visual stimulus representation, reduces the Stroop effect, whereas an increased role in driving the brain into hard to reach states exaggerates the Stroop effect.
Regional Controllers in the N-Back Working Memory Task
Finally, we examine the anatomical distribution of correlations between controllability values estimated from structural brain networks and true positives on a spatial n-back task (Fig. 6). We observe that modal controllability in the right cingulate isthmus and superior parietal lobe is positively associated with performance, whereas modal controllability in the left entorhinal and temporal pole is negatively associated with performance. Boundary controllability in the left precuneus and anterior cingulate cortex and right superior frontal gyrus is positively associated with performance, whereas boundary controllability in the left thalamus and medial orbitofrontal gyrus is negatively associated with performance.
Again, the anatomical loci of these control-cognition relationships are directly relevant to n-back memory performance and consistent with previous functional neuroimaging findings. For example, the superior parietal lobe contributes to the maintenance and manipulation of information in working memory , and contributes to spatial attention processing . The cingulate isthmus plays a more general role, being anatomically positioned to regulate activity in the default mode network and the frontoparietal network recruited for memory performance . Similar to the CPT, increasing boundary controllability in the left anterior cingulate cortex was advantageous for performance, potentially enabling intermodular communication of error-related signals. Finally, the superior frontal gyrus supports working memory and attention , and participates in the frontoparietal control system responsible for moment-to-moment response monitoring and set maintenance . Increases in this region’s modal control role may allow the brain to better maintain the difficult states necessary to facilitate monitoring among distractors.
Interestingly, in contrast to the supportive role of these areas, high control values in entorhinal cortex, temporal pole, and orbitofrontal gyrus is be detrimental to performance. To understand these results, we recall that the left entorhinal cortex and temporal pole contribute to familiarity for words  and binding perceptual and emotional information , respectively, and the orbitofrontal gyrus serves a key role in decision making, especially in the context of emotional and reward processing . Thus, our results suggest that dynamic brain trajectories driven more prominently by computations specialized for lexical recognition, emotional, and sensory processes may interfere with processes mediated through the cingulate and parietal systems.
It is interesting to consider these anatomical distributions of control-cognition relationships in the context of neuropsychiatric disorders. Specifically, in borderline personality disorder, a syndrome with substantial emotional processing and working memory deficits, orbitofrontal cortex dysfunction is associated with increased impulsive errors on attention demanding and working memory tasks  and altered temporal pole activity has been associated with social cognition and empathy dysfunction . An interesting future direction would be to test for alterations in structural controllability in these regions in clinical syndromes that present with emotion regulation and cognitive control dysfunction.
In this paper, we address the hypothesis that network control is a fundamental mechanism of cognitive control . We explicitly test this hypothesis by determining whether individual differences in putative control roles of brain regions are related to individual differences in performance on so-called cognitive control tasks. We construct structural brain networks from diffusion tensor imaging data acquired in 125 healthy adult individuals. We define a simplified model of brain dynamics and simulate network control to quantify modal and boundary controllability, which together describe complementary features of a region’s theoretically predicted preference to drive the brain into different cognitive states. We find that brain regions whose controllability values are associated with individual differences in cognitive control performance display non-trivial spatial clustering. Moreover, we find that distinct areas of the brain display controllability values that are strongly correlated with performance in a continuous performance test, a color-shape switching task, a Stroop task, and a n-back working memory task across individuals. Critically, these control-cognition relationships were two-sided: increases in controllability in some regions are advantageous for cognitive control performance, whereas increases in others are disadvantageous. These results provide the first evidence to support the hypothesis that human brains use network control to affect cognitive control.
Network Control Theory: a Framework to Study Cognition
Network control theory offers a novel framework for the study of human cognition. Unlike graph theory and its associated statistics [51, 52, 53], network control theory is built on a model of brain dynamics, which describes patterns of inter-region activity propagated along the wires of an underlying structural network . From such a model, one can mathematically derive theoretical predictions about how the structural network organization of the brain impacts on brain function . These tools can be tuned to study distinct types of control strategies. Here we focus on (i) modal controllability, which describes the role of a region in driving the brain into difficult-to-reach states, and (ii) boundary controllability, which describes the role of a region in driving the brain into states of network segregation or integration. This approach represents a substantial but complementary departure from traditional analyses of functional and diffusion imaging data at the voxel, tract, and graph levels. Instead, network controllability analysis can be used to better understand the role of each brain region in regulating whole-brain network function based explicitly on its structural fingerprint and a working model of node dynamics.
Relationship to Prior Theoretical Predictions
Prior work identified a spatial dissociation between different types of brain network controllers : modal controllers are prevalent in fronto-parietal cognitive control regions and boundary controllers are prevalent in dorsal and ventral attention regions. Based on this work, we anticipated that high levels of boundary controllability in anterior cingulate and opercular systems would be advantageous for attention. Further, we anticipated that high levels of modal controllability in frontoparietal systems would be associated with better working memory performance and task switching. Whereas we found strong support for the former hypothesis in the continuous performance test, we found weaker support for the latter hypothesis. Instead, we observed that individual differences in modal and boundary controllability in subcortical, default mode and lower association regions  were observed to positively relate to individual differences in cognitive control performance. These results suggest that while network controllers may be spatially clustered in certain brain subsystems , interindividual variability in regions that interact with frontoparietal systems through multiple pathways  may enhance brain-wide dynamics to support cognitive control functions.
Integrative Insights across Cognitive Control Tasks
In each cognitive task, we identified regions whose variation in controllability correlated with individual differences in performance. In some cases, these regions were consistent with those identified in activation-based studies as pertinent to the task, and in other cases, these regions were novel anatomical loci whose putative functions were supportive of the task. As an example of the former: the contributions of the anterior cingulate cortex in the CPT and N-Back tasks are consistent with classically defined roles for these regions in cognitive neuroscience. Our work extends this body of knowledge by suggesting that this region also plays a role in regulating the integration and segregation of network-wide dynamics, supporting its error monitoring, shifting and updating. As an example of the latter: the contributions of the basal ganglia in the color/shape task differ from the common emphasis on fronto-parietal and cingulo-opercular systems. Our results suggest that greater controllability in the putamen and pallidum (modal) as well as the nucleus accumbens (boundary) is supportive of better performance, potentially identifying a key role in regulating the brain via pathways involving cortico-basal ganglia-thalamo-cortical loops .
Competing Effects of Controllability on Cognition
Numerous functional tradeoffs exist in the brain . Cognitively, this includes the well-known oppositions between speed and accuracy , between value and time , and between exploration and exploitation . The approach we take in this work offers a novel dimension of dissociable effects on cognition: namely, that for some brain regions higher controllability values are benefitial for cognitive performance, while for others they are detrimental. While we do not yet know the fundamental mechanisms of these dissociations, it is interesting to speculate based on the theoretical intuitions behind the analytical approach. Specifically, in each task, the increasing controllability of regions may indicate dynamic interference from regions that do not contribute to the primary processing of stimuli or to the execution of network control that supports performance. It is also possible that higher controllability values in some regions may be detrimental to one cognitive function but supportive of another cognitive function, a possibility that could be examined directly in future empirical work.
Linking Cognitive Control and Network Control
From an engineering perspective, network control is a process in which a system’s dynamics are shifted or guided into new trajectories to support specific system goals. Such a process forms a natural description of cognitive control, a transient process putatively deployed to shift or guide the brain into new mental states in support of higher-order cognition. Prior work has reported that different brain regions are theoretically predicted to be effective at different control strategies: (i) regions in fronto-parietal cortex are suggested to be effective as modal controllers, moving the brain into difficult-to-reach states on an energy landscape, and (ii) regions in attentional circuits are suggested to be effective as boundary controllers, moving the brain into states in which cognitive subsystems or modules are either coupled or decoupled from one another . Our results provide the first empirical evidence that individual differences in structural controllability of single brain regions can explain significant variance in individual differences in cognitive control abilities, as measured by a battery of cognitive tests. These results support the possibility that network control might form a fundamental mechanism of cognitive control. This perspective complements other computational models of cognitive control  by highlighting the role of white matter architecture on large-scale functional dynamics.
Diffusion tensor imaging may undersample some white matter fibers, particularly those linking hemispheres or those that cross paths with other fibers . Future efforts could use diffusion spectrum imaging to improve estimates of structural network architecture. In addition, the cognitive organization of cognitive control remains contested, and it is possible that behavioral measures other than those studied here might elucidate other aspects of the relationship between cognitive control and network controllability. Furthermore, we note that the inter-measure behavioral correlations here are somewhat lower than those measured in prior factor analyses . While this allows us to speak to differences in cognitive control domains of function, future studies should investigate the generalizability of both the inter-measure relationships among cognitive control measures and robustness of associations with network controllability in samples of various socioeconomic backgrounds and ages. Finally, while we focus on modal and boundary controllability as an extension of previous research, other control strategies may also be relevant for the study of cognitive control in humans [14, 8, 63].
Our results link classic findings in cognitive neuroscience to the emerging field of network control theory. Moreover, they extend classical views of cognitive control by explicitly acknowledging the role of white matter organization in regulating the dynamic patterns of activity evolving on top of it. Future work could address the potential utility of these conceptual constructs in understanding deficits in cognitive control observed in neurological disease and psychiatric disorders.
Methods and any associated references are available in the online version of the paper.
Previous analyses involving diffusion tensor imaging tractography have found significant relationships with cognition in samples less than 20 individuals . However, the relationship between white matter network controllability and cognition has not been previously examined. To examine this relationship, we included the largest available sample of individuals with good quality diffusion images to examine our hypotheses. One hundred sixty six individuals completed a testing battery prior to a cognitive training intervention (study funding provided by NCI grant R01-CA170297) during daylight hours. Of these individuals, 145 had DTI images and a T1 acquired as part of a longer neuroimaging protocol. Using a criterion independently validated on the same scanner used to acquire these data, we excluded 17 individuals due to a signal to noise ratio of less than 6.47 . We removed a further three individuals for mean relative displacement greater than in a separate resting BOLD time series to reduce the likelihood that our results were influenced by subject motion. The final sample included 125 individuals (mean age , St.D. , range , females). All participants volunteered with written informed consent in accordance with the Institutional Review Board/Human Subjects Committee, University of Pennsylvania.
Data acquisition and structural network construction were performed as in previous work . DTI scans sampled 30 directions using a Q5 half-shell acquisition scheme with a maximum -value of 1,000 and an isotropic voxel size of 2.0 mm. We utilized an axial acquisition with the following parameters: repetition time (TR) = 8 s, echo time (TE)= 82 ms, 70 slices, field of view (FoV) (230, 230, 140 mm).
DTI data were reconstructed in DTI Studio (www.dsi-studio.labsolver.org) using -space diffeomorphic reconstruction (QSDR) . QSDR first reconstructs diffusion-weighted images in native space and computes the quantitative anisotropy (QA) in each voxel. These QA values are used to warp the brain to a template QA volume in Montreal Neurological Institute (MNI) space using the statistical parametric mapping (SPM) nonlinear registration algorithm. Once in MNI space, spin density functions were again reconstructed with a mean diffusion distance of 1.25 mm using three fiber orientations per voxel. Fiber tracking was performed in DSI Studio with an angular cutoff of 55, step size of 1.0 mm, minimum length of 10 mm, spin density function smoothing of 0.0, maximum length of 400 mm and a QA threshold determined by DWI signal in the colony-stimulating factor. Deterministic fiber tracking using a modified FACT algorithm was performed until 1,000,000 streamlines were reconstructed for each individual.
Anatomical scans were segmented using FreeSurfer  and parcellated using the connectome mapping toolkit . A parcellation scheme including regions was registered to the B0 volume from each subject’s DTI data. The B0 to MNI voxel mapping produced via QSDR was used to map region labels from native space to MNI coordinates. To extend region labels through the grey-white matter interface, the atlas was dilated by 4 mm . Dilation was accomplished by filling non-labelled voxels with the statistical mode of their neighbors’ labels. In the event of a tie, one of the modes was arbitrarily selected. Each streamline was labelled according to its terminal region pair. From these data, we constructed a structural connectivity matrix, whose element represented the number of streamlines connecting different regions, divided by the sum of volumes for regions and .
All cognitive measures were collected as part of a larger testing battery designed to assess cognitive control and decision making prior to a multi-week cognitive training program. We selected four measures based on (1) their statistical independence from one another and (2) their utility in predicting functional outcomes in health and disease. Specifically, measures included the continuous performance test number correct , color shape median switch cost (response times) , median Stroop effect (response times) , and a fractal (visuo-spatial) n-back [71, 19, 72]. All tasks were administered to participants seated at a desktop computer approximately two feet from the monitor. All stimuli were presented on a black background.
Continuous performance test.
The continuous performance test consisted of 360 trials of sequentially presented configurations of 2 to 7 red bars representing a digital display. Each trial was 1s in length with a fixed 700ms intertrial interval before the presentation of the next stimulus. One third of the stimuli (120 out of 360) formed numbers or letters. Participants were instructed to respond with a right index finger press of the spacebar on a standard keyboard when a number or letter was observed. No response was required for non-target trials. Responses made less than 100ms after the presentation of a stimulus were counted as incorrect due to inability to know if such responses were a reaction to the previous trial or impulsive response to the current trial .
Color shape task.
The color shape task consisted of two blocks of 48 trials (96 total) of sequentially presented figures with either a circle or triangle in the center flanked by either a red or green square surrounding the shape. Trial cues were presented 150ms before the stimulus, and both the cue and the stimulus remained on the screen until the participant responded. The cue “C” indicated that the participant should respond to the color (the “z” key with left index finger = green, the backslash key with right index finger = red), whereas an “S” indicated that the participant should respond to the shape (the “z” key with left index finger = circle, the backslash key with right index finger = triangle). Stimuli were separated by a 350ms blank screen and remained on screen until the participant responded .
The Stroop task consisted of 96 words printed on the screen sequentially in either red, green, or blue. Trials lasted until either the participant responded or 3.5 seconds elapsed and were separated by a fixed 100ms intertrial interval with a blank screen. In congruent trials, the printed word matched the color of the word. In incongruent trials, the color did not match the printed word. Participants were asked to respond only to the color of the word by pressing a colored button with one of the first three fingers on their right hand that matched the word color on each trial.
During the n-back, participants were instructed to remember the location of a stimulus, a grey circle (approximately 5 cm in diameter), as it appeared randomly in 8 possible locations around the perimeter of a computer screen. The stimulus appeared for 200 ms, followed by an interstimulus interval (ISI) of 2800 ms. A cross hair remained visible during the stimulus presentation to cue participants to look at the center of the screen so that all stimuli appearing around the perimeter of the screen could be seen clearly. The n-back task includes 4 conditions of 50 trials each of varying difficulty levels (200 trials total): the 0-back, 1-back, 2-back, and 3-back. During the 0-back, participants were asked to respond if the stimulus appeared in the upper left corner of the screen. During the 1-back condition, participants respond if the image is identical to the one preceding it. In the 2-back condition, they respond if the stimulus is identical to the one two trials before. In the 3-back condition they respond if the stimulus is identical to the one three trials before. Participants were asked to respond only to targets (30% of stimuli; 60 trials) by pressing the spacebar [71, 19, 72]. The primary outcome is the total number of true positives (correctly identifying a target) summed across all conditions.
See Table 1 for the mean and standard deviation of performance for the four cognitive measures across subjects. See Figure 7 for the inter-task correlations for the four measures.
|CPT(# correct)||Switch cost(ms)||Stroop(ms)||N-Back(#correct)|
|110.7 (9.6)||322.5 (253.5)||76.9 (89.5)||52.4 (4.9)|
Cognitive task performance is reported as mean (standard deviation) across all 125 participants. CPT maximum # correct = 120. N-Back maximum number correct = 60.
Network Control Theory
We follow a previous application of network control theory in diffusion weighted imaging data . While we briefly describe the mathematical approach here, we also point the reader to Ref.  for a more comprehensive description, and to Refs. [73, 12, 13] for additional mathematical details.
Our ability to understand neural systems is fundamentally related to our theoretical ability to control them . Network control theory is a branch of traditional control theory in engineering that examines how to control a system based on the links between its components; here, control means perturbing a system to reach a desired state. To apply a network control perspective, we require (i) knowledge of the network connectivity linking system components, and (ii) knowledge regarding how system components act, i.e., their dynamics. In contrast to traditional graph theory, network control theory offers mechanistic predictors of network dynamics. The use of mechanistic models offers the possibility for us to move from descriptive approaches in the human connectome toward an understanding of the structural and functional bases of cognition .
Mathematically, we can study the controllability of a networked system by defining a network represented by the graph , where and are the vertex and edge sets, respectively. Let be the weight associated with the edge , and define the weighted adjacency matrix of as , where whenever . We associate a real value (state) with each node, collect the node states into a vector (network state), and define the map to describe the evolution (network dynamics) of the network state over time. Given the network and node dynamics, we can use network control theory to quantitatively examine how the network structure predicts the types of control that nodes can exert.
Network Control Theory Applied to Neuroimaging
We begin with an analogous approach to prior work . We define structural brain networks by subdividing the entire brain into anatomically distinct brain areas (network nodes) in a commonly used anatomical atlas . Consistent with prior work [76, 77, 78, 8], we connect nodes by the number of white matter streamlines identified by a commonly used deterministic tractography algorithm (see  and above description). This procedure results in sparse, weighted, undirected structural brain networks for each subject (N = 125). The definition of structural brain networks based on tractography data in humans follows from our primary hypothesis that cognitive control processes are in part determined by the structural organization of the brain’s white matter tracts.
To define the dynamics of neural processes, we draw on prior models linking structural brain networks to resting state functional dynamics [80, 81, 82]. Although neural activity evolves through neural circuits as a collection of nonlinear dynamic processes, these prior studies have demonstrated that a significant amount of variance in neural dynamics as measured by fMRI can be predicted from simplified linear models [80, 83]. Based on this literature, we employ a simplified noise-free linear discrete-time and time-invariant network model:
where describes the state (e.g., a measure of the electrical charge, oxygen level, or firing rate) of brain regions over time, and is a symmetric and weighted adjacency matrix. In this case, we construct a weighted adjacency matrix whose elements indicate the number of white matter streamlines connecting two different brain regions – denoted here as and – and we stabilize this matrix by dividing by the mean edge weight. The diagonal elements of the matrix satisfy . The input matrix identifies the control points in the brain, where and
and denotes the -th canonical vector of dimension . The input denotes the control strategy.
To study the ability of a certain brain region to influence other regions in arbitrary ways we adopt the control theoretic notion of controllability. Controllability of a dynamical system refers to the possibility of driving the state of a dynamical system to a specific target state by means of an external control input . Classic results in control theory ensure that controllability of the network (1) from the set of network nodes is equivalent to the controllability Gramian being invertible, where
We utilize this framework to choose control nodes one at a time, and thus the input matrix in fact reduces to a one-dimensional vector.
Besides ensuring controllability, the eigenvalues of the controllability Gramian are a quantitative measure of the magnitude of the control input that drives a network to a desired target state , and the structure of the Gramian itself provides systematic guidelines for the selection of control areas that can theoretically optimize cognitive functions. While the magnitude of the control input may not be the unique feature to take into account when controlling brain dynamics , it allows us to better understand the relationship between the structural organization of the brain and its dynamics, and opens the door to the development of novel diagnostics. Here, this allows us to isolate the control role of each region separately in the context of the individual’s network, and then to associate variability in this property at each node with cognition.
Network Controllability Statistics
In this work, we study two types of controllability: modal controllability, and boundary controllability. Intuitively, modal controllability refers to the ability of a node to control each evolutionary mode of a dynamical network , and can be used to identify network configurations that are least controllable. Modal controllability is computed from the eigenvector matrix of the network adjacency matrix . By extension from the PBH test , if the entry is small, then the -th mode is poorly controllable from node . Following , we define as a scaled measure of the controllability of all modes from the brain region . Regions with high modal controllability are able to control all the dynamic modes of the network, and hence to drive the dynamics towards hard-to-reach configurations. Modal control values were assigned a ranked value to linearize their distribution .
To complement modal controllability, we also study boundary controllability, a metric developed in network control theory to quantify the role of a network node in controlling dynamics between modules in hierarchical networks. Boundary controllability identifies brain areas that can steer the system into states where different cognitive systems are either coupled or decoupled. Here, we apply a similar approach to that taken in  to quantify boundary controllability in our diffusion tractography networks and associate controllability variability with cognitive performance.
As boundary controllability describes a node’s role in controlling dynamics across modular network architecture, an initial identification of brain modules is required . To quantify the initial partition for the structural brain networks, we define the first level of the network empirically by maximizing the modularity quality function  using a Louvain-like  locally greedy algorithm . Our choice is based on extensive recent literature demonstrating that the brain is composed of many subnetworks (not just 2) [92, 11], which can be extracted using modularity maximization approaches [93, 94, 95, 96], and which correspond to sets of brain areas performing related functions [97, 94, 11].
The modularity quality function provides an estimate of the quality of a hard partition of the adjacency matrix into network communities (whereby each brain region is assigned to exactly one network community) [98, 99, 89, 100, 101]
where brain region is assigned to community , brain region is assigned to community , if and it equals otherwise, is a structural resolution parameter, and is the expected weight of the edge connecting node and node under a specified null model. Maximization of yields a hard partition of a network into communities such that the total edge weight inside of communities is as large as possible (relative to the null model and subject to the limitations of the employed computational heuristics, as optimizing is NP-hard [100, 101, 102]).
Because the modularity quality function has many near-degeneracies, it is important to perform the optimization algorithm multiple times . We perform 100 optimizations of the Louvain-like locally greedy algorithm  applied to the average structural adjacency matrix. To distill a single representative partition, we create a consensus partition from these 100 optimizations based on statistical comparison to an appropriate null model .
Importantly, we choose a value for the structural resolution parameter that produces a robust consensus partition over individuals. The choice is a common default decision, but it is important to consider multiple values of to examine the multiscale community structure of the A matrix [104, 100, 105]. Previous work has demonstrated that in some networks, a structural resolution parameter value that accurately captures the true community structure can be identified by the value at which the 100 optimizations produce similar partitions . To quantitatively estimate similarity in partitions, we adopt the -score of the Rand coefficient . For each pair of partitions and , we calculate the Rand -score in terms of the total number of pairs of nodes in the network , the number of pairs that are in the same community in partition , the number of pairs that are in the same community in partition , and the number of pairs of nodes that are assigned to the same community both in partition and in partition . The -score of the Rand coefficient comparing these two partitions is
where is the standard deviation of . Let the mean partition similarity denote the mean value of over all possible partition pairs for . Let the variance of the partition similarity denote the variance of over all possible partition pairs for .
Empirically, we calculated a group adjacency matrix by averaging the structural A matrices of all subjects. We optimized the modularity quality function 100 times and we computed the mean and variance of the partition similarity for a range of values. We observed that the mean partition similarity was high and the variance of the partition similarity was low for a value of at , which is within the range of stable partitions found in our prior analyses in diffusion spectrum imaging data . We therefore used the consens partition at for the remainder of the analysis in this study. To examine the fit of the group average consensus partition to each individual, we conducted a permutation test for each individual where the community assignments were randomly permuted 10,000 times per individual. Then, we computed the mean and variability of the -score of the Rand coefficient across all subjects and found that the subject similarity to the consensus partition was superior to randomly permuted assignments for all subjects (mean -Rand score = 55.4, standard deviation = 3.4 over all permutations for all subjects).
Boundary Point Criteria
The first modification concerns the definitions of the first level subnetworks for which we compute a two-partition based on the Fiedler eigenvector. Consistent with our prior work [66, 88] within the initial community partition we compute the Fiedler eigenvector of the adjacency matrix to create first level subnetworks defined by a two-partition. After calculting the two community partition, we must identify “boundary points”, which are nodes that contain connections to both communities. Following , we set a threshold ratio to identify boundary points. Considering the adaptivity to the local measure, we set a threshold ratio instead of a global threshold value. In detail, for a network with partition , a node is called a boundary node if
where is the adjacency matrix. Here, can be replaced with other statistics and needs to be chosen carefully. If is too small, there will be no effect of differences in edge weights within the subnetwork and the algorithm tends to add the total subnetwork as the set of boundary points. If is too large, there will be only a few points recognized as the boundary points.
In the results described in the main manuscript, we computed boundary controllability ranked over a range of to in increments of . The boundary control values are highly similar across choices of (See Fig. 8; minimum Pearson correlation approximately , corresponding to a , indicating that our results are robust to small variation in the boundary point criteria threshold. To use a stable measure of boundary controllability for association with cognitive variability, we compute boundary controllability for each node at each value of and take the average value of boundary controllability over this range.
Final Algorithm The final algorithm used in the calculation of boundary controllability in this paper can be summarized as follows. We begin with the application of a community detection method to the brain network to extract a partition of brain regions into network communities. We then recursively apply a Fiedler bipartition to add boundary nodes within communities, with the goal of improving the local controllability of the network. At each stage of the algorithm, we define the boundary nodes of the network as the nodes that maintain edges to nodes in other communities. Algorithmically, we can write:
[H] \SetAlgoLined\KwDataNetwork with adjacency matrix , Number of control nodes , threshold ratio ; \KwResultControl Nodes Index Set ; Define an empty set of control nodes Initialize the partition with the result of a community detection algorithm and initialize the boundary nodes set Add the boundary points of the initial partition \While Select least controllable community Compute Fiedler two-partition of -th community Compute boundary nodes of with the given threshold ratio Update partition with Update control nodes with boundary nodes \Return.
Associating Controllability with Cognition
To associate network controllability with measured cognitive performance, we computed the partial correlation between each cognitive variable and ranked controllability values controlling for age, sex, and DTI-scan total signal to noise values for each participant. Partial correlation coefficients were Fisher -to- transformed for visualization in the final maps. For each map, we tested for significant associations between the strength of controllability and geometric distance between regions. For each node pair in the adjacency matrix , we computed the Pearson’s correlation between the mean absolute value for the -transformed partial correlation coefficients and the Euclidean distance between region centers (two-sided test). The distance for region centers were transformed by empirically calculating the cumulative density function of observed distances and transforming this distribution to normal. Within the maps, nodes surpassing a Fisher’s value of +/- 2.0 were interpreted in the main manuscript.
Code used to calculate the controllability statistics described here can be found at: http://commdetect.weebly.com.
Data are available on request to the authors contingent upon IRB authorization.
JDM acknowledges support from the Office of the Director at the National Institutes of Health through grant number 1-DP5-OD-021352-01. DSB acknowledges support from the John D. and Catherine T. MacArthur Foundation, the Army Research Laboratory and the Army Research Office through contract numbers W911NF-10-2-0022 and W911NF-14-1-0679, the National Institute of Mental Health (2-R01-DC-009209-11), the National Institute of Child Health and Human Development (1R01HD086888-01), the Office of Naval Research, and the National Science Foundation (#BCS-1441502, #BCS-1430087, and #PHY-1554488). SG acknowledges support from the Applied Mathematics and Computational Science Graduate Program at the University of Pennsylvania. FP acknowledges support from the National Science Foundation award #BCS-1430279. CL acknowledges support from the National Cancer Institute (#R35-CA197461 and #R01-CA170297). The content is solely the responsibility of the authors and does not necessarily represent the official views of any of the funding agencies.
- S. Normandeau, F. Guay, Journal of Educational Psychology 90, 111 (1998).
- A. Diamond, W. S. Barnett, J. Thomas, S. Munro, Science (New York, NY) 318, 1387 (2007).
- M. W. Cole, et al., Nat Neurosci 16, 1348 (2013).
- J. D. Power, B. L. Schlaggar, C. N. Lessov-Schlaggar, S. E. Petersen, Neuron 79, 798 (2013).
- B. Voytek, R. T. Knight, Biological psychiatry 77, 1089 (2015).
- J. D. Medaglia, M. E. Lynall, D. S. Bassett, J Cogn Neurosci 27, 1471 (2015).
- L. Cocchi, A. Zalesky, A. Fornito, J. B. Mattingley, Trends in cognitive sciences 17, 493 (2013).
- S. Gu, et al., Nature Communications 6, 8414 (2015).
- M. G. Mattar, M. W. Cole, S. L. Thompson-Schill, D. S. Bassett, PLoS computational biology 11 (2015).
- U. Braun, et al., Proceedings of the National Academy of Sciences 112, 11678 (2015x).
- D. S. Bassett, M. Yang, N. F. Wymbs, S. T. Grafton, Nat Neurosci 18, 744 (2015).
- J. Ruths, D. Ruths, Science 343, 1373 (2014).
- F. Pasqualetti, S. Zampieri, F. Bullo, Control of Network Systems, IEEE Transactions on 1, 40 (2014).
- R. F. Betzel, S. Gu, J. D. Medaglia, F. Pasqualetti, D. S. Bassett, arXiv preprint arXiv:1603.05261 (2016).
- A. Miyake, et al., Cognitive psychology 41, 49 (2000).
- M. M. Kurtz, J. D. Ragland, W. Bilker, R. C. Gur, R. E. Gur, Schizophrenia research 48, 307 (2001).
- A. Miyake, M. J. Emerson, F. Padilla, J.-c. Ahn, Acta psychologica 115, 123 (2004).
- J. R. Stroop, Journal of experimental psychology 18, 643 (1935).
- A. Green, et al., Pharmacology Biochemistry and Behavior 81, 575 (2005).
- R. J. Ogg, et al., Magnetic resonance imaging 26, 504 (2008).
- R. Leech, D. J. Sharp, Brain 137, 12 (2014).
- A. W. Sali, S. M. Courtney, S. Yantis, The Journal of Neuroscience 36, 445 (2016).
- M. M. Botvinick, J. D. Cohen, C. S. Carter, Trends in cognitive sciences 8, 539 (2004).
- Y. Nagahama, et al., Neuroimage 10, 193 (1999).
- A. B. Smith, E. Taylor, M. Brammer, B. Toone, K. Rubia, American Journal of Psychiatry (2006).
- P. C. Holland, M. Gallagher, Trends in cognitive sciences 3, 65 (1999).
- S. Shanmugan, et al., American Journal of Psychiatry (2016).
- M. Isoda, O. Hikosaka, Nature neuroscience 10, 240 (2007).
- O. Hikosaka, M. Isoda, Trends in cognitive sciences 14, 154 (2010).
- S. B. Floresco, S. Ghods-Sharifi, C. Vexelman, O. Magyar, The Journal of neuroscience 26, 2449 (2006).
- K. Garner, P. E. Dux, Proceedings of the National Academy of Sciences 112, 14372 (2015).
- A. R. Aron, et al., The Journal of Neuroscience 27, 11860 (2007).
- A. Parent, L.-N. Hazrati, Brain Research Reviews 20, 91 (1995).
- R. Cools, R. A. Barker, B. J. Sahakian, T. W. Robbins, Brain 124, 2503 (2001).
- N. Derakshan, S. Smyth, M. W. Eysenck, Psychonomic Bulletin & Review 16, 1112 (2009).
- C. Gross, R. Hen, Nature Reviews Neuroscience 5, 545 (2004).
- A. Puce, T. Allison, M. Asgari, J. C. Gore, G. McCarthy, The Journal of Neuroscience 16, 5205 (1996).
- P. Ludersdorfer, M. Kronbichler, H. Wimmer, Human brain mapping 36, 1393 (2015).
- S. H. Creem, D. R. Proffitt, Acta psychologica 107, 43 (2001).
- M. T. Banich, et al., Cognitive Neuroscience, Journal of 12, 988 (2000).
- M. Koenigs, A. K. Barbey, B. R. Postle, J. Grafman, The Journal of Neuroscience 29, 14980 (2009).
- R. Vandenberghe, P. Molenberghs, C. R. Gillebert, Neuropsychologia 50, 1092 (2012).
- R. M. Braga, D. J. Sharp, C. Leeson, R. J. Wise, R. Leech, The Journal of Neuroscience 33, 14031 (2013).
- M. Petrides, Executive Control and the Frontal Lobe: Current Issues (Springer, 2000), pp. 44–54.
- J. D. Power, et al., Neuron 72, 665 (2011).
- K. R. Brandt, M. W. Eysenck, M. K. Nielsen, T. J. von Oertzen, Brain and cognition 104, 82 (2016).
- I. R. Olson, A. Plotzker, Y. Ezzyat, Brain 130, 1718 (2007).
- A. Bechara, H. Damasio, A. R. Damasio, Cerebral cortex 10, 295 (2000).
- H. A. Berlin, E. T. Rolls, S. D. Iversen, American journal of psychiatry 162, 2360 (2005).
- I. Dziobek, et al., Neuroimage 57, 539 (2011).
- E. Bullmore, O. Sporns, Nature Review Neuroscience 10, 186 (2009).
- M. Rubinov, O. Sporns, NeuroImage 52, 1059 (2010).
- E. T. Bullmore, D. S. Bassett, Annual Reviews of Clinical Psychology 7, 113 (2011).
- J. Sepulcre, Neuroscience Letters 567, 68 (2014).
- P. J. Hellyer, et al., Journal of Neuroscience 34, 451 (2014).
- S.-H. Wang, Brain, behavior and evolution 72, 159 (2008).
- L. Chittka, P. Skorupski, N. E. Raine, Trends in Ecology & Evolution 24, 400 (2009).
- S. Pinker, Annals of the New York Academy of Sciences 882, 119 (1999).
- J. D. Cohen, S. M. McClure, J. Y. Angela, Philosophical Transactions of the Royal Society of London B: Biological Sciences 362, 933 (2007).
- M. M. Botvinick, J. D. Cohen, Cognitive science 38, 1249 (2014).
- V. J. Wedeen, et al., Neuroimage 41, 1267 (2008).
- A. Miyake, N. P. Friedman, D. A. Rettinger, P. Shah, M. Hegarty, Journal of experimental psychology: General 130, 621 (2001).
- A. E. Motter, Chaos: An Interdisciplinary Journal of Nonlinear Science 25, 097621 (2015).
- M. F. Kraus, et al., Brain 130, 2508 (2007).
- D. R. Roalf, et al., NeuroImage 125, 903 (2016).
- S. Gu, et al., Nature communications 6 (2015).
- F.-C. Yeh, V. J. Wedeen, W.-Y. I. Tseng, Neuroimage 55, 1054 (2011).
- B. Fischl, Neuroimage 62, 774 (2012).
- L. Cammoun, et al., Journal of neuroscience methods 203, 386 (2012).
- M. Cieslak, S. Grafton, Brain imaging and behavior 8, 292 (2014).
- A.-C. Ehlis, C. G. Bähne, C. P. Jacob, M. J. Herrmann, A. J. Fallgatter, Journal of psychiatric research 42, 1060 (2008).
- A. M. Owen, K. M. McMillan, A. R. Laird, E. Bullmore, Human brain mapping 25, 46 (2005).
- Y.-Y. Liu, J.-J. Slotine, A.-L. Barabási, Nature 473, 167 (2011).
- S. J. Schiff, Neural Control Engineering (The MIT Press, 2012).
- P. Hagmann, et al., PLoS Biology 6, e159 (2008).
- D. S. Bassett, J. A. Brown, V. Deshpande, J. M. Carlson, S. T. Grafton, Neuroimage 54, 1262 (2011).
- A. M. Hermundstad, et al., Proc Natl Acad Sci U S A 110, 6169 (2013).
- A. M. Hermundstad, et al., PLoS Comput Biol 10, e1003591 (2014).
- M. Cieslak, S. T. Grafton, Brain Imaging Behav 8, 292 (2014).
- C. Honey, et al., Proceedings of the National Academy of Sciences 106, 2035 (2009).
- C. J. Honey, J.-P. Thivierge, O. Sporns, Neuroimage 52, 766 (2010).
- F. Abdelnour, H. U. Voss, A. Raj, Neuroimage 90, 335 (2014).
- R. Fernández Galán, PLoS One 3, e2148 (2008).
- R. E. Kalman, Y. C. Ho, S. K. Narendra, Contributions to Differential Equations 1, 189 (1963).
- T. Kailath, Linear Systems (Prentice-Hall, 1980).
- G. Kumar, D. Menolascino, S. Ching, arXiv preprint arXiv:1411.5892 (2014).
- A. M. A. Hamdan, A. H. Nayfeh, AIAA Journal of Guidance, Control, and Dynamics 12, 421 (1989).
- F. Pasqualetti, S. Zampieri, F. Bullo, IEEE Transactions on Control of Network Systems 1, 40 (2014).
- M. E. Newman, Proc Natl Acad Sci U S A 103, 8577 (2006).
- V. D. Blondel, J.-L. Guillaume, R. Lambiotte, E. Lefebvre, Journal of Statistical Mechanics: Theory and Experiment 10, P1000 (2008).
- I. S. Jutla, L. G. S. Jeub, P. J. Mucha, A generalized Louvain method for community detection implemented in MATLAB (2011–2012).
- D. Meunier, R. Lambiotte, E. T. Bullmore, Front Neurosci 4, 200 (2010).
- D. Meunier, S. Achard, A. Morcom, E. Bullmore, Neuroimage 44, 715 (2009).
- Z. J. Chen, Y. He, P. Rosa-Neto, J. Germann, A. C. Evans, Cereb Cortex 18, 2374 (2008).
- D. S. Bassett, F. Siebenhuhner, In Multiscale Analysis and Nonlinear Dynamics: From Genes to the Brain (Wiley, 2013), chap. Multiscale network organization in the human brain.
- D. S. Bassett, et al., Proc Natl Acad Sci U S A 108, 7641 (2011).
- J. D. Power, et al., Neuron 72, 665 (2011).
- M. E. Newman, M. Girvan, Phys Rev E 69, 026113 (2004).
- M. E. Newman, Phys Rev E 69, 066133 (2004).
- M. A. Porter, J.-P. Onnela, P. J. Mucha, Notices of the American Mathematical Society 56, 1082 (2009).
- S. Fortunato, Phys Rep 486, 75 (2010).
- U. Brandes, et al., IEEE Trans on Knowl Data Eng 20, 172 (2008).
- B. H. Good, Y. A. de Montjoye, A. Clauset, Phys Rev E 81, 046106 (2010).
- J. Reichardt, S. Bornholdt, Phys Rev E 74, 016110 (2006).
- J.-P. Onnela, et al., Phys Rev E 86, 036104 (2012).
- A. L. Traud, E. D. Kelsic, P. J. Mucha, M. A. Porter, SIAM Review 53, 526 (2011).