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.     

Loss of 'small-world' networks in Alzheimer's disease: graph analysis of FMRI resting-state functional connectivity

Background

Local network connectivity disruptions in Alzheimer''s disease patients have been found using graph analysis in BOLD fMRI. Other studies using MEG and cortical thickness measures, however, show more global long distance connectivity changes, both in functional and structural imaging data. The form and role of functional connectivity changes thus remains ambiguous. The current study shows more conclusive data on connectivity changes in early AD using graph analysis on resting-state condition fMRI data.

Methodology/Principal Findings

18 mild AD patients and 21 healthy age-matched control subjects without memory complaints were investigated in resting-state condition with MRI at 1.5 Tesla. Functional coupling between brain regions was calculated on the basis of pair-wise synchronizations between regional time-series. Local (cluster coefficient) and global (path length) network measures were quantitatively defined. Compared to controls, the characteristic path length of AD functional networks is closer to the theoretical values of random networks, while no significant differences were found in cluster coefficient. The whole-brain average synchronization does not differ between Alzheimer and healthy control groups. Post-hoc analysis of the regional synchronization reveals increased AD synchronization involving the frontal cortices and generalized decreases located at the parietal and occipital regions. This effectively translates in a global reduction of functional long-distance links between frontal and caudal brain regions.

Conclusions/Significance

We present evidence of AD-induced changes in global brain functional connectivity specifically affecting long-distance connectivity. This finding is highly relevant for it supports the anterior-posterior disconnection theory and its role in AD. Our results can be interpreted as reflecting the randomization of the brain functional networks in AD, further suggesting a loss of global information integration in disease.     

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.     

The graph theoretical analysis of the SSVEP harmonic response networks

Steady-state visually evoked potentials (SSVEP) have been widely used in the neural engineering and cognitive neuroscience researches. Previous studies have indicated that the SSVEP fundamental frequency responses are correlated with the topological properties of the functional networks entrained by the periodic stimuli. Given the different spatial and functional roles of the fundamental frequency and harmonic responses, in this study we further investigated the relation between the harmonic responses and the corresponding functional networks, using the graph theoretical analysis. We found that the second harmonic responses were positively correlated to the mean functional connectivity, clustering coefficient, and global and local efficiencies, while negatively correlated with the characteristic path lengths of the corresponding networks. In addition, similar pattern occurred with the lowest stimulus frequency (6.25 Hz) at the third harmonic responses. These findings demonstrate that more efficient brain networks are related to larger SSVEP responses. Furthermore, we showed that the main connection pattern of the SSVEP harmonic response networks originates from the interactions between the frontal and parietal–occipital regions. Overall, this study may bring new insights into the understanding of the brain mechanisms underlying SSVEP.     

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.     

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.     

Disrupted small-world brain networks in moderate Alzheimer's disease: a resting-state FMRI study

The small-world organization has been hypothesized to reflect a balance between local processing and global integration in the human brain. Previous multimodal imaging studies have consistently demonstrated that the topological architecture of the brain network is disrupted in Alzheimer's disease (AD). However, these studies have reported inconsistent results regarding the topological properties of brain alterations in AD. One potential explanation for these inconsistent results lies with the diverse homogeneity and distinct progressive stages of the AD involved in these studies, which are thought to be critical factors that might affect the results. We investigated the topological properties of brain functional networks derived from resting functional magnetic resonance imaging (fMRI) of carefully selected moderate AD patients and normal controls (NCs). Our results showed that the topological properties were found to be disrupted in AD patients, which showing increased local efficiency but decreased global efficiency. We found that the altered brain regions are mainly located in the default mode network, the temporal lobe and certain subcortical regions that are closely associated with the neuropathological changes in AD. Of note, our exploratory study revealed that the ApoE genotype modulates brain network properties, especially in AD patients.     

Graph analysis of functional brain networks: practical issues in translational neuroscience

The brain can be regarded as a network: a connected system where nodes, or units, represent different specialized regions and links, or connections, represent communication pathways. From a functional perspective, communication is coded by temporal dependence between the activities of different brain areas. In the last decade, the abstract representation of the brain as a graph has allowed to visualize functional brain networks and describe their non-trivial topological properties in a compact and objective way. Nowadays, the use of graph analysis in translational neuroscience has become essential to quantify brain dysfunctions in terms of aberrant reconfiguration of functional brain networks. Despite its evident impact, graph analysis of functional brain networks is not a simple toolbox that can be blindly applied to brain signals. On the one hand, it requires the know-how of all the methodological steps of the pipeline that manipulate the input brain signals and extract the functional network properties. On the other hand, knowledge of the neural phenomenon under study is required to perform physiologically relevant analysis. The aim of this review is to provide practical indications to make sense of brain network analysis and contrast counterproductive attitudes.     

Musical training induces functional plasticity in perceptual and motor networks: insights from resting-state FMRI

A number of previous studies have examined music-related plasticity in terms of multi-sensory and motor integration but little is known about the functional and effective connectivity patterns of spontaneous intrinsic activity in these systems during the resting state in musicians. Using functional connectivity and Granger causal analysis, functional and effective connectivity among the motor and multi-sensory (visual, auditory and somatosensory) cortices were evaluated using resting-state functional magnetic resonance imaging (fMRI) in musicians and non-musicians. The results revealed that functional connectivity was significantly increased in the motor and multi-sensory cortices of musicians. Moreover, the Granger causality results demonstrated a significant increase outflow-inflow degree in the auditory cortex with the strongest causal outflow pattern of effective connectivity being found in musicians. These resting state fMRI findings indicate enhanced functional integration among the lower-level perceptual and motor networks in musicians, and may reflect functional consolidation (plasticity) resulting from long-term musical training, involving both multi-sensory and motor functional integration.     

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.     

Dysregulated daily rhythmicity of neuronal resting-state networks in MCI patients

In young healthy participants, the degree of daily rhythmicity largely varies across different neuronal resting-state networks (RSNs), while it is to date unknown whether this temporal pattern of activity is conserved in healthy and pathological aging. Twelve healthy elderly (mean age?=?65.1?±?5.7 years) and 12 patients with amnestic mild cognitive impairment (aMCI; mean age?=?69.6?±?6.2 years) underwent four resting-state functional magnetic resonance imaging scans at fixed 2.5?h intervals throughout a day. Time courses of a RSN were extracted by a connectivity strength and a spatial extent approach performed individually for each participant. Highly rhythmic RSNs included a sensorimotor, a cerebellar and a visual network in healthy elderly; the least rhythmic RSNs in this group included a network associated with executive control and an orbitofrontal network. The degree of daily rhythmicity in aMCI patients was reduced and dysregulated. For healthy elderly, the findings are in accordance with results reported for young healthy participants suggesting a comparable distribution of daily rhythmicity across RSNs during healthy aging. In contrast, the reduction and dysregulation of daily rhythmicity observed in aMCI patients is presumably indicative of underlying neurodegenerative processes in this group.     

Identify schizophrenia using resting-state functional connectivity: an exploratory research and analysis

ABSTRACT: BACKGROUND: Schizophrenia is a severe mental illness associated with the symptoms such as hallucination and delusion. The objective of this study was to investigate the abnormal resting-state functional connectivity patterns of schizophrenic patients which could identify furthest patients from healthy controls. METHODS: The whole-brain resting-state fMRI was performed on patients diagnosed with schizophrenia (n=22) and on age- and gender-matched, healthy control subjects (n=22). To differentiate schizophrenic individuals from healthy controls, the multivariate classification analysis was employed. The weighted brain regions were got by reconstruction arithmetic to extract highly discriminative functional connectivity information. RESULTS: The results showed that 93.2% (p<0.001) of the subjects were correctly classified via the leave-one-out cross-validation method. And most of the altered functional connections identified located within the visual cortical-, default-mode-, and sensorimotor network. Furthermore, in reconstruction arithmetic, the fusiform gyrus exhibited the greatest amount of weight. CONCLUSIONS: This study demonstrates that schizophrenic patients may be successfully differentiated from healthy subjects by using whole-brain resting-state fMRI, and the fusiform gyrus may play an important functional role in the physiological symptoms manifested by schizophrenic patients. The brain region of great weight may be the problematic region of information exchange in schizophrenia. Thus, our result may provide insights into the identification of potentially effective biomarkers for the clinical diagnosis of schizophrenia.     

Auditory resting-state network connectivity in tinnitus: a functional MRI study

The underlying functional neuroanatomy of tinnitus remains poorly understood. Few studies have focused on functional cerebral connectivity changes in tinnitus patients. The aim of this study was to test if functional MRI "resting-state" connectivity patterns in auditory network differ between tinnitus patients and normal controls. Thirteen chronic tinnitus subjects and fifteen age-matched healthy controls were studied on a 3 tesla MRI. Connectivity was investigated using independent component analysis and an automated component selection approach taking into account the spatial and temporal properties of each component. Connectivity in extra-auditory regions such as brainstem, basal ganglia/NAc, cerebellum, parahippocampal, right prefrontal, parietal, and sensorimotor areas was found to be increased in tinnitus subjects. The right primary auditory cortex, left prefrontal, left fusiform gyrus, and bilateral occipital regions showed a decreased connectivity in tinnitus. These results show that there is a modification of cortical and subcortical functional connectivity in tinnitus encompassing attentional, mnemonic, and emotional networks. Our data corroborate the hypothesized implication of non-auditory regions in tinnitus physiopathology and suggest that various regions of the brain seem involved in the persistent awareness of the phenomenon as well as in the development of the associated distress leading to disabling chronic tinnitus.     

Altered small-world,functional brain networks in patients with lower back pain

In this study, we aimed to investigate the functional network changes that occur in patients with lower back pain(LBP). We also investigated the link between LBP and the small-world properties of functional networks within the brain. Functional MRI(fMRI) was performed on 20 individuals with LBP and 17 age and gender-matched normal controls during the resting state. The severity of the pain in the individuals with LBP ranged from 5 to 8 on a 0–10 scale, with 0 indicating no pain. Network-based statistics were performed to investigate the differences between the brain networks of individuals with LBP and those of normal controls. Several small-world parameters of brain networks were calculated, including the clustering coefficient, characteristic path length, local efficiency, and global efficiency. These criteria reflect the overall network efficiency. The brain networks in the individuals with LBP due to herniation of a lumbar disc demonstrated a significantly longer characteristic path length as well as a lower clustering coefficient, global efficiency, and local efficiency compared to those in control subjects. We found that LBP patients tended to have unstable and inefficient brain networks when compared with healthy controls. In addition, LBP individuals showed significantly decreased functional connectivity in the anterior cingulate cortex, middle cingulate cortex, post cingulate cortex, inferior frontal gyrus, middle temporal gyrus, occipital gyrus, postcentral gyrus, precentral gyrus, supplementary motor area, thalamus, fusiform, caudate, and cerebellum. We believe that these regions may be involved in the pathophysiology of lower back pain.     

Efficiency and cost of economical brain functional networks

Brain anatomical networks are sparse, complex, and have economical small-world properties. We investigated the efficiency and cost of human brain functional networks measured using functional magnetic resonance imaging (fMRI) in a factorial design: two groups of healthy old (N = 11; mean age = 66.5 years) and healthy young (N = 15; mean age = 24.7 years) volunteers were each scanned twice in a no-task or “resting” state following placebo or a single dose of a dopamine receptor antagonist (sulpiride 400 mg). Functional connectivity between 90 cortical and subcortical regions was estimated by wavelet correlation analysis, in the frequency interval 0.06–0.11 Hz, and thresholded to construct undirected graphs. These brain functional networks were small-world and economical in the sense of providing high global and local efficiency of parallel information processing for low connection cost. Efficiency was reduced disproportionately to cost in older people, and the detrimental effects of age on efficiency were localised to frontal and temporal cortical and subcortical regions. Dopamine antagonism also impaired global and local efficiency of the network, but this effect was differentially localised and did not interact with the effect of age. Brain functional networks have economical small-world properties—supporting efficient parallel information transfer at relatively low cost—which are differently impaired by normal aging and pharmacological blockade of dopamine transmission.     

Downregulated expression of HSP27 in human low-grade glioma tissues discovered by a quantitative proteomic analysis

Background  

Heat shock proteins (HSPs), including mainly HSP110, HSP90, HSP70, HSP60 and small HSP families, are evolutionary conserved proteins involved in various cellular processes. Abnormal expression of HSPs has been detected in several tumor types, which indicates that specific HSPs have different prognostic significance for different tumors. In the current studies, the expression profiling of HSPs in human low-grade glioma tissues (HGTs) were investigated using a sensitive, accurate SILAC (stable isotope labeling with amino acids in cell culture)-based quantitative proteomic strategy.     

Comparing brain networks of different size and connectivity density using graph theory

Graph theory is a valuable framework to study the organization of functional and anatomical connections in the brain. Its use for comparing network topologies, however, is not without difficulties. Graph measures may be influenced by the number of nodes (N) and the average degree (k) of the network. The explicit form of that influence depends on the type of network topology, which is usually unknown for experimental data. Direct comparisons of graph measures between empirical networks with different N and/or k can therefore yield spurious results. We list benefits and pitfalls of various approaches that intend to overcome these difficulties. We discuss the initial graph definition of unweighted graphs via fixed thresholds, average degrees or edge densities, and the use of weighted graphs. For instance, choosing a threshold to fix N and k does eliminate size and density effects but may lead to modifications of the network by enforcing (ignoring) non-significant (significant) connections. Opposed to fixing N and k, graph measures are often normalized via random surrogates but, in fact, this may even increase the sensitivity to differences in N and k for the commonly used clustering coefficient and small-world index. To avoid such a bias we tried to estimate the N,k-dependence for empirical networks, which can serve to correct for size effects, if successful. We also add a number of methods used in social sciences that build on statistics of local network structures including exponential random graph models and motif counting. We show that none of the here-investigated methods allows for a reliable and fully unbiased comparison, but some perform better than others.     

Clusters in alpha/beta barrel proteins: implications for protein structure, function, and folding: a graph theoretical approach

The alpha/beta barrel fold is adopted by most enzymes performing a variety of catalytic reactions, but with very low sequence similarity. In order to understand the stabilizing interactions important in maintaining the alpha/beta barrel fold, we have identified residue clusters in a dataset of 36 alpha/beta barrel proteins that have less than 10% sequence identity within themselves. A graph theoretical algorithm is used to identify backbone clusters. This approach uses the global information of the nonbonded interaction in the alpha/beta barrel fold for the clustering procedure. The nonbonded interactions are represented mathematically in the form of an adjacency matrix. On diagonalizing the adjacency matrix, clusters and cluster centers are obtained from the highest eigenvalue and its corresponding vector components. Residue clusters are identified in the strand regions forming the beta barrel and are topologically conserved in all 36 proteins studied. The residues forming the cluster in each of the alpha/beta protein are also conserved among the sequences belonging to the same family. The cluster centers are found to occur in the middle of the strands or in the C-terminal of the strands. In most cases, the residues forming the clusters are part of the active site or are located close to the active site. The folding nucleus of the alpha/beta fold is predicted based on hydrophobicity index evaluation of residues and identification of cluster centers. The predicted nucleation sites are found to occur mostly in the middle of the strands. Proteins 2001;43:103-112.