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.     

Abnormal reorganization of functional cortical small-world networks in focal hand dystonia

We investigated the large-scale functional cortical connectivity network in focal hand dystonia (FHD) patients using graph theoretic measures to assess efficiency. High-resolution EEGs were recorded in 15 FHD patients and 15 healthy volunteers at rest and during a simple sequential finger tapping task. Mutual information (MI) values of wavelet coefficients were estimated to create an association matrix between EEG electrodes, and to produce a series of adjacency matrices or graphs, G, by thresholding with network cost. Efficiency measures of small-world networks were assessed. As a result, we found that FHD patients have economical small-world properties in their brain functional networks in the alpha and beta bands. During a motor task, in the beta band network, FHD patients have decreased efficiency of small-world networks, whereas healthy volunteers increase efficiency. Reduced efficient beta band network in FHD patients during the task was consistently observed in global efficiency, cost-efficiency, and maximum cost-efficiency. This suggests that the beta band functional cortical network of FHD patients is reorganized even during a task that does not induce dystonic symptoms, representing a loss of long-range communication and abnormal functional integration in large-scale brain functional cortical networks. Moreover, negative correlations between efficiency measures and duration of disease were found, indicating that the longer duration of disease, the less efficient the beta band network in FHD patients. In regional efficiency analysis, FHD patients at rest have high regional efficiency at supplementary motor cortex (SMA) compared with healthy volunteers; however, it is diminished during the motor task, possibly reflecting abnormal inhibition in FHD patients. The present study provides the first evidence with graph theory for abnormal reconfiguration of brain functional networks in FHD during motor task.     

EEG-Based Functional Brain Networks: Does the Network Size Matter?

Functional connectivity in human brain can be represented as a network using electroencephalography (EEG) signals. These networks--whose nodes can vary from tens to hundreds--are characterized by neurobiologically meaningful graph theory metrics. This study investigates the degree to which various graph metrics depend upon the network size. To this end, EEGs from 32 normal subjects were recorded and functional networks of three different sizes were extracted. A state-space based method was used to calculate cross-correlation matrices between different brain regions. These correlation matrices were used to construct binary adjacency connectomes, which were assessed with regards to a number of graph metrics such as clustering coefficient, modularity, efficiency, economic efficiency, and assortativity. We showed that the estimates of these metrics significantly differ depending on the network size. Larger networks had higher efficiency, higher assortativity and lower modularity compared to those with smaller size and the same density. These findings indicate that the network size should be considered in any comparison of networks across studies.     

Graph theoretical analysis of complex networks in the brain

Since the discovery of small-world and scale-free networks the study of complex systems from a network perspective has taken an enormous flight. In recent years many important properties of complex networks have been delineated. In particular, significant progress has been made in understanding the relationship between the structural properties of networks and the nature of dynamics taking place on these networks. For instance, the 'synchronizability' of complex networks of coupled oscillators can be determined by graph spectral analysis. These developments in the theory of complex networks have inspired new applications in the field of neuroscience. Graph analysis has been used in the study of models of neural networks, anatomical connectivity, and functional connectivity based upon fMRI, EEG and MEG. These studies suggest that the human brain can be modelled as a complex network, and may have a small-world structure both at the level of anatomical as well as functional connectivity. This small-world structure is hypothesized to reflect an optimal situation associated with rapid synchronization and information transfer, minimal wiring costs, as well as a balance between local processing and global integration. The topological structure of functional networks is probably restrained by genetic and anatomical factors, but can be modified during tasks. There is also increasing evidence that various types of brain disease such as Alzheimer's disease, schizophrenia, brain tumours and epilepsy may be associated with deviations of the functional network topology from the optimal small-world pattern.     

EEGNET: An Open Source Tool for Analyzing and Visualizing M/EEG Connectome

The brain is a large-scale complex network often referred to as the “connectome”. Exploring the dynamic behavior of the connectome is a challenging issue as both excellent time and space resolution is required. In this context Magneto/Electroencephalography (M/EEG) are effective neuroimaging techniques allowing for analysis of the dynamics of functional brain networks at scalp level and/or at reconstructed sources. However, a tool that can cover all the processing steps of identifying brain networks from M/EEG data is still missing. In this paper, we report a novel software package, called EEGNET, running under MATLAB (Math works, inc), and allowing for analysis and visualization of functional brain networks from M/EEG recordings. EEGNET is developed to analyze networks either at the level of scalp electrodes or at the level of reconstructed cortical sources. It includes i) Basic steps in preprocessing M/EEG signals, ii) the solution of the inverse problem to localize / reconstruct the cortical sources, iii) the computation of functional connectivity among signals collected at surface electrodes or/and time courses of reconstructed sources and iv) the computation of the network measures based on graph theory analysis. EEGNET is the unique tool that combines the M/EEG functional connectivity analysis and the computation of network measures derived from the graph theory. The first version of EEGNET is easy to use, flexible and user friendly. EEGNET is an open source tool and can be freely downloaded from this webpage: https://sites.google.com/site/eegnetworks/.     

Graph analysis of functional brain network topology using minimum spanning tree in driver drowsiness

A large number of traffic accidents due to driver drowsiness have been under more attention of many countries. The organization of the functional brain network is associated with drowsiness, but little is known about the brain network topology that is modulated by drowsiness. To clarify this problem, in this study, we introduce a novel approach to detect driver drowsiness. Electroencephalogram (EEG) signals have been measured during a simulated driving task, in which participants are recruited to undergo both alert and drowsy states. The filtered EEG signals are then decomposed into multiple frequency bands by wavelet packet transform. Functional connectivity between all pairs of channels for multiple frequency bands is assessed using the phase lag index (PLI). Based on this, PLI-weighted networks are subsequently calculated, from which minimum spanning trees are constructed—a graph method that corrects for comparison bias. Statistical analyses are performed on graph-derived metrics as well as on the PLI connectivity values. The major finding is that significant differences in the delta frequency band for three graph metrics and in the theta frequency band for five graph metrics suggesting network integration and communication between network nodes are increased from alertness to drowsiness. Together, our findings also suggest a more line-like configuration in alert states and a more star-like topology in drowsy states. Collectively, our findings point to a more proficient configuration in drowsy state for lower frequency bands. Graph metrics relate to the intrinsic organization of functional brain networks, and these graph metrics may provide additional insights on driver drowsiness detection for reducing and preventing traffic accidents and further understanding the neural mechanisms of driver drowsiness.     

Whole-brain Functional Networks in Cognitively Normal,Mild Cognitive Impairment,and Alzheimer’s Disease

The conceptual significance of understanding functional brain alterations and cognitive deficits associated with Alzheimer’s disease (AD) process has been widely established. However, the whole-brain functional networks of AD and its prodromal stage, mild cognitive impairment (MCI), are not well clarified yet. In this study, we compared the characteristics of the whole-brain functional networks among cognitively normal (CN), MCI, and AD individuals by applying graph theoretical analyses to [¹⁸F] fluorodeoxyglucose positron emission tomography (FDG-PET) data. Ninety-four CN elderly, 183 with MCI, and 216 with AD underwent clinical evaluation and FDG-PET scan. The overall small-world property as seen in the CN whole-brain network was preserved in MCI and AD. In contrast, individual parameters of the network were altered with the following patterns of changes: local clustering of networks was lower in both MCI and AD compared to CN, while path length was not different among the three groups. Then, MCI had a lower level of local clustering than AD. Subgroup analyses for AD also revealed that very mild AD had lower local clustering and shorter path length compared to mild AD. Regarding the local properties of the whole-brain networks, MCI and AD had significantly decreased normalized betweenness centrality in several hubs regionally associated with the default mode network compared to CN. Our results suggest that the functional integration in whole-brain network progressively declines due to the AD process. On the other hand, functional relatedness between neighboring brain regions may not gradually decrease, but be the most severely altered in MCI stage and gradually re-increase in clinical AD stages.     

GAT: a graph-theoretical analysis toolbox for analyzing between-group differences in large-scale structural and functional brain networks

In recent years, graph theoretical analyses of neuroimaging data have increased our understanding of the organization of large-scale structural and functional brain networks. However, tools for pipeline application of graph theory for analyzing topology of brain networks is still lacking. In this report, we describe the development of a graph-analysis toolbox (GAT) that facilitates analysis and comparison of structural and functional network brain networks. GAT provides a graphical user interface (GUI) that facilitates construction and analysis of brain networks, comparison of regional and global topological properties between networks, analysis of network hub and modules, and analysis of resilience of the networks to random failure and targeted attacks. Area under a curve (AUC) and functional data analyses (FDA), in conjunction with permutation testing, is employed for testing the differences in network topologies; analyses that are less sensitive to the thresholding process. We demonstrated the capabilities of GAT by investigating the differences in the organization of regional gray-matter correlation networks in survivors of acute lymphoblastic leukemia (ALL) and healthy matched Controls (CON). The results revealed an alteration in small-world characteristics of the brain networks in the ALL survivors; an observation that confirm our hypothesis suggesting widespread neurobiological injury in ALL survivors. Along with demonstration of the capabilities of the GAT, this is the first report of altered large-scale structural brain networks in ALL survivors.     

Extracting Labeled Topological Patterns from Samples of Networks

An advanced graph theoretical approach is introduced that enables a higher level of functional interpretation of samples of directed networks with identical fixed pairwise different vertex labels that are drawn from a particular population. Compared to the analysis of single networks, their investigation promises to yield more detailed information about the represented system. Often patterns of directed edges in sample element networks are too intractable for a direct evaluation and interpretation. The new approach addresses the problem of simplifying topological information and characterizes such a sample of networks by finding its locatable characteristic topological patterns. These patterns, essentially sample-specific network motifs with vertex labeling, might represent the essence of the intricate topological information contained in all sample element networks and provides as well a means of differentiating network samples. Central to the accurateness of this approach is the null model and its properties, which is needed to assign significance to topological patterns. As a proof of principle the proposed approach has been applied to the analysis of networks that represent brain connectivity before and during painful stimulation in patients with major depression and in healthy subjects. The accomplished reduction of topological information enables a cautious functional interpretation of the altered neuronal processing of pain in both groups.     

Integrative network biology: graph prototyping for co-expression cancer networks

Network-based analysis has been proven useful in biologically-oriented areas, e.g., to explore the dynamics and complexity of biological networks. Investigating a set of networks allows deriving general knowledge about the underlying topological and functional properties. The integrative analysis of networks typically combines networks from different studies that investigate the same or similar research questions. In order to perform an integrative analysis it is often necessary to compare the properties of matching edges across the data set. This identification of common edges is often burdensome and computational intensive. Here, we present an approach that is different from inferring a new network based on common features. Instead, we select one network as a graph prototype, which then represents a set of comparable network objects, as it has the least average distance to all other networks in the same set. We demonstrate the usefulness of the graph prototyping approach on a set of prostate cancer networks and a set of corresponding benign networks. We further show that the distances within the cancer group and the benign group are statistically different depending on the utilized distance measure.     

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.     

Brain functional network modeling and analysis based on fMRI: a systematic review

In recent years, the number of patients with neurodegenerative diseases (i.e., Alzheimer’s disease, Parkinson’s disease, mild cognitive impairment) and mental disorders (i.e., depression, anxiety and schizophrenia) have increased dramatically. Researchers have found that complex network analysis can reveal the topology of brain functional networks, such as small-world, scale-free, etc. In the study of brain diseases, it has been found that these topologies have undergoed abnormal changes in different degrees. Therefore, the research of brain functional networks can not only provide a new perspective for understanding the pathological mechanism of neurological and psychiatric diseases, but also provide assistance for the early diagnosis. Focusing on the study of human brain functional networks, this paper reviews the research results in recent years. First, this paper introduces the background of the study of brain functional networks under complex network theory and the important role of topological properties in the study of brain diseases. Second, the paper describes how to construct a brain functional network using neural image data. Third, the common methods of functional network analysis, including network structure analysis and disease classification, are introduced. Fourth, the role of brain functional networks in pathological study, analysis and diagnosis of brain functional diseases is studied. Finally, the paper summarizes the existing studies of brain functional networks and points out the problems and future research directions.     

Analysing Local Sparseness in the Macaque Brain Network

Understanding the network structure of long distance pathways in the brain is a necessary step towards developing an insight into the brain’s function, organization and evolution. Dense global subnetworks of these pathways have often been studied, primarily due to their functional implications. Instead we study sparse local subnetworks of the pathways to establish the role of a brain area in enabling shortest path communication between its non-adjacent topological neighbours. We propose a novel metric to measure the topological communication load on a vertex due to its immediate neighbourhood, and show that in terms of distribution of this local communication load, a network of Macaque long distance pathways is substantially different from other real world networks and random graph models. Macaque network contains the entire range of local subnetworks, from star-like networks to clique-like networks, while other networks tend to contain a relatively small range of subnetworks. Further, sparse local subnetworks in the Macaque network are not only found across topographical super-areas, e.g., lobes, but also within a super-area, arguing that there is conservation of even relatively short-distance pathways. To establish the communication role of a vertex we borrow the concept of brokerage from social science, and present the different types of brokerage roles that brain areas play, highlighting that not only the thalamus, but also cingulate gyrus and insula often act as “relays” for areas in the neocortex. These and other analysis of communication load and roles of the sparse subnetworks of the Macaque brain provide new insights into the organisation of its pathways.     

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.     

Towards the "baby connectome": mapping the structural connectivity of the newborn brain

Defining the structural and functional connectivity of the human brain (the human "connectome") is a basic challenge in neuroscience. Recently, techniques for noninvasively characterizing structural connectivity networks in the adult brain have been developed using diffusion and high-resolution anatomic MRI. The purpose of this study was to establish a framework for assessing structural connectivity in the newborn brain at any stage of development and to show how network properties can be derived in a clinical cohort of six-month old infants sustaining perinatal hypoxic ischemic encephalopathy (HIE). Two different anatomically unconstrained parcellation schemes were proposed and the resulting network metrics were correlated with neurological outcome at 6 months. Elimination and correction of unreliable data, automated parcellation of the cortical surface, and assembling the large-scale baby connectome allowed an unbiased study of the network properties of the newborn brain using graph theoretic analysis. In the application to infants with HIE, a trend to declining brain network integration and segregation was observed with increasing neuromotor deficit scores.     

A linear model for characterization of synchronization frequencies of neural networks

The synchronization frequency of neural networks and its dynamics have important roles in deciphering the working mechanisms of the brain. It has been widely recognized that the properties of functional network synchronization and its dynamics are jointly determined by network topology, network connection strength, i.e., the connection strength of different edges in the network, and external input signals, among other factors. However, mathematical and computational characterization of the relationships between network synchronization frequency and these three important factors are still lacking. This paper presents a novel computational simulation framework to quantitatively characterize the relationships between neural network synchronization frequency and network attributes and input signals. Specifically, we constructed a series of neural networks including simulated small-world networks, real functional working memory network derived from functional magnetic resonance imaging, and real large-scale structural brain networks derived from diffusion tensor imaging, and performed synchronization simulations on these networks via the Izhikevich neuron spiking model. Our experiments demonstrate that both of the network synchronization strength and synchronization frequency change according to the combination of input signal frequency and network self-synchronization frequency. In particular, our extensive experiments show that the network synchronization frequency can be represented via a linear combination of the network self-synchronization frequency and the input signal frequency. This finding could be attributed to an intrinsically-preserved principle in different types of neural systems, offering novel insights into the working mechanism of neural systems.     

Effects of different correlation metrics and preprocessing factors on small-world brain functional networks: a resting-state functional MRI study

Graph theoretical analysis of brain networks based on resting-state functional MRI (R-fMRI) has attracted a great deal of attention in recent years. These analyses often involve the selection of correlation metrics and specific preprocessing steps. However, the influence of these factors on the topological properties of functional brain networks has not been systematically examined. Here, we investigated the influences of correlation metric choice (Pearson's correlation versus partial correlation), global signal presence (regressed or not) and frequency band selection [slow-5 (0.01-0.027 Hz) versus slow-4 (0.027-0.073 Hz)] on the topological properties of both binary and weighted brain networks derived from them, and we employed test-retest (TRT) analyses for further guidance on how to choose the "best" network modeling strategy from the reliability perspective. Our results show significant differences in global network metrics associated with both correlation metrics and global signals. Analysis of nodal degree revealed differing hub distributions for brain networks derived from Pearson's correlation versus partial correlation. TRT analysis revealed that the reliability of both global and local topological properties are modulated by correlation metrics and the global signal, with the highest reliability observed for Pearson's-correlation-based brain networks without global signal removal (WOGR-PEAR). The nodal reliability exhibited a spatially heterogeneous distribution wherein regions in association and limbic/paralimbic cortices showed moderate TRT reliability in Pearson's-correlation-based brain networks. Moreover, we found that there were significant frequency-related differences in topological properties of WOGR-PEAR networks, and brain networks derived in the 0.027-0.073 Hz band exhibited greater reliability than those in the 0.01-0.027 Hz band. Taken together, our results provide direct evidence regarding the influences of correlation metrics and specific preprocessing choices on both the global and nodal topological properties of functional brain networks. This study also has important implications for how to choose reliable analytical schemes in brain network studies.     

Complex networks and deep learning for EEG signal analysis

Electroencephalogram (EEG) signals acquired from brain can provide an effective representation of the human’s physiological and pathological states. Up to now, much work has been conducted to study and analyze the EEG signals, aiming at spying the current states or the evolution characteristics of the complex brain system. Considering the complex interactions between different structural and functional brain regions, brain network has received a lot of attention and has made great progress in brain mechanism research. In addition, characterized by autonomous, multi-layer and diversified feature extraction, deep learning has provided an effective and feasible solution for solving complex classification problems in many fields, including brain state research. Both of them show strong ability in EEG signal analysis, but the combination of these two theories to solve the difficult classification problems based on EEG signals is still in its infancy. We here review the application of these two theories in EEG signal research, mainly involving brain–computer interface, neurological disorders and cognitive analysis. Furthermore, we also develop a framework combining recurrence plots and convolutional neural network to achieve fatigue driving recognition. The results demonstrate that complex networks and deep learning can effectively implement functional complementarity for better feature extraction and classification, especially in EEG signal analysis.     

Interplay between Graph Topology and Correlations of Third Order in Spiking Neuronal Networks

The study of processes evolving on networks has recently become a very popular research field, not only because of the rich mathematical theory that underpins it, but also because of its many possible applications, a number of them in the field of biology. Indeed, molecular signaling pathways, gene regulation, predator-prey interactions and the communication between neurons in the brain can be seen as examples of networks with complex dynamics. The properties of such dynamics depend largely on the topology of the underlying network graph. In this work, we want to answer the following question: Knowing network connectivity, what can be said about the level of third-order correlations that will characterize the network dynamics? We consider a linear point process as a model for pulse-coded, or spiking activity in a neuronal network. Using recent results from theory of such processes, we study third-order correlations between spike trains in such a system and explain which features of the network graph (i.e. which topological motifs) are responsible for their emergence. Comparing two different models of network topology—random networks of Erdős-Rényi type and networks with highly interconnected hubs—we find that, in random networks, the average measure of third-order correlations does not depend on the local connectivity properties, but rather on global parameters, such as the connection probability. This, however, ceases to be the case in networks with a geometric out-degree distribution, where topological specificities have a strong impact on average correlations.