Latency analysis of resting-state BOLD-fMRI reveals traveling waves in visual cortex linking task-positive and task-negative networks.

Due to the low temporal resolution of BOLD-fMRI, imaging studies on human brain function have almost exclusively focused on instantaneous correlations within the data. Developments in hardware and acquisition protocols, however, are offering data with higher sampling rates that allow investigating the latency structure of BOLD-fMRI data. In this study we describe a method for analyzing the latency structure within BOLD-fMRI data and apply it to resting-state data of 94 participants from the Human Connectome Project. The method shows that task-positive and task-negative networks are integrated through traveling BOLD waves within early visual cortex. The waves are initiated at the periphery of the visual field and propagate towards the fovea. This observation suggests a mechanism for the functional integration of task-positive and task-negative networks, argues for an eccentricity-based view on visual information processing, and contributes to the emerging view that resting-state BOLD-fMRI fluctuations are superpositions of inherently spatiotemporal patterns.

Source-reconstruction of the sensorimotor network from resting-state macaque electrocorticography.

The discovery of hemodynamic (BOLD-fMRI) resting-state networks (RSNs) has brought about a fundamental shift in our thinking about the role of intrinsic brain activity. The electrophysiological underpinnings of RSNs remain largely elusive and it has been shown only recently that electric cortical rhythms are organized into the same RSNs as hemodynamic signals. Most electrophysiological studies into RSNs use magnetoencephalography (MEG) or scalp electroencephalography (EEG), which limits the spatial resolution with which electrophysiological RSNs can be observed. Due to their close proximity to the cortical surface, electrocorticographic (ECoG) recordings can potentially provide a more detailed picture of the functional organization of resting-state cortical rhythms, albeit at the expense of spatial coverage. In this study we propose using source-space spatial independent component analysis (spatial ICA) for identifying generators of resting-state cortical rhythms as recorded with ECoG and for reconstructing their functional connectivity. Network structure is assessed by two kinds of connectivity measures: instantaneous correlations between band-limited amplitude envelopes and oscillatory phase-locking. By simulating rhythmic cortical generators, we find that the reconstruction of oscillatory phase-locking is more challenging than that of amplitude correlations, particularly for low signal-to-noise levels. Specifically, phase-lags can both be over- and underestimated, which troubles the interpretation of lag-based connectivity measures. We illustrate the methodology on somatosensory beta rhythms recorded from a macaque monkey using ECoG. The methodology decomposes the resting-state sensorimotor network into three cortical generators, distributed across primary somatosensory and primary and higher-order motor areas. The generators display significant and reproducible amplitude correlations and phase-locking values with non-zero lags. Our findings illustrate the level of spatial detail attainable with source-projected ECoG and motivates wider use of the methodology for studying resting-state as well as event-related cortical dynamics in macaque and human.

Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI?

During the last several years, the focus of research on resting-state functional magnetic resonance imaging (fMRI) has shifted from the analysis of functional connectivity averaged over the duration of scanning sessions to the analysis of changes of functional connectivity within sessions. Although several studies have reported the presence of dynamic functional connectivity (dFC), statistical assessment of the results is not always carried out in a sound way and, in some studies, is even omitted. In this study, we explain why appropriate statistical tests are needed to detect dFC, we describe how they can be carried out and how to assess the performance of dFC measures, and we illustrate the methodology using spontaneous blood-oxygen level-dependent (BOLD) fMRI recordings of macaque monkeys under general anesthesia and in human subjects under resting-state conditions. We mainly focus on sliding-window correlations since these are most widely used in assessing dFC, but also consider a recently proposed non-linear measure. The simulations and methodology, however, are general and can be applied to any measure. The results are twofold. First, through simulations, we show that in typical resting-state sessions of 10 min, it is almost impossible to detect dFC using sliding-window correlations. This prediction is validated by both the macaque and the human data: in none of the individual recording sessions was evidence for dFC found. Second, detection power can be considerably increased by session- or subject-averaging of the measures. In doing so, we found that most of the functional connections are in fact dynamic. With this study, we hope to raise awareness of the statistical pitfalls in the assessment of dFC and how they can be avoided by using appropriate statistical methods.

Role of white-matter pathways in coordinating alpha oscillations in resting visual cortex.

In the absence of cognitive tasks and external stimuli, strong rhythmic fluctuations with a frequency ≈ 10 Hz emerge from posterior regions of human neocortex. These posterior α-oscillations can be recorded throughout the visual cortex and are particularly strong in the calcarine sulcus, where the primary visual cortex is located. The mechanisms and anatomical pathways through which local \alpha-oscillations are coordinated however, are not fully understood. In this study, we used a combination of magnetoencephalography (MEG), diffusion tensor imaging (DTI), and biophysical modeling to assess the role of white-matter pathways in coordinating cortical α-oscillations. Our findings suggest that primary visual cortex plays a special role in coordinating α-oscillations in higher-order visual regions. Specifically, the amplitudes of α-sources throughout visual cortex could be explained by propagation of α-oscillations from primary visual cortex through white-matter pathways. In particular, α-amplitudes within visual cortex correlated with both the anatomical and functional connection strengths to primary visual cortex. These findings reinforce the notion of posterior α-oscillations as intrinsic oscillations of the visual system. We speculate that they might reflect a default-mode of the visual system during which higher-order visual regions are rhythmically primed for expected visual stimuli by α-oscillations in primary visual cortex.

Unbiased estimation of Langevin dynamics from time series with application to hippocampal field potentials in vitro.

The last decade showed an increased interest in Langevin equations for modeling time series recorded from complex dynamical systems. These equations allow to discriminate between deterministic (drift) and stochastic (diffusion) components of the recorded time series. In practice, the estimation of drift and diffusion is often based on approximations of the models' dynamics that are valid only for high sampling frequencies. Also, model assessment is not or only indirectly performed, potentially leading to false claims. In this study we compare the performance of an asymptotically unbiased estimation method with a generally used approximate method, demonstrating the necessity of using (asymptotically) unbiased estimators. Furthermore, we describe how confidence intervals for the unknown parameters can be constructed and how model assessment can be carried out. We apply the methodology to local field potentials recorded in vitro from mouse hippocampus from eight genetically different strains. The recorded field potentials turn out to be well described by linearly damped Langevin equations with parabolic diffusion. The modeling enables a dynamical interpretation of the spectral power of the field potentials. It reveals that observed spectral power differences in the field potentials across hippocampal regions are associated with differences in the deterministic component of the system, and it reveals transiently active current dipoles, which are not detectable by conventional methods. Also, all estimated parameters have significant heritabilities, which suggests that the Langevin equations capture biological relevant aspects of electrical hippocampal activity.

Dynamics underlying spontaneous human alpha oscillations: a data-driven approach.

Although the cognitive and clinical correlates of spontaneous human alpha oscillations as recorded with electroencephalography (EEG) or magnetoencephalography (MEG) are well documented, the dynamics underlying these oscillations is still a matter of debate. This study proposes a data-driven method to reveal the dynamics of these oscillations. It demonstrates that spontaneous human alpha oscillations as recorded with MEG can be viewed as noise-perturbed damped harmonic oscillations. This provides evidence for the hypothesis that these oscillations reflect filtered noise and hence do not possess limit-cycle dynamics. To illustrate the use of the model, we apply it to two data-sets in which a decrease in alpha power can be observed across conditions. The associated differences in the estimated model parameters show that observed decreases in alpha power are associated with different kinds of changes in the dynamics. Thus, the model parameters are useful dynamical biomarkers for spontaneous human alpha oscillations.