Continuous EEG-fMRI in patients with partial epilepsy and focal interictal slow-wave discharges on EEG

Purpose

To verify whether in patients with partial epilepsy and routine electroenecephalogram (EEG) showing focal interictal slow-wave discharges without spikes combined EEG–functional magnetic resonance imaging (fMRI) would localize the corresponding epileptogenic focus, thus providing reliable information on the epileptic source.

Methods

Eight patients with partial epileptic seizures whose routine scalp EEG recordings on presentation showed focal interictal slow-wave activity underwent EEG–fMRI. EEG data were continuously recorded for 24 min (four concatenated sessions) from 18 scalp electrodes, while fMRI scans were simultaneously acquired with a 1.5-Tesla magnetic resonance imaging (MRI) scanner. After recording sessions and MRI artefact removal, EEG data were analyzed offline. We compared blood oxygen level-dependent (BOLD) signal changes on fMRI with EEG recordings obtained at rest and during activation (with and without focal interictal slow-wave discharges).

Results

In all patients, when the EEG tracing showed the onset of focal slow-wave discharges on a few lateralized electrodes, BOLD-fMRI activation in the corresponding brain area significantly increased. We detected significant concordance between focal EEG interictal slow-wave discharges and focal BOLD activation on fMRI. In patients with lesional epilepsy, the epileptogenic area corresponded to the sites of increased focal BOLD signal.

Conclusions

Even in patients with partial epilepsy whose standard EEGs show focal interictal slow-wave discharges without spikes, EEG–fMRI can visualize related focal BOLD activation thus providing useful information for pre-surgical planning.     

Controlling the false positive rate in fuzzy clustering using randomization: application to fMRI activation detection

Despite its potential advantages for fMRI analysis, fuzzy C-means (FCM) clustering suffers from limitations such as the need for a priori knowledge of the number of clusters, and unknown statistical significance and instability of the results. We propose a randomization-based method to control the false-positive rate and estimate statistical significance of the FCM results. Using this novel approach, we develop an fMRI activation detection method. The ability of the method in controlling the false-positive rate is shown by analysis of false positives in activation maps of resting-state fMRI data. Controlling the false-positive rate in FCM allows comparison of different fuzzy clustering methods, using different feature spaces, to other fMRI detection methods. In this article, using simulation and real fMRI data, we compare a novel feature space that takes the variability of the hemodynamic response function into account (HRF-based feature space) to the conventional cross-correlation analysis and FCM using the cross-correlation feature space. In both cases, the HRF-based feature space provides a greater sensitivity compared to the cross-correlation feature space and conventional cross-correlation analysis. Application of the proposed method to finger-tapping fMRI data, using HRF-based feature space, detected activation in sub-cortical regions, whereas both of the FCM with cross-correlation feature space and the conventional cross-correlation method failed to detect them.     

Comprehensive mathematical simulation of functional magnetic resonance imaging time series including motion-related image distortion and spin saturation effect

There has been vast interest in determining the feasibility of functional magnetic resonance imaging (fMRI) as an accurate method of imaging brain function for patient evaluations. The assessment of fMRI as an accurate tool for activation localization largely depends on the software used to process the time series data. The performance evaluation of different analysis tools is not reliable unless truths in motion and activation are known. Lack of valid truths has been the limiting factor for comparisons of different algorithms. Until now, currently available phantom data do not include comprehensive accounts of head motion. While most fMRI studies assume no interslice motion during the time series acquisition in fMRI data acquired using a multislice and single-shot echo-planar imaging sequence, each slice is subject to a different set of motion parameters. In this study, in addition to known three-dimensional motion parameters applied to each slice, included in the time series computation are geometric distortion from field inhomogeneity and spin saturation effect as a result of out-of-plane head motion. We investigated the effect of these head motion-related artifacts and present a validation of the mapping slice-to-volume (MSV) algorithm for motion correction and activation detection against the known truths. MSV was evaluated, and showed better performance in comparison with other widely used fMRI data processing software, which corrects for head motion with a volume-to-volume realignment method. Furthermore, improvement in signal detection was observed with the implementation of the geometric distortion correction and spin saturation effect compensation features in MSV.     

A quality control method for detecting and suppressing uncorrected residual motion in fMRI studies

Motion correction is an important step in the functional magnetic resonance imaging (fMRI) analysis pipeline. While many studies simply exclude subjects who are estimated to have moved beyond an arbitrary threshold, there exists no objective method for determining an appropriate threshold. Furthermore, any criterion based only upon motion estimation ignores the potential for proper realignment. The method proposed here uses unsupervised learning (specifically k-means clustering) on features derived from the mean square derivative (MSD) of the signal before and after realignment to identify problem data. These classifications are refined through analysis of correlation between subject activation maps and the mean activation map, as well as the relationship between tasking and motion as measured through regression of the canonical hemodynamic response functions to fit both estimated motion parameters and MSD. The MSD is further used to identify specific scans containing residual motion, data which is suppressed by adding nuisance regressors to the general linear model; this statistical suppression is performed for identified problem subjects, but has potential for use over all subjects. For problem subjects, our results show increased hemodynamic activity more consistent with group results; that is, the addition of nuisance regressors resulted in a doubling of the correlation between the activation map for the problem subjects and the activation map for all subjects. The proposed method should be useful in helping fMRI researchers make more efficient use of their data by reducing the need to exclude entire subjects from studies and thus collect new data to replace excluded subjects.     

Evaluation and optimization of fMRI single-subject processing pipelines with NPAIRS and second-level CVA

In functional magnetic resonance imaging (fMRI) analysis, although the univariate general linear model (GLM) is currently the dominant approach to brain activation detection, there is growing interest in multivariate approaches such as principal component analysis, canonical variate analysis (CVA), independent component analysis and cluster analysis, which have the potential to reveal neural networks and functional connectivity in the brain. To understand the effect of processing options on performance of multivariate model-based fMRI processing pipelines with real fMRI data, we investigated the impact of commonly used fMRI preprocessing steps and optimized the associated multivariate CVA-based, single-subject processing pipelines with the NPAIRS (nonparametric prediction, activation, influence and reproducibility resampling) performance metrics [prediction accuracy and statistical parametric image (SPI) reproducibility] on the Fiswidgets platform. We also compared the single-subject SPIs of univariate GLM with multivariate CVA-based processing pipelines from SPM, FSL.FEAT, NPAIRS.GLM and NPAIRS.CVA software packages (or modules) using a novel second-level CVA. We found that for the block-design data, (a) slice timing correction and global intensity normalization have little consistent impact on the fMRI processing pipeline, but spatial smoothing, temporal detrending or high-pass filtering, and motion correction significantly improved pipeline performance across all subjects; (b) the combined optimization of spatial smoothing, temporal detrending and CVA model parameters on average improved between-subject reproducibility; and (c) the most important pipeline choices include univariate or multivariate statistical models and spatial smoothing. This study suggests that considering options other than simply using GLM with a fixed spatial filter may be of critical importance in determining activation patterns in BOLD fMRI studies.     

Mapping of spikes, slow waves, and motor tasks in a patient with malformation of cortical development using simultaneous EEG and fMRI

We report on the simultaneous and continuous acquisition of EEG and functional MRI data in a patient with a left hemiparesis and focal epilepsy secondary to malformation of cortical development in the right hemisphere. EEG-triggered fMRI localization was previously demonstrated in this patient. In the experiments reported here, 322 spikes maximum at electrode C4 and 126 focal slow waves were identified offline. A hierarchy of models was explored in order to assess the relative contributions of each type of EEG event. Modeling the BOLD response to C4 spikes alone showed an area of activation within the large malformation, adjacent to the area of infolding cortex. However, also modeling slow-waves gave rise to a broader and stronger activation, suggesting that the generators overlap. Motor mapping of the right hand showed activation in the left sensorimotor cortex; left-hand tapping led to a more diffuse area of activation, displaced superiorly into the superior frontal gyrus, and a small area of activation within the lesion. In conclusion, continuous EEG-fMRI is useful to compare the functional mapping of epileptiform activity and eloquent cortices in individual patients.     

High-resolution, multiple gradient-echo functional MRI at 1.5 T

A multiple gradient echo, high resolution imaging method is proposed to better visualize different sources of activation in functional magnetic resonance imaging (fMRI) experiments. Eight echoes are collected from 30 ms to 205 ms with an echo spacing of 25 ms. All echoes show significant activation, but each echo reveals its own pattern of activation. From this variability, it appears that large vessel contributions can be separated from small vessel contributions using a fuzzy cluster analysis across echo times. The results demonstrate the importance of a multiple gradient echo data acquisition approach in localizing various vascular contributions to brain activation in fMRI.     

Flexible statistical modelling detects clinical functional magnetic resonance imaging activation in partially compliant subjects

Clinical functional magnetic resonance imaging (fMRI) occasionally fails to detect significant activation, often due to variability in task performance. The present study seeks to test whether a more flexible statistical analysis can better detect activation, by accounting for variance associated with variable compliance to the task over time. Experimental results and simulated data both confirm that even at 80% compliance to the task, such a flexible model outperforms standard statistical analysis when assessed using the extent of activation (experimental data), goodness of fit (experimental data), and area under the operator characteristic curve (simulated data). Furthermore, retrospective examination of 14 clinical fMRI examinations reveals that in patients where the standard statistical approach yields activation, there is a measurable gain in model performance in adopting the flexible statistical model, with little or no penalty in lost sensitivity. This indicates that a flexible model should be considered, particularly for clinical patients who may have difficulty complying fully with the study task.     

A kernel machine-based fMRI physiological noise removal method

Functional magnetic resonance imaging (fMRI) technique with blood oxygenation level dependent (BOLD) contrast is a powerful tool for noninvasive mapping of brain function under task and resting states. The removal of cardiac- and respiration-induced physiological noise in fMRI data has been a significant challenge as fMRI studies seek to achieve higher spatial resolutions and characterize more subtle neuronal changes. The low temporal sampling rate of most multi-slice fMRI experiments often causes aliasing of physiological noise into the frequency range of BOLD activation signal. In addition, changes of heartbeat and respiration patterns also generate physiological fluctuations that have similar frequencies with BOLD activation. Most existing physiological noise-removal methods either place restrictive limitations on image acquisition or utilize filtering or regression based post-processing algorithms, which cannot distinguish the frequency-overlapping BOLD activation and the physiological noise. In this work, we address the challenge of physiological noise removal via the kernel machine technique, where a nonlinear kernel machine technique, kernel principal component analysis, is used with a specifically identified kernel function to differentiate BOLD signal from the physiological noise of the frequency. The proposed method was evaluated in human fMRI data acquired from multiple task-related and resting state fMRI experiments. A comparison study was also performed with an existing adaptive filtering method. The results indicate that the proposed method can effectively identify and reduce the physiological noise in fMRI data. The comparison study shows that the proposed method can provide comparable or better noise removal performance than the adaptive filtering approach.     

Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study

Signal fluctuations in functional magnetic resonance imaging (fMRI) can result from a number of sources that may have a neuronal, physiologic or instrumental origin. To determine the relative contribution of these sources, we recorded physiological (respiration and cardiac) signals simultaneously with fMRI in human volunteers at rest with their eyes closed. State-of-the-art technology was used including high magnetic field (7 T), a multichannel detector array and high-resolution (3 mm³) echo-planar imaging. We investigated the relative contribution of thermal noise and other sources of variance to the observed fMRI signal fluctuations both in the visual cortex and in the whole brain gray matter. The following sources of variance were evaluated separately: low-frequency drifts due to scanner instability, effects correlated with respiratory and cardiac cycles, effects due to variability in the respiratory flow rate and cardiac rate, and other sources, tentatively attributed to spontaneous neuronal activity. We found that low-frequency drifts are the most significant source of fMRI signal fluctuations (3.0% signal change in the visual cortex, TE=32 ms), followed by spontaneous neuronal activity (2.9%), thermal noise (2.1%), effects due to variability in physiological rates (respiration 0.9%, heartbeat 0.9%), and correlated with physiological cycles (0.6%). We suggest the selection and use of four lagged physiological noise regressors as an effective model to explain the variance related to fluctuations in the rates of respiration volume change and cardiac pulsation. Our results also indicate that, compared to the whole brain gray matter, the visual cortex has higher sensitivity to changes in both the rate of respiration and the spontaneous resting-state activity. Under the conditions of this study, spontaneous neuronal activity is one of the major contributors to the measured fMRI signal fluctuations, increasing almost twofold relative to earlier experiments under similar conditions at 3 T.     

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.     

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.     

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.     

A data post-processing protocol for dynamic MRI data to discriminate brain activity from global physiological effects

All fMRI techniques measure stimulus induced focal metabolic and physiological changes in activated brain areas. During the entire fMRI experiment it is necessary to maintain the general physiological condition of the subject as stable as possible. This is not always an easy task. The typical block design in standard fMRI experiments minimizes most of the problems related with these general physiological changes. However in some fMRI experiments, like pharmacological MRI, the experimental setup makes the use of a blocked design impossible. Therefore signal correction algorithms have been developed to correct for these physiological signal instabilities. These algorithms often require elaborate calculation efforts and the data interpretation is often very difficult if no prior knowledge on the nature of the changes exists.In this work we present an algorithm, which has the advantage of being low in calculation effort and the resulting data after correction are easy to interpret. It makes use of a datafit between the general physiological and focal activation related signal changes to eliminate the generalized effects. This algorithm has been tested on simulated and experimentally obtained signal traces suffering both from substantial general signal changes overwhelming the smaller focal activation induced signal changes.     

Empirical optimization of ASL data analysis using an ASL data processing toolbox: ASLtbx

Arterial spin labeling (ASL) perfusion fMRI data differ in important respects from the more familiar blood oxygen level-dependent (BOLD) fMRI data and require specific processing strategies. In this paper, we examined several factors that may influence ASL data analysis, including data storage bit resolution, motion correction, preprocessing for cerebral blood flow (CBF) calculations and nuisance covariate modeling. Continuous ASL data were collected at 3 T from 10 subjects while they performed a simple sensorimotor task with an epoch length of 48 s. These data were then analyzed using systematic variations of the factors listed above to identify the approach that yielded optimal signal detection for task activation. Improvements in statistical power were found for use of at least 10 bits for data storage at 3 T. No significant difference was found in motor cortex regarding using simple subtraction or sinc subtraction, but the former presented minor but significantly (P<.024) larger peak t value in visual cortex. While artifactual head motion patterns were observed in synthetic data and background-suppressed ASL data when label/control images were realigned to a common target, independent realignment of label and control images did not yield significant improvements in activation in the sensorimotor data. It was also found that CBF calculations should be performed prior to spatial normalization and that modeling of global fluctuations yielded significantly increased peak t value in motor cortex. The implementation of all ASL data processing approaches is easily accomplished within an open-source toolbox, ASLtbx, and is advocated for most perfusion fMRI data sets.     

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.     

Regional brain activity after prolonged cholinergic enhancement in early Alzheimer's disease

This functional magnetic resonance imaging (fMRI) study examined changes in brain activation after prolonged (20 weeks) and stabilized treatment with the cholinesterase inhibitor galantamine in a small group of patients with very mild Alzheimer's disease (AD). Two cognitive activation paradigms were chosen: one requiring semantic association and the other relying on attention and requiring target detection. A group of age- and education-matched healthy controls was also scanned for comparison. A modest (but not statistically significant) improvement in behavioral scores after treatment was observed in both fMRI tasks. There were brain activation increases in the semantic association task after treatment, and the differences in brain activation present in the comparison of AD patients' baseline images with those of controls were not detectable after treatment. In the target detection task, regions that were activated in the elderly controls but not in the baseline images of the AD group also showed significant activation after treatment. Overall, however, the increases were modest and might reflect the heterogeneity of clinical response to treatment in this small group. Future pharmacological fMRI studies should include clinical response as a factor in the analysis of cholinergic enhancement effects in AD patients.     

Diffusion tensor MRI in temporal lobe epilepsy

The purpose of this study was to investigate the diffusion characteristics of white matter in patients with focal temporal lobe epilepsy (TLE). Diffusion tensor imaging (DTI) was applied to patients and normal controls. Rotationally invariant mean diffusivity and diffusion anisotropy maps were calculated for all subjects. Comparisons between the two groups were performed for several white matter structures. Mean diffusivity and diffusion anisotropy of each selected structure were tested for correlations with age at onset and duration of epilepsy. Significantly lower diffusion anisotropy, and higher diffusivity in directions perpendicular to the axons, was detected in several white matter structures of the patients when compared to the controls. These structures were not located in the temporal lobes. No significant difference in mean diffusivity was detected between the selected structures from the two groups. Diffusion anisotropy was significantly correlated with age at onset of epilepsy in the posterior corpus callosum. Duration of epilepsy was not significantly correlated with the diffusion indices from any of the selected structures. The results of this study suggest that diffusion anisotropy may reveal abnormalities in patients with focal TLE. In addition, these abnormal changes are not necessarily restricted to the temporal lobes but might extend in other brain regions as well. Furthermore, the age at onset of epilepsy may be an important factor in determining the extent of the effect of epilepsy on white matter.     

Functional MR imaging in the awake monkey: effects of motion on dynamic off-resonance and processing strategies

Functional MR imaging of the alert, behaving monkey is being used more and more often to detect activation patterns and guide electrophysiological research investigating the neural basis of behavior. Several labs have reported fMRI data from the awake monkey, but none of them has studied and systematically corrected the effects of monkeys' motion on fMRI time series. In this study, a significant refinement of acquisition and correction strategies is reported that can be used to minimize magnetic susceptibility artifacts induced by respiration and by jaw and body movement. Real-time acquisition of sensor signals (e.g., signals induced by jaw and body movement) and MR navigator data were combined to optimize fMRI signal-correction strategies. Within trials, the artifact-induced off-resonance changes were small and mainly reflected the effects of respiration; between trials, movements caused major changes of global frequency and shim (>20 Hz/cm). Several methods were used to assess the stability of the fMRI series: k-space analysis ('dynamic intensity and off-resonance changes in k-space', dubbed DICK and DORK) and image analysis using a Laplace operator and a center-of-mass metric. The variability between trials made it essential to correct for inter-trial variations. On the other hand, images were sufficiently stable with our approach to perform fMRI evaluations on single trials before averaging of trials. Different motion correction strategies were compared: DORK, McFLIRT (rigid body model with three translations and three rotations) and 2D image alignment based on a center-of-mass detection (in-plane translation). The latter yielded the best results and proved to be fast and robust for intra- and inter-trial alignment. Finally, fMRI in the behaving monkey was tested for spatial and temporal reproducibility on a trial-to-trial basis. Highly activated voxels also displayed good reproducibility between trials. On average, the BOLD amplitude response to a short 3-s visual stimulus was close to 2%.