Resting-state background features demonstrate multidien cycles in long-term EEG device recordings

BACKGROUND

Longitudinal EEG recorded by implanted devices is critical for understanding and managing epilepsy. Recent research reports patient-specific, multi-day cycles in device-detected epileptiform events that coincide with increased likelihood of clinical seizures. Understanding these cycles could elucidate mechanisms generating seizures and advance drug and neurostimulation therapies.

OBJECTIVE/HYPOTHESIS

We hypothesize that seizure-correlated cycles are present in background neural activity, independent of interictal epileptiform spikes, and that neurostimulation may temporarily interrupt these cycles.

METHODS

We analyzed regularly-recorded seizure-free data epochs from 20 patients implanted with a responsive neurostimulation (RNS) device for at least 1.5 years, to explore the relationship between cycles in device-detected interictal epileptiform activity (dIEA), clinician-validated interictal spikes, background EEG features, and neurostimulation.

RESULTS

Background EEG features tracked the cycle phase of dIEA in all patients (AUC: 0.63 [0.56-0.67]) with a greater effect size compared to clinically annotated spike rate alone (AUC: 0.55 [0.53-0.61], p < 0.01). After accounting for circadian variation and spike rate, we observed significant population trends in elevated theta and beta band power and theta and alpha connectivity features at the cycle peaks (sign test, p < 0.05). In the period directly after stimulation we observe a decreased association between cycle phase and EEG features compared to background recordings (AUC: 0.58 [0.55-0.64]).

CONCLUSIONS

Our findings suggest that seizure-correlated dIEA cycles are not solely due to epileptiform discharges but are associated with background measures of brain state; and that neurostimulation may temporarily interrupt these cycles. These results may help elucidate mechanisms underlying seizure generation, provide new biomarkers for seizure risk, and facilitate monitoring, treating, and managing epilepsy with implantable devices.

Resting-state background features demonstrate multidien cycles in long-term EEG device recordings

Background

Longitudinal EEG recorded by implanted devices is critical for understanding and managing epilepsy. Recent research reports patient-specific, multi-day cycles in device-detected epileptiform events that coincide with increased likelihood of clinical seizures. Understanding these cycles could elucidate mechanisms generating seizures and advance drug and neurostimulation therapies.

Objective/Hypothesis

We hypothesize that seizure-correlated cycles are present in background neural activity, independent of interictal epileptiform spikes, and that neurostimulation may disrupt these cycles.

Methods

We analyzed regularly-recorded seizure-free data epochs from 20 patients implanted with a responsive neurostimulation (RNS) device for at least 1.5 years, to explore the relationship between cycles in device-detected interictal epileptiform activity (dIEA), clinician-validated interictal spikes, background EEG features, and neurostimulation.

Results

Background EEG features tracked the cycle phase of dIEA in all patients (AUC: 0.63 [0.56 - 0.67]) with a greater effect size compared to clinically annotated spike rate alone (AUC: 0.55 [0.53-0.61], p < 0.01). After accounting for circadian variation and spike rate, we observed significant population trends in elevated theta and beta band power and theta and alpha connectivity features at the cycle peaks (sign test, p < 0.05). In the period directly after stimulation we observe a decreased association between cycle phase and EEG features compared to background recordings (AUC: 0.58 [0.55-0.64]).

Conclusions

Our findings suggest that seizure-correlated dIEA cycles are not solely due to epileptiform discharges but are associated with background measures of brain state; and that neurostimulation may disrupt these cycles. These results may help elucidate mechanisms underlying seizure generation, provide new biomarkers for seizure risk, and facilitate monitoring, treating, and managing epilepsy with implantable devices.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Thermal effects on neurons during stimulation of the brain

All electric and magnetic stimulation of the brain deposits thermal energy in the brain. This occurs through either Joule heating of the conductors carrying current through electrodes and magnetic coils, or through dissipation of energy in the conductive brain.Objective.Although electrical interaction with brain tissue is inseparable from thermal effects when electrodes are used, magnetic induction enables us to separate Joule heating from induction effects by contrasting AC and DC driving of magnetic coils using the same energy deposition within the conductors. Since mammalian cortical neurons have no known sensitivity to static magnetic fields, and if there is no evidence of effect on spike timing to oscillating magnetic fields, we can presume that the induced electrical currents within the brain are below the molecular shot noise where any interaction with tissue is purely thermal.Approach.In this study, we examined a range of frequencies produced from micromagnetic coils operating below the molecular shot noise threshold for electrical interaction with single neurons.Main results.We found that small temperature increases and decreases of 1∘C caused consistent transient suppression and excitation of neurons during temperature change. Numerical modeling of the biophysics demonstrated that the Na-K pump, and to a lesser extent the Nernst potential, could account for these transient effects. Such effects are dependent upon compartmental ion fluxes and the rate of temperature change.Significance.A new bifurcation is described in the model dynamics that accounts for the transient suppression and excitation; in addition, we note the remarkable similarity of this bifurcation's rate dependency with other thermal rate-dependent tipping points in planetary warming dynamics. These experimental and theoretical findings demonstrate that stimulation of the brain must take into account small thermal effects that are ubiquitously present in electrical and magnetic stimulation. More sophisticated models of electrical current interaction with neurons combined with thermal effects will lead to more accurate modulation of neuronal activity.

External drivers of BOLD signal's non-stationarity

A fundamental challenge in neuroscience is to uncover the principles governing how the brain interacts with the external environment. However, assumptions about external stimuli fundamentally constrain current computational models. We show in silico that unknown external stimulation can produce error in the estimated linear time-invariant dynamical system. To address these limitations, we propose an approach to retrieve the external (unknown) input parameters and demonstrate that the estimated system parameters during external input quiescence uncover spatiotemporal profiles of external inputs over external stimulation periods more accurately. Finally, we unveil the expected (and unexpected) sensory and task-related extra-cortical input profiles using functional magnetic resonance imaging data acquired from 96 subjects (Human Connectome Project) during the resting-state and task scans. This dynamical systems model of the brain offers information on the structure and dimensionality of the BOLD signal's external drivers and shines a light on the likely external sources contributing to the BOLD signal's non-stationarity. Our findings show the role of exogenous inputs in the BOLD dynamics and highlight the importance of accounting for external inputs to unravel the brain's time-varying functional dynamics.

Model-based stationarity filtering of long-term memory data applied to resting-state blood-oxygen-level-dependent signal

Resting-state blood-oxygen-level-dependent (BOLD) signal acquired through functional magnetic resonance imaging is a proxy of neural activity and a key mechanism for assessing neurological conditions. Therefore, practical tools to filter out artefacts that can compromise the assessment are required. On the one hand, a variety of tailored methods to preprocess the data to deal with identified sources of noise (e.g., head motion, heart beating, and breathing, just to mention a few) are in place. But, on the other hand, there might be unknown sources of unstructured noise present in the data. Therefore, to mitigate the effects of such unstructured noises, we propose a model-based filter that explores the statistical properties of the underlying signal (i.e., long-term memory). Specifically, we consider autoregressive fractional integrative process filters. Remarkably, we provide evidence that such processes can model the signals at different regions of interest to attain stationarity. Furthermore, we use a principled analysis where a ground-truth signal with statistical properties similar to the BOLD signal under the injection of noise is retrieved using the proposed filters. Next, we considered preprocessed (i.e., the identified sources of noise removed) resting-state BOLD data of 98 subjects from the Human Connectome Project. Our results demonstrate that the proposed filters decrease the power in the higher frequencies. However, unlike the low-pass filters, the proposed filters do not remove all high-frequency information, instead they preserve process-related higher frequency information. Additionally, we considered four different metrics (power spectrum, functional connectivity using the Pearson’s correlation, coherence, and eigenbrains) to infer the impact of such filter. We provided evidence that whereas the first three keep most of the features of interest from a neuroscience perspective unchanged, the latter exhibits some variations that could be due to the sporadic activity filtered out.

Pairwise maximum entropy model explains the role of white matter structure in shaping emergent co-activation states

A major challenge in neuroscience is determining a quantitative relationship between the brain's white matter structural connectivity and emergent activity. We seek to uncover the intrinsic relationship among brain regions fundamental to their functional activity by constructing a pairwise maximum entropy model (MEM) of the inter-ictal activation patterns of five patients with medically refractory epilepsy over an average of ~14 hours of band-passed intracranial EEG (iEEG) recordings per patient. We find that the pairwise MEM accurately predicts iEEG electrodes' activation patterns' probability and their pairwise correlations. We demonstrate that the estimated pairwise MEM's interaction weights predict structural connectivity and its strength over several frequencies significantly beyond what is expected based solely on sampled regions' distance in most patients. Together, the pairwise MEM offers a framework for explaining iEEG functional connectivity and provides insight into how the brain's structural connectome gives rise to large-scale activation patterns by promoting co-activation between connected structures.

Time-evolving controllability of effective connectivity networks during seizure progression

Over one third of the estimated 3 million people with epilepsy in the United States are medication resistant. Responsive neurostimulation from chronically implanted electrodes provides a promising treatment alternative to resective surgery. However, determining optimal personalized stimulation parameters, including when and where to intervene to guarantee a positive patient outcome, is a major open challenge. Network neuroscience and control theory offer useful tools that may guide improvements in parameter selection for control of anomalous neural activity. Here we use a method to characterize dynamic controllability across consecutive effective connectivity (EC) networks based on regularized partial correlations between implanted electrodes during the onset, propagation, and termination regimes of 34 seizures. We estimate regularized partial correlation adjacency matrices from 1-s time windows of intracranial electrocorticography recordings using the Graphical Least Absolute Shrinkage and Selection Operator (GLASSO). Average and modal controllability metrics calculated from each resulting EC network track the time-varying controllability of the brain on an evolving landscape of conditionally dependent network interactions. We show that average controllability increases throughout a seizure and is negatively correlated with modal controllability throughout. Our results support the hypothesis that the energy required to drive the brain to a seizure-free state from an ictal state is smallest during seizure onset, yet we find that applying control energy at electrodes in the seizure onset zone may not always be energetically favorable. Our work suggests that a low-complexity model of time-evolving controllability may offer insights for developing and improving control strategies targeting seizure suppression.

Multimodal in vivo recording using transparent graphene microelectrodes illuminates spatiotemporal seizure dynamics at the microscale

Neurological disorders such as epilepsy arise from disrupted brain networks. Our capacity to treat these disorders is limited by our inability to map these networks at sufficient temporal and spatial scales to target interventions. Current best techniques either sample broad areas at low temporal resolution (e.g. calcium imaging) or record from discrete regions at high temporal resolution (e.g. electrophysiology). This limitation hampers our ability to understand and intervene in aberrations of network dynamics. Here we present a technique to map the onset and spatiotemporal spread of acute epileptic seizures in vivo by simultaneously recording high bandwidth microelectrocorticography and calcium fluorescence using transparent graphene microelectrode arrays. We integrate dynamic data features from both modalities using non-negative matrix factorization to identify sequential spatiotemporal patterns of seizure onset and evolution, revealing how the temporal progression of ictal electrophysiology is linked to the spatial evolution of the recruited seizure core. This integrated analysis of multimodal data reveals otherwise hidden state transitions in the spatial and temporal progression of acute seizures. The techniques demonstrated here may enable future targeted therapeutic interventions and novel spatially embedded models of local circuit dynamics during seizure onset and evolution.

Fractional-order model predictive control as a framework for electrical neurostimulation in epilepsy

Electrical neurostimulation is an increasingly adopted therapeutic methodology for neurological conditions such as epilepsy. Electrical neurostimulation devices are commonly characterized by their limited sensing, actuating, and computational capabilities. However, the sensing mechanisms are often used only for their detection potential (e.g., to detect seizures), which automatically and dynamically trigger the actuation capabilities, but ultimately deploy prespecified stimulation doses that resulted from a period of manual (and empirical) calibration. The potential information contained in the measurements acquired by the sensing mechanisms is, therefore, considerably underutilized, given that this type of stimulation strategy only entails an event-triggered relationship between the sensors and actuators of the device. Such stimulation strategies are suboptimal in general and lack theoretical guarantees regarding their performance. In order to leverage the aforementioned information, harvested during normal sensing-actuating operation, we must consider a real-time feedback (closed-loop) strategy. More precisely, the stimulation signal itself should automatically adapt based upon the state of the neurophysiological system at hand, estimated from data collected in real-time through sensors in the device. In this work, we propose a model-based approach for (real-time) closed-loop electrical neurostimulation, in which the evolution of the system is captured by a fractional-order system (FOS). More precisely, we propose a model predictive control (MPC) approach with an underlying FOS predictive model, due to the ability of fractional-order dynamics to more accurately capture the long-term dependence present in biological systems, compared to the standard linear time-invariant models. Furthermore, MPC offers, by design, an additional layer of robustness to compensate for system-model mismatch, which the more traditional strategies lack. To establish the potential of our framework, we focus on epileptic seizure mitigation by computational simulation of our proposed strategy upon seizure-like events. Lastly, we provide evidence of the effectiveness of our method on seizures simulated by commonly adopted models in the neuroscience and medical community present in the literature, as well as real seizure data as obtained from subjects with epilepsy.

Reconfigurations within resonating communities of brain regions following TMS reveal different scales of processing

An overarching goal of neuroscience research is to understand how heterogeneous neuronal ensembles cohere into networks of coordinated activity to support cognition. To investigate how local activity harmonizes with global signals, we measured electroencephalography (EEG) while single pulses of transcranial magnetic stimulation (TMS) perturbed occipital and parietal cortices. We estimate the rapid network reconfigurations in dynamic network communities within specific frequency bands of the EEG, and characterize two distinct features of network reconfiguration, flexibility and allegiance, among spatially distributed neural sources following TMS. Using distance from the stimulation site to infer local and global effects, we find that alpha activity (8–12 Hz) reflects concurrent local and global effects on network dynamics. Pairwise allegiance of brain regions to communities on average increased near the stimulation site, whereas TMS-induced changes to flexibility were generally invariant to distance and stimulation site. In contrast, communities within the beta (13–20 Hz) band demonstrated a high level of spatial specificity, particularly within a cluster comprising paracentral areas. Together, these results suggest that focal magnetic neurostimulation to distinct cortical sites can help identify both local and global effects on brain network dynamics, and highlight fundamental differences in the manifestation of network reconfigurations within alpha and beta frequency bands.

TMS may be used to probe the causal link between local regional activity and global brain dynamics. Using simultaneous TMS-EEG and dynamic community detection, we introduce what we call “resonating communities” or frequency band-specific clusters in the brain, as a way to index local and global processing. These resonating communities within the alpha and beta bands display both global (or integrating) behavior and local specificity, highlighting fundamental differences in the manifestation of network reconfigurations.

Impact of childhood adversity on network reconfiguration dynamics during working memory in hypogonadal women

Many women with no history of cognitive difficulties experience executive dysfunction during menopause. Significant adversity during childhood negatively impacts executive function into adulthood and may be an indicator of women at risk of a mid-life cognitive decline. Previous studies have indicated that alterations in functional network connectivity underlie these negative effects of childhood adversity. There is growing evidence that functional brain networks are not static during executive tasks; instead, such networks reconfigure over time. Optimal dynamics are necessary for efficient executive function; while too little reconfiguration is insufficient for peak performance, too much reconfiguration (supra-optimal reconfiguration) is also maladaptive and associated with poorer performance. Here we examined the impact of adverse childhood experiences (ACEs) on network flexibility, a measure of dynamic reconfiguration, during a letter n-back task within three networks that support executive function: frontoparietal, salience, and default mode networks. Several animal and human subject studies have suggested that childhood adversity exerts lasting effects on executive function via serotonergic mechanisms. Tryptophan depletion (TD) was used to examine whether serotonin function drives ACE effects on network flexibility. We hypothesized that ACE would be associated with higher flexibility (supra-optimal flexibility) and that TD would further increase this measure. Forty women underwent functional imaging at two time points in this double-blind, placebo controlled, crossover study. Participants also completed the Penn Conditional Exclusion Test, a task assessing abstraction and mental flexibility. The effects of ACE and TD were evaluated using generalized estimating equations. ACE was associated with higher flexibility across networks (frontoparietal β = 0.00748, D = 2.79, p = 0.005; salience β = 0.00679, D = 3.02, p = 0.003; and default mode β = 0.00910, D = 3.53, p = 0.0004). While there was no interaction between ACE and TD, active TD increased network flexibility in both ACE groups in comparison to sham depletion (frontoparietal β = 0.00489, D = 2.15, p = 0.03; salience β = 0.00393, D = 1.91, p = 0.06; default mode β = 0.00334, D = 1.73, p = 0.08). These results suggest that childhood adversity has lasting impacts on dynamic reconfiguration of functional brain networks supporting executive function and that decreasing serotonin levels may exacerbate these effects.

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.

Model-based design for seizure control by stimulation

OBJECTIVE

Current brain stimulation paradigms are largely empirical rather than theoretical. An opportunity exists to improve upon their modest effectiveness in closed-loop control strategies with the development of theoretically grounded, model-based designs.

APPROACH

Inspired by this need, here we couple experimental data and mathematical modeling with a control-theoretic strategy for seizure termination. We begin by exercising a dynamical systems approach to model seizures (n = 94) recorded using intracranial EEG (iEEG) from 21 patients with medication-resistant, localization-related epilepsy.

MAIN RESULTS

Although each patient's seizures displayed unique spatial and temporal patterns, their evolution can be parsimoniously characterized by the same model form. Idiosyncracies of the model can inform individualized intervention strategies, specifically in iEEG samples with well-localized seizure onset zones. Temporal fluctuations in the spatial profiles of the oscillatory modes show that seizure onset marks a transition into a regime in which the underlying system supports prolonged rhythmic and focal activity. Based on these observations, we propose a control-theoretic strategy that aims to stabilize ictal activity using static output feedback for linear time-invariant switching systems. Finally, we demonstrate in silico that our proposed strategy allows us to dampen the emerging focal oscillatory sources using only a small set of electrodes.

SIGNIFICANCE

Our integrative study informs the development of modulation and control algorithms for neurostimulation that could improve the effectiveness of implantable, closed-loop anti-epileptic devices.

Characterizing the role of the structural connectome in seizure dynamics

How does the human brain's structural scaffold give rise to its intricate functional dynamics? This is a central question in translational neuroscience that is particularly relevant to epilepsy, a disorder affecting over 50 million subjects worldwide. Treatment for medication-resistant focal epilepsy is often structural-through surgery or laser ablation-but structural targets, particularly in patients without clear lesions, are largely based on functional mapping via intracranial EEG. Unfortunately, the relationship between structural and functional connectivity in the seizing brain is poorly understood. In this study, we quantify structure-function coupling, specifically between white matter connections and intracranial EEG, across pre-ictal and ictal periods in 45 seizures from nine patients with unilateral drug-resistant focal epilepsy. We use high angular resolution diffusion imaging (HARDI) tractography to construct structural connectivity networks and correlate these networks with time-varying broadband and frequency-specific functional networks derived from coregistered intracranial EEG. Across all frequency bands, we find significant increases in structure-function coupling from pre-ictal to ictal periods. We demonstrate that short-range structural connections are primarily responsible for this increase in coupling. Finally, we find that spatiotemporal patterns of structure-function coupling are highly stereotyped for each patient. These results suggest that seizures harness the underlying structural connectome as they propagate. Mapping the relationship between structural and functional connectivity in epilepsy may inform new therapies to halt seizure spread, and pave the way for targeted patient-specific interventions.

High interictal connectivity within the resection zone is associated with favorable post-surgical outcomes in focal epilepsy patients

Patients with drug-resistant focal epilepsy are often candidates for invasive surgical therapies. In these patients, it is necessary to accurately localize seizure generators to ensure seizure freedom following intervention. While intracranial electroencephalography (iEEG) is the gold standard for mapping networks for surgery, this approach requires inducing and recording seizures, which may cause patient morbidity. The goal of this study is to evaluate the utility of mapping interictal (non-seizure) iEEG networks to identify targets for surgical treatment. We analyze interictal iEEG recordings and neuroimaging from 27 focal epilepsy patients treated via surgical resection. We generate interictal functional networks by calculating pairwise correlation of iEEG signals across different frequency bands. Using image coregistration and segmentation, we identify electrodes falling within surgically resected tissue (i.e. the resection zone), and compute node-level and edge-level synchrony in relation to the resection zone. We further associate these metrics with post-surgical outcomes. Greater overlap between resected electrodes and highly synchronous electrodes is associated with favorable post-surgical outcomes. Additionally, good-outcome patients have significantly higher connectivity localized within the resection zone compared to those with poorer postoperative seizure control. This finding persists following normalization by a spatially-constrained null model. This study suggests that spatially-informed interictal network synchrony measures can distinguish between good and poor post-surgical outcomes. By capturing clinically-relevant information during interictal periods, our method may ultimately reduce the need for prolonged invasive implants and provide insights into the pathophysiology of an epileptic brain. We discuss next steps for translating these findings into a prospectively useful clinical tool.

•

We analyze interictal iEEG recordings and neuroimaging from epilepsy patients

•

We find that high interictal strength selectivity is associated with better outcomes

•

This effect appears to be driven largely by connectivity within the resection zone

•

Interictal recordings can guide identification of seizure-generating networks

Active spectroscopy measurements of the deuterium temperature, rotation, and density from the core to scrape off layer on the DIII-D tokamak (invited)

Main-ion charge exchange recombination spectroscopy (MICER) uses the neutral beam induced D _α spectrum to measure the local deuterium ion (D⁺) temperature, rotation, and density, as well as parameters related to the neutral beams, fast ions, and magnetic field. An edge MICER system consisting of 16 densely packed chords was recently installed on DIII-D, extending the MICER technique from the core to the pedestal and steep gradient region of H-mode plasmas where the D⁺ and commonly measured impurity ion properties can differ significantly. A combination of iterative collisional radiative modeling techniques and greatly accelerated spectral fitting allowed the extension of this diagnostic technique to the plasma edge where the steep gradients introduce significant diagnostic challenges. The importance of including the fast ion D _α emission in the fit to the spectrum for the edge system is investigated showing that it is typically not important except for cases which can have significant fast ion fractions near the plasma edge such as QH-mode. Example profiles from an Ohmic L-mode and a high power ITER baseline case show large differences in the toroidal rotation of the two species near the separatrix including a strong co-current D⁺ edge rotation. The measurements and analysis demonstrate the state of the art in active spectroscopy and integrated modeling for diagnosing fusion plasmas and the importance of direct main ion measurements.

Applications of community detection techniques to brain graphs: Algorithmic considerations and implications for neural function

The human brain can be represented as a graph in which neural units such as cells or small volumes of tissue are heterogeneously connected to one another through structural or functional links. Brain graphs are parsimonious representations of neural systems that have begun to offer fundamental insights into healthy human cognition, as well as its alteration in disease. A critical open question in network neuroscience lies in how neural units cluster into densely interconnected groups that can provide the coordinated activity that is characteristic of perception, action, and adaptive behaviors. Tools that have proven particularly useful for addressing this question are community detection approaches, which can identify communities or modules: groups of neural units that are densely interconnected with other units in their own group but sparsely interconnected with units in other groups. In this paper, we describe a common community detection algorithm known as modularity maximization, and we detail its applications to brain graphs constructed from neuroimaging data. We pay particular attention to important algorithmic considerations, especially in recent extensions of these techniques to graphs that evolve in time. After recounting a few fundamental insights that these techniques have provided into brain function, we highlight potential avenues of methodological advancements for future studies seeking to better characterize the patterns of coordinated activity in the brain that accompany human behavior. This tutorial provides a naive reader with an introduction to theoretical considerations pertinent to the generation of brain graphs, an understanding of modularity maximization for community detection, a resource of statistical measures that can be used to characterize community structure, and an appreciation of the usefulness of these approaches in uncovering behaviorally-relevant network dynamics in neuroimaging data.

Cohesive network reconfiguration accompanies extended training

Human behavior is supported by flexible neurophysiological processes that enable the fine‐scale manipulation of information across distributed neural circuits. Yet, approaches for understanding the dynamics of these circuit interactions have been limited. One promising avenue for quantifying and describing these dynamics lies in multilayer network models. Here, networks are composed of nodes (which represent brain regions) and time‐dependent edges (which represent statistical similarities in activity time series). We use this approach to examine functional connectivity measured by non‐invasive neuroimaging techniques. These multilayer network models facilitate the examination of changes in the pattern of statistical interactions between large‐scale brain regions that might facilitate behavior. In this study, we define and exercise two novel measures of network reconfiguration, and demonstrate their utility in neuroimaging data acquired as healthy adult human subjects learn a new motor skill. In particular, we identify putative functional modules in multilayer networks and characterize the degree to which nodes switch between modules. Next, we define cohesive switches, in which a set of nodes moves between modules together as a group, and we define disjoint switches, in which a single node moves between modules independently from other nodes. Together, these two concepts offer complementary yet distinct insights into the changes in functional connectivity that accompany motor learning. More generally, our work offers statistical tools that other researchers can use to better understand the reconfiguration patterns of functional connectivity over time. Hum Brain Mapp 38:4744–4759, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

The energy landscape underpinning module dynamics in the human brain connectome

Human brain dynamics can be viewed through the lens of statistical mechanics, where neurophysiological activity evolves around and between local attractors representing mental states. Many physically-inspired models of these dynamics define brain states based on instantaneous measurements of regional activity. Yet, recent work in network neuroscience has provided evidence that the brain might also be well-characterized by time-varying states composed of locally coherent activity or functional modules. We study this network-based notion of brain state to understand how functional modules dynamically interact with one another to perform cognitive functions. We estimate the functional relationships between regions of interest (ROIs) by fitting a pair-wise maximum entropy model to each ROI's pattern of allegiance to functional modules. This process uses an information theoretic notion of energy (as opposed to a metabolic one) to produce an energy landscape in which local minima represent attractor states characterized by specific patterns of modular structure. The clustering of local minima highlights three classes of ROIs with similar patterns of allegiance to community states. Visual, attention, sensorimotor, and subcortical ROIs are well-characterized by a single functional community. The remaining ROIs affiliate with a putative executive control community or a putative default mode and salience community. We simulate the brain's dynamic transitions between these community states using a random walk process. We observe that simulated transition probabilities between basins are statistically consistent with empirically observed transitions in resting state fMRI data. These results offer a view of the brain as a dynamical system that transitions between basins of attraction characterized by coherent activity in groups of brain regions, and that the strength of these attractors depends on the ongoing cognitive computations.