M/EEG+FMRI BIB

The successful use of functional magnetic resonance imaging (fMRI) as a way of visualizing cortical function depends largely on the important relationships between the signal observed and the underlying neuronal activity that it is believed to represent. Currently, a relatively direct correlation seems to be favoured between fMRI signals and population synaptic activity (including inhibitory and excitatory activity), with a secondary and potentially more variable correlation with cellular action potentials.

Keywords: Action Potentials/physiology ; Animals ; Cerebral Cortex/*physiology ; Cerebrovascular Circulation/*physiology ; Excitatory Postsynaptic Potentials/physiology ; Human ; *Magnetic Resonance Imaging ; Neural Inhibition/physiology ; Neurons/*physiology ; Support, Non-U.S. Gov't ; Synaptic Transmission/*physiology

The interpretation of task-induced functional imaging of the brain is critically dependent on understanding the relationship between observed haemodynamic responses and the underlying neural changes. However, the precise nature of this neurovascular coupling relationship remains unknown. In particular, it is unclear which measure of functional magnetic resonance imaging blood oxygen level dependent (fMRI BOLD) activity is the best correlate of neural activity. We measured the somatosensory evoked potential (SEP) amplitude at the scalp, and fMRI BOLD signal to increases in intensity of contralateral median nerve electrical stimulation in healthy non-anaesthetised subjects. We compared correlation analyses between SEP amplitude and both peak voxel fMRI BOLD percentage signal change and mean voxel fMRI BOLD percentage signal change across a somatosensory cluster, and we also performed a voxel-by-voxel correlation between fMRI BOLD activity and SEP amplitude. We found that fMRI BOLD changes in primary somatosensory cortex correlate significantly with SEP amplitudes, suggesting a linear neurovascular coupling relationship under the conditions investigated. We also found that mean changes across a cluster correlate less well with SEP amplitude than peak voxel levels. This suggests that the area of haemodynamic activity correlating with SEP amplitude is smaller than the entire cluster observed.

Keywords: Adult ; Brain Mapping ; Comparative Study ; Electric Stimulation ; Electrophysiology/*methods ; Evoked Potentials, Somatosensory/*physiology ; Female ; Hemodynamic Processes/physiology ; Human ; *Magnetic Resonance Imaging ; Male ; Nerve Net/physiology ; Oxygen/metabolism ; Somatosensory Cortex/*physiology ; Support, Non-U.S. Gov't

OBJECTIVE: To introduce a new technique for co-registration of Magnetoencephalography (MEG) with magnetic resonance imaging (MRI). We compare the accuracy of a new bite-bar with fixed fiducials to a previous technique whereby fiducial coils were attached proximal to landmarks on the skull. METHODS: A bite-bar with fixed fiducial coils is used to determine the position of the head in the MEG co-ordinate system. Co-registration is performed by a surface-matching technique. The advantage of fixing the coils is that the co-ordinate system is not based upon arbitrary and operator dependent fiducial points that are attached to landmarks (e.g. nasion and the preauricular points), but rather on those that are permanently fixed in relation to the skull. RESULTS: As a consequence of minimizing coil movement during digitization, errors in localization of the coils are significantly reduced, as shown by a randomization test. Displacement of the bite-bar caused by removal and repositioning between MEG recordings is minimal ( approximately 0.5 mm), and dipole localization accuracy of a somatosensory mapping paradigm shows a repeatability of approximately 5 mm. The overall accuracy of the new procedure is greatly improved compared to the previous technique. CONCLUSIONS: The test-retest reliability and accuracy of target localization with the new design is superior to techniques that incorporate anatomical-based fiducial points or coils placed on the circumference of the head.

Keywords: Brain/anatomy & histology ; Comparative Study ; Data Collection ; Equipment Design ; Head ; Human ; *Image Processing, Computer-Assisted ; *Magnetic Resonance Imaging ; *Magnetoencephalography ; Monte Carlo Method ; Posture ; Reproducibility of Results ; Stereotaxic Techniques/*instrumentation/standards

Independent component analysis (ICA) of functional magnetic resonance imaging (fMRI) data is commonly carried out under the assumption that each source may be represented as a spatially fixed pattern of activation, which leads to the instantaneous mixing model. To allow modeling patterns of spatio-temporal dynamics, in particular, the flow of oxygenated blood, we have developed a convolutive ICA approach: spatial complex ICA applied to frequency-domain fMRI data. In several frequency-bands, we identify components pertaining to activity in primary visual cortex (V1) and blood supply vessels. One such component, obtained in the 0.10 Hz band, is analyzed in detail and found to likely reflect flow of oxygenated blood in V1.

The haemodynamic responses to neural activity that underlie the blood-oxygen-level-dependent (BOLD) signal used in functional magnetic resonance imaging (fMRI) of the brain are often assumed to be driven by energy use, particularly in presynaptic terminals or glia. However, recent work has suggested that most brain energy is used to power postsynaptic currents and action potentials rather than presynaptic or glial activity and, furthermore, that haemodynamic responses are driven by neurotransmitter-related signalling and not directly by the local energy needs of the brain. A firm understanding of the BOLD response will require investigation to be focussed on the neural signalling mechanisms controlling blood flow rather than on the locus of energy use.

Keywords: Action Potentials/physiology ; Astrocytes/physiology ; Brain/*blood supply/physiology ; Brain Mapping ; Cerebrovascular Circulation/*physiology ; Energy Metabolism/*physiology ; Human ; Magnetic Resonance Imaging ; Neural Inhibition/physiology ; Presynaptic Terminals/physiology ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S.

Blood oxygenation dependent contrast (BOLD) fMRI is used increasingly to probe connectivity based on temporal correlations between signals from different brain regions. This approach assumes that there is constant local coupling of neuronal activity to the associated BOLD response. Here we test the alternative hypothesis that there is not a fixed relationship between these by determining whether attention modulates apparent neurovascular coupling. Electrical stimulation of the median nerve was applied with and without a concurrent distractor task (serial subtraction). Increasing stimulation intensity increased discomfort ratings ( p<0.001) and was associated with a significant increase in both somatosensory evoked potential (SEP) N20-P25 amplitude and BOLD fMRI response in the contralateral primary (SI) and bilaterally in the secondary somatosensory cortices. Attention to stimulation was reduced during distractor task performance and resulted in an overall trend for reduction in discomfort ( p=0.056), which was significant at the highest stimulation level ( p<0.05). A volume of interest analysis confined to SI confirmed a reduction in BOLD response with distraction ( p<0.001). However, distraction did not measurably affect SEP magnitude. The quantitative relationship between the BOLD fMRI response and the local field potential measured by the early SEP response therefore varies with attentional context. This may be a consequence of differences in either local spatial or temporal signal summation for the two methods. Either interpretation suggests caution in assuming a simple, fixed relationship between local BOLD changes and related electrophysiological activity.

Combined EEG/fMRI recording has been used to localize the generators of EEG events and to identify subject state in cognitive studies and is of increasing interest. However, the large EEG artifacts induced during fMRI have precluded simultaneous EEG and fMRI recording, restricting study design. Removing this artifact is difficult, as it normally exceeds EEG significantly and contains components in the EEG frequency range. We have developed a recording system and an artifact reduction method that reduce this artifact effectively. The recording system has large dynamic range to capture both low-amplitude EEG and large imaging artifact without distortion (resolution 2 microV, range 33.3 mV), 5-kHz sampling, and low-pass filtering prior to the main gain stage. Imaging artifact is reduced by subtracting an averaged artifact waveform, followed by adaptive noise cancellation to reduce any residual artifact. This method was validated in recordings from five subjects using periodic and continuous fMRI sequences. Spectral analysis revealed differences of only 10 to 18% between EEG recorded in the scanner without fMRI and the corrected EEG. Ninety-nine percent of spike waves (median 74 microV) added to the recordings were identified in the corrected EEG compared to 12% in the uncorrected EEG. The median noise after artifact reduction was 8 microV. All these measures indicate that most of the artifact was removed, with minimal EEG distortion. Using this recording system and artifact reduction method, we have demonstrated that simultaneous EEG/fMRI studies are for the first time possible, extending the scope of EEG/fMRI studies considerably.

Keywords: Adult ; Algorithms ; *Artifacts ; Electroencephalography/*methods/statistics & numerical data ; Female ; Human ; Image Processing, Computer-Assisted/*methods/statistics & numerical data ; Magnetic Resonance Imaging/*methods/statistics & numerical data ; Male ; Reproducibility of Results ; Signal Processing, Computer-Assisted

Ballistocardiogram and imaging artifacts cause major interference with simultaneous electroencephalogram (EEG) and functional magnetic resonance imaging (fMRI) recording. In particular, the large amplitude of the imaging artifact precludes easy retrieval of EEG signals during fMRI scanning. Recording with 20,000-Hz digitization rate combined with 3000-Hz low-pass filter revealed the real waveform of the imaging artifact, in which it was elucidated that each artifact peak precisely corresponded to each gradient component and actually had differential waveforms of the original gradient pulses. Based on this finding, to retrieve EEG signal during fMRI acquisition, a blip-type echo planar sequence was modified so that EEG sampling might be performed at every 1000 micros (digitization rate 1000 Hz) exclusively in the period in which the artifact resided around the baseline level. This method, called stepping stone sampling, substantially attenuated the amplitude of the imaging artifact. The remnant of the artifact was subtracted from the averaged artifact waveform. In human studies, alpha activity was successfully retrieved by inspection, and its attenuation/augmentation was observed during eyes open/closed periods. Fast Fourier transform analysis further revealed that even from DC up to 120 Hz, retrieved EEG data during scanning had very similar power distributions to the data retrieved during no scanning, implying the availability of the high-frequency band of the retrieved EEG signals, including even the gamma band.

Keywords: Adult ; Alpha Rhythm ; *Artifacts ; Ballistocardiography/methods ; Brain Mapping/methods ; Cerebral Cortex/*physiology ; Echo-Planar Imaging/methods ; Electroencephalography/*methods ; Female ; Fourier Analysis ; Human ; Image Interpretation, Computer-Assisted/*methods ; Magnetic Resonance Imaging/*methods ; Male ; Phantoms, Imaging ; Reference Values ; Support, Non-U.S. Gov't

The purpose of this study was to investigate the changes in specific absorption rate (SAR) in human-head tissues while using nonmagnetic metallic electroencephalography (EEG) electrodes and leads during magnetic resonance imaging (MRI). A realistic, high resolution (1 mm(3)) head model from individual MRI data was adopted to describe accurately thin tissues, such as bone marrow and skin. The RF power dissipated in the human head was evaluated using the FDTD algorithm. Both surface and bird cage coils were used. The following numbers of EEG electrodes/leads were considered: 16, 31, 62, and 124. Simulations were performed at 128 and 300 MHz. The difference in SAR between the electrodes/leads and no-electrodes conditions was greater with the bird cage coil than with the surface coil. The peak 1 g averaged SAR values were highest at 124 electrodes, increasing to as much as two orders of magnitude (x172.3) at 300 MHz compared to the original value. At 300 MHz, there was a fourfold (x3.6) increase of SAR averaged over the bone marrow, and a sevenfold (x7.4) increase in the skin. At 128 MHz, there was a fivefold (x5.6) increase of whole head SAR. Head models were obtained from two different subjects, with an inter-subject whole head SAR variability of 3%. .

Keywords: Adult ; *Electrodes ; Electroencephalography/*instrumentation ; Human ; Magnetic Resonance Imaging/*instrumentation ; Male ; Support, Non-U.S. Gov't

Magneto- and electroencephalography (MEG/EEG) and functional magnetic resonance imaging (fMRI) provide complementary information about the functional organization of the human brain. An important advantage of MEG/EEG is the millisecond time resolution in detecting electrical activity in the cerebral cortex. The interpretation of MEG/EEG signals, however, is limited by the difficulty of determining the spatial distribution of the neural activity. Functional MRI can help in the MEG/EEG source analysis by suggesting likely locations of activity. We present a geometric interpretation of fMRI-guided inverse solutions in which the MEG/EEG source estimate minimizes a distance to a subspace defined by the fMRI data. In this subspace regularization (SSR) approach, the fMRI bias does not assume preferred amplitudes for MEG/EEG sources, only locations. Characteristic dependence of the source estimates on the regularization parameters is illustrated with simulations. When the fMRI locations match the true MEG/EEG source locations, they serve to bias the underdetermined MEG/EEG inverse solution toward the fMRI loci. Importantly, when the fMRI loci do not match the true MEG/EEG loci, the solution is insensitive to those fMRI loci.

A sudden change in the direction of motion is a particularly salient and relevant feature of visual information. Extensive research has identified cortical areas responsive to visual motion and characterized their sensitivity to different features of motion, such as directional specificity. However, relatively little is known about responses to sudden changes in direction. Electrophysiological data from animals and functional imaging data from humans suggest a number of brain areas responsive to motion, presumably working as a network. Temporal patterns of activity allow the same network to process information in different ways. The present study in humans sought to determine which motion-sensitive areas are involved in processing changes in the direction of motion and to characterize the temporal patterns of processing within this network of brain regions. To accomplish this, we used both magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI). The fMRI data were used as supplementary information in the localization of MEG sources. The change in the direction of visual motion was found to activate a number of areas, each displaying a different temporal behavior. The fMRI revealed motion-related activity in areas MT+ (the human homologue of monkey middle temporal area and possibly also other motion sensitive areas next to MT), a region near the posterior end of the superior temporal sulcus (pSTS), V3A, and V1/V2. The MEG data suggested additional frontal sources. An equivalent dipole model for the generators of MEG signals indicated activity in MT+, starting at 130 ms and peaking at 170 ms after the reversal of the direction of motion, and then again at approximately 260 ms. Frontal activity began 0-20 ms later than in MT+, and peaked approximately 180 ms. Both pSTS and FEF+ showed long-duration activity continuing over the latency range of 200-400 ms. MEG responses in the region of V3A and V1/V2 were relatively small, and peaked at longer latencies than the initial peak in MT+. These data revealed characteristic patterns of activity in this cortical network for processing sudden changes in the direction of visual motion.

Keywords: Adult ; *Brain Mapping ; Cerebral Cortex/*physiology ; *Evoked Potentials, Visual ; Human ; Magnetic Resonance Imaging/*methods ; Magnetoencephalography/*methods ; Male ; Middle Aged ; Motion Perception/*physiology ; Nerve Net/physiology ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S.

Cerebral hemodynamic responses to brief periods of neural activity are delayed and dispersed in time. The specific shape of these responses is of some importance to the design and analysis of blood oxygenation level-dependent (BOLD), functional magnetic resonance imaging (fMRI) experiments. Using fMRI scanning, we examine here the characteristics and variability of hemodynamic responses from the central sulcus in human subjects during an event-related, simple reaction time task. Specifically, we determine the contribution of subject, day, and scanning session (within a day) to variability in the shape of evoked hemodynamic response. We find that while there is significant and substantial variability in the shape of responses collected across subjects, responses collected during multiple scans within a single subject are less variable. The results are discussed in terms of the impact of response variability upon sensitivity and specificity of analyses of event-related fMRI designs.

Keywords: Adult ; Brain/anatomy & histology ; Cerebrovascular Circulation/*physiology ; Female ; Hemodynamic Processes/*physiology ; Human ; Image Processing, Computer-Assisted/*methods ; Magnetic Resonance Imaging ; Male ; Models, Neurological ; Oxygen/*blood ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S.

Combined EEG-fMRI has recently been used to explore the BOLD responses to interictal epileptiform discharges. This study examines whether misspecification of the form of the haemodynamic response function (HRF) results in significant fMRI responses being missed in the statistical analysis. EEG-fMRI data from 31 patients with focal epilepsy were analysed with four HRFs peaking from 3 to 9 sec after each interictal event, in addition to a standard HRF that peaked after 5.4 sec. In four patients, fMRI responses were correlated with gadolinium-enhanced MR angiograms and with EEG data from intracranial electrodes. In an attempt to understand the absence of BOLD responses in a significant group of patients, the degree of signal loss occurring as a result of magnetic field inhomogeneities was compared with the detected fMRI responses in ten patients with temporal lobe spikes. Using multiple HRFs resulted in an increased percentage of data sets with significant fMRI activations, from 45% when using the standard HRF alone, to 62.5%. The standard HRF was good at detecting positive BOLD responses, but less appropriate for negative BOLD responses, the majority of which were more accurately modelled by an HRF that peaked later than the standard. Co-registration of statistical maps with gadolinium-enhanced MRIs suggested that the detected fMRI responses were not in general related to large veins. Signal loss in the temporal lobes seemed to be an important factor in 7 of 12 patients who did not show fMRI activations with any of the HRFs.

We describe two new methods for the inverse problem of electrocardiography. Both employ regularization with multiple constraints, rather than the standard single-constraint regularization. In one method, multiple constraints on the spatial behavior of the solution are used simultaneously. In the other, spatial constraints are used simultaneously with constraints on the temporal behavior of the solution. The specific cases of two spatial constraints and one spatial and one temporal constraint are considered in detail. A new method, the L-Surface, is presented to guide the choice of the required pairs of regularization parameters. In the case when both spatial and temporal regularization are used simultaneously, there is an increased computational burden, and two methods are presented to compute solutions efficiently. The methods are verified by simulations using both dipole sources and measured canine epicardial data.

Keywords: Animals ; Dogs ; Electrocardiography/*methods ; Mathematics ; *Models, Cardiovascular ; *Signal Processing, Computer-Assisted ; Support, U.S. Gov't, Non-P.H.S.

A novel method based on selective detection of rapidly changing DeltaB(0) magnetic fields and suppression of slowly changing DeltaB(0) fields is presented. The ultimate goal of this work is to present a method that may allow detection of transient and subtle changes in B(0) in cortical tissue associated with electrical currents produced by neuronal activity. The method involves the detection of NMR phase changes that occur during a single-shot spin-echo (SE) echo-planar sequence (EPI) echo time. SE EPI effectively rephases all changes in B(0) that occur on a time scale longer than the echo time (TE) and amplifies all DeltaB(0) changes that occur during TE/2. The method was tested on a phantom that contains wires in which current can be modulated. The sensitivity and flexibility of the technique was demonstrated by modulation of the temporal position and duration of the stimuli-evoked transient magnetic field relative to the 180 RF pulse in the imaging sequence-requiring precise stimulus timing. Currently, with this method magnetic field changes as small as 2 x 10(-10) T (200 pT) and lasting for 40 msec can be detected. Implications for direct mapping of brain neuronal activity with MRI are discussed.

Keywords: Brain Mapping/*instrumentation/methods ; Electromagnetic Fields ; Human ; Image Processing, Computer-Assisted ; Magnetic Resonance Imaging/*methods ; Neurons/*physiology ; *Phantoms, Imaging

OBJECTIVES: In this paper, we employed advanced methods for the modeling of human cortical activity related to voluntary right one-digit movements from combined high-resolution electroencepholography (EEG) and functional magnetic resonance imaging (fMRI). METHODS: Multimodal integration between EEG and fMRI data was performed by using realistic head models, a large number of scalp electrodes (128) and the estimation of current density strengths by linear inverse estimation. RESULTS: Increasing of spatial details of the estimated cortical density distributions has been detected by using the proposed integration method with respect to the estimation using EEG data alone. CONCLUSION: The proposed method of multimodal EEG-fMRI data is useful to increase spatial resolution of movement-related potentials and can also be applied to other kinds of event-related potentials.

Keywords: Brain Mapping/methods ; Cerebral Cortex/*physiology ; Cortical Synchronization ; Electrodes ; Electroencephalography/*methods ; Human ; Magnetoencephalography/*methods ; Motor Activity/*physiology ; Nerve Net ; Signal Processing, Computer-Assisted ; *Systems Integration

Previous simulation studies have stressed the importance of the use of fMRI priors in the estimation of cortical current density. However, no systematic variations of signal-to-noise ratio (SNR) and number of electrodes were explicitly taken into account in the estimation process. In this simulation study we considered the utility of including information as estimated from fMRI. This was done by using as the dependent variable both the correlation coefficient and the relative error between the imposed and the estimated waveforms at the level of cortical region of interests (ROI). A realistic head and cortical surface model was used. Factors used in the simulations were the different values of SNR of the scalp-generated data, the different inverse operators used to estimated the cortical source activity, the strengths of the fMRI priors in the fMRI-based inverse operators, and the number of scalp electrodes used in the analysis. Analysis of variance results suggested that all the considered factors significantly afflict the correlation and the relative error between the estimated and the simulated cortical activity. For the ROIs analyzed with simulated fMRI hot spots, it was observed that the best estimation of cortical source currents was performed with the inverse operators that used fMRI information. When the ROIs analyzed do not present fMRI hot spots, both standard (i.e., minimum norm) and fMRI-based inverse operators returned statistically equivalent correlation and relative error values.

Keywords: Analysis of Variance ; Brain Mapping ; Cerebral Cortex/*physiology ; *Computer Simulation ; *Electroencephalography ; Electrophysiology ; Human ; *Magnetic Resonance Imaging ; *Models, Neurological

The widely used technique of functional magnetic resonance imaging (fMRI) based on the blood oxygenation level-dependent (BOLD) effect is a tool for the investigation of changes in local brain activity upon stimulation. The principle of measurement is based on the assumption that there is a strong coupling between changes in neural activity, metabolism, vascular response and oxygen extraction in the area under investigation. As fMRI is on the way to become a routine tool in clinical examinations, we wanted to investigate whether, generally and under a variety of conditions, there is a strong link between the BOLD signal and neural activity. For clinical and experimental application of the method, it is crucial, whether the absence of changes in BOLD signal intensity upon stimulation can always be interpreted as an absence of changes in brain activity. We approached this question by inhibiting the nitric oxide mediated 'neurovascular coupling' via application of 7 nitroindazole. Before and after inhibition of this neurovascular coupling, we acquired evoked potentials and performed fMRI during somatosensory stimulation in rats. Cerebral blood flow response as well as BOLD signal intensity changes following electrical stimulation were abolished within 10 min after application of 7 nitroindazole, whereas somatosensory-evoked potentials were only slightly affected but still clearly detectable. Even 1 h after injection of 7 nitroindazole, there was still remaining electrical activity. Thus, we observed an uncoupling between electrical, i.e., neural activity and the BOLD signal. According to our results, the absence of BOLD signal changes did not permit the conclusion that there was no neural activity in the area under investigation. Our findings are especially relevant for the clinical application of fMRI in patients suffering from cerebrovascular and other brain diseases.

Keywords: Animals ; Blood Gas Analysis ; Cerebrovascular Circulation/*physiology ; Evoked Potentials, Somatosensory/physiology ; Heart Rate ; Magnetic Resonance Imaging ; Male ; Neurons/*physiology ; Oxygen/*blood ; Rats ; Rats, Sprague-Dawley ; Somatosensory Cortex/blood supply/physiology

Nowadays, several types of brain imaging device are available to provide images of the functional activity of the cerebral cortex 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 communicate with each other. In this study, advanced methods for the estimation of cortical connectivity from combined high-resolution electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data are presented. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual magnetic resonance images, multidipole source model, and regularized linear inverse source estimates of cortical current density. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) fMRI. We estimate functional cortical connectivity by computing the directed transfer function (DTF) on the estimated cortical current density waveforms in regions of interest (ROIs) on the modeled cortical mantle. The proposed method was able to unveil the direction of the information flow between the cortical regions of interest, as it is directional in nature. Furthermore, this method allows to detect changes in the time course of information flow between cortical regions in different frequency bands. The reliability of these techniques was further demonstrated by elaboration of high-resolution EEG and fMRI signals collected during visually triggered finger movements in four healthy subjects. Connectivity patterns estimated for this task reveal an involvement of right parietal and bilateral premotor and prefrontal cortical areas. This cortical region involvement resembles that revealed in previous studies where visually triggered finger movements were analyzed with the use of separate EEG or fMRI measurements.

Event-related fMRI (efMRI) has been repeatedly used to seek the neural sources of endogenous event-related potentials (ERP). However, significant discrepancies exist between the efMRI data and the results of previously published intracranial ERP studies of oddball task. To evaluate the capacity of efMRI to define the sources of the P3 component of ERP within the human brain, both efMRI and intracerebral ERP recordings were performed in eight patients with intractable epilepsy (five males and three females) during their preoperative invasive video-EEG monitoring. An identical auditory oddball task with frequent and target stimuli was completed in two sessions. A total of 606 intracerebral sites were electrophysiologically investigated by means of depth electrodes. In accordance with the finding of multiple intracerebral generators of P3 potential, the target stimuli evoked MRI signal increase in multiple brain regions. However, regions with evident hemodynamic and electrophysiological responses overlapped only partially. P3 generators were always found within hemodynamic-active sites, if these sites were investigated by means of depth electrodes. On the other hand, unequivocal local sources of P3 potential were apparently also located outside the regions with a significant hemodynamic response (typically in mesiotemporal regions). Both methods should thus be viewed as mutually complementary in investigations of the spatial distribution of cortical and subcortical activation during oddball task.

Keywords: Adult ; Auditory Cortex/physiology ; Auditory Perception/*physiology ; Brain/*physiology ; Cerebrovascular Circulation/physiology ; Electrodes, Implanted ; Electroencephalography ; Electrophysiology ; Epilepsy/physiopathology ; Evoked Potentials, Auditory/physiology ; Female ; Humans ; Magnetic Resonance Imaging ; Male ; Oxygen/blood ; Research Support, Non-U.S. Gov't

The linear transform model of functional magnetic resonance imaging (fMRI) hypothesizes that fMRI responses are proportional to local average neural activity averaged over a period of time. This work reports results from three empirical tests that support this hypothesis. First, fMRI responses in human primary visual cortex (V1) depend separably on stimulus timing and stimulus contrast. Second, responses to long-duration stimuli can be predicted from responses to shorter duration stimuli. Third, the noise in the fMRI data is independent of stimulus contrast and temporal period. Although these tests can not prove the correctness of the linear transform model, they might have been used to reject the model. Because the linear transform model is consistent with our data, we proceeded to estimate the temporal fMRI impulse-response function and the underlying (presumably neural) contrast-response function of human V1.

Keywords: Artifacts ; Human ; *Magnetic Resonance Imaging ; Models, Neurological ; Noise ; Photic Stimulation ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S. ; Time Factors ; Visual Cortex/*physiology

Keywords: ICP

The validity of functional magnetic resonance imaging (FMRI) brain maps with respect to the sites of neuronal activation is still unknown. One source of localization error may be pixels with large signal amplitudes, since such pixels may be expected to overlie large vessels, running remote from the centre of neuronal activation. In this study, magnetoencephalography was used to determine the centre of neuronal activation in a simple finger tapping task. The localization accuracy of conventional FMRI depending on FMRI signal enhancement was investigated relative to the magnetoencephalography reference. The results show a deterioration of FMRI localization with increasing signal amplitude related to increased contributions from large vessels. We conclude that FMRI data analysis should exclude large signal amplitudes and that magnetoencephalography may help to improve FMRI brain mapping results in a multimethod approach.

Keywords: Adult ; Brain/*physiology ; *Brain Mapping ; Human ; Magnetic Resonance Imaging/*methods ; *Magnetoencephalography ; Support, Non-U.S. Gov't

A general mathematical model for the delivery of O2 to the brain is presented, based on the assumptions that all of the brain capillaries are perfused at rest and that all of the oxygen extracted from the capillaries is metabolized. The model predicts that disproportionately large changes in blood flow are required in order to support small changes in the O2 metabolic rate. Interpreted in terms of this model, previous positron emission tomography (PET) studies of the human brain during neural stimulation demonstrating that cerebral blood flow (CBF) increases much more than the oxygen metabolic rate are consistent with tight coupling of flow and oxidative metabolism. The model provides a basis for the quantitative interpretation of functional magnetic resonance imaging (fMRI) studies in terms of changes in local CBF.

Keywords: Brain/physiology ; *Cerebrovascular Circulation ; Human ; *Models, Neurological ; *Oxygen Consumption ; *Regional Blood Flow ; Tomography, Emission-Computed

Recently, we introduced a new 'GLM-beamformer' technique for MEG analysis that enables accurate localisation of both phase-locked and non-phase-locked neuromagnetic effects, and their representation as statistical parametric maps (SPMs). This provides a useful framework for comparison of the full range of MEG responses with fMRI BOLD results. This paper reports a 'proof of principle' study using a simple visual paradigm (static checkerboard). The five subjects each underwent both MEG and fMRI paradigms. We demonstrate, for the first time, the presence of a sustained (DC) field in the visual cortex, and its co-localisation with the visual BOLD response. The GLM-beamformer analysis method is also used to investigate the main non-phase-locked oscillatory effects: an event-related desynchronisation (ERD) in the alpha band (8-13 Hz) and an event-related synchronisation (ERS) in the gamma band (55-70 Hz). We show, using SPMs and virtual electrode traces, the spatio-temporal covariance of these effects with the visual BOLD response. Comparisons between MEG and fMRI data sets generally focus on the relationship between the BOLD response and the transient evoked response. Here, we show that the stationary field and changes in oscillatory power are also important contributors to the BOLD response, and should be included in future studies on the relationship between neuronal activation and the haemodynamic response.

Keywords: Adult ; *Alpha Rhythm ; Brain Mapping ; Cerebrovascular Circulation ; Cortical Synchronization ; Female ; Humans ; Linear Models ; Magnetic Resonance Imaging/*methods ; Magnetoencephalography ; Male ; Oxygen/*blood ; Research Support, Non-U.S. Gov't ; Visual Cortex/*physiology

EEG-fMRI and EEG dipole source localisation are two non-invasive imaging methods that can be applied to the study of the haemodynamic and electrical consequences of epileptic discharges. Using them in combination has the potential to allow imaging with the spatial resolution of fMRI and the temporal resolution of EEG. However, although considerable data are available concerning their concordance in studies involving event-related potentials (ERPs), less is known about how well they agree in epilepsy. To this end, 17 patients were selected from a database of 57 who had undergone an EEG-fMRI scanning session followed by a separate EEG session outside of the scanner. Spatiotemporal dipole modelling was compared with the peak and closest EEG-fMRI activations and deactivations. On average, the dipoles were 58.5 mm from the voxel with the highest positive t value and 32.5 mm from the nearest activated voxel. For deactivations, the corresponding values were 60.8 and 34.0 mm. These values are considerably higher than is generally observed with ERPs, probably as a result of the relatively widespread field, which can lead to artificially deep dipoles, and the occurrence of EEG-fMRI responses remote from the presumed focus of the epileptic activity. The results suggest that EEG and MEG inverse solutions for equivalent current dipole approaches should not be strongly constrained by EEG-fMRI results in epilepsy, and that the use of distributed source modelling will be a more appropriate way of combining EEG-fMRI results with source localisation techniques.

Knowledge of regional cerebral hemodynamics has widespread application for both physiological research and clinical assessment because of the well-established interrelation between physiological function, energy metabolism, and localized blood supply. A magnetic resonance technique was developed for quantitative imaging of cerebral hemodynamics, allowing for measurement of regional cerebral blood volume during resting and activated cognitive states. This technique was used to generate the first functional magnetic resonance maps of human task activation, by using a visual stimulus paradigm. During photic stimulation, localized increases in blood volume (32 +/- 10 percent, n = 7 subjects) were detected in the primary visual cortex. Center-of-mass coordinates and linear extents of brain activation within the plane of the calcarine fissure are reported.

Keywords: Blood Volume ; *Brain Mapping ; Humans ; Magnetic Resonance Imaging/methods ; Magnetic Resonance Spectroscopy/methods ; Regional Blood Flow ; Research Support, Non-U.S. Gov't ; Research Support, U.S. Gov't, P.H.S. ; Visual Cortex/anatomy & histology/blood supply/*physiology

Although functional magnetic resonance imaging (fMRI) is a widely used and powerful tool for studying brain function, the quantitative interpretation of fMRI measurements for basic neuroscience and clinical studies can be complicated by inter-subject and inter-session variability arising from modulation of the baseline vascular state by disease, aging, diet, and pharmacological agents. In particular, recent studies have shown that the temporal dynamics of the cerebral blood flow (CBF) and the blood oxygenation level dependent (BOLD) responses to stimulus are modulated by changes in baseline CBF induced by various vasoactive agents and by decreases in vascular compliance associated with aging. These effects are not readily explained using current models of the CBF and BOLD responses. We present here a second-order nonlinear feedback model of the evoked CBF response in which neural activity modulates the compliance of arteriolar smooth muscle. Within this model framework, the baseline vascular state affects the dynamic response by changing the relative contributions of an active smooth muscle component and a passive connective tissue component to the overall vessel compliance. Baseline dependencies of the BOLD signal are studied by coupling the arteriolar compliance model with a previously described balloon model of the venous compartment. Numerical simulations show that the proposed model describes to first order the observed dependence of CBF and BOLD responses on the baseline vascular state.

Keywords: Aging/physiology ; Algorithms ; Arterioles/anatomy & histology/physiology ; Brain Chemistry/physiology ; Carbon Dioxide/physiology ; Cerebrovascular Circulation/*physiology ; Compliance ; Elasticity ; Hemoglobins/metabolism ; Humans ; Magnetic Resonance Imaging ; Models, Neurological ; Models, Statistical ; Muscle Contraction/physiology ; Muscle, Smooth, Vascular/anatomy & histology/*physiology ; Nonlinear Dynamics ; Oxygen/blood ; Research Support, Non-U.S. Gov't ; Viscosity

The acquisition of ERPs concurrently with fMRI in cognitive paradigms is appealing, but technically challenging. Little is known about the effects of the fMRI environment on the time-course and topography of previously documented ERP effects. We examined the replicability of ERP differences in the scanner at the level of individual subjects, using two cognitive paradigms and two statistical procedures. ERP P3 differences found outside the scanner in both paradigms were also robustly detected in the ERPs acquired during fMRI scanning. These P3 effects had equivalent time-courses and scalp topographies inside and outside the scanner. This replication at the level of individual data-sets has implications for the clinical applicability of ERP-fMRI and, more generally, for the quality of scanner recorded ERPs.

Alterations in gyral contour made it difficult to identify the motor cortex thought to be near an arteriovenous malformation (AVM) in a 24-year-old man considered for stereotactic radiosurgery. Functional imaging in three modalities was performed preoperatively to compare the reliability of localization using functional magnetic resonance imaging (fMRI) on a conventional scanner with positron emission tomography (PET) and magnetoencephalography (MEG). Similar tasks were used for each imaging modality in an attempt to activate and identify the sensory and motor cortex. Data from all three modalities converged for the sensory task, and fMRI and PET data converged for the motor task. The right hemisphere motor strip was localized adjacent and anterior to the AVM. These data were used in planning the radiosurgery isodose configuration to the AVM in order to reduce the irradiation of motor cortex parenchyma. A postoperative fMRI study was also performed using newer techniques to reduce head motion artifact and to improve signal-to-noise ratio. The data confirmed the conclusions derived from the preoperative evaluations. This study demonstrates how conventional MRI scanners can be used for functional studies of use in surgical planning.

Keywords: Adult ; Comparative Study ; Human ; Intracranial Arteriovenous Malformations/*pathology/radionuclide imaging/surgery ; *Magnetic Resonance Imaging ; *Magnetoencephalography ; Male ; Motor Cortex/*pathology/radionuclide imaging ; Radiosurgery ; Somatosensory Cortex/pathology/radionuclide imaging ; Stereotaxic Techniques ; *Tomography, Emission-Computed

An integrated model for magnetoencephalography (MEG) and functional Magnetic Resonance Imaging (fMRI) is proposed. In the model, the neural activity is related to the Post Synaptic Potentials (PSPs) which is common link between MEG and fMRI. Each PSP is modeled by the direction and strength of its current flow which are treated as random variables. The overall neural activity in each voxel is used for equivalent current dipole in MEG and as input of extended Balloon model in fMRI. The proposed model shows the possibility of detecting activation by fMRI in a voxel while the voxel is silent for MEG and vice versa. Parameters of the model can illustrate situations like closed field due to non-pyramidal cells, canceling effect of inhibitory PSP on excitatory PSP, and effect of synchronicity. In addition, the model shows that the crosstalk from neural activities of the adjacent voxels in fMRI may result in the detection of activations in these voxels that contain no neural activities. The proposed model is instrumental in evaluating and comparing different analysis methods of MEG and fMRI. It is also useful in characterizing the upcoming combined methods for simultaneous analysis of MEG and fMRI.

Keywords: Algorithms ; Excitatory Postsynaptic Potentials/physiology ; Humans ; Image Processing, Computer-Assisted/*statistics & numerical data ; Linear Models ; Magnetic Resonance Imaging/*statistics & numerical data ; Magnetoencephalography/*statistics & numerical data ; Models, Statistical ; Oxygen/blood

Artifacts generated by motion (e.g., ballistocardiac) of the head inside a high magnetic field corrupt recordings of EEG and EPs. This paper introduces a method for motion artifact cancellation. This method is based on adaptive filtering and takes advantage of piezoelectric motion sensor information to estimate the motion artifact noise. This filter estimates the mapping between motion sensor and EEG space, subtracting the motion-related noise from the raw EEG signal. Due to possible subject motion and changes in electrode impedance, a time-varying mapping of the motion versus EEG is required. We show that this filter is capable of removing both ballistocardiogram and gross motion artifacts, restoring EEG alpha waves (8-13 Hz), and visual evoked potentials (VEPs). This adaptive filter outperforms the simple band-pass filter for alpha detection because it is also capable of reducing noise within the frequency band of interest. In addition, this filter also removes the transient responses normally visible in the EEG window after echo planar image acquisition, observed during interleaved EEG/fMRI recordings. Our adaptive filter approach can be implemented in real-time to allow for continuous monitoring of EEG and fMRI during clinical and cognitive studies.

Keywords: Adult ; Alpha Rhythm ; *Artifacts ; Ballistocardiography ; Brain/*physiology ; *Electroencephalography ; *Evoked Potentials, Visual ; Female ; Human ; *Magnetic Resonance Imaging ; Male ; Motion ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S.

A component based method (CompCor) for the reduction of noise in both blood oxygenation level-dependent (BOLD) and perfusion-based functional magnetic resonance imaging (fMRI) data is presented. In the proposed method, significant principal components are derived from noise regions-of-interest (ROI) in which the time series data are unlikely to be modulated by neural activity. These components are then included as nuisance parameters within general linear models for BOLD and perfusion-based fMRI time series data. Two approaches for the determination of the noise ROI are considered. The first method uses high-resolution anatomical data to define a region of interest composed primarily of white matter and cerebrospinal fluid, while the second method defines a region based upon the temporal standard deviation of the time series data. With the application of CompCor, the temporal standard deviation of resting-state perfusion and BOLD data in gray matter regions was significantly reduced as compared to either no correction or the application of a previously described retrospective image based correction scheme (RETROICOR). For both functional perfusion and BOLD data, the application of CompCor significantly increased the number of activated voxels as compared to no correction. In addition, for functional BOLD data, there were significantly more activated voxels detected with CompCor as compared to RETROICOR. In comparison to RETROICOR, CompCor has the advantage of not requiring external monitoring of physiological fluctuations.

We used a real-skull phantom head to investigate the performances of representative methods for EEG source localization when considering various head models. We describe several experiments using a montage with current sources located at multiple positions and orientations inside a human skull filled with a conductive medium. The robustness of selected methods based on distributed source models is evaluated as various solutions to the forward problem (from the sphere to the finite element method) are considered. Experimental results indicate that inverse methods using appropriate cortex-based source models are almost always able to locate the active source with excellent precision, with little or no spurious activity in close or distant regions, even when two sources are simultaneously active. Superior regularization schemes for solving the inverse problem can dramatically help the estimation of sparse and focal active zones, despite significant approximation of the head geometry and the conductivity properties of the head tissues. Realistic head models are necessary, though, to fit the data with a reasonable level of residual variance.

Keywords: Electroencephalography/*methods ; Head/*radiation effects ; Human ; Models, Theoretical ; Phantoms, Imaging ; Reproducibility of Results ; Skull/*radiation effects ; Time Factors

We introduce a bottom-up model for integrating electroencephalography (EEG) or magnetoencephalography (MEG) with functional magnetic resonance imaging (fMRI). An extended neural mass model is proposed based on the physiological principles of cortical minicolumns and their connections. The fMRI signal is extracted from the proposed neural mass model by introducing a relationship between the stimulus and the neural activity and using the resultant neural activity as input of the extended Balloon model. The proposed model, validated using simulations, is instrumental in evaluating the upcoming combined methods for simultaneous analysis of MEG/EEG and fMRI.

Keywords: Algorithms ; Brain/*physiology ; Brain Mapping/*methods ; Computer Simulation ; Diagnosis, Computer-Assisted/*methods ; Electroencephalography/*methods ; Evoked Potentials/physiology ; Humans ; Magnetic Resonance Imaging/*methods ; Magnetoencephalography/*methods ; *Models, Neurological ; Nerve Net/physiology ; Systems Integration

Combined analysis of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) has the potential to provide higher spatiotemporal resolution than either method alone. In some situations, in which the activity of interest cannot be reliably reproduced (e.g., epilepsy, learning, sleep states), accurate combined analysis requires simultaneous acquisition of EEG and fMRI. Simultaneous measurements ensure that the EEG and fMRI recordings reflect the exact same brain activity state. We took advantage of the spatial filtering properties of the bipolar montage to allow recording of very short (125-250 ms) visual-evoked potentials (VEPs) during fMRI. These EEG and fMRI measurements are of sufficient quality to allow source localization of the cortical generators. In addition, our source localization approach provides a combined EEG/fMRI analysis that does not require any manual selection of fMRI activations or placement of source dipoles. The source of the VEP was found to be located in the occipital cortex. Separate analysis of EEG and fMRI data demonstrated good spatial overlap of the observed activated sites. As expected, the combined EEG/fMRI analysis provided better spatiotemporal resolution than either approach alone. The resulting spatiotemporal movie allows for the millisecond-to-millisecond display of changes in cortical activity caused by visual stimulation. These data reveal two peaks in activity corresponding to the N75 and the P100 components. This type of simultaneous acquisition and analysis allows for the accurate characterization of the location and timing of neurophysiological activity in the human brain.

Keywords: Adult ; *Brain Mapping ; Computer Graphics ; Data Display ; Dominance, Cerebral/physiology ; *Electroencephalography ; Evoked Potentials, Visual/*physiology ; Female ; Human ; *Image Enhancement ; *Image Processing, Computer-Assisted ; Imaging, Three-Dimensional ; *Magnetic Resonance Imaging ; Male ; Occipital Lobe/*physiology ; Photic Stimulation ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S.

An issue in blood oxygenation level dependent contrast-based functional MRI is the accurate interpretation of the activation-induced signal changes. Hemodynamic factors other than activation-induced changes in blood oxygenation are known to contribute to the signal change magnitudes and dynamics, and therefore need to be accounted for or removed. In this paper, a general method for removal of effects other than activation-induced blood oxygenation changes from fMRI brain activation maps by the use of hypercapnic stress normalization is introduced. First, the effects of resting blood volume distribution across voxels on activation-induced BOLD-based fMRI signal changes are shown to be significant. Second, the effects of hypercapnia and hypoxia on resting and activation-induced signal changes are demonstrated. These results suggest that global hemodynamic stresses may be useful for non-invasive mapping of blood volume. Third, the normalization technique is demonstrated.

Keywords: Brain/*anatomy & histology/blood supply/*physiology ; Brain Mapping/*methods ; Carbon Dioxide/*blood ; Humans ; Image Processing, Computer-Assisted/methods ; Magnetic Resonance Imaging/*methods ; Oxygen/blood

Integrating electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data may help to optimize anatomical and temporal resolution in the investigation of cortical function. Successful removal of fMRI scanning artifacts from continuous EEG in simultaneous recordings has been reported. We assessed the feasibility of recording reliable visual evoked potentials (VEPs) during fMRI scanning using available artifact removing procedures. EEG during administration of visual stimuli was recorded using MRI-compatible 32-channel equipment in nine normal subjects (mean age, 23.9 +/- 2.5 years), with and without fMRI acquisition. fMRI scanning and cardioballistographic artifacts were removed after subtraction of averaged artifact waveforms. Consistency between VEPs waveforms and of P1 and N1 peak latencies and amplitudes in the two conditions was assessed. Good correlation was found between VEP waveforms (Pearson's correlation coefficient: r(P) between 0.76-0.94 across subjects; P < 0.0001) and between latency or amplitude of P1 and N1 peaks (latencies: r = 0.7, P < 0.035; amplitudes: r > 0.65, P < 0.05; Spearman rank correlation coefficient) in the two recording conditions. No significant differences were found between P1 and N1 parameters in the two conditions (Wilcoxon signed rank test). Consistent VEP waveforms, latencies, and amplitudes with and without fMRI scanning indicate that reliable VEPs may be obtained simultaneously with fMRI recording. This possibility might be helpful by shortening recording times and reducing variability from learning, habituation, and fatigue phenomena from separate recordings for the integration of event-related EEG and fMRI data.

Keywords: Adult ; *Artifacts ; Brain/*physiology ; *Brain Mapping ; Electroencephalography ; Evoked Potentials, Visual/*physiology ; Female ; Humans ; *Magnetic Resonance Imaging ; Male ; Photic Stimulation ; Research Support, Non-U.S. Gov't

The coupling between neurotransmitter-induced changes in neuronal activity and the resultant hemodynamic response is central to the interpretation of neuroimaging techniques. In the present study, MRI experiments showed that dopamine transporter blockers such as cocaine and dopamine releasers such as amphetamine and D1 receptor agonists induced large positive increases in relative cerebral blood volume (rCBV) that were not sensitive to nitric oxide synthase inhibition. However, D1/D5 receptor antagonism with SCH-23390 prevented or blocked the hemodynamic response without any concomitant effect on dopamine release. Dopamine D2/D3 receptor agonists, in contrast, induced negative changes in rCBV in brain regions corresponding largely to those endowed with these receptors. D1 and D5 receptor mRNAs were expressed in microvessels of responsive brain areas, while D2 and D3 receptors were not consistently associated with the microvascular bed. D3 receptors had an astroglial localization. Together, these experiments show that direct effects of dopamine upon the vasculature cannot be ignored in measuring the hemodynamic coupling associated with dopaminergic drugs. These results further suggest that this coupling is partially mediated through D1/D5 receptors on the microvasculature leading to increased rCBV and through astroglial D3 receptors leading to decreased rCBV. These data provide additional support for the role of local post-synaptic events in neurovascular coupling and emphasize that the interpretation of fMRI signals exclusively in terms of neuronal activity may be incomplete.

Functional neuroimaging in humans is used widely to study brain function in relation to human disease and cognition. The neural basis of neuroimaging signals is probably synaptic activity, but the effect of context, defined as the interaction between synaptic inhibition, excitation, and the electroresponsive properties of the targeted neurons, is not well understood. We examined here the effect of interaction of synaptic excitation and net inhibition on the relationship between electrical activity and vascular signals in the cerebellar cortex. We show that stimulation of the net inhibitory parallel fibers simultaneously with stimulation of the excitatory climbing fibers leads to a further rise in total local field potentials (LFP) and cerebral blood flow (CBF) amplitudes, not a decrease, as predicted from theoretical studies. However, the combined stimulation of the parallel and climbing fiber systems produced changes in CBF and LFP that were smaller than their algebraic sum evoked by separate stimulation of either system. This finding was independent of the starting condition, i.e., whether inhibition was superimposed on a state of excitation or vice versa. The attenuation of the increases in LFP and CBF amplitudes was similar, suggesting that synaptic activity and CBF were coupled under these conditions. The result might be explained by a relative neuronal refractoriness that relates to the intrinsic membrane properties of Purkinje cells, which determine the recovery time of these cells. Our work implies that neuronal and vascular signals are context-sensitive and that their amplitudes are modulated by the electroresponsive properties of the targeted neurons.

Keywords: Animals ; Blood Flow Velocity/*physiology ; Brain/blood supply/*physiology ; Cerebrovascular Circulation/*physiology ; Electric Stimulation ; Humans ; Laser-Doppler Flowmetry ; Male ; Membrane Potentials/physiology ; Microelectrodes ; Olivary Nucleus/blood supply/physiology ; Purkinje Cells/physiology ; Rats ; Rats, Wistar ; Research Support, Non-U.S. Gov't ; Synapses/physiology

This paper deals with the estimation of the blood oxygen level-dependent response to a stimulus, as measured in functional magnetic resonance imaging (fMRI) data. A precise estimation is essential for a better understanding of cerebral activations. The most recent works have used a nonparametric framework for this estimation, considering each brain region as a system characterized by its impulse response, the so-called hemodynamic response function (HRF). However, the use of these techniques has remained limited since they are not well-adapted to real fMRI data. Here, we develop a threefold extension to previous works. We consider asynchronous event-related paradigms, account for different trial types and integrate several fMRI sessions into the estimation. These generalizations are simultaneously addressed through a badly conditioned observation model. Bayesian formalism is used to model temporal prior information of the underlying physiological process of the brain hemodynamic response. By this way, the HRF estimate results from a tradeoff between information brought by the data and by our prior knowledge. This tradeoff is modeled with hyperparameters that are set to the maximum-likelihood estimate using an expectation conditional maximization algorithm. The proposed unsupervised approach is validated on both synthetic and real fMRI data, the latter originating from a speech perception experiment.

Keywords: *Algorithms ; Brain/*blood supply/*physiology ; Brain Mapping/*methods ; Cerebrovascular Circulation/physiology ; Comparative Study ; Computer Simulation ; Hemodynamic Processes/physiology ; Human ; Image Interpretation, Computer-Assisted/*methods ; Imaging, Three-Dimensional/*methods ; Likelihood Functions ; Magnetic Resonance Imaging/*methods ; *Models, Cardiovascular ; Models, Statistical ; Reproducibility of Results ; Sensitivity and Specificity ; Speech Perception/physiology ; Support, Non-U.S. Gov't

All perfusion-based imaging modalities depend on the relationship between neuronal and vascular activity. However, the relationship between stimulus and response was never fully characterized. With the use of optical imaging (intrinsic signals and intravascular fluorescent dyes) during repetitive stimulation paradigms, we observed reduced responses with temporally close stimuli. Cortical evoked potentials, however, did not produce the same reduced responsiveness. We therefore termed these intervals of reduced responsiveness refractory periods. During these refractory periods an ability to respond was retained, but at a near 60% reduction in the initial magnitude. Although increasing the initial stimulus duration lengthened the observed refractory periods, significantly novel or temporally spaced stimuli overcame them. We observed this phenomenon in both rodent and human subjects in somatosensory and auditory cortices. These results have significant implications for understanding the capacities, mechanisms, and distributions of neurovascular coupling and thereby possess relevance to all perfusion-dependent functional imaging techniques.

Keywords: Acoustic Stimulation ; Animals ; Evoked Potentials, Somatosensory/physiology ; Fluorescent Dyes ; Humans ; Image Processing, Computer-Assisted ; Male ; Optics ; Proprioception/physiology ; Rats ; Rats, Sprague-Dawley ; Refractory Period, Neurologic/*physiology ; Research Support, U.S. Gov't, P.H.S. ; Somatosensory Cortex/*physiology ; Temporal Lobe/physiology

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) measurements were performed in six healthy subjects to determine the representation of stimulation of the right thumb in somatosensory cortex. In all subjects EEG-based dipole locations could be determined in primary (S1) and secondary (S2) somatosensory cortex. The stimulus-induced blood oxygenation level dependent response of the fMRI showed deviations of 23.5 mm (standard deviation, SD = 6.9) for S1 and 15.7 mm (SD = 3.5) for S2 cortex. fMRI constrained dipole searches lead to higher residual variances. The data show that simultaneous EEG and fMRI measurements of somatosensory activity are feasible and yield reliable and valid results.

Keywords: Adult ; Analysis of Variance ; Electric Stimulation/methods ; Electroencephalography/*methods/statistics & numerical data ; Female ; Humans ; Linear Models ; Magnetic Resonance Imaging/*methods/statistics & numerical data ; Male ; Research Support, Non-U.S. Gov't ; Somatosensory Cortex/*physiology

The blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal response to neural stimulation is influenced by many factors that are unrelated to the stimulus. These factors are physiological, such as the resting venous cerebral blood volume (CBV(v)) and vessel size, as well as experimental, such as pulse sequence and static magnetic field strength (B(0)). Thus, it is difficult to compare task-induced fMRI signals across subjects, field strengths, and pulse sequences. This problem can be overcome by normalizing the neural activity-induced BOLD fMRI response by a global hypercapnia-induced BOLD signal. To demonstrate the effectiveness of the BOLD normalization approach, gradient-echo BOLD fMRI at 1.5, 4, and 7 T and spin-echo BOLD fMRI at 4 T were performed in human subjects. For neural stimulation, subjects performed sequential finger movements at 2 Hz, while for global stimulation, subjects breathed a 5 CO(2) gas mixture. Under all conditions, voxels containing primarily large veins and those containing primarily active tissue (i.e., capillaries and small veins) showed distinguishable behavior after hypercapnic normalization. This allowed functional activity to be more accurately localized and quantified based on changes in venous blood oxygenation alone. The normalized BOLD signal induced by the motor task was consistent across different magnetic fields and pulse sequences, and corresponded well with cerebral blood flow measurements. Our data suggest that the hypercapnic normalization approach can improve the spatial specificity and interpretation of BOLD signals, allowing comparison of BOLD signals across subjects, field strengths, and pulse sequences. A theoretical framework for this method is provided.

Keywords: Adult ; Algorithms ; Brain Mapping ; Cerebrovascular Circulation ; Comparative Study ; Echo-Planar Imaging ; Female ; Humans ; Hypercapnia/*blood ; Image Interpretation, Computer-Assisted ; Magnetic Resonance Imaging/*methods ; Male ; Oxygen/*blood ; Psychomotor Performance/physiology ; Reference Values ; Research Support, Non-U.S. Gov't ; Research Support, U.S. Gov't, P.H.S.

We mapped ocular dominance columns (ODCs) in normal human subjects using high-field (4 T) functional magnetic resonance imaging (fMRI) with a segmented echo planar imaging technique and an in-plane resolution of 0.47 x 0.47 mm(2). The differential responses to left or right eye stimulation could be reliably resolved in anatomically well-defined sections of V1. The orientation and width ( approximately 1 mm) of mapped ODC stripes conformed to those previously revealed in postmortem brains stained with cytochrome oxidase. In addition, we showed that mapped ODC patterns could be largely reproduced in different experiments conducted within the same experimental session or over different sessions. Our results demonstrate that high-field fMRI can be used for studying the functions of human brains at columnar spatial resolution.

Keywords: Adult ; *Brain Mapping ; Dominance, Cerebral/*physiology ; Humans ; *Magnetic Resonance Imaging ; Male ; *Ocular Physiology ; Photic Stimulation ; Research Support, Non-U.S. Gov't ; Visual Cortex/physiology

Using a model of the functional MRI (fMRI) impulse response based on published data, we have demonstrated that the form of the fMRI response to stimuli of freely varied timing can be modeled well by convolution of the impulse response with the behavioral stimulus. The amplitudes of the responses as a function of parametrically varied behavioral conditions are fitted well using a piecewise linear approximation. Use of the combined model, in conjunction with correlation analysis, results in an increase in sensitivity for the MRI study. This approach, based on the well-established methods of linear systems analysis, also allows a quantitative comparison of the response amplitudes across subjects to a broad range of behavioral conditions. Fit parameters, derived from the amplitude data, are relatively insensitive to a variety of MRI-related artifacts and yield results that are compared readily across subjects.

Keywords: Adult ; Brain/anatomy & histology/*physiology ; Brain Mapping ; Cerebrovascular Circulation/physiology ; Human ; Linear Models ; Magnetic Resonance Imaging/*statistics & numerical data ; Photic Stimulation ; Psychomotor Performance/physiology

Recent advances in brain imaging techniques, including functional magnetic resonance imaging (fMRI), offer great promise for noninvasive mapping of brain function. However, the indirect nature of the imaging signals to the underlying neural activity limits the interpretation of the resulting maps. The present report represents the first systematic study with sufficient statistical power to quantitatively characterize the relationship between changes in blood oxygen content and the neural spiking and synaptic activity. Using two-dimensional optical measurements of hemodynamic signals, simultaneous recordings of neural activity, and an event-related stimulus paradigm, we demonstrate that (1) there is a strongly nonlinear relationship between electrophysiological measures of neuronal activity and the hemodynamic response, (2) the hemodynamic response continues to grow beyond the saturation of electrical activity, and (3) the initial increase in deoxyhemoglobin that precedes an increase in blood volume is counterbalanced by an equal initial decrease in oxyhemoglobin.

Keywords: Animals ; Brain Mapping ; Comparative Study ; Computer Simulation ; Demography ; Electric Stimulation ; Electrophysiology/methods ; Evoked Potentials, Somatosensory/physiology ; Hemodynamic Processes/physiology ; Hemoglobins/*metabolism ; Magnetic Resonance Imaging/methods ; Neurons/*physiology ; Nonlinear Dynamics ; Oxygen/*metabolism ; Rats ; Somatosensory Cortex/blood supply/cytology/*metabolism ; Spectrum Analysis/methods ; Support, Non-U.S. Gov't ; Support, U.S. Gov't, P.H.S. ; Time Factors

To assess the degree of fine-scale somatotopy within the hand area of the human primary motor cortex (M1), functional mapping of individual movements of all fingers was performed in healthy young subjects (n = 7) using MRI at 0.8 x 0.8 mm2 resolution and 4 mm section thickness. The experimental design comprised both a direct paradigm contrasting single digit movements vs. motor rest and multiple differential paradigms contrasting single digit movements vs. the movement of another digit. Direct mapping resulted in largely overlapping activations. A somatotopic arrangement was only recognizable when considering the mean center-of-mass coordinates of individual digit representations averaged across subjects. In contrast, differential paradigms revealed more segregated and somatotopically ordered activations in single subjects. The use of center-of-mass coordinates yielded inter-digit distances ranging from 2.0 to 16.8 mm, which reached statistical significance for pairs of more distant digits. For the middle fingers, the functional somatotopy obtained by differential mapping was dependent on the choice of the digit used for control. These results confirm previous concepts that finger somatotopy in the human M1 hand area emerges as a functional predominance of individual digit representations sharing common areas in a distributed though ordered network.

Keywords: Adult ; Analysis of Variance ; Brain Mapping/*methods ; Female ; Fingers/*physiology ; Human ; Least-Squares Analysis ; Male ; Motor Cortex/*physiology

There is an increasing interest in using physiologically plausible models in fMRI analysis. These models do raise new mathematical problems in terms of parameter estimation and interpretation of the measured data. In this paper, we show how to use physiological models to map and analyze brain activity from fMRI data. We describe a maximum likelihood parameter estimation algorithm and a statistical test that allow the following two actions: selecting the most statistically significant hemodynamic model for the measured data and deriving activation maps based on such model. Furthermore, as parameter estimation may leave much incertitude on the exact values of parameters, model identifiability characterization is a particular focus of our work. We applied these methods to different variations of the Balloon Model (Buxton, R.B., Wang, E.C., and Frank, L.R. 1998. Dynamics of blood flow and oxygenation changes during brain activation: the balloon model. Magn. Reson. Med. 39: 855-864; Buxton, R.B., Uludag, K., Dubowitz, D.J., and Liu, T.T. 2004. Modelling the hemodynamic response to brain activation. NeuroImage 23: 220-233; Friston, K. J., Mechelli, A., Turner, R., and Price, C. J. 2000. Nonlinear responses in fMRI: the balloon model, volterra kernels, and other hemodynamics. NeuroImage 12: 466-477) in a visual perception checkerboard experiment. Our model selection proved that hemodynamic models better explain the BOLD response than linear convolution, in particular because they are able to capture some features like poststimulus undershoot or nonlinear effects. On the other hand, nonlinear and linear models are comparable when signals get noisier, which explains that activation maps obtained in both frameworks are comparable. The tools we have developed prove that statistical inference methods used in the framework of the General Linear Model might be generalized to nonlinear models.

Several properties of the cerebral cortex, including its columnar and laminar organization, as well as the topographic organization of cortical areas, can only be properly understood in the context of the intrinsic two-dimensional structure of the cortical surface. In order to study such cortical properties in humans, it is necessary to obtain an accurate and explicit representation of the cortical surface in individual subjects. Here we describe a set of automated procedures for obtaining accurate reconstructions of the cortical surface, which have been applied to data from more than 100 subjects, requiring little or no manual intervention. Automated routines for unfolding and flattening the cortical surface are described in a companion paper. These procedures allow for the routine use of cortical surface-based analysis and visualization methods in functional brain imaging.

Keywords: Brain Mapping/instrumentation ; Cerebral Cortex/*anatomy & histology ; Human ; Image Processing, Computer-Assisted/*instrumentation ; Magnetic Resonance Imaging/*instrumentation ; Reference Values ; Software

We studied MEG and fMRI responses to electric median and tibial nerve stimulation in five healthy volunteers. The aim was to compare the results with those of a previous study using only fMRI on the primary and secondary somatosensory cortices in which the somatotopic organization of SII was observed wit