Fast 3D functional magnetic resonance imaging at 1.5 T with spiral acquisition

Mono-planar T-Hex: Speed and flexibility for high-resolution 3D imaging

Purpose

The aim of this work is the reconciliation of high spatial and temporal resolution for MRI. For this purpose, a novel sampling strategy for 3D encoding is proposed, which provides flexible k‐space segmentation along with uniform sampling density and benign filtering effects related to signal decay.

Methods

For time‐critical MRI applications such as functional MRI (fMRI), 3D k‐space is usually sampled by stacking together 2D trajectories such as echo planar imaging (EPI) or spiral readouts, where each shot covers one k‐space plane. For very high temporal and medium to low spatial resolution, tilted hexagonal sampling (T‐Hex) was recently proposed, which allows the acquisition of a larger k‐space volume per excitation than can be covered with a planar readout. Here, T‐Hex is described in a modified version where it instead acquires a smaller k‐space volume per shot for use with medium temporal and high spatial resolution.

Results

Mono‐planar T‐Hex sampling provides flexibility in the choice of speed, signal‐to‐noise ratio (SNR), and contrast for rapid MRI acquisitions. For use with a conventional gradient system, it offers the greatest benefit in a regime of high in‐plane resolution <1 mm. The sampling scheme is combined with spirals for high sampling speed as well as with more conventional EPI trajectories.

Conclusion

Mono‐planar T‐Hex sampling combines fast 3D encoding with SNR efficiency and favorable depiction characteristics regarding noise amplification and filtering effects from T2∗ decay, thereby providing flexibility in the choice of imaging parameters. It is attractive both for high‐resolution time series such as fMRI and for applications that require rapid anatomical imaging.

The feasibility of quantitative MRI of perivascular spaces at 7T

BACKGROUND

Dilated brain perivascular spaces (PVSs) are found to be associated with many conditions, including aging, dementia, and Alzheimer's disease (AD). Conventionally, PVS assessment is mainly based on subjective observations of the number, size and shape of PVSs in MR images collected at clinical field strengths (≤3T). This study tests the feasibility of imaging and quantifying brain PVS with an ultra-high 7T whole-body MRI scanner.

NEW METHOD

3D high resolution T2-weighted brain images from healthy subjects (n=3) and AD patients (n=5) were acquired on a 7T whole-body MRI scanner. To automatically segment the small hyperintensive fluid-filling PVS structures, we also developed a quantitative program based on algorithms for spatial gradient, component connectivity, edge-detection, k-means clustering, etc., producing quantitative results of white matter PVS volume densities.

RESULTS

The 3D maps of automatically segmented PVS show an apparent increase in PVS density in AD patients compared to age-matched healthy controls due to the PVS dilation (8.0±2.1 v/v% in AD vs. 4.9±1.3 v/v% in controls, p<0.05).

COMPARISON WITH EXISTING METHOD

We demonstrated that 7T provides sufficient SNR and resolution for quantitatively measuring PVSs in deep white matter that is challenging with clinical MRI systems (≤3T). Compared to the conventional visual counting and rating for the PVS assessment, the quantitation method we developed is automatic and objective.

CONCLUSIONS

Quantitative PVS MRI at 7T may serve as a non-invasive and endogenous imaging biomarker for diseases with PVS dilation.

Simultaneous functional MRI acquisition of distributed brain regions with high temporal resolution using a 2D-selective radiofrequency excitation

PURPOSE

To perform simultaneous functional MRI of multiple, distributed brain regions at high temporal resolution using a 2D-selective radiofrequency (2DRF) excitation.

METHODS

A tailored 2DRF excitation is used to excite several, small regions-of-interest distributed in the brain. They are acquired in a single projection image with an appropriately chosen orientation such that the different regions-of-interest can be discriminated by their position in the projection plane. Thus, they are excited and acquired simultaneously with a temporal resolution comparable to that of a single-slice measurement. The feasibility of this approach for functional neuroimaging (in-plane resolution 2 × 2 mm(2) ) at high temporal resolution (80 ms) is demonstrated in healthy volunteers for regions-of-interest in the visual and motor system using checkerboard and finger tapping block-design paradigms.

RESULTS

Task-related brain activation could be observed in both the visual and the motor system simultaneously with a high temporal resolution. For an onset shift of 240 ms for half of the checkerboard, a delay of the hemodynamic response in the corresponding hemisphere of the visual cortex could be detected.

CONCLUSION

Limiting the excited magnetization to the desired target regions with a 2DRF excitation reduces the imaging sampling requirements which can improve the temporal resolution significantly.

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

Quantification of pain-induced changes in cerebral blood flow by perfusion MRI

The purpose of this study was to assess if the functional activation caused by painful stimuli could be detected with arterial spin labeling (ASL), which is a non-invasive magnetic resonance imaging (MRI) technique for measuring cerebral blood flow (CBF). Because ASL directly measures blood flow, it is well suited to pain conditions that are difficult to assess with current functional MRI, such as chronic pain. However, the use of ASL in neuroimaging has been hampered by its low sensitivity. Recent improvements in MRI technology, namely increased magnetic field strengths and phased array receiver coils, should enable ASL to measure the small changes in CBF associated with pain. In this study, healthy volunteers underwent two ASL imaging sessions, during which a painful thermal stimulus was applied to the left hand. The results demonstrated that the ASL technique measured changes in regional CBF in brain regions that have been previously identified with pain perception. These included bilateral CBF changes in the insula, secondary somatosensory, and cingulate cortices, as well as the supplementary motor area (SMA). Also observed were contralateral primary somatosensory and ipsilateral thalamic CBF changes. The average change in CBF for all regions of interest was 3.68ml/100g/min, ranging from 2.97ml/100g/min in ipsilateral thalamus to 4.91ml/100g/min in contralateral insula. The average resting global CBF was 54+/-9.7ml/100g/min, and there was no change in global CBF due to the noxious thermal stimulus.

Reconstruction of undersampled dynamic images by modeling the motion of object elements

Dynamic MRI is restricted due to the time required to obtain enough data to reconstruct the image sequence. Several undersampled reconstruction techniques have been proposed to reduce the acquisition time. In most of these techniques the nonacquired data are recovered by modeling the temporal information as varying pixel intensities represented in time or in temporal frequencies. Here we propose a new approach that recovers the missing data through a motion estimation of the object elements ("obels," or pieces of tissue) of the image. This method assumes that an obel displacement through the sequence has lower bandwidth than fluctuations in pixel intensities caused by the motion, and thus it can be modeled with fewer parameters. Preliminary results show that this technique can effectively reconstruct (with root mean square (RMS) errors below 4%) cardiac images and joints with undersampling factors of 8 and 4, respectively. Moreover, in the reconstruction process an approximation of the motion vectors is obtained for each obel, which can be used to quantify dynamic information. In this method the motion need not be confined to a part of the field of view (FOV) or to a portion of the temporal frequency. It is appropriate for dynamic studies in which the obels' motion model has fewer parameters than the number of acquired samples.

Partial-k-space acquisition method for improved SNR efficiency and temporal resolution in 3D fMRI

Previous studies have shown the relative importance of physiological noise and thermal noise in 2D MR images. Since physiological noise is proportional to the signal, it can be the dominant component at the center of k-space. In this study we demonstrate that the signal-to-noise ratio (SNR) efficiency and temporal resolution for 3D functional MRI (fMRI) are increased by the use of a partial-k-space acquisition method. In partial-k-space methods, the high-spatial-frequency components are doubled in amplitude during reconstruction, resulting in twice as much noise from those components. However, in sum these contributions are relatively small compared to those at the low spatial frequencies, where physiological noise is dominant. Therefore, the effect on the final MR images is almost negligible due to the square summation rule. Thus, the partial-k-space 3D method sacrifices much less SNR than is expected from the thermal noise model, and the SNR efficiency is increased compared to a full-k-space acquisition since more time frames can be collected for the same scan time. Accordingly, the temporal resolution can be increased in 3D acquisitions because only partial coverage of k-space is necessary. Experimental results confirm that more activation with a higher average t-score is detected by this method.

Neurophysiological correlates of age-related changes in working memory capacity

Cognitive abilities such as working memory (WM) capacity decrease with age. To determine the neurophysiological correlates of age-related reduction in working memory capacity, we studied 10 young subjects (<35 years of age; mean age=29) and twelve older subjects (>55 years of age; mean age=59) with whole brain blood oxygenation-level dependent (BOLD) fMRI on a 1.5 T GE MR scanner using a SPIRAL FLASH pulse sequence (TE=24 ms, TR=56 ms, FA=60 degrees , voxel dimensions=3.75 mm(3)). Subjects performed a modified version of the "n" back working memory task at different levels of increasing working memory load (1-Back, 2-Back and 3-Back). Older subjects performed as well as the younger subjects at 1-Back (p=0.4), but performed worse than the younger subjects at 2-Back (p<0.01) and 3-Back (p=0.06). Older subjects had significantly longer reaction time (RT) than younger subjects (p<0.04) at all levels of task difficulty. Image analysis using SPM 99 revealed a similar distribution of cortical activity between younger and older subjects at all task levels. However, an analysis of variance revealed a significant group x task interaction in the prefrontal cortex bilaterally; within working memory capacity, as in 1-Back when the older subjects performed as well as the younger subjects, they showed greater prefrontal cortical (BA 9) activity bilaterally. At higher working memory loads, however, when they performed worse then the younger subjects, the older subjects showed relatively reduced activity in these prefrontal regions. These data suggest that, within capacity, compensatory mechanisms such as additional prefrontal cortical activity are called upon to maintain proficiency in task performance. As cognitive demand increases, however, they are pushed past a threshold beyond which physiological compensation cannot be made and, a decline in performance occurs.

Noise reduction in multi-slice arterial spin tagging imaging

Attenuating the static signal in arterial spin tagging (ASSIST) was initially developed for 3D imaging of cerebral blood flow. To enable the simultaneous collection of cerebral blood flow and BOLD data, a multi-slice version of ASSIST is proposed. As with the 3D version, this sequence uses multiple inversion pulses during the tagging period to suppress the static signal. To maintain background suppression in all slices, the multi-slice sequence applies additional inversion pulses between slice acquisitions. The utility of the sequence was demonstrated by simultaneously acquiring ASSIST and BOLD data during a functional task and by collecting resting-state ASSIST data over a large number of slices. In addition, the temporal stability of the perfusion signal was found to be 60% greater at 3 T compared to 1.5 T, which was attributed to the insensitivity of ASSIST to physiologic noise.

Comparison of Arterial Spin-Labeling Techniques and Dynamic Susceptibility-Weighted Contrast-Enhanced MRI in Perfusion Imaging of Normal Brain Tissue

OBJECTIVES

To evaluate relative cerebral blood flow (rCBF) in normal brain tissue using arterial spin-labeling (ASL) methods and first-pass dynamic susceptibility-weighted contrast-enhanced (DSC) magnetic resonance imaging (MRI).

METHODS

Sixty-two patients with brain metastases were examined on a 1.5 T-system up to 6 times during routine follow-up after stereotactic radiosurgery. Perfusion values in normal gray and white matter were measured using the ASL techniques ITS-FAIR in 38 patients, Q2TIPS in 62 patients, and the first-pass DSC echo-planar (EPI) MRI after bolus administration of gadopentetate dimeglumine in 42 patients. Precision of the ASL sequences was tested in follow-up examinations in 10 healthy volunteers.

RESULTS

Perfusion values in normal brain tissue obtained by all sequences correlated well by calculating Pearson's correlation coefficients (P < 0.0001) and remained unchanged after stereotactic radiosurgery as shown by analysis of variance (P > 0.05).

CONCLUSION

Both ASL and DSC EPI MRI yield highly comparable perfusion values in normal brain tissue.

Measuring the effects of indomethacin on changes in cerebral oxidative metabolism and cerebral blood flow during sensorimotor activation

The work presented here uses combined blood oxygenation level-dependent (BOLD) and arterial spin tagging (AST) approaches to study the effect of indomethacin on cerebral blood flow (CBF) and oxygen consumption (CMRO(2)) increases during motor activation. While indomethacin reduced the CBF increase during activation, it did not significantly affect the CMRO(2) increase during activation. The ratio of the activation-induced CBF increase in the presence and absence of indomethacin was 0.54 +/- 0.08 (+/-SEM, n = 8, P < 0.001), while the ratio of the CMRO(2) increase in the presence and absence of the drug was 1.02 +/- 0.08 (+/-SEM, N = 8, ns). Potential difficulties in estimating CMRO(2) changes from combined BOLD/AST data are discussed.

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.

UNFOLD using a temporal subtraction and spectral energy comparison technique

In dynamic MRI, several methods have been demonstrated to increase acquisition speed by decreasing the number of sequential phase encodings. The UNFOLD technique interleaves the measurements of k-space, reconstructs aliased images from each k-space interleaf, and applies a temporal low-pass filter to obtain the nonaliased images. However, low-pass filter resolution of the nonaliased images fails if there is overlap between the spatially aliased temporal spectra. In this study a subtraction method was used to remove the static portion of the image. The aliased and nonaliased dynamic portions are then resolved by comparing the temporal energy of bands in the power spectrum. This method was combined with the 3D 2 x 2 UNFOLD (a factor of 2 interleaves in two directions) technique. The combination resulted in a factor of 4 improvement in acquisition speed. Application of this method to a time-resolved, contrast-enhanced flow phantom study is presented.

Dopaminergic modulation of cortical function in patients with Parkinson's disease

Patients with idiopathic Parkinson's disease suffer not only from classic motor symptoms, but from deficits in cognitive function, primarily those subserved by the prefrontal cortex as well. The aim of the current study was to investigate the modulatory effects of dopaminergic therapy on neural systems subserving working memory and motor function in patients with Parkinson's disease. Ten patients with stage I and II Parkinson's disease were studied with functional magnetic resonance imaging, during a relatively hypodopaminergic state (ie, 12 hours after a last dose of dopamimetic treatment), and again during a dopamine-replete state. Functional magnetic resonance imaging was performed under three conditions: a working memory task, a cued sensorimotor task and rest. Consistent with prior data, the cortical motor regions activated during the motor task showed greater activation during the dopamine-replete state; however, the cortical regions subserving working memory displayed greater activation during the hypodopaminergic state. Interestingly, the increase in cortical activation during the working memory task in the hypodopaminergic state positively correlated with errors in task performance, and the increased activation in the cortical motor regions during the dopamine-replete state was positively correlated with improvement in motor function. These results support evidence from basic research that dopamine modulates cortical networks subserving working memory and motor function via two distinct mechanisms: nigrostriatal projections facilitate motor function indirectly via thalamic projections to motor cortices, whereas the mesocortical dopaminergic system facilitates working memory function via direct inputs to prefrontal cortex. The results are also consistent with evidence that the hypodopaminergic state is associated with decreased efficiency of prefrontal cortical information processing and that dopaminergic therapy improves the physiological efficiency of this region.

Effect of COMT Val108/158 Met genotype on frontal lobe function and risk for schizophrenia

Abnormalities of prefrontal cortical function are prominent features of schizophrenia and have been associated with genetic risk, suggesting that susceptibility genes for schizophrenia may impact on the molecular mechanisms of prefrontal function. A potential susceptibility mechanism involves regulation of prefrontal dopamine, which modulates the response of prefrontal neurons during working memory. We examined the relationship of a common functional polymorphism (Val(108/158) Met) in the catechol-O-methyltransferase (COMT) gene, which accounts for a 4-fold variation in enzyme activity and dopamine catabolism, with both prefrontally mediated cognition and prefrontal cortical physiology. In 175 patients with schizophrenia, 219 unaffected siblings, and 55 controls, COMT genotype was related in allele dosage fashion to performance on the Wisconsin Card Sorting Test of executive cognition and explained 4% of variance (P = 0.001) in frequency of perseverative errors. Consistent with other evidence that dopamine enhances prefrontal neuronal function, the load of the low-activity Met allele predicted enhanced cognitive performance. We then examined the effect of COMT genotype on prefrontal physiology during a working memory task in three separate subgroups (n = 11-16) assayed with functional MRI. Met allele load consistently predicted a more efficient physiological response in prefrontal cortex. Finally, in a family-based association analysis of 104 trios, we found a significant increase in transmission of the Val allele to the schizophrenic offspring. These data suggest that the COMT Val allele, because it increases prefrontal dopamine catabolism, impairs prefrontal cognition and physiology, and by this mechanism slightly increases risk for schizophrenia.

Temporal frequency analysis of dynamic MRI techniques

Dynamic imaging strategies often involve updating certain areas of k-space (i.e., the low spatial frequencies) more frequently than others. However, important dynamic signal changes may occur anywhere in k-space. In this study, a dynamic k-space sampling analysis method was developed to determine the energy error associated with specific dynamic sampling strategies. The method uses the temporal power spectrum of k-space signals to determine the level and k-space locations of sampling errors. The proposed method was used to compare two dynamic sampling strategies (full sequential and keyhole) for a dynamic first-pass bolus simulation and a continuous heart imaging study. The error analysis agreed well with the errors in the reconstructed images. The technique can be used to determine the minimum sampling frequency for any location in the k-space, and may ultimately be used to optimize dynamic sampling strategies. Magn Reson Med 45:550-556, 2001.

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.

H(2)(15)O PET validation of steady-state arterial spin tagging cerebral blood flow measurements in humans

Steady-state arterial spin tagging approaches can provide quantitative images of CBF, but have not been validated in humans. The work presented here compared CBF values measured using steady-state arterial spin tagging with CBF values measured in the same group of human subjects using the H(2)(15)O IV bolus PET method. Blood flow values determined by H(2)(15)O PET were corrected for the known effects of incomplete extraction of water across the blood brain barrier. For a cortical strip ROI, blood flow values determined using arterial spin tagging (64+/-12 cc/100 g/min) were not statistically different from corrected blood flow values determined using H(2)(15)O PET (67+/-13 cc/100 g/min). However, for a central white matter ROI, blood flow values determined using arterial spin tagging were significantly underestimated compared to corrected blood flow values determined using H(2)(15)O PET. This underestimation could be caused by an underestimation of the arterial transit time for white matter regions.

Noise reduction in 3D perfusion imaging by attenuating the static signal in arterial spin tagging (ASSIST)

Phase-encoded multishot SPIRAL approaches were used to acquire true 3D cerebral blood flow images of the human head using arterial spin tagging approaches. Multiple-inversion background suppression techniques, which suppress phase noise due to interacquisition fluctuations in the static magnetic field, reduced the temporal standard deviation of true 3D delta M images acquired using arterial spin tagging approaches by approximately 50%. Background suppressed arterial spin tagging (ASSIST) approaches were used to obtain high-resolution isotropic true 3D cerebral blood flow images, and to obtain true 3D activation images during cognitive (working memory) tasks. Magn Reson Med 44:92-100, 2000. Published 2000 Wiley-Liss, Inc.

Unaliasing by fourier-encoding the overlaps using the temporal dimension (UNFOLD), applied to cardiac imaging and fMRI

In several applications, MRI is used to monitor the time behavior of the signal in an organ of interest; e.g., signal evolution because of physiological motion, activation, or contrast-agent accumulation. Dynamic applications involve acquiring data in a k-t space, which contains both temporal and spatial information. It is shown here that in some dynamic applications, the t axis of k-t space is not densely filled with information. A method is introduced that can transfer information from the k axes to the t axis, allowing a denser, smaller k-t space to be acquired, and leading to significant reductions in the acquisition time of the temporal frames. Results are presented for cardiac-triggered imaging and functional MRI (fMRI), and are compared with data obtained in a conventional way. The temporal resolution was increased by nearly a factor of two in the cardiac-triggered study, and by as much as a factor of eight in the fMRI study. This increase allowed the acquisition of fMRI activation maps, even when the acquisition time for a single full time frame was actually longer than the paradigm cycle period itself. The new method can be used to significantly reduce the acquisition time of the individual temporal frames in certain dynamic studies. This can be used, for example, to increase the temporal or spatial resolution, increase the spatial coverage, decrease the total imaging time, or alter sequence parameters e.g., repetition time (TR) and echo time (TE) and thereby alter contrast. Magn Reson Med 42:813-828, 1999.

Simple but reliable solutions for spiral MRI gradient design

The efficiency of gradient design in MRI is limited by the simple fact that the gradient coil current and slew rate cannot exceed hardware threshold values. In spiral MRI, which requires gradients to be very rapidly switched between positive and negative values, minimization of the acquisition time is achieved by maintaining the current and slew rate as high as possible during the entire measurement. Since the current and slew rate compete against each other, an efficient gradient design consists of two parts in which current and slew rate are pushed alternatively to their limits. Values for these types of gradients can be obtained by solving numerically the equation of motion for the spiral trajectory. This paper shows that simple but reasonable mathematical approximations deliver reliable analytical solutions. Images obtained using these analytical solutions do not show evident distortions when compared with images obtained with numerical solutions.

Quantitation of regional cerebral blood flow increases during motor activation: A multislice, steady-state, arterial spin tagging study

Steady-state arterial spin tagging approaches were used to construct multislice images of relative cerebral blood flow changes during finger-tapping tasks. Statistically significant increases in cerebral blood flow were observed in primary sensorimotor cortex in all seven subjects. The mean volume of the activated region in the contralateral primary sensorimotor cortex was 0.9 cm(3), and the mean increase in cerebral blood flow in the activated area was 54% +/- 11%. Although the extended spatial coverage is advantageous for activation studies, the intrinsic sensitivity of the multislice approach is smaller than the intrinsic sensitivity of the single-slice, arterial spin tagging approach. Magn Reson Med 42:404-407, 1999. Published 1999 Wiley-Liss, Inc.

3D z-shim method for reduction of susceptibility effects in BOLD fMRI

Susceptibility-induced magnetic field gradients (SFGs) perpendicular to the slice plane often result in signal dropout in blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) experiments. Two-dimensional (2D) z-shim methods reduce these effects by acquiring multiple images with different slice refocusing gradient areas. In this work a 3D z-shim method is introduced as a more efficient alternative. The technique augments the k-space data for a conventional 3D phase encoding acquisition with N additional lines that extend the k(z) coverage sufficiently to sample k-space fully in regions with SFGs. Multiple subsets of these data are reconstructed using a sliding window that provides N +1 z-shim images. Fewer total acquisitions are required than with the 2D method for the same coverage, and finer z-shim steps are obtained. The technique is demonstrated with a motor task using intentionally introduced SFGs and compared with the 2D method. The results confirm increased BOLD SNR and activation with the new method in good agreement with theory. Magn Reson Med 42:290-299, 1999.

Imaging of acute cerebral ischemia

Until recently, there was no efficacious treatment for acute cerebral ischemia. As a result, the role of neuroimaging and the radiologist was peripheral in the diagnosis and management of this disease. The demonstration of efficacy using thrombolysis has redefined this role, with the success of intervention becoming increasingly dependent on timely imaging and accurate interpretation. The potential benefits of intervention have only begun to be realized. In this State-of-the-Art review of imaging of acute stroke, the role of imaging in the current and future management of stroke is presented. The role of computed tomography is emphasized in that it is currently the most utilized technique, and its value has been demonstrated in prospective clinical trials. Magnetic resonance techniques are equally emphasized in that they have the potential to provide a single modality evaluation of tissue viability and vessel patency in an increasingly rapid evaluation.

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.

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.

Organization of the human motor system as studied by functional magnetic resonance imaging

Blood oxygenation level dependent functional magnetic resonance imaging (BOLD fMRI), because of its superior resolution and unlimited repeatability, can be particularly useful in studying functional aspects of the human motor system, especially plasticity, and somatotopic and temporal organization. In this survey, while describing studies that have reliably used BOLD fMRI to examine these aspects of the motor system, we also discuss studies that investigate the neural substrates underlying motor skill acquisition, motor imagery, production of motor sequences; effect of rate and force of movement on brain activation and hemispheric control of motor function. In the clinical realm, in addition to the presurgical evaluation of neurosurgical patients, BOLD fMRI has been used to explore the mechanisms underlying motor abnormalities in patients with neuropsychiatric disorders and the mechanisms underlying reorganization or plasticity of the motor system following a cerebral insult.

3D spiral cardiac/respiratory ordered fMRI data acquisition at 3 Tesla

Three-dimensional (3D), multi-shot functional magnetic resonance imaging (fMRI) data acquisitions are desirable because of higher resolution and reduced susceptibility artifacts, due to shorter readouts and thinner slices. However, 3D multi-shot techniques are more susceptible to physiological noise, which can increase inter-image variance and lead to inaccurate assessment of activation. This work presents a 3D spiral fMRI data acquisition method at 3 T in which the acquisition of views was ordered to match the phase of either the respiratory or the cardiac cycle. For the acquisition timing parameters used in this work, cardiac ordering was found to reduce inter-image variance by 19%. Cardiac ordered data acquisitions showed the same reduction in variance as sequentially ordered data with cardiac contributions estimated and removed using an externally acquired reference prior to reconstruction. Respiratory ordering showed no reduction in fluctuation noise due to poor alignment of views to the respiratory phase.

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.

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.

Functional burst imaging

A quiet magnetic resonance (MR) imaging technique for detecting changes in cerebral activity functions is presented. This single-shot method, functional Burst imaging (FBI), combines elements of Burst imaging with an offset technique known as asymmetric spin echo (ASE). The FBI sequence has the unique feature of being nearly silent, because of the low number of gradient switching steps involved. Furthermore, this approach has the key advantage that it can be implemented on conventional MR systems. Established auditory and visual paradigms were used to evaluate whether FBI can detect changes in cerebral activity using a 1.5 Tesla MR system. In a second set of experiments, the FBI technique was used to evaluate cerebral activity changes during different sleep stages in humans. The results obtained demonstrate that the FBI sequence provides an alternative approach for functional imaging of brain activity in primary and secondary sensory areas of the human brain. Furthermore, in using this quiet MR technique, it was possible to scan continuously during different stages of human sleep without acoustic noise perturbation.

Partial Fourier reconstruction for three-dimensional gradient echo functional MRI: comparison of phase correction methods

Partial Fourier (PF) methods take advantage of data symmetry to allow for either faster image acquisition or increased image resolution. Faster acquisition and increased spatial resolution are advantageous for fMRI because of increased temporal resolution and/or reduced partial volume effects, respectively. Standard PF methods, which use a phase reference obtained from a low resolution image, are adequate for the reconstruction of time-stationary images acquired using either spin echoes or short TE gradient echoes. In fMRI, however, multiple images are acquired using long TE gradient echoes, which introduces possible phase drifts in the fMRI data and high spatial frequencies in the phase reference. This work investigates several techniques developed to reconstruct fMRI data obtained with PF acquisitions. PF methods that account for both high-frequency spatial variations and time-dependent drifts in the phase reference are discussed and are quantitatively evaluated using receiver operator characteristic curve analysis.

Quantitation of regional cerebral blood flow increases in prefrontal cortex during a working memory task: a steady-state arterial spin-tagging study

Steady-state arterial spin-tagging MRI approaches were used to quantitate regional cerebral blood flow increases in prefrontal cortex during a working memory ("two-back") task in six normal subjects. Statistically significant increases in cerebral blood flow in prefrontal cortex were observed in all six subjects: the average increase in cerebral blood flow in activated prefrontal cortex regions was 22 +/- 5 cc/100 g/min (23 +/- 7%). The results demonstrate that spin-tagging approaches can be used to follow focal activation in prefrontal cortex during cognitive tasks.

Three-dimensional spiral fMRI technique: a comparison with 2D spiral acquisition

A comparison study was performed of 3D and 2D spiral k-space fMRI techniques using BOLD contrast with a 4.25-min finger-tapping task paradigm. The 3D sequence uses the conventional 2D spiral technique in conjunction with Fourier phase encoding in the slice select direction. Characteristics that were compared included image appearances, signal-to-noise ratio (SNR), fluctuation noise, functional contrast-to-noise ratio (fCNR), detected activation areas, inflow effect, and large vessel involvement. The results showed that, with constant total scan time, the 3D spiral acquisition has higher SNR and fCNR despite its slightly higher fluctuation noise. The capability for thin-slice incoherent averaging is a further advantage of 3D over 2D. With these advantages as well as its intrinsic ability for contiguous slices and image reformatting, the 3D spiral fMRI technique may be superior to its 2D counterpart.

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

To evaluate the sensitivity of fast, gradient-echo MR scan techniques in their ability to detect blood oxygenation level dependent (BOLD) signal changes in task activation studies, three dedicated fast scan techniques, each with whole-brain coverage, were compared during a 3-min finger tapping paradigm on nine normal volunteers on a clinical 1.5 T scanner. Multislice (2D) single-shot spiral, 3D spiral, and multislice (2D) single-shot EPI scan techniques were done with similar temporal and spatial resolutions on each of the volunteers in random order. After image registration and statistical analysis, the sensitivity to detect activation was evaluated for the techniques by calculating t scores and number of activated voxels in predetermined regions of interest, including the contralateral primary sensorimotor cortex, the premotor region, the parietal region, the supplementary motor area, and the ipsilateral cerebellum. Baseline images acquired with the three techniques were qualitatively comparable and had a similar effective spatial resolution of around 5 x 5 x 5 mm3, as determined from autocorrelation analysis. The anatomical coverage was somewhat reduced (4 less slices per volume) with EPI at the identical temporal resolution of 1.76 s for all techniques. The use of multislice 2D spiral scan for motor cortex fMRI experiments provided for a superior overall temporal stability, and an increased sensitivity compared with multislice 2D EPI, and 3D spiral scan. The difference in sensitivity between multislice 2D spiral and EPI scans was small, in particular in the case of a ramp-sampled version of EPI. The difference in performance is attributed mainly to the difference in scan-to-scan stability.

Fast spiral magnetic resonance imaging with trapezoidal gradients

A modified spiral imaging technique is presented, in which the conventional sinusoidal gradient waveforms are replaced by trapezoidal ones. In addition to allowing a reduced data acquisition time, the new waveforms circumvent specific hardware restrictions on the minimum scan repetition time.

Quantitation of regional cerebral blood flow increases during motor activation: a steady-state arterial spin tagging study

Steady-state arterial spin tagging MRI approaches were used to quantitate regional cerebral blood flow increases during finger tapping tasks in seven normal subjects. Statistically significant increases in cerebral blood flow were observed in the contralateral primary sensorimotor cortex in all seven subjects and in the supplementary motor area in five subjects. The intrinsic spatial resolution of the cerebral blood flow images was approximately 4 mm. If no spatial filtering was applied, the average increase in cerebral blood flow in the activated primary sensorimotor cortex was 60 +/- 10 cc/100 g/min (91 +/- 32%). If the images were filtered to a spatial resolution of 15 mm, the average increase in cerebral blood flow in the activated primary sensorimotor cortex was 23 +/- 7 cc/100 g/min (42 +/- 15%), in agreement with previously reported 133Xe and PET results.