Large-amplitude, spatially correlated fluctuations in BOLD fMRI signals during extended rest and early sleep stages

A number of recent studies of human brain activity using blood-oxygen-level-dependent (BOLD) fMRI and EEG have reported the presence of spatiotemporal patterns of correlated activity in the absence of external stimuli. Although these patterns have been hypothesized to contain important information about brain architecture, little is known about their origin or about their relationship to active cognitive processes such as conscious awareness and monitoring of the environment. In this study, we have investigated the amplitude and spatiotemporal characteristics of resting-state activity patterns and their dependence on the subjects' alertness. For this purpose, BOLD fMRI was performed at 3.0 T on 12 normal subjects using a visual stimulation protocol, followed by a 27 min rest period, during which subjects were allowed to fall asleep. In subjects who were asleep at the end of the scan, we found (a) a higher amplitude of BOLD signal fluctuation during rest compared with subjects who were awake at the end of the scan; (b) spatially independent patterns of correlated activity that involve all of gray matter, including deep brain nuclei; (c) many patterns that were consistent across subjects; (d) that average percentage levels of fluctuation in visual cortex (VC) and whole brain were higher in subjects who were asleep (up to 1.71% and 1.16%, respectively) than in those who were awake (up to 1.15% and 0.96%) at the end of the scan and were comparable with those levels evoked by intense visual stimulation (up to 1.85% and 0.76% for two subject groups); (e) no confirmation of correlation, positive or negative, between thalamus and VC found in earlier studies. These findings suggest that resting-state activity continues during sleep and does not require active cognitive processes or conscious awareness.     

Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

Functional neural network analysis in frontotemporal dementia and Alzheimer's disease using EEG and graph theory

Background  

Although a large body of knowledge about both brain structure and function has been gathered over the last decades, we still have a poor understanding of their exact relationship. Graph theory provides a method to study the relation between network structure and function, and its application to neuroscientific data is an emerging research field. We investigated topological changes in large-scale functional brain networks in patients with Alzheimer's disease (AD) and frontotemporal lobar degeneration (FTLD) by means of graph theoretical analysis of resting-state EEG recordings. EEGs of 20 patients with mild to moderate AD, 15 FTLD patients, and 23 non-demented individuals were recorded in an eyes-closed resting-state. The synchronization likelihood (SL), a measure of functional connectivity, was calculated for each sensor pair in 0.5–4 Hz, 4–8 Hz, 8–10 Hz, 10–13 Hz, 13–30 Hz and 30–45 Hz frequency bands. The resulting connectivity matrices were converted to unweighted graphs, whose structure was characterized with several measures: mean clustering coefficient (local connectivity), characteristic path length (global connectivity) and degree correlation (network 'assortativity'). All results were normalized for network size and compared with random control networks.     

Stochastic qualifiers of epileptic brain dynamics

We evaluate the capability of reconstructing Fokker-Planck equations for an improved characterization of electroencephalographic (EEG) recordings from epilepsy patients. We derive stochastic qualifiers of brain dynamics that are based on specific characteristics of the Kramers-Moyal coefficients estimated from the EEG. Analyzing long-lasting multichannel EEG recordings from eight patients suffering from focal epilepsies we show that particularly the stochastic part of the dynamics can yield valuable information for diagnostic purposes.     

Estimation of the effective and functional human cortical connectivity with structural equation modeling and directed transfer function applied to high-resolution EEG

Different brain imaging devices are presently available to provide images of the human functional cortical activity, based on hemodynamic, metabolic or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions are interconnected. The concept of brain connectivity plays a central role in the neuroscience, and different definitions of connectivity, functional and effective, have been adopted in literature. While the functional connectivity is defined as the temporal coherence among the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally among cortical sites. The structural equation modeling (SEM) is the most used method to estimate effective connectivity in neuroscience, and its typical application is on data related to brain hemodynamic behavior tested by functional magnetic resonance imaging (fMRI), whereas the directed transfer function (DTF) method is a frequency-domain approach based on both a multivariate autoregressive (MVAR) modeling of time series and on the concept of Granger causality.

This study presents advanced methods for the estimation of cortical connectivity by applying SEM and DTF on the cortical signals estimated from high-resolution electroencephalography (EEG) recordings, since these signals exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. To estimate correctly the cortical signals, we used a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual MRI, a distributed source model and a regularized linear inverse source estimates of cortical current density. Before the application of SEM and DTF methodology to the cortical waveforms estimated from high-resolution EEG data, we performed a simulation study, in which different main factors (signal-to-noise ratio, SNR, and simulated cortical activity duration, LENGTH) were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). The statistical analysis returned that during simulations, both SEM and DTF estimators were able to correctly estimate the imposed connectivity patterns under reasonable operative conditions, that is, when data exhibit an SNR of at least 3 and a LENGTH of at least 75 s of nonconsecutive EEG recordings at 64 Hz of sampling rate.

Hence, effective and functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high-resolution EEG techniques and linear inverse estimation with SEM or DTF methods. We conclude that the estimation of cortical connectivity can be performed not only with hemodynamic measurements, but also with EEG signals treated with advanced computational techniques.     

A Comparative Study of Functional Connectivity Measures for Brain Network Analysis in the Context of AD Detection with EEG

This work addresses brain network analysis considering different clinical severity stages of cognitive dysfunction, based on resting-state electroencephalography (EEG). We use a cohort acquired in real-life clinical conditions, which contains EEG data of subjective cognitive impairment (SCI) patients, mild cognitive impairment (MCI) patients, and Alzheimer’s disease (AD) patients. We propose to exploit an epoch-based entropy measure to quantify the connectivity links in the networks. This entropy measure relies on a refined statistical modeling of EEG signals with Hidden Markov Models, which allow a better estimation of the spatiotemporal characteristics of EEG signals. We also propose to conduct a comparative study by considering three other measures largely used in the literature: phase lag index, coherence, and mutual information. We calculated such measures at different frequency bands and computed different local graph parameters considering different proportional threshold values for a binary network analysis. After applying a feature selection procedure to determine the most relevant features for classification performance with a linear Support Vector Machine algorithm, our study demonstrates the effectiveness of the statistical entropy measure for analyzing the brain network in patients with different stages of cognitive dysfunction.     

静息态脑电信号动态功能连接分析

静息态脑功能连接分析是近年来脑研究的一个热点问题, 对于某些脑疾病的诊断及成因理解具有重要意义. 已有的脑功能连接研究基本上都假设功能连接网络在一段时间内是稳定不变的, 但越来越多的证据表明它应该是随时间动态变化的. 对25名被试睁眼和闭眼状态的64电极脑电生理信号, 采用独立成分分析、滑动时间窗、低分辨率脑电断层溯源、图论等方法和技术进行动态功能连接分析, 展现了睁眼和闭眼两种基线状态下视觉网络、默认网络等功能连接网络随时间的动态变化, 并对动态连接矩阵进行主成分分析得到了在整个时间段内具有代表意义的功能连接模式. 该结论支持和补充了传统稳态脑功能连接的研究, 也将为相关实验设计以及脑电信号临床研究提供基线选择依据.     

Nonlinear statistical properties of fMRI signals in human visual cortex during resting-state

The purpose of the current work is to study nonlinear dynamics in neuronal activity within human brain visual cortex based on blood-oxygen-level dependent (BOLD) contrast imaging. In particular, based on functional magnetic resonance imaging (fMRI) signals, measures of fractality, complexity, and state disorder are estimated from central and peripheral eccentricity bands across three visual areas. Statistical results from analysis of 48750 resting-state fMRI signals show evidence that nonlinear dynamics of neuronal activity in resting-state in central and peripheral eccentricity bands of human visual cortex are persistent. However, they exhibit heterogeneous variability across eccentricity bands and visual areas. Also, information content in first visual area is more ordered than in the second one, whilst information content in the third visual area is the least ordered. These interesting nonlinear statistical properties are a further step toward understanding neuronal activity and nonlinear dynamics in human brain visual cortex.     

Preoperative localization of the sensorimotor area using independent component analysis of resting-state fMRI

Analysis of resting-state functional magnetic resonance imaging (fMRI) data is based on detecting low-frequency signal fluctuations in functionally connected brain areas. These synchronous fluctuations in resting-state networks have been observed in several studies with healthy subjects. In this study, we explored if independent component analysis (ICA) can be used to localize the sensorimotor area from resting-state fMRI data in patients with brain tumors. Finger-tapping activation task and resting-state blood-oxygenation-level-dependent fMRI data were acquired from 8 patients with brain tumors and 10 healthy volunteers. Sensorimotor task independent components (IC_task) were used to verify resting-state independent components (IC_rest) individually. In addition, sensorimotor IC_rests were compared between the groups and no significant differences were detected in volume, spatial correlation or temporal correlation. These results show that it is possible to localize a sensorimotor area from resting-state data using ICA in patients with brain tumors. This offers a complementary method for assessing the sensorimotor area in subjects with brain tumors who have difficulties in performing motor paradigms.     

Sleep disturbances in highly stress reactive mice: Modeling endophenotypes of major depression

Background  

Neuronal mechanisms underlying affective disorders such as major depression (MD) are still poorly understood. By selectively breeding mice for high (HR), intermediate (IR), or low (LR) reactivity of the hypothalamic-pituitary-adrenocortical (HPA) axis, we recently established a new genetic animal model of extremes in stress reactivity (SR). Studies characterizing this SR mouse model on the behavioral, endocrine, and neurobiological levels revealed several similarities with key endophenotypes observed in MD patients. HR mice were shown to have changes in rhythmicity and sleep measures such as rapid eye movement sleep (REMS) and non-REM sleep (NREMS) as well as in slow wave activity, indicative of reduced sleep efficacy and increased REMS. In the present study we were interested in how far a detailed spectral analysis of several electroencephalogram (EEG) parameters, including relevant frequency bands, could reveal further alterations of sleep architecture in this animal model. Eight adult males of each of the three breeding lines were equipped with epidural EEG and intramuscular electromyogram (EMG) electrodes. After recovery, EEG and EMG recordings were performed for two days.     

Simultaneous EEG and fMRI in the macaque monkey at 4.7 Tesla

Simultaneous electroencephalography (EEG)/functional magnetic resonance imaging (fMRI) acquisition can identify the brain networks involved in generating specific EEG patterns. Yet, the combination of these methodologies is hampered by strong artifacts that arise due to electromagnetic interference during magnetic resonance (MR) image acquisition. Here, we report corrections of the gradient-induced artifact in phantom measurements and in experiments with an awake behaving macaque monkey during fMRI acquisition at a magnetic field strength of 4.7 T. Ninety-one percent of the amplitude of a 10 microV, 10 Hz phantom signal could successfully be recovered without phase distortions. Using this method, we were able to extract the monkey EEG from scalp recordings obtained during MR image acquisition. Visual evoked potentials could also be reliably identified. In conclusion, simultaneous EEG/fMRI acquisition is feasible in the macaque monkey preparation at 4.7 T and holds promise for investigating the neural processes that give rise to particular EEG patterns.     

Combination of BOLD-fMRI and VEP recordings for spin-echo MRI detection of primary magnetic effects caused by neuronal currents

In the present paper, for the first time, the feasibility to detect primary magnetic field changes caused by neuronal activity in vivo by spin-echo (SE) magnetic resonance imaging (MRI) is investigated. The detection of effects more directly linked to brain activity than secondary hemodynamic–metabolic changes would enable the study of brain function with improved specificity. However, the detection of neuronal currents by MRI is hampered by such accompanying hemodynamic changes. Therefore, SE image acquisition, rather than gradient-echo (GE) image acquisition, was preferred in the present work since the detection of primary neuronal and not blood oxygenation level-dependent (BOLD)-related effects may be facilitated by this approach. First of all, a precise spatiotemporal synchronization of image acquisition with the neuronal event had to be performed to avoid refocusing of the dephasing phenomenon during the course of the SE sequence. At this aim, we propose the combined use of visual evoked potential (VEP) recordings and BOLD-fMRI measurements prior to SE MRI scanning. Moreover, we exemplify by theory and experimentation how the control of artefactual signal changes due to BOLD and movement effects may be further improved by the experimental design. Finally, results from a pilot study using the proposed combination of VEP recordings and MRI techniques are reported, suggesting the feasibility of this method.     

Approximate Entropy of Brain Network in the Study of Hemispheric Differences

Human brain, a dynamic complex system, can be studied with different approaches, including linear and nonlinear ones. One of the nonlinear approaches widely used in electroencephalographic (EEG) analyses is the entropy, the measurement of disorder in a system. The present study investigates brain networks applying approximate entropy (ApEn) measure for assessing the hemispheric EEG differences; reproducibility and stability of ApEn data across separate recording sessions were evaluated. Twenty healthy adult volunteers were submitted to eyes-closed resting EEG recordings, for 80 recordings. Significant differences in the occipital region, with higher values of entropy in the left hemisphere than in the right one, show that the hemispheres become active with different intensities according to the performed function. Besides, the present methodology proved to be reproducible and stable, when carried out on relatively brief EEG epochs but also at a 1-week distance in a group of 36 subjects. Nonlinear approaches represent an interesting probe to study the dynamics of brain networks. ApEn technique might provide more insight into the pathophysiological processes underlying age-related brain disconnection as well as for monitoring the impact of pharmacological and rehabilitation treatments.     

Modulation of electroencephalograph activity by manual acupuncture stimulation in healthy subjects:An autoregressive spectral analysis

To investigate whether and how manual acupuncture(MA) modulates brain activities,we design an experiment where acupuncture at acupoint ST36 of the right leg is used to obtain electroencephalograph(EEG) signals in healthy subjects.We adopt the autoregressive(AR) Burg method to estimate the power spectrum of EEG signals and analyze the relative powers in delta(0 Hz-4 Hz),theta(4 Hz-8 Hz),alpha(8 Hz-13 Hz),and beta(13 Hz-30 Hz) bands.Our results show that MA at ST36 can significantly increase the EEG slow wave relative power(delta band) and reduce the fast wave relative powers(alpha and beta bands),while there are no statistical differences in theta band relative power between different acupuncture states.In order to quantify the ratio of slow to fast wave EEG activity,we compute the power ratio index.It is found that the MA can significantly increase the power ratio index,especially in frontal and central lobes.All the results highlight the modulation of brain activities with MA and may provide potential help for the clinical use of acupuncture.The proposed quantitative method of acupuncture signals may be further used to make MA more standardized.     

基于自适应模板法的脑电信号转移熵分析

脑电信号是一种产生机理相当复杂且非常微弱的随机信号, 综合反映了大脑组织的脑电活动及大脑的功能状态. 由于脑电信号的微弱性, 传统的基本模板方法在脑电信号分析上得到了良好的应用. 为进一步提升分析脑电信号的性能, 提出了一种新的基于自适应模板的转移熵方法并分析了青少年脑电与成年人脑电信号. 结果表明: 对于青少年脑电还是成年人脑电, 与基本模板法相比, 基于自适应模板法的转移熵可以更显著地表示脑电信号的耦合作用, 并且具有更好的区分度, 这将能更好地捕捉到信号中的动态信息、系统动力学复杂性的改变. 同时, 该方法将更有利于医学临床诊断的辅助检测, 对脑电信号是否处于病理状态的诊断提供了新的更好的判断依据.     

Comparison of analytical strategies for EEG-correlated fMRI data in patients with epilepsy

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.     

EEG Fractal Analysis Reflects Brain Impairment after Stroke

Stroke is the commonest cause of disability. Novel treatments require an improved understanding of the underlying mechanisms of recovery. Fractal approaches have demonstrated that a single metric can describe the complexity of seemingly random fluctuations of physiological signals. We hypothesize that fractal algorithms applied to electroencephalographic (EEG) signals may track brain impairment after stroke. Sixteen stroke survivors were studied in the hyperacute (<48 h) and in the acute phase (∼1 week after stroke), and 35 stroke survivors during the early subacute phase (from 8 days to 32 days and after ∼2 months after stroke): We compared resting-state EEG fractal changes using fractal measures (i.e., Higuchi Index, Tortuosity) with 11 healthy controls. Both Higuchi index and Tortuosity values were significantly lower after a stroke throughout the acute and early subacute stage compared to healthy subjects, reflecting a brain activity which is significantly less complex. These indices may be promising metrics to track behavioral changes in the very early stage after stroke. Our findings might contribute to the neurorehabilitation quest in identifying reliable biomarkers for a better tailoring of rehabilitation pathways.     

Extreme events in epileptic EEG of rodents after ischemic stroke

Extreme events are observed on electroencephalographic (EEG) recordings of rodents with induced ischemic stroke. These events represent themselves as post-stroke epileptic seizures in the form of spontaneous high-amplitude oscillations, which appear during the first 2–3 h after induced focal cerebral ischemia. The analysis of the EEG time-frequency structure reveals these extreme events as a sharp sudden growth of the wavelet energy in a particular frequency band, while the energy in the resting part of the power spectrum remains normal. The distinguished features of the extreme events are used for detection and quantification of the pathological brain activity.     

Distributed analysis of simultaneous EEG-fMRI time-series: modeling and interpretation issues

Functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) represent brain activity in terms of a reliable anatomical localization and a detailed temporal evolution of neural signals. Simultaneous EEG-fMRI recordings offer the possibility to greatly enrich the significance and the interpretation of the single modality results because the same neural processes are observed from the same brain at the same time. Nonetheless, the different physical nature of the measured signals by the two techniques renders the coupling not always straightforward, especially in cognitive experiments where spatially localized and distributed effects coexist and evolve temporally at different temporal scales.     

Automatic localization of EEG electrode markers within 3D MR data

The electrical activity of the brain can be monitored using ElectroEncephaloGraphy (EEG). From the positions of the EEG electrodes, it is possible to localize focal brain activity. Thereby, the accuracy of the localization strongly depends on the accuracy with which the positions of the electrodes can be determined. In this work, we present an automatic, simple, and accurate scheme that detects EEG electrode markers from 3D MR data of the human head.