Author Summary Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. We constructed a multilayer network model of neuron-astrocyte connectivity based on calcium activity in mixed cortical cultures, and used this model to evaluate the effect of glutamatergic inhibition and mechanical injury on network topology. We found that injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in injured networks pretreated with a glutamate receptor antagonist. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations. Abstract Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. Despite recent advances in understanding neuron-astrocyte signaling, little is known about astrocytic modulation of neuronal activity at the population level, particularly in disease or following injury. We used high-speed calcium imaging of mixed cortical cultures in vitro to determine how population activity changes after disruption of glutamatergic signaling and mechanical injury. We constructed a multilayer network model of neuron-astrocyte connectivity, which captured distinct topology and response behavior from single-cell-type networks. mGluR 5 inhibition decreased neuronal activity, but did not on its own disrupt functional connectivity or network topology. In contrast, injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in networks pretreated with mGluR 5 inhibition. Comparison of spatial and functional connectivity revealed that functional connectivity is largely independent of spatial proximity at the microscale, but mechanical injury increased the spatial-functional correlation. Finally, we found that astrocyte segments of the same cell often belong to separate functional communities based on neuronal connectivity, suggesting that astrocyte segments function as independent entities. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations.

Recent advances in computational models of signal propagation and routing in the human brain have underscored the critical role of white-matter structure. A complementary approach has utilized the framework of network control theory to better understand how white matter constrains the manner in which a region or set of regions can direct or control the activity of other regions. Despite the potential for both of these approaches to enhance our understanding of the role of network structure in brain function, little work has sought to understand the relations between them. Here, we seek to explicitly bridge computational models of communication and principles of network control in a conceptual review of the current literature. By drawing comparisons between communication and control models in terms of the level of abstraction, the dynamical complexity, the dependence on network attributes, and the interplay of multiple spatiotemporal scales, we highlight the convergence of and distinctions between the two frameworks. Based on the understanding of the intertwined nature of communication and control in human brain networks, this work provides an integrative perspective for the field and outlines exciting directions for future work. Author Summary Models of communication in brain networks have been essential in building a quantitative understanding of the relationship between structure and function. More recently, control-theoretic models have also been applied to brain networks to quantify the response of brain networks to exogenous and endogenous perturbations. Mechanistically, both of these frameworks investigate the role of interregional communication in determining the behavior and response of the brain. Theoretically, both of these frameworks share common features, indicating the possibility of combining the two approaches. Drawing on a large body of past and ongoing works, this review presents a discussion of convergence and distinctions between the two approaches, and argues for the development of integrated models at the confluence of the two frameworks, with potential applications to various topics in neuroscience.

Author Summary A central question in network neuroscience is how the structure of the brain constrains the patterns of communication dynamics that underlie function. At the mesoscale of network organization, this question has been examined through the lens of modularity. Recent work has demonstrated a diversity in the mesoscale architecture of the human connectome. Further diversity in the characterization of structural brain networks is introduced by the fact that the distribution of edge weights in a network depends on the precise empirical measurement whose value is assigned to an edge. This paper explores network controllability in light of the variety of community interaction motifs and edge weight distributions that may be used to characterize structural brain networks. Abstract The human brain displays rich communication dynamics that are thought to be particularly well-reflected in its marked community structure. Yet, the precise relationship between community structure in structural brain networks and the communication dynamics that can emerge therefrom is not well understood. In addition to offering insight into the structure-function relationship of networked systems, such an understanding is a critical step toward the ability to manipulate the brain’s large-scale dynamical activity in a targeted manner. We investigate the role of community structure in the controllability of structural brain networks. At the region level, we find that certain network measures of community structure are sometimes statistically correlated with measures of linear controllability. However, we then demonstrate that this relationship depends on the distribution of network edge weights. We highlight the complexity of the relationship between community structure and controllability by performing numerical simulations using canonical graph models with varying mesoscale architectures and edge weight distributions. Finally, we demonstrate that weighted subgraph centrality , a measure rooted in the graph spectrum, and which captures higher order graph architecture, is a stronger and more consistent predictor of controllability. Our study contributes to an understanding of how the brain’s diverse mesoscale structure supports transient communication dynamics.

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. Author Summary 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.

Author Summary Network neuroscience applied to epileptic brains seeks to identify pathological neural connections that promote and maintain seizures, and holds promise to guide surgical planning in patients with intractable epilepsy. However, sampling of the epileptic network in intracranial EEG recording is limited by the choice of where to place intracranial electrodes, which is highly variable within and between epilepsy centers. The effect of incomplete spatial sampling generated by sparse electrode placement on network statistics is unknown. Here, we determine the sensitivity of several network statistics to incomplete spatial sampling, and we propose a method using electrode subsampling to determine patient-specific confidence intervals in network model predictions. Abstract Network neuroscience applied to epilepsy holds promise to map pathological networks, localize seizure generators, and inform targeted interventions to control seizures. However, incomplete sampling of the epileptic brain because of sparse placement of intracranial electrodes may affect model results. In this study, we evaluate the sensitivity of several published network measures to incomplete spatial sampling and propose an algorithm using network subsampling to determine confidence in model results. We retrospectively evaluated intracranial EEG data from 28 patients implanted with grid, strip, and depth electrodes during evaluation for epilepsy surgery. We recalculated global and local network metrics after randomly and systematically removing subsets of intracranial EEG electrode contacts. We found that sensitivity to incomplete sampling varied significantly across network metrics. This sensitivity was largely independent of whether seizure onset zone contacts were targeted or spared from removal. We present an algorithm using random subsampling to compute patient-specific confidence intervals for network localizations. Our findings highlight the difference in robustness between commonly used network metrics and provide tools to assess confidence in intracranial network localization. We present these techniques as an important step toward translating personalized network models of seizures into rigorous, quantitative approaches to invasive therapy.

The brain is a complex, multiscale dynamical system composed of many interacting regions. Knowledge of the spatiotemporal organization of these interactions is critical for establishing a solid understanding of the brain’s functional architecture and the relationship between neural dynamics and cognition in health and disease. The possibility of studying these dynamics through careful analysis of neuroimaging data has catalyzed substantial interest in methods that estimate time-resolved fluctuations in functional connectivity (often referred to as “dynamic” or time-varying functional connectivity; TVFC). At the same time, debates have emerged regarding the application of TVFC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive and behavioral relevance of resting TVFC. These and other unresolved issues complicate interpretation of resting TVFC findings and limit the insights that can be gained from this promising new research area. This article brings together scientists with a variety of perspectives on resting TVFC to review the current literature in light of these issues. We introduce core concepts, define key terms, summarize controversies and open questions, and present a forward-looking perspective on how resting TVFC analyses can be rigorously and productively applied to investigate a wide range of questions in cognitive and systems neuroscience.

Data analysis techniques from network science have fundamentally improved our understanding of neural systems and the complex behaviors that they support. Yet the restriction of network techniques to the study of pairwise interactions prevents us from taking into account intrinsic topological features such as cavities that may be crucial for system function. To detect and quantify these topological features, we must turn to algebro-topological methods that encode data as a simplicial complex built from sets of interacting nodes called simplices. We then use the relations between simplices to expose cavities within the complex, thereby summarizing its topological features. Here we provide an introduction to persistent homology, a fundamental method from applied topology that builds a global descriptor of system structure by chronicling the evolution of cavities as we move through a combinatorial object such as a weighted network. We detail the mathematics and perform demonstrative calculations on the mouse structural connectome, synapses in C. elegans , and genomic interaction data. Finally, we suggest avenues for future work and highlight new advances in mathematics ready for use in neural systems. Author Summary For the network neuroscientist, this exposition aims to communicate both the mathematics and the advantages of using tools from applied topology for the study of neural systems. Using data from the mouse connectome, electrical and chemical synapses in C. elegans , and chromatin interaction data, we offer example computations and applications to further demonstrate the power of topological data analysis in neuroscience. Finally, we expose the reader to novel developments in applied topology and relate these developments to current questions and methodological difficulties in network neuroscience.

Author Summary Brain stimulation devices capable of perturbing the physiological state of neural systems are rapidly gaining popularity for their potential to treat neurological and psychiatric disease. A root problem is that underlying dysfunction spans a large-scale network of brain regions, requiring the ability to control the complex interactions between multiple brain areas. Here, we use tools from network control theory to examine the dynamic reconfiguration of functionally interacting neuronal ensembles during targeted neurostimulation of cortical and subcortical brain structures. We demonstrate the ability to predictably reconfigure patterns of interactions between functional brain areas by modulating the strength and location of stimulation. Our findings have high significance for designing stimulation protocols capable of modulating distributed neural circuits in the human brain. Abstract Chronically implantable neurostimulation devices are becoming a clinically viable option for treating patients with neurological disease and psychiatric disorders. Neurostimulation offers the ability to probe and manipulate distributed networks of interacting brain areas in dysfunctional circuits. Here, we use tools from network control theory to examine the dynamic reconfiguration of functionally interacting neuronal ensembles during targeted neurostimulation of cortical and subcortical brain structures. By integrating multimodal intracranial recordings and diffusion-weighted imaging from patients with drug-resistant epilepsy, we test hypothesized structural and functional rules that predict altered patterns of synchronized local field potentials. We demonstrate the ability to predictably reconfigure functional interactions depending on stimulation strength and location. Stimulation of areas with structurally weak connections largely modulates the functional hubness of downstream areas and concurrently propels the brain towards more difficult-to-reach dynamical states. By using focal perturbations to bridge large-scale structure, function, and markers of behavior, our findings suggest that stimulation may be tuned to influence different scales of network interactions driving cognition.

Dynamic functional connectivity reflects the spatiotemporal organization of spontaneous brain activity in health and disease. Dynamic functional connectivity may be susceptible to artifacts induced by participant motion. This report provides a systematic evaluation of 12 commonly used participant-level confound regression strategies designed to mitigate the effects of micromovements in a sample of 393 youths (ages 8–22 years). Each strategy was evaluated according to a number of benchmarks, including (a) the residual association between participant motion and edge dispersion, (b) distance-dependent effects of motion on edge dispersion, (c) the degree to which functional subnetworks could be identified by multilayer modularity maximization, and (d) measures of module reconfiguration, including node flexibility and node promiscuity. Results indicate variability in the effectiveness of the evaluated pipelines across benchmarks. Methods that included global signal regression were the most consistently effective de-noising strategies. Author Summary Dynamic functional connectivity reflects the spatiotemporal organization of spontaneous brain activity in health and disease, but it can be susceptible to motion artifacts. Here we provide a systematic evaluation of 12 commonly used participant-level confound regression strategies designed to mitigate the effects of micromovements in a sample of 393 youths (ages 8–22 years). Each strategy was evaluated according to the residual association between participant motion and edge dispersion, distance-dependent effects of motion on edge dispersion, the degree to which functional subnetworks could be identified by multilayer modularity maximization, and measures of module reconfiguration, including node flexibility and node promiscuity. Results indicate variability in the effectiveness of the evaluated pipelines across benchmarks. Methods that included global signal regression were the most consistently effective de-noising strategies.

Connectomics is an integral part of network neuroscience. The field has undergone rapid expansion over recent years and increasingly involves a blend of experimental and computational approaches to brain connectivity. This Focus Feature on “New Trends in Connectomics” aims to track the progress of the field and its many applications across different neurobiological systems and species.

Value guides behavior. With knowledge of stimulus values and action consequences, behaviors that maximize expected reward can be selected. Prior work has identified several brain structures critical for representing both stimuli and their values. Yet, it remains unclear how these structures interact with one another and with other regions of the brain to support the dynamic acquisition of value-related knowledge. Here, we use a network neuroscience approach to examine how BOLD functional networks change as 20 healthy human subjects learn the values of novel visual stimuli over the course of four consecutive days. We show that connections between regions of the visual, frontal, and cingulate cortices become stronger as learning progresses, with some of these changes being specific to the type of feedback received during learning. These results demonstrate that functional networks dynamically track behavioral improvement in value judgments, and that interactions between network communities form predictive biomarkers of learning. Author Summary Rational human behavior is the pursuit of actions that maximize expected reward. These rewards can be understood as stimulus-value contingencies, learned by experience throughout our lives. Various structures have been recognized to participate in these learning processes. Yet, an understanding of how these structures interact with one another and with other brain regions remains vastly unexplored. Here, we propose a novel analytical framework utilizing and extending techniques from the dynamic network neuroscience to ask “How do our brains change when we learn values?” We find that interactions between sensory and fronto-cingulate structures grow stronger as learning progresses, bringing together several isolated findings in the cognitive neuroscience of value-based behavior and extending our understanding of human learning in general.

AUTHOR SUMMARY The human brain is characterized by a complex pattern of anatomical wiring, in the form of white-matter tracts that link large volumes of neural tissue. The organization of this pattern is likely driven by many factors, including evolutionary adaptability, robustness to perturbations, and a separation of the timescales necessary to produce a diverse repertoire of neural dynamics. In this study, we sought to disentangle two such factors—the drive to decrease the cost of wiring, and the putative drive to increase the efficiency of the network topology—and we explored the impacts of these factors on the brain’s modular organization. The contributions of this work include a new algorithmic approach to community detection and novel insights into the role of modules in human brain function. Abstract Brain networks are expected to be modular. However, existing techniques for estimating a network’s modules make it difficult to assess the influence of organizational principles such as wiring cost reduction on the detected modules. Here we present a modification of an existing module detection algorithm that allowed us to focus on connections that are unexpected under a cost-reduction wiring rule and to identify modules from among these connections. We applied this technique to anatomical brain networks and showed that the modules we detected differ from those detected using the standard technique. We demonstrated that these novel modules are spatially distributed, exhibit unique functional fingerprints, and overlap considerably with rich clubs, giving rise to an alternative and complementary interpretation of the functional roles of specific brain regions. Finally, we demonstrated that, using the modified module detection approach, we can detect modules in a developmental dataset that track normative patterns of maturation. Collectively, these findings support the hypothesis that brain networks are composed of modules and provide additional insight into the function of those modules.

AUTHOR SUMMARY Our ability to thoughtfully engage with the world around us changes appreciably as we transition from childhood to adulthood. Yet, how our brains develop to enable that change remains far from understood. Here we used network science—traditionally applied to the study of social networks like Facebook or Twitter—and machine learning to show that growing cognitive abilities are accompanied by greater flexibility of brain regions within distributed networks. This flexibility is greatest in the executive system, which is critical for higher-order cognitive functions and increases in expression and flexibility from childhood to young adulthood. These results suggest that healthy development is facilitated by an increasing precedence of executive networks and a greater switching of the regions and interactions subserving these networks. Abstract Cognitive function evolves significantly over development, enabling flexible control of human behavior. Yet, how these functions are instantiated in spatially distributed and dynamically interacting networks, or graphs, that change in structure from childhood to adolescence is far from understood. Here we applied a novel machine-learning method to track continuously overlapping and time-varying subgraphs in the brain at rest within a sample of 200 healthy youth (ages 8–11 and 19–22) drawn from the Philadelphia Neurodevelopmental Cohort. We uncovered a set of subgraphs that capture surprisingly integrated and dynamically changing interactions among known cognitive systems. We observed that subgraphs that were highly expressed were especially transient, flexibly switching between high and low expression over time. This transience was particularly salient in a subgraph predominantly linking frontoparietal regions of the executive system, which increases in both expression and flexibility from childhood to young adulthood. Collectively, these results suggest that healthy development is accompanied by an increasing precedence of executive networks and a greater switching of the regions and interactions subserving these networks.