Analysis of Dynamic Brain Imaging Data

Modern imaging techniques for probing brain function, including functional Magnetic Resonance Imaging, intrinsic and extrinsic contrast optical imaging, and magnetoencephalography, generate large data sets with complex content. In this paper we develop appropriate techniques of analysis and visualization of such imaging data, in order to separate the signal from the noise, as well as to characterize the signal. The techniques developed fall into the general category of multivariate time series analysis, and in particular we extensively use the multitaper framework of spectral analysis. We develop specific protocols for the analysis of fMRI, optical imaging and MEG data, and illustrate the techniques by applications to real data sets generated by these imaging modalities. In general, the analysis protocols involve two distinct stages: `noise' characterization and suppression, and `signal' characterization and visualization. An important general conclusion of our study is the utility of a frequency-based representation, with short, moving analysis windows to account for non-stationarity in the data. Of particular note are (a) the development of a decomposition technique (`space-frequency singular value decomposition') that is shown to be a useful means of characterizing the image data, and (b) the development of an algorithm, based on multitaper methods, for the removal of approximately periodic physiological artifacts arising from cardiac and respiratory sources.Comment: 40 pages; 26 figures with subparts including 3 figures as .gif files. Originally submitted to the neuro-sys archive which was never publicly announced (was 9804003

Gearing up for action: attentive tracking dynamically tunes sensory and motor oscillations in the alpha and beta band

Allocation of attention during goal-directed behavior entails simultaneous processing of relevant and attenuation of irrelevant information. How the brain delegates such processes when confronted with dynamic (biological motion) stimuli and harnesses relevant sensory information for sculpting prospective responses remains unclear. We analyzed neuromagnetic signals that were recorded while participants attentively tracked an actor’s pointing movement that ended at the location where subsequently the response-cue indicated the required response. We found the observers’ spatial allocation of attention to be dynamically reflected in lateralized parieto-occipital alpha (8-12Hz) activity and to have a lasting influence on motor preparation. Specifically, beta (16-25Hz) power modulation reflected observers’ tendency to selectively prepare for a spatially compatible response even before knowing the required one. We discuss the observed frequency-specific and temporally evolving neural activity within a framework of integrated visuomotor processing and point towards possible implications about the mechanisms involved in action observation

MEG sensor and source measures of visually induced gamma-band oscillations are highly reliable

High frequency brain oscillations are associated with numerous cognitive and behavioral processes. Non-invasive measurements using electro-/magnetoencephalography (EEG/MEG) have revealed that high frequency neural signals are heritable and manifest changes with age as well as in neuropsychiatric illnesses. Despite the extensive use of EEG/MEG-measured neural oscillations in basic and clinical research, studies demonstrating test–retest reliability of power and frequency measures of neural signals remain scarce. Here, we evaluated the test–retest reliability of visually induced gamma (30–100 Hz) oscillations derived from sensor and source signals acquired over two MEG sessions. The study required participants (N = 13) to detect the randomly occurring stimulus acceleration while viewing a moving concentric grating. Sensor and source MEG measures of gamma-band activity yielded comparably strong reliability (average intraclass correlation, ICC = 0.861). Peak stimulus-induced gamma frequency (53–72 Hz) yielded the highest measures of stability (ICCsensor = 0.940; ICCsource = 0.966) followed by spectral signal change (ICCsensor = 0.890; ICCsource = 0.893) and peak frequency bandwidth (ICCsensor = 0.856; ICCsource = 0.622). Furthermore, source-reconstruction significantly improved signal-to-noise for spectral amplitude of gamma activity compared to sensor estimates. Our assessments highlight that both sensor and source derived estimates of visually induced gamma-band oscillations from MEG signals are characterized by high test–retest reliability, with source derived oscillatory measures conferring an improvement in the stability of peak-frequency estimates. Importantly, our finding of high test–retest reliability supports the feasibility of pharma-MEG studies and longitudinal aging or clinical studies

Frequency and Phase Synchronization in Neuromagnetic Cortical Responses to Flickering-Color Stimuli

In our earlier study dealing with the analysis of neuromagnetic responses (magnetoencephalograms - MEG) to flickering-color stimuli for a group of control human subjects (9 volunteers) and a patient with photosensitive epilepsy (a 12-year old girl), it was shown that Flicker-Noise Spectroscopy (FNS) was able to identify specific differences in the responses of each organism. The high specificity of individual MEG responses manifested itself in the values of FNS parameters for both chaotic and resonant components of the original signal. The present study applies the FNS cross-correlation function to the analysis of correlations between the MEG responses simultaneously measured at spatially separated points of the human cortex processing the red-blue flickering color stimulus. It is shown that the cross-correlations for control (healthy) subjects are characterized by frequency and phase synchronization at different points of the cortex, with the dynamics of neuromagnetic responses being determined by the low-frequency processes that correspond to normal physiological rhythms. But for the patient, the frequency and phase synchronization breaks down, which is associated with the suppression of cortical regulatory functions when the flickering-color stimulus is applied, and higher frequencies start playing the dominating role. This suggests that the disruption of correlations in the MEG responses is the indicator of pathological changes leading to photosensitive epilepsy, which can be used for developing a method of diagnosing the disease based on the analysis with the FNS cross-correlation function.Comment: 21 pages, 14 figures; submitted to "Laser Physics", 2010, 2

Perceptually relevant speech tracking in auditory and motor cortex reflects distinct linguistic features

During online speech processing, our brain tracks the acoustic fluctuations in speech at different timescales. Previous research has focused on generic timescales (for example, delta or theta bands) that are assumed to map onto linguistic features such as prosody or syllables. However, given the high intersubject variability in speaking patterns, such a generic association between the timescales of brain activity and speech properties can be ambiguous. Here, we analyse speech tracking in source-localised magnetoencephalographic data by directly focusing on timescales extracted from statistical regularities in our speech material. This revealed widespread significant tracking at the timescales of phrases (0.6–1.3 Hz), words (1.8–3 Hz), syllables (2.8–4.8 Hz), and phonemes (8–12.4 Hz). Importantly, when examining its perceptual relevance, we found stronger tracking for correctly comprehended trials in the left premotor (PM) cortex at the phrasal scale as well as in left middle temporal cortex at the word scale. Control analyses using generic bands confirmed that these effects were specific to the speech regularities in our stimuli. Furthermore, we found that the phase at the phrasal timescale coupled to power at beta frequency (13–30 Hz) in motor areas. This cross-frequency coupling presumably reflects top-down temporal prediction in ongoing speech perception. Together, our results reveal specific functional and perceptually relevant roles of distinct tracking and cross-frequency processes along the auditory–motor pathway

Efficient transfer entropy analysis of non-stationary neural time series

Information theory allows us to investigate information processing in neural systems in terms of information transfer, storage and modification. Especially the measure of information transfer, transfer entropy, has seen a dramatic surge of interest in neuroscience. Estimating transfer entropy from two processes requires the observation of multiple realizations of these processes to estimate associated probability density functions. To obtain these observations, available estimators assume stationarity of processes to allow pooling of observations over time. This assumption however, is a major obstacle to the application of these estimators in neuroscience as observed processes are often non-stationary. As a solution, Gomez-Herrero and colleagues theoretically showed that the stationarity assumption may be avoided by estimating transfer entropy from an ensemble of realizations. Such an ensemble is often readily available in neuroscience experiments in the form of experimental trials. Thus, in this work we combine the ensemble method with a recently proposed transfer entropy estimator to make transfer entropy estimation applicable to non-stationary time series. We present an efficient implementation of the approach that deals with the increased computational demand of the ensemble method's practical application. In particular, we use a massively parallel implementation for a graphics processing unit to handle the computationally most heavy aspects of the ensemble method. We test the performance and robustness of our implementation on data from simulated stochastic processes and demonstrate the method's applicability to magnetoencephalographic data. While we mainly evaluate the proposed method for neuroscientific data, we expect it to be applicable in a variety of fields that are concerned with the analysis of information transfer in complex biological, social, and artificial systems.Comment: 27 pages, 7 figures, submitted to PLOS ON