Macroscopic resting-state brain dynamics are best described by linear models.

It is typically assumed that large networks of neurons exhibit a large repertoire of nonlinear behaviours. Here we challenge this assumption by leveraging mathematical models derived from measurements of local field potentials via intracranial electroencephalography and of whole-brain blood-oxygen-level-dependent brain activity via functional magnetic resonance imaging. We used state-of-the-art linear and nonlinear families of models to describe spontaneous resting-state activity of 700 participants in the Human Connectome Project and 122 participants in the Restoring Active Memory project. We found that linear autoregressive models provide the best fit across both data types and three performance metrics: predictive power, computational complexity and the extent of the residual dynamics unexplained by the model. To explain this observation, we show that microscopic nonlinear dynamics can be counteracted or masked by four factors associated with macroscopic dynamics: averaging over space and over time, which are inherent to aggregated macroscopic brain activity, and observation noise and limited data samples, which stem from technological limitations. We therefore argue that easier-to-interpret linear models can faithfully describe macroscopic brain dynamics during resting-state conditions.

Improved Seizure Onset-Zone Lateralization in Temporal Lobe Epilepsy using 7T Resting-State fMRI: A Direct Comparison with 3T.

Objective: Resting-state functional magnetic resonance imaging (rs-fMRI) at ultra high-field strengths (≥7T) is known to provide superior signal-to-noise and statistical power than comparable acquisitions at lower field strengths. In this study, we aim to provide a direct comparison of the seizure onset-zone (SOZ) lateralizing ability of 7T rs-fMRI and 3T rs-fMRI. Methods: We investigated a cohort of 70 temporal lobe epilepsy (TLE) patients. A paired cohort of 19 patients had 3T and 7T rs-fMRI acquisitions for direct comparison between the two field strengths. Forty-three patients had only 3T, and 8 patients had only 7T rs-fMRI acquisitions. We quantified the functional connectivity between the hippocampus and other nodes within the default mode network (DMN) using seed-to-voxel connectivity, and measured how hippocampo-DMN connectivity could inform SOZ lateralization at 7T and 3T field strengths. Results: Differences between hippocampo-DMN connectivity ipsilateral and contralateral to the SOZ were significantly higher at 7T (pFDR=0.008) than at 3T (pFDR=0.80) when measured in the same subjects. We found that our ability to lateralize the SOZ, by distinguishing subjects with left TLE from subjects with right TLE, was superior at 7T (AUC = 0.97) than 3T (AUC = 0.68). Our findings were reproduced in extended cohorts of subjects scanned at either 3T or 7T. Our rs-fMRI findings at 7T, but not 3T, are consistent and highly correlated (Spearman Rho=0.65) with clinical FDG-PET lateralizing hypometabolism. Significance: We show superior SOZ lateralization in TLE patients when using 7T relative to 3T rs-fMRI, supporting the adoption of high-field strength functional imaging in the epilepsy presurgical evaluation.

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study.

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

Subcortical functional connectivity gradients in temporal lobe epilepsy.

BACKGROUND AND MOTIVATION: Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE. METHODS: In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 24 R-TLE patients and 31 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures. RESULTS: We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE. CONCLUSION: Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

Resting state functional connectivity demonstrates increased segregation in bilateral temporal lobe epilepsy.

OBJECTIVE: Temporal lobe epilepsy (TLE) is the most common type of focal epilepsy. An increasingly identified subset of patients with TLE consists of those who show bilaterally independent temporal lobe seizures. The purpose of this study was to leverage network neuroscience to better understand the interictal whole brain network of bilateral TLE (BiTLE). METHODS: In this study, using a multicenter resting state functional magnetic resonance imaging (rs-fMRI) data set, we constructed whole-brain functional networks of 19 patients with BiTLE, and compared them to those of 75 patients with unilateral TLE (UTLE). We quantified resting-state, whole-brain topological properties using metrics derived from network theory, including clustering coefficient, global efficiency, participation coefficient, and modularity. For each metric, we computed an average across all brain regions, and iterated this process across network densities. Curves of network density vs each network metric were compared between groups. Finally, we derived a combined metric, which we term the "integration-segregation axis," by combining whole-brain average clustering coefficient and global efficiency curves, and applying principal component analysis (PCA)-based dimensionality reduction. RESULTS: Compared to UTLE, BiTLE had decreased global efficiency (p = .031), and decreased whole brain average participation coefficient across a range of network densities (p = .019). Modularity maximization yielded a larger number of smaller communities in BiTLE than in UTLE (p = .020). Differences in network properties separate BiTLE and UTLE along the integration-segregation axis, with regions within the axis having a specificity of up to 0.87 for BiTLE. Along the integration-segregation axis, UTLE patients with poor surgical outcomes were distributed in the same regions as BiTLE, and network metrics confirmed similar patterns of increased segregation in both BiTLE and poor outcome UTLE. SIGNIFICANCE: Increased interictal whole-brain network segregation, as measured by rs-fMRI, is specific to BiTLE, as well as poor surgical outcome UTLE, and may assist in non-invasively identifying this patient population prior to intracranial electroencephalography or device implantation.

Quantifying trial-by-trial variability during cortico-cortical evoked potential mapping of epileptogenic tissue.

OBJECTIVE: Measuring cortico-cortical evoked potentials (CCEPs) is a promising tool for mapping epileptic networks, but it is not known how variability in brain state and stimulation technique might impact the use of CCEPs for epilepsy localization. We test the hypotheses that (1) CCEPs demonstrate systematic variability across trials and (2) CCEP amplitudes depend on the timing of stimulation with respect to endogenous, low-frequency oscillations. METHODS: We studied 11 patients who underwent CCEP mapping after stereo-electroencephalography electrode implantation for surgical evaluation of drug-resistant epilepsy. Evoked potentials were measured from all electrodes after each pulse of a 30 s, 1 Hz bipolar stimulation train. We quantified monotonic trends, phase dependence, and standard deviation (SD) of N1 (15-50 ms post-stimulation) and N2 (50-300 ms post-stimulation) amplitudes across the 30 stimulation trials for each patient. We used linear regression to quantify the relationship between measures of CCEP variability and the clinical seizure-onset zone (SOZ) or interictal spike rates. RESULTS: We found that N1 and N2 waveforms exhibited both positive and negative monotonic trends in amplitude across trials. SOZ electrodes and electrodes with high interictal spike rates had lower N1 and N2 amplitudes with higher SD across trials. Monotonic trends of N1 and N2 amplitude were more positive when stimulating from an area with higher interictal spike rate. We also found intermittent synchronization of trial-level N1 amplitude with low-frequency phase in the hippocampus, which did not localize the SOZ. SIGNIFICANCE: These findings suggest that standard approaches for CCEP mapping, which involve computing a trial-averaged response over a .2-1 Hz stimulation train, may be masking inter-trial variability that localizes to epileptogenic tissue. We also found that CCEP N1 amplitudes synchronize with ongoing low-frequency oscillations in the hippocampus. Further targeted experiments are needed to determine whether phase-locked stimulation could have a role in localizing epileptogenic tissue.

Mindful attention promotes control of brain network dynamics for self-regulation and discontinues the past from the present.

Mindful attention is characterized by acknowledging the present experience as a transient mental event. Early stages of mindfulness practice may require greater neural effort for later efficiency. Early effort may self-regulate behavior and focalize the present, but this understanding lacks a computational explanation. Here we used network control theory as a model of how external control inputs-operationalizing effort-distribute changes in neural activity evoked during mindful attention across the white matter network. We hypothesized that individuals with greater network controllability, thereby efficiently distributing control inputs, effectively self-regulate behavior. We further hypothesized that brain regions that utilize greater control input exhibit shorter intrinsic timescales of neural activity. Shorter timescales characterize quickly discontinuing past processing to focalize the present. We tested these hypotheses in a randomized controlled study that primed participants to either mindfully respond or naturally react to alcohol cues during fMRI and administered text reminders and measurements of alcohol consumption during 4 wk postscan. We found that participants with greater network controllability moderated alcohol consumption. Mindful regulation of alcohol cues, compared to one's own natural reactions, reduced craving, but craving did not differ from the baseline group. Mindful regulation of alcohol cues, compared to the natural reactions of the baseline group, involved more-effortful control of neural dynamics across cognitive control and attention subnetworks. This effort persisted in the natural reactions of the mindful group compared to the baseline group. More-effortful neural states had shorter timescales than less effortful states, offering an explanation for how mindful attention promotes being present.

Subcortical Functional Connectivity Gradients in Temporal Lobe Epilepsy.

Background and Motivation: Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE. Methods: In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 23 R-TLE patients and 32 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures. Results: We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE. Conclusion: Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

Information content of brain states is explained by structural constraints on state energetics.

Signal propagation along the structural connectome of the brain induces changes in the patterns of activity. These activity patterns define global brain states and contain information in accordance with their expected probability of occurrence. Being the physical substrate upon which information propagates, the structural connectome, in conjunction with the dynamics, determines the set of possible brain states and constrains the transition between accessible states. Yet, precisely how these structural constraints on state transitions relate to their information content remains unexplored. To address this gap in knowledge, we defined the information content as a function of the activation distribution, where statistically rare values of activation correspond to high information content. With this numerical definition in hand, we studied the spatiotemporal distribution of information content in functional magnetic resonance imaging (fMRI) data from the Human Connectome Project during different tasks, and report four key findings. First, information content strongly depends on cognitive context; its absolute level and spatial distribution depend on the cognitive task. Second, while information content shows similarities to other measures of brain activity, it is distinct from both Neurosynth maps and task contrast maps generated by a general linear model applied to the fMRI data. Third, the brain's structural wiring constrains the cost to control its state, where the cost to transition into high information content states is larger than that to transition into low information content states. Finally, all state transitions-especially those to high information content states-are less costly than expected from random network null models, thereby indicating the brains marked efficiency. Taken together, our findings establish an explanatory link between the information contained in a brain state and the energetic cost of attaining that state, thereby laying important groundwork for our understanding of large-scale cognitive computations.

Altered functional brain dynamics in chromosome 22q11.2 deletion syndrome during facial affect processing.

Chromosome 22q11.2 deletion syndrome (22q11.2DS) is a multisystem disorder associated with multiple congenital anomalies, variable medical features, and neurodevelopmental differences resulting in diverse psychiatric phenotypes, including marked deficits in facial memory and social cognition. Neuroimaging in individuals with 22q11.2DS has revealed differences relative to matched controls in BOLD fMRI activation during facial affect processing tasks. However, time-varying interactions between brain areas during facial affect processing have not yet been studied with BOLD fMRI in 22q11.2DS. We applied constrained principal component analysis to identify temporally overlapping brain activation patterns from BOLD fMRI data acquired during an emotion identification task from 58 individuals with 22q11.2DS and 58 age-, race-, and sex-matched healthy controls. Delayed frontal-motor feedback signals were diminished in individuals with 22q11.2DS, as were delayed emotional memory signals engaging amygdala, hippocampus, and entorhinal cortex. Early task-related engagement of motor and visual cortices and salience-related insular activation were relatively preserved in 22q11.2DS. Insular activation was associated with task performance within the 22q11.2DS sample. Differences in cortical surface area, but not cortical thickness, showed spatial alignment with an activation pattern associated with face processing. These findings suggest that relative to matched controls, primary visual processing and insular function are relatively intact in individuals with 22q11.22DS, while motor feedback, face processing, and emotional memory processes are more affected. Such insights may help inform potential interventional targets and enhance the specificity of neuroimaging indices of cognitive dysfunction in 22q11.2DS.

Broken detailed balance and entropy production in the human brain.

Living systems break detailed balance at small scales, consuming energy and producing entropy in the environment to perform molecular and cellular functions. However, it remains unclear how broken detailed balance manifests at macroscopic scales and how such dynamics support higher-order biological functions. Here we present a framework to quantify broken detailed balance by measuring entropy production in macroscopic systems. We apply our method to the human brain, an organ whose immense metabolic consumption drives a diverse range of cognitive functions. Using whole-brain imaging data, we demonstrate that the brain nearly obeys detailed balance when at rest, but strongly breaks detailed balance when performing physically and cognitively demanding tasks. Using a dynamic Ising model, we show that these large-scale violations of detailed balance can emerge from fine-scale asymmetries in the interactions between elements, a known feature of neural systems. Together, these results suggest that violations of detailed balance are vital for cognition and provide a general tool for quantifying entropy production in macroscopic systems.

Computational modeling of tau pathology spread reveals patterns of regional vulnerability and the impact of a genetic risk factor.

Neuropathological staging studies have suggested that tau pathology spreads through the brain in Alzheimer's disease (AD) and other tauopathies, but it is unclear how neuroanatomical connections, spatial proximity, and regional vulnerability contribute. In this study, we seed tau pathology in the brains of nontransgenic mice with AD tau and quantify pathology development over 9 months in 134 brain regions. Network modeling of pathology progression shows that diffusion through the connectome is the best predictor of tau pathology patterns. Further, deviations from pure neuroanatomical spread are used to estimate regional vulnerability to tau pathology and identify related gene expression patterns. Last, we show that pathology spread is altered in mice harboring a mutation in leucine-rich repeat kinase 2. While tau pathology spread is still constrained by anatomical connectivity in these mice, it spreads preferentially in a retrograde direction. This study provides a framework for understanding neuropathological progression in tauopathies.

Modeling brain, symptom, and behavior in the winds of change.

Neuropsychopharmacology addresses pressing questions in the study of three intertwined complex systems: the brain, human behavior, and symptoms of illness. The field seeks to understand the perturbations that impinge upon those systems, either driving greater health or illness. In the pursuit of this aim, investigators often perform analyses that make certain assumptions about the nature of the systems that are being perturbed. Those assumptions can be encoded in powerful computational models that serve to bridge the wide gulf between a descriptive analysis and a formal theory of a system's response. Here we review a set of three such models along a continuum of complexity, moving from a local treatment to a network treatment: one commonly applied form of the general linear model, impulse response models, and network control models. For each, we describe the model's basic form, review its use in the field, and provide a frank assessment of its relative strengths and weaknesses. The discussion naturally motivates future efforts to interlink data analysis, computational modeling, and formal theory. Our goal is to inspire practitioners to consider the assumptions implicit in their analytical approach, align those assumptions to the complexity of the systems under study, and take advantage of exciting recent advances in modeling the relations between perturbations and system function.

Defining and predicting transdiagnostic categories of neurodegenerative disease.

The prevalence of concomitant proteinopathies and heterogeneous clinical symptoms in neurodegenerative diseases hinders the identification of individuals who might be candidates for a particular intervention. Here, by applying an unsupervised clustering algorithm to post-mortem histopathological data from 895 patients with degeneration in the central nervous system, we show that six non-overlapping disease clusters can simultaneously account for tau neurofibrillary tangles, α-synuclein inclusions, neuritic plaques, inclusions of the transcriptional repressor TDP-43, angiopathy, neuron loss and gliosis. We also show that membership to the six transdiagnostic disease clusters, which explains more variance in cognitive phenotypes than can be explained by individual diagnoses, can be accurately predicted from scores of the Mini-Mental Status Exam, protein levels in cerebrospinal fluid, and genotype at the APOE and MAPT loci, via cross-validated multiple logistic regression. This combination of unsupervised and supervised data-driven tools provides a framework that could be used to identify latent disease subtypes in other areas of medicine.

Gene coexpression patterns predict opiate-induced brain-state transitions.

Opioid addiction is a chronic, relapsing disorder associated with persistent changes in brain plasticity. Reconfiguration of neuronal connectivity may explain heightened abuse liability in individuals with a history of chronic drug exposure. To characterize network-level changes in neuronal activity induced by chronic opiate exposure, we compared FOS expression in mice that are morphine-naïve, morphine-dependent, or have undergone 4 wk of withdrawal from chronic morphine exposure, relative to saline-exposed controls. Pairwise interregional correlations in FOS expression data were used to construct network models that reveal a persistent reduction in connectivity strength following opiate dependence. Further, we demonstrate that basal gene expression patterns are predictive of changes in FOS correlation networks in the morphine-dependent state. Finally, we determine that regions of the hippocampus, striatum, and midbrain are most influential in driving transitions between opiate-naïve and opiate-dependent brain states using a control theoretic approach. This study provides a framework for predicting the influence of specific therapeutic interventions on the state of the opiate-dependent brain.

Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands.

A diverse set of white matter connections supports seamless transitions between cognitive states. However, it remains unclear how these connections guide the temporal progression of large-scale brain activity patterns in different cognitive states. Here, we analyze the brain's trajectories across a set of single time point activity patterns from functional magnetic resonance imaging data acquired during the resting state and an n-back working memory task. We find that specific temporal sequences of brain activity are modulated by cognitive load, associated with age, and related to task performance. Using diffusion-weighted imaging acquired from the same subjects, we apply tools from network control theory to show that linear spread of activity along white matter connections constrains the probabilities of these sequences at rest, while stimulus-driven visual inputs explain the sequences observed during the n-back task. Overall, these results elucidate the structural underpinnings of cognitively and developmentally relevant spatiotemporal brain dynamics.

Glucocerebrosidase Activity Modulates Neuronal Susceptibility to Pathological α-Synuclein Insult.

Mutations in the GBA1 gene are the most common genetic risk factor for Parkinson's disease (PD) and dementia with Lewy bodies (DLB). GBA1 encodes the lysosomal lipid hydrolase glucocerebrosidase (GCase), and its activity has been linked to accumulation of α-synuclein. The current study systematically examines the relationship between GCase activity and both pathogenic and non-pathogenic forms of α-synuclein in primary hippocampal, cortical, and midbrain neuron and astrocyte cultures, as well as in transgenic mice and a non-transgenic mouse model of PD. We find that reduced GCase activity does not result in aggregation of α-synuclein. However, in the context of extant misfolded α-synuclein, GCase activity modulates neuronal susceptibility to pathology. Furthermore, this modulation does not depend on neuron type but rather is driven by the level of pathological α-synuclein seeds. This study has implications for understanding how GBA1 mutations influence PD pathogenesis and provides a platform for testing novel therapeutics.

Spread of α-synuclein pathology through the brain connectome is modulated by selective vulnerability and predicted by network analysis.

Studies of patients afflicted by neurodegenerative diseases suggest that misfolded proteins spread through the brain along anatomically connected networks, prompting progressive decline. Recently, mouse models have recapitulated the cell-to-cell transmission of pathogenic proteins and neuron death observed in patients. However, the factors regulating the spread of pathogenic proteins remain a matter of debate due to an incomplete understanding of how vulnerability functions in the context of spread. Here we use quantitative pathology mapping in the mouse brain, combined with network modeling to understand the spatiotemporal pattern of spread. Patterns of α-synuclein pathology are well described by a network model that is based on two factors: anatomical connectivity and endogenous α-synuclein expression. The map and model allow the assessment of selective vulnerability to α-synuclein pathology development and neuron death. Finally, we use quantitative pathology to understand how the G2019S LRRK2 genetic risk factor affects the spread and toxicity of α-synuclein pathology.

Harnessing networks and machine learning in neuropsychiatric care.

The development of next-generation therapies for neuropsychiatric illness will likely rely on a precise and accurate understanding of human brain dynamics. Toward this end, researchers have focused on collecting large quantities of neuroimaging data. For simplicity, we will refer to large cross-sectional neuroimaging studies as broad studies and to intensive longitudinal studies as deep studies. Recent progress in identifying illness subtypes and predicting treatment response in neuropsychiatry has been supported by these study designs, along with methods bridging machine learning and network science. Such methods combine analytic power, interpretability, and direct connection to underlying theory in cognitive neuroscience. Ultimately, we propose a general framework for the treatment of neuropsychiatric illness relying on the findings from broad and deep studies combined with basic cognitive and physiologic measurements.

Sex differences in network controllability as a predictor of executive function in youth.

Executive function is a quintessential human capacity that emerges late in development and displays different developmental trends in males and females. Sex differences in executive function in youth have been linked to vulnerability to psychopathology as well as to behaviors that impinge on health, wellbeing, and longevity. Yet, the neurobiological basis of these differences is not well understood, in part due to the spatiotemporal complexity inherent in patterns of brain network maturation supporting executive function. Here we test the hypothesis that sex differences in impulsivity in youth stem from sex differences in the controllability of structural brain networks as they rewire over development. Combining methods from network neuroscience and network control theory, we characterize the network control properties of structural brain networks estimated from diffusion imaging data acquired in males and females in a sample of 879 youth aged 8-22 years. We summarize the control properties of these networks by estimating average and modal controllability, two statistics that probe the ease with which brain areas can drive the network towards easy versus difficult-to-reach states. We find that females have higher modal controllability in frontal, parietal, and subcortical regions while males have higher average controllability in frontal and subcortical regions. Furthermore, controllability profiles in males are negatively related to the false positive rate on a continuous performance task, a common measure of impulsivity. Finally, we find associations between average controllability and individual differences in activation during an n-back working memory task. Taken together, our findings support the notion that sex differences in the controllability of structural brain networks can partially explain sex differences in executive function. Controllability of structural brain networks also predicts features of task-relevant activation, suggesting the potential for controllability to represent context-specific constraints on network state more generally.