A brainstem integrator for self-location memory and positional homeostasis in zebrafish.

To track and control self-location, animals integrate their movements through space. Representations of self-location are observed in the mammalian hippocampal formation, but it is unknown if positional representations exist in more ancient brain regions, how they arise from integrated self-motion, and by what pathways they control locomotion. Here, in a head-fixed, fictive-swimming, virtual-reality preparation, we exposed larval zebrafish to a variety of involuntary displacements. They tracked these displacements and, many seconds later, moved toward their earlier location through corrective swimming ("positional homeostasis"). Whole-brain functional imaging revealed a network in the medulla that stores a memory of location and induces an error signal in the inferior olive to drive future corrective swimming. Optogenetically manipulating medullary integrator cells evoked displacement-memory behavior. Ablating them, or downstream olivary neurons, abolished displacement corrections. These results reveal a multiregional hindbrain circuit in vertebrates that integrates self-motion and stores self-location to control locomotor behavior.

Glia Accumulate Evidence that Actions Are Futile and Suppress Unsuccessful Behavior.

When a behavior repeatedly fails to achieve its goal, animals often give up and become passive, which can be strategic for preserving energy or regrouping between attempts. It is unknown how the brain identifies behavioral failures and mediates this behavioral-state switch. In larval zebrafish swimming in virtual reality, visual feedback can be withheld so that swim attempts fail to trigger expected visual flow. After tens of seconds of such motor futility, animals became passive for similar durations. Whole-brain calcium imaging revealed noradrenergic neurons that responded specifically to failed swim attempts and radial astrocytes whose calcium levels accumulated with increasing numbers of failed attempts. Using cell ablation and optogenetic or chemogenetic activation, we found that noradrenergic neurons progressively activated brainstem radial astrocytes, which then suppressed swimming. Thus, radial astrocytes perform a computation critical for behavior: they accumulate evidence that current actions are ineffective and consequently drive changes in behavioral states. VIDEO ABSTRACT.

Unbiased and efficient sampling of timeseries reveals redundancy of brain network and gradient structure.

Many studies in human neuroscience seek to understand the structure of brain networks and gradients. Few studies, however, have tested the redundancy between these outwardly distinct features. Here, we developed methods to directly enable such tests. We built on insights from linear algebra to develop methods for unbiased and efficient sampling of timeseries with network or gradient constraints. We used these methods to show considerable redundancy between popular definitions of network and gradient structure in functional MRI data. On the one hand, we found that network constraints largely accounted for the structure of three major gradients. On the other hand, we found that gradient constraints largely accounted for the structure of seven major networks. Our results imply that some networks and gradients may denote discrete and continuous representations of the same aspects of functional MRI data. We suggest that integrated explanations can reduce redundancy by avoiding the attribution of independent existence or function to these features.

Wiring cost and topological participation of the mouse brain connectome.

Brain connectomes are topologically complex systems, anatomically embedded in 3D space. Anatomical conservation of "wiring cost" explains many but not all aspects of these networks. Here, we examined the relationship between topology and wiring cost in the mouse connectome by using data from 461 systematically acquired anterograde-tracer injections into the right cortical and subcortical regions of the mouse brain. We estimated brain-wide weights, distances, and wiring costs of axonal projections and performed a multiscale topological and spatial analysis of the resulting weighted and directed mouse brain connectome. Our analysis showed that the mouse connectome has small-world properties, a hierarchical modular structure, and greater-than-minimal wiring costs. High-participation hubs of this connectome mediated communication between functionally specialized and anatomically localized modules, had especially high wiring costs, and closely corresponded to regions of the default mode network. Analyses of independently acquired histological and gene-expression data showed that nodal participation colocalized with low neuronal density and high expression of genes enriched for cognition, learning and memory, and behavior. The mouse connectome contains high-participation hubs, which are not explained by wiring-cost minimization but instead reflect competitive selection pressures for integrated network topology as a basis for higher cognitive and behavioral functions.

Correspondence between evoked and intrinsic functional brain network configurations.

Much of the literature exploring differences between intrinsic and task-evoked brain architectures has examined changes in functional connectivity patterns between specific brain regions. While informative, this approach overlooks important overall functional changes in hub organization and network topology that may provide insights about differences in integration between intrinsic and task-evoked states. Examination of changes in overall network organization, such as a change in the concentration of hub nodes or a quantitative change in network organization, is important for understanding the underlying processes that differ between intrinsic and task-evoked brain architectures. The present study used graph-theoretical techniques applied to publicly available neuroimaging data collected from a large sample of individuals (N = 202), and a within-subject design where resting-state and several task scans were collected from each participant as part of the Human Connectome Project. We demonstrate that differences between intrinsic and task-evoked brain networks are characterized by a task-general shift in high-connectivity hubs from primarily sensorimotor/auditory processing areas during the intrinsic state to executive control/salience network areas during task performance. In addition, we demonstrate that differences between intrinsic and task-evoked architectures are associated with changes in overall network organization, such as increases in network clustering, global efficiency and integration between modules. These findings offer a new perspective on the principles guiding functional brain organization by identifying unique and divergent properties of overall network organization between the resting-state and task performance. Hum Brain Mapp 38:1992-2007, 2017. © 2017 Wiley Periodicals, Inc.

State and Trait Components of Functional Connectivity: Individual Differences Vary with Mental State.

Resting-state functional connectivity, as measured by functional magnetic resonance imaging (fMRI), is often treated as a trait, used, for example, to draw inferences about individual differences in cognitive function, or differences between healthy or diseased populations. However, functional connectivity can also depend on the individual's mental state. In the present study, we examined the relative contribution of state and trait components in shaping an individual's functional architecture. We used fMRI data from a large, population-based human sample (N = 587, age 18-88 years), as part of the Cambridge Centre for Aging and Neuroscience (Cam-CAN), which were collected in three mental states: resting, performing a sensorimotor task, and watching a movie. Whereas previous studies have shown commonalities across mental states in the average functional connectivity across individuals, we focused on the effects of states on the pattern of individual differences in functional connectivity. We found that state effects were as important as trait effects in shaping individual functional connectivity patterns, each explaining an approximately equal amount of variance. This was true when we looked at aging, as one specific dimension of individual differences, as well as when we looked at generic aspects of individual variation. These results show that individual differences in functional connectivity consist of state-dependent aspects, as well as more stable, trait-like characteristics. Studying individual differences in functional connectivity across a wider range of mental states will therefore provide a more complete picture of the mechanisms underlying factors such as cognitive ability, aging, and disease. SIGNIFICANCE STATEMENT: The brain's functional architecture is remarkably similar across different individuals and across different mental states, which is why many studies use functional connectivity as a trait measure. Despite these trait-like aspects, functional connectivity varies over time and with changes in cognitive state. We measured connectivity in three different states to quantify the size of the trait-like component of functional connectivity, compared with the state-dependent component. Our results show that studying individual differences within one state (such as resting) uncovers only part of the relevant individual differences in brain function, and that the study of functional connectivity under multiple mental states is essential to disentangle connectivity differences that are transient versus those that represent more stable, trait-like characteristics of an individual.

Modern Methods for Interrogating the Human Connectome.

OBJECTIVES: Connectionist theories of brain function took hold with the seminal contributions of Norman Geschwind a half century ago. Modern neuroimaging techniques have expanded the scientific interest in the study of brain connectivity to include the intact as well as disordered brain. METHODS: In this review, we describe the most common techniques used to measure functional and structural connectivity, including resting state functional MRI, diffusion MRI, and electroencephalography and magnetoencephalography coherence. We also review the most common analytical approaches used for examining brain interconnectivity associated with these various imaging methods. RESULTS: This review presents a critical analysis of the assumptions, as well as methodological limitations, of each imaging and analysis approach. CONCLUSIONS: The overall goal of this review is to provide the reader with an introduction to evaluating the scientific methods underlying investigations that probe the human connectome.

Network topology and functional connectivity disturbances precede the onset of Huntington's disease.

Cognitive, motor and psychiatric changes in prodromal Huntington's disease have nurtured the emergent need for early interventions. Preventive clinical trials for Huntington's disease, however, are limited by a shortage of suitable measures that could serve as surrogate outcomes. Measures of intrinsic functional connectivity from resting-state functional magnetic resonance imaging are of keen interest. Yet recent studies suggest circumscribed abnormalities in resting-state functional magnetic resonance imaging connectivity in prodromal Huntington's disease, despite the spectrum of behavioural changes preceding a manifest diagnosis. The present study used two complementary analytical approaches to examine whole-brain resting-state functional magnetic resonance imaging connectivity in prodromal Huntington's disease. Network topology was studied using graph theory and simple functional connectivity amongst brain regions was explored using the network-based statistic. Participants consisted of gene-negative controls (n = 16) and prodromal Huntington's disease individuals (n = 48) with various stages of disease progression to examine the influence of disease burden on intrinsic connectivity. Graph theory analyses showed that global network interconnectivity approximated a random network topology as proximity to diagnosis neared and this was associated with decreased connectivity amongst highly-connected rich-club network hubs, which integrate processing from diverse brain regions. However, functional segregation within the global network (average clustering) was preserved. Functional segregation was also largely maintained at the local level, except for the notable decrease in the diversity of anterior insula intermodular-interconnections (participation coefficient), irrespective of disease burden. In contrast, network-based statistic analyses revealed patterns of weakened frontostriatal connections and strengthened frontal-posterior connections that evolved as disease burden increased. These disturbances were often related to long-range connections involving peripheral nodes and interhemispheric connections. A strong association was found between weaker connectivity and decreased rich-club organization, indicating that whole-brain simple connectivity partially expressed disturbances in the communication of highly-connected hubs. However, network topology and network-based statistic connectivity metrics did not correlate with key markers of executive dysfunction (Stroop Test, Trail Making Test) in prodromal Huntington's disease, which instead were related to whole-brain connectivity disturbances in nodes (right inferior parietal, right thalamus, left anterior cingulate) that exhibited multiple aberrant connections and that mediate executive control. Altogether, our results show for the first time a largely disease burden-dependent functional reorganization of whole-brain networks in prodromal Huntington's disease. Both analytic approaches provided a unique window into brain reorganization that was not related to brain atrophy or motor symptoms. Longitudinal studies currently in progress will chart the course of functional changes to determine the most sensitive markers of disease progression.

A wavelet method for modeling and despiking motion artifacts from resting-state fMRI time series.

The impact of in-scanner head movement on functional magnetic resonance imaging (fMRI) signals has long been established as undesirable. These effects have been traditionally corrected by methods such as linear regression of head movement parameters. However, a number of recent independent studies have demonstrated that these techniques are insufficient to remove motion confounds, and that even small movements can spuriously bias estimates of functional connectivity. Here we propose a new data-driven, spatially-adaptive, wavelet-based method for identifying, modeling, and removing non-stationary events in fMRI time series, caused by head movement, without the need for data scrubbing. This method involves the addition of just one extra step, the Wavelet Despike, in standard pre-processing pipelines. With this method, we demonstrate robust removal of a range of different motion artifacts and motion-related biases including distance-dependent connectivity artifacts, at a group and single-subject level, using a range of previously published and new diagnostic measures. The Wavelet Despike is able to accommodate the substantial spatial and temporal heterogeneity of motion artifacts and can consequently remove a range of high and low frequency artifacts from fMRI time series, that may be linearly or non-linearly related to physical movements. Our methods are demonstrated by the analysis of three cohorts of resting-state fMRI data, including two high-motion datasets: a previously published dataset on children (N=22) and a new dataset on adults with stimulant drug dependence (N=40). We conclude that there is a real risk of motion-related bias in connectivity analysis of fMRI data, but that this risk is generally manageable, by effective time series denoising strategies designed to attenuate synchronized signal transients induced by abrupt head movements. The Wavelet Despiking software described in this article is freely available for download at www.brainwavelet.org.

Dynamic change of global and local information processing in propofol-induced loss and recovery of consciousness.

Whether unique to humans or not, consciousness is a central aspect of our experience of the world. The neural fingerprint of this experience, however, remains one of the least understood aspects of the human brain. In this paper we employ graph-theoretic measures and support vector machine classification to assess, in 12 healthy volunteers, the dynamic reconfiguration of functional connectivity during wakefulness, propofol-induced sedation and loss of consciousness, and the recovery of wakefulness. Our main findings, based on resting-state fMRI, are three-fold. First, we find that propofol-induced anesthesia does not bear differently on long-range versus short-range connections. Second, our multi-stage design dissociated an initial phase of thalamo-cortical and cortico-cortical hyperconnectivity, present during sedation, from a phase of cortico-cortical hypoconnectivity, apparent during loss of consciousness. Finally, we show that while clustering is increased during loss of consciousness, as recently suggested, it also remains significantly elevated during wakefulness recovery. Conversely, the characteristic path length of brain networks (i.e., the average functional distance between any two regions of the brain) appears significantly increased only during loss of consciousness, marking a decrease of global information-processing efficiency uniquely associated with unconsciousness. These findings suggest that propofol-induced loss of consciousness is mainly tied to cortico-cortical and not thalamo-cortical mechanisms, and that decreased efficiency of information flow is the main feature differentiating the conscious from the unconscious brain.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena.

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

Neurobiologically realistic determinants of self-organized criticality in networks of spiking neurons.

Self-organized criticality refers to the spontaneous emergence of self-similar dynamics in complex systems poised between order and randomness. The presence of self-organized critical dynamics in the brain is theoretically appealing and is supported by recent neurophysiological studies. Despite this, the neurobiological determinants of these dynamics have not been previously sought. Here, we systematically examined the influence of such determinants in hierarchically modular networks of leaky integrate-and-fire neurons with spike-timing-dependent synaptic plasticity and axonal conduction delays. We characterized emergent dynamics in our networks by distributions of active neuronal ensemble modules (neuronal avalanches) and rigorously assessed these distributions for power-law scaling. We found that spike-timing-dependent synaptic plasticity enabled a rapid phase transition from random subcritical dynamics to ordered supercritical dynamics. Importantly, modular connectivity and low wiring cost broadened this transition, and enabled a regime indicative of self-organized criticality. The regime only occurred when modular connectivity, low wiring cost and synaptic plasticity were simultaneously present, and the regime was most evident when between-module connection density scaled as a power-law. The regime was robust to variations in other neurobiologically relevant parameters and favored systems with low external drive and strong internal interactions. Increases in system size and connectivity facilitated internal interactions, permitting reductions in external drive and facilitating convergence of postsynaptic-response magnitude and synaptic-plasticity learning rate parameter values towards neurobiologically realistic levels. We hence infer a novel association between self-organized critical neuronal dynamics and several neurobiologically realistic features of structural connectivity. The central role of these features in our model may reflect their importance for neuronal information processing.

Network connectivity in primary generalized tonic-clonic seizures.

OBJECTIVE: To determine EEG spatiospectral activation and connectivity in the generalized tonic-clonic seizure (GTCS) semiological subtypes. METHODS: 39 patients with genetic generalized epilepsy (GGE) who had GTCS (n = 58) during video-EEG monitoring were identified in the Vanderbilt Epilepsy database. GTCSs were classified as absence tonic-clonic, myoclonic tonic-clonic, or tonic-clonic. Patient characteristics and semiological features were compared. Spectral power and node degree, a network measure of connectivity, were calculated at two seizure epochs, electrographic and tonic-start. RESULTS: Different GTCS subtypes occurred within individual patients. At electrographic-onset, all subtypes activated midline frontal cortex at delta/theta and beta frequencies but differed in network connectivity. In all subtypes, GTCS evolution from electrographic to tonic-start associated with preserved beta frequency spectral power, but reduced connectivity and delta/theta power. CONCLUSIONS: Our findings suggest that at GTCS onset, the subtypes activate similar cortical regions and their different initial semiologies relate to their distinct onset long-range connectivity. Upon transition to the tonic-start epoch, the ictal activity is predominantly conveyed by ß frequency activity and connectivity. SIGNIFICANCE: Future neurostimulation therapies for medically intractable GTCSs may target the same brain regions for all GTCS subtypes and may be most effective prior to the tonic-start epoch.

Weight-conserving characterization of complex functional brain networks.

Complex functional brain networks are large networks of brain regions and functional brain connections. Statistical characterizations of these networks aim to quantify global and local properties of brain activity with a small number of network measures. Important functional network measures include measures of modularity (measures of the goodness with which a network is optimally partitioned into functional subgroups) and measures of centrality (measures of the functional influence of individual brain regions). Characterizations of functional networks are increasing in popularity, but are associated with several important methodological problems. These problems include the inability to characterize densely connected and weighted functional networks, the neglect of degenerate topologically distinct high-modularity partitions of these networks, and the absence of a network null model for testing hypotheses of association between observed nontrivial network properties and simple weighted connectivity properties. In this study we describe a set of methods to overcome these problems. Specifically, we generalize measures of modularity and centrality to fully connected and weighted complex networks, describe the detection of degenerate high-modularity partitions of these networks, and introduce a weighted-connectivity null model of these networks. We illustrate our methods by demonstrating degenerate high-modularity partitions and strong correlations between two complementary measures of centrality in resting-state functional magnetic resonance imaging (MRI) networks from the 1000 Functional Connectomes Project, an open-access repository of resting-state functional MRI datasets. Our methods may allow more sound and reliable characterizations and comparisons of functional brain networks across conditions and subjects.

Role of the Nucleus Basalis as a Key Network Node in Temporal Lobe Epilepsy.

OBJECTIVE: To determine whether the nucleus basalis of Meynert (NBM) may be a key network structure of altered functional connectivity in temporal lobe epilepsy (TLE), we examined fMRI with network-based analyses. METHODS: We acquired resting-state fMRI in 40 adults with TLE and 40 matched healthy control participants. We calculated functional connectivity of NBM and used multiple complementary network-based analyses to explore the importance of NBM in TLE networks without biasing our results by our approach. We compared patients to controls and examined associations of network properties with disease metrics and neurocognitive testing. RESULTS: We observed marked decreases in connectivity between NBM and the rest of the brain in patients with TLE (0.91 ± 0.88, mean ± SD) vs controls (1.96 ± 1.13, p < 0.001, t test). Larger decreases in connectivity between NBM and fronto-parietal-insular regions were associated with higher frequency of consciousness-impairing seizures (r = -0.41, p = 0.008, Pearson). A core network of altered nodes in TLE included NBM ipsilateral to the epileptogenic side and bilateral limbic structures. Furthermore, normal community affiliation of ipsilateral NBM was lost in patients, and this structure displayed the most altered clustering coefficient of any node examined (3.46 ± 1.17 in controls vs 2.23 ± 0.93 in patients). Abnormal connectivity between NBM and subcortical arousal community was associated with modest neurocognitive deficits. Finally, a logistic regression model incorporating connectivity properties of ipsilateral NBM successfully distinguished patients from control datasets with moderately high accuracy (78%). CONCLUSIONS: These results suggest that while NBM is rarely studied in epilepsy, it may be one of the most perturbed network nodes in TLE, contributing to widespread neural effects in this disabling disorder.

Complex network measures of brain connectivity: uses and interpretations.

Brain connectivity datasets comprise networks of brain regions connected by anatomical tracts or by functional associations. Complex network analysis-a new multidisciplinary approach to the study of complex systems-aims to characterize these brain networks with a small number of neurobiologically meaningful and easily computable measures. In this article, we discuss construction of brain networks from connectivity data and describe the most commonly used network measures of structural and functional connectivity. We describe measures that variously detect functional integration and segregation, quantify centrality of individual brain regions or pathways, characterize patterns of local anatomical circuitry, and test resilience of networks to insult. We discuss the issues surrounding comparison of structural and functional network connectivity, as well as comparison of networks across subjects. Finally, we describe a Matlab toolbox (http://www.brain-connectivity-toolbox.net) accompanying this article and containing a collection of complex network measures and large-scale neuroanatomical connectivity datasets.

Small-world properties of nonlinear brain activity in schizophrenia.

A disturbance in the interactions between distributed cortical regions may underlie the cognitive and perceptual dysfunction associated with schizophrenia. In this article, nonlinear measures of cortical interactions and graph-theoretical metrics of network topography are combined to investigate this schizophrenia "disconnection hypothesis." This is achieved by analyzing the spatiotemporal structure of resting state scalp EEG data previously acquired from 40 young subjects with a recent first episode of schizophrenia and 40 healthy matched controls. In each subject, a method of mapping the topography of nonlinear interactions between cortical regions was applied to a widely distributed array of these data. The resulting nonlinear correlation matrices were converted to weighted graphs. The path length (a measure of large-scale network integration), clustering coefficient (a measure of "cliquishness"), and hub structure of these graphs were used as metrics of the underlying brain network activity. The graphs of both groups exhibited high levels of local clustering combined with comparatively short path lengths--features consistent with a "small-world" topology--as well as the presence of strong, central hubs. The graphs in the schizophrenia group displayed lower clustering and shorter path lengths in comparison to the healthy group. Whilst still "small-world," these effects are consistent with a subtle randomization in the underlying network architecture--likely associated with a greater number of links connecting disparate clusters. This randomization may underlie the cognitive disturbances characteristic of schizophrenia.

Symbiotic relationship between brain structure and dynamics.

BACKGROUND: Brain structure and dynamics are interdependent through processes such as activity-dependent neuroplasticity. In this study, we aim to theoretically examine this interdependence in a model of spontaneous cortical activity. To this end, we simulate spontaneous brain dynamics on structural connectivity networks, using coupled nonlinear maps. On slow time scales structural connectivity is gradually adjusted towards the resulting functional patterns via an unsupervised, activity-dependent rewiring rule. The present model has been previously shown to generate cortical-like, modular small-world structural topology from initially random connectivity. We provide further biophysical justification for this model and quantitatively characterize the relationship between structure, function and dynamics that accompanies the ensuing self-organization. RESULTS: We show that coupled chaotic dynamics generate ordered and modular functional patterns, even on a random underlying structural connectivity. Consequently, structural connectivity becomes more modular as it rewires towards these functional patterns. Functional networks reflect the underlying structural networks on slow time scales, but significantly less so on faster time scales. In spite of ordered functional topology, structural networks remain robustly interconnected--and therefore small-world--due to the presence of central, inter-modular hub nodes. The noisy dynamics of these hubs enable them to persist despite ongoing rewiring and despite their comparative absence in functional networks. CONCLUSION: Our results outline a theoretical mechanism by which brain dynamics may facilitate neuroanatomical self-organization. We find time scale dependent differences between structural and functional networks. These differences are likely to arise from the distinct dynamics of central structural nodes.

Simulation of neuronal death and network recovery in a computational model of distributed cortical activity.

OBJECTIVES: The authors utilize a model of activity-dependent neuronal plasticity to study the interplay between synaptogenesis, neuronal death, and neurogenesis on the resulting pattern of neuronal connectivity. DESIGN: A mathematical model of neuronal network activity was employed, with plasticity instantiated by an activity-dependent rewiring rule. In particular, the authors modeled a neural system as a collection of "nodes" (neural subsystems) connected by "links" (anatomical connectivity). Neuronal damage was simulated by deletion of nodes in this evolving network through either random or targeted attack. Neurogenesis was likewise simulated by insertion of new nodes with random connections. MEASUREMENTS: Local and global structural network properties were characterized using the metrics of local and global "efficiency," and network "reachability." RESULTS: Activity-dependent plasticity yields a network that is robust to random node deletion, with preservation of a "small-world" architecture, characterized by high local and global efficiency. In contrast, targeted deletion of central nodes leads to a drop in reachability and global efficiency, with a consequent loss of small-world properties. Simulated neurogenesis is able to compensate for this targeted cell loss even when rates of new cell formation are considerably slower than that of simulated cell death. CONCLUSIONS: The rapid growth of computational neuroscience enables to study the interplay between neuronal plasticity and cell death in computational models of brain network activity. Although the current simulations lack much of the rich physiology of real neuronal systems, they nevertheless allow us to make tentative hypotheses of the effects of neuronal lesions on the resulting neuroanatomical connectivity networks.