Connectome spectrum electromagnetic tomography: A method to reconstruct electrical brain source networks at high-spatial resolution.

Connectome spectrum electromagnetic tomography (CSET) combines diffusion MRI-derived structural connectivity data with well-established graph signal processing tools to solve the M/EEG inverse problem. Using simulated EEG signals from fMRI responses, and two EEG datasets on visual-evoked potentials, we provide evidence supporting that (i) CSET captures realistic neurophysiological patterns with better accuracy than state-of-the-art methods, (ii) CSET can reconstruct brain responses more accurately and with more robustness to intrinsic noise in the EEG signal. These results demonstrate that CSET offers high spatio-temporal accuracy, enabling neuroscientists to extend their research beyond the current limitations of low sampling frequency in functional MRI and the poor spatial resolution of M/EEG.

Computational Models in Electroencephalography.

Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses in silico and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by "computational model" is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses in silico. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.

The connectome spectrum as a canonical basis for a sparse representation of fast brain activity.

The functional organization of neural processes is constrained by the brain's intrinsic structural connectivity, i.e., the connectome. Here, we explore how structural connectivity can improve the representation of brain activity signals and their dynamics. Using a multi-modal imaging dataset (electroencephalography, structural MRI, and diffusion MRI), we represent electrical brain activity at the cortical surface as a time-varying composition of harmonic modes of structural connectivity. These harmonic modes are known as connectome harmonics. Here we describe brain activity signal as a time-varying combination of connectome harmonics. We term this description as the connectome spectrum of the signal. We found that: first, the brain activity signal is represented more compactly by the connectome spectrum than by the traditional area-based representation; second, the connectome spectrum characterizes fast brain dynamics in terms of signal broadcasting profile, revealing different temporal regimes of integration and segregation that are consistent across participants. And last, the connectome spectrum characterizes fast brain dynamics with fewer degrees of freedom than area-based signal representations. Specifically, we show that a smaller number of dimensions capture the differences between low-level and high-level visual processing in the connectome spectrum. Also, we demonstrate that connectome harmonics capture more sensitively the topological properties of brain activity. In summary, this work provides statistical, functional, and topological evidence indicating that the description of brain activity in terms of structural connectivity fosters a more comprehensive understanding of large-scale dynamic neural functioning.

Functional harmonics reveal multi-dimensional basis functions underlying cortical organization.

The human brain consists of specialized areas that flexibly interact to form a multitude of functional networks. Complementary to this notion of modular organization, brain function has been shown to vary along a smooth continuum across the whole cortex. We demonstrate a mathematical framework that accounts for both of these perspectives: harmonic modes. We calculate the harmonic modes of the brain's functional connectivity graph, called "functional harmonics," revealing a multi-dimensional, frequency-ordered set of basis functions. Functional harmonics link characteristics of cortical organization across several spatial scales, capturing aspects of intra-areal organizational features (retinotopy, somatotopy), delineating brain areas, and explaining macroscopic functional networks as well as global cortical gradients. Furthermore, we show how the activity patterns elicited by seven different tasks are reconstructed from a very small subset of functional harmonics. Our results suggest that the principle of harmonicity, ubiquitous in nature, also underlies functional cortical organization in the human brain.

Using structural connectivity to augment community structure in EEG functional connectivity.

Recently, EEG recording techniques and source analysis have improved, making it feasible to tap into fast network dynamics. Yet, analyzing whole-cortex EEG signals in source space is not standard, partly because EEG suffers from volume conduction: Functional connectivity (FC) reflecting genuine functional relationships is impossible to disentangle from spurious FC introduced by volume conduction. Here, we investigate the relationship between white matter structural connectivity (SC) and large-scale network structure encoded in EEG-FC. We start by confirming that FC (power envelope correlations) is predicted by SC beyond the impact of Euclidean distance, in line with the assumption that SC mediates genuine FC. We then use information from white matter structural connectivity in order to smooth the EEG signal in the space spanned by graphs derived from SC. Thereby, FC between nearby, structurally connected brain regions increases while FC between nonconnected regions remains unchanged, resulting in an increase in genuine, SC-mediated FC. We analyze the induced changes in FC, assessing the resemblance between EEG-FC and volume-conduction- free fMRI-FC, and find that smoothing increases resemblance in terms of overall correlation and community structure. This result suggests that our method boosts genuine FC, an outcome that is of interest for many EEG network neuroscience questions.

Connectome spectral analysis to track EEG task dynamics on a subsecond scale.

We present an approach for tracking fast spatiotemporal cortical dynamics in which we combine white matter connectivity data with source-projected electroencephalographic (EEG) data. We employ the mathematical framework of graph signal processing in order to derive the Fourier modes of the brain structural connectivity graph, or "network harmonics". These network harmonics are naturally ordered by smoothness. Smoothness in this context can be understood as the amount of variation along the cortex, leading to a multi-scale representation of brain connectivity. We demonstrate that network harmonics provide a sparse representation of the EEG signal, where, at certain times, the smoothest 15 network harmonics capture 90% of the signal power. This suggests that network harmonics are functionally meaningful, which we demonstrate by using them as a basis for the functional EEG data recorded from a face detection task. There, only 13 network harmonics are sufficient to track the large-scale cortical activity during the processing of the stimuli with a 50 â€‹ms resolution, reproducing well-known activity in the fusiform face area as well as revealing co-activation patterns in somatosensory/motor and frontal cortices that an unconstrained ROI-by-ROI analysis fails to capture. The proposed approach is simple and fast, provides a means of integration of multimodal datasets, and is tied to a theoretical framework in mathematics and physics. Thus, network harmonics point towards promising research directions both theoretically - for example in exploring the relationship between structure and function in the brain - and practically - for example for network tracking in different tasks and groups of individuals, such as patients.

Neural mechanisms of vibrotactile categorization.

The grouping of sensory stimuli into categories is fundamental to cognition. Previous research in the visual and auditory systems supports a two-stage processing hierarchy that underlies perceptual categorization: (a) a "bottom-up" perceptual stage in sensory cortices where neurons show selectivity for stimulus features and (b) a "top-down" second stage in higher level cortical areas that categorizes the stimulus-selective input from the first stage. In order to test the hypothesis that the two-stage model applies to the somatosensory system, 14 human participants were trained to categorize vibrotactile stimuli presented to their right forearm. Then, during an fMRI scan, participants actively categorized the stimuli. Representational similarity analysis revealed stimulus selectivity in areas including the left precentral and postcentral gyri, the supramarginal gyrus, and the posterior middle temporal gyrus. Crucially, we identified a single category-selective region in the left ventral precentral gyrus. Furthermore, an estimation of directed functional connectivity delivered evidence for robust top-down connectivity from the second to first stage. These results support the validity of the two-stage model of perceptual categorization for the somatosensory system, suggesting common computational principles and a unified theory of perceptual categorization across the visual, auditory, and somatosensory systems.

Stereotypical modulations in dynamic functional connectivity explained by changes in BOLD variance.

Spontaneous activity measured in human subject under the absence of any task exhibits complex patterns of correlation that largely correspond to large-scale functional topographies obtained with a wide variety of cognitive and perceptual tasks. These "resting state networks" (RSNs) fluctuate over time, forming and dissolving on the scale of seconds to minutes. While these fluctuations, most prominently those of the default mode network, have been linked to cognitive function, it remains unclear whether they result from random noise or whether they index a nonstationary process which could be described as state switching. In this study, we use a sliding windows-approach to relate temporal dynamics of RSNs to global modulations in correlation and BOLD variance. We compare empirical data, phase-randomized surrogate data, and data simulated with a stationary model. We find that RSN time courses exhibit a large amount of coactivation in all three cases, and that the modulations in their activity are closely linked to global dynamics of the underlying BOLD signal. We find that many properties of the observed fluctuations in FC and BOLD, including their ranges and their correlations amongst each other, are explained by fluctuations around the average FC structure. However, we also report some interesting characteristics that clearly support nonstationary features in the data. In particular, we find that the brain spends more time in the troughs of modulations than can be expected from stationary dynamics.

Resting state networks in empirical and simulated dynamic functional connectivity.

It is well-established that patterns of functional connectivity (FC) - measures of correlated activity between pairs of voxels or regions observed in the human brain using neuroimaging - are robustly expressed in spontaneous activity during rest. These patterns are not static, but exhibit complex spatio-temporal dynamics. Over the last years, a multitude of methods have been proposed to reveal these dynamics on the level of the whole brain. One finding is that the brain transitions through different FC configurations over time, and substantial effort has been put into characterizing these configurations. However, the dynamics governing these transitions are more elusive, specifically, the contribution of stationary vs. non-stationary dynamics is an active field of inquiry. In this study, we use a whole-brain approach, considering FC dynamics between 66 ROIs covering the entire cortex. We combine an innovative dimensionality reduction technique, tensor decomposition, with a mean field model which possesses stationary dynamics. It has been shown to explain resting state FC averaged over time and multiple subjects, however, this average FC summarizes the spatial distribution of correlations while hiding their temporal dynamics. First, we apply tensor decomposition to resting state scans from 24 healthy controls in order to characterize spatio-temporal dynamics present in the data. We simultaneously utilize temporal and spatial information by creating tensors that are subsequently decomposed into sets of brain regions ("communities") that share similar temporal dynamics, and their associated time courses. The tensors contain pairwise FC computed inside of overlapping sliding windows. Communities are discovered by clustering features pooled from all subjects, thereby ensuring that they generalize. We find that, on the group level, the data give rise to four distinct communities that resemble known resting state networks (RSNs): default mode network, visual network, control networks, and somatomotor network. Second, we simulate data with our stationary mean field model whose nodes are connected according to results from DTI and fiber tracking. In this model, all spatio-temporal structure is due to noisy fluctuations around the average FC. We analyze the simulated data in the same way as the empirical data in order to determine whether stationary dynamics can explain the emergence of distinct FC patterns (RSNs) which have their own time courses. We find that this is the case for all four networks using the spatio-temporal information revealed by tensor decomposition if nodes in the simulation are connected according to model-based effective connectivity. Furthermore, we find that these results require only a small part of the FC values, namely the highest values that occur across time and ROI pair. Our findings show that stationary dynamics can account for the emergence of RSNs. We provide an innovative method that does not make strong assumptions about the underlying data and is generally applicable to resting state or task data from different subject populations.