Beyond modularity: Fine-scale mechanisms and rules for brain network reconfiguration

Extended nonnegative matrix factorization for dynamic functional connectivity analysis of fMRI data

Dynamic functional connectivity (DFC) analysis using functional magnetic resonance imaging (fMRI) technology has attracted increasing attention in revealing brain dynamics in recent years. Although the nonnegative matrix factorization (NMF) method was applied to dynamic subgraph analysis to reveal brain dynamics, its application in DFC analysis was largely limited due to its nonnegative constraint on the input data. This study proposed the extended NMF (eNMF) method that allowed the input matrix and decomposed basis matrix to have negative values without altering the NMF algorithm. The eNMF method was applied to DFC analysis of both simulated and real resting fMRI data. The simulated data demonstrated that eNMF successfully decomposed the mixed-sign matrix into one positive matrix and one mixed-sign matrix. In contrast to K-means, eNMF extracted more accurate brain state patterns in all cases and estimated better DFC temporal properties for uneven brain state distribution. The real resting-fMRI data demonstrated that eNMF can provide more temporal measures of DFC and was more sensitive to detect intergroup differences of DFC than K-means. Results of eNMF revealed that the female group possibly showed worse relaxation and produced stronger spontaneous cognitive processes although they tended to spend more time in relaxation state and less time in states relevant to cognitive processes in contrast to the male group.

The Profile of Network Spontaneous Activity and Functional Organization Interplay in Hierarchically Connected Modular Neural Networks In Vitro

Modern microtechnology methods are widely used to create neural networks on a chip with a connection architecture demonstrating properties of modularity and hierarchy similar to brain networks. Such in vitro networks serve as a valuable model for studying the interplay of functional architecture within modules, their activity, and the effectiveness of inter-module interaction. In this study, we use a two-chamber microfluidic platform to investigate functional connectivity and global activity in hierarchically connected modular neural networks. We found that the strength of functional connections within the module and the profile of network spontaneous activity determine the effectiveness of inter-modular interaction and integration activity in the network. The direction of intermodular activity propagation configures the different densities of inhibitory synapses in the network. The developed microfluidic platform holds the potential to explore function-structure relationships and efficient information processing in two- or multilayer neural networks, in both healthy and pathological states.

A dynamic brain network decomposition method discovers effective brain hemodynamic sub-networks for Parkinson's disease

Objective.Dopaminergic treatment is effective for Parkinson's disease (PD). Nevertheless, the conventional treatment assessment mainly focuses on human-administered behavior examination while the underlying functional improvements have not been well explored. This paper aims to investigate brain functional variations of PD patients after dopaminergic therapy.Approach.This paper proposed a dynamic brain network decomposition method and discovered brain hemodynamic sub-networks that well characterized the efficacy of dopaminergic treatment in PD. Firstly, a clinical walking procedure with functional near-infrared spectroscopy was developed, and brain activations during the procedure from fifty PD patients under the OFF and ON states (without and with dopaminergic medication) were captured. Then, dynamic brain networks were constructed with sliding-window analysis of phase lag index and integrated time-varying functional networks across all patients. Afterwards, an aggregated network decomposition algorithm was formulated based on aggregated effectiveness optimization of functional networks in spanning network topology and cross-validation network variations, and utilized to unveil effective brain hemodynamic sub-networks for PD patients. Further, dynamic sub-network features were constructed to characterize the brain flexibility and dynamics according to the temporal switching and activation variations of discovered sub-networks, and their correlations with differential treatment-induced gait alterations were analyzed.Results.The results demonstrated that PD patients exhibited significantly enhanced flexibility after dopaminergic therapy within a sub-network related to the improvement of motor functions. Other sub-networks were significantly correlated with trunk-related axial symptoms and exhibited no significant treatment-induced dynamic interactions.Significance.The proposed method promises a quantified and objective approach for dopaminergic treatment evaluation. Moreover, the findings suggest that the gait of PD patients comprises distinct motor domains, and the corresponding neural controls are selectively responsive to dopaminergic treatment.

Longitudinal regression of covariance matrix outcomes

In this study, a longitudinal regression model for covariance matrix outcomes is introduced. The proposal considers a multilevel generalized linear model for regressing covariance matrices on (time-varying) predictors. This model simultaneously identifies covariate-associated components from covariance matrices, estimates regression coefficients, and captures the within-subject variation in the covariance matrices. Optimal estimators are proposed for both low-dimensional and high-dimensional cases by maximizing the (approximated) hierarchical-likelihood function. These estimators are proved to be asymptotically consistent, where the proposed covariance matrix estimator is the most efficient under the low-dimensional case and achieves the uniformly minimum quadratic loss among all linear combinations of the identity matrix and the sample covariance matrix under the high-dimensional case. Through extensive simulation studies, the proposed approach achieves good performance in identifying the covariate-related components and estimating the model parameters. Applying to a longitudinal resting-state functional magnetic resonance imaging data set from the Alzheimer's Disease (AD) Neuroimaging Initiative, the proposed approach identifies brain networks that demonstrate the difference between males and females at different disease stages. The findings are in line with existing knowledge of AD and the method improves the statistical power over the analysis of cross-sectional data.

Estimating Functional Brain Networks by Low-Rank Representation With Local Constraint

The functional architecture undergoes alterations during the preclinical phase of Alzheimer's disease. Consequently, the primary research focus has shifted towards identifying Alzheimer's disease and its early stages by constructing a functional connectivity network based on resting-state fMRI data. Recent investigations show that as Alzheimer's Disease (AD) progresses, modular tissue and connections in the core brain areas of AD patients diminish. Sparse learning methods are powerful tools for understanding Functional Brain Networks (FBNs) with Regions of Interest (ROIs) and a connectivity matrix measuring functional coherence between them. However, these tools often focus exclusively on functional connectivity measures, neglecting the brain network's modularity. Modularity orchestrates dynamic activities within the FBN to execute intricate cognitive tasks. To provide a comprehensive delineation of the FBN, we propose a local similarity-constrained low-rank sparse representation (LSLRSR) method that encodes modularity information under a manifold-regularized network learning framework and further formulate it as a low-rank sparse graph learning problem, which can be solved by an efficient optimization algorithm. Specifically, for each modularity structure, the Schatten p-norm regularizer reduces the reconstruction error and provides a better approximation of the low-rank constraint. Furthermore, we adopt a manifold-regularized local similarity prior to infer the intricate relationship between subnetwork similarity and modularity, guiding the modeling of FBN. Additionally, the proximal average method approximates the joint solution's proximal map, and the resulting nonconvex optimization problems are solved using the alternating direction multiplier method (ADMM). Compared to state-of-the-art methods for constructing FBNs, our algorithm generates a more modular FBN. This lays the groundwork for further research into alterations in brain network modularity resulting from diseases.

Network targets for therapeutic brain stimulation: towards personalized therapy for pain

Precision neuromodulation of central brain circuits is a promising emerging therapeutic modality for a variety of neuropsychiatric disorders. Reliably identifying in whom, where, and in what context to provide brain stimulation for optimal pain relief are fundamental challenges limiting the widespread implementation of central neuromodulation treatments for chronic pain. Current approaches to brain stimulation target empirically derived regions of interest to the disorder or targets with strong connections to these regions. However, complex, multidimensional experiences like chronic pain are more closely linked to patterns of coordinated activity across distributed large-scale functional networks. Recent advances in precision network neuroscience indicate that these networks are highly variable in their neuroanatomical organization across individuals. Here we review accumulating evidence that variable central representations of pain will likely pose a major barrier to implementation of population-derived analgesic brain stimulation targets. We propose network-level estimates as a more valid, robust, and reliable way to stratify personalized candidate regions. Finally, we review key background, methods, and implications for developing network topology-informed brain stimulation targets for chronic pain.

Evolution of brain network dynamics in early Parkinson's disease with mild cognitive impairment

How mild cognitive impairment (MCI) is instantiated in dynamically interacting and spatially distributed functional brain networks remains an unexplored mystery in early Parkinson's disease (PD). We applied a machine-learning technology based on personalized sliding-window algorithm to track continuously time-varying and overlapping subnetworks under the functional brain networks calculated form resting state electroencephalogram data within a sample of 33 early PD patients (13 early PD patients with MCI and 20 early PD patients without MCI). We decoded a set of subnetworks that captured surprisingly dynamically varying and integrated interactions among certain brain lobes. We observed that the master expressed subnetworks were particularly transient, and flexibly switching between high and low expression during integration into a dynamic brain network. This transience was particularly salient in a subnetwork predominantly linking temporal-parietal-occipital lobes, which decreases in both expression and flexibility in early PD patients with MCI and expresses their degree of cognitive impairment. Moreover, MCI induced a regularly interrupted, slow evolution of subnetworks in functional brain network dynamics in early PD at the individual level, and the dynamic expression characteristics of subnetworks also reflected the degree of cognitive impairment in patients with early PD. Collectively, these results provide novel and deeper insights regarding MCI-induced abnormal dynamical interaction and large-scale changes in functional brain network of early PD.

Investigating cognitive neuroscience theories of human intelligence: A connectome-based predictive modeling approach

Central to modern neuroscientific theories of human intelligence is the notion that general intelligence depends on a primary brain region or network, engaging spatially localized (rather than global) neural representations. Recent findings in network neuroscience, however, challenge this assumption, providing evidence that general intelligence may depend on system-wide network mechanisms, suggesting that local representations are necessary but not sufficient to account for the neural architecture of human intelligence. Despite the importance of this key theoretical distinction, prior research has not systematically investigated the role of local versus global neural representations in predicting general intelligence. We conducted a large-scale connectome-based predictive modeling study (N = 297), administering resting-state fMRI and a comprehensive cognitive battery to evaluate the efficacy of modern neuroscientific theories of human intelligence, including spatially localized theories (Lateral Prefrontal Cortex Theory, Parieto-Frontal Integration Theory, and Multiple Demand Theory) and recent global accounts (Process Overlap Theory and Network Neuroscience Theory). The results of our study demonstrate that general intelligence can be predicted by local functional connectivity profiles but is most robustly explained by global profiles of whole-brain connectivity. Our findings further suggest that the improved efficacy of global theories is not reducible to a greater strength or number of connections, but instead results from considering both strong and weak connections that provide the basis for intelligence (as predicted by the Network Neuroscience Theory). Our results highlight the importance of considering local neural representations in the context of a global information-processing architecture, suggesting future directions for theory-driven research on system-wide network mechanisms underlying general intelligence.

An application of dynamical directed connectivity of ictal intracranial EEG recordings in seizure onset zone localization

BACKGROUND

Identification of the seizure onset zone (SOZ) is a challenging task in epilepsy surgery. Patients with epilepsy have an altered brain network, allowing connectivity-based analyses to have a great potential in SOZ identification. We investigated a dynamical directed connectivity analysis utilizing ictal intracranial electroencephalographic (iEEG) recordings and proposed an algorithm for SOZ identification based on grouping iEEG contacts.

NEW METHODS

Granger Causality was used for directed connectivity analysis in this study. The intracranial contacts were grouped into visually detected contacts (VDCs), which were identified as SOZ by epileptologists, and non-resected contacts (NRCs). The intragroup and intergroup directed connectivity for VDCs and NRCs were calculated around seizure onset. We then proposed an algorithm for SOZ identification based on the cross-correlation of intragroup outflow and inflow of SOZ candidate contacts.

RESULTS

Our results revealed that the intragroup connectivity of VDCs (VDC→VDC) was significantly larger than the intragroup connectivity of NRCs (NRC→NRC) and the intergroup connectivity between NRCs and VDCs (NRC→VDC) around seizure onset. We found that the proposed algorithm had 90.1 % accuracy for SOZ identification in the seizure-free patients.

COMPARISON WITH EXISTING METHODS

The existing connectivity-based methods for SOZ identification often use either outflow or inflow. In this study, SOZ contacts were identified by integrating outflow and inflow based on the cross correlation between these two measures.

CONCLUSIONS

The proposed group-based dynamical connectivity analysis in this study can aid our understanding of underlying seizure network and may be used to assist in identifying the SOZ contacts before epilepsy surgery.

Plasticity impairment alters community structure but permits successful pattern separation in a hippocampal network model

Patients who suffer from traumatic brain injury (TBI) often complain of learning and memory problems. Their symptoms are principally mediated by the hippocampus and the ability to adapt to stimulus, also known as neural plasticity. Therefore, one plausible injury mechanism is plasticity impairment, which currently lacks comprehensive investigation across TBI research. For these studies, we used a computational network model of the hippocampus that includes the dentate gyrus, CA3, and CA1 with neuron-scale resolution. We simulated mild injury through weakened spike-timing-dependent plasticity (STDP), which modulates synaptic weights according to causal spike timing. In preliminary work, we found functional deficits consisting of decreased firing rate and broadband power in areas CA3 and CA1 after STDP impairment. To address structural changes with these studies, we applied modularity analysis to evaluate how STDP impairment modifies community structure in the hippocampal network. We also studied the emergent function of network-based learning and found that impaired networks could acquire conditioned responses after training, but the magnitude of the response was significantly lower. Furthermore, we examined pattern separation, a prerequisite of learning, by entraining two overlapping patterns. Contrary to our initial hypothesis, impaired networks did not exhibit deficits in pattern separation with either population- or rate-based coding. Collectively, these results demonstrate how a mechanism of injury that operates at the synapse regulates circuit function.

Evaluating the reliability, validity, and utility of overlapping networks: Implications for network theories of cognition

Brain network definitions typically assume nonoverlap or minimal overlap, ignoring regions' connections to multiple networks. However, new methods are emerging that emphasize network overlap. Here, we investigated the reliability and validity of one assignment method, the mixed membership algorithm, and explored its potential utility for identifying gaps in existing network models of cognition. We first assessed between-sample reliability of overlapping assignments with a split-half design; a bootstrapped Dice similarity analysis demonstrated good agreement between the networks from the two subgroups. Next, we assessed whether overlapping networks captured expected nonoverlapping topographies; overlapping networks captured portions of one to three nonoverlapping topographies, which aligned with canonical network definitions. Following this, a relative entropy analysis showed that a majority of regions participated in more than one network, as is seen biologically, and many regions did not show preferential connection to any one network. Finally, we explored overlapping network membership in regions of the dual-networks model of cognitive control, showing that almost every region was a member of multiple networks. Thus, the mixed membership algorithm produces consistent and biologically plausible networks, which presumably will allow for the development of more complete network models of cognition.

A spatio-temporal decomposition framework for dynamic functional connectivity in the human brain

Much recent attention has been directed toward investigating the spatial and temporal organization of brain dynamics, but the rules which constrain the variation of spatio-temporal organization in functional connectivity under different brain states remain unclear. Here, we developed a novel computational approach based on tensor decomposition and regularization to represent dynamic functional connectivity as a linear combination of dynamic modules and time-varying weights. In this approach, dynamic modules represent co-activating functional connectivity patterns, and time-varying weights represent the temporal expression of dynamic modules. We applied this dynamic decomposition model (DDM) on a resting-state fMRI dataset and found that whole-brain dynamic functional connectivity can be decomposed as a linear combination of eight dynamic modules which we summarize as 'high order modules' and 'primary-high order modules', according to their spatial attributes and correspondence with existing intrinsic functional brain networks. By clustering the time-varying weights, we identified five brain states including three major states and two minor states. We found that state transitions mainly occurred between the three major states, and that temporal variation of dynamic modules may contribute to brain state transitions. We then conceptualized the variability of weights as the flexibility of the corresponding dynamic modules and found that different dynamic modules exhibit different amounts of flexibility and contribute to different cognitive measures. Finally, we applied DDM to a schizophrenia resting-state fMRI dataset and found that atypical flexibility of dynamic modules correlates with impaired cognitive flexibility in schizophrenia. Overall, this work provides a quantitative framework that characterizes temporal variation in the topology of dynamic functional connectivity.

Intra- and inter-regional dynamics in cortical-striatal-tegmental networks

It is increasingly recognized that networks of brain areas work together to accomplish computational goals. However, functional connectivity networks are not often compared between different behavioral states and across different frequencies of electrical oscillatory signals. In addition, connectivity is always defined as the strength of signal relatedness between two atlas-based anatomical locations. Here, we performed an exploratory analysis using data collectected from high density arrays in the prefrontal cortex (PFC), striatum (STR), and ventral tegmental area (VTA) of male rats. These areas have all been implicated in a wide range of different tasks and computations including various types of memory as well as reward valuation, habit formation and execution, and skill learning. Novel intra-regional clustering analyses identified patterns of spatially restricted, temporally coherent, and frequency specific signals that were reproducible across days and were modulated by behavioral states. Multiple clusters were identified within each anatomical region, indicating a mesoscopic scale of organization. Generalized eigendecomposition (GED) was used to dimension-reduce each cluster to a single component time series. Dense inter-cluster connectivity was modulated by behavioral state, with connectivity becoming reduced when the animals were exposed to a novel object, compared to a baseline condition. Behavior-modulated connectivity changes were seen across the spectrum, with delta, theta, and gamma all being modulated. These results demonstrate the brain's ability to reorganize functionally at both the intra- and inter-regional levels during different behavioral states.

Long-term brain network reorganization predicts responsive neurostimulation outcomes for focal epilepsy

Responsive neurostimulation (RNS) devices, able to detect imminent seizures and to rapidly deliver electrical stimulation to the brain, are effective in reducing seizures in some patients with focal epilepsy. However, therapeutic response to RNS is often slow, is highly variable, and defies prognostication based on clinical factors. A prevailing view holds that RNS efficacy is primarily mediated by acute seizure termination; yet, stimulations greatly outnumber seizures and occur mostly in the interictal state, suggesting chronic modulation of brain networks that generate seizures. Here, using years-long intracranial neural recordings collected during RNS therapy, we found that patients with the greatest therapeutic benefit undergo progressive, frequency-dependent reorganization of interictal functional connectivity. The extent of this reorganization scales directly with seizure reduction and emerges within the first year of RNS treatment, enabling potential early prediction of therapeutic response. Our findings reveal a mechanism for RNS that involves network plasticity and may inform development of next-generation devices for epilepsy.

Nodal Memberships to Communities of Functional Brain Networks Reveal Functional Flexibility and Individualized Connectome

Human brain network is organized as interconnected communities for supporting cognition and behavior. Despite studies on the nonoverlapping communities of brain network, overlapping community structure and its relationship to brain function remain largely unknown. With this consideration, we employed the Bayesian nonnegative matrix factorization to decompose the functional brain networks constructed from resting-state fMRI data into overlapping communities with interdigitated mapping to functional subnetworks. By examining the heterogeneous nodal membership to communities, we classified nodes into three classes: Most nodes in somatomotor and limbic subnetworks were affiliated with one dominant community and classified as unimodule nodes; most nodes in attention and frontoparietal subnetworks were affiliated with more than two communities and classified as multimodule nodes; and the remaining nodes affiliated with two communities were classified as bimodule nodes. This three-class paradigm was highly reproducible across sessions and subjects. Furthermore, the more likely a node was classified as multimodule node, the more flexible it will be engaged in multiple tasks. Finally, the FC feature vector associated with multimodule nodes could serve as connectome "fingerprinting" to gain high subject discriminability. Together, our findings offer new insights on the flexible spatial overlapping communities that related to task-based functional flexibility and individual connectome "fingerprinting."

A single mode of population covariation associates brain networks structure and behavior and predicts individual subjects' age

Multiple human behaviors improve early in life, peaking in young adulthood, and declining thereafter. Several properties of brain structure and function progress similarly across the lifespan. Cognitive and neuroscience research has approached aging primarily using associations between a few behaviors, brain functions, and structures. Because of this, the multivariate, global factors relating brain and behavior across the lifespan are not well understood. We investigated the global patterns of associations between 334 behavioral and clinical measures and 376 brain structural connections in 594 individuals across the lifespan. A single-axis associated changes in multiple behavioral domains and brain structural connections (r = 0.5808). Individual variability within the single association axis well predicted the age of the subject (r = 0.6275). Representational similarity analysis evidenced global patterns of interactions across multiple brain network systems and behavioral domains. Results show that global processes of human aging can be well captured by a multivariate data fusion approach.

The modulation of brain network integration and arousal during exploration

There is growing interest in how neuromodulators shape brain networks. Recent neuroimaging studies provide evidence that brainstem arousal systems, such as the locus coeruleus-norepinephrine system (LC-NE), influence functional connectivity and brain network topology, suggesting they have a role in flexibly reconfiguring brain networks in order to adapt behavior and cognition to environmental demands. To date, however, the relationship between brainstem arousal systems and functional connectivity has not been assessed within the context of a task with an established relationship between arousal and behavior, with most prior studies relying on incidental variations in arousal or pharmacological manipulation and static brain networks constructed over long periods of time. These factors have likely contributed to a heterogeneity of effects across studies. To address these issues, we took advantage of the association between LC-NE-linked arousal and exploration to probe the relationships between exploratory choice, arousal—as measured indirectly via pupil diameter—and brain network dynamics. Exploration in a bandit task was associated with a shift toward fewer, more weakly connected modules that were more segregated in terms of connectivity and topology but more integrated with respect to the diversity of cognitive systems represented in each module. Functional connectivity strength decreased, and changes in connectivity were correlated with changes in pupil diameter, in line with the hypothesis that brainstem arousal systems influence the dynamic reorganization of brain networks. More broadly, we argue that carefully aligning dynamic network analyses with task designs can increase the temporal resolution at which behaviorally- and cognitively-relevant modulations can be identified, and offer these results as a proof of concept of this approach.

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.

Measurement reliability for individual differences in multilayer network dynamics: Cautions and considerations

Multilayer network models have been proposed as an effective means of capturing the dynamic configuration of distributed neural circuits and quantitatively describing how communities vary over time. Beyond general insights into brain function, a growing number of studies have begun to employ these methods for the study of individual differences. However, test-retest reliabilities for multilayer network measures have yet to be fully quantified or optimized, potentially limiting their utility for individual difference studies. Here, we systematically evaluated the impact of multilayer community detection algorithms, selection of network parameters, scan duration, and task condition on test-retest reliabilities of multilayer network measures (i.e., flexibility, integration, and recruitment). A key finding was that the default method used for community detection by the popular generalized Louvain algorithm can generate erroneous results. Although available, an updated algorithm addressing this issue is yet to be broadly adopted in the neuroimaging literature. Beyond the algorithm, the present work identified parameter selection as a key determinant of test-retest reliability; however, optimization of these parameters and expected reliabilities appeared to be dataset-specific. Once parameters were optimized, consistent with findings from the static functional connectivity literature, scan duration was a much stronger determinant of reliability than scan condition. When the parameters were optimized and scan duration was sufficient, both passive (i.e., resting state, Inscapes, and movie) and active (i.e., flanker) tasks were reliable, although reliability in the movie watching condition was significantly higher than in the other three tasks. The minimal data requirement for achieving reliable measures for the movie watching condition was 20 min, and 30 min for the other three tasks. Our results caution the field against the use of default parameters without optimization based on the specific datasets to be employed - a process likely to be limited for most due to the lack of test-retest samples to enable parameter optimization.

Transcranial direct current stimulation modulates brain functional connectivity in autism

Autism spectrum disorder (ASD) is characterized by deficits in social interactions, impairments in language and communication, and highly restricted behavioral interests. Transcranial direct current stimulation (tDCS) is a widely used form of noninvasive stimulation and may have therapeutic potential for ASD. So far, despite the widespread use of this technique in the neuroscience field, its effects on network-level neural activity and the underlying mechanisms of any effects are still unclear. In the present study, we used electroencephalography (EEG) to investigate tDCS induced brain network changes in children with ASD before and after active and sham stimulation. We recorded 5 min of resting state EEG before and after a single session of tDCS (of approximately 20 min) over dorsolateral prefrontal cortex (DLPFC). Two network-based methods were applied to investigate tDCS modulation on brain networks: 1) temporal network dynamics were analyzed by comparing "flexibility" changes before vs after stimulation, and 2) frequency specific network changes were identified using non-negative matrix factorization (NMF). We found 1) an increase in network flexibility following tDCS (rapid network configuration of dynamic network communities), 2) specific increase in interhemispheric connectivity within the alpha frequency band following tDCS. Together, these results demonstrate that tDCS could help modify both local and global brain network dynamics, and highlight stimulation-induced differences in the manifestation of network reconfiguration. Meanwhile, frequency-specific subnetworks, as a way to index local and global information processing, highlight the core modulatory effects of tDCS on the modular architecture of the functional connectivity patterns within higher frequency bands.

Roles of the prefrontal cortex in learning to time the onset of pre-existing motor programs

The prefrontal cortex (PFC) is involved in cognitive control of motor activities and timing of future intensions. This study investigated the cognitive control of balance recovery in response to unpredictable gait perturbations and the role of PFC subregions in learning by repetition. Bilateral dorsolateral (DLPFC), ventrolateral (VLPFC), frontopolar (FPFC) and orbitofrontal (OFC) cortex hemodynamic changes induced by unpredictable slips were analyzed as a function of successive trials in ten healthy young adults. Slips were induced by the acceleration of one belt as the participant walked on a split-belt treadmill. A portable functional near-infrared spectroscope monitored PFC activities quantified by oxyhemoglobin (ΔO₂Hb) and deoxyhemoglobin (ΔHbR) during the consecutive trial phases: standing, walking, slip-recovery. During the first 3 trials, the average oxyhemoglobin (ΔO₂Hb_avg) in the DLPFC, VLPFC, FPFC, and OFC cortex was significantly higher during slip-recovery than unperturbed walking or the standing baseline. Then, ΔO₂Hb_avg decreased progressively from trial-to-trial in the DLPFC, VLPFC, and FPFC, but increased and then remained constant in the OFC. The average deoxyhemoglobin (ΔHbR_avg) presented mirror patterns. These changes after the third trial were paralleled by the progressive improvement of recovery revealed by kinematic variables. The results corroborate our previous hypothesis that only timing of the onset of a “good enough recovery motor program” is learned with practice. They also strongly support the assumption that the PFC contributes to the recall of pre-existing motor programs whose onset timing is adjusted by the OFC. Hence, learning is clearly divided into two steps delineated by the switch in activity of the OFC. Additionally, motor processes appear to share the working memory as well as decisional and predictive resources of the cognitive system.

Learning in brain-computer interface control evidenced by joint decomposition of brain and behavior

OBJECTIVE

Motor imagery-based brain-computer interfaces (BCIs) use an individual's ability to volitionally modulate localized brain activity, often as a therapy for motor dysfunction or to probe causal relations between brain activity and behavior. However, many individuals cannot learn to successfully modulate their brain activity, greatly limiting the efficacy of BCI for therapy and for basic scientific inquiry. Formal experiments designed to probe the nature of BCI learning have offered initial evidence that coherent activity across spatially distributed and functionally diverse cognitive systems is a hallmark of individuals who can successfully learn to control the BCI. However, little is known about how these distributed networks interact through time to support learning.

APPROACH

Here, we address this gap in knowledge by constructing and applying a multimodal network approach to decipher brain-behavior relations in motor imagery-based brain-computer interface learning using magnetoencephalography. Specifically, we employ a minimally constrained matrix decomposition method - non-negative matrix factorization - to simultaneously identify regularized, covarying subgraphs of functional connectivity, to assess their similarity to task performance, and to detect their time-varying expression.

MAIN RESULTS

We find that learning is marked by diffuse brain-behavior relations: good learners displayed many subgraphs whose temporal expression tracked performance. Individuals also displayed marked variation in the spatial properties of subgraphs such as the connectivity between the frontal lobe and the rest of the brain, and in the temporal properties of subgraphs such as the stage of learning at which they reached maximum expression. From these observations, we posit a conceptual model in which certain subgraphs support learning by modulating brain activity in sensors near regions important for sustaining attention. To test this model, we use tools that stipulate regional dynamics on a networked system (network control theory), and find that good learners display a single subgraph whose temporal expression tracked performance and whose architecture supports easy modulation of sensors located near brain regions important for attention.

SIGNIFICANCE

The nature of our contribution to the neuroscience of BCI learning is therefore both computational and theoretical; we first use a minimally-constrained, individual specific method of identifying mesoscale structure in dynamic brain activity to show how global connectivity and interactions between distributed networks supports BCI learning, and then we use a formal network model of control to lend theoretical support to the hypothesis that these identified subgraphs are well suited to modulate attention.

Functional brain network reconfiguration during learning in a dynamic environment

When learning about dynamic and uncertain environments, people should update their beliefs most strongly when new evidence is most informative, such as when the environment undergoes a surprising change or existing beliefs are highly uncertain. Here we show that modulations of surprise and uncertainty are encoded in a particular, temporally dynamic pattern of whole-brain functional connectivity, and this encoding is enhanced in individuals that adapt their learning dynamics more appropriately in response to these factors. The key feature of this whole-brain pattern of functional connectivity is stronger connectivity, or functional integration, between the fronto-parietal and other functional systems. Our results provide new insights regarding the association between dynamic adjustments in learning and dynamic, large-scale changes in functional connectivity across the brain.

Early childhood developmental functional connectivity of autistic brains with non-negative matrix factorization

Autism spectrum disorder (ASD) is associated with altered patterns of over- and under-connectivity of neural circuits. Age-related changes in neural connectivities remain unclear for autistic children as compared with normal children. In this study, a parts-based network-decomposition technique, known as non-negative matrix factorization (NMF), was applied to identify a set of possible abnormal connectivity patterns in brains affected by ASD, using resting-state electroencephalographic (EEG) data. Age-related changes in connectivities in both inter- and intra-hemispheric areas were studied in a total of 256 children (3-6 years old), both with and without ASD. The findings included the following: (1) the brains of children affected by ASD were characterized by a general trend toward long-range under-connectivity, particularly in interhemispheric connections, combined with short-range over-connectivity; (2) long-range connections were often associated with slower rhythms (δ and θ), whereas synchronization of short-range networks tended to be associated with faster frequencies (α and β); and (3) the α-band specific patterns of interhemispheric connections in ASD could be the most prominent during early childhood neurodevelopment. Therefore, NMF would be useful for further exploring the early childhood developmental functional connectivity of children aged 3-6 with ASD as well as with typical development. Additionally, long-range interhemispheric alterations in connectivity may represent a potential biomarker for the identification of ASD.

LRRK2 and GBA Variants Exert Distinct Influences on Parkinson's Disease-Specific Metabolic Networks

The natural history of idiopathic Parkinson's disease (PD) varies considerably across patients. While PD is generally sporadic, there are known genetic influences: the two most common, mutations in the LRRK2 or GBA1 gene, are associated with slower and more aggressive progression, respectively. Here, we applied graph theory to metabolic brain imaging to understand the effects of genotype on the organization of previously established PD-specific networks. We found that closely matched PD patient groups with the LRRK2-G2019S mutation (PD-LRRK2) or GBA1 variants (PD-GBA) expressed the same disease networks as sporadic disease (sPD), but PD-LRRK2 and PD-GBA patients exhibited abnormal increases in network connectivity that were not present in sPD. Using a community detection strategy, we found that the location and modular distribution of these connections differed strikingly across genotypes. In PD-LRRK2, connections were gained within the network core, with the formation of distinct functional pathways linking the cerebellum and putamen. In PD-GBA, by contrast, the majority of functional connections were formed outside the core, involving corticocortical pathways at the network periphery. Strategically localized connections within the core in PD-LRRK2 may maintain PD network activity at lower levels than in PD-GBA, resulting in a less aggressive clinical course.

Linking personality types to depressive symptoms: A prospective typology based on neuroticism, extraversion and conscientiousness

From a traditional variable-centered perspective, the personality traits specifically linked to depressive symptoms are neuroticism, extraversion, and conscientiousness (NEC). Few studies have considered the interactive effects of personality traits within individuals from the taxonomic perspective. We explored novel NEC personality subtypes in general and subthreshold depressive subjects by using graph-based classification algorithms and multiple external validators. Personality and depressive symptoms were assessed in 1055 healthy subjects (150 with subthreshold depression) using the NEO-PI-R and BDI at baseline. NEC personality data were used to identify subtypes in healthy and subthreshold depressive samples, and external validators, including current and longitudinal depressive symptoms and seven subcortical gray matter volumes, were examined. Four novel NEC personality types in the general sample showed superior current and longitudinal behavioral validation of depressive symptoms as well as some discrimination in biological indicators (putamen, nucleus accumbens, and amygdala). Four profiles defined in subthreshold depression selectively exhibited meaningful differences in longitudinal depressive symptoms. In both samples, some types adhere to the principles previously described NEC three-way interaction. The resulting typology, especially the four types in the general population, linked depressive symptoms provided a superior description of within-person organization of neuroticism, extraversion, and conscientiousness.

A Minimum Free Energy Model of Motor Learning

Even highly trained behaviors demonstrate variability, which is correlated with performance on current and future tasks. An objective of motor learning that is general enough to explain these phenomena has not been precisely formulated. In this six-week longitudinal learning study, participants practiced a set of motor sequences each day, and neuroimaging data were collected on days 1, 14, 28, and 42 to capture the neural correlates of the learning process. In our analysis, we first modeled the underlying neural and behavioral dynamics during learning. Our results demonstrate that the densities of whole-brain response, task-active regional response, and behavioral performance evolve according to a Fokker-Planck equation during the acquisition of a motor skill. We show that this implies that the brain concurrently optimizes the entropy of a joint density over neural response and behavior (as measured by sampling over multiple trials and subjects) and the expected performance under this density; we call this formulation of learning minimum free energy learning (MFEL). This model provides an explanation as to how behavioral variability can be tuned while simultaneously improving performance during learning. We then develop a novel variant of inverse reinforcement learning to retrieve the cost function optimized by the brain during the learning process, as well as the parameter used to tune variability. We show that this population-level analysis can be used to derive a learning objective that each subject optimizes during his or her study. In this way, MFEL effectively acts as a unifying principle, allowing users to precisely formulate learning objectives and infer their structure.

Memory systems 2018 - Towards a new paradigm

The multiple memory systems theory (MMS) postulates that the brain stores information based on the independent and parallel activity of a number of modules, each with distinct properties, dynamics, and neural basis. Much of the evidence for this theory comes from dissociation studies indicating that damage to restricted brain areas cause selective types of memory deficits. MMS has been the prevalent paradigm in memory research for more than thirty years, even as it has been adjusted several times to accommodate new data. However, recent empirical results indicating that the memory systems are not always dissociable constitute a challenge to fundamental tenets of the current theory because they suggest that representations formed by individual memory systems can contribute to more than one type of memory-driven behavioral strategy. This problem can be addressed by applying a dynamic network perspective to memory architecture. According to this view, memory networks can reconfigure or transiently couple in response to environmental demands. Within this context, the neural network underlying a specific memory system can act as an independent unit or as an integrated component of a higher order meta-network. This dynamic network model proposes a way in which empirical evidence that challenges the idea of distinct memory systems can be incorporated within a modular memory architecture. The model also provides a framework to account for the complex interactions among memory systems demonstrated at the behavioral level. Advances in the study of dynamic networks can generate new ideas to experimentally manipulate and control memory in basic or clinical research.

Modeling and interpreting mesoscale network dynamics

Recent advances in brain imaging techniques, measurement approaches, and storage capacities have provided an unprecedented supply of high temporal resolution neural data. These data present a remarkable opportunity to gain a mechanistic understanding not just of circuit structure, but also of circuit dynamics, and its role in cognition and disease. Such understanding necessitates a description of the raw observations, and a delineation of computational models and mathematical theories that accurately capture fundamental principles behind the observations. Here we review recent advances in a range of modeling approaches that embrace the temporally-evolving interconnected structure of the brain and summarize that structure in a dynamic graph. We describe recent efforts to model dynamic patterns of connectivity, dynamic patterns of activity, and patterns of activity atop connectivity. In the context of these models, we review important considerations in statistical testing, including parametric and non-parametric approaches. Finally, we offer thoughts on careful and accurate interpretation of dynamic graph architecture, and outline important future directions for method development.

Network Approaches to Understand Individual Differences in Brain Connectivity: Opportunities for Personality Neuroscience

Over the past decade, advances in the interdisciplinary field of network science have provided a framework for understanding the intrinsic structure and function of human brain networks. A particularly fruitful area of this work has focused on patterns of functional connectivity derived from non-invasive neuroimaging techniques such as functional magnetic resonance imaging (fMRI). An important subset of these efforts has bridged the computational approaches of network science with the rich empirical data and biological hypotheses of neuroscience, and this research has begun to identify features of brain networks that explain individual differences in social, emotional, and cognitive functioning. The most common approach estimates connections assuming a single configuration of edges that is stable across the experimental session. In the literature, this is referred to as a static network approach, and researchers measure static brain networks while a subject is either at rest or performing a cognitively demanding task. Research on social and emotional functioning has primarily focused on linking static brain networks with individual differences, but recent advances have extended this work to examine temporal fluctuations in dynamic brain networks. Mounting evidence suggests that both the strength and flexibility of time-evolving brain networks influence individual differences in executive function, attention, working memory, and learning. In this review, we first examine the current evidence for brain networks involved in cognitive functioning. Then we review some preliminary evidence linking static network properties to individual differences in social and emotional functioning. We then discuss the applicability of emerging dynamic network methods for examining individual differences in social and emotional functioning. We close with an outline of important frontiers at the intersection between network science and neuroscience that will enhance our understanding of the neurobiological underpinnings of social behavior.

Subgraphs of functional brain networks identify dynamical constraints of cognitive control

Brain anatomy and physiology support the human ability to navigate a complex space of perceptions and actions. To maneuver across an ever-changing landscape of mental states, the brain invokes cognitive control-a set of dynamic processes that engage and disengage different groups of brain regions to modulate attention, switch between tasks, and inhibit prepotent responses. Current theory posits that correlated and anticorrelated brain activity may signify cooperative and competitive interactions between brain areas that subserve adaptive behavior. In this study, we use a quantitative approach to identify distinct topological motifs of functional interactions and examine how their expression relates to cognitive control processes and behavior. In particular, we acquire fMRI BOLD signal in twenty-eight healthy subjects as they perform two cognitive control tasks-a Stroop interference task and a local-global perception switching task using Navon figures-each with low and high cognitive control demand conditions. Based on these data, we construct dynamic functional brain networks and use a parts-based, network decomposition technique called non-negative matrix factorization to identify putative cognitive control subgraphs whose temporal expression captures distributed network structures involved in different phases of cooperative and competitive control processes. Our results demonstrate that temporal expression of the subgraphs fluctuate alongside changes in cognitive demand and are associated with individual differences in task performance. These findings offer insight into how coordinated changes in the cooperative and competitive roles of cognitive systems map trajectories between cognitively demanding brain states.