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.

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.

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.

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.

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.

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.

Methods and Considerations for Dynamic Analysis of Functional MR Imaging Data.

Functional MR imaging (fMR imaging) studies have recently begun to examine spontaneous changes in interregional interactions (functional connectivity) over seconds to minutes, and their relation to natural shifts in cognitive and physiologic states. This practice opens the potential for uncovering structured, transient configurations of coordinated brain activity whose features may provide novel cognitive and clinical biomarkers. However, analysis of these time-varying phenomena requires careful differentiation between neural and nonneural contributions to the fMR imaging signal and thorough validation and statistical testing. In this article, the authors present an overview of methodological and interpretational considerations in this emerging field.

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.

Constraints and spandrels of interareal connectomes.

Interareal connectomes are whole-brain wiring diagrams of white-matter pathways. Recent studies have identified modules, hubs, module hierarchies and rich clubs as structural hallmarks of these wiring diagrams. An influential current theory postulates that connectome modules are adequately explained by evolutionary pressures for wiring economy, but that the other hallmarks are not explained by such pressures and are therefore less trivial. Here, we use constraint network models to test these postulates in current gold-standard vertebrate and invertebrate interareal-connectome reconstructions. We show that empirical wiring-cost constraints inadequately explain connectome module organization, and that simultaneous module and hub constraints induce the structural byproducts of hierarchies and rich clubs. These byproducts, known as spandrels in evolutionary biology, include the structural substrate of the default-mode network. Our results imply that currently standard connectome characterizations are based on circular analyses or double dipping, and we emphasize an integrative approach to future connectome analyses for avoiding such pitfalls.

Regional expression of the MAPT gene is associated with loss of hubs in brain networks and cognitive impairment in Parkinson disease and progressive supranuclear palsy.

Abnormalities of tau protein are central to the pathogenesis of progressive supranuclear palsy, whereas haplotype variation of the tau gene MAPT influences the risk of Parkinson disease and Parkinson's disease dementia. We assessed whether regional MAPT expression might be associated with selective vulnerability of global brain networks to neurodegenerative pathology. Using task-free functional magnetic resonance imaging in progressive supranuclear palsy, Parkinson disease, and healthy subjects (n = 128), we examined functional brain networks and measured the connection strength between 471 gray matter regions. We obtained MAPT and SNCA microarray expression data in healthy subjects from the Allen brain atlas. Regional connectivity varied according to the normal expression of MAPT. The regional expression of MAPT correlated with the proportionate loss of regional connectivity in Parkinson's disease. Executive cognition was impaired in proportion to the loss of hub connectivity. These effects were not seen with SNCA, suggesting that alpha-synuclein pathology is not mediated through global network properties. The results establish a link between regional MAPT expression and selective vulnerability of functional brain networks to neurodegeneration.

Default Mode Hypoconnectivity Underlies a Sex-Related Autism Spectrum.

BACKGROUND: Females and males differ significantly in the prevalence and presentation of autism spectrum conditions. One theory of this effect postulates that autistic traits lie on a sex-related continuum in the general population, and autism represents the extreme male end of this spectrum. This theory predicts that any feature of autism in males should 1) be present in autistic females, 2) differentiate between the sexes in the typical population, and 3) correlate with autistic traits. We tested these three predictions for default mode network (DMN) hypoconnectivity during the resting state, one of the most robustly found neurobiological differences in autism. METHODS: We analyzed a primary dataset of adolescents (N = 121, 12-18 years of age) containing a relatively large number of females and a replication multisite dataset including children, adolescents, and adults (N = 980, 6-58 years of age). We quantified the average connectivity between DMN regions and tested for group differences and correlation with behavioral performance using robust regression. RESULTS: We found significant differences in DMN intraconnectivity between female controls and females with autism (p = .001 in the primary dataset; p = .009 in the replication dataset), and between female controls and male controls (p = .036 in the primary dataset; p = .002 in the replication dataset). We also found a significant correlation between DMN intraconnectivity and performance on a mentalizing task (p = .001) in the primary dataset. CONCLUSIONS: Collectively, these findings provide the first evidence for DMN hypoconnectivity as a behaviorally relevant neuroimaging phenotype of the sex-related spectrum of autistic traits, of which autism represents the extreme case.

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.

Altered Hub Functioning and Compensatory Activations in the Connectome: A Meta-Analysis of Functional Neuroimaging Studies in Schizophrenia.

BACKGROUND: Functional neuroimaging studies of schizophrenia have identified abnormal activations in many brain regions. In an effort to interpret these findings from a network perspective, we carried out a meta-analysis of this literature, mapping anatomical locations of under- and over-activation to the topology of a normative human functional connectome. METHODS: We included 314 task-based functional neuroimaging studies including more than 5000 patients with schizophrenia and over 5000 controls. Coordinates of significant under- or over-activations in patients relative to controls were mapped to nodes of a normative connectome defined by a prior meta-analysis of 1641 functional neuroimaging studies of task-related activation in healthy volunteers. RESULTS: Under-activations and over-activations were reported in a wide diversity of brain regions. Both under- and over-activations were significantly more likely to be located in hub nodes that constitute the "rich club" or core of the normative connectome. In a subset of 121 studies that reported both under- and over-activations in the same patients, we found that, in network terms, these abnormalities were located in close topological proximity to each other. Under-activation in a peripheral node was more frequently associated specifically with over-activation of core nodes than with over-activation of another peripheral node. CONCLUSIONS: Although schizophrenia is associated with altered brain functional activation in a wide variety of regions, abnormal responses are concentrated in hubs of the normative connectome. Task-specific under-activation in schizophrenia is accompanied by over-activation of topologically central, less functionally specialized network nodes, which may represent a compensatory response.

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.

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.

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 unifying framework for measuring weighted rich clubs.

Network analysis can help uncover meaningful regularities in the organization of complex systems. Among these, rich clubs are a functionally important property of a variety of social, technological and biological networks. Rich clubs emerge when nodes that are somehow prominent or 'rich' (e.g., highly connected) interact preferentially with one another. The identification of rich clubs is non-trivial, especially in weighted networks, and to this end multiple distinct metrics have been proposed. Here we describe a unifying framework for detecting rich clubs which intuitively generalizes various metrics into a single integrated method. This generalization rests upon the explicit incorporation of randomized control networks into the measurement process. We apply this framework to real-life examples, and show that, depending on the selection of randomized controls, different kinds of rich-club structures can be detected, such as topological and weighted rich clubs.

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.

Fledgling pathoconnectomics of psychiatric disorders.

Pathoconnectomics, the mapping of abnormal brain networks, is a popular current framework for the study of brain dysfunction in psychiatric disorders. In this review we evaluate the conceptual foundations of this framework, describe the construction and analysis of empirical models of brain networks or connectomes, and summarize recent reports of the large-scale whole-brain connectome organization of two candidate brain-network disorders, schizophrenia and autism. We consider the evidence for the abnormal brain-network nature of psychiatric disorders and find it inconclusive. For instance, although there is some evidence for more random whole-brain network organization in schizophrenia and autism, future studies need to determine if these and other observed brain-network abnormalities represent sufficient phenotypes of psychiatric disorders, in order to validate pathoconnectomics as a scientific and clinical framework.