The human brain is a complex system with high metabolic demands and extensive connectivity that requires control to balance energy consumption and functional efficiency over time. How this control is manifested on a whole-brain scale is largely unexplored, particularly what the associated costs are. Using the network control theory, here, we introduce a novel concept, time-averaged control energy (TCE), to quantify the cost of controlling human brain dynamics at rest, as measured from functional and diffusion MRI. Importantly, TCE spatially correlates with oxygen metabolism measures from the positron emission tomography, providing insight into the bioenergetic footing of resting-state control. Examining the temporal dimension of control costs, we find that brain state transitions along a hierarchical axis from sensory to association areas are more efficient in terms of control costs and more frequent within hierarchical groups than between. This inverse correlation between temporal control costs and state visits suggests a mechanism for maintaining functional diversity while minimizing energy expenditure. By unpacking the temporal dimension of control costs, we contribute to the neuroscientific understanding of how the brain governs its functionality while managing energy expenses.

Understanding how the brain balances functional efficiency with energy conservation is a central question in neuroscience. The network control theory (NCT) views this question from a network perspective where the brain manages signal propagations along its structural connections to transition across desired activity states. Our study thus presents a novel framework based on the NCT to analyze the costs associated with transitioning across resting states, revealing that regions with high control costs on average are also metabolically demanding in terms of oxygen use. Our findings further show that transitions between sensory and association states are infrequent due to high control costs, while transitions within these states are more common. This suggests that the brain employs a mechanism to preserve functional diversity while minimizing energy costs.

The intricate nature of the human brain, characterized by its extensive connectivity (Bassett & Sporns, 2017; Bullmore & Sporns, 2012; Thiebaut de Schotten & Forkel, 2022) and correspondingly high metabolic demands (Padamsey & Rochefort, 2023; Raichle & Gusnard, 2002), has led to the assumption that its inherent organization represents a delicate balance between energy consumption and functional efficiency (Shine & Poldrack, 2018). This notion underscores the role of the nervous system as a regulatory entity (Friston, 2010), orchestrating its functions to fulfill cognitive demands while simultaneously managing energy-related constraints. While prior research has shed light on how the brain manages its energy expenditures at a cellular level (Gautron, Elmquist, & Williams, 2015; Levin, 2006; Padamsey & Rochefort, 2023), the broader exploration of how these mechanisms translate to the regulation of functional costs at a whole-brain scale remains relatively uncharted.

Drawing inspiration from engineering principles, network control theory (NCT) offers a novel perspective on this problem by conceptualizing the brain as a networked control system in order to explain its dynamics (Gu et al., 2015). In its most basic form, NCT considers brain dynamics as a composite outcome of a region’s connectivity profile and the necessary control inputs to guide a neural activity toward a desired state (Gu et al., 2015; Karrer et al., 2020; Parkes et al., 2024; Figure 1A). The former aspect delves into the constant anatomical interactions between various brain regions, while the latter presents an adaptable measure to optimally transition between states within the confines of energy constraints (Figure 1B). These constraints are quantified as control energy, which measures the costs entailed in controlling the brain across diverse states (Figure 1C).

Figure 1.

Network Control Theory. (A) The brain as a network of connected regions fluctuates through different modes of activity across time. We can model these transitions throughout time using the NCT. The NCT is a framework that uses information about the anatomical connections in the brain, which are fixed, and regional control inputs, which are optimized, as bases to model the temporal evolution of brain activity (highlighted by the changing node color patterns on the right). In other words, the changes in brain activity of a particular region are modeled as depending on its connectivity profile and an “activation” that spreads out along such profile, but also spreads in from other regions. Such activation is referred to control input and is fitted to represent an optimal bridge between the brain’s current and next state (Parkes et al., 2024). (B) Before employing a control model, we define what a brain state is, that is, how maps of brain activity are represented. In previous literature, states have been defined as functional intrinsic networks (Cornblath et al., 2020; Cui et al., 2020; He et al., 2022), cytoarchitectonical hierarchy levels (Parkes et al., 2022), or cognitive maps (Luppi, Singleton, et al., 2024), to name a few. These states are represented as single points in an N-dimensional state space, with N being the number of brain regions that we control. We seek to find the shortest possible trajectory between two points in this space and generate activations such that the brain is steered, or “controlled,” along a previously optimized trajectory. (C) We define the current and next state of brain activity to be related by the structural connectome of an individual (see the Methods section). The connectome imposes a resistance, or flow, into the state space such that it is easier to make activity changes between states in one direction than another. Think of walking on a hill where the walk uphill will be more strenuous than the downhill walk. In analogy to this, it costs the brain more energy, called “control energy,” to control its activity between particular states than others. This asymmetry has been the subject of study in previous literature (Luppi, Singleton, et al., 2024; Parkes et al., 2022).

Figure 1.

Network Control Theory. (A) The brain as a network of connected regions fluctuates through different modes of activity across time. We can model these transitions throughout time using the NCT. The NCT is a framework that uses information about the anatomical connections in the brain, which are fixed, and regional control inputs, which are optimized, as bases to model the temporal evolution of brain activity (highlighted by the changing node color patterns on the right). In other words, the changes in brain activity of a particular region are modeled as depending on its connectivity profile and an “activation” that spreads out along such profile, but also spreads in from other regions. Such activation is referred to control input and is fitted to represent an optimal bridge between the brain’s current and next state (Parkes et al., 2024). (B) Before employing a control model, we define what a brain state is, that is, how maps of brain activity are represented. In previous literature, states have been defined as functional intrinsic networks (Cornblath et al., 2020; Cui et al., 2020; He et al., 2022), cytoarchitectonical hierarchy levels (Parkes et al., 2022), or cognitive maps (Luppi, Singleton, et al., 2024), to name a few. These states are represented as single points in an N-dimensional state space, with N being the number of brain regions that we control. We seek to find the shortest possible trajectory between two points in this space and generate activations such that the brain is steered, or “controlled,” along a previously optimized trajectory. (C) We define the current and next state of brain activity to be related by the structural connectome of an individual (see the Methods section). The connectome imposes a resistance, or flow, into the state space such that it is easier to make activity changes between states in one direction than another. Think of walking on a hill where the walk uphill will be more strenuous than the downhill walk. In analogy to this, it costs the brain more energy, called “control energy,” to control its activity between particular states than others. This asymmetry has been the subject of study in previous literature (Luppi, Singleton, et al., 2024; Parkes et al., 2022).

Close modal

Control energy has been shown to be a powerful tool in various domains of neuroscience, providing an explanation of how the enhancement of executive functions is supported by structural changes across development (Cui et al., 2020), explaining the effects of psychedelics on brain dynamics (Singleton et al., 2022) or predicting cortical responses to electrical stimulation (Stiso et al., 2019). Nevertheless, studies employing NCT traditionally focus on the control energy between single pairs of states without accounting for how such costs organize in time. Furthermore, control energy remains a purely statistical notion (Srivastava et al., 2020), dissociated from the actual energetic currency employed by the brain (e.g., glucose and oxygen metabolism), with only recent efforts attempting to establish a connection between control costs and metabolism by showing that relative hemispheric differences in glucose uptake are mirrored by differences in control energy for temporal lobe epilepsy (TLE) patients (He et al., 2022).

Despite these advances, the costs associated with controlling a neurotypical resting-state activity as well as their energetic signature remain elusive. In light of these considerations, we propose a novel framework to examine the expenses incurred in regulating resting-state dynamics, that is, activity over time. Building upon the findings of previous reports (He et al., 2022), we further investigate whether theoretically derived control costs manifest in tangible biological metrics of energy expenditure by comparing them with normative measures of energy metabolism. Moreover, we extend our methodology to investigate organizational principles in brain dynamics through control costs. Ultimately, our efforts seek to establish a foundation for comprehending how the brain governs its functionality and the associated costs entailed in this intricate control.

Here, we present our methodology to study the control costs of human brain dynamics, that is, activity changes over time, using NCT. First, we begin by describing our analysis pipeline and the modeling choices made to simulate the control costs of resting-state dynamics. Second, we ask how our new measure of control costs compares with measurements of energy metabolism based on the positron emission tomography (PET). Finally, we conclude by showing how our framework can be leveraged to compute the transition costs between functional hierarchical levels in the brain. For all our main results, we parcellate our data into 400 cortical regions according to the delineation by Schaefer et al. (2018). Further, our analysis is primarily based on human anatomical and functional data from diffusion and functional magnetic resonance imaging (dMRI and fMRI) recordings of n = 327 participants, readily available in the Human Connectome Project Young Adult dataset (Van Essen et al., 2013). Further replications with other parcellations and datasets are presented at the end of the Results section.

Simulating the Control Costs of Resting-State Dynamics

We begin our analysis by defining brain states for our resting-state data so that we can later simulate the costs to transition between them. While the definition of a brain state is subject to debate (Kringelbach & Deco, 2020), brain states in the context of brain network control have been typically defined as task activations, for example, based on statistical maps of parameter estimates (Braun et al., 2021) or meta-analytic maps (Luppi, Singleton, et al., 2024). In resting-state studies, brain states have been defined using two different approaches. On the one hand, states have been defined according to a priori specified parcellations, where areas that fall into a given intrinsic network are represented with ones and all other areas with zeros (Cui et al., 2020; He et al., 2022). On the other, resting states have been derived from clustering techniques and then categorized based on their similarity to binary intrinsic network maps like in the previous approach (Cornblath et al., 2020; Singleton et al., 2022).

While both approaches offer a respective advantage, they also come with their drawbacks. The former approach profits from the interpretability of previously established brain networks with defined neurocognitive functions. However, such binary, discontinuous state representations are equally defined for all individuals and do not convey any temporal order of the state sequence, thereby losing unique information about the subject. The latter approach circumvents this by using data-driven clustering techniques to derive states with a temporal order from individual data. However, given the possibility of multiple solutions in clustering (Eickhoff, Thirion, Varoquaux, & Bzdok, 2015; Jain, 2010; Kleinberg, 2002), there is no guarantee that any representation of intrinsic networks are present, thereby complicating the interpretability and replicability of findings across studies. We thus propose a new approach that preserves the intelligibility of intrinsic networks as states, is nonbinary, and presents a temporal sequence tailored to the individual data.

We begin with a denoised time series activity from 400 cortical regions defined in the Schaefer atlas (Schaefer et al., 2018) that are z-scored across time. Each region is uniquely assigned to one of seven canonical intrinsic networks as proposed by Yeo-Krienen (Yeo et al., 2011). We therefore proceed by averaging regional time series data per network, resulting in seven time series, each for a Yeo-Krienen intrinsic network (Figure 2A). Intrinsic networks are then compared at each time point based on their average activity, and the network with the highest activity is designated the dominating network of that time point. Note that in all steps hereafter, as well as the rest of our analyses, we work with the 400 region representation of the data and not their network aggregate. Moreover, to control for noise, we discard time points where no network activity surpasses a fixed threshold set at half a standard deviation in an average activity. We redo our analysis with a lower threshold and also leave out the network modeling altogether to validate our approach (see the Sensitivity and Robustness Analysis section).

Figure 2.

Estimating the costs of brain dynamics. (A) We first average the activity of parcels into their respective Yeo-Krienen intrinsic networks (Yeo et al., 2011) in the Schaefer-400 parcellation (Schaefer et al., 2018), thereby aggregating the data from 400 to 7 time series (left). Then, we select the network with the highest amplitude at each time point, yielding a sequence of network dominance across time (right). (B) We proceed by averaging the original 400-parcel BOLD activity of each network across time points it dominated. (C) This yields seven representative network state maps, reflecting the average magnitude and sign of cortical activity during each network dominance. (D) Following the sequence of network dominance, we simulate the optimal control strategy in which the brain transitions between pairs of network states. Thus, each transition yields a map of momentary control costs to move between states. We then average all control maps across transitions. (E) The previous steps culminate in a map of TCE for each individual, where the value of each region represents the amount of control energy integration a region saw throughout the recording on average. VIS = visual network; SMN = somatomotor network; DAN = dorsal attention network; SAL = salience ventral network; LIM = limbic network; FPN = frontoparietal network; DMN = default mode network.

Figure 2.

Estimating the costs of brain dynamics. (A) We first average the activity of parcels into their respective Yeo-Krienen intrinsic networks (Yeo et al., 2011) in the Schaefer-400 parcellation (Schaefer et al., 2018), thereby aggregating the data from 400 to 7 time series (left). Then, we select the network with the highest amplitude at each time point, yielding a sequence of network dominance across time (right). (B) We proceed by averaging the original 400-parcel BOLD activity of each network across time points it dominated. (C) This yields seven representative network state maps, reflecting the average magnitude and sign of cortical activity during each network dominance. (D) Following the sequence of network dominance, we simulate the optimal control strategy in which the brain transitions between pairs of network states. Thus, each transition yields a map of momentary control costs to move between states. We then average all control maps across transitions. (E) The previous steps culminate in a map of TCE for each individual, where the value of each region represents the amount of control energy integration a region saw throughout the recording on average. VIS = visual network; SMN = somatomotor network; DAN = dorsal attention network; SAL = salience ventral network; LIM = limbic network; FPN = frontoparietal network; DMN = default mode network.

Close modal

With the identified time points of dominance for each intrinsic network in hand, we proceed to compute the actual state representations of networks. The dominance labels serve as temporal markers to define which time points belong together. We average the activity within such markers to yield a centroid-like representation that resembles an intrinsic network, with a continuous magnitude and sign that is captured by the activity of each individual (Figure 2B). These state coefficient maps represent each intrinsic network in the a priori defined parcellation, with additional information about regions outside the network that cofluctuate with it (Figure 2C). Essentially, we are defining network representations for each individual by tailoring the network weights to their brain activity.

To simulate the control energy between each network coefficient map, we leverage structural information about the anatomical connections of each individual, derived from dMRI. As referenced earlier, network control methods require state maps to represent initial and target positions from where we simulate the control costs to move between them. Such costs are subject to anatomical constrains imposed by the individual’s connectome. In other words, we quantify the effort required in each region to change its current activity to a desired state, with other connected regions also contributing to its change (Gu et al., 2015; Karrer et al., 2020; Parkes et al., 2024). To maintain biological plausibility, we set our simulated path to follow an optimal trajectory that balances the energy required to enable a transition with the proximity to the desired target state (Gu et al., 2015). More simply, we simulated transitions that focus on energy minimization while imposing constraints to prevent them from deviating into implausible, nonphysiological pathways that may otherwise provide a low-energy pathway.

The results of our modeling of pairwise state transitions manifest in an asymmetric matrix representing the energy cost required to transition the entire brain from one network state to another (Supporting Information Figure S2). The inherent asymmetry in these costs arises from the directional flow imposed by the structural connectivity within the state space. In another context, the states we have modeled could be thought of lying on peaks and valleys, where access is significantly easier in one direction than the other—a phenomenon analogous to the contrast between ascending and descending a hill, with the former requiring more energy.

Having captured the control cost of the transition between all pairs of resting states, we now reconcile this information with the temporal sequence of the dominant networks. This means that we use the network sequence to simulate the energy required for the transitions between states as they occurred during the recording (Figure 2D). We aggregate the transition energy over all pairs of states in the recording sequence and then normalize the cumulative energy by dividing it by the number of transitions, that is, we average the overall energy across transitions. This normalization ensures that our measurement remains comparable across different recording durations. Ultimately, our approach culminates in the derivation of a metric we refer to as time-averaged control energy, or TCE (Figure 2E).

TCE is an estimate of the temporal control costs associated with idiosyncratic intrinsic network activations across time (Figure 3A). It is a spatially heterogeneous measure, but we find that its distribution in the left and right hemispheres is statistically indistinguishable (Mann–Whitney U = 76,390, p = 0.27; Figure 3B). Moreover, previous evidence indicates that control energy is modulated by age (Cui et al., 2020; Tang et al., 2017). We find that stratifying by age (using bins of 22–25, 26–30, 31–35, and 36+, with 247, 527, 418, and 14 individuals, respectively) does not lead to any statistically significant group differences (Levene W = 108, p = 0.955; ANOVA (analysis of variance), F(3, 1596) = 0.133, p = 0.94; Figure 3C). We thus show that the control costs of human brain dynamics are region specific, consistent between hemispheres, and conserved across ages studied here.

Figure 3.

TCE. (A) Averaging all TCE maps across individuals, we observe that the costs to control resting-state dynamics are regionally heterogeneous. (B) Regions in the dorsal attention network require the most control with limbic regions costing the least. In general, temporal control costs display no preference for a hemisphere, with comparable values evident in both the left and right hemispheres (Mann–Whitney U = 76,390, p = 0.27). (C) When stratifying our analysis by age, we see no statistical difference in temporal control costs across ages groups (ANOVA F[3, 1596] = 0.133, p = 0.94). Each dot represents the subject-averaged value for a region of interest. Bars around the mean indicate the standard deviation of the data.

Figure 3.

TCE. (A) Averaging all TCE maps across individuals, we observe that the costs to control resting-state dynamics are regionally heterogeneous. (B) Regions in the dorsal attention network require the most control with limbic regions costing the least. In general, temporal control costs display no preference for a hemisphere, with comparable values evident in both the left and right hemispheres (Mann–Whitney U = 76,390, p = 0.27). (C) When stratifying our analysis by age, we see no statistical difference in temporal control costs across ages groups (ANOVA F[3, 1596] = 0.133, p = 0.94). Each dot represents the subject-averaged value for a region of interest. Bars around the mean indicate the standard deviation of the data.

Close modal

Recapitulating, our framework first begins by identifying the dominant intrinsic networks at each time point (Figure 2A). This information is used to estimate the average activity across time points where each network dominates (Figure 2B), resulting in state maps for each intrinsic network (Figure 2C). We then simulate the control energy to transition between each state in order to map these energetic costs to the sequence of dominating networks from our first step. This allows us to track the amount of control energy required to transition between time points (Figure 2D). Finally, we average across all transitions to determine the TCE of an individual, which is an estimate of the average costs to control their resting-state activity over time (Figure 2E). Our measure offers a new insight into the energy required to control resting-state activity, which we proceed to validate with actual measures of energetic costs, such as glucose and oxygen metabolism.

Grounding Control Costs in Metabolism

Previously, NCT metrics were shown to correlate with physiological markers of energy metabolism (He et al., 2022), thereby establishing an initial bridge between theoretical and tangible indicators of energy. As such, we contextualize our results with neurobiological maps to look for energetic correlates of temporal control costs. To accomplish this, we scrutinize the relationship between temporal control costs and two pivotal substrates of metabolic energy: oxygen and glucose metabolism (Camandola & Mattson, 2017; Raichle & Gusnard, 2002; Figure 4).

Figure 4.

Relationship to metabolism. (A) Left: We compared how TCE spatially covaries with glucose metabolism from an internally acquired PET dataset (Castrillon et al., 2023). Glucose uptake has been previously shown to scale with momentary control energy (He et al., 2022). We, however, do not find a similar significant relationship to control energy when considering its temporal dimension (Spearman ρ = 0.144, PSMASH = 0.347). Shades around the regression line represent 95% confidence interval with bootstrap samples. Right: Distribution of spatial nulls. Correlation results accounted for spatial autocorrelation in the data by using variogram-matched nulls (Burt et al., 2020; Markello & Misic, 2021; Váša & Mišić, 2022). (B) Autocorrelation in the data by using variogram-matched nulls (Burt et al., 2020; Markello & Misic, 2021; Váša & Mišić, 2022).

Figure 4.

Relationship to metabolism. (A) Left: We compared how TCE spatially covaries with glucose metabolism from an internally acquired PET dataset (Castrillon et al., 2023). Glucose uptake has been previously shown to scale with momentary control energy (He et al., 2022). We, however, do not find a similar significant relationship to control energy when considering its temporal dimension (Spearman ρ = 0.144, PSMASH = 0.347). Shades around the regression line represent 95% confidence interval with bootstrap samples. Right: Distribution of spatial nulls. Correlation results accounted for spatial autocorrelation in the data by using variogram-matched nulls (Burt et al., 2020; Markello & Misic, 2021; Váša & Mišić, 2022). (B) Autocorrelation in the data by using variogram-matched nulls (Burt et al., 2020; Markello & Misic, 2021; Váša & Mišić, 2022).

Close modal

Our initial investigation involved comparing TCE with the regional variations in cerebral glucose metabolism (CMRglc) across the cortex. To achieve this, we employed a group-average template derived from an internally acquired dataset (Castrillon et al., 2023) comprising n = 20 subjects who underwent PET scans, utilizing 18F-fluorodeoxyglucose (FDG) as a tracer to monitor their glucose uptake (see the Methods section).

Given the inherent spatial smoothness of our PET map, we employ rigorous measures to validate our results by contrasting them against 10,000 null maps that replicate the variogram of the data (Burt et al., 2020; see the Methods section). These null maps randomize the values in our PET map while accounting for the spatial autocorrelation inherent in the recordings, enabling us to distinguish true effects from spatially induced artifacts using a modified p value (PSMASH; Burt et al., 2020; Markello & Misic, 2021). Our results showed a weak correlation between TCE and CMRglc, which, however, did not exhibit statistical significance when the spatial smoothness of the dataset was taken into account (Spearman ρ = 0.144, p = 0.003, PSMASH = 0.347).

Subsequently, we investigate the relationship between TCE and cerebral oxygen uptake within the cortex. Here, we employed a normative map of cerebral oxygen metabolism (CMRO2) derived from Vaishnavi et al. (2010). This map was computed from n = 33 participants who underwent PET scans to estimate their cerebral oxygen consumption, utilizing 15O-labeled oxygen (see the Methods section). We note that the map is not quantified into cerebral metabolic units; however, we retain the name for oxygen uptake as CMRO2, in accordance with the original authors’ nomenclature for their activity concentration.

Given the same technical limitations as with the CMRglc map, we generate 10,000 null maps that preserve the spatial autocorrelation of the original data using the same variogram-matching algorithm (Burt et al., 2020). Our results show that TCE significantly correlates with CMRO2 and that this relationship is significant even after accounting for the spatial autocorrelation in the data (Spearman ρ = 0.35, p < 0.001, PSMASH = 0.027). Altogether, we find that the regional cortical pattern of temporal control costs resembles the spatial pattern of oxygen uptake, establishing further evidence about the link between biological and control energy.

Temporal Organization of Control Costs

We now turn to a new perspective on brain states by reframing our analysis from the perspective of functional hierarchical processing as defined by Mesulam (1998). Previous work has shown that control energy differs across cortical hierarchies (Parkes et al., 2022). In this context, we examine the control costs of switching within and between hierarchical levels across time, ultimately looking at their dynamical organization in time.

We begin by mapping the sequence of dominant networks to their respective levels within the functional hierarchy of the brain (Figure 5A). Visual and sensorimotor networks are hereafter referred to as “unimodal,” a distinction that arises from their specialized processing that assumes specific inputs (Margulies et al., 2016; Mesulam, 1998). In turn, we refer to all other networks as “heteromodal” (associative) as they integrate multiple functional streams (Mesulam, 1998). Bimodal mapping ensures consistency in our subsequent validation with a different parcellation that defines intrinsic networks differently (Supporting Information Figure S5). Moreover, we show that adopting this new perspective provides a fresh view of the dynamics of the brain, uncovering a relationship between TCE and state-switching behavior. We follow the same framework to compute TCE, but now separately average transitions within unimodal networks, that is, from unimodal to unimodal, and within heteromodal networks, as well as transitions between unimodal and heteromodal levels.

Figure 5.

Temporal organization of control costs. (A) We map our sequence of dominating networks to their respective hierarchical level. This means sensory networks are jointly labeled as unimodal networks (U), with the remaining networks defined as heteromodal, association networks (H). (B) Transitions within hierarchical levels represent the majority of transitions (86.5%) at rest. (C) Averaging control costs across the whole brain, we observed that transitions within a hierarchical level required significantly less TCE than between. Asterisks indicate high significance (Mann–Whitney U = 376,072, p < 0.001). Each dot per transition type represents a subject. (D) We relate whole-brain control costs to transition frequency using repeated measures correlation (Bakdash & Marusich, 2017) to assess common intra-individual relationships and find that the number of transitions per transition type is negatively correlated to the cost needed to perform such switch (rrm = −0.619, p < 0.001).

Figure 5.

Temporal organization of control costs. (A) We map our sequence of dominating networks to their respective hierarchical level. This means sensory networks are jointly labeled as unimodal networks (U), with the remaining networks defined as heteromodal, association networks (H). (B) Transitions within hierarchical levels represent the majority of transitions (86.5%) at rest. (C) Averaging control costs across the whole brain, we observed that transitions within a hierarchical level required significantly less TCE than between. Asterisks indicate high significance (Mann–Whitney U = 376,072, p < 0.001). Each dot per transition type represents a subject. (D) We relate whole-brain control costs to transition frequency using repeated measures correlation (Bakdash & Marusich, 2017) to assess common intra-individual relationships and find that the number of transitions per transition type is negatively correlated to the cost needed to perform such switch (rrm = −0.619, p < 0.001).

Close modal

After mapping the sequence of dominating networks, we observe that, on average for individuals, transitions occur predominantly within a hierarchy rather than between hierarchies (Figure 5B). Analogously, the control costs to transition between hierarchies are more demanding than those within a hierarchical level (Mann–Whitney U = 376,072, p < 0.001; Figure 5C). Using repeated-measures correlation (Bakdash & Marusich, 2017; see the Methods section) to account for repeated measures of control costs for each participant as well as their individual traits, we examine the commonalities in the relationship between participants’ control costs and transition frequency across hierarchies. We find a robust negative correlation throughout the cohort (rrm = −0.619, p < 0.001), indicating that more frequent transitions within a hierarchy are associated with lower average control energy costs, in stark contrast to infrequent, costly transitions that bridge different hierarchical domains. We show that this relationship is not only an artifact of the network partition in our parcellation but is also evident when using alternative partitions (see the Sensitivity and Robustness Analysis section). Altogether, our findings suggests a dynamic temporal principle within the brain that strives to minimize its control costs by keeping expensive transitions sparse.

Sensitivity and Robustness Analysis

To ensure the robustness of our results, we replicate our analyses on an independent dataset (Castrillon et al., 2023). This additional dataset has lower spatial and temporal resolution and is subjected to different preprocessing protocols (see the Methods section). Indeed, even in the face of these differences, our results remain unchanged (Supporting Information Figure S6).

We then repeat our analysis using an alternative cortical parcellation (Gordon et al., 2016). This approach seeks to validate the reproducibility of our results with a different definition of intrinsic networks. We find that this approach leads to consistent results (Supporting Information Figure S5). Additionally, we investigate the impact of a coarser partition by dividing the cortex into 200 regions instead of the original 400 regions. We find that this modification does not affect our results (Supporting Information Figure S4).

Further, we test the impact of our state modeling step by leaving it out and instead simulate the control costs to transition between time points of raw activity, that is, our state maps are represented as the BOLD activity at each time point. We form a group-average template like in our main results and compare these results with the temporal signal-to-noise ratio (tSNR) of each region. tSNR is a metric to quantify the quality of a BOLD signal, calculated as the mean of the signal prior to demeaning and divided by the signal standard deviation (Welvaert & Rosseel, 2013). We compute the tSNR for each participant run in the Human Connectome Project (HCP) dataset, rescale it from 0 to 1 using min-max normalization, and finally average it across runs.

We find that, when leaving out the state modeling step, regions with high TCE are also regions that exhibit a low tSNR (Spearman ρ = 0.36, Pspin = 0.042). In contrast, using our approach, we observe that the group-level results are unlikely related to the group-level regional tSNR (Spearman ρ = 0.234, PSMASH = 0.232). Following our main approach, we discard time points with an average activity below 0.5 but observe that, using a lower threshold of 0 instead, our group-level results again show no significant relationship to the regional signal quality (Spearman ρ = 0.274, PSMASH = 0.099). Moreover, we observed on an individual level that a larger extent of participants’ TCE maps were significantly related to their tSNR when leaving out our modeling approach (Supporting Information Figure S1). Taken together, this suggests that modeling states using our approach mitigates the influence of noise in our analysis.

We conclude by examining whether temporal control costs are related to the occurrence of resting states. This would mean that the patterning of individual TCE maps are a reflection of the most frequently occurring states. To test this, we create a representative map of the occurring intrinsic networks by counting how frequent, that is, how many time points, each state occurs. Dividing this number by the total number of time points yields an estimate of how much each network was present in the recording, for example, the visual network dominated 20% of the time. We multiply this ratio with its corresponding state coefficient map, thereby creating time-weighted maps, which we can sum together to create a map that represents the frequency of each network and its spatial distribution. Interestingly, we find no significant association between this map and our TCE map (Spearman ρ = 0.27, PSMASH = 0.092; Supporting Information Figure S3). Altogether, we conclude that our findings result from a convergence of characteristic sequences of resting states and individual anatomical constraints, which together form the basis to simulate the control costs of resting-state dynamics.

In the present study, we introduce a metric for assessing the control costs inherent in human brain dynamics, referred to as TCE. TCE serves as a composite measure quantifying the average control exerted by an individual to sustain a dynamic neural activity over time. Notably, our findings reveal a statistically significant spatial resemblance between the distribution of TCE and cortical oxygen metabolism, a major energetic substrate crucial for compensating a neural activity (Harris, Jolivet, & Attwell, 2012). Our results further show that transitions within uni- or heteromodal networks entail distinct temporal control costs as opposed to transitions between modalities, thereby offering insight into the temporal organization of brain dynamics. Importantly, we provide an open-access code to facilitate the replication and extension of our results in subsequent studies.

Temporal Control Is a Confluence of Anatomical Constraints and State Sequences

NCT is a powerful tool to investigate the emergence of dynamics from the structural scaffold of the brain (Parkes et al., 2024). Its application in neuroscience has yielded valuable insights, for example, into individual differences in psychiatric conditions (Luppi, Singleton, et al., 2024; Parkes et al., 2021) or developmental changes in youth (Cui et al., 2020; Tang et al., 2017). Notably, controllability is reported to increase asymptotically and plateau with adolescence (Tang et al., 2017), potentially explaining the absence of age effects observed in our results.

We find that the distribution of temporal control costs is spatially distinct, with certain brain regions requiring higher TCE than others (Figure 3). The extent of engagement in this control dynamic may be influenced by the manner in which control inputs propagate, reminiscent of the role of connectivity profiles in facilitating diverse spreading dynamics across brain regions (Mišić et al., 2015). This spreading process may parallel the flow of control throughout the connectome (Srivastava et al., 2020).

Control costs are further influenced by the specific states the brain traverses. Our results corroborate previous findings (Luppi, Singleton, et al., 2024; Parkes et al., 2022), indicating that control energy is dependent on both initial and goal states (Figure 5; Supporting Information Figure S2). This underscores the significance of our state-modeling approach that determines the sequence of network dominance based on intrinsic network activations, thereby establishing the order of state transitions.

Our methodology demonstrates robustness across different partitions of intrinsic networks (Supporting Information Figures S4 and S5), as well as separate datasets (Supporting Information Figure S6). Nonetheless, our state-modeling approach trades off robustness for temporal granularity, thereby losing the ability to capture temporal shifts in brain activity. Indeed, brain dynamics are hypothesized to be temporally organized such that high-amplitude modes of activity are counterbalanced by transient “off” states (Saggar, Shine, Liégeois, Dosenbach, & Fair, 2022; Vidaurre, Smith, & Woolrich, 2017; Zamani Esfahlani et al., 2020). We prompt future investigations to experiment with alternative state-modeling methods, like sliding-window approaches, for improved temporal precision. Appropriate tests should, however, be developed in tandem to distill real fluctuations from spurious ones (Leonardi & Van De Ville, 2015; Zalesky & Breakspear, 2015).

Finally, we note that our network control model is linear, therefore assuming that linear dynamics govern brain transitions (Karrer et al., 2020; Parkes et al., 2024). Previous studies demonstrated how local nonlinear fluctuations in the cortex support the emergence of intrinsic resting-state networks (Cabral, Hugues, Sporns, & Deco, 2011; Mantini, Perrucci, Del Gratta, Romani, & Corbetta, 2007) and, as such, denoted the utility of nonlinear models to simulate neural dynamics (Breakspear, 2017). However, recent evidence has shown that models based on first-order linear approximations can approximate macroscale activity (Nozari et al., 2023). Indeed, previous work showed that linear control model predictions align with nonlinear biophysical model simulations for personalized stimulation protocols (Muldoon et al., 2016). Moreover, compared with nonlinear control models, a desirable feature of linear control models is that they are more interpretable, which can facilitate incorporation of additional biological properties, such as receptor availability (Singleton et al., 2022, 2023) or cortical abnormalities in patient populations (Luppi, Singleton, et al., 2024). Consequently, we see the potential for future efforts to expand on our work in order to investigate how temporal control costs are modulated by incorporating additional biological insights or clinical distinctions.

Control Costs Reflect Metabolic Features of the Brain

Our results show that temporal control costs colocalize with oxygen metabolism, but not glucose metabolism (Figure 4). A previous study from He et al. (2022) showed that control energy and glucose uptake are similarly put out of balance between hemispheres in TLE patients. The authors, however, only show this relationship in the limbic system, and, moreover, their methodology does not consider the costs of temporal switching but seeks to link relative hemispheric asymmetries in momentary control costs and glucose uptake to explain structural changes in TLE patients. In contrast, our study relates control and metabolism on an absolute whole-brain basis and considers the temporal dimension of control costs.

Moreover, while both glucose and oxygen constitute fundamental components for meeting the energetic demands of neural activity (Camandola & Mattson, 2017; Harris et al., 2012; Li & Sheng, 2022; Mergenthaler, Lindauer, Dienel, & Meisel, 2013; Vaishnavi et al., 2010; Watts, Pocock, & Claudianos, 2018), they exhibit distinct spatial patterns. The predominant energy-producing pathways in the brain encompass glycolysis and oxidative metabolism, where glycolysis relies on glucose and oxidative metabolism relies on glucose-derived pyruvate and oxygen (Camandola & Mattson, 2017; Dienel, 2012). The disparate spatial distribution of glucose and oxygen metabolism is attributed to their dissimilar upscaling capacities for energy generation (Díaz-García et al., 2017; Dienel, 2019; Zala et al., 2013). Specifically, certain regions are posited to favor glycolysis beyond basal levels of oxygen and glucose consumption at an equable ratio in order to afford greater flexibility in energy demand necessary for supporting human cognition (Luppi, Rosas, et al., 2024; Vaishnavi et al., 2010). In light of these considerations, our association of temporal control costs with the spatial distribution of oxygen, rather than glucose metabolism, supports the interpretation that the control of human brain dynamics is grounded in the baseline metabolic activity constituent of oxidative energy generation.

We note, however, that our conclusion is based on group-level comparisons of control costs and metabolic maps. Further research is warranted to investigate the nuanced relationship between control costs and energy metabolism at an individual level. Advances in MRI now enable the simultaneous acquisition of BOLD and oxygen metabolism data from the same subjects (Epp et al., 2024), providing an avenue for an in-depth exploration of this relationship.

State Switching Is Grounded by Control Costs

We observed that the costs required to transition within and between hierarchical levels, that is, unimodal or heteromodal, were significantly different (Figures 5C and 5D). Similarly, this phenomenon was mirrored in the frequency with which the brain either transitions out or persists in a hierarchical level (Figures 5B and 5D). This means that costly transitions between hierarchical levels were rare, with more efficient transitions within hierarchical levels happening more frequently.

In order to enable information flow, the human brain is posited to switch between modes of integration and segregation over time, where segregation refers to regions clustering together to form functionally distinct modules and integration represents a form of global intercommunication between modules (Deco, Tononi, Boly, & Kringelbach, 2015; Finc et al., 2020; Shine & Poldrack, 2018). Building upon this, we propose that the transitions between functional hierarchical levels may reflect an integration-segregation dynamic. In this sense, switches between hierarchical levels might represent modes of communication among functional hierarchies, counterbalanced by switches within hierarchical levels as modes of segregated activity. As previous studies have shown the disruption of integration-segregation dynamics by pathology or altered states of consciousness (Lord, Stevner, Deco, & Kringelbach, 2017; Luppi et al., 2019), future works could investigate whether such conditions equally disrupt state-switching dynamics.

Furthermore, we hypothesize that the inverse correlation between control costs and state visits represents a mechanism where the brain seeks to maintain functional diversity, while minimizing its energy expenditure. This is similar to the hypothesis put forward by Zalesky, Fornito, Cocchi, Gollo, and Breakspear (2014), where the authors posit that the brain transitions through costly, intermittent states in order to facilitate global integration while keeping energy demand minimal. As such, we encourage future efforts to further elucidate our understanding of the intricacies governing state-switching dynamics and control costs.

In conclusion, our work presents a new perspective of control energy as a measure of the costs to control human brain dynamics. We are able to show that temporal control costs are spatially related to oxygen metabolism as an energy substrate and uncover an organizational principle that relates temporal dynamics to their costs in the brain. As such, we envision that our work can help advance research on the costs of regulating brain activity.

All preprocessed data are available at https://osf.io/nw9zt. The code and additional data used to perform the analyses are available at https://github.com/NeuroenergeticsLab/control_costs.

Data

HCP.

The main data used for this study consisted of resting-state time series from fMRI and structural connectomes from dMRI taken from the Human Connectome Project S900 Young Adult release (Van Essen et al., 2013). Scans from 327 unrelated participants (mean age = 28.6 ± 3.73 years, 55% females) were used to ensure that familial factors do not confound our analysis (Miranda-Dominguez et al., 2018). Informed consent was obtained for all subjects (the protocol was approved by the Washington University Institutional Review Board as part of the HCP). The participants were scanned in the HCP’s custom Siemens 3T “Connectome Skyra” scanner, and the acquisition protocol included four 15-min resting-state fMRI sessions and a high angular resolution diffusion imaging (HARDI) sequence. The resting-state fMRI data were acquired using a gradient-echo EPI sequence (TR = 720 ms; TE = 33.1 ms; FOV = 208 × 180 mm2; voxel size = 2 mm3; number of slices = 72; and number of volumes = 1,200). The dMRI data were acquired with a spin-echo EPI sequence (TR = 5520 ms; TE = 89.5 ms; FOV = 210 × 180 mm2; voxel size = 1.25 mm3; b value = three different shells, i.e., 1000, 2,000, and 3,000 s/mm2; number of diffusion directions = 270; and number of b0 images = 18). Additional information regarding the acquisition protocol can be found under Van Essen et al. (2013).

To process the functional data, each run of each subject’s resting-state fMRI recording was preprocessed in terms of gradient distortion correction, motion correction, and spatial normalization according to S. M. Smith et al. (2013) and Glasser et al. (2013). Artifacts were then removed using ICA-FIX (Salimi-Khorshidi et al., 2014). Intersubject registration of the cerebral cortex was carried out using areal feature-based alignment and the multimodal surface matching algorithm (Robinson et al., 2014). The clean time series were then parcellated into 400 and 200 cortical region time series according to the Schaefer functional atlas (Schaefer et al., 2018).

To process the diffusion data, structural connectomes were reconstructed from the dMRI data using the MRtrix3 package (Tournier et al., 2019). Gray matter was parcellated into 400 and 200 cortical regions according to the Schaefer functional atlas (Schaefer et al., 2018), and fiber orientation distributions were generated using a multishell, multitissue constrained spherical deconvolution algorithm (Dhollander, Raffelt, & Connelly, 2016; Jeurissen, Tournier, Dhollander, Connelly, & Sijbers, 2014). The initial tractogram was generated with 40 million streamlines, with a maximum tract length of 250 and a fractional anisotropy cutoff of 0.06. Spherical-deconvolution informed filtering of tractograms (SIFT2) was used to reconstruct whole-brain streamlines weighted by cross-section multipliers (R. E. Smith, Tournier, Calamante, & Connelly, 2015). More information regarding the individual network reconstructions is available in Park et al. (2021).

TU Munich dataset.

The replication data consisted of simultaneously measured FDG-PET and fMRI, with subsequent dMRI recordings, all while the participants kept their eyes open. Scans from 20 participants (mean age = 34.2 ± 5.99 years, 50% females) were used. Informed consent was obtained for all subjects (the protocol was approved by the Technical University of Munich Review Board). The participants were scanned in an integrated PET/MR (3T) Siemens Biograph mMR scanner (Siemens, Erlangen, Germany) and used a 12-channel phase-array head coil for the MRI acquisition. The PET data were collected in a list-mode format with an average intravenous bolus injection of 184 MBq (SD = 12 MBq) of [18F]FDG. In parallel to the PET measurement, automatic arterial blood samples were taken from the radial artery every second to measure blood radioactivity using a Twilite blood sampler (Swisstrace, Zurich, Switzerland).

The fMRI data were acquired during a 10-min time interval using a single-shot EPI sequence (300 volumes; 35 slices; repetition time, TR = 2,000 ms; echo time, TE = 30 ms; flip angle [FA] = 90°; FOV = 192 × 192 mm2; matrix size = 64 × 64; voxel size = 3 × 3 × 3.6 mm3). Diffusion-weighted images were acquired using a single-shot EPI sequence (60 slices; 30 noncolinear gradient directions; b value = 800 s/mm2 and one b = 0 s/mm2 image; TR = 10,800 ms, TE = 82 ms; FA = 90°; FOV = 260 × 264 mm2; matrix size = 130 × 132; voxel size = 2 × 2 × 2 mm3). Additional information regarding the acquisition protocol can be found under Castrillon et al. (2023).

To process the functional data, each subject’s resting-state fMRI recording was preprocessed using the Configurable Pipeline for the Analysis of Connectomes (CPAC v1.4.0; Craddock et al., 2013). This included slice-timing correction, motion correction, spatial normalization, quadratic and linear detrending of scanner drift, anatomical CompCor regression of white matter and cerebrospinal fluid activity (Behzadi, Restom, Liau, & Liu, 2007), and subsequent bandpass filtering (0.01–0.1 Hz). More details can be found in Castrillon et al. (2023).

To process the diffusion data, structural connectomes were generated using the MRtrix3_connectome BIDS App (Gorgolewski et al., 2017), which operates principally using tools provided in the MRtrix3 package (Tournier et al., 2019). This included denoising (Veraart et al., 2016), Gibbs ringing removal (Kellner, Dhital, Kiselev, & Reisert, 2016), preprocessing (Andersson, Skare, & Ashburner, 2003; Andersson & Sotiropoulos, 2016), bias field correction (Tustison et al., 2010), intermodal registration (Bhushan et al., 2015), brain extraction (S. M. Smith, 2002), T1 tissue segmentation (Patenaude, Smith, Kennedy, & Jenkinson, 2011; R. E. Smith, Tournier, Calamante, & Connelly, 2012; S. M. Smith, 2002; Zhang, Brady, & Smith, 2001), spherical deconvolution (Jeurissen et al., 2014; Tournier, Calamante, Gadian, & Connelly, 2004), and probabilistic tractography (Tournier, Calamante, & Connelly, 2011) utilizing Anatomically Constrained Tractography (R. E. Smith et al., 2012) and dynamic seeding (R. E. Smith et al., 2015). The resulting fiber track files were subsequently converted into streamline counts by counting the number of streamlines that passed through one of the 400 cortical regions according to the Schaefer functional atlas (Schaefer et al., 2018). In order to compensate for the bias toward longer fibers inherent in the tractography procedure, as well as differences in region size, we normalized the streamline count by the average length of the streamlines and average surface area of the two regions (Hagmann et al., 2008).

To process the PET data, the first 45 min of the PET acquisition were reconstructed offline using the NiftyPET library (Markiewicz et al., 2018) based on the ordered subsets expectation maximization (OSEM) algorithm with 14 subsets, four iterations, and divided into 33 dynamic frames: 10 × 12 s, 8 × 30 s, 8 × 60 s, 2 × 180 s, and 5 × 300 s. The attenuation correction was based on the T1-derived pseudo-CT images (Burgos et al., 2014). All reconstructed PET images were motion corrected and spatially smoothed (Gaussian filter, full width at half maximum [FWHM] = 6 mm). The net uptake rate constant (Ki) was calculated using the Patlak plot model (Patlak & Blasberg, 1985) based on the last five frames of the preprocessed PET images (frames between 20 and 45 min) and the arterial input function derived from the preprocessed arterial blood samples. The cerebral metabolic rate of glucose (CMRglc) was calculated by multiplying the Ki map with the concentration of glucose in the plasma of every participant, divided by a lumped constant of 0.65 (Wu et al., 2003). Finally, the CMRglc maps were partial volume corrected using the grey matter, white matter, and cerebrospinal fluid masks derived from the T1 images using the iterative Yang method (Yang et al., 2017) and finally registered to the MNI152NLin6ASym 3 mm template through the anatomical image.

Additional data.

In addition to the metabolic map of glucose uptake from the TUM dataset, we complemented our analysis with a normative, group-average map of oxygen metabolism (CMRO2) from Vaishnavi et al. (2010), acquired through neuromaps (Markello et al., 2022). The brain map was derived from 33 healthy, right-handed, neurologically normal participants (mean age = 25.4 ± 2.6 years, 58% females) that were recruited from the Washington University community. The recording was performed using a Siemens model 961 ECAT EXACT HR 47 PET scanner (Siemens/CTI) with 47 slices encompassing an axial field of view of 15 cm. Transverse resolution was 3.8–5.0 mm FWHM, and axial resolution was 4.7- to 5.7-mm FWHM. Attenuation data were obtained using 6868GeGa rotating rod sources to enable quantitative reconstruction of subsequent emission scans. Emission data were obtained in the 2D mode (interslice septa extended). The PET data were reconstructed using a ramp filter (6-mmm FWHM) and then blurred to a 12-mm FWHM. Distribution of CMRO2 was measured with a 40-s emission scan (derived from a 120-s dynamic scan) after brief inhalation of 60 mCi of [15O]oxygen in room air. More details on the recording parameters and processing can be found in Vaishnavi et al. (2010).

Network Control Theory (NCT)

Network Control Theory (NCT) is a framework that provides the approach to computing metrics related to the control of brain activity.

Fundamental to NCT is the definition of the brain as a networked system with nodes connected by edges and defined by an adjacency matrix A. Biologically, the matrix A represents the thick bundles of myelinated axonal fibers that run through large-scale connections of brain regions (Sporns, 2011). These fibers are thought to play a critical role in coupling the activity of distant brain regions (Avena-Koenigsberger et al., 2017).

By accounting for connections of the brain-region level, an equation relating the brain structure to time-evolving brain activity x(t) can be formulated as:
(1)
where x˙ and x are of size N × 1 and AN × N, with N is the number of brain regions. In practice, the matrix A is normalized in order to avoid infinite growth, and as such defined as:
(2)
where λ(A)max is the largest eigenvalue of A and c a constant that determines if the system approaches the corresponding mode of A (c > 0) or decays toward 0 (c > 0). Here, c is set to 1. For simplicity, Anorm will be referred to as A from now on, assuming that all the following equations make use of a normalized matrix.

Looking at the assumptions made, Equation 1 imposes that the temporal evolution of the brain is a linear function that is described by its connectivity and its current state in time. In addition to linearity, it is important to note that this equation assumes that A is not changing, that is, the brain is a time-invariant system and that brain dynamics are noise free (Karrer et al., 2020; Srivastava et al., 2020).

From here, Equation 1 can be extended to account for controlled dynamics that steer the system away from its natural trajectories through an external input. Formally,
(3)
with u(t) representing the injected control input into the system at time t and the columns in B representing the total number of voxels or regions of interest to control. Depending on the number of regions m that the modeler decides to control, u(t) is of size m × 1 and, correspondingly, B is of size N × m. Each element in bij represents the influence that an external input uj(t) has on region i, where j = 1, …, m and i = 1, …, N. By including controlled dynamics in Equation 3, it is now possible to model optimal trajectories between brain states and calculate their optimal control energy.

Optimal control energy.

Optimal control energy provides a measure of the controllability of a system in terms of (a) the state trajectory traversed and (b) the work required to reach a state. It further specifies a time horizon T during which the control input u(t) is effective in moving the system from the initial state x(0) = x0 to the goal state x(T) = xT.

Formally, the problem of reaching xT starting from x0 while keeping (a) and (b) minimal is as follows:
(4)
with Δx(t) being the distance to the goal state at time t, that is,
(5)
Moreover, ρ determines the relative weighting between the costs associated with the length of the state trajectory and input energy. A usual choice is to set ρ = 1 as to weight both goals equally (Gu et al., 2017; Karrer et al., 2020). Once Equation 4 is minimized, an optimal control energy can be derived from u(t)*:
(6)

The optimal energy E* was computed for each spatial unit—in this case, regions from a parcellation. These values therefore represent the quadratic input across model time to optimally move toward a goal state given an initial state.

Readers interested in a more in-depth discussion about the mechanisms of NCT are invited to read the review from Srivastava et al. (2020) and Parkes et al. (2024). For more discussions on parameter sensibilities of the model, interested readers are referred to the exhaustive study by Karrer et al. (2020).

Time-averaged control energy.

We begin by averaging at each time point the BOLD activity of all parcels belonging to one of seven intrinsic networks in the Schaefer parcellation (Schaefer et al., 2018). The network with the highest average activity is consequently called the dominant network for that time point, which results in a sequence of network dominance across time. For each network, we proceed to average all the time points that were assigned to it, thereby estimating what the brain looks like when it is dominated by an intrinsic network. These cortical maps represent states in our control framework that are continuous, rather than binary (i.e., ones in regions assigned to an intrinsic network, and zeros everywhere else). We then set the control input in our model to be uniform across the cortex, meaning that B is set to be the identity matrix. Finally, using Equation 6, we simulate the control energy to transition between all possible pairs of dominating states, resulting in 49 different control energy maps. These maps, together with the sequence of network dominance, allows us to determine the TCE for each individual:
(7)

Here, E01* is defined as the optimal control energy to transition from time point 0 to 1, and P is the total number of transitions, which is equivalent to the total number of time points in a recording minus one, that is P = Trecording − 1. The resulting TCE value represents the average control energy exerted onto a region across a resting-state recording. As we show in our results, each region sees a different demand for control energy, resulting in a spatially heterogeneous map of control costs for each individual. Finally, we average each map of TCE across individuals to have a group representation of TCE, which we can statistically compare with other brain maps.

Spatial Autocorrelation Preserving Null Model

For all of our hypothesis tests that relied on comparing two brain maps, we employed BrainSMASH (Burt et al., 2020) as implemented in the neuromaps (Markello et al., 2022) library. Briefly, the methodology of BrainSMASH entails two steps: (a) random permutation of values within a designated brain map and (b) smoothing and rescaling to restore the spatial autocorrelation characteristic of the original data.

In the first step, values within the brain map undergo random permutation. Subsequently, the permuted data undergo a transformation to reintroduce spatial autocorrelation. This transformation is expressed as
(8)
where x′ represents the permuted data, zN(0, 1), and α and β are estimated via least-squares optimization, aligning variograms of the original and permuted data.

Repeated-Measures Correlation

Repeated-measures correlation (rrm) is a statistical technique developed by Bakdash and Marusich (2017) to examine the common association between variables within each subject, focusing on accounting for within-individual associations in paired measures. This method is particularly valuable when dealing with repeated measures on the same participants—in our case four control cost estimates for each hierarchical transition. Because this method only assesses the common slope across participants, we can investigate whether there is a consistent relationship trend throughout measurements while automatically correcting for individual confounding factors such as the participant’s gender or age.

We thank Filip Milisav, Justine Hansen, Vincent Bazinet, Asa Farahani, Zhen-Qi Liu, Moohebat Pourmajidian, Laura Suarez, and Markus Ploner for their helpful discussion. E.G.C. acknowledges support from the Molson Foundation, the IFI fellowship of the German Academic Exchange Service, and the Elite Network of Bavaria. A.I.L. was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC; funding reference number 202209BPF-489453-401636; Banting Postdoctoral Fellowship). V.R. has been funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (ERC Starting Grant, ID 759659).

Supporting information for this article is available at https://doi.org/10.1162/netn_a_00425.

Eric G. Ceballos: Conceptualization; Data curation; Formal analysis; Funding acquisition; Investigation; Methodology; Project administration; Resources; Software; Supervision; Validation; Visualization; Writing – original draft; Writing – review & editing. Andrea I. Luppi: Conceptualization; Investigation; Methodology; Supervision; Validation; Writing – original draft; Writing – review & editing. Gabriel Castrillon: Conceptualization; Data curation; Software; Validation. Manish Saggar: Conceptualization; Formal analysis; Funding acquisition; Investigation; Methodology; Project administration; Supervision; Validation; Writing – review & editing. Bratislav Misic: Conceptualization; Data curation; Funding acquisition; Investigation; Methodology; Project administration; Supervision; Validation; Visualization; Writing – review & editing. Valentin Riedl: Conceptualization; Data curation; Funding acquisition; Investigation; Methodology; Project administration; Supervision; Validation; Visualization; Writing – original draft; Writing – review & editing.

Valentin Riedl, HORIZON EUROPE European Research Council (https://dx.doi.org/10.13039/100019180), Award ID: 759659.

Network control theory:

Engineering framework used in neuroscience to depict the brain as a networked dynamical system to model its intrinsic activity.

Brain state:

Pattern of momentary activity that the brain can express.

Transition:

Switch from one pattern of momentary activity to another.

Control energy:

Energy to transition from one single brain state to another, modeled according to the network control theory.

Time-averaged control energy:

Average energy required to transition across a sequence of physiological states, in this case resting states.

Unimodal:

Processing of a single modality stream, for example, visual inputs.

Heteromodal:

Processing of multiple modality streams like combined visual, auditory, and proprioceptive cues.

Controllability:

Theoretical measure of a node’s capacity to spread an input throughout its network vicinity.

Intrinsic networks:

Patterns of co-activation in the cortex that pertain to resting-state activity.

Andersson
,
J. L. R.
,
Skare
,
S.
, &
Ashburner
,
J.
(
2003
).
How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging
.
NeuroImage
,
20
(
2
),
870
888
. ,
[PubMed]
Andersson
,
J. L. R.
, &
Sotiropoulos
,
S. N.
(
2016
).
An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging
.
NeuroImage
,
125
,
1063
1078
. ,
[PubMed]
Avena-Koenigsberger
,
A.
,
Mišić
,
B.
,
Hawkins
,
R. X. D.
,
Griffa
,
A.
,
Hagmann
,
P.
,
Goñi
,
J.
, &
Sporns
,
O.
(
2017
).
Path ensembles and a tradeoff between communication efficiency and resilience in the human connectome
.
Brain Structure and Function
,
222
(
1
),
603
618
. ,
[PubMed]
Bakdash
,
J. Z.
, &
Marusich
,
L. R.
(
2017
).
Repeated measures correlation
.
Frontiers in Psychology
,
8
,
456
. ,
[PubMed]
Bassett
,
D. S.
, &
Sporns
,
O.
(
2017
).
Network neuroscience
.
Nature Neuroscience
,
20
(
3
),
353
364
. ,
[PubMed]
Behzadi
,
Y.
,
Restom
,
K.
,
Liau
,
J.
, &
Liu
,
T. T.
(
2007
).
A component based noise correction method (CompCor) for BOLD and perfusion based fMRI
.
NeuroImage
,
37
(
1
),
90
101
. ,
[PubMed]
Bhushan
,
C.
,
Haldar
,
J. P.
,
Choi
,
S.
,
Joshi
,
A. A.
,
Shattuck
,
D. W.
, &
Leahy
,
R. M.
(
2015
).
Co-registration and distortion correction of diffusion and anatomical images based on inverse contrast normalization
.
NeuroImage
,
115
,
269
280
. ,
[PubMed]
Braun
,
U.
,
Harneit
,
A.
,
Pergola
,
G.
,
Menara
,
T.
,
Schäfer
,
A.
,
Betzel
,
R. F.
, …
Tost
,
H.
(
2021
).
Brain network dynamics during working memory are modulated by dopamine and diminished in schizophrenia
.
Nature Communications
,
12
(
1
),
3478
. ,
[PubMed]
Breakspear
,
M.
(
2017
).
Dynamic models of large-scale brain activity
.
Nature Neuroscience
,
20
(
3
),
340
352
. ,
[PubMed]
Bullmore
,
E.
, &
Sporns
,
O.
(
2012
).
The economy of brain network organization
.
Nature Reviews Neuroscience
,
13
(
5
),
336
349
. ,
[PubMed]
Burgos
,
N.
,
Cardoso
,
M. J.
,
Thielemans
,
K.
,
Modat
,
M.
,
Pedemonte
,
S.
,
Dickson
,
J.
, …
Ourselin
,
S.
(
2014
).
Attenuation correction synthesis for hybrid PET-MR scanners: Application to brain studies
.
IEEE Transactions on Medical Imaging
,
33
(
12
),
2332
2341
. ,
[PubMed]
Burt
,
J. B.
,
Helmer
,
M.
,
Shinn
,
M.
,
Anticevic
,
A.
, &
Murray
,
J. D.
(
2020
).
Generative modeling of brain maps with spatial autocorrelation
.
NeuroImage
,
220
,
117038
. ,
[PubMed]
Cabral
,
J.
,
Hugues
,
E.
,
Sporns
,
O.
, &
Deco
,
G.
(
2011
).
Role of local network oscillations in resting-state functional connectivity
.
NeuroImage
,
57
(
1
),
130
139
. ,
[PubMed]
Camandola
,
S.
, &
Mattson
,
M. P.
(
2017
).
Brain metabolism in health, aging, and neurodegeneration
.
EMBO Journal
,
36
(
11
),
1474
1492
. ,
[PubMed]
Castrillon
,
G.
,
Epp
,
S.
,
Bose
,
A.
,
Fraticelli
,
L.
,
Hechler
,
A.
,
Belenya
,
R.
, …
Riedl
,
V.
(
2023
).
An energy costly architecture of neuromodulators for human brain evolution and cognition
.
Science Advances
,
9
(
50
),
eadi7632
. ,
[PubMed]
Cornblath
,
E. J.
,
Ashourvan
,
A.
,
Kim
,
J. Z.
,
Betzel
,
R. F.
,
Ciric
,
R.
,
Adebimpe
,
A.
, …
Bassett
,
D. S.
(
2020
).
Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands
.
Communications Biology
,
3
(
1
),
261
. ,
[PubMed]
Craddock
,
C.
,
Sikka
,
S.
,
Cheung
,
B.
,
Khanuja
,
R.
,
Ghosh
,
S. S.
,
Yan
,
C.
, …
Milham
,
M.
(
2013
).
Towards automated analysis of connectomes: The Configurable Pipeline for the Analysis of Connectomes (C-PAC)
.
Frontiers in Neuroinformatics
.
Cui
,
Z.
,
Stiso
,
J.
,
Baum
,
G. L.
,
Kim
,
J. Z.
,
Roalf
,
D. R.
,
Betzel
,
R. F.
, …
Satterthwaite
,
T. D.
(
2020
).
Optimization of energy state transition trajectory supports the development of executive function during youth
.
eLife
,
9
,
e53060
. ,
[PubMed]
Deco
,
G.
,
Tononi
,
G.
,
Boly
,
M.
, &
Kringelbach
,
M. L.
(
2015
).
Rethinking segregation and integration: Contributions of whole-brain modelling
.
Nature Reviews Neuroscience
,
16
(
7
),
430
439
. ,
[PubMed]
Dhollander
,
T.
,
Raffelt
,
D.
, &
Connelly
,
A.
(
2016
).
Unsupervised 3-tissue response function estimation from single-shell or multi-shell diffusion MR data without a co-registered T1 image
. In
ISMRM Workshop on Breaking the Barriers of Diffusion MRI
(
Vol. 5
).
Díaz-García
,
C. M.
,
Mongeon
,
R.
,
Lahmann
,
C.
,
Koveal
,
D.
,
Zucker
,
H.
, &
Yellen
,
G.
(
2017
).
Neuronal stimulation triggers neuronal glycolysis and not lactate uptake
.
Cell Metabolism
,
26
(
2
),
361
374
. ,
[PubMed]
Dienel
,
G. A.
(
2012
).
Fueling and imaging brain activation
.
ASN Neuro
,
4
(
5
),
e00093
. ,
[PubMed]
Dienel
,
G. A.
(
2019
).
Brain glucose metabolism: Integration of energetics with function
.
Physiological Reviews
,
99
(
1
),
949
1045
. ,
[PubMed]
Eickhoff
,
S. B.
,
Thirion
,
B.
,
Varoquaux
,
G.
, &
Bzdok
,
D.
(
2015
).
Connectivity-based parcellation: Critique and implications
.
Human Brain Mapping
,
36
(
12
),
4771
4792
. ,
[PubMed]
Epp
,
S. M.
,
Castrillón
,
G.
,
Yuan
,
B.
,
Andrews-Hanna
,
J.
,
Preibisch
,
C.
, &
Riedl
,
V.
(
2024
).
Two distinct modes of hemodynamic responses across the human cortex
.
bioRxiv
.
Finc
,
K.
,
Bonna
,
K.
,
He
,
X.
,
Lydon-Staley
,
D. M.
,
Kühn
,
S.
,
Duch
,
W.
, &
Bassett
,
D. S.
(
2020
).
Dynamic reconfiguration of functional brain networks during working memory training
.
Nature Communications
,
11
(
1
),
2435
. ,
[PubMed]
Friston
,
K.
(
2010
).
The free-energy principle: A unified brain theory?
Nature Reviews Neuroscience
,
11
(
2
),
127
138
. ,
[PubMed]
Gautron
,
L.
,
Elmquist
,
J. K.
, &
Williams
,
K. W.
(
2015
).
Neural control of energy balance: Translating circuits to therapies
.
Cell
,
161
(
1
),
133
145
. ,
[PubMed]
Glasser
,
M. F.
,
Sotiropoulos
,
S. N.
,
Wilson
,
J. A.
,
Coalson
,
T. S.
,
Fischl
,
B.
,
Andersson
,
J. L.
, …
WU-Minn HCP Consortium
. (
2013
).
The minimal preprocessing pipelines for the Human Connectome Project
.
NeuroImage
,
80
,
105
124
. ,
[PubMed]
Gordon
,
E. M.
,
Laumann
,
T. O.
,
Adeyemo
,
B.
,
Huckins
,
J. F.
,
Kelley
,
W. M.
, &
Petersen
,
S. E.
(
2016
).
Generation and evaluation of a cortical area parcellation from resting-state correlations
.
Cerebral Cortex
,
26
(
1
),
288
303
. ,
[PubMed]
Gorgolewski
,
K. J.
,
Alfaro-Almagro
,
F.
,
Auer
,
T.
,
Bellec
,
P.
,
Capotă
,
M.
,
Chakravarty
,
M. M.
, …
Poldrack
,
R. A.
(
2017
).
BIDS apps: Improving ease of use, accessibility, and reproducibility of neuroimaging data analysis methods
.
PLOS Computational Biology
,
13
(
3
),
e1005209
. ,
[PubMed]
Gu
,
S.
,
Betzel
,
R. F.
,
Mattar
,
M. G.
,
Cieslak
,
M.
,
Delio
,
P. R.
,
Grafton
,
S. T.
, …
Bassett
,
D. S.
(
2017
).
Optimal trajectories of brain state transitions
.
NeuroImage
,
148
,
305
317
. ,
[PubMed]
Gu
,
S.
,
Pasqualetti
,
F.
,
Cieslak
,
M.
,
Telesford
,
Q. K.
,
Yu
,
A. B.
,
Kahn
,
A. E.
, …
Bassett
,
D. S.
(
2015
).
Controllability of structural brain networks
.
Nature Communications
,
6
,
8414
. ,
[PubMed]
Hagmann
,
P.
,
Cammoun
,
L.
,
Gigandet
,
X.
,
Meuli
,
R.
,
Honey
,
C. J.
,
Wedeen
,
V. J.
, &
Sporns
,
O.
(
2008
).
Mapping the structural core of human cerebral cortex
.
PLoS Biology
,
6
(
7
),
e159
. ,
[PubMed]
Harris
,
J. J.
,
Jolivet
,
R.
, &
Attwell
,
D.
(
2012
).
Synaptic energy use and supply
.
Neuron
,
75
(
5
),
762
777
. ,
[PubMed]
He
,
X.
,
Caciagli
,
L.
,
Parkes
,
L.
,
Stiso
,
J.
,
Karrer
,
T. M.
,
Kim
,
J. Z.
, …
Bassett
,
D. S.
(
2022
).
Uncovering the biological basis of control energy: Structural and metabolic correlates of energy inefficiency in temporal lobe epilepsy
.
Science Advances
,
8
(
45
),
eabn2293
. ,
[PubMed]
Jain
,
A. K.
(
2010
).
Data clustering: 50 years beyond K-means
.
Pattern Recognition Letters
,
31
(
8
),
651
666
.
Jeurissen
,
B.
,
Tournier
,
J.-D.
,
Dhollander
,
T.
,
Connelly
,
A.
, &
Sijbers
,
J.
(
2014
).
Multi-tissue constrained spherical deconvolution for improved analysis of multi-shell diffusion MRI data
.
NeuroImage
,
103
,
411
426
. ,
[PubMed]
Karrer
,
T. M.
,
Kim
,
J. Z.
,
Stiso
,
J.
,
Kahn
,
A. E.
,
Pasqualetti
,
F.
,
Habel
,
U.
, &
Bassett
,
D. S.
(
2020
).
A practical guide to methodological considerations in the controllability of structural brain networks
.
Journal of Neural Engineering
,
17
(
2
),
026031
. ,
[PubMed]
Kellner
,
E.
,
Dhital
,
B.
,
Kiselev
,
V. G.
, &
Reisert
,
M.
(
2016
).
Gibbs-ringing artifact removal based on local subvoxel-shifts
.
Magnetic Resonance in Medicine
,
76
(
5
),
1574
1581
. ,
[PubMed]
Kleinberg
,
J.
(
2002
).
An impossibility theorem for clustering
. In
Proceedings of the 16th International Conference on Neural Information Processing Systems
(pp.
463
470
).
MIT Press
.
Kringelbach
,
M. L.
, &
Deco
,
G.
(
2020
).
Brain states and transitions: Insights from computational neuroscience
.
Cell Reports
,
32
(
10
),
108128
. ,
[PubMed]
Leonardi
,
N.
, &
Van De Ville
,
D.
(
2015
).
On spurious and real fluctuations of dynamic functional connectivity during rest
.
NeuroImage
,
104
,
430
436
. ,
[PubMed]
Levin
,
B. E.
(
2006
).
Metabolic sensing neurons and the control of energy homeostasis
.
Physiology & Behavior
,
89
(
4
),
486
489
. ,
[PubMed]
Li
,
S.
, &
Sheng
,
Z.-H.
(
2022
).
Energy matters: Presynaptic metabolism and the maintenance of synaptic transmission
.
Nature Reviews Neuroscience
,
23
(
1
),
4
22
. ,
[PubMed]
Lord
,
L.-D.
,
Stevner
,
A. B.
,
Deco
,
G.
, &
Kringelbach
,
M. L.
(
2017
).
Understanding principles of integration and segregation using whole-brain computational connectomics: Implications for neuropsychiatric disorders
.
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
,
375
(
2096
),
20160283
. ,
[PubMed]
Luppi
,
A. I.
,
Craig
,
M. M.
,
Pappas
,
I.
,
Finoia
,
P.
,
Williams
,
G. B.
,
Allanson
,
J.
, …
Stamatakis
,
E. A.
(
2019
).
Consciousness-specific dynamic interactions of brain integration and functional diversity
.
Nature Communications
,
10
(
1
),
4616
. ,
[PubMed]
Luppi
,
A. I.
,
Rosas
,
F. E.
,
Noonan
,
M. P.
,
Mediano
,
P. A. M.
,
Kringelbach
,
M. L.
,
Carhart-Harris
,
R. L.
, …
Turkheimer
,
F. E.
(
2024
).
Oxygen and the spark of human brain evolution: Complex interactions of metabolism and cortical expansion across development and evolution
.
Neuroscientist
,
30
(
2
),
173
198
. ,
[PubMed]
Luppi
,
A. I.
,
Singleton
,
S. P.
,
Hansen
,
J. Y.
,
Jamison
,
K. W.
,
Bzdok
,
D.
,
Kuceyeski
,
A.
, …
Misic
,
B.
(
2024
).
Contributions of network structure, chemoarchitecture and diagnostic categories to transitions between cognitive topographies
.
Nature Biomedical Engineering
,
8
(
9
),
1142
1161
. ,
[PubMed]
Mantini
,
D.
,
Perrucci
,
M. G.
,
Del Gratta
,
C.
,
Romani
,
G. L.
, &
Corbetta
,
M.
(
2007
).
Electrophysiological signatures of resting state networks in the human brain
.
Proceedings of the National Academy of Sciences
,
104
(
32
),
13170
13175
. ,
[PubMed]
Margulies
,
D. S.
,
Ghosh
,
S. S.
,
Goulas
,
A.
,
Falkiewicz
,
M.
,
Huntenburg
,
J. M.
,
Langs
,
G.
, …
Smallwood
,
J.
(
2016
).
Situating the default-mode network along a principal gradient of macroscale cortical organization
.
Proceedings of the National Academy of Sciences
,
113
(
44
),
12574
12579
. ,
[PubMed]
Markello
,
R. D.
,
Hansen
,
J. Y.
,
Liu
,
Z.-Q.
,
Bazinet
,
V.
,
Shafiei
,
G.
,
Suárez
,
L. E.
, &
Misic
,
B.
(
2022
).
Neuromaps: Structural and functional interpretation of brain maps
.
Nature Methods
,
19
(
11
),
1472
1479
. ,
[PubMed]
Markello
,
R. D.
, &
Misic
,
B.
(
2021
).
Comparing spatial null models for brain maps
.
NeuroImage
,
236
,
118052
. ,
[PubMed]
Markiewicz
,
P. J.
,
Ehrhardt
,
M. J.
,
Erlandsson
,
K.
,
Noonan
,
P. J.
,
Barnes
,
A.
,
Schott
,
J. M.
, …
Ourselin
,
S.
(
2018
).
NiftyPET: A high-throughput software platform for high quantitative accuracy and precision PET imaging and analysis
.
Neuroinformatics
,
16
(
1
),
95
115
. ,
[PubMed]
Mergenthaler
,
P.
,
Lindauer
,
U.
,
Dienel
,
G. A.
, &
Meisel
,
A.
(
2013
).
Sugar for the brain: The role of glucose in physiological and pathological brain function
.
Trends in Neurosciences
,
36
(
10
),
587
597
. ,
[PubMed]
Mesulam
,
M. M.
(
1998
).
From sensation to cognition
.
Brain
,
121
(
6
),
1013
1052
. ,
[PubMed]
Miranda-Dominguez
,
O.
,
Feczko
,
E.
,
Grayson
,
D. S.
,
Walum
,
H.
,
Nigg
,
J. T.
, &
Fair
,
D. A.
(
2018
).
Heritability of the human connectome: A connectotyping study
.
Network Neuroscience
,
2
(
2
),
175
199
. ,
[PubMed]
Mišić
,
B.
,
Betzel
,
R. F.
,
Nematzadeh
,
A.
,
Goñi
,
J.
,
Griffa
,
A.
,
Hagmann
,
P.
, …
Sporns
,
O.
(
2015
).
Cooperative and competitive spreading dynamics on the human connectome
.
Neuron
,
86
(
6
),
1518
1529
. ,
[PubMed]
Muldoon
,
S. F.
,
Pasqualetti
,
F.
,
Gu
,
S.
,
Cieslak
,
M.
,
Grafton
,
S. T.
,
Vettel
,
J. M.
, &
Bassett
,
D. S.
(
2016
).
Stimulation-based control of dynamic brain networks
.
PLoS Computational Biology
,
12
(
9
),
e1005076
. ,
[PubMed]
Nozari
,
E.
,
Bertolero
,
M. A.
,
Stiso
,
J.
,
Caciagli
,
L.
,
Cornblath
,
E. J.
,
He
,
X.
, …
Bassett
,
D. S.
(
2023
).
Macroscopic resting-state brain dynamics are best described by linear models
.
Nature Biomedical Engineering
,
8
(
1
),
68
84
. ,
[PubMed]
Padamsey
,
Z.
, &
Rochefort
,
N. L.
(
2023
).
Paying the brain’s energy bill
.
Current Opinion in Neurobiology
,
78
,
102668
. ,
[PubMed]
Park
,
B.-Y.
,
Vos de Wael
,
R.
,
Paquola
,
C.
,
Larivière
,
S.
,
Benkarim
,
O.
,
Royer
,
J.
, …
Bernhardt
,
B. C.
(
2021
).
Signal diffusion along connectome gradients and inter-hub routing differentially contribute to dynamic human brain function
.
NeuroImage
,
224
,
117429
. ,
[PubMed]
Parkes
,
L.
,
Kim
,
J. Z.
,
Stiso
,
J.
,
Brynildsen
,
J. K.
,
Cieslak
,
M.
,
Covitz
,
S.
, …
Bassett
,
D. S.
(
2024
).
A network control theory pipeline for studying the dynamics of the structural connectome
.
Nature Protocols
,
19
(
12
),
3721
3749
. ,
[PubMed]
Parkes
,
L.
,
Kim
,
J. Z.
,
Stiso
,
J.
,
Calkins
,
M. E.
,
Cieslak
,
M.
,
Gur
,
R. E.
, …
Bassett
,
D. S.
(
2022
).
Asymmetric signaling across the hierarchy of cytoarchitecture within the human connectome
.
Science Advances
,
8
(
50
),
eadd2185
. ,
[PubMed]
Parkes
,
L.
,
Moore
,
T. M.
,
Calkins
,
M. E.
,
Cieslak
,
M.
,
Roalf
,
D. R.
,
Wolf
,
D. H.
, …
Bassett
,
D. S.
(
2021
).
Network controllability in transmodal cortex predicts positive psychosis spectrum symptoms
.
Biological Psychiatry
,
90
(
6
),
409
418
. ,
[PubMed]
Patenaude
,
B.
,
Smith
,
S. M.
,
Kennedy
,
D. N.
, &
Jenkinson
,
M.
(
2011
).
A Bayesian model of shape and appearance for subcortical brain segmentation
.
NeuroImage
,
56
(
3
),
907
922
. ,
[PubMed]
Patlak
,
C. S.
, &
Blasberg
,
R. G.
(
1985
).
Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. Generalizations
.
Journal of Cerebral Blood Flow & Metabolism
,
5
(
4
),
584
590
. ,
[PubMed]
Raichle
,
M. E.
, &
Gusnard
,
D. A.
(
2002
).
Appraising the brain’s energy budget
.
Proceedings of the National Academy of Sciences
,
99
(
16
),
10237
10239
. ,
[PubMed]
Robinson
,
E. C.
,
Jbabdi
,
S.
,
Glasser
,
M. F.
,
Andersson
,
J.
,
Burgess
,
G. C.
,
Harms
,
M. P.
, …
Jenkinson
,
M.
(
2014
).
MSM: A new flexible framework for multimodal surface matching
.
NeuroImage
,
100
,
414
426
. ,
[PubMed]
Saggar
,
M.
,
Shine
,
J. M.
,
Liégeois
,
R.
,
Dosenbach
,
N. U. F.
, &
Fair
,
D.
(
2022
).
Precision dynamical mapping using topological data analysis reveals a hub-like transition state at rest
.
Nature Communications
,
13
(
1
),
4791
. ,
[PubMed]
Salimi-Khorshidi
,
G.
,
Douaud
,
G.
,
Beckmann
,
C. F.
,
Glasser
,
M. F.
,
Griffanti
,
L.
, &
Smith
,
S. M.
(
2014
).
Automatic denoising of functional MRI data: Combining independent component analysis and hierarchical fusion of classifiers
.
NeuroImage
,
90
,
449
468
. ,
[PubMed]
Schaefer
,
A.
,
Kong
,
R.
,
Gordon
,
E. M.
,
Laumann
,
T. O.
,
Zuo
,
X.-N.
,
Holmes
,
A. J.
, …
Yeo
,
B. T. T.
(
2018
).
Local–global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI
.
Cerebral Cortex
,
28
(
9
),
3095
3114
. ,
[PubMed]
Shine
,
J. M.
, &
Poldrack
,
R. A.
(
2018
).
Principles of dynamic network reconfiguration across diverse brain states
.
NeuroImage
,
180
,
396
405
. ,
[PubMed]
Singleton
,
S. P.
,
Luppi
,
A. I.
,
Carhart-Harris
,
R. L.
,
Cruzat
,
J.
,
Roseman
,
L.
,
Nutt
,
D. J.
, …
Kuceyeski
,
A.
(
2022
).
Receptor-informed network control theory links LSD and psilocybin to a flattening of the brain’s control energy landscape
.
Nature Communications
,
13
(
1
),
5812
. ,
[PubMed]
Singleton
,
S. P.
,
Timmermann
,
C.
,
Luppi
,
A. I.
,
Eckernäs
,
E.
,
Roseman
,
L.
,
Carhart-Harris
,
R. L.
, &
Kuceyeski
,
A.
(
2023
).
Time-resolved network control analysis links reduced control energy under DMT with the serotonin 2a receptor, signal diversity, and subjective experience
.
bioRxiv
. ,
[PubMed]
Smith
,
R. E.
,
Tournier
,
J.-D.
,
Calamante
,
F.
, &
Connelly
,
A.
(
2012
).
Anatomically-constrained tractography: Improved diffusion MRI streamlines tractography through effective use of anatomical information
.
NeuroImage
,
62
(
3
),
1924
1938
. ,
[PubMed]
Smith
,
R. E.
,
Tournier
,
J.-D.
,
Calamante
,
F.
, &
Connelly
,
A.
(
2015
).
SIFT2: Enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography
.
NeuroImage
,
119
,
338
351
. ,
[PubMed]
Smith
,
S. M.
(
2002
).
Fast robust automated brain extraction
.
Human Brain Mapping
,
17
(
3
),
143
155
. ,
[PubMed]
Smith
,
S. M.
,
Beckmann
,
C. F.
,
Andersson
,
J.
,
Auerbach
,
E. J.
,
Bijsterbosch
,
J.
,
Douaud
,
G.
, …
WU-Minn HCP Consortium
. (
2013
).
Resting-state fMRI in the Human Connectome Project
.
NeuroImage
,
80
,
144
168
. ,
[PubMed]
Sporns
,
O.
(
2011
).
The human connectome: A complex network
.
Annals of the New York Academy of Sciences
,
1224
(
1
),
109
125
. ,
[PubMed]
Srivastava
,
P.
,
Nozari
,
E.
,
Kim
,
J. Z.
,
Ju
,
H.
,
Zhou
,
D.
,
Becker
,
C.
, …
Bassett
,
D. S.
(
2020
).
Models of communication and control for brain networks: Distinctions, convergence, and future outlook
.
Network Neuroscience
,
4
(
4
),
1122
1159
. ,
[PubMed]
Stiso
,
J.
,
Khambhati
,
A. N.
,
Menara
,
T.
,
Kahn
,
A. E.
,
Stein
,
J. M.
,
Das
,
S. R.
, …
Bassett
,
D. S.
(
2019
).
White matter network architecture guides direct electrical stimulation through optimal state transitions
.
Cell Reports
,
28
(
10
),
2554
2566
. ,
[PubMed]
Tang
,
E.
,
Giusti
,
C.
,
Baum
,
G. L.
,
Gu
,
S.
,
Pollock
,
E.
,
Kahn
,
A. E.
, …
Bassett
,
D. S.
(
2017
).
Developmental increases in white matter network controllability support a growing diversity of brain dynamics
.
Nature Communications
,
8
(
1
),
1252
. ,
[PubMed]
Thiebaut de Schotten
,
M.
, &
Forkel
,
S. J.
(
2022
).
The emergent properties of the connected brain
.
Science
,
378
(
6619
),
505
510
. ,
[PubMed]
Tournier
,
J.-D.
,
Calamante
,
F.
, &
Connelly
,
A.
(
2011
).
Effect of step size on probabilistic streamlines: Implications for the interpretation of connectivity analysis
. In
Proceedings of the International Society for Magnetic Resonance in Medicine
(
Vol. 19
, p.
2019
).
Tournier
,
J.-D.
,
Calamante
,
F.
,
Gadian
,
D. G.
, &
Connelly
,
A.
(
2004
).
Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution
.
NeuroImage
,
23
(
3
),
1176
1185
. ,
[PubMed]
Tournier
,
J.-D.
,
Smith
,
R.
,
Raffelt
,
D.
,
Tabbara
,
R.
,
Dhollander
,
T.
,
Pietsch
,
M.
, …
Connelly
,
A.
(
2019
).
MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation
.
NeuroImage
,
202
,
116137
. ,
[PubMed]
Tustison
,
N. J.
,
Avants
,
B. B.
,
Cook
,
P. A.
,
Zheng
,
Y.
,
Egan
,
A.
,
Yushkevich
,
P. A.
, &
Gee
,
J. C.
(
2010
).
N4ITK: Improved N3 bias correction
.
IEEE Transactions on Medical Imaging
,
29
(
6
),
1310
1320
. ,
[PubMed]
Vaishnavi
,
S. N.
,
Vlassenko
,
A. G.
,
Rundle
,
M. M.
,
Snyder
,
A. Z.
,
Mintun
,
M. A.
, &
Raichle
,
M. E.
(
2010
).
Regional aerobic glycolysis in the human brain
.
Proceedings of the National Academy of Sciences
,
107
(
41
),
17757
17762
. ,
[PubMed]
Van Essen
,
D. C.
,
Smith
,
S. M.
,
Barch
,
D. M.
,
Behrens
,
T. E. J.
,
Yacoub
,
E.
, &
Ugurbil
,
K.
(
2013
).
The WU-Minn Human Connectome Project: An overview
.
NeuroImage
,
80
,
62
79
. ,
[PubMed]
Váša
,
F.
, &
Mišić
,
B.
(
2022
).
Null models in network neuroscience
.
Nature Reviews Neuroscience
,
23
(
8
),
493
504
. ,
[PubMed]
Veraart
,
J.
,
Novikov
,
D. S.
,
Christiaens
,
D.
,
Ades-aron
,
B.
,
Sijbers
,
J.
, &
Fieremans
,
E.
(
2016
).
Denoising of diffusion MRI using random matrix theory
.
NeuroImage
,
142
,
394
406
. ,
[PubMed]
Vidaurre
,
D.
,
Smith
,
S. M.
, &
Woolrich
,
M. W.
(
2017
).
Brain network dynamics are hierarchically organized in time
.
Proceedings of the National Academy of Sciences
,
114
(
48
),
12827
12832
. ,
[PubMed]
Watts
,
M. E.
,
Pocock
,
R.
, &
Claudianos
,
C.
(
2018
).
Brain energy and oxygen metabolism: Emerging role in normal function and disease
.
Frontiers in Molecular Neuroscience
,
11
,
216
. ,
[PubMed]
Welvaert
,
M.
, &
Rosseel
,
Y.
(
2013
).
On the definition of signal-to-noise ratio and contrast-to-noise ratio for fMRI data
.
PLOS ONE
,
8
(
11
),
e77089
. ,
[PubMed]
Wu
,
H.-M.
,
Bergsneider
,
M.
,
Glenn
,
T. C.
,
Yeh
,
E.
,
Hovda
,
D. A.
,
Phelps
,
M. E.
, &
Huang
,
S.-C.
(
2003
).
Measurement of the global lumped constant for 2-deoxy-2-[18F]fluoro-D-glucose in normal human brain using [15O]water and 2-deoxy-2-[18F]fluoro-D-glucose positron emission tomography imaging: A method with validation based on multiple methodologies
.
Molecular Imaging & Biology
,
5
(
1
),
32
41
. ,
[PubMed]
Yang
,
J.
,
Hu
,
C.
,
Guo
,
N.
,
Dutta
,
J.
,
Vaina
,
L. M.
,
Johnson
,
K. A.
, …
Li
,
Q.
(
2017
).
Partial volume correction for PET quantification and its impact on brain network in Alzheimer’s disease
.
Scientific Reports
,
7
(
1
),
13035
. ,
[PubMed]
Yeo
,
B. T. T.
,
Krienen
,
F. M.
,
Sepulcre
,
J.
,
Sabuncu
,
M. R.
,
Lashkari
,
D.
,
Hollinshead
,
M.
, …
Buckner
,
R. L.
(
2011
).
The organization of the human cerebral cortex estimated by intrinsic functional connectivity
.
Journal of Neurophysiology
,
106
(
3
),
1125
1165
. ,
[PubMed]
Zala
,
D.
,
Hinckelmann
,
M.-V.
,
Yu
,
H.
,
Lyra da Cunha
,
M. M.
,
Liot
,
G.
,
Cordelières
,
F. P.
, …
Saudou
,
F.
(
2013
).
Vesicular glycolysis provides on-board energy for fast axonal transport
.
Cell
,
152
(
3
),
479
491
. ,
[PubMed]
Zalesky
,
A.
, &
Breakspear
,
M.
(
2015
).
Towards a statistical test for functional connectivity dynamics
.
NeuroImage
,
114
,
466
470
. ,
[PubMed]
Zalesky
,
A.
,
Fornito
,
A.
,
Cocchi
,
L.
,
Gollo
,
L. L.
, &
Breakspear
,
M.
(
2014
).
Time-resolved resting-state brain networks
.
Proceedings of the National Academy of Sciences
,
111
(
28
),
10341
10346
. ,
[PubMed]
Zamani Esfahlani
,
F.
,
Jo
,
Y.
,
Faskowitz
,
J.
,
Byrge
,
L.
,
Kennedy
,
D. P.
,
Sporns
,
O.
, &
Betzel
,
R. F.
(
2020
).
High-amplitude cofluctuations in cortical activity drive functional connectivity
.
Proceedings of the National Academy of Sciences
,
117
(
45
),
28393
28401
. ,
[PubMed]
Zhang
,
Y.
,
Brady
,
M.
, &
Smith
,
S.
(
2001
).
Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm
.
IEEE Transactions on Medical Imaging
,
20
(
1
),
45
57
. ,
[PubMed]

Competing Interests

Competing Interests: The authors have declared that no competing interests exist.

Author notes

Handling Editor: Claus Hilgetag

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a full description of the license, please visit https://creativecommons.org/licenses/by/4.0/legalcode.

Supplementary data