Analytical Operations Relate Structural and Functional Connectivity in the Brain

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

The Hopf whole-brain model and its linear approximation

Whole-brain models have proven to be useful to understand the emergence of collective activity among neural populations or brain regions. These models combine connectivity matrices, or connectomes, with local node dynamics, noise, and, eventually, transmission delays. Multiple choices for the local dynamics have been proposed. Among them, nonlinear oscillators corresponding to a supercritical Hopf bifurcation have been used to link brain connectivity and collective phase and amplitude dynamics in different brain states. Here, we studied the linear fluctuations of this model to estimate its stationary statistics, i.e., the instantaneous and lagged covariances and the power spectral densities. This linear approximation-that holds in the case of heterogeneous parameters and time-delays-allows analytical estimation of the statistics and it can be used for fast parameter explorations to study changes in brain state, changes in brain activity due to alterations in structural connectivity, and modulations of parameter due to non-equilibrium dynamics.

Alzheimer's Disease: Insights from Large-Scale Brain Dynamics Models

Alzheimer's disease (AD) is a degenerative brain disease, and the condition is difficult to assess. In the past, numerous brain dynamics models have made remarkable contributions to neuroscience and the brain from the microcosmic to the macroscopic scale. Recently, large-scale brain dynamics models have been developed based on dual-driven multimodal neuroimaging data and neurodynamics theory. These models bridge the gap between anatomical structure and functional dynamics and have played an important role in assisting the understanding of the brain mechanism. Large-scale brain dynamics have been widely used to explain how macroscale neuroimaging biomarkers emerge from potential neuronal population level disturbances associated with AD. In this review, we describe this emerging approach to studying AD that utilizes a biophysically large-scale brain dynamics model. In particular, we focus on the application of the model to AD and discuss important directions for the future development and analysis of AD models. This will facilitate the development of virtual brain models in the field of AD diagnosis and treatment and add new opportunities for advancing clinical neuroscience.

A joint subspace mapping between structural and functional brain connectomes

Understanding the connection between the brain's structural connectivity and its functional connectivity is of immense interest in computational neuroscience. Although some studies have suggested that whole brain functional connectivity is shaped by the underlying structure, the rule by which anatomy constraints brain dynamics remains an open question. In this work, we introduce a computational framework that identifies a joint subspace of eigenmodes for both functional and structural connectomes. We found that a small number of those eigenmodes are sufficient to reconstruct functional connectivity from the structural connectome, thus serving as low-dimensional basis function set. We then develop an algorithm that can estimate the functional eigen spectrum in this joint space from the structural eigen spectrum. By concurrently estimating the joint eigenmodes and the functional eigen spectrum, we can reconstruct a given subject's functional connectivity from their structural connectome. We perform elaborate experiments and demonstrate that the proposed algorithm for estimating functional connectivity from the structural connectome using joint space eigenmodes gives competitive performance as compared to the existing benchmark methods with better interpretability.

Data-driven discovery of canonical large-scale brain dynamics

Human behavior and cognitive function correlate with complex patterns of spatio-temporal brain dynamics, which can be simulated using computational models with different degrees of biophysical realism. We used a data-driven optimization algorithm to determine and classify the types of local dynamics that enable the reproduction of different observables derived from functional magnetic resonance recordings. The phase space analysis of the resulting equations revealed a predominance of stable spiral attractors, which optimized the similarity to the empirical data in terms of the synchronization, metastability, and functional connectivity dynamics. For stable limit cycles, departures from harmonic oscillations improved the fit in terms of functional connectivity dynamics. Eigenvalue analyses showed that proximity to a bifurcation improved the accuracy of the simulation for wakefulness, whereas deep sleep was associated with increased stability. Our results provide testable predictions that constrain the landscape of suitable biophysical models, while supporting noise-driven dynamics close to a bifurcation as a canonical mechanism underlying the complex fluctuations that characterize endogenous brain activity.

Reliability and subject specificity of personalized whole-brain dynamical models

Dynamical whole-brain models were developed to link structural (SC) and functional connectivity (FC) together into one framework. Nowadays, they are used to investigate the dynamical regimes of the brain and how these relate to behavioral, clinical and demographic traits. However, there is no comprehensive investigation on how reliable and subject specific the modeling results are given the variability of the empirical FC. In this study, we show that the parameters of these models can be fitted with a "poor" to "good" reliability depending on the exact implementation of the modeling paradigm. We find, as a general rule of thumb, that enhanced model personalization leads to increasingly reliable model parameters. In addition, we observe no clear effect of the model complexity evaluated by separately sampling results for linear, phase oscillator and neural mass network models. In fact, the most complex neural mass model often yields modeling results with "poor" reliability comparable to the simple linear model, but demonstrates an enhanced subject specificity of the model similarity maps. Subsequently, we show that the FC simulated by these models can outperform the empirical FC in terms of both reliability and subject specificity. For the structure-function relationship, simulated FC of individual subjects may be identified from the correlations with the empirical SC with an accuracy up to 70%, but not vice versa for non-linear models. We sample all our findings for 8 distinct brain parcellations and 6 modeling conditions and show that the parcellation-induced effect is much more pronounced for the modeling results than for the empirical data. In sum, this study provides an exploratory account on the reliability and subject specificity of dynamical whole-brain models and may be relevant for their further development and application. In particular, our findings suggest that the application of the dynamical whole-brain modeling should be tightly connected with an estimate of the reliability of the results.

Disentangling the critical signatures of neural activity

The critical brain hypothesis has emerged as an attractive framework to understand neuronal activity, but it is still widely debated. In this work, we analyze data from a multi-electrodes array in the rat's cortex and we find that power-law neuronal avalanches satisfying the crackling-noise relation coexist with spatial correlations that display typical features of critical systems. In order to shed a light on the underlying mechanisms at the origin of these signatures of criticality, we introduce a paradigmatic framework with a common stochastic modulation and pairwise linear interactions inferred from our data. We show that in such models power-law avalanches that satisfy the crackling-noise relation emerge as a consequence of the extrinsic modulation, whereas scale-free correlations are solely determined by internal interactions. Moreover, this disentangling is fully captured by the mutual information in the system. Finally, we show that analogous power-law avalanches are found in more realistic models of neural activity as well, suggesting that extrinsic modulation might be a broad mechanism for their generation.

A Riemannian Revisiting of Structure-Function Mapping Based on Eigenmodes

Understanding the link between brain structure and function may not only improve our knowledge of brain organization, but also lead to better quantification of pathology. To quantify this link, recent studies have attempted to predict the brain's functional connectivity from its structural connectivity. However, functional connectivity matrices live in the Riemannian manifold of the symmetric positive definite space and a specific attention must be paid to operate on this appropriate space. In this work we investigated the implications of using a distance based on an affine invariant Riemannian metric in the context of structure-function mapping. Specifically, we revisit previously proposed structure-function mappings based on eigendecomposition and test them on 100 healthy subjects from the Human Connectome Project using this adapted notion of distance. First, we show that using this Riemannian distance significantly alters the notion of similarity between subjects from a functional point of view. We also show that using this distance improves the correlation between the structural and functional similarity of different subjects. Finally, by using a distance appropriate to this manifold, we demonstrate the importance of mapping function from structure under the Riemannian manifold and show in particular that it is possible to outperform the group average and the so-called glass ceiling on the performance of mappings based on eigenmodes.

Predicting Functional Connectivity From Observed and Latent Structural Connectivity via Eigenvalue Mapping

Understanding how complex dynamic activity propagates over a static structural network is an overarching question in the field of neuroscience. Previous work has demonstrated that linear graph-theoretic models perform as well as non-linear neural simulations in predicting functional connectivity with the added benefits of low dimensionality and a closed-form solution which make them far less computationally expensive. Here we show a simple model relating the eigenvalues of the structural connectivity and functional networks using the Gamma function, producing a reliable prediction of functional connectivity with a single model parameter. We also investigate the impact of local activity diffusion and long-range interhemispheric connectivity on the structure-function model and show an improvement in functional connectivity prediction when accounting for such latent variables which are often excluded from traditional diffusion tensor imaging (DTI) methods.

High-resolution Virtual Brain Modeling Personalizes Deep Brain Stimulation for Treatment-Resistant Depression: Spatiotemporal Response Characteristics Following Stimulation of Neural Fiber Pathways

Over the past 15 years, deep brain stimulation (DBS) has been actively investigated as a groundbreaking therapy for patients with treatment-resistant depression (TRD); nevertheless, outcomes have varied from patient to patient, with an average response rate of ∼50%. The engagement of specific fiber tracts at the stimulation site has been hypothesized to be an important factor in determining outcomes, however, the resulting individual network effects at the whole-brain scale remain largely unknown. Here we provide a computational framework that can explore each individual's brain response characteristics elicited by selective stimulation of fiber tracts. We use a novel personalized in-silico approach, the Virtual Big Brain, which makes use of high-resolution virtual brain models at a mm-scale and explicitly reconstructs more than 100 000 fiber tracts for each individual. Each fiber tract is active and can be selectively stimulated. Simulation results demonstrate distinct stimulus-induced event-related potentials as a function of stimulation location, parametrized by the contact positions of the electrodes implanted in each patient, even though validation against empirical patient data reveals some limitations (i.e., the need for individual parameter adjustment, and differential accuracy across stimulation locations). This study provides evidence for the capacity of personalized high-resolution virtual brain models to investigate individual network effects in DBS for patients with TRD and opens up novel avenues in the personalized optimization of brain stimulation.

Network Controllability in Transmodal Cortex Predicts Positive Psychosis Spectrum Symptoms

BACKGROUND

The psychosis spectrum (PS) is associated with structural dysconnectivity concentrated in transmodal cortex. However, understanding of this pathophysiology has been limited by an overreliance on examining direct interregional connectivity. Using network control theory, we measured variation in both direct and indirect connectivity to a region to gain new insights into the pathophysiology of the PS.

METHODS

We used psychosis symptom data and structural connectivity in 1068 individuals from the Philadelphia Neurodevelopmental Cohort. Applying a network control theory metric called average controllability, we estimated each brain region's capacity to leverage its direct and indirect structural connections to control linear brain dynamics. Using nonlinear regression, we determined the accuracy with which average controllability could predict PS symptoms in out-of-sample testing. We also examined the predictive performance of regional strength, which indexes only direct connections to a region, as well as several graph-theoretic measures of centrality that index indirect connectivity. Finally, we assessed how the prediction performance for PS symptoms varied over the functional hierarchy spanning unimodal to transmodal cortex.

RESULTS

Average controllability outperformed all other connectivity features at predicting positive PS symptoms and was the only feature to yield above-chance predictive performance. Improved prediction for average controllability was concentrated in transmodal cortex, whereas prediction performance for strength was uniform across the cortex, suggesting that indexing indirect connections through average controllability is crucial in association cortex.

CONCLUSIONS

Examining interindividual variation in direct and indirect structural connections to transmodal cortex is crucial for accurate prediction of positive PS symptoms.

ADHD symptoms map onto noise-driven structure-function decoupling between hub and peripheral brain regions

Adults with childhood-onset attention-deficit hyperactivity disorder (ADHD) show altered whole-brain connectivity. However, the relationship between structural and functional brain abnormalities, the implications for the development of life-long debilitating symptoms, and the underlying mechanisms remain uncharted. We recruited a unique sample of 80 medication-naive adults with a clinical diagnosis of childhood-onset ADHD without psychiatric comorbidities, and 123 age-, sex-, and intelligence-matched healthy controls. Structural and functional connectivity matrices were derived from diffusion spectrum imaging and multi-echo resting-state functional MRI data. Hub, feeder, and local connections were defined using diffusion data. Individual-level measures of structural connectivity and structure-function coupling were used to contrast groups and link behavior to brain abnormalities. Computational modeling was used to test possible neural mechanisms underpinning observed group differences in the structure-function coupling. Structural connectivity did not significantly differ between groups but, relative to controls, ADHD showed a reduction in structure-function coupling in feeder connections linking hubs with peripheral regions. This abnormality involved connections linking fronto-parietal control systems with sensory networks. Crucially, lower structure-function coupling was associated with higher ADHD symptoms. Results from our computational model further suggest that the observed structure-function decoupling in ADHD is driven by heterogeneity in neural noise variability across brain regions. By highlighting a neural cause of a clinically meaningful breakdown in the structure-function relationship, our work provides novel information on the nature of chronic ADHD. The current results encourage future work assessing the genetic and neurobiological underpinnings of neural noise in ADHD, particularly in brain regions encompassed by fronto-parietal systems.

Virtual Connectomic Datasets in Alzheimer's Disease and Aging Using Whole-Brain Network Dynamics Modelling

Large neuroimaging datasets, including information about structural connectivity (SC) and functional connectivity (FC), play an increasingly important role in clinical research, where they guide the design of algorithms for automated stratification, diagnosis or prediction. A major obstacle is, however, the problem of missing features [e.g., lack of concurrent DTI SC and resting-state functional magnetic resonance imaging (rsfMRI) FC measurements for many of the subjects]. We propose here to address the missing connectivity features problem by introducing strategies based on computational whole-brain network modeling. Using two datasets, the Alzheimer’s Disease Neuroimaging Initiative (ADNI) dataset and a healthy aging dataset, for proof-of-concept, we demonstrate the feasibility of virtual data completion (i.e., inferring “virtual FC” from empirical SC or “virtual SC” from empirical FC), by using self-consistent simulations of linear and nonlinear brain network models. Furthermore, by performing machine learning classification (to separate age classes or control from patient subjects), we show that algorithms trained on virtual connectomes achieve discrimination performance comparable to when trained on actual empirical data; similarly, algorithms trained on virtual connectomes can be used to successfully classify novel empirical connectomes. Completion algorithms can be combined and reiterated to generate realistic surrogate connectivity matrices in arbitrarily large number, opening the way to the generation of virtual connectomic datasets with network connectivity information comparable to the one of the original data.

Brain Connectivity Studies on Structure-Function Relationships: A Short Survey with an Emphasis on Machine Learning

This short survey reviews the recent literature on the relationship between the brain structure and its functional dynamics. Imaging techniques such as diffusion tensor imaging (DTI) make it possible to reconstruct axonal fiber tracks and describe the structural connectivity (SC) between brain regions. By measuring fluctuations in neuronal activity, functional magnetic resonance imaging (fMRI) provides insights into the dynamics within this structural network. One key for a better understanding of brain mechanisms is to investigate how these fast dynamics emerge on a relatively stable structural backbone. So far, computational simulations and methods from graph theory have been mainly used for modeling this relationship. Machine learning techniques have already been established in neuroimaging for identifying functionally independent brain networks and classifying pathological brain states. This survey focuses on methods from machine learning, which contribute to our understanding of functional interactions between brain regions and their relation to the underlying anatomical substrate.

A graph neural network framework for causal inference in brain networks

A central question in neuroscience is how self-organizing dynamic interactions in the brain emerge on their relatively static structural backbone. Due to the complexity of spatial and temporal dependencies between different brain areas, fully comprehending the interplay between structure and function is still challenging and an area of intense research. In this paper we present a graph neural network (GNN) framework, to describe functional interactions based on the structural anatomical layout. A GNN allows us to process graph-structured spatio-temporal signals, providing a possibility to combine structural information derived from diffusion tensor imaging (DTI) with temporal neural activity profiles, like that observed in functional magnetic resonance imaging (fMRI). Moreover, dynamic interactions between different brain regions discovered by this data-driven approach can provide a multi-modal measure of causal connectivity strength. We assess the proposed model's accuracy by evaluating its capabilities to replicate empirically observed neural activation profiles, and compare the performance to those of a vector auto regression (VAR), like that typically used in Granger causality. We show that GNNs are able to capture long-term dependencies in data and also computationally scale up to the analysis of large-scale networks. Finally we confirm that features learned by a GNN can generalize across MRI scanner types and acquisition protocols, by demonstrating that the performance on small datasets can be improved by pre-training the GNN on data from an earlier study. We conclude that the proposed multi-modal GNN framework can provide a novel perspective on the structure-function relationship in the brain. Accordingly this approach appears to be promising for the characterization of the information flow in brain networks.

Robustness of connectome harmonics to local gray matter and long-range white matter connectivity changes

Recently, it has been proposed that the harmonic patterns emerging from the brain's structural connectivity underlie the resting state networks of the human brain. These harmonic patterns, termed connectome harmonics, are estimated as the Laplace eigenfunctions of the combined gray and white matters connectivity matrices and yield a connectome-specific extension of the well-known Fourier basis. However, it remains unclear how topological properties of the combined connectomes constrain the precise shape of the connectome harmonics and their relationships to the resting state networks. Here, we systematically study how alterations of the local and long-range connectivity matrices affect the spatial patterns of connectome harmonics. Specifically, the proportion of local gray matter homogeneous connectivity versus long-range white-matter heterogeneous connectivity is varied by means of weight-based matrix thresholding, distance-based matrix trimming, and several types of matrix randomizations. We demonstrate that the proportion of local gray matter connections plays a crucial role for the emergence of wide-spread, functionally meaningful, and originally published connectome harmonic patterns. This finding is robust for several different cortical surface templates, mesh resolutions, or widths of the local diffusion kernel. Finally, using the connectome harmonic framework, we also provide a proof-of-concept for how targeted structural changes such as the atrophy of inter-hemispheric callosal fibers and gray matter alterations may predict functional deficits associated with neurodegenerative conditions.

Whole-Brain Models to Explore Altered States of Consciousness from the Bottom Up

The scope of human consciousness includes states departing from what most of us experience as ordinary wakefulness. These altered states of consciousness constitute a prime opportunity to study how global changes in brain activity relate to different varieties of subjective experience. We consider the problem of explaining how global signatures of altered consciousness arise from the interplay between large-scale connectivity and local dynamical rules that can be traced to known properties of neural tissue. For this purpose, we advocate a research program aimed at bridging the gap between bottom-up generative models of whole-brain activity and the top-down signatures proposed by theories of consciousness. Throughout this paper, we define altered states of consciousness, discuss relevant signatures of consciousness observed in brain activity, and introduce whole-brain models to explore the biophysics of altered consciousness from the bottom-up. We discuss the potential of our proposal in view of the current state of the art, give specific examples of how this research agenda might play out, and emphasize how a systematic investigation of altered states of consciousness via bottom-up modeling may help us better understand the biophysical, informational, and dynamical underpinnings of consciousness.

A unified framework for multimodal structure-function mapping based on eigenmodes

Characterizing the connection between brain structure and brain function is essential for understanding how behaviour emerges from the underlying anatomy. A number of studies have shown that the network structure of the white matter shapes functional connectivity. Therefore, it should be possible to predict, at least partially, functional connectivity given the structural network. Many structure-function mappings have been proposed in the literature, including several direct mappings between the structural and functional connectivity matrices. However, the current literature is fragmented and does not provide a uniform treatment of current methods based on eigendecompositions. In particular, existing methods have never been compared to each other and their relationship explicitly derived in the context of brain structure-function mapping. In this work, we propose a unified computational framework that generalizes recently proposed structure-function mappings based on eigenmodes. Using this unified framework, we highlight the link between existing models and show how they can be obtained by specific choices of the parameters of our framework. By applying our framework to 50 subjects of the Human Connectome Project, we reproduce 6 recently published results, devise two new models and provide a direct comparison between all mappings. Finally, we show that a glass ceiling on the performance of mappings based on eigenmodes seems to be reached and conclude with possible approaches to break this performance limit.

The role of node dynamics in shaping emergent functional connectivity patterns in the brain

The contribution of structural connectivity to functional brain states remains poorly understood. We present a mathematical and computational study suited to assess the structure–function issue, treating a system of Jansen–Rit neural mass nodes with heterogeneous structural connections estimated from diffusion MRI data provided by the Human Connectome Project. Via direct simulations we determine the similarity of functional (inferred from correlated activity between nodes) and structural connectivity matrices under variation of the parameters controlling single-node dynamics, highlighting a nontrivial structure–function relationship in regimes that support limit cycle oscillations. To determine their relationship, we firstly calculate network instabilities giving rise to oscillations, and the so-called ‘false bifurcations’ (for which a significant qualitative change in the orbit is observed, without a change of stability) occurring beyond this onset. We highlight that functional connectivity (FC) is inherited robustly from structure when node dynamics are poised near a Hopf bifurcation, whilst near false bifurcations, and structure only weakly influences FC. Secondly, we develop a weakly coupled oscillator description to analyse oscillatory phase-locked states and, furthermore, show how the modular structure of FC matrices can be predicted via linear stability analysis. This study thereby emphasises the substantial role that local dynamics can have in shaping large-scale functional brain states.

Patterns of oscillation across the brain arise because of structural connections between brain regions. However, the type of oscillation at a site may also play a contributory role. We focus on an idealised model of a neural mass network, coupled using estimates of structural connections obtained via tractography on Human Connectome Project MRI data. Using a mixture of computational and mathematical techniques, we show that functional connectivity is inherited most strongly from structural connectivity when the network nodes are poised at a Hopf bifurcation. However, beyond the onset of this oscillatory instability a phase-locked network state can undergo a false bifurcation, and structural connectivity only weakly influences functional connectivity. This highlights the important effect that local dynamics can have on large-scale brain states.

Mapping functional brain networks from the structural connectome: Relating the series expansion and eigenmode approaches

Functional brain networks are shaped and constrained by the underlying structural network. However, functional networks are not merely a one-to-one reflection of the structural network. Several theories have been put forward to understand the relationship between structural and functional networks. However, it remains unclear how these theories can be unified. Two existing recent theories state that 1) functional networks can be explained by all possible walks in the structural network, which we will refer to as the series expansion approach, and 2) functional networks can be explained by a weighted combination of the eigenmodes of the structural network, the so-called eigenmode approach. To elucidate the unique or common explanatory power of these approaches to estimate functional networks from the structural network, we analysed the relationship between these two existing views. Using linear algebra, we first show that the eigenmode approach can be written in terms of the series expansion approach, i.e., walks on the structural network associated with different hop counts correspond to different weightings of the eigenvectors of this network. Second, we provide explicit expressions for the coefficients for both the eigenmode and series expansion approach. These theoretical results were verified by empirical data from Diffusion Tensor Imaging (DTI) and functional Magnetic Resonance Imaging (fMRI), demonstrating a strong correlation between the mappings based on both approaches. Third, we analytically and empirically demonstrate that the fit of the eigenmode approach to measured functional data is always at least as good as the fit of the series expansion approach, and that errors in the structural data lead to large errors of the estimated coefficients for the series expansion approach. Therefore, we argue that the eigenmode approach should be preferred over the series expansion approach. Results hold for eigenmodes of the weighted adjacency matrices as well as eigenmodes of the graph Laplacian. ​Taken together, these results provide an important step towards unification of existing theories regarding the structure-function relationships in brain networks.

Deriving pairwise transfer entropy from network structure and motifs

Transfer entropy (TE) is an established method for quantifying directed statistical dependencies in neuroimaging and complex systems datasets. The pairwise (or bivariate) TE from a source to a target node in a network does not depend solely on the local source-target link weight, but on the wider network structure that the link is embedded in. This relationship is studied using a discrete-time linearly coupled Gaussian model, which allows us to derive the TE for each link from the network topology. It is shown analytically that the dependence on the directed link weight is only a first approximation, valid for weak coupling. More generally, the TE increases with the in-degree of the source and decreases with the in-degree of the target, indicating an asymmetry of information transfer between hubs and low-degree nodes. In addition, the TE is directly proportional to weighted motif counts involving common parents or multiple walks from the source to the target, which are more abundant in networks with a high clustering coefficient than in random networks. Our findings also apply to Granger causality, which is equivalent to TE for Gaussian variables. Moreover, similar empirical results on random Boolean networks suggest that the dependence of the TE on the in-degree extends to nonlinear dynamics.

The role of node dynamics in shaping emergent functional connectivity patterns in the brain

The contribution of structural connectivity to functional brain states remains poorly understood. We present a mathematical and computational study suited to assess the structure-function issue, treating a system of Jansen-Rit neural mass nodes with heterogeneous structural connections estimated from diffusion MRI data provided by the Human Connectome Project. Via direct simulations we determine the similarity of functional (inferred from correlated activity between nodes) and structural connectivity matrices under variation of the parameters controlling single-node dynamics, highlighting a nontrivial structure-function relationship in regimes that support limit cycle oscillations. To determine their relationship, we firstly calculate network instabilities giving rise to oscillations, and the so-called 'false bifurcations' (for which a significant qualitative change in the orbit is observed, without a change of stability) occurring beyond this onset. We highlight that functional connectivity (FC) is inherited robustly from structure when node dynamics are poised near a Hopf bifurcation, whilst near false bifurcations, and structure only weakly influences FC. Secondly, we develop a weakly coupled oscillator description to analyse oscillatory phase-locked states and, furthermore, show how the modular structure of FC matrices can be predicted via linear stability analysis. This study thereby emphasises the substantial role that local dynamics can have in shaping large-scale functional brain states.

Individual structural features constrain the mouse functional connectome

Whole brain dynamics intuitively depend upon the internal wiring of the brain; but to which extent the individual structural connectome constrains the corresponding functional connectome is unknown, even though its importance is uncontested. After acquiring structural data from individual mice, we virtualized their brain networks and simulated in silico functional MRI data. Theoretical results were validated against empirical awake functional MRI data obtained from the same mice. We demonstrate that individual structural connectomes predict the functional organization of individual brains. Using a virtual mouse brain derived from the Allen Mouse Brain Connectivity Atlas, we further show that the dominant predictors of individual structure-function relations are the asymmetry and the weights of the structural links. Model predictions were validated experimentally using tracer injections, identifying which missing connections (not measurable with diffusion MRI) are important for whole brain dynamics in the mouse. Individual variations thus define a specific structural fingerprint with direct impact upon the functional organization of individual brains, a key feature for personalized medicine.

How do spatially distinct frequency specific MEG networks emerge from one underlying structural connectome? The role of the structural eigenmodes

Functional networks obtained from magnetoencephalography (MEG) from different frequency bands show distinct spatial patterns. It remains to be elucidated how distinct spatial patterns in MEG networks emerge given a single underlying structural network. Recent work has suggested that the eigenmodes of the structural network might serve as a basis set for functional network patterns in the case of functional MRI. Here, we take this notion further in the context of frequency band specific MEG networks. We show that a selected set of eigenmodes of the structural network can predict different frequency band specific networks in the resting state, ranging from delta (1-4 Hz) to the high gamma band (40-70 Hz). These predictions outperform predictions based from surrogate data, suggesting a genuine relationship between eigenmodes of the structural network and frequency specific MEG networks. We then show that the relevant set of eigenmodes can be excited in a network of neural mass models using linear stability analysis only by including delays. Excitation of an eigenmode in this context refers to a dynamic instability of a network steady state to a spatial pattern with a corresponding coherent temporal oscillation. Simulations verify the results from linear stability analysis and suggest that theta, alpha and beta band networks emerge very near to the bifurcation. The delta and gamma bands in the resting state emerges further away from the bifurcation. These results show for the first time how delayed interactions can excite the relevant set of eigenmodes that give rise to frequency specific functional connectivity patterns.

Homeostatic plasticity and emergence of functional networks in a whole-brain model at criticality

Understanding the relationship between large-scale structural and functional brain networks remains a crucial issue in modern neuroscience. Recently, there has been growing interest in investigating the role of homeostatic plasticity mechanisms, across different spatiotemporal scales, in regulating network activity and brain functioning against a wide range of environmental conditions and brain states (e.g., during learning, development, ageing, neurological diseases). In the present study, we investigate how the inclusion of homeostatic plasticity in a stochastic whole-brain model, implemented as a normalization of the incoming node's excitatory input, affects the macroscopic activity during rest and the formation of functional networks. Importantly, we address the structure-function relationship both at the group and individual-based levels. In this work, we show that normalization of the node's excitatory input improves the correspondence between simulated neural patterns of the model and various brain functional data. Indeed, we find that the best match is achieved when the model control parameter is in its critical value and that normalization minimizes both the variability of the critical points and neuronal activity patterns among subjects. Therefore, our results suggest that the inclusion of homeostatic principles lead to more realistic brain activity consistent with the hallmarks of criticality. Our theoretical framework open new perspectives in personalized brain modeling with potential applications to investigate the deviation from criticality due to structural lesions (e.g. stroke) or brain disorders.

Dynamics of large-scale electrophysiological networks: A technical review

For several years it has been argued that neural synchronisation is crucial for cognition. The idea that synchronised temporal patterns between different neural groups carries information above and beyond the isolated activity of these groups has inspired a shift in focus in the field of functional neuroimaging. Specifically, investigation into the activation elicited within certain regions by some stimulus or task has, in part, given way to analysis of patterns of co-activation or functional connectivity between distal regions. Recently, the functional connectivity community has been looking beyond the assumptions of stationarity that earlier work was based on, and has introduced methods to incorporate temporal dynamics into the analysis of connectivity. In particular, non-invasive electrophysiological data (magnetoencephalography/electroencephalography (MEG/EEG)), which provides direct measurement of whole-brain activity and rich temporal information, offers an exceptional window into such (potentially fast) brain dynamics. In this review, we discuss challenges, solutions, and a collection of analysis tools that have been developed in recent years to facilitate the investigation of dynamic functional connectivity using these imaging modalities. Further, we discuss the applications of these approaches in the study of cognition and neuropsychiatric disorders. Finally, we review some existing developments that, by using realistic computational models, pursue a deeper understanding of the underlying causes of non-stationary connectivity.

Warnings and caveats in brain controllability

A recent article by Gu et al. (Nat. Commun. 6, 2015) proposed to characterize brain networks, quantified using anatomical diffusion imaging, in terms of their "controllability", drawing on concepts and methods of control theory. They reported that brain activity is controllable from a single node, and that the topology of brain networks provides an explanation for the types of control roles that different regions play in the brain. In this work, we first briefly review the framework of control theory applied to complex networks. We then show contrasting results on brain controllability through the analysis of five different datasets and numerical simulations. We find that brain networks are not controllable (in a statistical significant way) by one single region. Additionally, we show that random null models, with no biological resemblance to brain network architecture, produce the same type of relationship observed by Gu et al. between the average/modal controllability and weighted degree. Finally, we find that resting state networks defined with fMRI cannot be attributed specific control roles. In summary, our study highlights some warning and caveats in the brain controllability framework.

Toward a theory of coactivation patterns in excitable neural networks

The relationship between the structural connectivity (SC) and functional connectivity (FC) of neural systems is of central importance in brain network science. It is an open question, however, how the SC-FC relationship depends on specific topological features of brain networks or the models used for describing neural dynamics. Using a basic but general model of discrete excitable units that follow a susceptible—excited—refractory activity cycle (SER model), we here analyze how the network activity patterns underlying functional connectivity are shaped by the characteristic topological features of the network. We develop an analytical framework for describing the contribution of essential topological elements, such as common inputs and pacemakers, to the coactivation of nodes, and demonstrate the validity of the approach by comparison of the analytical predictions with numerical simulations of various exemplar networks. The present analytic framework may serve as an initial step for the mechanistic understanding of the contributions of brain network topology to brain dynamics.

Functional connectivity, as reflected in the statistical dependencies of distributed activity, is widely used to probe the organization of complex systems such as the brain. While this measure has been helpful for characterizing brain states and highlighting alterations of brain dynamics in various diseases, the mechanisms underlying the generation of FC patterns remain poorly understood. One prominent factor shaping FC is the underlying neural network structure. Using a minimalist model of excitation, we investigate how the topology of excitable neural networks contributes to FC. Specifically, we show that FC can be analytically predicted from the way in which the nodes are embedded in the network and how they are related to basic self-organizing units of excitable dynamics, particularly, short pacemaker cycles. These insights are a step towards a mechanistic understanding of the activation patterns of complex neural networks.

Multiple Kernel Learning Model for Relating Structural and Functional Connectivity in the Brain

A challenging problem in cognitive neuroscience is to relate the structural connectivity (SC) to the functional connectivity (FC) to better understand how large-scale network dynamics underlying human cognition emerges from the relatively fixed SC architecture. Recent modeling attempts point to the possibility of a single diffusion kernel giving a good estimate of the FC. We highlight the shortcomings of the single-diffusion-kernel model (SDK) and propose a multi-scale diffusion scheme. Our multi-scale model is formulated as a reaction-diffusion system giving rise to spatio-temporal patterns on a fixed topology. We hypothesize the presence of inter-regional co-activations (latent parameters) that combine diffusion kernels at multiple scales to characterize how FC could arise from SC. We formulated a multiple kernel learning (MKL) scheme to estimate the latent parameters from training data. Our model is analytically tractable and complex enough to capture the details of the underlying biological phenomena. The parameters learned by the MKL model lead to highly accurate predictions of subject-specific FCs from test datasets at a rate of 71%, surpassing the performance of the existing linear and non-linear models. We provide an example of how these latent parameters could be used to characterize age-specific reorganization in the brain structure and function.

How structure sculpts function: Unveiling the contribution of anatomical connectivity to the brain's spontaneous correlation structure

Intrinsic brain activity is characterized by highly organized co-activations between different regions, forming clustered spatial patterns referred to as resting-state networks. The observed co-activation patterns are sustained by the intricate fabric of millions of interconnected neurons constituting the brain's wiring diagram. However, as for other real networks, the relationship between the connectional structure and the emergent collective dynamics still evades complete understanding. Here, we show that it is possible to estimate the expected pair-wise correlations that a network tends to generate thanks to the underlying path structure. We start from the assumption that in order for two nodes to exhibit correlated activity, they must be exposed to similar input patterns from the entire network. We then acknowledge that information rarely spreads only along a unique route but rather travels along all possible paths. In real networks, the strength of local perturbations tends to decay as they propagate away from the sources, leading to a progressive attenuation of the original information content and, thus, of their influence. Accordingly, we define a novel graph measure, topological similarity, which quantifies the propensity of two nodes to dynamically correlate as a function of the resemblance of the overall influences they are expected to receive due to the underlying structure of the network. Applied to the human brain, we find that the similarity of whole-network inputs, estimated from the topology of the anatomical connectome, plays an important role in sculpting the backbone pattern of time-average correlations observed at rest.

Mapping complementary features of cross-species structural connectivity to construct realistic "Virtual Brains"

Modern systems neuroscience increasingly leans on large-scale multi-lab neuroinformatics initiatives to provide necessary capacity for biologically realistic modeling of primate whole-brain activity. Here, we present a framework to assemble primate brain's biologically plausible anatomical backbone for such modeling initiatives. In this framework, structural connectivity is determined by adding complementary information from invasive macaque axonal tract tracing and non-invasive human diffusion tensor imaging. Both modalities are combined by means of available interspecies registration tools and a newly developed Bayesian probabilistic modeling approach to extract common connectivity evidence. We demonstrate how this novel framework is embedded in the whole-brain simulation platform called The Virtual Brain (TVB). Hum Brain Mapp 38:2080-2093, 2017. © 2017 Wiley Periodicals, Inc.

Cognitive decline in metabolic syndrome is linked to microstructural white matter abnormalities

Subjects with metabolic syndrome (MetS) often show worse cognitive performance compared with the healthy population. We investigated whether microstructural white matter abnormalities are associated with cognitive performance in adults with MetS using diffusion tensor MR imaging. A total of 32 subjects with MetS (age 64.8 ± 7.8, 56.25 % female) and 23 age-, gender-, and education-matched healthy controls completed a battery of neuropsychological tests and diffusion tensor imaging (DTI) at 3-T MRI. Brain global and regional volumes, white matter fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (LD) were calculated. The least-square models adjusted for age, sex, HbA1c, hypertension, body mass index, hyperlipidemia, and white matter hyperintensities were used to evaluate the relationship between cognitive function and DTI. The MetS group had worse performance in verbal fluency (VF) and learning and memory function (total VF: T score (p = 0.01), VF: animals T score (p = 0.0001), Hopkins Verbal Learning Test (HVLT): Total recall T score (p = 0.0001), and HVLT: delayed recall T score (p = 0.002), as compared with controls. In the MetS group, abnormalities in diffusivity measures were associated with worse cognitive performance [VF: animals T score and left post-central gyrus-LD (p = 0.0007, r _adj 0.4), R angular gyrus-RD (p = 0.0008, r _adj 0.3), L supra-marginal gyrus-RD (p = 0.009, r _adj 0.2) after adjusting for age, sex, HbA1c, 24 h mean BP, presence of hyperlipidemia, and global white matter hyperintensities]. Microstructural white matter abnormalities in the MetS group might be the underlying mechanisms of worse verbal learning and memory performance.