Altered Network Topologies and Hub Organization in Adults with Autism: A Resting-State fMRI Study

Recent functional magnetic resonance imaging (fMRI) studies on autism spectrum condition (ASC) have identified dysfunctions in specific brain networks involved in social and non-social cognition that persist into adulthood. Although increasing numbers of fMRI studies have revealed atypical functional connectivity in the adult ASC brain, such functional alterations at the network level have not yet been fully characterized within the recently developed graph-theoretical framework. Here, we applied a graph-theoretical analysis to resting-state fMRI data acquired from 46 adults with ASC and 46 age- and gender-matched controls, to investigate the topological properties and organization of autistic brain network. Analyses of global metrics revealed that, relative to the controls, participants with ASC exhibited significant decreases in clustering coefficient and characteristic path length, indicating a shift towards randomized organization. Furthermore, analyses of local metrics revealed a significantly altered organization of the hub nodes in ASC, as shown by analyses of hub disruption indices using multiple local metrics and by a loss of “hubness” in several nodes (e.g., the bilateral superior temporal sulcus, right dorsolateral prefrontal cortex, and precuneus) that are critical for social and non-social cognitive functions. In particular, local metrics of the anterior cingulate cortex consistently showed significant negative correlations with the Autism-Spectrum Quotient score. Our results demonstrate altered patterns of global and local topological properties that may underlie impaired social and non-social cognition in ASC.     

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.     

Epidemic threshold and network structure: The interplay of probability of transmission and of persistence in small-size directed networks

A growing number of studies are investigating the effect of contact structure on the dynamics of epidemics in large-scale complex networks. Whether findings thus obtained apply also to networks of small size, and thus to many real-world biological applications, is still an open question. We use numerical simulations of disease spread in directed networks of 100 individual nodes with a constant number of links. We show that, no matter the type of network structure (local, small-world, random and scale-free), there is a linear threshold determined by the probability of infection transmission between connected nodes and the probability of infection persistence in an infected node. The threshold is significantly lower for scale-free networks compared to local, random and small-world ones only if super-connected nodes have a higher number of links both to and from other nodes. The starting point, the node at which the epidemic starts, does not affect the threshold conditions, but has a marked influence on the final size of the epidemic in all kinds of network. There is evidence that contact structure has an influence on the average final size of an epidemic across all starting nodes, with significantly lower values in scale-free networks at equilibrium. Simulations in scale-free networks show a distinctive time-series pattern, which, if found in a real epidemic, can be used to infer the underlying network structure. The findings have relevance also for meta-population ecology and species conservation.     

fMRI Resting Slow Fluctuations Correlate with the Activity of Fast Cortico-Cortical Physiological Connections

Recording of slow spontaneous fluctuations at rest using functional magnetic resonance imaging (fMRI) allows distinct long-range cortical networks to be identified. The neuronal basis of connectivity as assessed by resting-state fMRI still needs to be fully clarified, considering that these signals are an indirect measure of neuronal activity, reflecting slow local variations in de-oxyhaemoglobin concentration. Here, we combined fMRI with multifocal transcranial magnetic stimulation (TMS), a technique that allows the investigation of the causal neurophysiological interactions occurring in specific cortico-cortical connections. We investigated whether the physiological properties of parieto-frontal circuits mapped with short-latency multifocal TMS at rest may have some relationship with the resting-state fMRI measures of specific resting-state functional networks (RSNs). Results showed that the activity of fast cortico-cortical physiological interactions occurring in the millisecond range correlated selectively with the coupling of fMRI slow oscillations within the same cortical areas that form part of the dorsal attention network, i.e., the attention system believed to be involved in reorientation of attention. We conclude that resting-state fMRI ongoing slow fluctuations likely reflect the interaction of underlying physiological cortico-cortical connections.     

Topological fractionation of resting-state networks

Exploring topological properties of human brain network has become an exciting topic in neuroscience research. Large-scale structural and functional brain networks both exhibit a small-world topology, which is evidence for global and local parallel information processing. Meanwhile, resting state networks (RSNs) underlying specific biological functions have provided insights into how intrinsic functional architecture influences cognitive and perceptual information processing. However, topological properties of single RSNs remain poorly understood. Here, we have two hypotheses: i) each RSN also has optimized small-world architecture; ii) topological properties of RSNs related to perceptual and higher cognitive processes are different. To test these hypotheses, we investigated the topological properties of the default-mode, dorsal attention, central-executive, somato-motor, visual and auditory networks derived from resting-state functional magnetic resonance imaging (fMRI). We found small-world topology in each RSN. Furthermore, small-world properties of cognitive networks were higher than those of perceptual networks. Our findings are the first to demonstrate a topological fractionation between perceptual and higher cognitive networks. Our approach may be useful for clinical research, especially for diseases that show selective abnormal connectivity in specific brain networks.     

Disease spread in small-size directed networks: Epidemic threshold, correlation between links to and from nodes, and clustering

Network epidemiology has mainly focused on large-scale complex networks. It is unclear whether findings of these investigations also apply to networks of small size. This knowledge gap is of relevance for many biological applications, including meta-communities, plant–pollinator interactions and the spread of the oomycete pathogen Phytophthora ramorum in networks of plant nurseries. Moreover, many small-size biological networks are inherently asymmetrical and thus cannot be realistically modelled with undirected networks. We modelled disease spread and establishment in directed networks of 100 and 500 nodes at four levels of connectance in six network structures (local, small-world, random, one-way, uncorrelated, and two-way scale-free networks). The model was based on the probability of infection persistence in a node and of infection transmission between connected nodes. Regardless of the size of the network, the epidemic threshold did not depend on the starting node of infection but was negatively related to the correlation coefficient between in- and out-degree for all structures, unless networks were sparsely connected. In this case clustering played a significant role. For small-size scale-free directed networks to have a lower epidemic threshold than other network structures, there needs to be a positive correlation between number of links to and from nodes. When this correlation is negative (one-way scale-free networks), the epidemic threshold for small-size networks can be higher than in non-scale-free networks. Clustering does not necessarily have an influence on the epidemic threshold if connectance is kept constant. Analyses of the influence of the clustering on the epidemic threshold in directed networks can also be spurious if they do not consider simultaneously the effect of the correlation coefficient between in- and out-degree.     

Scalability of Asynchronous Networks Is Limited by One-to-One Mapping between Effective Connectivity and Correlations

Network models are routinely downscaled compared to nature in terms of numbers of nodes or edges because of a lack of computational resources, often without explicit mention of the limitations this entails. While reliable methods have long existed to adjust parameters such that the first-order statistics of network dynamics are conserved, here we show that limitations already arise if also second-order statistics are to be maintained. The temporal structure of pairwise averaged correlations in the activity of recurrent networks is determined by the effective population-level connectivity. We first show that in general the converse is also true and explicitly mention degenerate cases when this one-to-one relationship does not hold. The one-to-one correspondence between effective connectivity and the temporal structure of pairwise averaged correlations implies that network scalings should preserve the effective connectivity if pairwise averaged correlations are to be held constant. Changes in effective connectivity can even push a network from a linearly stable to an unstable, oscillatory regime and vice versa. On this basis, we derive conditions for the preservation of both mean population-averaged activities and pairwise averaged correlations under a change in numbers of neurons or synapses in the asynchronous regime typical of cortical networks. We find that mean activities and correlation structure can be maintained by an appropriate scaling of the synaptic weights, but only over a range of numbers of synapses that is limited by the variance of external inputs to the network. Our results therefore show that the reducibility of asynchronous networks is fundamentally limited.     

Influence of Choice of Null Network on Small-World Parameters of Structural Correlation Networks

In recent years, coordinated variations in brain morphology (e.g., volume, thickness) have been employed as a measure of structural association between brain regions to infer large-scale structural correlation networks. Recent evidence suggests that brain networks constructed in this manner are inherently more clustered than random networks of the same size and degree. Thus, null networks constructed by randomizing topology are not a good choice for benchmarking small-world parameters of these networks. In the present report, we investigated the influence of choice of null networks on small-world parameters of gray matter correlation networks in healthy individuals and survivors of acute lymphoblastic leukemia. Three types of null networks were studied: 1) networks constructed by topology randomization (TOP), 2) networks matched to the distributional properties of the observed covariance matrix (HQS), and 3) networks generated from correlation of randomized input data (COR). The results revealed that the choice of null network not only influences the estimated small-world parameters, it also influences the results of between-group differences in small-world parameters. In addition, at higher network densities, the choice of null network influences the direction of group differences in network measures. Our data suggest that the choice of null network is quite crucial for interpretation of group differences in small-world parameters of structural correlation networks. We argue that none of the available null models is perfect for estimation of small-world parameters for correlation networks and the relative strengths and weaknesses of the selected model should be carefully considered with respect to obtained network measures.     

Altered topological properties of functional network connectivity in schizophrenia during resting state: a small-world brain network study

Aberrant topological properties of small-world human brain networks in patients with schizophrenia (SZ) have been documented in previous neuroimaging studies. Aberrant functional network connectivity (FNC, temporal relationships among independent component time courses) has also been found in SZ by a previous resting state functional magnetic resonance imaging (fMRI) study. However, no study has yet determined if topological properties of FNC are also altered in SZ. In this study, small-world network metrics of FNC during the resting state were examined in both healthy controls (HCs) and SZ subjects. FMRI data were obtained from 19 HCs and 19 SZ. Brain images were decomposed into independent components (ICs) by group independent component analysis (ICA). FNC maps were constructed via a partial correlation analysis of ICA time courses. A set of undirected graphs were built by thresholding the FNC maps and the small-world network metrics of these maps were evaluated. Our results demonstrated significantly altered topological properties of FNC in SZ relative to controls. In addition, topological measures of many ICs involving frontal, parietal, occipital and cerebellar areas were altered in SZ relative to controls. Specifically, topological measures of whole network and specific components in SZ were correlated with scores on the negative symptom scale of the Positive and Negative Symptom Scale (PANSS). These findings suggest that aberrant architecture of small-world brain topology in SZ consists of ICA temporally coherent brain networks.     

Functional Connectivity in MRI Is Driven by Spontaneous BOLD Events

Functional brain signals are frequently decomposed into a relatively small set of large scale, distributed cortical networks that are associated with different cognitive functions. It is generally assumed that the connectivity of these networks is static in time and constant over the whole network, although there is increasing evidence that this view is too simplistic. This work proposes novel techniques to investigate the contribution of spontaneous BOLD events to the temporal dynamics of functional connectivity as assessed by ultra-high field functional magnetic resonance imaging (fMRI). The results show that: 1) spontaneous events in recognised brain networks contribute significantly to network connectivity estimates; 2) these spontaneous events do not necessarily involve whole networks or nodes, but clusters of voxels which act in concert, forming transiently synchronising sub-networks and 3) a task can significantly alter the number of localised spontaneous events that are detected within a single network. These findings support the notion that spontaneous events are the main driver of the large scale networks that are commonly detected by seed-based correlation and ICA. Furthermore, we found that large scale networks are manifestations of smaller, transiently synchronising sub-networks acting dynamically in concert, corresponding to spontaneous events, and which do not necessarily involve all voxels within the network nodes oscillating in unison.     

Altered small-world brain networks in schizophrenia patients during working memory performance

Impairment of working memory (WM) performance in schizophrenia patients (SZ) is well-established. Compared to healthy controls (HC), SZ patients show aberrant blood oxygen level dependent (BOLD) activations and disrupted functional connectivity during WM performance. In this study, we examined the small-world network metrics computed from functional magnetic resonance imaging (fMRI) data collected as 35 HC and 35 SZ performed a Sternberg Item Recognition Paradigm (SIRP) at three WM load levels. Functional connectivity networks were built by calculating the partial correlation on preprocessed time courses of BOLD signal between task-related brain regions of interest (ROIs) defined by group independent component analysis (ICA). The networks were then thresholded within the small-world regime, resulting in undirected binarized small-world networks at different working memory loads. Our results showed: 1) at the medium WM load level, the networks in SZ showed a lower clustering coefficient and less local efficiency compared with HC; 2) in SZ, most network measures altered significantly as the WM load level increased from low to medium and from medium to high, while the network metrics were relatively stable in HC at different WM loads; and 3) the altered structure at medium WM load in SZ was related to their performance during the task, with longer reaction time related to lower clustering coefficient and lower local efficiency. These findings suggest brain connectivity in patients with SZ was more diffuse and less strongly linked locally in functional network at intermediate level of WM when compared to HC. SZ show distinctly inefficient and variable network structures in response to WM load increase, comparing to stable highly clustered network topologies in HC.     

Large-Scale Functional Networks Identified from Resting-State EEG Using Spatial ICA

Several methods have been applied to EEG or MEG signals to detect functional networks. In recent works using MEG/EEG and fMRI data, temporal ICA analysis has been used to extract spatial maps of resting-state networks with or without an atlas-based parcellation of the cortex. Since the links between the fMRI signal and the electromagnetic signals are not fully established, and to avoid any bias, we examined whether EEG alone was able to derive the spatial distribution and temporal characteristics of functional networks. To do so, we propose a two-step original method: 1) An individual multi-frequency data analysis including EEG-based source localisation and spatial independent component analysis, which allowed us to characterize the resting-state networks. 2) A group-level analysis involving a hierarchical clustering procedure to identify reproducible large-scale networks across the population. Compared with large-scale resting-state networks obtained with fMRI, the proposed EEG-based analysis revealed smaller independent networks thanks to the high temporal resolution of EEG, hence hierarchical organization of networks. The comparison showed a substantial overlap between EEG and fMRI networks in motor, premotor, sensory, frontal, and parietal areas. However, there were mismatches between EEG-based and fMRI-based networks in temporal areas, presumably resulting from a poor sensitivity of fMRI in these regions or artefacts in the EEG signals. The proposed method opens the way for studying the high temporal dynamics of networks at the source level thanks to the high temporal resolution of EEG. It would then become possible to study detailed measures of the dynamics of connectivity.     

Emerging concepts for the dynamical organization of resting-state activity in the brain

A broad body of experimental work has demonstrated that apparently spontaneous brain activity is not random. At the level of large-scale neural systems, as measured with functional MRI (fMRI), this ongoing activity reflects the organization of a series of highly coherent functional networks. These so-called resting-state networks (RSNs) closely relate to the underlying anatomical connectivity but cannot be understood in those terms alone. Here we review three large-scale neural system models of primate neocortex that emphasize the key contributions of local dynamics, signal transmission delays and noise to the emerging RSNs. We propose that the formation and dissolution of resting-state patterns reflects the exploration of possible functional network configurations around a stable anatomical skeleton.     

Consistency of network modules in resting-state FMRI connectome data

At rest, spontaneous brain activity measured by fMRI is summarized by a number of distinct resting state networks (RSNs) following similar temporal time courses. Such networks have been consistently identified across subjects using spatial ICA (independent component analysis). Moreover, graph theory-based network analyses have also been applied to resting-state fMRI data, identifying similar RSNs, although typically at a coarser spatial resolution. In this work, we examined resting-state fMRI networks from 194 subjects at a voxel-level resolution, and examined the consistency of RSNs across subjects using a metric called scaled inclusivity (SI), which summarizes consistency of modular partitions across networks. Our SI analyses indicated that some RSNs are robust across subjects, comparable to the corresponding RSNs identified by ICA. We also found that some commonly reported RSNs are less consistent across subjects. This is the first direct comparison of RSNs between ICAs and graph-based network analyses at a comparable resolution.     

Disturbed Small-World Networks and Neurocognitive Function in Frontal Lobe Low-Grade Glioma Patients

Background

Brain tumor patients often associated with losses of the small-world configuration and neurocognitive functions before operations. However, few studies were performed on the impairments of frontal lobe low-grade gliomas (LGG) after tumor resection using small-world network features.

Methodology/Principal Findings

To detect differences in the whole brain topology among LGG patients before and after operation, a combined study of neurocognitive assessment and graph theoretical network analysis of fMRI data was performed. We collected resting-state fMRI data of 12 carefully selected frontal lobe LGG patients before and after operation. We calculated the topological properties of brain functional networks in the 12 LGG, and compared with 12 healthy controls (HCs). We also applied Montreal Cognitive Assessment (MoCA) in a subset of patients (n = 12, including before and after operation groups) and HCs (n = 12). The resulting functional connectivity matrices were constructed for all 12 patients, and binary network analysis was performed. In the range of , the functional networks in preoperative LGG and postoperative one both fitted the definition of small-worldness. We proposed as small-world network interval, and the results showed that the topological properties were found to be disrupted in the two LGG groups, meanwhile the global efficiency increased and the local efficiency decreased. in the two LGG groups both were longer than HCs. in the LGG groups were smaller than HCs. Compared with the Hcs, MoCA in the two LGG groups were lower than HCs with significant difference, and the disturbed networks in the LGG were negatively related to worse MoCA scores.

Conclusions

Disturbed small-worldness preperty in the two LGG groups was found and widely spread in the strength and spatial organization of brain networks, and the alterated small-world network may be responsible for cognitive dysfunction in frontal lobe LGG patients.     

A Study of Brain Networks Associated with Swallowing Using Graph-Theoretical Approaches

Functional connectivity between brain regions during swallowing tasks is still not well understood. Understanding these complex interactions is of great interest from both a scientific and a clinical perspective. In this study, functional magnetic resonance imaging (fMRI) was utilized to study brain functional networks during voluntary saliva swallowing in twenty-two adult healthy subjects (all females, years of age). To construct these functional connections, we computed mean partial correlation matrices over ninety brain regions for each participant. Two regions were determined to be functionally connected if their correlation was above a certain threshold. These correlation matrices were then analyzed using graph-theoretical approaches. In particular, we considered several network measures for the whole brain and for swallowing-related brain regions. The results have shown that significant pairwise functional connections were, mostly, either local and intra-hemispheric or symmetrically inter-hemispheric. Furthermore, we showed that all human brain functional network, although varying in some degree, had typical small-world properties as compared to regular networks and random networks. These properties allow information transfer within the network at a relatively high efficiency. Swallowing-related brain regions also had higher values for some of the network measures in comparison to when these measures were calculated for the whole brain. The current results warrant further investigation of graph-theoretical approaches as a potential tool for understanding the neural basis of dysphagia.     

Spatial vs. Temporal Features in ICA of Resting-State fMRI – A Quantitative and Qualitative Investigation in the Context of Response Inhibition

Independent component analysis (ICA) can identify covarying functional networks in the resting brain. Despite its relatively widespread use, the potential of the temporal information (unlike spatial information) obtained by ICA from resting state fMRI (RS-fMRI) data is not always fully utilized. In this study, we systematically investigated which features in ICA of resting-state fMRI relate to behaviour, with stop signal reaction time (SSRT) in a stop-signal task taken as a test case. We did this by correlating SSRT with the following three kinds of measure obtained from RS-fMRI data: (1) the amplitude of each resting state network (RSN) (evaluated by the standard deviation of the RSN timeseries), (2) the temporal correlation between every pair of RSN timeseries, and (3) the spatial map of each RSN. For multiple networks, we found significant correlations not only between SSRT and spatial maps, but also between SSRT and network activity amplitude. Most of these correlations are of functional interpretability. The temporal correlations between RSN pairs were of functional significance, but these correlations did not appear to be very sensitive to finding SSRT correlations. In addition, we also investigated the effects of the decomposition dimension, spatial smoothing and Z-transformation of the spatial maps, as well as the techniques for evaluating the temporal correlation between RSN timeseries. Overall, the temporal information acquired by ICA enabled us to investigate brain function from a complementary perspective to the information provided by spatial maps.     

Topological Organization of Functional Brain Networks in Healthy Children: Differences in Relation to Age,Sex, and Intelligence

Recent studies have demonstrated developmental changes of functional brain networks derived from functional connectivity using graph theoretical analysis, which has been rapidly translated to studies of brain network organization. However, little is known about sex- and IQ-related differences in the topological organization of functional brain networks during development. In this study, resting-state fMRI (rs-fMRI) was used to map the functional brain networks in 51 healthy children. We then investigated the effects of age, sex, and IQ on economic small-world properties and regional nodal properties of the functional brain networks. At a global level of whole networks, we found significant age-related increases in the small-worldness and local efficiency, significant higher values of the global efficiency in boys compared with girls, and no significant IQ-related difference. Age-related increases in the regional nodal properties were found predominately in the frontal brain regions, whereas the parietal, temporal, and occipital brain regions showed age-related decreases. Significant sex-related differences in the regional nodal properties were found in various brain regions, primarily related to the default mode, language, and vision systems. Positive correlations between IQ and the regional nodal properties were found in several brain regions related to the attention system, whereas negative correlations were found in various brain regions primarily involved in the default mode, emotion, and language systems. Together, our findings of the network topology of the functional brain networks in healthy children and its relationship with age, sex, and IQ bring new insights into the understanding of brain maturation and cognitive development during childhood and adolescence.