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 conditional Granger causality model approach for group analysis in functional magnetic resonance imaging

Granger causality model (GCM) derived from multivariate vector autoregressive models of data has been employed to identify effective connectivity in the human brain with functional magnetic resonance imaging (fMRI) and to reveal complex temporal and spatial dynamics underlying a variety of cognitive processes. In the most recent fMRI effective connectivity measures, pair-wise GCM has commonly been applied based on single-voxel values or average values from special brain areas at the group level. Although a few novel conditional GCM methods have been proposed to quantify the connections between brain areas, our study is the first to propose a viable standardized approach for group analysis of fMRI data with GCM. To compare the effectiveness of our approach with traditional pair-wise GCM models, we applied a well-established conditional GCM to preselected time series of brain regions resulting from general linear model (GLM) and group spatial kernel independent component analysis of an fMRI data set in the temporal domain. Data sets consisting of one task-related and one resting-state fMRI were used to investigate connections among brain areas with the conditional GCM method. With the GLM-detected brain activation regions in the emotion-related cortex during the block design paradigm, the conditional GCM method was proposed to study the causality of the habituation between the left amygdala and pregenual cingulate cortex during emotion processing. For the resting-state data set, it is possible to calculate not only the effective connectivity between networks but also the heterogeneity within a single network. Our results have further shown a particular interacting pattern of default mode network that can be characterized as both afferent and efferent influences on the medial prefrontal cortex and posterior cingulate cortex. These results suggest that the conditional GCM approach based on a linear multivariate vector autoregressive model can achieve greater accuracy in detecting network connectivity than the widely used pair-wise GCM, and this group analysis methodology can be quite useful to extend the information obtainable in fMRI.     

Cortex-based independent component analysis of fMRI time series

The cerebral cortex is the main target of analysis in many functional magnetic resonance imaging (fMRI) studies. Since only about 20% of the voxels of a typical fMRI data set lie within the cortex, statistical analysis can be restricted to the subset of the voxels obtained after cortex segmentation. While such restriction does not influence conventional univariate statistical tests, it may have a substantial effect on the performance of multivariate methods.

Here, we describe a novel approach for data-driven analysis of single-subject fMRI time series that combines techniques for the segmentation and reconstruction of the cortical surface of the brain and the spatial independent component analysis (sICA) of the functional time courses (TCs). We use the mesh of the white matter/gray matter boundary, automatically reconstructed from high-spatial-resolution anatomical MR images, to limit the sICA decomposition of a coregistered functional time series to those voxels which are within a specified region with respect to the cortical sheet (cortex-based ICA, or cbICA). We illustrate our analysis method in the context of fMRI blocked and event-related experimental designs and in an fMRI experiment with perceptually ambiguous stimulation, in which an a priori specification of the stimulation protocol is not possible.

A comparison between cbICA and conventional hypothesis-driven statistical methods shows that cortical surface maps and component TCs blindly obtained with cbICA reliably reflect task-related spatiotemporal activation patterns. Furthermore, the advantages of using cbICA when the specification of a temporal model of the expected hemodynamic response is not straightforward are illustrated and discussed. A comparison between cbICA and anatomically unconstrained ICA reveals that — beside reducing computational demand — the cortex-based approach improves the fitting of the ICA model in the gray matter voxels, the separation of cortical components and the estimation of their TCs, particularly in the case of fMRI data sets with a complex spatiotemporal statistical structure.     

The effect of sensorimotor activation on functional connectivity mapping with MRI

The correlations in the fluctuations in the blood oxygenation level-dependent (BOLD) MRI signal between anatomically distinct regions of the cortex that are known components of functional systems have been previously studied as possible indicators of functional connectivity. The objective of this study was to examine the effect of sensorimotor brain activity, as assessed by task-based functional magnetic resonance imaging (fMRI), on functional connectivity indices in the same region. Regions of activation for sequential finger motion were determined using a task-based, block-design fMRI study. Functional connectivity measurements based on interregional correlations were acquired at rest and during continuous, sequential finger motion. Connectivity indices were determined using normalized mean correlations within and between three regions of interest activated for the finger motion task. Connectivity indices were also determined for a control region that was not activated for the task. Continuous motor tasks performed during BOLD measurements did not significantly affect the functional connectivity as compared to the connectivity at rest within or between regions known to be activated by the task. However, there appeared to be a trend suggesting a slight reduction in connectivity indices during the motor task. The connectivity within and between those areas not activated for the task remained unchanged between conditions. These results suggest that in the motor system investigated, the recruitment of neurons to perform a specific task may moderately reduce the degree of hemodynamic coupling within and between regions.     

Investigating the neural basis for fMRI-based functional connectivity in a blocked design: application to interregional correlations and psycho-physiological interactions

This paper investigates how well different kinds of fMRI functional connectivity analysis reflect the underlying interregional neural interactions. This is hard to evaluate using real experimental data where such relationships are unknown. Rather, we use a biologically realistic neural model to simulate both neuronal activities and multiregional fMRI data from a blocked design. Because we know how every element in the model is related to every other element, we can compare functional connectivity measurements across different spatial and temporal scales. We focus on (1) psycho-physiological interaction (PPI) analysis, which is a simple brain connectivity method that characterizes the activity in one brain region by the interaction between another region's activity and a psychological factor, and (2) interregional correlation analysis. We investigated the neurobiological underpinnings of PPI using simulated neural activities and fMRI signals generated by a large-scale neural model that performs a visual delayed match-to-sample task. Simulated fMRI data are generated by convolving integrated synaptic activities (ISAs) with a hemodynamic response function. The simulation was done under three task conditions: high-attention, low-attention and a control task ('passive viewing'). We investigated how biological and scanning parameters affect PPI and compared these with functional connectivity measures obtained using correlation analysis. We performed correlational and PPI analyses with three types of time-series data: ISA, fMRI and deconvolved fMRI (which yields estimated neural signals) obtained using a deconvolution algorithm. The simulated ISA can be considered as the 'gold standard' because it represents the underlying neural activity. Our main findings show (1) that evaluating the change in an interregional functional connection using the difference in regression coefficients (as is essentially done in the PPI method) produces results that better reflect the underlying changes in neural interrelationships than does evaluating the functional connectivity difference as a change in correlation coefficient; (2) that using fMRI and deconvolved fMRI data led to similar conclusions in the PPI-based functional connectivity results, and these generally agreed with the nature of the underlying neural interactions; and (3) the functional connectivity correlation measures often led to different conclusions regarding significance for different scanning and hemodynamic parameters, but the significances of the PPI regression parameters were relatively robust. These results highlight the way in which neural modeling can be used to help validate the inferences one can make about functional connectivity based on fMRI data.     

Nonlinear responses of cerebral blood volume, blood flow and blood oxygenation signals during visual stimulation

Hemodynamic-based functional magnetic resonance imaging (fMRI) techniques provide a great utility for noninvasive functional brain mapping. However, because the hemodynamic signals reflect underlying neural activity indirectly, characterization of these signals following brain activation is essential for experimental design and data interpretation. In this report, the linear (or nonlinear) responses to neuronal activation of three hemodynamic parameters based primarily on changes of cerebral blood volume, blood flow and blood oxygenation were investigated by testing these hemodynamic responses' additivity property. Using a recently developed fMRI technique that acquires vascular space occupancy (VASO), arterial spin labeling (ASL) perfusion and blood oxygenation level-dependent (BOLD) signals simultaneously, the additivity property of the three hemodynamic responses in human visual cortex was assessed using various visual stimulus durations. Experiments on healthy volunteers showed that all three hemodynamic-weighted signals responded nonlinearly to stimulus durations less than 4 s, with the degree of nonlinearity becoming more severe as the stimulus duration decreased. Vascular space occupancy and ASL perfusion signals showed similar nonlinearity properties, whereas the BOLD signal was the most nonlinear. These data suggest that caution should be taken in the interpretation of hemodynamic-based signals in fMRI.     

A comparison of Gamma and Gaussian dynamic convolution models of the fMRI BOLD response

Blood oxygenation level-dependent (BOLD) contrast-based functional magnetic resonance imaging (fMRI) has been widely utilized to detect brain neural activities and great efforts are now stressed on the hemodynamic processes of different brain regions activated by a stimulus. The focus of this paper is the comparison of Gamma and Gaussian dynamic convolution models of the fMRI BOLD response. The convolutions are between the perfusion function of the neural response to a stimulus and a Gaussian or Gamma function. The parameters of the two models are estimated by a nonlinear least-squares optimal algorithm for the fMRI data of eight subjects collected in a visual stimulus experiment. The results show that the Gaussian model is better than the Gamma model in fitting the data. The model parameters are different in the left and right occipital regions, which indicate that the dynamic processes seem different in various cerebral functional regions.     

Biophysical and neural basis of resting state functional connectivity: Evidence from non-human primates

Functional MRI (fMRI) has evolved from simple observations of regional changes in MRI signals caused by cortical activity induced by a task or stimulus, to task-free acquisitions of images in a resting state. Such resting state signals contain low frequency fluctuations which may be correlated between voxels, and strongly correlated regions are deemed to reflect functional connectivity within synchronized circuits. Resting state functional connectivity (rsFC) measures have been widely adopted by the neuroscience community, and are being used and interpreted as indicators of intrinsic neural circuits and their functional states in a broad range of applications, both basic and clinical. However, there has been relatively little work reported that validates whether inter-regional correlations in resting state fluctuations of fMRI (rsfMRI) signals actually measure functional connectivity between brain regions, or to establish how MRI data correlate with other metrics of functional connectivity. In this mini-review, we summarize recent studies of rsFC within mesoscopic scale cortical networks (100 μm–10 mm) within a well defined functional region of primary somatosensory cortex (S1), as well as spinal cord and brain white matter in non-human primates, in which we have measured spatial patterns of resting state correlations and validated their interpretation with electrophysiological signals and anatomic connections. Moreover, we emphasize that low frequency correlations are a general feature of neural systems, as evidenced by their presence in the spinal cord as well as white matter. These studies demonstrate the valuable role of high field MRI and invasive measurements in an animal model to inform the interpretation of human imaging studies.     

Integration of fMRI, NIROT and ERP for studies of human brain function

Different methods of assessing human brain function possess specific advantages and disadvantages compared to others, but it is believed that combining different approaches will provide greater information than can be obtained from each alone. For example, functional magnetic resonance imaging (fMRI) has good spatial resolution but poor temporal resolution, whereas the converse is true for electrophysiological recordings (event-related potentials or ERPs). In this review of recent work, we highlight a novel approach to combining these modalities in a manner designed to increase information on the origins and locations of the generators of specific ERPs and the relationship between fMRI and ERP signals. Near infrared imaging techniques have also been studied as alternatives to fMRI and can be readily integrated with simultaneous electrophysiological recordings. Each of these modalities may in principle be also used in so-called steady-state acquisitions in which the correlational structure of signals from the brain may be analyzed to provide new insights into brain function.     

Neurophysiological and functional MRI evidence of reorganization of cortical motor areas in cerebral arteriovenous malformation

Functional magnetic resonance imaging (fMRI) research has shown that brain arteriovenous malformations (AVMs) lead to reorganization of cortical motor areas. Since it is known that blood oxygenation level-dependent signal in fMRI may be influenced by the hemodynamic perturbation associated with the presence of the AVM, in the present study, a combined exploration with fMRI and transcranial magnetic stimulation was performed in a patient with a right rolandic AVM in order to explore the relationship between neuronal and hemodynamic activity. The combined protocol of investigation adopted in this study was able to provide significant information regarding neuronal activity of the different cortical areas that partake to post-lesional reorganization.     

Dissecting cognitive stages with time-resolved fMRI data: a comparison of fuzzy clustering and independent component analysis

In combination with cognitive tasks entailing sequences of sensory and cognitive processes, event-related acquisition schemes allow using functional MRI to examine not only the topography but also the temporal sequence of cortical activation across brain regions (time-resolved fMRI). In this study, we compared two data-driven methods--fuzzy clustering method (FCM) and independent component analysis (ICA)--in the context of time-resolved fMRI data collected during the performance of a newly devised visual imagery task. We analyzed a multisubject fMRI data set using both methods and compared their results in terms of within- and between-subject consistency and spatial and temporal correspondence of obtained maps and time courses. Both FCM and spatial ICA allowed discriminating the contribution of distinct networks of brain regions to the main cognitive stages of the task (auditory perception, mental imagery and behavioural response), with good agreement across methods. Whereas ICA worked optimally on the original time series, averaging with respect to the task onset (and thus introducing some a priori information on the stimulation protocol) was found to be indispensable in the case of FCM. On averaged time series, FCM led to a richer decomposition of the spatio-temporal patterns of activation and allowed a finer separation of the neurocognitive processes subserving the mental imagery task. This study confirms the efficacy of the two examined methods in the data-driven estimation of hemodynamic responses in time-resolved fMRI studies and provides empirical guidelines to their use.     

Functional Magnetic Resonance Imaging Data Analysis by Autoregressive Estimator

The autoregressive (AR) estimator, a non-parametric method, is used to analyze functional magnetic resonance imaging (fMRI) data. The same method has been used, with success, in several other time series data analysis. It uses exclusively the available experimental data points to estimate the most plausible power spectra compatible with the experimental data and there is no need to make any assumption about non-measured points. The time series, obtained from fMRI block paradigm data, is analyzed by the AR method to determine the brain active regions involved in the processing of a given stimulus. This method is considerably more reliable than the fast Fourier transform or the parametric methods. The time series corresponding to each image pixel is analyzed using the AR estimator and the corresponding poles are obtained. The pole distribution gives the shape of power spectra, and the pixels with poles at the stimulation frequency are considered as the active regions. The method was applied in simulated and real data, its superiority is shown by the receiver operating characteristic curves which were obtained using the simulated data.     

静息状态下脑功能连接的磁共振成像研究

静息状态下脑功能连接的磁共振成像研究近年来取得了迅猛发展. 通过对fMRI信号低频涨落成分的同步性分析,可以得到大脑静息态任意脑区的功能连接和多套网络系统,其中“默认网络”的发现可能为人脑固有网络的研究提供新的思路. 而静息态网络与解剖连接之间可能存在的对应,以及在神经精神疾病患者脑中性质和连接的异常改变,使其具有重要的研究和临床应用价值. 该文总结了静息状态功能磁共振成像的主要研究成果,对静息状态脑功能网络的发现和发展、研究方法、各网络及其特点以及在临床方面的应用进行简单的介绍和分析.     

Relationship between neural and hemodynamic signals during spontaneous activity studied with temporal kernel CCA

Functional magnetic resonance imaging (fMRI) based on the so-called blood oxygen level-dependent (BOLD) contrast is a powerful tool for studying brain function not only locally but also on the large scale. Most studies assume a simple relationship between neural and BOLD activity, in spite of the fact that it is important to elucidate how the “when” and “what” components of neural activity are correlated to the “where” of fMRI data. Here we conducted simultaneous recordings of neural and BOLD signal fluctuations in primary visual (V1) cortex of anesthetized monkeys. We explored the neurovascular relationship during periods of spontaneous activity by using temporal kernel canonical correlation analysis (tkCCA). tkCCA is a multivariate method that can take into account any features in the signals that univariate analysis cannot. The method detects filters in voxel space (for fMRI data) and in frequency–time space (for neural data) that maximize the neurovascular correlation without any assumption of a hemodynamic response function (HRF). Our results showed a positive neurovascular coupling with a lag of 4–5 s and a larger contribution from local field potentials (LFPs) in the γ range than from low-frequency LFPs or spiking activity. The method also detected a higher correlation around the recording site in the concurrent spatial map, even though the pattern covered most of the occipital part of V1. These results are consistent with those of previous studies and represent the first multivariate analysis of intracranial electrophysiology and high-resolution fMRI.     

Component structure of event-related fMRI responses in the different neurovascular compartments

In most functional magnetic resonance imaging (fMRI) studies, brain activity is localized by observing changes in the blood oxygenation level-dependent (BOLD) signal that are believed to arise from capillaries, venules and veins in and around the active neuronal population. However, the contribution from veins can be relatively far downstream from active neurons, thereby limiting the ability of BOLD imaging methods to precisely pinpoint neural generators. Hemodynamic measures based on apparent diffusion coefficients (ADCs) have recently been used to identify more upstream functional blood flow changes in the capillaries, arterioles and arteries. In particular, we recently showed that, due to the complementary vascular sensitivities of ADC and BOLD signals, the voxels conjointly activated by both measures may identify the capillary networks of the active neuronal areas. In this study, we first used simultaneously acquired ADC and BOLD functional imaging signals to identify brain voxels activated by ADC only, by both ADC and BOLD and by BOLD only, thereby delineating voxels relatively dominated by the arterial, capillary, and draining venous neurovascular compartments, respectively. We then examined the event-related fMRI BOLD responses in each of these delineated neurovascular compartments, hypothesizing that their event-related responses would show different temporal componentries. In the regions activated by both the BOLD and ADC contrasts, but not in the BOLD-only areas, we observed an initial transient signal reduction (an initial dip), consistent with the local production of deoxyhemoglobin by the active neuronal population. In addition, the BOLD-ADC overlap areas and the BOLD-only areas showed a clear poststimulus undershoot, whereas the compartment activated by only ADC did not show this component. These results indicate that using ADC contrast in conjunction with BOLD imaging can help delineate the various neurovascular compartments, improve the localization of active neural populations, and provide insight into the physiological mechanisms underlying the hemodynamic signals.     

On the use of information theory for the analysis of the relationship between neural and imaging signals

Functional magnetic resonance imaging (fMRI) is a widely used method for studying the neural basis of cognition and of sensory function. A potential problem in the interpretation of fMRI data is that fMRI measures neural activity only indirectly, as a local change of deoxyhemoglobin concentration due to the metabolic demands of neural function. To build correct sensory and cognitive maps in the human brain, it is thus crucial to understand whether fMRI and neural activity convey the same type of information about external correlates. While a substantial experimental effort has been devoted to the simultaneous recordings of hemodynamic and neural signals, so far, the development of analysis methods that elucidate how neural and hemodynamic signals represent sensory information has received less attention. In this article, we critically review why the analytical framework of information theory, the mathematical theory of communication, is ideally suited to this purpose. We review the principles of information theory and explain how they could be applied to the analysis of fMRI and neural signals. We show that a critical advantage of information theory over more traditional analysis paradigms commonly used in the fMRI literature is that it can elucidate, within a single framework, whether an empirically observed correlation between neural and fMRI signals reflects either a similar stimulus tuning or a common source of variability unrelated to the external stimuli. In addition, information theory determines the extent to which these shared sources of stimulus signal and of variability lead fMRI and neural signals to convey similar information about external correlates. We then illustrate the formalism by applying it to the analysis of the information carried by different bands of the local field potential. We conclude by discussing the current methodological challenges that need to be addressed to make the information-theoretic approach more robustly applicable to the simultaneous recordings of neural and imaging data.     

The BOLD response and the gamma oscillations respond differently than evoked potentials: an interleaved EEG-fMRI study

Background  

The integration of EEG and fMRI is attractive because of their complementary precision regarding time and space. But the relationship between the indirect hemodynamic fMRI signal and the more direct EEG signal is uncertain. Event-related EEG responses can be analyzed in two different ways, reflecting two different kinds of brain activity: evoked, i.e. phase-locked to the stimulus, such as evoked potentials, or induced, i.e. non phase-locked to the stimulus such as event-related oscillations. In order to determine which kind of EEG activity was more closely related with fMRI, EEG and fMRI signals were acquired together, while subjects were presented with two kinds of rare events intermingled with frequent distractors. Target events had to be signaled by pressing a button and Novel events had to be ignored.     

On the nature of the BOLD fMRI contrast mechanism

Since its development about 15 years ago, functional magnetic resonance imaging (fMRI) has become the leading research tool for mapping brain activity. The technique works by detecting the levels of oxygen in the blood, point by point, throughout the brain. In other words, it relies on a surrogate signal, resulting from changes in oxygenation, blood volume and flow, and does not directly measure neural activity. Although a relationship between changes in brain activity and blood flow has long been speculated, indirectly examined and suggested and surely anticipated and expected, the neural basis of the fMRI signal was only recently demonstrated directly in experiments using combined imaging and intracortical recordings. In the present paper, we discuss the results obtained from such combined experiments. We also discuss our current knowledge of the extracellularly measured signals of the neural processes that they represent and of the structural and functional neurovascular coupling, which links such processes with the hemodynamic changes that offer the surrogate signal that we use to map brain activity. We conclude by considering applications of invasive MRI, including injections of paramagnetic tracers for the study of connectivity in the living animal and simultaneous imaging and electrical microstimulation.     

A signal subspace approach for modeling the hemodynamic response function in fMRI

Many fMRI analysis methods use a model for the hemodynamic response function (HRF). Common models of the HRF, such as the Gaussian or Gamma functions, have parameters that are usually selected a priori by the data analyst. A new method is presented that characterizes the HRF over a wide range of parameters via three basis signals derived using principal component analysis (PCA). Covering the HRF variability, these three basis signals together with the stimulation pattern define signal subspaces which are applicable to both linear and nonlinear modeling and identification of the HRF and for various activation detection strategies. Analysis of simulated fMRI data using the proposed signal subspace showed increased detection sensitivity compared to the case of using a previously proposed trigonometric subspace. The methodology was also applied to activation detection in both event-related and block design experimental fMRI data using both linear and nonlinear modeling of the HRF. The activated regions were consistent with previous studies, indicating the ability of the proposed approach in detecting brain activation without a priori assumptions about the shape parameters of the HRF. The utility of the proposed basis functions in identifying the HRF is demonstrated by estimating the HRF in different activated regions.     

Impaired functionality of reperfused brain tissue following short transient focal ischemia in rats

Functional magnetic resonance imaging (fMRI) has been applied to study the consequences of transient focal ischemia on neuronal excitability in the rat brain. The experimental paradigm consisted of measuring the changes in local cerebral blood volume (CBV) induced by systemic infusion of the GABA(A) antagonist bicuculline after occlusion of the middle cerebral artery (MCA) for durations of 5, 15, 30 and 60 min using the intraluminal thread model. fMRI studies were carried out 60 min after successful reperfusion of the ischemic territory. Bicuculline-induced dynamic changes in local CBV were assessed in three brain regions: Parietal cortex, caudate putamen and thalamus. The measured CBV response was negatively correlated with the ischemia duration. Additionally, the three regions showed different vulnerability to the transient MCA occlusion, caudate being the most susceptible followed by parietal cortex and thalamus. The fMRI signals weakly correlated with basal CBF and CBV following reperfusion. Our results indicate that fMRI is a sensitive method to assess functional integrity of the brain. Activation maps allow to quantitatively assess the functionally compromized territory at an early stage following the ischemic event prior to the manifestation of pathomorphological changes.