A comparison of fast MR scan techniques for cerebral activation studies at 1.5 tesla

Feasibility of spiral fMRI based on an LTI gradient model

Spiral imaging is very well suited for functional MRI, however its use has been limited by the fact that artifacts caused by gradient imperfections and B0 inhomogeneity are more difficult to correct compared to EPI. Effective correction requires accurate knowledge of the traversed k-space trajectory. With the goal of making spiral fMRI more accessible, we have evaluated image reconstruction using trajectories predicted by the gradient impulse response function (GIRF), which can be determined in a one-time calibration step. GIRF-predicted reconstruction was tested for high-resolution (0.8 mm) fMRI at 7T. Image quality and functional results of the reconstructions using GIRF-prediction were compared to reconstructions using the nominal trajectory and concurrent field monitoring. The reconstructions using nominal spiral trajectories contain substantial artifacts and the activation maps contain misplaced activation. Image artifacts are substantially reduced when using the GIRF-predicted reconstruction, and the activation maps for the GIRF-predicted and monitored reconstructions largely overlap. The GIRF reconstruction provides a large increase in the spatial specificity of the activation compared to the nominal reconstruction. The GIRF-reconstruction generates image quality and fMRI results similar to using a concurrently monitored trajectory. The presented approach does not prolong or complicate the fMRI acquisition. Using GIRF-predicted trajectories has the potential to enable high-quality spiral fMRI in situations where concurrent trajectory monitoring is not available.

Pneumatic artificial muscle-based stimulator for passive functional magnetic resonance imaging sensorimotor mapping in patients with brain tumours

BACKGROUND

Two concerns with respect to pre-operative task-based motor functional magnetic resonance imaging (fMRI) in patients with brain tumours are inadequate performance due to patients' impaired motor function and head motion artefacts.

NEW METHOD

In the present study we validate the use of a stimulator based on a pneumatic artificial muscle (PAM) for fMRI mapping of the primary sensorimotor (SM1) cortex in twenty patients with rolandic or perirolandic brain tumours. All patients underwent both active and passive motor block-design fMRI paradigms, performing comparable active and passive PAM-induced flexion-extensions of the icontralesional index finger.

RESULTS

PAM-induced movements resulted in a significant BOLD signal increase in contralateral primary motor (M1) and somatosensory (S1) cortices in 18/20 and 19/20 (p<.05 FWE corrected in 16/18 and 18/19) patients, versus 18/20 and 16/20 (p<.05 FWE corrected) during active movements. The two patients in whom the PAM-based stimulator failed to induce any significant BOLD signal change in the contralateral M1 cortex differed from the two in whom active motion was conversely ineffective. At the group level, no significant difference in contrast magnitude was observed within the contralateral SM1 cortex when comparing active with passive movements. During passive movements, head motion was significantly reduced. Comparison with existing method(s) As compared to the several robotic devices for passive motion that were introduced in the past decades, our PAM-based stimulator appears smaller, handier, and easier to use.

CONCLUSION

The use of PAM-based stimulators should be included in routine pre-operative fMRI protocols along with active paradigms in such patients' population.

TE-dependent spatial and spectral specificity of functional connectivity

Previous studies suggest that spontaneous fluctuations in the resting-state fMRI (RS-fMRI) signal may reflect fluctuations in transverse relaxation time (T(2)(*)) rather than spin density (S(0)). However, such S(0) and T(2)(*) features have not been well characterized. In this study, spatial and spectral characteristics of functional connectivity on sensorimotor, default-mode, dorsal attention, and primary visual systems were examined using a multiple gradient-echo sequence at 3T. In the spatial domain, we found broad, local correlations at short echo times (TE ≤ 14 ms) due to dominant S(0) contribution, whereas long-range connections mediated by T(2)(*) became explicit at TEs longer than 22 ms. In the frequency domain, compared with the flat spectrum of S(0), spectral power of the T(2)(*)-weighted signal elevated significantly with increasing TE, particularly in the frequency ranges of 0.008-0.023 Hz and 0.037-0.043 Hz. Using the S(0) spectrum as a reference, we propose two indices to measure spectral signal change (SSC) and spectral contrast-to-noise ratio (SCNR), respectively, for quantifying the RS-fMRI signal. These indices demonstrated TE dependency of connectivity-related fluctuation strength, resembling functional contrasts in activation-based fMRI. These findings further confirm that large-scale functional circuit connectivity based on BOLD contrast may be constrained within specific frequency ranges in every brain network, and the spectral features of S(0) and T(2)(*) could be valuable for interpreting and quantifying RS-fMRI data.

High-field FMRI for human applications: an overview of spatial resolution and signal specificity

In the last decade, dozens of 7 Tesla scanners have been purchased or installed around the world, while 3 Tesla systems have become a standard. This increased interest in higher field strengths is driven by a demonstrated advantage of high fields for available signal-to-noise ratio (SNR) in the magnetic resonance signal. Functional imaging studies have additional advantages of increases in both the contrast and the spatial specificity of the susceptibility based BOLD signal. One use of this resultant increase in the contrast to noise ratio (CNR) for functional MRI studies at high field is increased image resolution. However, there are many factors to consider in predicting exactly what kind of resolution gains might be made at high fields, and what the opportunity costs might be. The first part of this article discusses both hardware and image quality considerations for higher resolution functional imaging. The second part draws distinctions between image resolution, spatial specificity, and functional specificity of the fMRI signals that can be acquired at high fields, suggesting practical limitations for attainable resolutions of fMRI experiments at a given field, given the current state of the art in imaging techniques. Finally, practical resolution limitations and pulse sequence options for studies in human subjects are considered.

Spiral imaging in fMRI

T2*-weighted Blood Oxygen Level Dependent (BOLD) functional magnetic resonance imaging (fMRI) requires efficient acquisition methods in order to fully sample the brain in a several second time period. The most widely used approach is Echo Planar Imaging (EPI), which utilizes a Cartesian trajectory to cover k-space. This trajectory is subject to ghosts from off-resonance and gradient imperfections and is intrinsically sensitive to cardiac-induced pulsatile motion from substantial first- and higher order moments of the gradient waveform near the k-space origin. In addition, only the readout direction gradient contributes significant energy to the trajectory. By contrast, the spiral method samples k-space with an Archimedean or similar trajectory that begins at the k-space center and spirals to the edge (spiral-out), or its reverse, ending at the origin (spiral-in). Spiral methods have reduced sensitivity to motion, shorter readout times, improved signal recovery in most frontal and parietal brain regions, and exhibit blurring artifacts instead of ghosts or geometric distortion. Methods combining spiral-in and spiral-out trajectories have further advantages in terms of diminished susceptibility-induced signal dropout and increased BOLD signal. In measurements of temporal signal to noise ratio measured in 8 subjects, spiral-in/out exhibited significant increases over EPI in voxel volumes recovered in frontal and whole brain regions (18% and 10%, respectively).

Ultrafast bold fMRI using single-shot spin-echo echo planar imaging

The choice of imaging parameters for functional MRI can have an impact on the accuracy of functional localization by affecting the image quality and the degree of blood oxygenation-dependent (BOLD) contrast achieved. By improving sampling efficiency, parallel acquisition techniques such as sensitivity encoding (SENSE) have been used to shorten readout trains in single-shot (SS) echo planar imaging (EPI). This has been applied to susceptibility artifact reduction and improving spatial resolution. SENSE together with single-shot spin-echo (SS-SE) imaging may also reduce off-resonance artifacts. The goal of this work was to investigate the BOLD response of a SENSE-adapted SE-EPI on a three Tesla scanner. Whole-brain fMRI studies of seven healthy right hand-dominant volunteers were carried out in a three Tesla scanner. fMRI was performed using an SS-SE EPI sequence with SENSE. The data was processed using statistical parametric mapping. Both, group and individual subject data analyses were performed. Individual average percentage and maximal percentage signal changes attributed to the BOLD effect in M1 were calculated for all the subjects as a function of echo time. Corresponding activation maps and the sizes of the activated clusters were also calculated. Our results show that susceptibility artifacts were reduced with the use of SENSE; and the acquired BOLD images were free of the typical quadrature artifacts of SS-EPI. Such measures are crucial at high field strengths. SS SE-EPI with SENSE offers further benefits in this regard and is more specific for oxygenation changes in the microvasculature bed. Functional brain activity can be investigated with the help of single-shot spin echo EPI using SENSE at high magnetic fields.

Static and dynamic characteristics of cerebral blood flow during the resting state

In this study, the static and dynamic characteristics of cerebral blood flow (CBF) in the resting state were investigated using an arterial spin labeling (ASL) perfusion imaging technique. Consistent with previous PET results, static CBF measured by ASL was significantly higher in the posterior cingulate cortex (PCC), thalamus, insula/superior temporal gyrus (STG) and medial prefrontal cortex (MPFC) than the average CBF of the brain. The dynamic measurement of CBF fluctuations showed high correlation (functional connectivity) between components in the default mode network. These brain regions also had high local temporal synchrony and high fluctuation amplitude, as measured by regional homogeneity (ReHo) and amplitude of low-frequency fluctuation (ALFF) analyses. The spatial pattern of the static CBF correlated well with that of the dynamic indices. The high static and dynamic activities in the PCC, MPFC, insula/STG and thalamus suggest that these regions play a vital role in maintaining and facilitating fundamental brain functions.

Comparison of spiral imaging and SENSE-EPI at 1.5 and 3.0 T using a controlled cerebrovascular challenge

PURPOSE

To quantitatively compare spiral imaging and sensitivity-encoded-echo-planar-imaging (SENSE-EPI) methods for blood oxygen level-dependent (BOLD) imaging using controlled changes in the end-tidal partial pressure of CO(2) (PetCO(2)) to provide a global BOLD response. Specifically, we examined susceptibility-field-gradient effects on the BOLD sensitivity throughout the brain.

MATERIALS AND METHODS

We quantified cerebrovascular reactivity (CVR) using the BOLD response to cyclic changes in PetCO(2) in five healthy volunteers at 1.5 and 3.0 T using spiral imaging and SENSE-EPI. We compared the two techniques with respect to susceptibility-induced signal dropout and CVR t-statistic.

RESULTS

Compared to spiral imaging, SENSE-EPI significantly reduced the volume of signal dropout by 32 +/- 18% at 3.0 T. In regions with large susceptibility gradients, SENSE-EPI demonstrated a trend for a greater t-statistic than spiral imaging, particularly at 3.0 T. However, no statistically significant between-technique differences existed.

CONCLUSION

The results at 3.0 T suggest that, compared with spiral imaging, SENSE-EPI reduces signal loss associated with susceptibility field gradients in affected regions without affecting BOLD sensitivity. This study also demonstrates a unique application of controlled PetCO(2) changes to quantitatively compare BOLD techniques, which may be useful for the design of future fMRI studies.

Laminar profiles of functional activity in the human brain

Functional magnetic resonance imaging (fMRI) data were obtained in human visual cortex using sub-millimeter voxels at a field strength of 3 T. Reliable functional signals were largely confined to the gray matter and these responses measure the retinotopic organization of visual cortex. Functional signals were further characterized with respect to their laminar position within the cortical gray matter. The laminar response profiles during our visuospatial attention task, normalized for cortical thickness, had a stereotypical shape, with a peak in the superficial gray matter and declining in the deeper layers. The thickness of the sheet producing functional signals was in excellent agreement with the estimated structural thickness of the gray matter throughout early visual cortex (error < 0.5 mm). Thickness measurements were highly repeatable from session-to-session (error < 0.4 mm). Hence, it is feasible and useful to use high-resolution fMRI to measure laminar activity profiles. The ability to distinguish signals arising in different lamina has significant potential scientific and clinical applications.

The impact of susceptibility gradients on cartesian and spiral EPI for BOLD fMRI

High sensitivity to magnetic susceptibility changes and accurate localization of functional activations are key requisites for pulse sequences used for BOLD fMRI. This paper seeks to develop a framework for analysing the performance of various k-space sampling techniques in this respect, with special emphasis on spiral EPI (spiral) and cartesian EPI (EPI) and their performance under influence of induced field gradients (SFGs) and stochastic noise. A numerical method for calculating synthetic MR images is developed and used to simulate BOLD fMRI experiments using EPI and spirals. The data is then examined for activation using a pixel-wise t test. Nine subjects are scanned with both techniques while performing a motor task. SPM99 is used for analysing the experimental data. The simulated spirals provide generally higher t scores at low SFGs but lose more strength than EPI at higher SFGs, where EPI activation is offset from the true position. In the primary motor area spirals provide significantly higher t scores (P < 0.0002). In-plane variation of EPI is higher in phase-encoding direction than in frequency-encoding direction (P < 0.003). In the low SFG areas spirals provide stronger activation than EPI and less spatial variability. Thus, spirals are recommended for fMRI in motor area and language areas.

Anterior medial temporal lobe activation during encoding of words: FMRI methods to optimize sensitivity

The existence of a rostrocaudal gradient of medial temporal lobe (MTL) activation during memory encoding has historically received support from positron emission tomography studies, but less so from functional MRI (FMRI) studies. More recently, FMRI studies have demonstrated that characteristics of the stimuli can affect the location of activation seen in the MTL when those stimuli are encoded. The current study tested the hypothesis that MTL activation during memory encoding is related to the modality of stimulus presentation. Subjects encoded auditorily or visually presented words in an FMRI novelty paradigm. Imaging and analysis parameters were optimized to minimize susceptibility artifact in the anterior MTL. Greater activation was observed in the anterior than posterior MTL for both modalities of stimulus presentation. The results indicate that anterior MTL activation occurred during encoding, independent of stimulus modality and provide support for the hypothesis that verbal-semantic memory processing occurs in anterior MTL. The authors suggest that technical factors are critical for observing the rostrocaudal gradient in MTL memory activation.

Preoperative and postoperative mapping of cerebrovascular reactivity in moyamoya disease by using blood oxygen level-dependent magnetic resonance imaging

OBJECT

The ability to map cerebrovascular reactivity (CVR) at the tissue level in patients with moyamoya disease could have considerable impact on patient management, especially in guiding surgical intervention and assessing the effectiveness of surgical revascularization. This paper introduces a new noninvasive magnetic resonance (MR) imaging-based method to map CVR. Preoperative and postoperative results are reported in three cases to demonstrate the efficacy of this technique in assessing vascular reserve at the microvascular level.

METHODS

Three patients with angiographically confirmed moyamoya disease were evaluated before and after surgical revascularization. Measurements of CVR were obtained by rapidly manipulating end-tidal PCO2 between hypercapnic and hypocapnic states during MR imaging. The CVR maps were then calculated by comparing the percentage of changes in MR signal with changes in end-tidal PCO2. Presurgical CVR maps showed distinct regions of positive and negative reactivity that correlated precisely with the vascular territories supplied by severely narrowed vessels. Postsurgical reactivity maps demonstrated improvement in the two patients with positive clinical outcome and no change in the patient in whom a failed superficial temporal artery-middle cerebral artery bypass was performed.

CONCLUSIONS

Magnetic imaging-based CVR mapping during rapid manipulation of end-tidal PCO2 is an exciting new method for determining the location and extent of abnormal vascular reactivity secondary to proximal large-vessel stenoses in moyamoya disease. Although the study group is small, there seems to be considerable potential for guiding preoperative decisions and monitoring efficacy of surgical revascularization.

Resolution and reproducibility of BOLD and perfusion functional MRI at 3.0 Tesla

Visual and somatosensory activation studies were performed on normal subjects to compare the spatial discrimination and reproducibility between functional MRI (fMRI) methods based on blood oxygen level-dependent (BOLD) and perfusion contrast. To allow simultaneous measurement of BOLD and perfusion contrast, a dedicated MRI acquisition technique was developed. Repeated experiments of sensory stimulation of single digits of the right hand showed an average variability of activation amplitude of 25% for BOLD data, and a significantly lower variability of 21% for perfusion data. No significant difference in the variability of the locus of activity was observed between the BOLD and perfusion data. In somatotopy experiments, digits II and V were subjected to passive sensory stimulation. Both the BOLD and perfusion data showed substantial overlap in the activation patterns from the two digits. In a retinotopy study, two stimuli were alternated to excite different patches of V1. Again there was substantial overlap between the activation patterns from both stimuli, although the perfusion performed somewhat better than the BOLD method. Particularly for the visual studies, the overlap in activation patterns was more than expected based on the fine-scale retinotopic mapping of cortical activity, suggesting that both BOLD and perfusion contrast mechanisms contribute substantially to the point-spread function (PSF).

Functional brain imaging using a blood oxygenation sensitive steady state

Blood oxygenation level dependent (BOLD) functional MRI (fMRI) is an important method for functional neuroimaging that is sensitive to changes in blood oxygenation related to brain activation. While BOLD imaging has good spatial coverage and resolution relative to other neuroimaging methods (such as positron emission tomography (PET)), it has significant limitations relative to other MRI techniques, including poor spatial resolution, low signal levels, limited contrast, and image artifacts. These limitations derive from the coupling of BOLD functional contrast to sources of image degradation. This work presents an alternative method for fMRI that may over-come these limitations by establishing a blood oxygenation sensitive steady-state (BOSS) that inverts the signal from deoxygenated blood relative to the water signal. BOSS fMRI allows the imaging parameters to be optimized independently of the functional contrast, resulting in fewer image artifacts and higher signal-to-noise ratio (SNR). In addition, BOSS fMRI has greater functional contrast than BOLD. BOSS fMRI requires careful shimming and multiple acquisitions to obtain a precise alignment of the magnetization to the SSFP frequency response.

A PRESTO-SENSE sequence with alternating partial-Fourier encoding for rapid susceptibility-weighted 3D MRI time series

A 3D sequence for dynamic susceptibility imaging is proposed which combines echo-shifting principles (such as PRESTO), sensitivity encoding (SENSE), and partial-Fourier acquisition. The method uses a moderate SENSE factor of 2 and takes advantage of an alternating partial k-space acquisition in the "slow" phase encode direction allowing an iterative reconstruction using high-resolution phase estimates. Offering an isotropic spatial resolution of 4 x 4 x 4 mm(3), the novel sequence covers the whole brain including parts of the cerebellum in 0.5 sec. Its temporal signal stability is comparable to that of a full-Fourier, full-FOV EPI sequence having the same dynamic scan time but much less brain coverage. Initial functional MRI experiments showed consistent activation in the motor cortex with an average signal change slightly less than that of EPI.

Single-shot MR imaging using trapezoidal-gradient-based Lissajous trajectories

A novel single-shot trapezoidal-gradient-based Lissajous trajectory is described and implemented on a 3-tesla magnetic resonance (MR) scanner. A feature of this trajectory is that its sampling points are located on a nonequidistant rectangular grid, which permits the usage of one-dimensional optimal algorithms to increase the robustness and speed of image reconstruction. Another advantage of the trajectory is that two images with different effective echo times can be obtained within a single excitation, which might be used for fast T2* mapping, in functional MR imaging scanning of brain activity associated with mental processes. Potential artifacts in reconstructed images were investigated and methods for suppressing these artifacts were developed. Experiments on normal subjects at rest and during brain activation were performed to demonstrate the feasibility of the new sequence.

Functional magnetic resonance imaging: a review of methodological aspects and clinical applications

This paper gives an overview of the recent literature on methodological developments of functional magnetic resonance imaging (fMRI) and recent trends in clinical applications. With the recent introduction of high-field systems and methodological developments leading to more robust signal behavior, fMRI is in a transition state from a research modality for use by experts to a standard procedure with useful applications in patient management. Compared to the use in neuroscientific research, which is often based on BOLD techniques alone, the application in patients is distinguished by a multiparametric characterization of the brain using a combination of several techniques. Neuronal fiber tracking based on diffusion anisotropy measurements, in particular, has already turned out to provide relevant supplementary information to the BOLD-based cortical activation maps.

EPI-BOLD fMRI of human motor cortex at 1.5 T and 3.0 T: Sensitivity dependence on echo time and acquisition bandwidth

PURPOSE

To investigate the sensitivity dependence of BOLD functional imaging on MRI acquisition parameters in motor stimulation experiments using a finger tapping paradigm.

MATERIALS AND METHODS

Gradient-echo echo-planar fMRI experiments were performed at 1.5 T and 3.0 T with varying acquisition echo time and bandwidth, and with a 4 mm isotropic voxel size. To analyze the BOLD sensitivity, the relative contributions of BOLD signal amplitude and thermal and physiologic noise sources were evaluated, and statistical t-scores were compared in the motor area.

RESULTS

At 1.5 T, the number of activated pixels and the average t-score showed a relatively broad optimum over a TE range of 60-160 msec. At 3.0 T, an optimum range was observed between TEs of 30-130 msec. Averaged over nine subjects, maxima in the number of pixels and t-score values were 59% and 18% higher at 3.0 T than at 1.5 T, respectively, an improvement that was lower than the observed 100% to 110% increase in signal-to-noise ratio at 3.0 T.

CONCLUSION

The somewhat disappointing increase in t-scores at 3.0 T was attributed to the increased contribution of physiologic noise at the higher field strength under the given experimental conditions. At both field strengths, reducing the effective image acquisition bandwidth from 35 to 17 Hz per pixel did not affect or only marginally affect the BOLD sensitivity.

Application of sensitivity-encoded echo-planar imaging for blood oxygen level-dependent functional brain imaging

The benefits of sensitivity-encoded (SENSE) echo-planar imaging (EPI) for functional MRI (fMRI) based on blood oxygen level-dependent (BOLD) contrast were quantitatively investigated at 1.5 T. For experiments with 3.4 x 3.4 x 4.0 mm(3) resolution, SENSE allowed the single-shot EPI image acquisition duration to be shortened from 24.1 to 12.4 ms, resulting in a reduced sensitivity to geometric distortions and T(*)(2) blurring. Finger-tapping fMRI experiments, performed on eight normal volunteers, showed an overall 18% loss in t-score in the activated area, which was substantially smaller than expected based on the image signal-to-noise ratio (SNR) and g-factor, but similar to the loss predicted by a model that takes physiologic noise into account.

Simultaneous perfusion and BOLD imaging using reverse spiral scanning at 3T: characterization of functional contrast and susceptibility artifacts

Reverse spiral scanning with arterial spin-labeling was developed at 3T to simultaneously detect perfusion and BOLD signals in the brain by subtracting or adding the control and labeled images, respectively, in the same dataset. BOLD contrast was improved with the longer effective echo time achieved in the reverse spiral scan compared to conventional forward spiral scans. Susceptibility artifacts near air-tissue interfaces in the brain were substantially reduced in the reverse spiral images due to their early data acquisition time and, hence, less signal attenuation. Brain activation experiments with the reverse spiral scan were performed on normal subjects and were compared to forward spiral imaging in the same subjects. The experiments demonstrated that reverse spiral imaging was able to detect perfusion and BOLD signals simultaneously and reliably, even in the brain regions with severe susceptibility-induced local gradients, while forward spiral scans were either not optimal for detecting the two functional signals at the same time or were vulnerable to susceptibility artifacts.

Effect of attention on central auditory processing: an fMRI study

Functional magnetic resonance imaging was used to investigate preattentive and attentional processing of auditory stimuli in 18 right-handed normal volunteers. Responses to trains of 1000-Hz pure tones and infrequent (15%) deviant 1300-Hz tones were characterized while subjects ignored all tones; listened for deviants in the left ear; or listened for deviants in the right ear. Preattentive detection of deviants, associated with the mismatch negativity in electrophysiology, was associated with bilateral temporal lobe activation, with a rightward predominance. Processing of deviant stimuli while attending to either ear produced a more robust and widespread activation of these temporal regions, again with a rightward predominance. Thus, preattentive tone processing appears to be linked to asymmetric activation of a core set of temporal regions in which activity is significantly amplified by selective attention. Extratemporal regions activated by attending to targets in either ear included the anterior cingulate cortex, supramarginal gyrus, and dorsolateral prefrontal cortex.

fMRI study comparing names versus pictures of objects

We performed an fMRI one-back recognition study aimed at distinguishing the semantic versus perceptual aspects of how objects and their written forms are processed. There were three types of visually presented items: pictures (schematic drawings of objects); words identifying these objects; and a mixed condition in which pictures were interleaved with words. A semantic decision about object identity was required when pictures were interleaved with words. This condition, contrasted with the other two, invoked a larger signal in multiple areas, including frontal cortex, bilateral occipitotemporal cortex, and the right middle temporal gyrus. We propose that the left occipitotemporal and right temporal activations are indicative of the neural substrate mediating picture-word conversions, whereas the frontal activations reflect the coordinating functions of the central executive.

Repetition time in echo planar functional MRI

To date, surprisingly little attention has been directed toward determining the optimum TR in a functional imaging experiment. A survey of the literature reveals a wide range of TRs, but little justification for a specific TR. Long-TR functional imaging experiments provide maximum signal-to-noise ratio (SNR) in the raw images; allow for the collection of a large number of slice locations; and decrease the size of the data set acquired, simplifying storage and handling. This work, however, demonstrates that long-TR imaging sacrifices statistical power when the paradigm timing is held fixed. That is, for a fixed-run duration consisting of multiple activation/control blocks, shorter TR acquisitions (on the order of 1000 ms) provide better discrimination between the activated and nonactivated brain tissue regions than do long-TR acquisitions (on the order of 4000 ms). Results are shown for modeling the functional imaging experiment and for three different paradigms performed on normal subjects.

Quantifying head motion associated with motor tasks used in fMRI

In functional magnetic resonance imaging (fMRI) studies, long experiment times and small intensity changes associated with brain activation frequently lead to image artifacts due to head motion. Methods to minimize and correct for head motion by restraint, fast imaging, and retrospective image registration are typically combined but do not completely solve the problem, particularly for specific patient populations. As an initial step toward optimizing future designs of head restraints and improving motion correction techniques, the head motion characteristics of groups of stroke subjects, age-matched controls, and young adults were investigated with the aid of an MR simulator and a highly accurate position tracking system. Position measurements were recorded during motor tasks involving either the hand or the foot. Head motion was strongly dependent on the subject group and less upon the task conditions based on ANOVA calculations (P < 0.05). The stroke subjects exhibited approximately twice the head motion compared to that of age-matched controls, and the latter's head motion was about twice that of young adults. Moreover, the range of head motion in stroke subjects over all tasks was approximately 2 +/- 1 mm, with the motion occurring predominantly as translation in the superior-inferior direction and pitch rotation (nodding). These results lead to several recommendations on the design of fMRI motor experiments and suggest that improved motion correction strategies are required to examine such patient populations comprehensively.

High-sensitivity single-shot perfusion-weighted fMRI

A method is presented for measurement of perfusion changes during brain activation using a single-shot pulsed spin labeling technique. By employing a double-inversion labeling strategy, stationary tissue (background) signal was suppressed while minimally affecting perfusion sensitivity. This allowed omission of the otherwise required reference scan, resulting in twofold-improved temporal resolution. The method was applied to visual and motor cortex activation studies in humans, and compared to standard FAIR-type perfusion labeling techniques. Experiments performed at 1.5T and 3.0T indicate a close to 90% suppression of background signal, at a cost of an 11% and 9%, respectively, reduction in perfusion signal. Combined with the twofold increase in signal averaging, and a reduction in background signal fluctuations, this resulted in a 64% (1.5T, N = 3) and a 128% (3T, N = 4) overall improvement in sensitivity for the detection of activation-related perfusion changes. Magn Reson Med 46:88-94, 2001. Published 2001 Wiley-Liss, Inc.

Single-shot T2(*) measurement to establish optimum echo time for fMRI: studies of the visual, motor, and auditory cortices at 3.0 T

The signal change in fMRI is dependent on the echo time and the rate of decay of transverse magnetization. The latter factor may vary across regions of the brain as a result of variations in field homogeneity. Previous measurements of the signal change with echo time have generally employed relatively slow multi-echo techniques, which may be sensitive to movement and habituation effects. Here a fast T(2)(*) measurement technique, involving the generation of six low-resolution echo planar images from a single FID, is described, and its use in the evaluation of the optimum echo time for visual, motor, and auditory fMRI experiments at 3.0 T is outlined.

Transit time, trailing time, and cerebral blood flow during brain activation: measurement using multislice, pulsed spin-labeling perfusion imaging

Transit time and trailing time in pulsed spin-labeling perfusion imaging are likely to be modulated by local blood flow changes, such as those accompanying brain activation. The majority of transit/trailing time is due to the passage of the tagged blood bolus through the arteriole/capillary regions, because of lower blood flow velocity in these regions. Changes of transit/trailing time during activation could affect the quantification of CBF in functional neuroimaging studies, and are therefore important to characterize. In this work, the measurement of transit and trailing times and CBF during sensorimotor activation using multislice perfusion imaging with pulsed arterial spin-labeling is described. While CBF elevated dramatically ( thick similar80.7%) during the sensorimotor activation, sizable reductions of transit time ( thick similar0.11 sec) and trailing time ( thick similar0.26 sec) were observed. Transit and trailing times were dependent on the distances from the leading and trailing edges of the tagged blood bolus to the location of the imaging slices. The effects of transit/trailing time changes on CBF quantification during brain activation were analyzed by simulation studies. Significant errors can be caused in the estimation of CBF if such changes of transit/trailing time are not taken into account.

A CBF-based event-related brain activation paradigm: characterization of impulse-response function and comparison to BOLD

A perfusion-based event-related functional MRI method for the study of brain activation is presented. In this method, cerebral blood flow (CBF) was measured using a recently developed multislice arterial spin-labeling (ASL) perfusion imaging method with rapid spiral scanning. Temporal resolution of the perfusion measurement was substantially improved by employing intertrial subtraction and stimulus-shifting schemes. Perfusion and blood oxygenation level-dependent (BOLD) signals were obtained simultaneously by subtracting or adding the control and labeled images, respectively, in the same data sets. The impulse response function (IRF) of perfusion during brain activation was characterized for multiple stimulus durations and compared to the simultaneously acquired BOLD response. The CBF response curve preceded the BOLD curve by 0.21 s in the rising phase and 0.64 s in the falling phase. Linear additivity of the CBF and BOLD responses was assessed with rapidly repeated stimulations within single trials, and departure from linearity was found in both responses, characterized as attenuated amplitude and delayed rising time. Event-related visual and sensorimotor activation experiments were successfully performed with the new perfusion technique.

PRESTO-SENSE: an ultrafast whole-brain fMRI technique

A new ultrafast MR imaging method is proposed and tested, which enables whole-brain fMRI with a true temporal resolution of 1 sec. The method combines a 3D PRESTO pulse sequence with the concept of sensitivity-encoding with multiple receiver coils (SENSE). The so-called PRESTO-SENSE technique is demonstrated on a set of functional block-type motor and visual experiments and compared with conventional functional imaging techniques, such as PRESTO and EPI. Comparable image quality and activation areas are found with all sequences. The noise characteristics of the proposed method are analyzed in detail and their implications for ultrafast fMRI studies are discussed. Magn Reson Med 43:779-786, 2000.

Functional MRI of the human motor cortex using single-shot, multiple gradient-echo spiral imaging

In this study, we combined the advantages of a fast multi-slice spiral imaging approach with a multiple gradient-echo sampling scheme at high magnetic field strength to improve quantification of BOLD and inflow effects and to estimate T2* relaxation times in functional brain imaging. Eight echoes are collected with echo time (TE) ranging from 5 to 180 ms. Acquisition time per slice and echo time is 25 ms for a nominal resolution of 4 x 4 x 4 mm3. Evaluation of parameter images during rest and stimulation yields no significant activation on the inflow sensitive spin-density images (rho or I0-maps) whereas clear activation patterns in primary human motor cortex (M1) and supplementary motor area (SMA) are detected on BOLD sensitive T2*-maps. The calculation of relaxation times and rates of the activated areas over all subjects yields an average T2* +/- standard deviation (SD) of 46.1+/-4.5 ms (R2* of 21.8+/-2.2 s(-1)) and an average increase (deltaT2* +/- SD) of 0.93+/-0.47 ms (deltaR2* of -0.4+/-0.14 s(-1)). Our findings demonstrate the usefulness of a multiple gradient echo data acquisition approach in separating various vascular contributions to brain activation in fMRI.

Functional magnetic resonance imaging: imaging techniques and contrast mechanisms

Functional magnetic resonance imaging (fMRI) is a widely used technique for generating images or maps of human brain activity. The applications of the technique are widespread in cognitive neuroscience and it is hoped they will eventually extend into clinical practice. The activation signal measured with fMRI is predicated on indirectly measuring changes in the concentration of deoxyhaemoglobin which arise from an increase in blood oxygenation in the vicinity of neuronal firing. The exact mechanisms of this blood oxygenation level dependent (BOLD) contrast are highly complex. The signal measured is dependent on both the underlying physiological events and the imaging physics. BOLD contrast, although sensitive, is not a quantifiable measure of neuronal activity. A number of different imaging techniques and parameters can be used for fMRI, the choice of which depends on the particular requirements of each functional imaging experiment. The high-speed MRI technique, echo-planar imaging provides the basis for most fMRI experiments. The problems inherent to this method and the ways in which these may be overcome are particularly important in the move towards performing functional studies on higher field MRI systems. Future developments in techniques and hardware are also likely to enhance the measurement of brain activity using MRI.

Gradient echo time dependence and quantitative parameter maps for somatosensory activation in rats at 7 T

The dependence of functional magnetic resonance imaging (MRI) contrast on the gradient echo time TE in T2*-weighted blood oxygenation level-dependent (BOLD) fast low-angle shot (FLASH) imaging has been studied at 7 T for electrical forepaw stimulation in alpha-chloralose anesthetized rats. The observed variation of both the activation signal intensity and spatial pattern with echo time TE, resulting from the regional heterogeneity of T2*, was assessed by the calculation of quantitative T2* and quantitative STE = 0 maps, the latter representing the back-extrapolated signal intensity for TE = 0. The subsequently determined T2* and STE = 0 activation maps allowed a pixelwise separation of true BOLD from inflow contributions to forepaw stimulation-induced signal change in the somatosensory cortex of rat brain. For functional activation experiments performed with one single echo time the prior measurement of a quantitative T2* map is recommended as minimum further information to judge the intensity and the regional pattern of the resulting activation maps.

Investigation of low frequency drift in fMRI signal

Low frequency drift (0.0-0.015 Hz) has often been reported in time series fMRI data. This drift has often been attributed to physiological noise or subject motion, but no studies have been done to test this assumption. Time series T*2-weighted volumes were acquired on two clinical 1.5 T MRI systems using spiral and EPI readout gradients from cadavers, a normal volunteer, and nonhomogeneous and homogeneous phantoms. The data were tested for significant differences (P = 0.001) from Gaussian noise in the frequency range 0.0-0.015 Hz. The percentage of voxels that were significant in data from the cadaver, normal volunteer, nonhomogeneous and homogeneous phantoms were 13.7-49.0%, 22.1-61.9%, 46.4-68.0%, and 1.10%, respectively. Low frequency drift was more pronounced in regions with high spatial intensity gradients. Significant drifting was present in data acquired from cadavers and nonhomogeneous phantoms and all pulse sequences tested, implying that scanner instabilities and not motion or physiological noise may be the major cause of the drift.

Comparison of 3D BOLD functional MRI with spiral acquisition at 1.5 and 4.0 T

In order to investigate the merit of high field strength for BOLD-contrast-based functional magnetic resonance imaging (fMRI) studies, multishot gradient-echo fMRI experiments during motor cortex activation were performed on 1.5- and 4.0-T scanners with equivalent hardware, on the same volunteers. In these studies, artifactual vascular enhancement related to inflow effects was minimized, and large brain areas were covered by using a 3D scan technique. Temporal signal stability was optimized by using spiral readout gradients. The sensitivity for detection of activated regions was assessed by measuring the number of "activated voxels" and their average t score in predefined regions of interest. When comparing fMRI experiments with the same total scan time, performed on six subjects, and with acquisition parameters optimized for each field strength separately, the 4.0-T scanner proved to give superior results, with a 70% greater number of activated voxels and a 20% higher average t score for the activated voxels.

The effect of slice order and thickness on fMRI activation data using multislice echo-planar imaging

Multislice echo-planar imaging (EPI) is a commonly used technique for fMRI studies. Brain activation images acquired using fMRI are sensitive to T2* changes, reflecting the level of blood oxygenation (BOLD contrast), and may also contain an element of T1 contrast which detects blood flow changes in large vessels. If slice inflow (T1) effects are significant in multislice EPI, then as the order in which the slices are acquired is changed, differences in the activation maps are predicted. However, in experiments presented here using visual stimulation, the data demonstrate that highly consistent results can be achieved for repetition times (TR) of 6.0, 3.0, and 1.5 s. This suggests that, for whole-brain multislice EPI, fMRI activation is dominated by T2*, BOLD contrast. The thickness of the imaging slice is also an important parameter in these studies, having implications for spatial resolution, sensitivity, and acquisition time. In separate visual cortex experiments the effect on the values of the fMRI Z scores and the number of activated voxels is investigated as a function of slice thickness (from 1 to 8 mm). The maximum Z scores in the data are similar for all slice thicknesses and, after resampling to allow a direct comparison to be made, the volume of visual cortex detected as significantly activated increases with slice thickness.

Submillimeter functional localization in human striate cortex using BOLD contrast at 4 Tesla: implications for the vascular point-spread function

Using multislice segmented echoplanar imaging at 4 T, we have measured an upper bound to the cortical vasculature point-spread function (PSF) using functional magnetic resonance imaging (fMRI) in humans. Our experiments demonstrate that cortical subunits that are approximately 700 microm apart can be resolved using the early part of the hyperoxygenation phase of the blood oxygenation level-dependent (BOLD) effect. This was accomplished using brief (4 sec) single trials of monocular and binocular stimulation of ocular dominance columns in human primary visual cortex. The data suggest that at even higher magnetic fields, the cortical vasculature PSF may be limited by the extent and nature of horizontal connections and not signal-to-noise ratio.

Optimization of fast acquisition methods for whole-brain relative cerebral blood volume (rCBV) mapping with susceptibility contrast agents

Fast gradient-echo magnetic resonance scan techniques with spiral and rectilinear (echoplanar) k-space trajectories were optimized to perform bolus-tracking studies of human brain. Cerebral hemodynamics were studied with full brain coverage, a spatial resolution of 4 mm, and a temporal resolution of 2 seconds. The sensitivity of the techniques to detect image signal-intensity changes during the first pass of the contrast agent was studied at a range of TEs using dedicated experiments. For single-shot versions of spiral scanning and echoplanar imaging techniques with a 0.1-mmol/kg injection of gadolinium diethylenetriamine pentaacetic acid using a mechanical injector at 10 mL/sec under 1.5 T, the maximum sensitivity was obtained at TEs between 35 and 45 msec. At TEs less than 35 msec, signal-intensity artifacts were observed in the images. Analysis of the point-spread function revealed that susceptibility changes induced by the contrast agent can result in signal shifts to neighboring voxels. These artifacts are attributed to susceptibility-related signal changes during the acquisition window.

Technical solution for an interactive functional MR imaging examination: application to a physiologic interview and the study of cerebral physiology

Studies with functional magnetic resonance (MR) imaging produce large unprocessed raw data sets in minutes. The analysis usually requires transferring of the data to an off-line workstation, and this process frequently occurs after the subject has left the MR unit. The authors describe a hardware configuration and processing software that captures whole-brain raw data files as they are being produced from the MR unit. It then performs the reconstruction, registration, and statistical analysis, and displays the results in seconds after completion of the MR image acquisition.

Multislice imaging of quantitative cerebral perfusion with pulsed arterial spin labeling

A method is presented for multislice measurements of quantitative cerebral perfusion based on magnetic labeling of arterial spins. The method combines a pulsed arterial inversion, known as the FAIR (Flow-sensitive Alternating Inversion Recovery) experiment, with a fast spiral scan image acquisition. The short duration (22 ms) of the spiral data collection allows simultaneous measurement of up to 10 slices per labeling period, thus dramatically increasing efficiency compared to current single slice acquisition protocols. Investigation of labeling efficiency, suppression of unwanted signals from stationary as well as intraarterial spins, and the FAIR signal change as a function of inversion delay are presented. The assessment of quantitative cerebral blood flow (CBF) with the new technique is demonstrated and shown to require measurement of arterial transit time as well as suppression of intraarterial spin signals. CBF values measured on normal volunteers are consistent with results obtained from H2O15 positron emission tomography (PET) studies and other radioactive tracer approaches. In addition, the new method allows detection of activation-related perfusion changes in a finger-tapping experiment, with locations of activation corresponding well to those observed with blood oxygen level dependent (BOLD) fMRI.