Ming Long Wu

Quantitative susceptibility mapping of human brain at 3T: a multisite reproducibility study.

BACKGROUND AND PURPOSE: Quantitative susceptibility mapping of the human brain has demonstrated strong potential in examining iron deposition, which may help in investigating possible brain pathology. This study assesses the reproducibility of quantitative susceptibility mapping across different imaging sites. MATERIALS AND METHODS: In this study, the susceptibility values of 5 regions of interest in the human brain were measured on 9 healthy subjects following calibration by using phantom experiments. Each of the subjects was imaged 5 times on 1 scanner with the same procedure repeated on 3 different 3T systems so that both within-site and cross-site quantitative susceptibility mapping precision levels could be assessed. Two quantitative susceptibility mapping algorithms, similar in principle, one by using iterative regularization (iterative quantitative susceptibility mapping) and the other with analytic optimal solutions (deterministic quantitative susceptibility mapping), were implemented, and their performances were compared. RESULTS: Results show that while deterministic quantitative susceptibility mapping had nearly 700 times faster computation speed, residual streaking artifacts seem to be more prominent compared with iterative quantitative susceptibility mapping. With quantitative susceptibility mapping, the putamen, globus pallidus, and caudate nucleus showed smaller imprecision on the order of 0.005 ppm, whereas the red nucleus and substantia nigra, closer to the skull base, had a somewhat larger imprecision of approximately 0.01 ppm. Cross-site errors were not significantly larger than within-site errors. Possible sources of estimation errors are discussed. CONCLUSIONS: The reproducibility of quantitative susceptibility mapping in the human brain in vivo is regionally dependent, and the precision levels achieved with quantitative susceptibility mapping should allow longitudinal and multisite studies such as aging-related changes in brain tissue magnetic susceptibility.

Frequency stabilization using infinite impulse response filtering for SSFP fMRI at 3T.

The steady‐state free precession (SSFP) method has been shown to exhibit strong potential for distortion‐free functional magnetic resonance imaging (fMRI). One major challenge of SSFP fMRI is that the frequency band corresponding to the highest functional sensitivity is extremely narrow, leading to substantial loss of functional contrast in the presence of magnetic field drifts. In this study we propose a frequency stabilization scheme whereby an RF pulse with small flip angle is applied before each image scan, and the initial phase of the free induction decay (FID) signals is extracted to reflect temporal field drifts. A simple infinite impulse response (IIR) filter is further employed to obtain a low‐pass‐filtered estimate of the central reference frequency for the upcoming scan. Experimental results suggest that the proposed scheme can stabilize the frequency settings in accordance with field drifts, with oscillation amplitudes of <0.5 Hz. Phantom studies showed that both slow drifts and fast fluctuations were prominently reduced, resulting in less than 5% signal variations. Visual fMRI at submillimeter in‐plane resolution further demonstrated 15% activation signals that were nicely registered in the microvessels within the sulci. It is concluded that the IIR‐filtered frequency stabilization is an effective technique for achieving reliable SSFP fMR images at high field strengths. Magn Reson Med 57:369–379, 2007.

Are the Local Blood Oxygen Level–Dependent (BOLD) Signals Caused by Neural Stimulation Response Dependent on Global BOLD Signals Induced by Hypercapnia in the Functional MR Imaging Experiment? Experiments of Long-Duration Hypercapnia and Multilevel Carbon Dioxide Concentration

BACKGROUND AND PURPOSE: The relationship between the local blood oxygen level–dependent (BOLD) signals caused by neural stimulation (fBOLD) and the global BOLD signals induced by hypercapnia (hBOLD) has not been fully investigated. In this study, we examine whether fBOLD is modulated by hBOLD signals, by means of experiments using a relatively wide range of inhaled carbon dioxide (CO2) for a long duration of 5 minutes. MATERIALS AND METHODS: Ten healthy volunteers were recruited, each undergoing 6 separate experiments by inhaling gas mixtures with different fractions of CO2 (room air, 3%–7%). Each experiment contained 3 phases, prehypercapnic, hypercapnic, and posthypercapnic, during which boxcar visual stimulus was given. The local fBOLD signals were measured from areas showing activation patterns highly correlated with the visual stimulus paradigm, whereas the global hBOLD signals were measured from areas showing no visual activations. Percentage changes in fBOLD during transient-state hypercapnia and steady-state hypercapnia were both investigated in response to varying degrees of hypercapnic perturbations. RESULTS: The hBOLD signals increased with increase of inhaled CO2 fractions. The duration for the hBOLD signals to reach steady state prolonged with increase of inhaled CO2 fractions. Normalized fBOLD ratio was inversely related to the inhaled CO2 during steady-state hypercapnia but showed positive association with hBOLD during transient-state hypercapnia. CONCLUSION: Our study concludes that the steady-state fBOLD signal intensity is dependent on and inversely related to the hBOLD signals. Previous reports documenting independent and additive relationships between hBOLD and fBOLD may likely be due to transient-state observations.

Tensor deflection (TEND) tractography with adaptive subvoxel stepping

To develop an adaptive subvoxel stepping scheme, as an adjunct to tensor deflection (TEND) tractography, that automatically adjusts the stepping size by considering the tensor linearity to properly trace fiber bundles in regions with different degrees of tensor anisotropy.

Application of model-free analysis in the MR assessment of pulmonary perfusion dynamics

Dynamic contrast‐enhanced (DCE) MRI has been used to quantitatively evaluate pulmonary perfusion based on the assumption of a gamma‐variate function and an arterial input function (AIF) for deconvolution. However, these assumptions may be too simplistic and may not be valid in pathological conditions, especially in patients with complex inflow patterns (such as in congenital heart disease). Exploratory data analysis methods make minimal assumptions on the data and could overcome these pitfalls. In this work, two temporal clustering methods—Kohonen clustering network (KCN) and Fuzzy C‐Means (FCM)—were concatenated to identify pixel time‐course patterns. The results from seven normal volunteers show that this technique is superior for discriminating vessels and compartments in the pulmonary circulation. Patient studies with five cases of acquired or congenital pulmonary perfusion disorders demonstrate that pathologies can be highlighted in a concise map that combines information of the mean transit time (MTT) and pulmonary blood volume (PBV). The method was found to provide greater insight into the perfusion dynamics that might be overlooked by current model‐based analyses, and may serve as a basis for optimal hemodynamic quantitative modeling of the interrogated perfusion compartments. Magn Reson Med 54:299–308, 2005.

Effects of RF profile on precision of quantitative T2 mapping using dual-echo steady-state acquisition

The dual echo steady-state (DESS) sequence has been shown successful in achieving fast T2 mapping with good precision. Under-estimation of T2, however, becomes increasingly prominent as the flip angle decreases. In 3D DESS imaging, therefore, the derived T2 values would become a function of the slice location in the presence of non-ideal slice profile of the excitation RF pulse. Furthermore, the pattern of slice-dependent variation in T2 estimates is dependent on the RF pulse waveform. Multi-slice 2D DESS imaging provides better inter-slice consistency, but the signal intensity is subject to integrated effects of within-slice distribution of the actual flip angle. Consequently, T2 measured using 2D DESS is prone to inaccuracy even at the designated flip angle of 90°. In this study, both phantom and human experiments demonstrate the above phenomena in good agreement with model prediction.

High spatial resolution brain functional MRI using submillimeter balanced steady-state free precession acquisition

PURPOSE One of the technical advantages of functional magnetic resonance imaging (fMRI) is its precise localization of changes from neuronal activities. While current practice of fMRI acquisition at voxel size around 3 × 3 × 3 mm(3) achieves satisfactory results in studies of basic brain functions, higher spatial resolution is required in order to resolve finer cortical structures. This study investigated spatial resolution effects on brain fMRI experiments using balanced steady-state free precession (bSSFP) imaging with 0.37 mm(3) voxel volume at 3.0 T. METHODS In fMRI experiments, full and unilateral visual field 5 Hz flashing checkerboard stimulations were given to healthy subjects. The bSSFP imaging experiments were performed at three different frequency offsets to widen the coverage, with functional activations in the primary visual cortex analyzed using the general linear model. Variations of the spatial resolution were achieved by removing outer k-space data components. RESULTS Results show that a reduction in voxel volume from 3.44 × 3.44 × 2 mm(3) to 0.43 × 0.43 × 2 mm(3) has resulted in an increase of the functional activation signals from (7.7 ± 1.7)% to (20.9 ± 2.0)% at 3.0 T, despite of the threefold SNR decreases in the original images, leading to nearly invariant functional contrast-to-noise ratios (fCNR) even at high spatial resolution. Activation signals aligning nicely with gray matter sulci at high spatial resolution would, on the other hand, have possibly been mistaken as noise at low spatial resolution. CONCLUSIONS It is concluded that the bSSFP sequence is a plausible technique for fMRI investigations at submillimeter voxel widths without compromising fCNR. The reduction of partial volume averaging with nonactivated brain tissues to retain fCNR is uniquely suitable for high spatial resolution applications such as the resolving of columnar organization in the brain.

Efficient imaging of midbrain nuclei using inverse double-echo steady-state acquisition

PURPOSE Imaging of midbrain nuclei using T2- or T2*-weighted MRI often entails long echo time, leading to long scan time. In this study, an inverse double-echo steady-state (iDESS) technique is proposed for efficiently depicting midbrain nuclei. METHODS Thirteen healthy subjects participated in this study. iDESS was performed along with two sets of T2*-weighted spoiled gradient-echo images (SPGR1, with scan time identical to iDESS and SPGR2, using clinical scanning parameters as a reference standard) for comparison. Generation of iDESS composite images combining two echo signals was optimized for maximal contrast-to-noise ratio (CNR) between the red nuclei and surrounding tissues. Signal-to-noise ratios (SNRs) were calculated from the occipital lobe. Comparison was also made using phase-enhanced images as in standard susceptibility-weighted imaging (SWI). RESULTS The iDESS images present significantly higher SNR efficiency (171.3) than SPGR1 (158.7, p = 0.013) and SPGR2 (95.5, p < 10(-8)). iDESS CNR efficiency (19.2) is also significantly greater than SPGR1 (6.9, p < 10(-6)) and SPGR2 (14.3, p = 0.0016). Compared with DESS, iDESS provides further advantage on enhanced phase information and hence improved contrast on SWI-processed images. CONCLUSIONS iDESS efficiently depicts midbrain nuclei with improved CNR efficiency, increased SNR efficiency, and reduced scan time and is less prone to susceptibility signal loss from air-tissue interfaces.

De-aliasing for signal restoration in Propeller MR imaging

PURPOSE Objects falling outside of the true elliptical field-of-view (FOV) in Propeller imaging show unique aliasing artifacts. This study proposes a de-aliasing approach to restore the signal intensities in Propeller images without extra data acquisition. MATERIALS AND METHODS Computer simulation was performed on the Shepp-Logan head phantom deliberately placed obliquely to examine the signal aliasing. In addition, phantom and human imaging experiments were performed using Propeller imaging with various readouts on a 3.0 Tesla MR scanner. De-aliasing using the proposed method was then performed, with the first low-resolution single-blade image used to find out the aliasing patterns in all the single-blade images, followed by standard Propeller reconstruction. The Propeller images without and with de-aliasing were compared. RESULTS Computer simulations showed signal loss at the image corners along with aliasing artifacts distributed along directions corresponding to the rotational blades, consistent with clinical observations. The proposed de-aliasing operation successfully restored the correct images in both phantom and human experiments. CONCLUSION The de-aliasing operation is an effective adjunct to Propeller MR image reconstruction for retrospective restoration of aliased signals.