The wiring economy principle: connectivity determines anatomy in the human brain

Minimization of the wiring cost of white matter fibers in the human brain appears to be an organizational principle. We investigate this aspect in the human brain using whole brain connectivity networks extracted from high resolution diffusion MRI data of 14 normal volunteers. We specifically address the question of whether brain anatomy determines its connectivity or vice versa. Unlike previous studies we use weighted networks, where connections between cortical nodes are real-valued rather than binary off-on connections. In one set of analyses we found that the connectivity structure of the brain has near optimal wiring cost compared to random networks with the same number of edges, degree distribution and edge weight distribution. A specifically designed minimization routine could not find cheaper wiring without significantly degrading network performance. In another set of analyses we kept the observed brain network topology and connectivity but allowed nodes to freely move on a 3D manifold topologically identical to the brain. An efficient minimization routine was written to find the lowest wiring cost configuration. We found that beginning from any random configuration, the nodes invariably arrange themselves in a configuration with a striking resemblance to the brain. This confirms the widely held but poorly tested claim that wiring economy is a driving principle of the brain. Intriguingly, our results also suggest that the brain mainly optimizes for the most desirable network connectivity, and the observed brain anatomy is merely a result of this optimization.     

Dynamics of excitable neural networks with heterogeneous connectivity

A central issue of neuroscience is to understand how neural units integrates internal and external signals to create coherent states. Recently, it has been shown that the sensitivity and dynamic range of neural assemblies are optimal at a critical coupling among its elements. Complex architectures of connections seem to play a constructive role on the reliable coordination of neural units. Here we show that, the synchronizability and sensitivity of excitable neural networks can be tuned by diversity in the connections strengths. We illustrate our findings for weighted networks with regular, random and complex topologies. Additional comparisons of real brain networks support previous studies suggesting that heterogeneity in the connectivity may play a constructive role on information processing. These findings provide insights into the relationship between structure and function of neural circuits.     

Analysis of Whole-Brain Resting-State fMRI Data Using Hierarchical Clustering Approach

Background

Previous studies using hierarchical clustering approach to analyze resting-state fMRI data were limited to a few slices or regions-of-interest (ROIs) after substantial data reduction.

Purpose

To develop a framework that can perform voxel-wise hierarchical clustering of whole-brain resting-state fMRI data from a group of subjects.

Materials and Methods

Resting-state fMRI measurements were conducted for 86 adult subjects using a single-shot echo-planar imaging (EPI) technique. After pre-processing and co-registration to a standard template, pair-wise cross-correlation coefficients (CC) were calculated for all voxels inside the brain and translated into absolute Pearson''s distances after imposing a threshold CC≥0.3. The group averages of the Pearson''s distances were then used to perform hierarchical clustering with the developed framework, which entails gray matter masking and an iterative scheme to analyze the dendrogram.

Results

With the hierarchical clustering framework, we identified most of the functional connectivity networks reported previously in the literature, such as the motor, sensory, visual, memory, and the default-mode functional networks (DMN). Furthermore, the DMN and visual system were split into their corresponding hierarchical sub-networks.

Conclusion

It is feasible to use the proposed hierarchical clustering scheme for voxel-wise analysis of whole-brain resting-state fMRI data. The hierarchical clustering result not only confirmed generally the finding in functional connectivity networks identified previously using other data processing techniques, such as ICA, but also revealed directly the hierarchical structure within the functional connectivity networks.     

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.     

An Adaptive Complex Network Model for Brain Functional Networks

Brain functional networks are graph representations of activity in the brain, where the vertices represent anatomical regions and the edges their functional connectivity. These networks present a robust small world topological structure, characterized by highly integrated modules connected sparsely by long range links. Recent studies showed that other topological properties such as the degree distribution and the presence (or absence) of a hierarchical structure are not robust, and show different intriguing behaviors. In order to understand the basic ingredients necessary for the emergence of these complex network structures we present an adaptive complex network model for human brain functional networks. The microscopic units of the model are dynamical nodes that represent active regions of the brain, whose interaction gives rise to complex network structures. The links between the nodes are chosen following an adaptive algorithm that establishes connections between dynamical elements with similar internal states. We show that the model is able to describe topological characteristics of human brain networks obtained from functional magnetic resonance imaging studies. In particular, when the dynamical rules of the model allow for integrated processing over the entire network scale-free non-hierarchical networks with well defined communities emerge. On the other hand, when the dynamical rules restrict the information to a local neighborhood, communities cluster together into larger ones, giving rise to a hierarchical structure, with a truncated power law degree distribution.     

Graph theoretical analysis of functional brain networks: test-retest evaluation on short- and long-term resting-state functional MRI data

Graph-based computational network analysis has proven a powerful tool to quantitatively characterize functional architectures of the brain. However, the test-retest (TRT) reliability of graph metrics of functional networks has not been systematically examined. Here, we investigated TRT reliability of topological metrics of functional brain networks derived from resting-state functional magnetic resonance imaging data. Specifically, we evaluated both short-term (<1 hour apart) and long-term (>5 months apart) TRT reliability for 12 global and 6 local nodal network metrics. We found that reliability of global network metrics was overall low, threshold-sensitive and dependent on several factors of scanning time interval (TI, long-term>short-term), network membership (NM, networks excluding negative correlations>networks including negative correlations) and network type (NT, binarized networks>weighted networks). The dependence was modulated by another factor of node definition (ND) strategy. The local nodal reliability exhibited large variability across nodal metrics and a spatially heterogeneous distribution. Nodal degree was the most reliable metric and varied the least across the factors above. Hub regions in association and limbic/paralimbic cortices showed moderate TRT reliability. Importantly, nodal reliability was robust to above-mentioned four factors. Simulation analysis revealed that global network metrics were extremely sensitive (but varying degrees) to noise in functional connectivity and weighted networks generated numerically more reliable results in compared with binarized networks. For nodal network metrics, they showed high resistance to noise in functional connectivity and no NT related differences were found in the resistance. These findings provide important implications on how to choose reliable analytical schemes and network metrics of interest.     

Scale-free behavior in protein domain networks

Several technical, social, and biological networks were recently found to demonstrate scale-free and small-world behavior instead of random graph characteristics. In this work, the topology of protein domain networks generated with data from the ProDom, Pfam, and Prosite domain databases was studied. It was found that these networks exhibited small-world and scale-free topologies with a high degree of local clustering accompanied by a few long-distance connections. Moreover, these observations apply not only to the complete databases, but also to the domain distributions in proteomes of different organisms. The extent of connectivity among domains reflects the evolutionary complexity of the organisms considered.     

Consensus between Pipelines in Structural Brain Networks

Structural brain networks may be reconstructed from diffusion MRI tractography data and have great potential to further our understanding of the topological organisation of brain structure in health and disease. Network reconstruction is complex and involves a series of processesing methods including anatomical parcellation, registration, fiber orientation estimation and whole-brain fiber tractography. Methodological choices at each stage can affect the anatomical accuracy and graph theoretical properties of the reconstructed networks, meaning applying different combinations in a network reconstruction pipeline may produce substantially different networks. Furthermore, the choice of which connections are considered important is unclear. In this study, we assessed the similarity between structural networks obtained using two independent state-of-the-art reconstruction pipelines. We aimed to quantify network similarity and identify the core connections emerging most robustly in both pipelines. Similarity of network connections was compared between pipelines employing different atlases by merging parcels to a common and equivalent node scale. We found a high agreement between the networks across a range of fiber density thresholds. In addition, we identified a robust core of highly connected regions coinciding with a peak in similarity across network density thresholds, and replicated these results with atlases at different node scales. The binary network properties of these core connections were similar between pipelines but showed some differences in atlases across node scales. This study demonstrates the utility of applying multiple structural network reconstrution pipelines to diffusion data in order to identify the most important connections for further study.     

Hierarchical Alteration of Brain Structural and Functional Networks in Female Migraine Sufferers

Background

Little is known about the changes of brain structural and functional connectivity networks underlying the pathophysiology in migraine. We aimed to investigate how the cortical network reorganization is altered by frequent cortical overstimulation associated with migraine.

Methodology/Principal Findings

Gray matter volumes and resting-state functional magnetic resonance imaging signal correlations were employed to construct structural and functional networks between brain regions in 43 female patients with migraine (PM) and 43 gender-matched healthy controls (HC) by using graph theory-based approaches. Compared with the HC group, the patients showed abnormal global topology in both structural and functional networks, characterized by higher mean clustering coefficients without significant change in the shortest absolute path length, which indicated that the PM lost optimal topological organization in their cortical networks. Brain hubs related to pain-processing revealed abnormal nodal centrality in both structural and functional networks, including the precentral gyrus, orbital part of the inferior frontal gyrus, parahippocampal gyrus, anterior cingulate gyrus, thalamus, temporal pole of the middle temporal gyrus and the inferior parietal gyrus. Negative correlations were found between migraine duration and regions with abnormal centrality. Furthermore, the dysfunctional connections in patients'' cortical networks formed into a connected component and three dysregulated modules were identified involving pain-related information processing and motion-processing visual networks.

Conclusions

Our results may reflect brain alteration dynamics resulting from migraine and suggest that long-term and high-frequency headache attacks may cause both structural and functional connectivity network reorganization. The disrupted information exchange between brain areas in migraine may be reshaped into a hierarchical modular structure progressively.     

A novel index of functional connectivity: phase lag based on Wilcoxon signed rank test

Phase synchronization has been an effective measurement of functional connectivity, detecting similar dynamics over time among distinct brain regions. However, traditional phase synchronization-based functional connectivity indices have been proved to have some drawbacks. For example, the phase locking value (PLV) index is sensitive to volume conduction, while the phase lag index (PLI) and the weighted phase lag index (wPLI) are easily affected by noise perturbations. In addition, thresholds need to be applied to these indices to obtain the binary adjacency matrix that determines the connections. However, the selection of the thresholds is generally arbitrary. To address these issues, in this paper we propose a novel index of functional connectivity, named the phase lag based on the Wilcoxon signed-rank test (PLWT). Specifically, it characterizes the functional connectivity based on the phase lag with a weighting procedure to reduce the influence of volume conduction and noise. Besides, it automatically identifies the important connections without relying on thresholds, by taking advantage of the framework of the Wilcoxon signed-rank test. The performance of the proposed PLWT index is evaluated on simulated electroencephalograph (EEG) datasets, as well as on two resting-state EEG datasets. The experimental results on the simulated EEG data show that the PLWT index is robust to volume conduction and noise. Furthermore, the brain functional networks derived by PLWT on the real EEG data exhibit a reasonable scale-free characteristic and high test–retest (TRT) reliability of graph measures. We believe that the proposed PLWT index provides a useful and reliable tool to identify the underlying neural interactions, while effectively diminishing the influence of volume conduction and noise.     

Resting-state quantitative electroencephalography reveals increased neurophysiologic connectivity in depression

Symptoms of Major Depressive Disorder (MDD) are hypothesized to arise from dysfunction in brain networks linking the limbic system and cortical regions. Alterations in brain functional cortical connectivity in resting-state networks have been detected with functional imaging techniques, but neurophysiologic connectivity measures have not been systematically examined. We used weighted network analysis to examine resting state functional connectivity as measured by quantitative electroencephalographic (qEEG) coherence in 121 unmedicated subjects with MDD and 37 healthy controls. Subjects with MDD had significantly higher overall coherence as compared to controls in the delta (0.5-4 Hz), theta (4-8 Hz), alpha (8-12 Hz), and beta (12-20 Hz) frequency bands. The frontopolar region contained the greatest number of "hub nodes" (surface recording locations) with high connectivity. MDD subjects expressed higher theta and alpha coherence primarily in longer distance connections between frontopolar and temporal or parietooccipital regions, and higher beta coherence primarily in connections within and between electrodes overlying the dorsolateral prefrontal cortical (DLPFC) or temporal regions. Nearest centroid analysis indicated that MDD subjects were best characterized by six alpha band connections primarily involving the prefrontal region. The present findings indicate a loss of selectivity in resting functional connectivity in MDD. The overall greater coherence observed in depressed subjects establishes a new context for the interpretation of previous studies showing differences in frontal alpha power and synchrony between subjects with MDD and normal controls. These results can inform the development of qEEG state and trait biomarkers for MDD.     

Effective Connectivity Extracted from Resting-State fMRI Images Using Transfer Entropy

Background and objectiveBased on magnetic resonance imaging (MRI), macroscopic structural and functional connectivity of human brain has been widely explored in the last decade. However, little work has been done on effective connectivity between individual brain parcels. In this preliminary study, we aim to investigate whole-brain effective connectivity networks from resting-state functional MRI (rs-fMRI) images.Material and methodsAfter the functional connectivity networks of 26 healthy subjects (aged from 25 to 35 years old) from Human Connectome Project database were derived from rs-fMRI images with dynamic time warping, proportional thresholding (PT) was performed on the functional connectivity matrices by retaining the PT% strongest functional connections. PT% ranges from 40% to 10% in steps of 5%. Then, effective connections corresponding to the PT% strongest functional connections, both bi-directional and unidirectional, were estimated with Renyi's 2-order transfer entropy (TE) method. Topological metrics of the built functional and effective connectivity networks were further characterized, including clustering coefficient, transitivity, and modularity.ResultsIt is found that the effective connectivity networks exhibit small world attributes, and that the networks contain a subset of highly interactive regions, including right frontal pole (in-degree 6), left middle frontal gyrus (in-degree 8, out-degree 1), right precentral gyrus (out-degree 9), left precentral gyrus (out-degree 7), right posterior division of supramarginal gyrus (in-degree 2, out-degree 3), left angular gyrus (out-degree 6), left inferior division of lateral occipital cortex (out-degree 6), right occipital pole (in-degree 5), right cerebellum 7b parcel (in-degree 15), and right cerebellum 8 parcel (in-degree 7, out-degree 1).ConclusionsThe observations in this study provide information about the casual interactions among brain parcels in resting state, helping reveal how different subregions of large-scale distributed neural networks are coupled together in performing cognitive functions.     

Resolving Anatomical and Functional Structure in Human Brain Organization: Identifying Mesoscale Organization in Weighted Network Representations

Human brain anatomy and function display a combination of modular and hierarchical organization, suggesting the importance of both cohesive structures and variable resolutions in the facilitation of healthy cognitive processes. However, tools to simultaneously probe these features of brain architecture require further development. We propose and apply a set of methods to extract cohesive structures in network representations of brain connectivity using multi-resolution techniques. We employ a combination of soft thresholding, windowed thresholding, and resolution in community detection, that enable us to identify and isolate structures associated with different weights. One such mesoscale structure is bipartivity, which quantifies the extent to which the brain is divided into two partitions with high connectivity between partitions and low connectivity within partitions. A second, complementary mesoscale structure is modularity, which quantifies the extent to which the brain is divided into multiple communities with strong connectivity within each community and weak connectivity between communities. Our methods lead to multi-resolution curves of these network diagnostics over a range of spatial, geometric, and structural scales. For statistical comparison, we contrast our results with those obtained for several benchmark null models. Our work demonstrates that multi-resolution diagnostic curves capture complex organizational profiles in weighted graphs. We apply these methods to the identification of resolution-specific characteristics of healthy weighted graph architecture and altered connectivity profiles in psychiatric disease.     

Disrupted functional brain connectivity in partial epilepsy: a resting-state fMRI study

Examining the spontaneous activity to understand the neural mechanism of brain disorder is a focus in recent resting-state fMRI. In the current study, to investigate the alteration of brain functional connectivity in partial epilepsy in a systematical way, two levels of analyses (functional connectivity analysis within resting state networks (RSNs) and functional network connectivity (FNC) analysis) were carried out on resting-state fMRI data acquired from the 30 participants including 14 healthy controls(HC) and 16 partial epilepsy patients. According to the etiology, all patients are subdivided into temporal lobe epilepsy group (TLE, included 7 patients) and mixed partial epilepsy group (MPE, 9 patients). Using group independent component analysis, eight RSNs were identified, and selected to evaluate functional connectivity and FNC between groups. Compared with the controls, decreased functional connectivity within all RSNs was found in both TLE and MPE. However, dissociating patterns were observed within the 8 RSNs between two patient groups, i.e, compared with TLE, we found decreased functional connectivity in 5 RSNs increased functional connectivity in 1 RSN, and no difference in the other 2 RSNs in MPE. Furthermore, the hierarchical disconnections of FNC was found in two patient groups, in which the intra-system connections were preserved for all three subsystems while the lost connections were confined to intersystem connections in patients with partial epilepsy. These findings may suggest that decreased resting state functional connectivity and disconnection of FNC are two remarkable characteristics of partial epilepsy. The selective impairment of FNC implicated that it is unsuitable to understand the partial epilepsy only from global or local perspective. We presumed that studying epilepsy in the multi-perspective based on RSNs may be a valuable means to assess the functional changes corresponding to specific RSN and may contribute to the understanding of the neuro-pathophysiological mechanism of epilepsy.     

Preserved Modular Network Organization in the Sedated Rat Brain

Translation of resting-state functional connectivity (FC) magnetic resonance imaging (rs-fMRI) applications from human to rodents has experienced growing interest, and bears a great potential in pre-clinical imaging as it enables assessing non-invasively the topological organization of complex FC networks (FCNs) in rodent models under normal and various pathophysiological conditions. However, to date, little is known about the organizational architecture of FCNs in rodents in a mentally healthy state, although an understanding of the same is of paramount importance before investigating networks under compromised states. In this study, we characterized the properties of resting-state FCN in an extensive number of Sprague-Dawley rats (n = 40) under medetomidine sedation by evaluating its modular organization and centrality of brain regions and tested for reproducibility. Fully-connected large-scale complex networks of positively and negatively weighted connections were constructed based on Pearson partial correlation analysis between the time courses of 36 brain regions encompassing almost the entire brain. Applying recently proposed complex network analysis measures, we show that the rat FCN exhibits a modular architecture, comprising six modules with a high between subject reproducibility. In addition, we identified network hubs with strong connections to diverse brain regions. Overall our results obtained under a straight medetomidine protocol show for the first time that the community structure of the rat brain is preserved under pharmacologically induced sedation with a network modularity contrasting from the one reported for deep anesthesia but closely resembles the organization described for the rat in conscious state.     

Network analysis of intrinsic functional brain connectivity in Alzheimer's disease

Functional brain networks detected in task-free (“resting-state”) functional magnetic resonance imaging (fMRI) have a small-world architecture that reflects a robust functional organization of the brain. Here, we examined whether this functional organization is disrupted in Alzheimer's disease (AD). Task-free fMRI data from 21 AD subjects and 18 age-matched controls were obtained. Wavelet analysis was applied to the fMRI data to compute frequency-dependent correlation matrices. Correlation matrices were thresholded to create 90-node undirected-graphs of functional brain networks. Small-world metrics (characteristic path length and clustering coefficient) were computed using graph analytical methods. In the low frequency interval 0.01 to 0.05 Hz, functional brain networks in controls showed small-world organization of brain activity, characterized by a high clustering coefficient and a low characteristic path length. In contrast, functional brain networks in AD showed loss of small-world properties, characterized by a significantly lower clustering coefficient (p<0.01), indicative of disrupted local connectivity. Clustering coefficients for the left and right hippocampus were significantly lower (p<0.01) in the AD group compared to the control group. Furthermore, the clustering coefficient distinguished AD participants from the controls with a sensitivity of 72% and specificity of 78%. Our study provides new evidence that there is disrupted organization of functional brain networks in AD. Small-world metrics can characterize the functional organization of the brain in AD, and our findings further suggest that these network measures may be useful as an imaging-based biomarker to distinguish AD from healthy aging.     

Comparing structural and functional graph theory features in the human brain using multimodal MRI

Recent advances in magnetic resonance imaging (MRI) are allowing neuroscientists to gain critical insights into the neural networks mediating a variety of cognitive processes. This work investigates structural and functional connectivity in the human brain under different experimental conditions through multimodal MRI acquisitions. To define the nodes of a full-brain network, a set of regions was identified from resting-state functional MRI (fMRI) data using spatial independent component analysis (sICA) and a hierarchical clustering technique. Diffusion-weighted imaging (DWI) data were acquired from the same subjects and a probabilistic fiber tracking method was used to estimate the structure of this network. Using features originating from graph theory, such as small-world properties and network efficiency, we characterized the structural and functional connectivities of the full-brain network and we compared them quantitatively. We showed that structural and functional networks shared some properties in terms of topology as measured by the distribution of the node degrees, hence supporting the existence of an underlying anatomical substrate for functional networks.     

Altered Network Communication Following a Neuroprotective Drug Treatment

Preconditioning is defined as a range of stimuli that allow cells to withstand subsequent anaerobic and other deleterious conditions. While cell protection under preconditioning is well established, this paper investigates the influence of neuroprotective preconditioning drugs, 4-aminopyridine and bicuculline (4-AP/bic), on synaptic communication across a broad network of in vitro rat cortical neurons. Using a permutation test, we evaluated cross-correlations of extracellular spiking activity across all pairs of recording electrodes on a 64-channel multielectrode array. The resulting functional connectivity maps were analyzed in terms of their graph-theoretic properties. A small-world effect was found, characterized by a functional network with high clustering coefficient and short average path length. Twenty-four hours after exposure to 4-AP/bic, small-world properties were comparable to control cultures that were not treated with the drug. Four hours following drug washout, however, the density of functional connections increased, while path length decreased and clustering coefficient increased. These alterations in functional connectivity were maintained at four days post-washout, suggesting that 4-AP/bic preconditioning leads to long-term effects on functional networks of cortical neurons. Because of their influence on communication efficiency in neuronal networks, alterations in small-world properties hold implications for information processing in brain systems. The observed relationship between density, path length, and clustering coefficient is captured by a phenomenological model where connections are added randomly within a spatially-embedded network. Taken together, results provide information regarding functional consequences of drug therapies that are overlooked in traditional viability studies and present the first investigation of functional networks under neuroprotective preconditioning.