Functional magnetic resonance imaging of awake monkeys: some approaches for improving imaging quality

Functional magnetic resonance imaging (fMRI) at high magnetic field strength can suffer from serious degradation of image quality because of motion and physiological noise, as well as spatial distortions and signal losses due to susceptibility effects. Overcoming such limitations is essential for sensitive detection and reliable interpretation of fMRI data. These issues are particularly problematic in studies of awake animals. As part of our initial efforts to study functional brain activations in awake, behaving monkeys using fMRI at 4.7 T, we have developed acquisition and analysis procedures to improve image quality with encouraging results.We evaluated the influence of two main variables on image quality. First, we show how important the level of behavioral training is for obtaining good data stability and high temporal signal-to-noise ratios. In initial sessions, our typical scan session lasted 1.5 h, partitioned into short (<10 min) runs. During reward periods and breaks between runs, the monkey exhibited movements resulting in considerable image misregistrations. After a few months of extensive behavioral training, we were able to increase the length of individual runs and the total length of each session. The monkey learned to wait until the end of a block for fluid reward, resulting in longer periods of continuous acquisition. Each additional 60 training sessions extended the duration of each session by 60 min, culminating, after about 140 training sessions, in sessions that last about 4 h. As a result, the average translational movement decreased from over 500 μm to less than 80 μm, a displacement close to that observed in anesthetized monkeys scanned in a 7-T horizontal scanner.Another major source of distortion at high fields arises from susceptibility variations. To reduce such artifacts, we used segmented gradient-echo echo-planar imaging (EPI) sequences. Increasing the number of segments significantly decreased susceptibility artifacts and image distortion. Comparisons of images from functional runs using four segments with those using a single-shot EPI sequence revealed a roughly twofold improvement in functional signal-to-noise-ratio and 50% decrease in distortion. These methods enabled reliable detection of neural activation and permitted blood-oxygenation-level-dependent-based mapping of early visual areas in monkeys using a volume coil.In summary, both extensive behavioral training of monkeys and application of segmented gradient-echo EPI sequence improved signal-to-noise ratio and image quality. Understanding the effects these factors have is important for the application of high field imaging methods to the detection of submillimeter functional structures in the awake monkey brain.     

Functional phantom for fMRI: a feasibility study

Functional magnetic resonance imaging (fMRI) reveals changes in blood oxygen level-dependent (BOLD) signal after considerable processing. This paper describes the implementation and testing of an fMRI phantom where electric current applied to a thin wire within a proton-rich medium substituted BOLD distortion of the magnetic field; the scanner detects these two distortions as practically identical signal changes. The magnitude of the change depended on the current strength. The phantom has a number of possible applications. Signal changes across sessions, days, instruments and individuals could be monitored. Placing the phantom close to a subject during an fMRI experiment could allow differentiating sensitivity changes in the scanner due to instrumentation from changes in the subject's state and performance during the experiment. The spatial extent of brain activations and effects of various changes in the chain of image formation could be analyzed using current-induced "activations". Furthermore, the phantom could expedite fMRI sequence development by reducing the need to scan human subjects, who introduce uncertainty to the signal. Thus, this fMRI phantom could be useful for both cognitive fMRI studies and scanner calibration.     

增强低场下脑功能磁共振成像显著性的研究

实现了基于低场0.35 T磁共振成像系统的大脑功能磁共振成像（functional Magnetic Resonance Imaging,fMRI）的研究. 基于质子密度加权的快速自旋回波(Turbo Spin Echo,TSE)图像,重点研究增强低场fMRI显著性的方法,目的在于提高低场fMRI的可用性. 结果表明：健康受试者在执行手动任务期间,大脑运动区的信号强度变化可以由基于血管外质子信号增强 (Signal Enhancement by Extravascular water Protons,EEP)的对比机制探测. 优化TSE序列参数能提高图像SNR和扫描速度,并在统计分析中增加外在屏蔽图像,可以有效地提高低场下fMRI研究结果的显著性.     

The effect of artifacts on dependence measurement in fMRI

The study of effective connectivity by means of neuroimaging depends on the measurement of similarity between activity patterns at different locations in the brain, without necessarily presupposing a particular model for this dependence. When these interactions are measured using functional magnetic resonance imaging (fMRI) techniques, however, imaging and physiological artifacts create patterns of dependence that may be unrelated to cortical activity. We demonstrate some of these effects through the measurement of short-range dependencies present in fMRI scans of the primary visual cortex (V1) in the anaesthetized macaque monkey. High-field (4.7 T) fMRI scans were conducted to measure responses based on the blood oxygen level-dependent contrast mechanism, during periods of no sensory stimulation and of visual stimulation with rotating polar-transformed checkerboard gratings. Dependence between the haemodynamic activity at different spatial locations (i.e., different voxels) was measured using correlation, mutual information and functional covariance. Particular attention was paid to understanding the sources of spurious dependence that may be observed during such investigations. Two main effects were detected: (a) short-range correlations introduced by the process of image reconstruction and (b) perturbations in the haemodynamic response caused by breathing. The image reconstruction artifacts were shown to create an artificially high short-range dependence in the readout direction of the scan, and the breathing artifacts caused enhanced short-range dependence in both the readout and phase-encode directions. Additional dependence in the phase-encode direction due to image-ghosting is also possible but will not be discussed in this report, as it can be alleviated by fine adjustment of preemphasis (elimination of eddy currents). A technique is described for removing breathing artifacts, and the effect of breathing on the apparent dependence between voxels is illustrated. The correlation of haemodynamic activity with the stimulus was found to be affected by breathing, although this effect can be neutralised by averaging the haemodynamic responses over many repetitions of the stimulus. Nonetheless, patterns of dependent activity between voxels may be lost in this averaging process, which makes the removal of breathing artifacts necessary if statistical dependence and the study of effective connectivity is the primary aim of an investigation.     

A two-compartment gel phantom for optimization and quality assurance in clinical BOLD fMRI

Clinical applications of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) depend heavily on robust paradigms, imaging methods and analysis procedures. In this work, as a means to optimize and perform quality assurance of the entire imaging and analysis chain, a phantom that provides a well known and reproducible signal change similar to a block type fMRI experiment is presented. It consists of two gel compartments with slightly different T2 that dynamically enter and leave the imaged volume. The homogeneous gel in combination with a cylindrical geometry results in a well-defined T*2 difference causing a signal difference between the two compartments in T*2-weighted MR images. From time series data obtained with the phantom, maps of percent signal change (PSC) and t-values are calculated. As an example of image parameter optimisation, the phantom is demonstrated to be useful for accurate determination of the influence of echo time (TE) on BOLD fMRI results, taking the t-value as a measure of sensitivity. In addition, the phantom is proposed as a tool for quality assurance (QA) since reproducible time series and t-maps are obtained in a series of independent repeat experiments. The phantom is relatively simple to build and can therefore be used by any clinical fMRI center.     

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%.     

Accurate and sensitive measurements of magnetic susceptibility using echo planar imaging

Susceptibility differences are common causes for artifacts in magnetic resonance (MR); therefore, it is important to choose phantom materials in a way that these artifacts are kept at a minimum. In this study, a previously proposed MR imaging (MRI) method [Beuf O, Briguet A, Lissac M, Davis R. Magnetic resonance imaging for the determination of magnetic susceptibility of materials. J Magn Reson 1996; Series B(112):111-118] was improved to facilitate sensitive in-house measurements of different phantom materials so that such artifacts can more easily be minimized. Using standard MRI protocols and distilled water as reference, we measured magnetic volume susceptibility differences with a clinical MR system. Two imaging techniques, echo planar imaging (EPI) and spin echo, were compared using liquid samples whose susceptibilities were verified by MR spectroscopy. The EPI sequence has a very narrow bandwidth in the phase-encoding direction, which gives an increased sensitivity to magnetic field inhomogeneities. All MRI measurements were evaluated in two ways: (1) manual image analysis and (2) model fitting. The narrow bandwidth of the EPI made it possible to detect very small susceptibility differences (equivalent susceptibility difference, Deltachi(e)> or =0.02 ppm), and even plastics could be measured. Model fitting yielded high accuracy and high sensitivity and was less sensitive to other image artifacts as compared with manual image analysis.     

Detection of pseudoperiodic patterns using partial acquisition of magnetic resonance images

Improving the resolution of magnetic resonance imaging (MRI), or, alternatively, reducing the acquisition time, can be quite beneficial for many applications. The main motivation of this work is the assumption that any information that is a priori available on the target image could be used to achieve this goal. In order to demonstrate this approach, we present a novel partial acquisition strategy and reconstruction algorithm, suitable for the special case of detection of pseudoperiodic patterns. Pseudoperiodic patterns are frequently encountered in the cerebral cortex due to its columnar functional organization (best exemplified by orientation columns and ocular dominance columns of the visual cortex). We present a new MRI research methodology, in which we seek an activity pattern, and a pattern-specific experiment is devised to detect it. Such specialized experiments extend the limits of conventional MRI experiments by substantially reducing the scan time. Using the fact that pseudoperiodic patterns are localized in the Fourier domain, we present an optimality criterion for partial acquisition of the MR signal and a strategy for obtaining the optimal discrete Fourier transform (DFT) coefficients. A by-product of this strategy is an optimal linear extrapolation estimate. We also present a nonlinear spectral extrapolation algorithm, based on projections onto convex sets (POCSs), used to perform the actual reconstruction. The proposed strategy was tested and analyzed on simulated signals and in MRI phantom experiments.     

Performance of a novel piezoelectric motor at 4.7 T: applications and initial tests

The focus of this report was to test the performance of a novel piezoelectric motor under high magnetic field strength conditions and to investigate its potential applications in small animal magnetic resonance imaging (MRI). The device is made entirely of nonferrous materials and consists of four piezoelectric ceramic plates connected to a threaded metal tube through which a screw migrates. Ultrasonic vibrations of the threads inherent to the tube result in rotational and translational motion of the screw. Potential applications of the piezoelectric motor were investigated at 4.7 T. Firstly, phantom studies showed the motor was capable of accurately delivering low injection volumes ( approximately 0.01 ml). Dynamic contrast-enhanced MRI (DCE-MRI) studies performed in vivo using serially acquired T1-weighted, spin-echo imaging demonstrated the ability of the motor to reliably administer MR contrast-enhancing agent into live tumor-bearing mice without the introduction of image artifacts. In a second set of experiments, the motor allowed for controlled, dynamic repositioning of an anatomic slice of interest in a live animal to magnetic field isocenter, which resulted in reduced geometric distortion and image artifact due to improved radiofrequency and gradient field homogeneity. In conclusion, piezoelectric motors are MR compatible and offer great potential for improving MRI efficiency and throughput, particularly in a preclinical setting. Further investigation into applications such as automated capacitor tuning and impedance matching for MR transceiver coils is warranted.     

Evaluating gradient artifact correction of EEG data acquired simultaneously with fMRI

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) has become a widely used application in spite of EEG perturbations due to electromagnetic interference in the MR environment. The most prominent and disturbing artifacts in the EEG are caused by the alternating magnetic fields (gradients) of the MR scanner. Different methods for gradient artifact correction have been developed. Here we propose an approach for the systematic evaluation and comparison of these gradient artifact correction methods. Exemplarily, we evaluate different algorithms all based on artifact template subtraction--the currently most established means of gradient artifact removal. We introduce indices for the degree of gradient artifact reduction and physiological signal preservation. The combination of both indices was used as a measure for the overall performance of gradient artifact removal and was shown to be useful in identifying problems during artifact removal. We demonstrate that the evaluation as proposed here allows to reveal frequency-band specific performance differences among the algorithms. This emphasizes the importance of carefully selecting the artifact correction method appropriate for the respective case.     

Functional magnetic resonance imaging reference phantom

Functional magnetic resonance imaging (fMRI) is widely used to pinpoint active brain areas. Changes in neuronal activity modulate the local blood oxygenation level, and the associated modulation of the magnetic field homogeneity can be detected with magnetic resonance imaging. Thus, the blood oxygenation level-dependent (BOLD) fMRI indirectly measures neuronal activity. Similar modulation of magnetic field homogeneity was here elicited by other means to generate a BOLD-like change in a new phantom constructed to provide reference activations during fMRI. Magnetic inhomogeneities were produced by applying current to coils located near the phantom containing 1.5 ml of Gd-doped water. The signal-to-noise ratio of the images, produced by gradient-recalled echo-planar imaging, varied between 104 and 107 at a selected voxel when the field was and was not inhomogenized, respectively. The contrast of signals between homogeneous and inhomogeneous conditions was generally stable, except in 3% of time points. During the periods of greatest deviations an observable change would have been detected in a simultaneously measured BOLD signal. Such changes could result from the imaging method or occur through glitches in hardware or alterations in the measurement environment. With identical measurement setups, the phantom could allow comparing intersession or intersubject brain activations.     

Improved functional mapping of the human amygdala using a standard functional magnetic resonance imaging sequence with simple modifications

As the amygdala is involved in various aspects of emotional processing, its characterization using neuroimaging modalities, such as functional magnetic resonance imaging (fMRI), is of great interest. However, in fMRI, the amygdala region suffers from susceptibility artifacts that are composed of signal dropouts and image distortions. Various technically demanding approaches to reduce these artifacts have been proposed, and most require alterations beyond a mere change of the acquisition parameters and cannot be easily implemented by the user without changing the MR sequence code. In the present study, we therefore evaluated the impact of simple alterations of the acquisition parameters of a standard gradient-echo echo-planar imaging technique at 3 T composed of echo times (TEs) of 27 and 36 ms as well as section thicknesses of 2 and 4 mm while retaining a section orientation parallel to the intercommissural plane and an in-plane resolution of 2x2 mm(2). In contrast to previous studies, we based our evaluation on the resulting activation maps using an emotional stimulation paradigm rather than on MR raw image quality only. Furthermore, we tested the effects of spatial smoothing of the functional raw data in the course of postprocessing using spatial filters of 4 and 8 mm. Regarding MR raw image quality, a TE of 27 ms and 2-mm sections resulted in the least susceptibility artifacts in the anteromedial aspect of the temporal lobe. The emotional stimulation paradigm resulted in robust bilateral amygdala activation for the approaches with 2-mm sections only -- but with larger activation volumes for a TE of 36 ms as compared with that of 27 ms. Moderate smoothing with a 4-mm spatial filter represented a good compromise between increased sensitivity and preserved specificity. In summary, we showed that rather than applying advanced modifications of the MR sequence, a simple increase in spatial resolution (i.e., the reduction of section thickness) is sufficient to improve the detectability of amygdala activation.     

Optimizing the imaging of the monkey auditory cortex: sparse vs. continuous fMRI

The noninvasive imaging of the monkey auditory system with functional magnetic resonance imaging (fMRI) can bridge the gap between electrophysiological studies in monkeys and imaging studies in humans. Some of the recent imaging of monkey auditory cortical and subcortical structures relies on a technique of “sparse imaging,” which was developed in human studies to sidestep the negative influence of scanner noise by adding periods of silence in between volume acquisition. Among the various aspects that have gone into the ongoing optimization of fMRI of the monkey auditory cortex, replacing the more common continuous-imaging paradigm with sparse imaging seemed to us to make the most obvious difference in the amount of activity that we could reliably obtain from awake or anesthetized animals. Here, we directly compare the sparse- and continuous-imaging paradigms in anesthetized animals. We document a strikingly greater auditory response with sparse imaging, both quantitatively and qualitatively, which includes a more expansive and robust tonotopic organization. There were instances where continuous imaging could better reveal organizational properties that sparse imaging missed, such as aspects of the hierarchical organization of auditory cortex. We consider the choice of imaging paradigm as a key component in optimizing the fMRI of the monkey auditory cortex.     

Adaptive data acquisition in MRI

We propose an adaptive data acquisition technique that depends on the object to be imaged in magnetic resonance (MR) imaging. In this paper, we employed a matching pursuit (MP) algorithm to achieve the adaptive data acquisition. Since the matching pursuit is a greedy algorithm to find RF and gradient waveforms which are the best match for an object-signal, the signal can be decomposed with a few iterations and thereby lead reduction of imaging time in MR. To adopt the matching pursuit algorithm to the adaptive data acquisition in MRI, we have designed a dictionary which contains a windowed Fourier basis set. Because the basis set is localized spatially, the image signal could be divided into segmented signals so that matching pursuit with the segmented signals could lead to effective and object-dependent data acquisition. To verify the proposed technique, computer simulations and experiments are performed with a 1.0 T whole body MRI system.     

Functional MR imaging of visual and motor cortex stimulation at high temporal resolution using a flash technique on a standard 1.5 tesla scanner

Functional magnetic resonance imaging (fMRI) was performed on a conventional 1.5 T scanner by means of a modified FLASH-technique at temporal resolutions of 80 and 320 ms. The method's stability was assessed by phantom measurements and by investigation of three volunteers resulting in a low amplitude (3%) periodic (4 s) signal modulation for the in vivo measurements, which was not observable in the phantom experiments. fMRI activation studies of motor and visual cortices of four adjacent slices were carried out on 12 healthy right-handed volunteers. Stimulation was performed by a triggered single white light flash or single finger-to-thumb opposition movement, respectively. Event-related response of visual and motor activation was traced over 10.24 s with a temporal resolution of 320 ms for the four slice measurements. Brain activation maps were calculated by correlation of measured signal time courses with a time-shifted boxcar function. Activation was quantified by calculation of percentual signal change in relation to the baseline. Observed signal magnitudes were about 5–7% in visual and about 8–12% in primary motor cortex. While photic response was delayed by about 2 s, motor stimulation showed an instantaneous increase of the MR signal. MR signal responses for both stimuli had decayed completely after about 5 s. Our results show that event-related fMRI enables mapping of brain function at sufficient spatial resolution with a temporal resolution of up to 80 ms on a conventional scanner.     

A brain phantom for motion-corrected PROPELLER showing image contrast and construction similar to those of in vivo MRI

PurposeA fast spin-echo sequence based on the Periodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction (PROPELLER) technique is a magnetic resonance (MR) imaging data acquisition and reconstruction method for correcting motion during scans. Previous studies attempted to verify the in vivo capabilities of motion-corrected PROPELLER in real clinical situations. However, such experiments are limited by repeated, stray head motion by research participants during the prescribed and precise head motion protocol of a PROPELLER acquisition. Therefore, our purpose was to develop a brain phantom set for motion-corrected PROPELLER.Materials and methodsThe profile curves of the signal intensities on the in vivo T₂-weighted image (T₂WI) and 3-D rapid prototyping technology were used to produce the phantom. In addition, we used a homemade driver system to achieve in-plane motion at the intended timing. We calculated the Pearson's correlation coefficient (R²) between the signal intensities of the in vivo T₂WI and the phantom T₂WI and clarified the rotation precision of the driver system. In addition, we used the phantom set to perform initial experiments to show the rotational angle and frequency dependences of PROPELLER.ResultsThe in vivo and phantom T₂WIs were visually congruent, with a significant correlation (R²) of 0.955 (p < .001). The rotational precision of the driver system was within 1 degree of tolerance. The experiment on the rotational angle dependency showed image discrepancies between the rotational angles. The experiment on the rotational frequency dependency showed that the reconstructed images became increasingly blurred by the corruption of the blades as the number of motions increased.ConclusionsIn this study, we developed a phantom that showed image contrasts and construction similar to the in vivo T₂WI. In addition, our homemade driver system achieved precise in-plane motion at the intended timing. Our proposed phantom set could perform systematic experiments with a real clinical MR image, which to date has not been possible in in vivo studies. Further investigation should focus on the improvement of the motion-correction algorithm in PROPELLER using our phantom set for what would traditionally be considered problematic patients (children, emergency patients, elderly, those with dementia, and so on).     

Fine-scale functional connectivity in somatosensory cortex revealed by high-resolution fMRI

High-resolution functional magnetic resonance imaging (fMRI) at high field (9.4 T) has been used to measure functional connectivity between subregions within the primary somatosensory (SI) cortex of the squirrel monkey brain. The hand-face region within the SI cortex of the squirrel monkey has been previously well mapped with functional imaging and electrophysiological and anatomical methods, and the orderly topographic map of the hand region is characterized by a lateral to medial representation of individual digits in four subregions of areas 3a, 3b, 1 and 2. With submillimeter resolution, we are able to detect not only the separate islands of activation corresponding to vibrotactile stimulations of single digits but also, in subsequent acquisitions, the degree of correlation between voxels within the SI cortex in the resting state. The results suggest that connectivity patterns are very similar to stimulus-driven distributions of activity and that connectivity varies on the scale of millimeters within the same primary region. Connectivity strength is not a reflection of global larger-scale changes in blood flow and is not directly dependent on distance between regions. Preliminary electrophysiological recordings agree well with the fMRI data. In human studies at 7 T, high-resolution fMRI may also be used to identify the same subregions and assess responses to sensory as well as painful stimuli, and to measure connectivity dynamically before and after such stimulations.     

Discussion on the choice of separated components in fMRI data analysis by spatial independent component analysis

By measuring the changes of magnetic resonance signals during a stimulation, the functional magnetic resonance imaging (fMRI) is able to localize the neural activation in the brain. In this report, we discuss the fMRI application of the spatial independent component analysis (spatial ICA), which maximizes statistical independence over spatial images. Included simulations show the possibility of the spatial ICA on discriminating asynchronous activations or different response patterns in an fMRI data set. An in vivo visual stimulation fMRI test was conducted, and the result shows a proper sum of the separated components as the final image is better than a single component, using fMRI data analysis by spatial ICA. Our result means that spatial ICA is a useful tool for the detection of different response activations and suggests that a proper sum of the separated independent components should be used for the imaging result of fMRI data processing.     

Signal profile measurements of single- and double-volume acquisitions with image-selected in vivo spectroscopy for 31P magnetic resonance spectroscopy

The volume-selection performance was studied for single- and double-volume-of-interest (VOI) acquisition with the volume-selection method image-selected in vivo spectroscopy for ³¹P magnetic resonance spectroscopy. High-resolution signal profiles were measured using a phantom simulating a brain. Inside the phantom there was a small, remotely controlled, movable signal source filled with ortho-phosphoric acid. Signal profiles of the VOI were measured in three perpendicular directions for 1VOI (single VOI) and 2VOI (double VOI) acquisition. The measured signal profiles for both acquisitions were very similar, but they showed a discrepancy with regard to the intended VOI (iVOI). The transition regions were on average 3.8 mm and the average full width at half maximum of the signal profile was 30 mm for an iVOI size of 30130130 (mm³). No displacement was observed in the signal profiles. To avoid overlapping signal profiles, the minimum separation between two iVOIs was found to be 10 mm in our magnetic resonance (MR) system. A substantial negative signal contribution from regions outside the iVOI was measured in the y-direction for 1VOI acquisition and one of the two VOIs in 2VOI acquisition. The other VOI in 2VOI acquisition exhibited only minor contamination. The measurements presented underline the importance of detailed knowledge on the volume selection performance in in vivo magnetic resonance spectroscopy.     

Real-time MR artifacts filtering during continuous EEG/fMRI acquisition

The purpose of this study was the development of a real-time filtering procedure of MRI artifacts in order to monitor the EEG activity during continuous EEG/fMRI acquisition. The development of a combined EEG and fMRI technique has increased in the past few years. Preliminary “spike-triggered” applications have been possible because in this method, EEG knowledge was only necessary to identify a trigger signal to start a delayed fMRI acquisition. In this way, the two methods were used together but in an interleaved manner. In real simultaneous applications, like event-related fMRI study, artifacts induced by MRI events on EEG traces represent a substantial obstacle for a right analysis. Up until now, the methods proposed to solve this problem are mainly based on procedures to remove post-processing artifacts without the possibility to control electrophysiological behavior of the patient during fMRI scan. Moreover, these methods are not characterized by a strong “prior knowledge” of the artifact, which is an imperative condition to avoid any loss of information on the physiological signals recovered after filtering. In this work, we present a new method to perform simultaneous EEG/fMRI study with real-time artifacts filtering characterized by a procedure based on a preliminary analytical study of EPI sequence parameters-related EEG-artifact shapes. Standard EEG equipment was modified in order to work properly during ultra-fast MRI acquisitions. Changes included: high-performance acquisition device; electrodes/cap/wires/cables materials and geometric design; shielding box for EEG signal receiver; optical fiber link; and software. The effects of the RF pulse and time-varying magnetic fields were minimized by using a correct head cap wires-locked environment montage and then removed during EEG/fMRI acquisition with a subtraction algorithm that takes in account the most significant EPI sequence parameters. The on-line method also allows a further post-processing utilization.