Endogenous electric fields may guide neocortical network activity

Local field potentials and the underlying endogenous electric fields (EFs) are traditionally considered to be epiphenomena of structured neuronal network activity. Recently, however, externally applied EFs have been shown to modulate pharmacologically evoked network activity in rodent hippocampus. In contrast, very little is known about the role of endogenous EFs during physiological activity states in neocortex. Here, we used the neocortical slow oscillation in vitro as a model system to show that weak sinusoidal and naturalistic EFs enhance and entrain physiological neocortical network activity with an amplitude threshold within the range of in vivo endogenous field strengths. Modulation of network activity by positive and negative feedback fields based on the network activity in real-time provide direct evidence for a feedback loop between neuronal activity and endogenous EF. This significant susceptibility of active networks to EFs that only cause small changes in membrane potential in individual neurons suggests that endogenous EFs could guide neocortical network activity.

Cartography of complex networks: modules and universal roles

Integrative approaches to the study of complex systems demand that one knows the manner in which the parts comprising the system are connected. The structure of the complex network defining the interactions provides insight into the function and evolution of the components of the system. Unfortunately, the large size and intricacy of these networks implies that such insight is usually difficult to extract. Here, we propose a method that allows one to systematically extract and display information contained in complex networks. Specifically, we demonstrate that one can (i) find modules in complex networks and (ii) classify nodes into universal roles according to their pattern of within- and between-module connections. The method thus yields a 'cartographic representation' of complex networks.

Transcranial magnetic stimulation modulates the brain's intrinsic activity in a frequency-dependent manner

Intrinsic activity in the brain is organized into networks. Although constrained by their anatomical connections, functional correlations between nodes of these networks reorganize dynamically. Dynamic organization implies that couplings between network nodes can be reconfigured to support processing demands. To explore such reconfigurations, we combined repetitive transcranial magnetic stimulation (rTMS) and functional connectivity MRI (fcMRI) to modulate cortical activity in one node of the default network, and assessed the effect of this upon functional correlations throughout the network. Two different frequencies of rTMS to the same default network node (the left posterior inferior parietal lobule, lpIPL) induced two topographically distinct changes in functional connectivity. High-frequency rTMS to lpIPL decreased functional correlations between cortical default network nodes, but not between these nodes and the hippocampal formation. In contrast, low frequency rTMS to lpIPL did not alter connectivity between cortical default network nodes, but increased functional correlations between lpIPL and the hippocampal formation. These results suggest that the default network is composed of (at least) two subsystems. More broadly, the finding that two rTMS stimulation regimens to the same default network node have distinct effects reveals that this node is embedded within a network that possesses multiple, functionally distinct relationships among its distributed partners.

Toward the dynamic interactome: it's about time

Dynamic molecular interactions play a central role in regulating the functioning of cells and organisms. The availability of experimentally determined large-scale cellular networks, along with other high-throughput experimental data sets that provide snapshots of biological systems at different times and conditions, is increasingly helpful in elucidating interaction dynamics. Here we review the beginnings of a new subfield within computational biology, one focused on the global inference and analysis of the dynamic interactome. This burgeoning research area, which entails a shift from static to dynamic network analysis, promises to be a major step forward in our ability to model and reason about cellular function and behavior.

Decomposition of spontaneous brain activity into distinct fMRI co-activation patterns

Recent fMRI studies have shown that analysis of the human brain's spontaneous activity may provide a powerful approach to reveal its functional organization. Dedicated methods have been proposed to investigate co-variation of signals from different brain regions, with the goal of revealing neuronal networks (NNs) that may serve specialized functions. However, these analysis methods generally do not take into account a potential non-stationary (variable) interaction between brain regions, and as a result have limited effectiveness. To address this, we propose a novel analysis method that uses clustering analysis to sort and selectively average fMRI activity time frames to produce a set of co-activation patterns. Compared to the established networks extracted with conventional analysis methods, these co-activation patterns demonstrate novel network features with apparent relevance to the brain's functional organization.

Quantifying environmental adaptation of metabolic pathways in metagenomics

Recently, approaches have been developed to sample the genetic content of heterogeneous environments (metagenomics). However, by what means these sequences link distinct environmental conditions with specific biological processes is not well understood. Thus, a major challenge is how the usage of particular pathways and subnetworks reflects the adaptation of microbial communities across environments and habitats-i.e., how network dynamics relates to environmental features. Previous research has treated environments as discrete, somewhat simplified classes (e.g., terrestrial vs. marine), and searched for obvious metabolic differences among them (i.e., treating the analysis as a typical classification problem). However, environmental differences result from combinations of many factors, which often vary only slightly. Therefore, we introduce an approach that employs correlation and regression to relate multiple, continuously varying factors defining an environment to the extent of particular microbial pathways present in a geographic site. Moreover, rather than looking only at individual correlations (one-to-one), we adapted canonical correlation analysis and related techniques to define an ensemble of weighted pathways that maximally covaries with a combination of environmental variables (many-to-many), which we term a metabolic footprint. Applied to available aquatic datasets, we identified footprints predictive of their environment that can potentially be used as biosensors. For example, we show a strong multivariate correlation between the energy-conversion strategies of a community and multiple environmental gradients (e.g., temperature). Moreover, we identified covariation in amino acid transport and cofactor synthesis, suggesting that limiting amounts of cofactor can (partially) explain increased import of amino acids in nutrient-limited conditions.

Neural correlates of tactile detection: a combined magnetoencephalography and biophysically based computational modeling study

Previous reports conflict as to the role of primary somatosensory neocortex (SI) in tactile detection. We addressed this question in normal human subjects using whole-head magnetoencephalography (MEG) recording. We found that the evoked signal (0-175 ms) showed a prominent equivalent current dipole that localized to the anterior bank of the postcentral gyrus, area 3b of SI. The magnitude and timing of peaks in the SI waveform were stimulus amplitude dependent and predicted perception beginning at approximately 70 ms after stimulus. To make a direct and principled connection between the SI waveform and underlying neural dynamics, we developed a biophysically realistic computational SI model that contained excitatory and inhibitory neurons in supragranular and infragranular layers. The SI evoked response was successfully reproduced from the intracellular currents in pyramidal neurons driven by a sequence of lamina-specific excitatory input, consisting of output from the granular layer (approximately 25 ms), exogenous input to the supragranular layers (approximately 70 ms), and a second wave of granular output (approximately 135 ms). The model also predicted that SI correlates of perception reflect stronger and shorter-latency supragranular and late granular drive during perceived trials. These findings strongly support the view that signatures of tactile detection are present in human SI and are mediated by local neural dynamics induced by lamina-specific synaptic drive. Furthermore, our model provides a biophysically realistic solution to the MEG signal and can predict the electrophysiological correlates of human perception.

Exploring the roles of noise in the eukaryotic cell cycle

The DNA replication-division cycle of eukaryotic cells is controlled by a complex network of regulatory proteins, called cyclin-dependent kinases, and their activators and inhibitors. Although comprehensive and accurate deterministic models of the control system are available for yeast cells, reliable stochastic simulations have not been carried out because the full reaction network has yet to be expressed in terms of elementary reaction steps. As a first step in this direction, we present a simplified version of the control system that is suitable for exact stochastic simulation of intrinsic noise caused by molecular fluctuations and extrinsic noise because of unequal division. The model is consistent with many characteristic features of noisy cell cycle progression in yeast populations, including the observation that mRNAs are present in very low abundance (approximately 1 mRNA molecule per cell for each expressed gene). For the control system to operate reliably at such low mRNA levels, some specific mRNAs in our model must have very short half-lives (<1 min). If these mRNA molecules are longer-lived (perhaps 2 min), then the intrinsic noise in our simulations is too large, and there must be some additional noise suppression mechanisms at work in cells.

Progression from a stem cell-like state to early differentiation in the C. elegans germ line

Controls of stem cell maintenance and early differentiation are known in several systems. However, the progression from stem cell self-renewal to overt signs of early differentiation is a poorly understood but important problem in stem cell biology. The Caenorhabditis elegans germ line provides a genetically defined model for studying that progression. In this system, a single-celled mesenchymal niche, the distal tip cell (DTC), employs GLP-1/Notch signaling and an RNA regulatory network to balance self-renewal and early differentiation within the "mitotic region," which continuously self-renews while generating new gametes. Here, we investigate germ cells in the mitotic region for their capacity to differentiate and their state of maturation. Two distinct pools emerge. The "distal pool" is maintained by the DTC in an essentially uniform and immature or "stem cell-like" state; the "proximal pool," by contrast, contains cells that are maturing toward early differentiation and are likely transit-amplifying cells. A rough estimate of pool sizes is 30-70 germ cells in the distal immature pool and approximately 150 in the proximal transit-amplifying pool. We present a simple model for how the network underlying the switch between self-renewal and early differentiation may be acting in these two pools. According to our model, the self-renewal mode of the network maintains the distal pool in an immature state, whereas the transition between self-renewal and early differentiation modes of the network underlies the graded maturation of germ cells in the proximal pool. We discuss implications of this model for controls of stem cells more broadly.

Dynamics of multistable states during ongoing and evoked cortical activity

Single-trial analyses of ensemble activity in alert animals demonstrate that cortical circuits dynamics evolve through temporal sequences of metastable states. Metastability has been studied for its potential role in sensory coding, memory, and decision-making. Yet, very little is known about the network mechanisms responsible for its genesis. It is often assumed that the onset of state sequences is triggered by an external stimulus. Here we show that state sequences can be observed also in the absence of overt sensory stimulation. Analysis of multielectrode recordings from the gustatory cortex of alert rats revealed ongoing sequences of states, where single neurons spontaneously attain several firing rates across different states. This single-neuron multistability represents a challenge to existing spiking network models, where typically each neuron is at most bistable. We present a recurrent spiking network model that accounts for both the spontaneous generation of state sequences and the multistability in single-neuron firing rates. Each state results from the activation of neural clusters with potentiated intracluster connections, with the firing rate in each cluster depending on the number of active clusters. Simulations show that the model's ensemble activity hops among the different states, reproducing the ongoing dynamics observed in the data. When probed with external stimuli, the model predicts the quenching of single-neuron multistability into bistability and the reduction of trial-by-trial variability. Both predictions were confirmed in the data. Together, these results provide a theoretical framework that captures both ongoing and evoked network dynamics in a single mechanistic model.

The Neuroanatomical Ultrastructure and Function of a Biological Ring Attractor

Neural representations of head direction (HD) have been discovered in many species. Theoretical work has proposed that the dynamics associated with these representations are generated, maintained, and updated by recurrent network structures called ring attractors. We evaluated this theorized structure-function relationship by performing electron-microscopy-based circuit reconstruction and RNA profiling of identified cell types in the HD system of Drosophila melanogaster. We identified motifs that have been hypothesized to maintain the HD representation in darkness, update it when the animal turns, and tether it to visual cues. Functional studies provided support for the proposed roles of individual excitatory or inhibitory circuit elements in shaping activity. We also discovered recurrent connections between neuronal arbors with mixed pre- and postsynaptic specializations. Our results confirm that the Drosophila HD network contains the core components of a ring attractor while also revealing unpredicted structural features that might enhance the network's computational power.

Topological effects of network structure on long-term social network dynamics in a wild mammal

Social structure influences ecological processes such as dispersal and invasion, and affects survival and reproductive success. Recent studies have used static snapshots of social networks, thus neglecting their temporal dynamics, and focused primarily on a limited number of variables that might be affecting social structure. Here, instead we modelled effects of multiple predictors of social network dynamics in the spotted hyena, using observational data collected during 20 years of continuous field research in Kenya. We tested the hypothesis that the current state of the social network affects its long-term dynamics. We employed stochastic agent-based models that allowed us to estimate the contribution of multiple factors to network changes. After controlling for environmental and individual effects, we found that network density and individual centrality affected network dynamics, but that social bond transitivity consistently had the strongest effects. Our results emphasise the significance of structural properties of networks in shaping social dynamics.

Integration or predictability? A further specification of the functional role of gamma oscillations in language comprehension

Gamma-band neuronal synchronization during sentence-level language comprehension has previously been linked with semantic unification. Here, we attempt to further narrow down the functional significance of gamma during language comprehension, by distinguishing between two aspects of semantic unification: successful integration of word meaning into the sentence context, and prediction of upcoming words. We computed event-related potentials (ERPs) and frequency band-specific electroencephalographic (EEG) power changes while participants read sentences that contained a critical word (CW) that was (1) both semantically congruent and predictable (high cloze, HC), (2) semantically congruent but unpredictable (low cloze, LC), or (3) semantically incongruent (and therefore also unpredictable; semantic violation, SV). The ERP analysis showed the expected parametric N400 modulation (HC < LC < SV). The time-frequency analysis showed qualitatively different results. In the gamma-frequency range, we observed a power increase in response to the CW in the HC condition, but not in the LC and the SV conditions. Additionally, in the theta frequency range we observed a power increase in the SV condition only. Our data provide evidence that gamma power increases are related to the predictability of an upcoming word based on the preceding sentence context, rather than to the integration of the incoming word's semantics into the preceding context. Further, our theta band data are compatible with the notion that theta band synchronization in sentence comprehension might be related to the detection of an error in the language input.

Time-Resolved Resting-State Functional Magnetic Resonance Imaging Analysis: Current Status, Challenges, and New Directions

Time-resolved analysis of resting-state functional magnetic resonance imaging (rs-fMRI) data allows researchers to extract more information about brain function than traditional functional connectivity analysis, yet a number of challenges in data analysis and interpretation remain. This article briefly summarizes common methods for time-resolved analysis and presents some of the pressing issues and opportunities in the field. From there, the discussion moves to interpretation of the network dynamics observed with rs-fMRI and the role that rs-fMRI can play in elucidating the large-scale organization of brain activity.

Breathing Rhythm and Pattern and Their Influence on Emotion

Breathing is a vital rhythmic motor behavior with a surprisingly broad influence on the brain and body. The apparent simplicity of breathing belies a complex neural control system, the breathing central pattern generator (bCPG), that exhibits diverse operational modes to regulate gas exchange and coordinate breathing with an array of behaviors. In this review, we focus on selected advances in our understanding of the bCPG. At the core of the bCPG is the preBötzinger complex (preBötC), which drives inspiratory rhythm via an unexpectedly sophisticated emergent mechanism. Synchronization dynamics underlying preBötC rhythmogenesis imbue the system with robustness and lability. These dynamics are modulated by inputs from throughout the brain and generate rhythmic, patterned activity that is widely distributed. The connectivity and an emerging literature support a link between breathing, emotion, and cognition that is becoming experimentally tractable. These advances bring great potential for elucidating function and dysfunction in breathing and other mammalian neural circuits.

Compromised hippocampus-striatum pathway as a potential imaging biomarker of mild-traumatic brain injury and posttraumatic stress disorder

OBJECTIVES

Military service members risk acquiring posttraumatic stress disorder (PTSD) and mild-traumatic brain injury (mTBI), with high comorbidity. Owing to overlapping symptomatology in chronic mTBI or postconcussion syndrome (PCS) and PTSD, it is difficult to assess the etiology of a patient's condition without objective measures. Using resting-state functional MRI in a novel framework, we tested the hypothesis that their neural signatures are characterized by functionally hyperconnected brain regions which are less variable over time. Additionally, we predicted that such connectivities possessed the highest ability in predicting the diagnostic membership of a novel subject (top-predictors) in addition to being statistically significant.

METHODS

U.S. Army Soldiers (N = 87) with PTSD and comorbid PCS + PTSD were recruited along with combat controls. Static and dynamic functional connectivities were evaluated. Group differences were obtained in accordance with our hypothesis. Machine learning classification (MLC) was employed to determine top predictors.

RESULTS

From whole-brain connectivity, we identified the hippocampus-striatum connectivity to be significantly altered in accordance with our hypothesis. Diffusion tractography revealed compromised white-matter integrity between aforementioned regions only in the PCS + PTSD group, suggesting a structural etiology for the PCS + PTSD group rather than being an extreme subset of PTSD. Employing MLC, connectivities provided worst-case accuracy of 84% (9% more than psychological measures). Additionally, the hippocampus-striatum connectivities were found to be top predictors and thus a potential biomarker of PTSD/mTBI.

CONCLUSIONS

PTSD/mTBI are associated with hippocampal-striatal hyperconnectivity from which it is difficult to disengage, leading to a habit-like response toward episodic traumatic memories, which fits well with behavioral manifestations of combat-related PTSD/mTBI. Hum Brain Mapp 38:2843-2864, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Linking the response properties of cells in auditory cortex with network architecture: cotuning versus lateral inhibition

The frequency-intensity receptive fields (RF) of neurons in primary auditory cortex (AI) are heterogeneous. Some neurons have V-shaped RFs, whereas others have enclosed ovoid RFs. Moreover, there is a wide range of temporal response profiles ranging from phasic to tonic firing. The mechanisms underlying this diversity of receptive field properties are yet unknown. Here we study the characteristics of thalamocortical (TC) and intracortical connectivity that give rise to the individual cell responses. Using a mouse auditory TC slice preparation, we found that the amplitude of synaptic responses in AI varies non-monotonically with the intensity of the stimulation in the medial geniculate nucleus (MGv). We constructed a network model of MGv and AI that was simulated using either rate model cells or in vitro neurons through an iterative procedure that used the recorded neural responses to reconstruct network activity. We compared the receptive fields and firing profiles obtained with networks configured to have either cotuned excitatory and inhibitory inputs or relatively broad, lateral inhibitory inputs. Each of these networks yielded distinct response properties consistent with those documented in vivo with natural stimuli. The cotuned network produced V-shaped RFs, phasic-tonic firing profiles, and predominantly monotonic rate-level functions. The lateral inhibitory network produced enclosed RFs with narrow frequency tuning, a variety of firing profiles, and robust non-monotonic rate-level functions. We conclude that both types of circuits must be present to account for the wide variety of responses observed in vivo.

Anesthesia differentially modulates spontaneous network dynamics by cortical area and layer

Anesthesia is widely used in medicine and research to achieve altered states of consciousness and cognition. Whereas changes to macroscopic cortical activity patterns by anesthesia measured at the spatial resolution of electroencephalography have been widely studied, modulation of mesoscopic and microscopic network dynamics by anesthesia remain poorly understood. To address this gap in knowledge, we recorded spontaneous mesoscopic (local field potential) and microscopic (multiunit activity) network dynamics in primary visual cortex (V1) and prefrontal cortex (PFC) of awake and isoflurane anesthetized ferrets (Mustela putoris furo). This approach allowed for examination of activity as a function of cortical area, cortical layer, and anesthetic depth with much higher spatial and temporal resolution than in previous studies. We hypothesized that a primary sensory area and an association cortical area would exhibit different patterns of network modulation by anesthesia due to their different functional roles. Indeed, we found effects specific to cortical area and cortical layer. V1 exhibited minimal changes in rhythmic structure with anesthesia but differential modulation of input layer IV. In contrast, anesthesia profoundly altered spectral power in PFC, with more uniform modulation across cortical layers. Our results demonstrate that anesthesia modulates spontaneous cortical activity in an area- and layer-specific manner. These finding provide the basis for 1) refining anesthesia monitoring algorithms, 2) reevaluating the large number of systems neuroscience studies performed in anesthetized animals, and 3) increasing our understanding of differential dynamics across cortical layers and areas.

A critical role for network structure in seizure onset: a computational modeling approach

Recent clinical work has implicated network structure as critically important in the initiation of seizures in people with idiopathic generalized epilepsies. In line with this idea, functional networks derived from the electroencephalogram (EEG) at rest have been shown to be significantly different in people with generalized epilepsy compared to controls. In particular, the mean node degree of networks from the epilepsy cohort was found to be statistically significantly higher than those of controls. However, the mechanisms by which these network differences can support recurrent transitions into seizures remain unclear. In this study, we use a computational model of the transition into seizure dynamics to explore the dynamic consequences of these differences in functional networks. We demonstrate that networks with higher mean node degree are more prone to generating seizure dynamics in the model and therefore suggest a mechanism by which increased mean node degree of brain networks can cause heightened ictogenicity.

Increased structural white and grey matter network connectivity compensates for functional decline in early multiple sclerosis

BACKGROUND

The pathology of multiple sclerosis (MS) consists of demyelination and neuronal injury, which occur early in the disease; yet, remission phases indicate repair. Whether and how the central nervous system (CNS) maintains homeostasis to counteract clinical impairment is not known.

OBJECTIVE

We analyse the structural connectivity of white matter (WM) and grey matter (GM) networks to understand the absence of clinical decline as the disease progresses.

METHODS

A total of 138 relapsing-remitting MS patients (classified into six groups by disease duration) and 32 healthy controls were investigated using 3-Tesla magnetic resonance imaging (MRI). Networks were analysed using graph theoretical approaches based on connectivity patterns derived from diffusion-tensor imaging with probabilistic tractography for WM and voxel-based morphometry and regional-volume-correlation matrix for GM.

RESULTS

In the first year after disease onset, WM networks evolved to a structure of increased modularity, strengthened local connectivity and increased local clustering while no clinical decline occurred. GM networks showed a similar dynamic of increasing modularity. This modified connectivity pattern mainly involved the cerebellum, cingulum and temporo-parietal regions. Clinical impairment was associated at later disease stages with a divergence of the network patterns.

CONCLUSION

Our findings suggest that network functionality in MS is maintained through structural adaptation towards increased local and modular connectivity, patterns linked to adaptability and homeostasis.

Integration in emerging social networks explains academic failure and success

Understanding the factors that explain academic failure and success of university students is a core interest of educational researchers, teachers, and managers. We demonstrate how the dynamic social networks that informally evolve between students can affect their academic performance. We closely followed the emergence of multiple social networks within a cohort of 226 undergraduate university students. They were strangers to each other on their first day at university, but developed densely knit social networks through time. We show that functional studying relationships tended to evolve from informal friendship relations. In a critical examination period after one year, these networks proved to be crucial: Socially isolated students had significantly lower examination grades and were more likely to drop out of university.

Academic success of students has been explained with a variety of individual and socioeconomic factors. Social networks that informally emerge within student communities can have an additional effect on their achievement. However, this effect of social ties is difficult to measure and quantify, because social networks are multidimensional and dynamically evolving within the educational context. We repeatedly surveyed a cohort of 226 engineering undergraduates between their first day at university and a crucial examination at the end of the academic year. We investigate how social networks emerge between previously unacquainted students and how integration in these networks explains academic success. Our study measures multiple important dimensions of social ties between students: their positive interactions, friendships, and studying relations. By using statistical models for dynamic network data, we are able to investigate the processes of social network formation in the cohort. We find that friendship ties informally evolve into studying relationships over the academic year. This process is crucial, as studying together with others, in turn, has a strong impact on students’ success at the examination. The results are robust to individual differences in socioeconomic background factors and to various indirect measures of cognitive abilities, such as prior academic achievement and being perceived as smart by other students. The findings underline the importance of understanding social network dynamics in educational settings. They call for the creation of university environments promoting the development of positive relationships in pursuit of academic success.

Parameters for burst detection

Bursts of action potentials within neurons and throughout networks are believed to serve roles in how neurons handle and store information, both in vivo and in vitro. Accurate detection of burst occurrences and durations are therefore crucial for many studies. A number of algorithms have been proposed to do so, but a standard method has not been adopted. This is due, in part, to many algorithms requiring the adjustment of multiple ad-hoc parameters and further post-hoc criteria in order to produce satisfactory results. Here, we broadly catalog existing approaches and present a new approach requiring the selection of only a single parameter: the number of spikes N comprising the smallest burst to consider. A burst was identified if N spikes occurred in less than T ms, where the threshold T was automatically determined from observing a probability distribution of inter-spike-intervals. Performance was compared vs. different classes of detectors on data gathered from in vitro neuronal networks grown over microelectrode arrays. Our approach offered a number of useful features including: a simple implementation, no need for ad-hoc or post-hoc criteria, and precise assignment of burst boundary time points. Unlike existing approaches, detection was not biased toward larger bursts, allowing identification and analysis of a greater range of neuronal and network dynamics.

Flexible Coordinator and Switcher Hubs for Adaptive Task Control

Functional connectivity (FC) studies have identified at least two large-scale neural systems that constitute cognitive control networks, the frontoparietal network (FPN) and cingulo-opercular network (CON). Control networks are thought to support goal-directed cognition and behavior. It was previously shown that the FPN flexibly shifts its global connectivity pattern according to task goal, consistent with a "flexible hub" mechanism for cognitive control. Our aim was to build on this finding to develop a functional cartography (a multimetric profile) of control networks in terms of dynamic network properties. We quantified network properties in (male and female) humans using a high-control-demand cognitive paradigm involving switching among 64 task sets. We hypothesized that cognitive control is enacted by the FPN and CON via distinct but complementary roles reflected in network dynamics. Consistent with a flexible "coordinator" mechanism, FPN connections were varied across tasks, while maintaining within-network connectivity to aid cross-region coordination. Consistent with a flexible "switcher" mechanism, CON regions switched to other networks in a task-dependent manner, driven primarily by reduced within-network connections to other CON regions. This pattern of results suggests FPN acts as a dynamic, global coordinator of goal-relevant information, while CON transiently disbands to lend processing resources to other goal-relevant networks. This cartography of network dynamics reveals a dissociation between two prominent cognitive control networks, suggesting complementary mechanisms underlying goal-directed cognition.SIGNIFICANCE STATEMENT Cognitive control supports a variety of behaviors requiring flexible cognition, such as rapidly switching between tasks. Furthermore, cognitive control is negatively impacted in a variety of mental illnesses. We used tools from network science to characterize the implementation of cognitive control by large-scale brain systems. This revealed that two systems, the frontoparietal (FPN) and cingulo-opercular (CON) networks, have distinct but complementary roles in controlling global network reconfigurations. The FPN exhibited properties of a flexible coordinator (orchestrating task changes), while CON acted as a flexible switcher (switching specific regions to other systems to lend processing resources). These findings reveal an underlying distinction in cognitive processes that may be applicable to clinical, educational, and machine learning work targeting cognitive flexibility.

Dynamic connectivity detection: an algorithm for determining functional connectivity change points in fMRI data

Recently there has been an increased interest in using fMRI data to study the dynamic nature of brain connectivity. In this setting, the activity in a set of regions of interest (ROIs) is often modeled using a multivariate Gaussian distribution, with a mean vector and covariance matrix that are allowed to vary as the experiment progresses, representing changing brain states. In this work, we introduce the Dynamic Connectivity Detection (DCD) algorithm, which is a data-driven technique to detect temporal change points in functional connectivity, and estimate a graph between ROIs for data within each segment defined by the change points. DCD builds upon the framework of the recently developed Dynamic Connectivity Regression (DCR) algorithm, which has proven efficient at detecting changes in connectivity for problems consisting of a small to medium (< 50) number of regions, but which runs into computational problems as the number of regions becomes large (>100). The newly proposed DCD method is faster, requires less user input, and is better able to handle high-dimensional data. It overcomes the shortcomings of DCR by adopting a simplified sparse matrix estimation approach and a different hypothesis testing procedure to determine change points. The application of DCD to simulated data, as well as fMRI data, illustrates the efficacy of the proposed method.