Petra E. Vértes

Cognitive relevance of the community structure of the human brain functional coactivation network

There is growing interest in the complex topology of human brain functional networks, often measured using resting-state functional MRI (fMRI). Here, we used a meta-analysis of the large primary literature that used fMRI or PET to measure task-related activation (>1,600 studies; 1985–2010). We estimated the similarity (Jaccard index) of the activation patterns across experimental tasks between each pair of 638 brain regions. This continuous coactivation matrix was used to build a weighted graph to characterize network topology. The coactivation network was modular, with occipital, central, and default-mode modules predominantly coactivated by specific cognitive domains (perception, action, and emotion, respectively). It also included a rich club of hub nodes, located in parietal and prefrontal cortex and often connected over long distances, which were coactivated by a diverse range of experimental tasks. Investigating the topological role of edges between a deactivated and an activated node, we found that such competitive interactions were most frequent between nodes in different modules or between an activated rich-club node and a deactivated peripheral node. Many aspects of the coactivation network were convergent with a connectivity network derived from resting state fMRI data (n = 27, healthy volunteers); although the connectivity network was more parsimoniously connected and differed in the anatomical locations of some hubs. We conclude that the community structure of human brain networks is relevant to cognitive function. Deactivations may play a role in flexible reconfiguration of the network according to cognitive demand, varying the integration between modules, and between the periphery and a central rich club.

Simple models of human brain functional networks

Human brain functional networks are embedded in anatomical space and have topological properties—small-worldness, modularity, fat-tailed degree distributions—that are comparable to many other complex networks. Although a sophisticated set of measures is available to describe the topology of brain networks, the selection pressures that drive their formation remain largely unknown. Here we consider generative models for the probability of a functional connection (an edge) between two cortical regions (nodes) separated by some Euclidean distance in anatomical space. In particular, we propose a model in which the embedded topology of brain networks emerges from two competing factors: a distance penalty based on the cost of maintaining long-range connections; and a topological term that favors links between regions sharing similar input. We show that, together, these two biologically plausible factors are sufficient to capture an impressive range of topological properties of functional brain networks. Model parameters estimated in one set of functional MRI (fMRI) data on normal volunteers provided a good fit to networks estimated in a second independent sample of fMRI data. Furthermore, slightly detuned model parameters also generated a reasonable simulation of the abnormal properties of brain functional networks in people with schizophrenia. We therefore anticipate that many aspects of brain network organization, in health and disease, may be parsimoniously explained by an economical clustering rule for the probability of functional connectivity between different brain areas.

The Rich Club of the C. elegans Neuronal Connectome

There is increasing interest in topological analysis of brain networks as complex systems, with researchers often using neuroimaging to represent the large-scale organization of nervous systems without precise cellular resolution. Here we used graph theory to investigate the neuronal connectome of the nematode worm Caenorhabditis elegans, which is defined anatomically at a cellular scale as 2287 synaptic connections between 279 neurons. We identified a small number of highly connected neurons as a rich club (N = 11) interconnected with high efficiency and high connection distance. Rich club neurons comprise almost exclusively the interneurons of the locomotor circuits, with known functional importance for coordinated movement. The rich club neurons are connector hubs, with high betweenness centrality, and many intermodular connections to nodes in different modules. On identifying the shortest topological paths (motifs) between pairs of peripheral neurons, the motifs that are found most frequently traverse the rich club. The rich club neurons are born early in development, before visible movement of the animal and before the main phase of developmental elongation of its body. We conclude that the high wiring cost of the globally integrative rich club of neurons in the C. elegans connectome is justified by the adaptive value of coordinated movement of the animal. The economical trade-off between physical cost and behavioral value of rich club organization in a cellular connectome confirms theoretical expectations and recapitulates comparable results from human neuroimaging on much larger scale networks, suggesting that this may be a general and scale-invariant principle of brain network organization.

Hubs of brain functional networks are radically reorganized in comatose patients.

Human brain networks have topological properties in common with many other complex systems, prompting the following question: what aspects of brain network organization are critical for distinctive functional properties of the brain, such as consciousness? To address this question, we used graph theoretical methods to explore brain network topology in resting state functional MRI data acquired from 17 patients with severely impaired consciousness and 20 healthy volunteers. We found that many global network properties were conserved in comatose patients. Specifically, there was no significant abnormality of global efficiency, clustering, small-worldness, modularity, or degree distribution in the patient group. However, in every patient, we found evidence for a radical reorganization of high degree or highly efficient “hub” nodes. Cortical regions that were hubs of healthy brain networks had typically become nonhubs of comatose brain networks and vice versa. These results indicate that global topological properties of complex brain networks may be homeostatically conserved under extremely different clinical conditions and that consciousness likely depends on the anatomical location of hub nodes in human brain networks.

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.

Integrated strategy for improving functional connectivity mapping using multiecho fMRI

Functional connectivity analysis of resting state blood oxygen level–dependent (BOLD) functional MRI is widely used for noninvasively studying brain functional networks. Recent findings have indicated, however, that even small (≤1 mm) amounts of head movement during scanning can disproportionately bias connectivity estimates, despite various preprocessing efforts. Further complications for interregional connectivity estimation from time domain signals include the unaccounted reduction in BOLD degrees of freedom related to sensitivity losses from high subject motion. To address these issues, we describe an integrated strategy for data acquisition, denoising, and connectivity estimation. This strategy builds on our previously published technique combining data acquisition with multiecho (ME) echo planar imaging and analysis with spatial independent component analysis (ICA), called ME-ICA, which distinguishes BOLD (neuronal) and non-BOLD (artifactual) components based on linear echo-time dependence of signals—a characteristic property of BOLD signal changes. Here we show for 32 control subjects that this method provides a physically principled and nearly operator-independent way of removing complex artifacts such as motion from resting state data. We then describe a robust estimator of functional connectivity based on interregional correlation of BOLD-independent component coefficients. This estimator, called independent components regression, considerably simplifies statistical inference for functional connectivity because degrees of freedom equals the number of independent coefficients. Compared with traditional connectivity estimation methods, the proposed strategy results in fourfold improvements in signal-to-noise ratio, functional connectivity analysis with improved specificity, and valid statistical inference with nominal control of type 1 error in contrasts of connectivity between groups with different levels of subject motion.

Annual Research Review: Growth connectomics – the organization and reorganization of brain networks during normal and abnormal development

Background We first give a brief introduction to graph theoretical analysis and its application to the study of brain network topology or connectomics. Within this framework, we review the existing empirical data on developmental changes in brain network organization across a range of experimental modalities (including structural and functional MRI, diffusion tensor imaging, magnetoencephalography and electroencephalography in humans). Synthesis We discuss preliminary evidence and current hypotheses for how the emergence of network properties correlates with concomitant cognitive and behavioural changes associated with development. We highlight some of the technical and conceptual challenges to be addressed by future developments in this rapidly moving field. Given the parallels previously discovered between neural systems across species and over a range of spatial scales, we also review some recent advances in developmental network studies at the cellular scale. We highlight the opportunities presented by such studies and how they may complement neuroimaging in advancing our understanding of brain development. Finally, we note that many brain and mind disorders are thought to be neurodevelopmental in origin and that charting the trajectory of brain network changes associated with healthy development also sets the stage for understanding abnormal network development. Conclusions We therefore briefly review the clinical relevance of network metrics as potential diagnostic markers and some recent efforts in computational modelling of brain networks which might contribute to a more mechanistic understanding of neurodevelopmental disorders in future.

Adolescence is associated with genomically patterned consolidation of the hubs of the human brain connectome

Significance Adolescence is a period of human brain growth and high incidence of mental health disorders. Here, we show consistently in two MRI cohorts that human brain changes in adolescence were concentrated on the more densely connected hubs of the connectome (i.e., association cortical regions that mediated efficient connectivity throughout the human brain structural network). Hubs were less myelinated at 14 y but had faster rates of myelination and cortical shrinkage in the 14- to 24-y period. This topologically focused process of cortical consolidation was associated with expression of genes enriched for normal synaptic and myelin-related processes and risk of schizophrenia. Consolidation of anatomical network hubs could be important for normal and clinically disordered adolescent brain development. How does human brain structure mature during adolescence? We used MRI to measure cortical thickness and intracortical myelination in 297 population volunteers aged 14–24 y old. We found and replicated that association cortical areas were thicker and less myelinated than primary cortical areas at 14 y. However, association cortex had faster rates of shrinkage and myelination over the course of adolescence. Age-related increases in cortical myelination were maximized approximately at the internal layer of projection neurons. Adolescent cortical myelination and shrinkage were coupled and specifically associated with a dorsoventrally patterned gene expression profile enriched for synaptic, oligodendroglial- and schizophrenia-related genes. Topologically efficient and biologically expensive hubs of the brain anatomical network had greater rates of shrinkage/myelination and were associated with overexpression of the same transcriptional profile as cortical consolidation. We conclude that normative human brain maturation involves a genetically patterned process of consolidating anatomical network hubs. We argue that developmental variation of this consolidation process may be relevant both to normal cognitive and behavioral changes and the high incidence of schizophrenia during human brain adolescence.

Volitional eyes opening perturbs brain dynamics and functional connectivity regardless of light input

The act of opening (or closing) ones eyes has long been demonstrated to impact on brain function. However, the eyes open condition is usually accompanied by visual input, and this effect may have been a significant confounding factor in previous studies. To clarify this situation, we extended the traditional eyes open/closed study to a two-factor balanced, repeated measures resting state fMRI (rs-fMRI) experiment, in which light on/off was also included as a factor. In 16 healthy participants, we estimated the univariate properties of the BOLD signal, as well as a bivariate measure of functional connectivity and multivariate network topology measures. Across all these measures, we demonstrate that human brain adopts a distinctive configuration when eyes are open (compared to when eyes are closed) independently of exogenous light input: (i) the eyes open states were associated with decreased BOLD signal variance (P-value=0.0004), decreased fractional amplitude of low frequency fluctuation (fALFF. P-value=0.0061), and decreased Hurst exponent (H. P-value=0.0321) mainly in the primary and secondary sensory cortical areas, the insula, and the thalamus. (ii) The strength of functional connectivity (FC) between the posterior cingulate cortex (PCC), a major component of the default mode network (DMN), and the bilateral perisylvian and perirolandic regions was also significantly decreased during eyes open states. (iii) On the other hand, the average network connection distance increased during eyes open states (P-value=0.0139). Additionally, the metrics of univariate, bivariate, and multivariate analyses in this study are significantly correlated. In short, we have shown that the marked effects on the dynamics and connectivity of fMRI time series brought by volitional eyes open or closed are simply endogenous and irrespective of exogenous visual stimulus. The state of eyes open (or closed) may thus be an important factor to control in design of rs-fMRI and even other cognitive block or event-related experiments.

The Multilayer Connectome of Caenorhabditis elegans

Connectomics has focused primarily on the mapping of synaptic links in the brain; yet it is well established that extrasynaptic volume transmission, especially via monoamines and neuropeptides, is also critical to brain function and occurs primarily outside the synaptic connectome. We have mapped the putative monoamine connections, as well as a subset of neuropeptide connections, in C. elegans based on new and published gene expression data. The monoamine and neuropeptide networks exhibit distinct topological properties, with the monoamine network displaying a highly disassortative star-like structure with a rich-club of interconnected broadcasting hubs, and the neuropeptide network showing a more recurrent, highly clustered topology. Despite the low degree of overlap between the extrasynaptic (or wireless) and synaptic (or wired) connectomes, we find highly significant multilink motifs of interaction, pinpointing locations in the network where aminergic and neuropeptide signalling modulate synaptic activity. Thus, the C. elegans connectome can be mapped as a multiplex network with synaptic, gap junction, and neuromodulator layers representing alternative modes of interaction between neurons. This provides a new topological plan for understanding how aminergic and peptidergic modulation of behaviour is achieved by specific motifs and loci of integration between hard-wired synaptic or junctional circuits and extrasynaptic signals wirelessly broadcast from a small number of modulatory neurons.