Fractal connectivity of long-memory networks

Resting-state fractal brain connectivity is associated with impaired cognitive performance in healthy aging

Aging affects cognitive functions even in the absence of ongoing pathologies. The neurophysiological basis of age-related cognitive decline (CD), however, is not completely understood. Alterations in both functional brain connectivity and in the fractal scaling of neuronal dynamics have been linked to aging and cognitive performance. Recently, fractal connectivity (FrC) has been proposed - combining the two concepts - for capturing long-term interactions among brain regions. FrC was shown to be influenced by increased mental workload; however, no prior studies investigated how resting-state FrC relates to cognitive performance and plausible CD in healthy aging. We recruited 19 healthy elderly (HE) and 24 young control (YC) participants, who underwent resting-state electroencephalography (EEG) measurements and comprehensive cognitive evaluation using 7 tests of the Cambridge Neurophysiological Test Automated Battery. FrC networks were reconstructed from EEG data using the recently introduced multiple-resampling cross-spectral analysis (MRCSA). Elderly individuals could be characterized with increased response latency and reduced performance in 4-4 tasks, respectively, with both reaction time and accuracy being affected in two tasks. Auto- and cross-spectral exponents - characterizing regional fractal dynamics and FrC, respectively, - were found reduced in HE when compared to YC over most of the cortex. Additionally, fractal scaling of frontoparietal connections expressed an inverse relationship with task performance in visual memory and sustained attention domains in elderly, but not in young individuals. Our results confirm that the fractal nature of brain connectivity - as captured by MRCSA - is affected in healthy aging. Furthermore, FrC appears as a sensitive neurophysiological marker of age-related CD.

Fingerprints of decreased cognitive performance on fractal connectivity dynamics in healthy aging

Analysis of brain functional connectivity (FC) could provide insight in how and why cognitive functions decline even in healthy aging (HA). Despite FC being established as fluctuating over time even in the resting state (RS), dynamic functional connectivity (DFC) studies involving healthy elderly individuals and assessing how these patterns relate to cognitive performance are yet scarce. In our recent study we showed that fractal temporal scaling of functional connections in RS is not only reduced in HA, but also predicts increased response latency and reduced task solving accuracy. However, in that work we did not address changes in the dynamics of fractal connectivity (FrC) strength itself and its plausible relationship with mental capabilities. Therefore, here we analyzed RS electroencephalography recordings of the same subject cohort as previously, consisting of 24 young and 19 healthy elderly individuals, who also completed 7 different cognitive tasks after data collection. Dynamic fractal connectivity (dFrC) analysis was carried out via sliding-window detrended cross-correlation analysis (DCCA). A machine learning method based on recursive feature elimination was employed to select the subset of connections most discriminative between the two age groups, identifying 56 connections that allowed for classifying participants with an accuracy surpassing 92%. Mean of DCCA was found generally increased, while temporal variability of FrC decreased in the elderly when compared to the young group. Finally, dFrC indices expressed an elaborate pattern of associations-assessed via Spearman correlation-with cognitive performance scores in both groups, linking fractal connectivity strength and variance to increased response latency and reduced accuracy in the elderly population. Our results provide further support for the relevance of FrC dynamics in understanding age-related cognitive decline and might help to identify potential targets for future intervention strategies.

Functional brain networks reflect spatial and temporal autocorrelation

High-throughput experimental methods in neuroscience have led to an explosion of techniques for measuring complex interactions and multi-dimensional patterns. However, whether sophisticated measures of emergent phenomena can be traced back to simpler, low-dimensional statistics is largely unknown. To explore this question, we examined resting-state functional magnetic resonance imaging (rs-fMRI) data using complex topology measures from network neuroscience. Here we show that spatial and temporal autocorrelation are reliable statistics that explain numerous measures of network topology. Surrogate time series with subject-matched spatial and temporal autocorrelation capture nearly all reliable individual and regional variation in these topology measures. Network topology changes during aging are driven by spatial autocorrelation, and multiple serotonergic drugs causally induce the same topographic change in temporal autocorrelation. This reductionistic interpretation of widely used complexity measures may help link them to neurobiology.

Brain Network Analysis: A Review on Multivariate Analytical Methods

Despite the explosive growth of neuroimaging studies aimed at analyzing the brain as a complex system, critical methodological gaps remain to be addressed. Most tools currently used for analyzing network data of the brain are univariate in nature and are based on assumptions borne out of previous techniques not directly related to the big and complex data of the brain. Although graph-based methods have shown great promise, the development of principled multivariate models to address inherent limitations of graph-based methods, such as their dependence on network size and degree distributions, and to allow assessing the effects of multiple phenotypes on the brain and simulating brain networks has largely lagged behind. Although some studies have been made in developing multivariate frameworks to fill this gap, in the absence of a "gold-standard" method or guidelines, choosing the most appropriate method for each study can be another critical challenge for investigators in this multidisciplinary field. Here, we briefly introduce important multivariate methods for brain network analyses in two main categories: data-driven and model-based methods. We discuss whether/how such methods are suited for examining connectivity (edge-level), topology (system-level), or both. This review will aid in choosing an appropriate multivariate method with respect to variables such as network type, number of subjects and brain regions included, and the interest in connectivity, topology, or both. This review is aimed to be accessible to investigators from different backgrounds, with a focus on applications in brain network studies, though the methods may be applicable in other areas too. Impact statement As the U.S. National Institute of Health notes, the rich biomedical data can greatly improve our knowledge of human health if new analytical tools are developed, and their applications are broadly disseminated. A major challenge in analyzing the brain as a complex system is about developing parsimonious multivariate methods, and particularly choosing the most appropriate one among the existing methods with respect to the study variables in this multidisciplinary field. This study provides a review on the most important multivariate methods to aid in helping the most appropriate ones with respect to the desired variables for each study.

Wavelet-Based Fractal Analysis of rs-fMRI for Classification of Alzheimer's Disease

The resting-state functional magnetic resonance imaging (rs-fMRI) modality has gained widespread acceptance as a promising method for analyzing a variety of neurological and psychiatric diseases. It is established that resting-state neuroimaging data exhibit fractal behavior, manifested in the form of slow-decaying auto-correlation and power-law scaling of the power spectrum across low-frequency components. With this property, the rs-fMRI signal can be broken down into fractal and nonfractal components. The fractal nature originates from several sources, such as cardiac fluctuations, respiration and system noise, and carries no information on the brain’s neuronal activities. As a result, the conventional correlation of rs-fMRI signals may not accurately reflect the functional dynamic of spontaneous neuronal activities. This problem can be solved by using a better representation of neuronal activities provided by the connectivity of nonfractal components. In this work, the nonfractal connectivity of rs-fMRI is used to distinguish Alzheimer’s patients from healthy controls. The automated anatomical labeling (AAL) atlas is used to extract the blood-oxygenation-level-dependent time series signals from 116 brain regions, yielding a 116 × 116 nonfractal connectivity matrix. From this matrix, significant connections evaluated using the p-value are selected as an input to a classifier for the classification of Alzheimer’s vs. normal controls. The nonfractal-based approach provides a good representation of the brain’s neuronal activity. It outperformed the fractal and Pearson-based connectivity approaches by 16.4% and 17.2%, respectively. The classification algorithm developed based on the nonfractal connectivity feature and support vector machine classifier has shown an excellent performance, with an accuracy of 90.3% and 83.3% for the XHSLF dataset and ADNI dataset, respectively. For further validation of our proposed work, we combined the two datasets (XHSLF+ADNI) and still received an accuracy of 90.2%. The proposed work outperformed the recently published work by a margin of 8.18% and 11.2%, respectively.

Multiple-Resampling Cross-Spectral Analysis: An Unbiased Tool for Estimating Fractal Connectivity With an Application to Neurophysiological Signals

Investigating scale-free (i.e., fractal) functional connectivity in the brain has recently attracted increasing attention. Although numerous methods have been developed to assess the fractal nature of functional coupling, these typically ignore that neurophysiological signals are assemblies of broadband, arrhythmic activities as well as oscillatory activities at characteristic frequencies such as the alpha waves. While contribution of such rhythmic components may bias estimates of fractal connectivity, they are also likely to represent neural activity and coupling emerging from distinct mechanisms. Irregular-resampling auto-spectral analysis (IRASA) was recently introduced as a tool to separate fractal and oscillatory components in the power spectrum of neurophysiological signals by statistically summarizing the power spectra obtained when resampling the original signal by several non-integer factors. Here we introduce multiple-resampling cross-spectral analysis (MRCSA) as an extension of IRASA from the univariate to the bivariate case, namely, to separate the fractal component of the cross-spectrum between two simultaneously recorded neural signals by applying the same principle. MRCSA does not only provide a theoretically unbiased estimate of the fractal cross-spectrum (and thus its spectral exponent) but also allows for computing the proportion of scale-free coupling between brain regions. As a demonstration, we apply MRCSA to human electroencephalographic recordings obtained in a word generation paradigm. We show that the cross-spectral exponent as well as the proportion of fractal coupling increases almost uniformly over the cortex during the rest-task transition, likely reflecting neural desynchronization. Our results indicate that MRCSA can be a valuable tool for scale-free connectivity studies in characterizing various cognitive states, while it also can be generalized to other applications outside the field of neuroscience.

Constructing Multi-view High-order Functional Connectivity Networks for Diagnosis of Autism Spectrum Disorder

Brain functional connectivity network (FCN) based on resting-state functional magnetic resonance imaging (rs-fMRI) has been widely used to identify neuropsychiatric disorders such as autism spectrum disorder (ASD). Most existing FCN-based methods only estimate the correlation between brain regions of in terest (ROIs), without exploring more informative higher-level inter actions among multiple ROIs which could be beneficial to disease diagnosis. To fully explore the discriminative information provided by different brain networks, a cluster-based multi-view high-order FCN (Ho-FCN) framework is proposed in this paper. Specifically, we first group the functional connectivity (FC) time series into different clusters and compute the multi-order central moment series for the FC time series in each cluster. Then we utilize the correlation of central moment series between different clusters to reveal the high-order FC relationships among multiple ROIs. In addition, to address the phase mismatch issue in conventional FCNs, we also adopt the central moments of the correlation time series as the temporal-invariance features to capture the dynamic characteristics of low-order dynamic FCN (Lo-D-FCN). Experimental results on the ABIDE dataset validate that: 1) the proposed multi-view Ho-FCNs is able to explore rich discriminative information for ASD diagnosis; 2) the phase mismatch issue can be well circumvented by using central moments; and 3) the combination of different types of FCNs can significantly improve the diagnostic accuracy of ASD (86.2%).

Multifractal Functional Connectivity Analysis of Electroencephalogram Reveals Reorganization of Brain Networks in a Visual Pattern Recognition Paradigm

The human brain consists of anatomically distant neuronal assemblies that are interconnected via a myriad of synapses. This anatomical network provides the neurophysiological wiring framework for functional connectivity (FC), which is essential for higher-order brain functions. While several studies have explored the scale-specific FC, the scale-free (i.e., multifractal) aspect of brain connectivity remains largely neglected. Here we examined the brain reorganization during a visual pattern recognition paradigm, using bivariate focus-based multifractal (BFMF) analysis. For this study, 58 young, healthy volunteers were recruited. Before the task, 3-3 min of resting EEG was recorded in eyes-closed (EC) and eyes-open (EO) states, respectively. The subsequent part of the measurement protocol consisted of 30 visual pattern recognition trials of 3 difficulty levels graded as Easy, Medium, and Hard. Multifractal FC was estimated with BFMF analysis of preprocessed EEG signals yielding two generalized Hurst exponent-based multifractal connectivity endpoint parameters, H(2) and ΔH ₁₅; with the former indicating the long-term cross-correlation between two brain regions, while the latter captures the degree of multifractality of their functional coupling. Accordingly, H(2) and ΔH ₁₅ networks were constructed for every participant and state, and they were characterized by their weighted local and global node degrees. Then, we investigated the between- and within-state variability of multifractal FC, as well as the relationship between global node degree and task performance captured in average success rate and reaction time. Multifractal FC increased when visual pattern recognition was administered with no differences regarding difficulty level. The observed regional heterogeneity was greater for ΔH ₁₅ networks compared to H(2) networks. These results show that reorganization of scale-free coupled dynamics takes place during visual pattern recognition independent of difficulty level. Additionally, the observed regional variability illustrates that multifractal FC is region-specific both during rest and task. Our findings indicate that investigating multifractal FC under various conditions - such as mental workload in healthy and potentially in diseased populations - is a promising direction for future research.

Constructing high-order functional connectivity network based on central moment features for diagnosis of autism spectrum disorder

The sliding-window-based dynamic functional connectivity network (D-FCN) has been becoming an increasingly useful tool for understanding the changes of brain connectivity patterns and the association of neurological diseases with these dynamic variations. However, conventional D-FCN is essentially low-order network, which only reflects the pairwise interaction pattern between brain regions and thus overlooking the high-order interactions among multiple brain regions. In addition, D-FCN is innate with temporal sensitivity issue, i.e., D-FCN is sensitive to the chronological order of its subnetworks. To deal with the above issues, we propose a novel high-order functional connectivity network framework based on the central moment feature of D-FCN. Specifically, we firstly adopt a central moment approach to extract multiple central moment feature matrices from D-FCN. Furthermore, we regard the matrices as the profiles to build multiple high-order functional connectivity networks which further capture the higher level and more complex interaction relationships among multiple brain regions. Finally, we use the voting strategy to combine the high-order networks with D-FCN for autism spectrum disorder diagnosis. Experimental results show that the combination of multiple functional connectivity networks achieves accuracy of 88.06%, and the best single network achieves accuracy of 79.5%.

Scale-Free Coupled Dynamics in Brain Networks Captured by Bivariate Focus-Based Multifractal Analysis

While most connectivity studies investigate functional connectivity (FC) in a scale-dependent manner, coupled neural processes may also exhibit broadband dynamics, manifesting as power-law scaling of their measures of interdependence. Here we introduce the bivariate focus-based multifractal (BFMF) analysis as a robust tool for capturing such scale-free relations and use resting-state electroencephalography (EEG) recordings of 12 subjects to demonstrate its performance in reconstructing physiological networks. BFMF was employed to characterize broadband FC between 62 cortical regions in a pairwise manner, with all investigated connections being tested for true bivariate multifractality. EEG channels were also grouped to represent the activity of six resting-state networks (RSNs) in the brain, thus allowing for the analysis of within- and between- RSNs connectivity, separately. Most connections featured true bivariate multifractality, which could be attributed to the genuine scale-free coupling of neural dynamics. Bivariate multifractality showed a characteristic topology over the cortex that was highly concordant among subjects. Long-term autocorrelation was higher in within-RSNs, while the degree of multifractality was generally found stronger in between-RSNs connections. These results offer statistical evidence of the bivariate multifractal nature of functional coupling in the brain and validate BFMF as a robust method to capture such scale-independent coupled dynamics.

Revisiting Functional Connectivity for Infraslow Scale-Free Brain Dynamics Using Complex Wavelets

The analysis of human brain functional networks is achieved by computing functional connectivity indices reflecting phase coupling and interactions between remote brain regions. In magneto- and electroencephalography, the most frequently used functional connectivity indices are constructed based on Fourier-based cross-spectral estimation applied to specific fast and band-limited oscillatory regimes. Recently, infraslow arrhythmic fluctuations (below the 1 Hz) were recognized as playing a leading role in spontaneous brain activity. The present work aims to propose to assess functional connectivity from fractal dynamics, thus extending the assessment of functional connectivity to the infraslow arrhythmic or scale-free temporal dynamics of M/EEG-quantified brain activity. Instead of being based on Fourier analysis, new Imaginary Coherence and weighted Phase Lag indices are constructed from complex-wavelet representations. Their performances are first assessed on synthetic data by means of Monte-Carlo simulations, and they are then compared favorably against the classical Fourier-based indices. These new assessments of functional connectivity indices are also applied to MEG data collected on 36 individuals both at rest and during the learning of a visual motion discrimination task. They demonstrate a higher statistical sensitivity, compared to their Fourier counterparts, in capturing significant and relevant functional interactions in the infraslow regime and modulations from rest to task. Notably, the consistent overall increase in functional connectivity assessed from fractal dynamics from rest to task correlated with a change in temporal dynamics as well as with improved performance in task completion, which suggests that the complex-wavelet weighted Phase Lag index is the sole index is able to capture brain plasticity in the infraslow scale-free regime.

Diagnosis of Autism Spectrum Disorder Using Central-Moment Features From Low- and High-Order Dynamic Resting-State Functional Connectivity Networks

The sliding-window-based dynamic functional connectivity networks (D-FCNs) derived from resting-state functional magnetic resonance imaging (rs-fMRI) are effective methods for diagnosing various neurological diseases, including autism spectrum disorder (ASD). However, traditional D-FCNs are low-order networks based on pairwise correlation between brain regions, thus overlooking high-level interactions across multiple regions of interest (ROIs). Moreover, D-FCNs suffer from the temporal mismatching issue, i.e., subnetworks in the same temporal window do not have temporal correspondence across different subjects. To address the above problems, we first construct a novel high-order D-FCNs based on the principle of “correlation’s correlation” to further explore the higher level and more complex interaction relationships among multiple ROIs. Furthermore, we propose to use a central-moment method to extract temporal-invariance properties contained in either low- or high-order D-FCNs. Finally, we design and train an ensemble classifier by fusing the features extracted from conventional FCN, low-order D-FCNs, and high-order D-FCNs for the diagnosis of ASD and normal control subjects. Our method achieved the best ASD classification accuracy (83%), and our results revealed the features extracted from different networks fingerprinting the autistic brain at different connectional levels.

More Flexible Integration of Functional Systems After Musical Training in Young Adults

Musical training, because it involves the interaction and integration of diverse functional systems, is an excellent model to investigate training-induced brain plasticity. The human brain functions in a network architecture in which dynamic modules and subgraphs are considered to enable efficient information communication. However, it remains largely unknown how the dynamic integration of functional systems changes with musical training, which may provide new insight into musical training-induced brain plasticity and further the use of music therapy for neuropsychiatric disease and brain injury. Here, 29 healthy young adult novices who received 24 weeks of piano training, and another 27 novices without any intervention were scanned at three time points-before and after musical training and 12 weeks after training. We used nonnegative matrix factorization to identify a set of subgraphs and their corresponding time-dependent coefficients from a concatenated functional network of all the subjects in sliding time windows. The energy and entropy of the time-dependent coefficients were computed to quantify the subgraph's dynamic changes in expression. The training group showed a significantly increased energy of the time-dependent coefficients of 3 subgraphs after training. Furthermore, one of the subgraphs, comprised of primary functional systems and cingulo-opercular task control and salience systems, showed significantly changed entropy in the training group after training. Our results suggest that the integration of functional systems undergoes increased flexibility in fine-scale dynamics after musical training, which reveals how brain functional systems engage in musical performance. The efficacy of musical training induced brain plasticity may provide new therapeutic strategies for brain injury and neuropsychiatric disorders.

Brain Networks Reveal the Effects of Antipsychotic Drugs on Schizophrenia Patients and Controls

The study of brain networks, including those derived from functional neuroimaging data, attracts a broad interest and represents a rapidly growing interdisciplinary field. Comparing networks of healthy volunteers with those of patients can potentially offer new, quantitative diagnostic methods and a framework for better understanding brain and mind disorders. We explore resting state functional Magnetic Resonance Imaging (fMRI) data through network measures. We construct networks representing 15 healthy individuals and 12 schizophrenia patients (males and females), all of whom are administered three drug treatments: i) a placebo; and two antipsychotic medications ii) aripiprazole and iii) sulpiride. We compare these resting state networks to a performance at an "N-back" working memory task. We demonstrate that not only is there a distinctive network architecture in the healthy brain that is disrupted in schizophrenia but also that both networks respond to antipsychotic medication. We first reproduce the established finding that brain networks of schizophrenia patients exhibit increased efficiency and reduced clustering compared with controls. Our data then reveal that the antipsychotic medications mitigate this effect, shifting the metrics toward those observed in healthy volunteers, with a marked difference in efficacy between the two drugs. Additionally, we find that aripiprazole considerably alters the network statistics of healthy controls. Examining the "N-back" working memory task, we establish that aripiprazole also adversely affects their performance. This suggests that changes to macroscopic brain network architecture result in measurable behavioral differences. This is one of the first studies to directly compare different medications using a whole-brain graph theoretical analysis with accompanying behavioral data. The small sample size is an inherent limitation and means a degree of caution is warranted in interpreting the findings. Our results lay the groundwork for an objective methodology with which to calculate and compare the efficacy of different treatments of mind and brain disorders.

Multi-scale detection of hierarchical community architecture in structural and functional brain networks

Community detection algorithms have been widely used to study the organization of complex networks like the brain. These techniques provide a partition of brain regions (or nodes) into clusters (or communities), where nodes within a community are densely interconnected with one another. In their simplest application, community detection algorithms are agnostic to the presence of community hierarchies: clusters embedded within clusters of other clusters. To address this limitation, we exercise a multi-scale extension of a common community detection technique, and we apply the tool to synthetic graphs and to graphs derived from human neuroimaging data, including structural and functional imaging data. Our multi-scale community detection algorithm links a graph to copies of itself across neighboring topological scales, thereby becoming sensitive to conserved community organization across levels of the hierarchy. We demonstrate that this method is sensitive to topological inhomogeneities of the graph's hierarchy by providing a local measure of community stability and inter-scale reliability across topological scales. We compare the brain's structural and functional network architectures, and we demonstrate that structural graphs display a more prominent hierarchical community organization than functional graphs. Finally, we build an explicitly multimodal multiplex graph that combines both structural and functional connectivity in a single model, and we identify the topological scales where resting state functional connectivity and underlying structural connectivity show similar versus unique hierarchical community architecture. Together, our results demonstrate the advantages of the multi-scale community detection algorithm in studying hierarchical community structure in brain graphs, and they illustrate its utility in modeling multimodal neuroimaging data.

Subnetwork mining on functional connectivity network for classification of minimal hepatic encephalopathy

Hepatic encephalopathy (HE), as a complication of cirrhosis, is a serious brain disease, which may lead to death. Accurate diagnosis of HE and its intermediate stage, i.e., minimal HE (MHE), is very important for possibly early diagnosis and treatment. Brain connectivity network, as a simple representation of brain interaction, has been widely used for the brain disease (e.g., HE and MHE) analysis. However, those studies mainly focus on finding disease-related abnormal connectivity between brain regions, although a large number of studies have indicated that some brain diseases are usually related to local structure of brain connectivity network (i.e., subnetwork), rather than solely on some single brain regions or connectivities. Also, mining such disease-related subnetwork is a challenging task because of the complexity of brain network. To address this problem, we proposed a novel frequent-subnetwork-based method to mine disease-related subnetworks for MHE classification. Specifically, we first mine frequent subnetworks from both groups, i.e., MHE patients and non-HE (NHE) patients, respectively. Then we used the graph-kernel based method to select the most discriminative subnetworks for subsequent classification. We evaluate our proposed method on a MHE dataset with 77 cirrhosis patients, including 38 MHE patients and 39 NHE patients. The results demonstrate that our proposed method can not only obtain the improved classification performance in comparison with state-of-the-art network-based methods, but also identify disease-related subnetworks which can help us better understand the pathology of the brain diseases.

Modular reconfiguration of an auditory control brain network supports adaptive listening behavior

How do brain networks shape our listening behavior? We here develop and test the hypothesis that, during challenging listening situations, intrinsic brain networks are reconfigured to adapt to the listening demands and thus, to enable successful listening. We find that, relative to a task-free resting state, networks of the listening brain show higher segregation of temporal auditory, ventral attention, and frontal control regions known to be involved in speech processing, sound localization, and effortful listening. Importantly, the relative change in modularity of this auditory control network predicts individuals’ listening success. Our findings shed light on how cortical communication dynamics tune selection and comprehension of speech in challenging listening situations and suggest modularity as the network principle of auditory attention.

Speech comprehension in noisy, multitalker situations poses a challenge. Successful behavioral adaptation to a listening challenge often requires stronger engagement of auditory spatial attention and context-dependent semantic predictions. Human listeners differ substantially in the degree to which they adapt behaviorally and can listen successfully under such circumstances. How cortical networks embody this adaptation, particularly at the individual level, is currently unknown. We here explain this adaptation from reconfiguration of brain networks for a challenging listening task (i.e., a linguistic variant of the Posner paradigm with concurrent speech) in an age-varying sample of n = 49 healthy adults undergoing resting-state and task fMRI. We here provide evidence for the hypothesis that more successful listeners exhibit stronger task-specific reconfiguration (hence, better adaptation) of brain networks. From rest to task, brain networks become reconfigured toward more localized cortical processing characterized by higher topological segregation. This reconfiguration is dominated by the functional division of an auditory and a cingulo-opercular module and the emergence of a conjoined auditory and ventral attention module along bilateral middle and posterior temporal cortices. Supporting our hypothesis, the degree to which modularity of this frontotemporal auditory control network is increased relative to resting state predicts individuals’ listening success in states of divided and selective attention. Our findings elucidate how fine-tuned cortical communication dynamics shape selection and comprehension of speech. Our results highlight modularity of the auditory control network as a key organizational principle in cortical implementation of auditory spatial attention in challenging listening situations.

Diagnosis of Autism Spectrum Disorders Using Multi-Level High-Order Functional Networks Derived From Resting-State Functional MRI

Functional brain networks derived from resting-state functional magnetic resonance imaging (rs-fMRI) have been widely used for Autism Spectrum Disorder (ASD) diagnosis. Typically, these networks are constructed by calculating functional connectivity (FC) between any pair of brain regions of interest (ROIs), i.e., using Pearson's correlation between rs-fMRI time series. However, this can only be called as a low-order representation of the functional interaction, because the relationship is investigated just between two ROIs. Brain disorders might not only affect low-order FC, but also high-order FC, i.e., the higher-level relationship among multiple brain regions, which might be more crucial for diagnosis. To comprehensively characterize such relationship for better diagnosis of ASD, we propose a multi-level, high-order FC network representation that can nicely capture complex interactions among brain regions. Then, we design a feature selection method to identify those discriminative multi-level, high-order FC features for ASD diagnosis. Finally, we design an ensemble classifier with multiple linear SVMs, each trained on a specific level of FC networks, for boosting the final classification accuracy. Experimental results show that the integration of both low-order and first-level high-order FC networks achieves the best ASD diagnostic accuracy (81%). We further investigated those selected discriminative low-order and high-order FC features and found that the high-order FC features can provide complementary information to the low-order FC features in the ASD diagnosis.

Time scale properties of task and resting-state functional connectivity: Detrended partial cross-correlation analysis

Functional connectivity analysis is an essential tool for understanding brain function. Previous studies showed that brain regions are functionally connected through low-frequency signals both within the default mode network (DMN) and task networks. However, no studies have directly compared the time scale (frequency) properties of network connectivity during task versus rest, or examined how they relate to task performance. Here, using fMRI data collected from sixty-eight subjects at rest and during a stop signal task, we addressed this issue with a novel functional connectivity measure based on detrended partial cross-correlation analysis (DPCCA). DPCCA has the advantage of quantifying correlations between two variables in different time scales while controlling for the influence of other variables. The results showed that the time scales of within-network connectivity of the DMN and task networks are modulated in opposite directions across rest and task, with the time scale increased during rest vs. task in the DMN and vice versa in task networks. In regions of interest analysis, the within-network connectivity time scale of the pre-supplementary motor area - a medial prefrontal cortical structure of the task network and critical to proactive inhibitory control - correlated inversely with Barratt impulsivity and stop signal reaction time. Together, these findings demonstrate that time scale properties of brain networks may vary across mental states and provide evidence in support of a role of low frequency fluctuations of BOLD signals in behavioral control.

Brain state flexibility accompanies motor-skill acquisition

Learning requires the traversal of inherently distinct cognitive states to produce behavioral adaptation. Yet, tools to explicitly measure these states with non-invasive imaging – and to assess their dynamics during learning – remain limited. Here, we describe an approach based on a distinct application of graph theory in which points in time are represented by network nodes, and similarities in brain states between two different time points are represented as network edges. We use a graph-based clustering technique to identify clusters of time points representing canonical brain states, and to assess the manner in which the brain moves from one state to another as learning progresses. We observe the presence of two primary states characterized by either high activation in sensorimotor cortex or high activation in a frontal-subcortical system. Flexible switching among these primary states and other less common states becomes more frequent as learning progresses, and is inversely correlated with individual differences in learning rate. These results are consistent with the notion that the development of automaticity is associated with a greater freedom to use cognitive resources for other processes. Taken together, our work offers new insights into the constrained, low dimensional nature of brain dynamics characteristic of early learning, which give way to less constrained, high-dimensional dynamics in later learning.

Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks

Human functional magnetic resonance imaging (fMRI) brain networks have a complex topology comprising integrative components, e.g. long-distance inter-modular edges, that are theoretically associated with higher biological cost. Here, we estimated intra-modular degree, inter-modular degree and connection distance for each of 285 cortical nodes in multi-echo fMRI data from 38 healthy adults. We used the multivariate technique of partial least squares (PLS) to reduce the dimensionality of the relationships between these three nodal network parameters and prior microarray data on regional expression of 20 737 genes. The first PLS component defined a transcriptional profile associated with high intra-modular degree and short connection distance, whereas the second PLS component was associated with high inter-modular degree and long connection distance. Nodes in superior and lateral cortex with high inter-modular degree and long connection distance had local transcriptional profiles enriched for oxidative metabolism and mitochondria, and for genes specific to supragranular layers of human cortex. In contrast, primary and secondary sensory cortical nodes in posterior cortex with high intra-modular degree and short connection distance had transcriptional profiles enriched for RNA translation and nuclear components. We conclude that, as predicted, topologically integrative hubs, mediating long-distance connections between modules, are more costly in terms of mitochondrial glucose metabolism.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Feature fusion via hierarchical supervised local CCA for diagnosis of autism spectrum disorder

Early diagnosis of autism spectrum disorder (ASD) is critical for timely medical intervention, for improving patient quality of life, and for reducing the financial burden borne by the society. A key issue in neuroimaging-based ASD diagnosis is the identification of discriminating features and then fusing them to produce accurate diagnosis. In this paper, we propose a novel framework for fusing complementary and discriminating features from different imaging modalities. Specifically, we integrate the Fisher discriminant criterion and local correlation information into the canonical correlation analysis (CCA) framework, giving a new feature fusion method, called Supervised Local CCA (SL-CCA), which caters specifically to local and global multimodal features. To alleviate the neighborhood selection problem associated with SL-CCA, we further propose a hierarchical SL-CCA (HSL-CCA), by performing SL-CCA with the gradually varying neighborhood sizes. Extensive experiments on the multimodal ABIDE database show that the proposed method achieves superior performance. In addition, based on feature weight analysis, we found that only a few specific brain regions play active roles in ASD diagnosis. These brain regions include the putamen, precuneus, and orbitofrontal cortex, which are highly associated with human emotional modulation and memory formation. These finding are consistent with the behavioral phenotype of ASD.

Functional hypergraph uncovers novel covariant structures over neurodevelopment

Brain development during adolescence is marked by substantial changes in brain structure and function, leading to a stable network topology in adulthood. However, most prior work has examined the data through the lens of brain areas connected to one another in large-scale functional networks. Here, we apply a recently developed hypergraph approach that treats network connections (edges) rather than brain regions as the unit of interest, allowing us to describe functional network topology from a fundamentally different perspective. Capitalizing on a sample of 780 youth imaged as part of the Philadelphia Neurodevelopmental Cohort, this hypergraph representation of resting-state functional MRI data reveals three distinct classes of subnetworks (hyperedges): clusters, bridges, and stars, which respectively represent homogeneously connected, bipartite, and focal architectures. Cluster hyperedges show a strong resemblance to previously-described functional modules of the brain including somatomotor, visual, default mode, and salience systems. In contrast, star hyperedges represent highly localized subnetworks centered on a small set of regions, and are distributed across the entire cortex. Finally, bridge hyperedges link clusters and stars in a core-periphery organization. Notably, developmental changes within hyperedges are ordered in a similar core-periphery fashion, with the greatest developmental effects occurring in networked hyperedges within the functional core. Taken together, these results reveal a novel decomposition of the network organization of human brain, and further provide a new perspective on the role of local structures that emerge across neurodevelopment. Hum Brain Mapp 38:3823-3835, 2017. © 2017 Wiley Periodicals, Inc.

Cohesive network reconfiguration accompanies extended training

Human behavior is supported by flexible neurophysiological processes that enable the fine‐scale manipulation of information across distributed neural circuits. Yet, approaches for understanding the dynamics of these circuit interactions have been limited. One promising avenue for quantifying and describing these dynamics lies in multilayer network models. Here, networks are composed of nodes (which represent brain regions) and time‐dependent edges (which represent statistical similarities in activity time series). We use this approach to examine functional connectivity measured by non‐invasive neuroimaging techniques. These multilayer network models facilitate the examination of changes in the pattern of statistical interactions between large‐scale brain regions that might facilitate behavior. In this study, we define and exercise two novel measures of network reconfiguration, and demonstrate their utility in neuroimaging data acquired as healthy adult human subjects learn a new motor skill. In particular, we identify putative functional modules in multilayer networks and characterize the degree to which nodes switch between modules. Next, we define cohesive switches, in which a set of nodes moves between modules together as a group, and we define disjoint switches, in which a single node moves between modules independently from other nodes. Together, these two concepts offer complementary yet distinct insights into the changes in functional connectivity that accompany motor learning. More generally, our work offers statistical tools that other researchers can use to better understand the reconfiguration patterns of functional connectivity over time. Hum Brain Mapp 38:4744–4759, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Self-similarity and quasi-idempotence in neural networks and related dynamical systems

Self-similarity across length scales is pervasively observed in natural systems. Here, we investigate topological self-similarity in complex networks representing diverse forms of connectivity in the brain and some related dynamical systems, by considering the correlation between edges directly connecting any two nodes in a network and indirect connection between the same via all triangles spanning the rest of the network. We note that this aspect of self-similarity, which is distinct from hierarchically nested connectivity (coarse-grain similarity), is closely related to idempotence of the matrix representing the graph. We introduce two measures, ι(1) and ι(∞), which represent the element-wise correlation coefficients between the initial matrix and the ones obtained after squaring it once or infinitely many times, and term the matrices which yield large values of these parameters "quasi-idempotent". These measures delineate qualitatively different forms of "shallow" and "deep" quasi-idempotence, which are influenced by nodal strength heterogeneity. A high degree of quasi-idempotence was observed for partially synchronized mean-field Kuramoto oscillators with noise, electronic chaotic oscillators, and cultures of dissociated neurons, wherein the expression of quasi-idempotence correlated strongly with network maturity. Quasi-idempotence was also detected for macro-scale brain networks representing axonal connectivity, synchronization of slow activity fluctuations during idleness, and co-activation across experimental tasks, and preliminary data indicated that quasi-idempotence of structural connectivity may decrease with ageing. This initial study highlights that the form of network self-similarity indexed by quasi-idempotence is detectable in diverse dynamical systems, and draws attention to it as a possible basis for measures representing network "collectivity" and pattern formation.

Interhemispheric connectivity in amyotrophic lateral sclerosis: A near-infrared spectroscopy and diffusion tensor imaging study

PURPOSE

Aim of the present study was to investigate potential impairment of non-motor areas in amyotrophic lateral sclerosis (ALS) using near-infrared spectroscopy (NIRS) and diffusion tensor imaging (DTI). In particular, we evaluated whether homotopic resting-state functional connectivity (rs-FC) of non-motor associated cortical areas correlates with clinical parameters and disease-specific degeneration of the corpus callosum (CC) in ALS.

MATERIAL AND METHODS

Interhemispheric homotopic rs-FC was assessed in 31 patients and 30 healthy controls (HCs) for 8 cortical sites, from prefrontal to occipital cortex, using NIRS. DTI was performed in a subgroup of 21 patients. All patients were evaluated for cognitive dysfunction in the executive, memory, and visuospatial domains.

RESULTS

ALS patients displayed an altered spatial pattern of correlation between homotopic rs-FC values when compared to HCs (p = 0.000013). In patients without executive dysfunction a strong correlation existed between the rate of motor decline and homotopic rs-FC of the anterior temporal lobes (ATLs) (ρ = - 0.85, p = 0.0004). Furthermore, antero-temporal homotopic rs-FC correlated with fractional anisotropy in the central corpus callosum (CC), corticospinal tracts (CSTs), and forceps minor as determined by DTI (p < 0.05).

CONCLUSIONS

The present study further supports involvement of non-motor areas in ALS. Our results render homotopic rs-FC as assessed by NIRS a potential clinical marker for disease progression rate in ALS patients without executive dysfunction and a potential anatomical marker for ALS-specific degeneration of the CC and CSTs.

Choosing Wavelet Methods, Filters, and Lengths for Functional Brain Network Construction

Wavelet methods are widely used to decompose fMRI, EEG, or MEG signals into time series representing neurophysiological activity in fixed frequency bands. Using these time series, one can estimate frequency-band specific functional connectivity between sensors or regions of interest, and thereby construct functional brain networks that can be examined from a graph theoretic perspective. Despite their common use, however, practical guidelines for the choice of wavelet method, filter, and length have remained largely undelineated. Here, we explicitly explore the effects of wavelet method (MODWT vs. DWT), wavelet filter (Daubechies Extremal Phase, Daubechies Least Asymmetric, and Coiflet families), and wavelet length (2 to 24)—each essential parameters in wavelet-based methods—on the estimated values of graph metrics and in their sensitivity to alterations in psychiatric disease. We observe that the MODWT method produces less variable estimates than the DWT method. We also observe that the length of the wavelet filter chosen has a greater impact on the estimated values of graph metrics than the type of wavelet chosen. Furthermore, wavelet length impacts the sensitivity of the method to detect differences between health and disease and tunes classification accuracy. Collectively, our results suggest that the choice of wavelet method and length significantly alters the reliability and sensitivity of these methods in estimating values of metrics drawn from graph theory. They furthermore demonstrate the importance of reporting the choices utilized in neuroimaging studies and support the utility of exploring wavelet parameters to maximize classification accuracy in the development of biomarkers of psychiatric disease and neurological disorders.

Robust Statistical Detection of Power-Law Cross-Correlation

We show that widely used approaches in statistical physics incorrectly indicate the existence of power-law cross-correlations between financial stock market fluctuations measured over several years and the neuronal activity of the human brain lasting for only a few minutes. While such cross-correlations are nonsensical, no current methodology allows them to be reliably discarded, leaving researchers at greater risk when the spurious nature of cross-correlations is not clear from the unrelated origin of the time series and rather requires careful statistical estimation. Here we propose a theory and method (PLCC-test) which allows us to rigorously and robustly test for power-law cross-correlations, correctly detecting genuine and discarding spurious cross-correlations, thus establishing meaningful relationships between processes in complex physical systems. Our method reveals for the first time the presence of power-law cross-correlations between amplitudes of the alpha and beta frequency ranges of the human electroencephalogram.

Multi-Dimensional Dynamics of Human Electromagnetic Brain Activity

Magnetoencephalography (MEG) and electroencephalography (EEG) are invaluable neuroscientific tools for unveiling human neural dynamics in three dimensions (space, time, and frequency), which are associated with a wide variety of perceptions, cognition, and actions. MEG/EEG also provides different categories of neuronal indices including activity magnitude, connectivity, and network properties along the three dimensions. In the last 20 years, interest has increased in inter-regional connectivity and complex network properties assessed by various sophisticated scientific analyses. We herein review the definition, computation, short history, and pros and cons of connectivity and complex network (graph-theory) analyses applied to MEG/EEG signals. We briefly describe recent developments in source reconstruction algorithms essential for source-space connectivity and network analyses. Furthermore, we discuss a relatively novel approach used in MEG/EEG studies to examine the complex dynamics represented by human brain activity. The correct and effective use of these neuronal metrics provides a new insight into the multi-dimensional dynamics of the neural representations of various functions in the complex human brain.

Exploring multifractal-based features for mild Alzheimer's disease classification

PURPOSE

Multifractal applications to resting state functional MRI (rs-fMRI) time series for diagnosing Alzheimer's disease (AD) are still limited. We aim to address two issues: (I) if and what multifractal features are sufficiently discriminative to detect AD from the healthy; (II) if AD classification could be further improved by combining multifractal features with traditional features in this field.

METHODS

Rs-fMRI data of 25 AD patients and 38 normal controls were analyzed. A set of multifractal features were systematically investigated. Traditional features in monofractal, linear, and network-based categories were also extracted for comparison and combination. Both support vector machines and multiple kernel learning (MKL) were used to perform classification with individual and combined features.

RESULTS

We identified a multifractal feature, Δf, which has the strongest discriminative power among all the features in our study. Moreover, we found that the classification accuracy could be significantly improved from 69% (by Δf only) to up to 76%, when nonsparse MKL is used to combine Δf with the monofractal feature, Hurst. Finally, we showed that incorporating other multifractal features, α(0), Δα and Pc, could also improve traditional-feature-based AD classification.

CONCLUSION

Our work demonstrated the potential usefulness of multifractal analysis for AD research, especially when combining with the traditional rs-fMRI features. It contributes to distinguishing AD from NC subjects. Magn Reson Med 76:259-269, 2016. © 2015 Wiley Periodicals, Inc.

Detrended partial cross-correlation analysis of two nonstationary time series influenced by common external forces

When common factors strongly influence two power-law cross-correlated time series recorded in complex natural or social systems, using detrended cross-correlation analysis (DCCA) without considering these common factors will bias the results. We use detrended partial cross-correlation analysis (DPXA) to uncover the intrinsic power-law cross correlations between two simultaneously recorded time series in the presence of nonstationarity after removing the effects of other time series acting as common forces. The DPXA method is a generalization of the detrended cross-correlation analysis that takes into account partial correlation analysis. We demonstrate the method by using bivariate fractional Brownian motions contaminated with a fractional Brownian motion. We find that the DPXA is able to recover the analytical cross Hurst indices, and thus the multiscale DPXA coefficients are a viable alternative to the conventional cross-correlation coefficient. We demonstrate the advantage of the DPXA coefficients over the DCCA coefficients by analyzing contaminated bivariate fractional Brownian motions. We calculate the DPXA coefficients and use them to extract the intrinsic cross correlation between crude oil and gold futures by taking into consideration the impact of the U.S. dollar index. We develop the multifractal DPXA (MF-DPXA) method in order to generalize the DPXA method and investigate multifractal time series. We analyze multifractal binomial measures masked with strong white noises and find that the MF-DPXA method quantifies the hidden multifractal nature while the multifractal DCCA method fails.

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.

Scale invariance properties of intracerebral EEG improve seizure prediction in mesial temporal lobe epilepsy

Although treatment for epilepsy is available and effective for nearly 70 percent of patients, many remain in need of new therapeutic approaches. Predicting the impending seizures in these patients could significantly enhance their quality of life if the prediction performance is clinically practical. In this study, we investigate the improvement of the performance of a seizure prediction algorithm in 17 patients with mesial temporal lobe epilepsy by means of a novel measure. Scale-free dynamics of the intracerebral EEG are quantified through robust estimates of the scaling exponents--the first cumulants--derived from a wavelet leader and bootstrap based multifractal analysis. The cumulants are investigated for the discriminability between preictal and interictal epochs. The performance of our recently published patient-specific seizure prediction algorithm is then out-of-sample tested on long-lasting data using combinations of cumulants and state similarity measures previously introduced. By using the first cumulant in combination with state similarity measures, up to 13 of 17 patients had seizures predicted above chance with clinically practical levels of sensitivity (80.5%) and specificity (25.1% of total time under warning) for prediction horizons above 25 min. These results indicate that the scale-free dynamics of the preictal state are different from those of the interictal state. Quantifiers of these dynamics may carry a predictive power that can be used to improve seizure prediction performance.

Learning-induced autonomy of sensorimotor systems

Distributed networks of brain areas interact with one another in a time-varying fashion to enable complex cognitive and sensorimotor functions. Here we used new network-analysis algorithms to test the recruitment and integration of large-scale functional neural circuitry during learning. Using functional magnetic resonance imaging data acquired from healthy human participants, we investigated changes in the architecture of functional connectivity patterns that promote learning from initial training through mastery of a simple motor skill. Our results show that learning induces an autonomy of sensorimotor systems and that the release of cognitive control hubs in frontal and cingulate cortices predicts individual differences in the rate of learning on other days of practice. Our general statistical approach is applicable across other cognitive domains and provides a key to understanding time-resolved interactions between distributed neural circuits that enable task performance.

Frequent and discriminative subnetwork mining for mild cognitive impairment classification

Recent studies on brain networks have suggested that many brain diseases, such as Alzheimer's disease and mild cognitive impairment (MCI), are related to a large-scale brain network, rather than individual brain regions. However, it is challenging to find such a network from the whole brain network due to the complexity of brain networks. In this article, the authors propose a novel method to mine the discriminative subnetworks for classifying MCI patients from healthy controls (HC). Specifically, the authors first extract a set of frequent subnetworks from each of the two groups (i.e., MCI and HC), respectively. Then, measure the discriminative ability of those frequent subnetworks using the graph kernel-based classification method and select the most discriminative subnetworks for subsequent classification. The results on the functional connectivity networks of 12 MCI and 25 HC show that this method can obtain competitive results compared with state-of-the-art methods on MCI classification.

Brain Network Organization in Focal Epilepsy: A Systematic Review and Meta-Analysis

Normal brain functioning is presumed to depend upon interacting regions within large-scale neuronal networks. Increasing evidence exists that interictal network alterations in focal epilepsy are associated with cognitive and behavioral deficits. Nevertheless, the reported network alterations are inconclusive and prone to low statistical power due to small sample sizes as well as modest effect sizes. We therefore systematically reviewed the existing literature and conducted a meta-analysis to characterize the changes in whole-brain interictal focal epilepsy networks at sufficient power levels. We focused on the two most commonly used metrics in whole-brain networks: average path length and average clustering coefficient. Twelve studies were included that reported whole-brain network average path length and average clustering coefficient characteristics in patients and controls. The overall group difference, quantified as the standardized mean average path length difference between epilepsy and control groups, corresponded to a significantly increased average path length of 0.29 (95% confidence interval (CI): 0.12 to 0.45, p = 0.0007) in the epilepsy group. This suggests a less integrated interictal whole-brain network. Similarly, a significantly increased standardized mean average clustering coefficient of 0.35 (CI: 0.05 to 0.65, p = 0.02) was found in the epilepsy group in comparison with controls, pointing towards a more segregated interictal network. Sub-analyses revealed similar results for functional and structural networks in terms of effect size and directionality for both metrics. In addition, we found individual network studies to be prone to low power due to the relatively small group differences in average path length and average clustering coefficient in combination with small sample sizes. The pooled network characteristics support the hypothesis that focal epilepsy has widespread detrimental effects, that is, reduced integration and increased segregation, on whole brain interictal network organization, which may relate to the co-morbid cognitive and behavioral impairments often reported in patients with focal epilepsy.

Resting-state fMRI activity predicts unsupervised learning and memory in an immersive virtual reality environment

In the real world, learning often proceeds in an unsupervised manner without explicit instructions or feedback. In this study, we employed an experimental paradigm in which subjects explored an immersive virtual reality environment on each of two days. On day 1, subjects implicitly learned the location of 39 objects in an unsupervised fashion. On day 2, the locations of some of the objects were changed, and object location recall performance was assessed and found to vary across subjects. As prior work had shown that functional magnetic resonance imaging (fMRI) measures of resting-state brain activity can predict various measures of brain performance across individuals, we examined whether resting-state fMRI measures could be used to predict object location recall performance. We found a significant correlation between performance and the variability of the resting-state fMRI signal in the basal ganglia, hippocampus, amygdala, thalamus, insula, and regions in the frontal and temporal lobes, regions important for spatial exploration, learning, memory, and decision making. In addition, performance was significantly correlated with resting-state fMRI connectivity between the left caudate and the right fusiform gyrus, lateral occipital complex, and superior temporal gyrus. Given the basal ganglia's role in exploration, these findings suggest that tighter integration of the brain systems responsible for exploration and visuospatial processing may be critical for learning in a complex environment.

A wavelet method for modeling and despiking motion artifacts from resting-state fMRI time series

The impact of in-scanner head movement on functional magnetic resonance imaging (fMRI) signals has long been established as undesirable. These effects have been traditionally corrected by methods such as linear regression of head movement parameters. However, a number of recent independent studies have demonstrated that these techniques are insufficient to remove motion confounds, and that even small movements can spuriously bias estimates of functional connectivity. Here we propose a new data-driven, spatially-adaptive, wavelet-based method for identifying, modeling, and removing non-stationary events in fMRI time series, caused by head movement, without the need for data scrubbing. This method involves the addition of just one extra step, the Wavelet Despike, in standard pre-processing pipelines. With this method, we demonstrate robust removal of a range of different motion artifacts and motion-related biases including distance-dependent connectivity artifacts, at a group and single-subject level, using a range of previously published and new diagnostic measures. The Wavelet Despike is able to accommodate the substantial spatial and temporal heterogeneity of motion artifacts and can consequently remove a range of high and low frequency artifacts from fMRI time series, that may be linearly or non-linearly related to physical movements. Our methods are demonstrated by the analysis of three cohorts of resting-state fMRI data, including two high-motion datasets: a previously published dataset on children (N=22) and a new dataset on adults with stimulant drug dependence (N=40). We conclude that there is a real risk of motion-related bias in connectivity analysis of fMRI data, but that this risk is generally manageable, by effective time series denoising strategies designed to attenuate synchronized signal transients induced by abrupt head movements. The Wavelet Despiking software described in this article is freely available for download at www.brainwavelet.org.

Effects of transient unilateral functional brain disruption on global neural network status in rats: a methods paper

Permanent focal brain damage can have critical effects on the function of nearby as well as remote brain regions. However, the effects of transient disturbances on global brain function are largely unknown. Our goal was to develop an experimental in vivo model to map the impact of transient functional brain impairment on large-scale neural networks in the absence of structural damage. We describe a new rat model of transient functional hemispheric disruption using unilateral focal anesthesia by intracarotid pentobarbital injection. The brain's functional status was assessed with resting-state fMRI (rs-fMRI) and electroencephalography (EEG). We performed network analysis to identify and quantify highly connected network hubs, i.e., "rich-club organization," in pre- and postbarbital functional networks. Perfusion MRI data demonstrated that the catheterized carotid artery predominantly supplied the ipsilateral hemisphere, allowing for selective hemispheric brain silencing. The prebarbital baseline network displayed strong functional connectivity (FC) within and between hemispheres. Following pentobarbital injection, the disrupted hemisphere revealed increased intrahemispheric FC with concomitant decrease of interhemispheric connectivity. The bilateral functional network was characterized by a strong positive rich-club effect, which was not affected by ipsilateral disruption. Nevertheless, the rich-club value was significantly decreased in the ipsilateral hemisphere and to a lesser extent contralaterally. Loss of interhemispheric EEG synchronization supported the rs-fMRI findings. Our data support the concept that densely connected rich-club regions play a central role in global brain communication, and show that network hub configurations can be significantly affected by focal temporary functional hemispheric disruption without structural neuronal damage. Further studies with this rat model will provide essential additional insights into network reorganization patterns in response to transient functional brain disruption.

Group-constrained sparse fMRI connectivity modeling for mild cognitive impairment identification

Emergence of advanced network analysis techniques utilizing resting-state functional magnetic resonance imaging (R-fMRI) has enabled a more comprehensive understanding of neurological disorders at a whole-brain level. However, inferring brain connectivity from R-fMRI is a challenging task, particularly when the ultimate goal is to achieve good control-patient classification performance, owing to perplexing noise effects, curse of dimensionality, and inter-subject variability. Incorporating sparsity into connectivity modeling may be a possible solution to partially remedy this problem since most biological networks are intrinsically sparse. Nevertheless, sparsity constraint, when applied at an individual level, will inevitably cause inter-subject variability and hence degrade classification performance. To this end, we formulate the R-fMRI time series of each region of interest (ROI) as a linear representation of time series of other ROIs to infer sparse connectivity networks that are topologically identical across individuals. This formulation allows simultaneous selection of a common set of ROIs across subjects so that their linear combination is best in estimating the time series of the considered ROI. Specifically, l 1-norm is imposed on each subject to filter out spurious or insignificant connections to produce sparse networks. A group-constraint is hence imposed via multi-task learning using a l 2-norm to encourage consistent non-zero connections across subjects. This group-constraint is crucial since the network topology is identical for all subjects while still preserving individual information via different connectivity values. We validated the proposed modeling in mild cognitive impairment identification and promising results achieved demonstrate its superiority in disease characterization, particularly greater sensitivity to early stage brain pathologies. The inferred group-constrained sparse network is found to be biologically plausible and is highly associated with the disease-associated anatomical anomalies. Furthermore, our proposed approach achieved similar classification performance when finer atlas was used to parcellate the brain space.

Integration of network topological and connectivity properties for neuroimaging classification

Rapid advances in neuroimaging techniques have provided an efficient and noninvasive way for exploring the structural and functional connectivity of the human brain. Quantitative measurement of abnormality of brain connectivity in patients with neurodegenerative diseases, such as mild cognitive impairment (MCI) and Alzheimer's disease (AD), have also been widely reported, especially at a group level. Recently, machine learning techniques have been applied to the study of AD and MCI, i.e., to identify the individuals with AD/MCI from the healthy controls (HCs). However, most existing methods focus on using only a single property of a connectivity network, although multiple network properties, such as local connectivity and global topological properties, can potentially be used. In this paper, by employing multikernel based approach, we propose a novel connectivity based framework to integrate multiple properties of connectivity network for improving the classification performance. Specifically, two different types of kernels (i.e., vector-based kernel and graph kernel) are used to quantify two different yet complementary properties of the network, i.e., local connectivity and global topological properties. Then, multikernel learning (MKL) technique is adopted to fuse these heterogeneous kernels for neuroimaging classification. We test the performance of our proposed method on two different data sets. First, we test it on the functional connectivity networks of 12 MCI and 25 HC subjects. The results show that our method achieves significant performance improvement over those using only one type of network property. Specifically, our method achieves a classification accuracy of 91.9%, which is 10.8% better than those by single network-property-based methods. Then, we test our method for gender classification on a large set of functional connectivity networks with 133 infants scanned at birth, 1 year, and 2 years, also demonstrating very promising results.

Inhibition of propofol anesthesia on functional connectivity between LFPs in PFC during rat working memory task

Working memory (WM) refers to the temporary storage and manipulation of information necessary for performance of complex cognitive tasks. There is a growing interest in whether and how propofol anesthesia inhibits WM function. The aim of this study is to investigate the possible inhibition mechanism of propofol anesthesia based on the functional connections of multi-local field potentials (LFPs) and behavior during WM tasks. Adult SD rats were randomly divided into 3 groups: pro group (0.5 mg·kg⁻¹·min⁻¹,2 h), PRO group (0.9 mg·kg⁻¹·min⁻¹, 2 h) and control group. The experimental data were 16-channel LFPs obtained at prefrontal cortex with implanted microelectrode array in SD rats during WM tasks in Y-maze at 24, 48, 72, 96, 120 hours (day 1-day 5) after propofol anesthesia, and the behavior results of WM were recoded at the same time. Directed transfer function (DTF) method was applied to analyze the connections among LFPs directly. Furthermore, the causal networks were identified by DTF. The clustering coefficient (C), network density (D) and global efficiency (E_global) were selected to describe the functional connectivity quantitatively. The results show that: comparing with the control group, the LFPs functional connectivity in pro group were no significantly difference (p>0.05); the connectivity in PRO group were significantly decreased (p<0.05 at 24 hours, p<0.05 at 48 hours), while no significant difference at 72, 96 and 120 hours for rats (p>0.05), which were consistent with the behavior results. These findings could lead to improved understanding the mechanism of inhibition of anesthesia on WM functions from the view of connections among LFPs.

Task-based core-periphery organization of human brain dynamics

As a person learns a new skill, distinct synapses, brain regions, and circuits are engaged and change over time. In this paper, we develop methods to examine patterns of correlated activity across a large set of brain regions. Our goal is to identify properties that enable robust learning of a motor skill. We measure brain activity during motor sequencing and characterize network properties based on coherent activity between brain regions. Using recently developed algorithms to detect time-evolving communities, we find that the complex reconfiguration patterns of the brain's putative functional modules that control learning can be described parsimoniously by the combined presence of a relatively stiff temporal core that is composed primarily of sensorimotor and visual regions whose connectivity changes little in time and a flexible temporal periphery that is composed primarily of multimodal association regions whose connectivity changes frequently. The separation between temporal core and periphery changes over the course of training and, importantly, is a good predictor of individual differences in learning success. The core of dynamically stiff regions exhibits dense connectivity, which is consistent with notions of core-periphery organization established previously in social networks. Our results demonstrate that core-periphery organization provides an insightful way to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

When someone learns a new skill, his/her brain dynamically alters individual synapses, regional activity, and larger-scale circuits. In this paper, we capture some of these dynamics by measuring and characterizing patterns of coherent brain activity during the learning of a motor skill. We extract time-evolving communities from these patterns and find that a temporal core that is composed primarily of primary sensorimotor and visual regions reconfigures little over time, whereas a periphery that is composed primarily of multimodal association regions reconfigures frequently. The core consists of densely connected nodes, and the periphery consists of sparsely connected nodes. Individual participants with a larger separation between core and periphery learn better in subsequent training sessions than individuals with a smaller separation. Conceptually, core-periphery organization provides a framework in which to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

Topological graph kernel on multiple thresholded functional connectivity networks for mild cognitive impairment classification

Recently, brain connectivity networks have been used for classification of Alzheimer's disease and mild cognitive impairment (MCI) from normal controls (NC). In typical connectivity-networks-based classification approaches, local measures of connectivity networks are first extracted from each region-of-interest as network features, which are then concatenated into a vector for subsequent feature selection and classification. However, some useful structural information of network, especially global topological information, may be lost in this type of approaches. To address this issue, in this article, we propose a connectivity-networks-based classification framework to identify accurately the MCI patients from NC. The core of the proposed method involves the use of a new graph-kernel-based approach to measure directly the topological similarity between connectivity networks. We evaluate our method on functional connectivity networks of 12 MCI and 25 NC subjects. The experimental results show that our proposed method achieves a classification accuracy of 91.9%, a sensitivity of 100.0%, a balanced accuracy of 94.0%, and an area under receiver operating characteristic curve of 0.94, demonstrating a great potential in MCI classification, based on connectivity networks. Further connectivity analysis indicates that the connectivity of the selected brain regions is different between MCI patients and NC, that is, MCI patients show reduced functional connectivity compared with NC, in line with the findings reported in the existing studies.