Functional Connectivity Graphs Research Articles

Complexity of weighted graph: A new technique to investigate structural complexity of brain activities with applications to aging and autism

In recent years complexity of the brain structure in healthy and disordered subjects has been studied increasingly. But to the best of the authors’ knowledge, researchers so far have investigated the structural complexity only in the context of two restricted networks known as Small-World and Scale-free networks; whereas other aspects of the structural complexity of brain activities may be affected by aging and neurodegenerative disorders such as the Alzheimer’s disease and autism spectrum disorder. In this study, two general complexity metrics of graphs, Graph Index Complexity and Offdiagonal Complexity are proposed as general measures of complexity, not restricted to SWN only. They are adopted to measure the structural complexity of the weighted graphs instead of the common binary graphs. Fuzzy Synchronization Likelihood is applied to the EEGs and their sub-bands, as a functional connectivity metric of the brain, to construct the functional connectivity graphs. Two applications are used to evaluate the efficacy of the complexity measures: diagnosis of autism and aging, both based on EEG. It was discovered that the Graph Index Complexity of gamma band is discriminative in distinguishing autistic children from non-autistic children. Also, Offdiagonal Complexity of theta band in young subjects was observed to be significantly different than old subjects. This study shows that changes in the structure of functional connectivity of brain in disorders and different healthy states can be revealed by unrestricted metrics of graph complexity. While the applications presented in this paper are based on EEG, the approach is general and can be used with other modalities such as fMRI, MEG, etc. Further, it can be used to study every other neurological and psychiatric disorder.

Large-scale DCMs for resting-state fMRI.

This paper considers the identification of large directed graphs for resting-state brain networks based on biophysical models of distributed neuronal activity, that is, effective connectivity. This identification can be contrasted with functional connectivity methods based on symmetric correlations that are ubiquitous in resting-state functional MRI (fMRI). We use spectral dynamic causal modeling (DCM) to invert large graphs comprising dozens of nodes or regions. The ensuing graphs are directed and weighted, hence providing a neurobiologically plausible characterization of connectivity in terms of excitatory and inhibitory coupling. Furthermore, we show that the use of Bayesian model reduction to discover the most likely sparse graph (or model) from a parent (e.g., fully connected) graph eschews the arbitrary thresholding often applied to large symmetric (functional connectivity) graphs. Using empirical fMRI data, we show that spectral DCM furnishes connectivity estimates on large graphs that correlate strongly with the estimates provided by stochastic DCM. Furthermore, we increase the efficiency of model inversion using functional connectivity modes to place prior constraints on effective connectivity. In other words, we use a small number of modes to finesse the potentially redundant parameterization of large DCMs. We show that spectral DCM—with functional connectivity priors—is ideally suited for directed graph theoretic analyses of resting-state fMRI. We envision that directed graphs will prove useful in understanding the psychopathology and pathophysiology of neurodegenerative and neurodevelopmental disorders. We will demonstrate the utility of large directed graphs in clinical populations in subsequent reports, using the procedures described in this paper.

Mining cross-frequency coupling microstates from resting state MEG: An application to mild traumatic brain injury.

Recent studies have investigated the possible role of dynamic functional connectivity and the role of cross-frequency coupling (CFC) to provide the substrate for reliable biomarkers of brain disorders. In this study, we analyzed time-varying CFC profiles from resting state Magnetoencephal-ographic recordings of 30 mild Traumatic Brain Injury (mTBI) patients and 50 normal controls. Interactions among sensors at specific pairs of frequency bands were computed via estimation of phase-to-amplitude couplings. We then computed time-varying functional connectivity graphs that were described in terms of segregation (local efficiency, LE) and integration (global efficiency, GE) and mapped those graphs to time series of GE/LE estimates. The resulting dynamic network revealed transitions between a limited number of microstates for mTBI subjects compared to controls. The significant differences in transition probability between the two groups, along with the limited repertoire of possible states, can form the basis for a robust dynamic connectomic biomarker for the diagnosis of mTBI.

Impaired functional default mode network in patients with mild neurological Wilson’s disease

Wilson’s disease (WD) is an autosomal recessive metabolic disorder characterized by cognitive, psychiatric and motor signs and symptoms that are associated with structural and pathological brain abnormalities, in addition to liver changes. However, functional brain connectivity pattern of WD patients remains largely unknown. In the present study, we investigated functional brain connectivity pattern of WD patients using resting state functional magnetic resonance imaging. Particularly, we studied default mode network (DMN) using posterior cingulate cortex (PCC) based seed functional connectivity analysis and graph theoretic functional brain network analysis tools, and investigated the relationship between the DMN’s functional connectivity pattern of WD patients and their attention functions examined using the attention network test (ANT). Our results demonstrated that WD patients had altered DMN’s functional connectivity and lower local and global network efficiency compared with normal controls (NCs). In addition, the functional connectivity between left inferior temporal cortex and right lateral parietal cortex was correlated with altering function, one of the attention functions, across WD and NC subjects. These findings indicated that the DMN’s functional connectivity was altered in WD patients, which might be correlated with their attention dysfunction.

Functional Connectivity of the Dorsal Striatum in Female Musicians.

The dorsal striatum (caudate/putamen) is a node of the cortico-striato-pallido-thalamo-cortical (CSPTC) motor circuit, which plays a central role in skilled motor learning, a critical feature of musical performance. The dorsal striatum receives input from a large part of the cerebral cortex, forming a hub in the cortical-subcortical network. This study sought to examine how the functional network of the dorsal striatum differs between musicians and nonmusicians. Resting state functional magnetic resonance imaging (fMRI) data were acquired from female university students majoring in music and nonmusic disciplines. The data were subjected to functional connectivity analysis and graph theoretical analysis. The functional connectivity analysis indicated that compared with nonmusicians, musicians had significantly decreased connectivity between the left putamen and bilateral frontal operculum (FO) and between the left caudate nucleus and cerebellum. The graph theoretical analysis of the entire brain revealed that the degrees, which represent the numbers of connections, of the bilateral putamen were significantly lower in musicians than in nonmusicians. In conclusion, compared with nonmusicians, female musicians have a smaller functional network of the dorsal striatum with decreased connectivity. These data are consistent with previous anatomical studies reporting a reduced volume of the dorsal striatum in musicians and ballet dancers, suggesting that long-term musical training reshapes the functional network of the dorsal striatum to be less extensive or selective.

ToolConnect: A Functional Connectivity Toolbox for In vitro Networks.

Nowadays, the use of in vitro reduced models of neuronal networks to investigate the interplay between structural-functional connectivity and the emerging collective dynamics is a widely accepted approach. In this respect, a relevant advance for this kind of studies has been given by the recent introduction of high-density large-scale Micro-Electrode Arrays (MEAs) which have favored the mapping of functional connections and the recordings of the neuronal electrical activity. Although, several toolboxes have been implemented to characterize network dynamics and derive functional links, no specifically dedicated software for the management of huge amount of data and direct estimation of functional connectivity maps has been developed. toolconnect offers the implementation of up to date algorithms and a user-friendly Graphical User Interface (GUI) to analyze recorded data from large scale networks. It has been specifically conceived as a computationally efficient open-source software tailored to infer functional connectivity by analyzing the spike trains acquired from in vitro networks coupled to MEAs. In the current version, toolconnect implements correlation- (cross-correlation, partial-correlation) and information theory (joint entropy, transfer entropy) based core algorithms, as well as useful and practical add-ons to visualize functional connectivity graphs and extract some topological features. In this work, we present the software, its main features and capabilities together with some demonstrative applications on hippocampal recordings.

Altered cross-frequency coupling in resting-state MEG after mild traumatic brain injury

Cross-frequency coupling (CFC) is thought to represent a basic mechanism of functional integration of neural networks across distant brain regions. In this study, we analyzed CFC profiles from resting state Magnetoencephalographic (MEG) recordings obtained from 30 mild traumatic brain injury (mTBI) patients and 50 controls. We used mutual information (MI) to quantify the phase-to-amplitude coupling (PAC) of activity among the recording sensors in six nonoverlapping frequency bands. After forming the CFC-based functional connectivity graphs, we employed a tensor representation and tensor subspace analysis to identify the optimal set of features for subject classification as mTBI or control. Our results showed that controls formed a dense network of stronger local and global connections indicating higher functional integration compared to mTBI patients. Furthermore, mTBI patients could be separated from controls with more than 90% classification accuracy. These findings indicate that analysis of brain networks computed from resting-state MEG with PAC and tensorial representation of connectivity profiles may provide a valuable biomarker for the diagnosis of mTBI.

Identification of Voxels Confounded by Venous Signals Using Resting-State fMRI Functional Connectivity Graph Community Identification.

Identifying venous voxels in fMRI datasets is important to increase the specificity of fMRI analyses to microvasculature in the vicinity of the neural processes triggering the BOLD response. This is, however, difficult to achieve in particular in typical studies where magnitude images of BOLD EPI are the only data available. In this study, voxelwise functional connectivity graphs were computed on minimally preprocessed low TR (333 ms) multiband resting-state fMRI data, using both high positive and negative correlations to define edges between nodes (voxels). A high correlation threshold for binarization ensures that most edges in the resulting sparse graph reflect the high coherence of signals in medium to large veins. Graph clustering based on the optimization of modularity was then employed to identify clusters of coherent voxels in this graph, and all clusters of 50 or more voxels were then interpreted as corresponding to medium to large veins. Indeed, a comparison with SWI reveals that 75.6±5.9% of voxels within these large clusters overlap with veins visible in the SWI image or lie outside the brain parenchyma. Some of the remaining differences between the two modalities can be explained by imperfect alignment or geometric distortions between the two images. Overall, the graph clustering based method for identifying venous voxels has a high specificity as well as the additional advantages of being computed in the same voxel grid as the fMRI dataset itself and not needing any additional data beyond what is usually acquired (and exported) in standard fMRI experiments.

Graph theoretical comparison of functional connectivity between cLTP treated and untreated microelectrode arrays

Analyzing graph properties of neural networks has recently gained much attention in attempts to understand how information is processed in the brain. Using in-vitro techniques to form neural networks has increased in popularity as it allows one to develop small, easy to record networks that maintain many of the graph properties of larger brain networks [1]. One widely recognized tool for studying in vitro networks is the Microelectrode Array (MEAs) on which neurons can be cultured and recorded simultaneously. MEAs can be used to grow neural networks from disassociated cells to understand how neurons spontaneously connect to create networks and how these networks then evolve over time. In addition, these cultures can be treated with pharmacological agents to study how these agents affect the networks as a whole [2,3]. To understand the network formation of MEA cultured neurons, we study the graph theoretical properties of two MEAs networks, the control MEA network and the MEA network treated with chemical Long Term Potentiation (cLTP). The data sets for each MEA network consists of recording from three days: baseline, 2 days past baseline and 5 days past baseline. Based on these data sets and the assumption that each electrode on the MEA records one neuron, we construct functional connectivity graphs of MEA networks for different days. Nodes in such a connectivity graph represent the electrodes (also neurons). To determine whether there is a connection (an edge on the graph) between two nodes, we carry out several steps of computations. We first filter the recorded spike trains with a Gaussian kernel, and perform cross-correlation analysis using the Pearson product moment correlation coefficient [4]. We set a correlation threshold by applying a shuffling method to the inter-spike intervals of a spike train. Using thresholded correlations, unweighted, undirected adjacency matrices, we create corresponding graphs for untreated (not shown) and treated MEA networks (shown baseline and 5 days past baseline in Figure 1). We find that the synchronization and average node degree increase dramatically for the cLTP treated networks while the untreated network shows no obvious change. Figure 1 Graph models. (a) cLTP treated MEA network at baseline; (b) cLTP treated MEA network at 5 days past baseline. To better understand the treated and untreated MEA network, we will evaluate the graph theoretic properties, such as degree distribution and clustering coefficient. We will determine how cLTP affects these properties. The graphical analysis will enable us to identify what type of network each is (such as a small-world or a scale free network) and determine whether cLTP has an effect on the network development or merely on the strength of connectivity. We conjecture that cLTP treated networks have more efficient and quicker communication between nodes. Therefore, the cLTP treated networks show greater clustering as well as shorter path length than the untreated networks. Information flow is another important aspect of such graph model. We intend to develop directed graphs using transfer entropy to study how information flow of the network may change during its development.

A wavelet-based estimator of the degrees of freedom in denoised fMRI time series for probabilistic testing of functional connectivity and brain graphs

Connectome mapping using techniques such as functional magnetic resonance imaging (fMRI) has become a focus of systems neuroscience. There remain many statistical challenges in analysis of functional connectivity and network architecture from BOLD fMRI multivariate time series. One key statistic for any time series is its (effective) degrees of freedom, df, which will generally be less than the number of time points (or nominal degrees of freedom, N). If we know the df, then probabilistic inference on other fMRI statistics, such as the correlation between two voxel or regional time series, is feasible. However, we currently lack good estimators of df in fMRI time series, especially after the degrees of freedom of the “raw” data have been modified substantially by denoising algorithms for head movement. Here, we used a wavelet-based method both to denoise fMRI data and to estimate the (effective) df of the denoised process. We show that seed voxel correlations corrected for locally variable df could be tested for false positive connectivity with better control over Type I error and greater specificity of anatomical mapping than probabilistic connectivity maps using the nominal degrees of freedom. We also show that wavelet despiked statistics can be used to estimate all pairwise correlations between a set of regional nodes, assign a P value to each edge, and then iteratively add edges to the graph in order of increasing P. These probabilistically thresholded graphs are likely more robust to regional variation in head movement effects than comparable graphs constructed by thresholding correlations. Finally, we show that time-windowed estimates of df can be used for probabilistic connectivity testing or dynamic network analysis so that apparent changes in the functional connectome are appropriately corrected for the effects of transient noise bursts. Wavelet despiking is both an algorithm for fMRI time series denoising and an estimator of the (effective) df of denoised fMRI time series. Accurate estimation of df offers many potential advantages for probabilistically thresholding functional connectivity and network statistics tested in the context of spatially variant and non-stationary noise. Code for wavelet despiking, seed correlational testing and probabilistic graph construction is freely available to download as part of the BrainWavelet Toolbox at www.brainwavelet.org.

On the danger of detecting network states in white noise.

On the danger of detecting network states in white noise.

Functional connectivity changes detected with magnetoencephalography after mild traumatic brain injury.

Mild traumatic brain injury (mTBI) may affect normal cognition and behavior by disrupting the functional connectivity networks that mediate efficient communication among brain regions. In this study, we analyzed brain connectivity profiles from resting state Magnetoencephalographic (MEG) recordings obtained from 31 mTBI patients and 55 normal controls. We used phase-locking value estimates to compute functional connectivity graphs to quantify frequency-specific couplings between sensors at various frequency bands. Overall, normal controls showed a dense network of strong local connections and a limited number of long-range connections that accounted for approximately 20% of all connections, whereas mTBI patients showed networks characterized by weak local connections and strong long-range connections that accounted for more than 60% of all connections. Comparison of the two distinct general patterns at different frequencies using a tensor representation for the connectivity graphs and tensor subspace analysis for optimal feature extraction showed that mTBI patients could be separated from normal controls with 100% classification accuracy in the alpha band. These encouraging findings support the hypothesis that MEG-based functional connectivity patterns may be used as biomarkers that can provide more accurate diagnoses, help guide treatment, and monitor effectiveness of intervention in mTBI.

Functional connectivity estimation over large networks at cellular resolution based on electrophysiological recordings and structural prior.

Despite many structural and functional aspects of the brain organization have been extensively studied in neuroscience, we are still far from a clear understanding of the intricate structure-function interactions occurring in the multi-layered brain architecture, where billions of different neurons are involved. Although structure and function can individually convey a large amount of information, only a combined study of these two aspects can probably shade light on how brain circuits develop and operate at the cellular scale. Here, we propose a novel approach for refining functional connectivity estimates within neuronal networks using the structural connectivity as prior. This is done at the mesoscale, dealing with thousands of neurons while reaching, at the microscale, an unprecedented cellular resolution. The High-Density Micro Electrode Array (HD-MEA) technology, combined with fluorescence microscopy, offers the unique opportunity to acquire structural and functional data from large neuronal cultures approaching the granularity of the single cell. In this work, an advanced method based on probabilistic directional features and heat propagation is introduced to estimate the structural connectivity from the fluorescence image while functional connectivity graphs are obtained from the cross-correlation analysis of the spiking activity. Structural and functional information are then integrated by reweighting the functional connectivity graph based on the structural prior. Results show that the resulting functional connectivity estimates are more coherent with the network topology, as compared to standard measures purely based on cross-correlations and spatio-temporal filters. We finally use the obtained results to gain some insights on which features of the functional activity are more relevant to characterize actual neuronal interactions.

Fast construction of voxel-level functional connectivity graphs.

BackgroundGraph-based analysis of fMRI data has recently emerged as a promising approach to study brain networks. Based on the assessment of synchronous fMRI activity at separate brain sites, functional connectivity graphs are constructed and analyzed using graph-theoretical concepts. Most previous studies investigated region-level graphs, which are computationally inexpensive, but bring along the problem of choosing sensible regions and involve blurring of more detailed information. In contrast, voxel-level graphs provide the finest granularity attainable from the data, enabling analyses at superior spatial resolution. They are, however, associated with considerable computational demands, which can render high-resolution analyses infeasible. In response, many existing studies investigating functional connectivity at the voxel-level reduced the computational burden by sacrificing spatial resolution.MethodsHere, a novel, time-efficient method for graph construction is presented that retains the original spatial resolution. Performance gains are instead achieved through data reduction in the temporal domain based on dichotomization of voxel time series combined with tetrachoric correlation estimation and efficient implementation.ResultsBy comparison with graph construction based on Pearson’s r, the technique used by the majority of previous studies, we find that the novel approach produces highly similar results an order of magnitude faster.ConclusionsIts demonstrated performance makes the proposed approach a sensible and efficient alternative to customary practice. An open source software package containing the created programs is freely available for download.

Subject-specific functional parcellation via Prior Based Eigenanatomy

We present a new framework for prior-constrained sparse decomposition of matrices derived from the neuroimaging data and apply this method to functional network analysis of a clinically relevant population. Matrix decomposition methods are powerful dimensionality reduction tools that have found widespread use in neuroimaging. However, the unconstrained nature of these totally data-driven techniques makes it difficult to interpret the results in a domain where network-specific hypotheses may exist.We propose a novel approach, Prior Based Eigenanatomy (p-Eigen), which seeks to identify a data-driven matrix decomposition but at the same time constrains the individual components by spatial anatomical priors (probabilistic ROIs). We formulate our novel solution in terms of prior-constrained ℓ1 penalized (sparse) principal component analysis. p-Eigen starts with a common functional parcellation for all the subjects and refines it with subject-specific information. This enables modeling of the inter-subject variability in the functional parcel boundaries and allows us to construct subject-specific networks with reduced sensitivity to ROI placement.We show that while still maintaining correspondence across subjects, p-Eigen extracts biologically-relevant and patient-specific functional parcels that facilitate hypothesis-driven network analysis. We construct default mode network (DMN) connectivity graphs using p-Eigen refined ROIs and use them in a classification paradigm. Our results show that the functional connectivity graphs derived from p-Eigen significantly aid classification of mild cognitive impairment (MCI) as well as the prediction of scores in a Delayed Recall memory task when compared to graph metrics derived from 1) standard registration-based seed ROI definitions, 2) totally data-driven ROIs, 3) a model based on standard demographics plus hippocampal volume as covariates, and 4) Ward Clustering based data-driven ROIs. In summary, p-Eigen incarnates a new class of prior-constrained dimensionality reduction tools that may improve our understanding of the relationship between MCI and functional connectivity.

Local Functional Overconnectivity in Posterior Brain Regions Is Associated with Symptom Severity in Autism Spectrum Disorders

Although growing evidence indicates atypical long-distance connectivity in autism spectrum disorder (ASD), much less is known about local connectivity, despite conjectures that local overconnectivity may be causally involved in the disorder. Using functional connectivity MRI and graph theory, we found that local functional connectivity was atypically increased in adolescents with ASD in temporo-occipital regions bilaterally. Posterior overconnectivity was found to be associated with higher ASD symptom severity, whereas an ASD subsample with low severity showed frontal underconnectivity. The findings suggest links between symptomatology and local connectivity, which vary within the autism spectrum.

Detection of preclinical neural dysfunction from functional connectivity graphs derived from task fMRI. An example from degeneration

The early, preferably pre-clinical, identification of neurodegenerative diseases is important as treatment will be most successful before substantial neuronal loss. Here, we reasoned that functional brain changes as measured using functional magnetic resonance imaging (fMRI) will precede neurodegeneration. Three independent cohorts of patients with the genetic mutation leading to Huntington's Disease (HD) but without any clinical symptoms and matched controls performed three different fMRI tasks: Sequential finger tapping engaged the motor system, which is primarily affected by HD, whereas a working-memory task and a task aiming to induce irritation represented the range of low- and high-level cognitive functions also affected by HD. Each diagnostic group of every cohort included 11–16 subjects. After segmentation into 76 cortical and 14 subcortical regions, we extracted functional connectivity patterns through pairwise correlation between the signals in the regions. The resulting coefficients were directly embedded as input to a pattern classifier aiming to separate controls from gene mutation carriers. Alternatively, graph-theory measures such as degree and clustering coefficient were used as features for the discrimination. Classification accuracy never outperformed the accuracy of a grouping based on parameter estimates from a general-linear model approach or a grouping based on features extracted from anatomical images as reported in a previous analysis. Despite good within-subject stability between two runs of the same task, a high between-subject variability led to chance-level accuracy. These results indicate that standard graph-metrics are insufficient to detect subtle disease related changes when within-group variability is high. Developing methods that reduce variability related to noise should be the focus of future studies.

Neurophysiology of juvenile myoclonic epilepsy: EEG-based network and graph analysis of the interictal and immediate preictal states

The neuronal mechanisms of enduring seizure propensity and seizure precipitation in juvenile myoclonic epilepsy (JME) are not known. We investigated these issues, within the framework of the "network concept" of epilepsy. Design1: 19, unmedicated JME patients were compared with nineteen, age-, and sex-matched normal control persons (NC). A total of 120s, artifact-free, paroxysm-free, eyes-closed, resting state EEG background activity was analyzed for each person. Design2: interictal and immediate preictal periods of the JME patients were compared in order to explore interictal-preictal network differences. For both comparison designs, statistically significant differences of EEG functional connectivity (EEGfC), nodal and global graph parameters were evaluated. Design1: maximum abnormalities were: increased delta, theta, alpha1 EEGfC and decreased alpha2 and beta EEGfC in the JME group as compared to the NC group, mainly among cortical areas that are involved in sensory-motor integration. Nodal degree and efficiency of three, medial, basal frontal nodes were greater in JME than in NC, in the alpha1 band. Design2: preictal delta EEGfC showed further increase in the above-mentioned areas, as compared to the interictal state. Increased EEGfC indicates a hypercoupled state among the specified cortical areas. This interictal abnormality further increases in the preictal state. Nodal graph statistics indicates abnormal neuronal dynamics in the cortical area that is the ictal onset zone in JME. Interictal and preictal neuronal dysfunction has been described in terms of network dynamics and topography in JME patients. Forthcoming investigations of seizure precipitation and therapeutic drug effects are encouraged on this basis.

Large-scale brain functional modularity is reflected in slow electroencephalographic rhythms across the human non-rapid eye movement sleep cycle

Large-scale brain functional networks (measured with functional magnetic resonance imaging, fMRI) are organized into separated but interacting modules, an architecture supporting the integration of distinct dynamical processes. In this work we study how the aforementioned modular architecture changes with the progressive loss of vigilance occurring in the descent to deep sleep and we examine the relationship between the ensuing slow electroencephalographic rhythms and large-scale network modularity as measured with fMRI. Graph theoretical methods are used to analyze functional connectivity graphs obtained from fifty-five participants at wakefulness, light and deep sleep. Network modularity (a measure of functional segregation) was found to increase during deeper sleep stages but not in light sleep. By endowing functional networks with dynamical properties, we found a direct link between increased electroencephalographic (EEG) delta power (1-4 Hz) and a breakdown of inter-modular connectivity. Both EEG slowing and increased network modularity were found to quickly decrease during awakenings from deep sleep to wakefulness, in a highly coordinated fashion. Studying the modular structure itself by means of a permutation test, we revealed different module memberships when deep sleep was compared to wakefulness. Analysis of node roles in the modular structure revealed an increase in the number of locally well-connected nodes and a decrease in the number of globally well-connected hubs, which hinders interactions between separated functional modules. Our results reveal a well-defined sequence of changes in brain modular organization occurring during the descent to sleep and establish a close parallel between modularity alterations in large-scale functional networks (accessible through whole brain fMRI recordings) and the slowing of scalp oscillations (visible on EEG). The observed re-arrangement of connectivity might play an important role in the processes underlying loss of vigilance and sensory awareness during deep sleep.

On the discovery of group-consistent graph substructure patterns from brain networks

Complex networks constitute a recurring issue in the analysis of neuroimaging data. Recently, network motifs have been identified as patterns of interconnections since they appear in a significantly higher number than in randomized networks, in a given ensemble of anatomical or functional connectivity graphs. The current approach for detecting and enumerating motifs in brain networks requires a predetermined motif repertoire and can operate only with motifs of small size (consisting of few nodes).There is a growing interest in methodologies for frequent graph-based pattern mining in large graph datasets that can facilitate adaptive design of motifs. The results presented in this paper are based on the graph-based Substructure pattern mining (gSpan) algorithm and introduce a manifold of ways to exploit it for data-driven motif extraction in connectomics research.Functional connectivity graphs from electroencephalographic (EEG) recordings during resting state and mental calculations are used to demonstrate our approach. Relying on either time-invariant or time-evolving graphs, characteristic motifs associated with various frequency bands were derived and compared. With a suitable manipulation, the gSpan discovers motifs which are specific to performing mental arithmetics. Finally, the subject-dependent temporal signatures of motifs’ appearance revealed the transient nature of the evolving functional connectivity (math-related motifs “come and go”).