Alternating repetition time balanced steady state free precession

Submillimeter balanced SSFP BOLD-functional MRI accelerated with 3D stack-of-spirals at 9.4 T

PURPOSE

This work aims to improve the speed of balanced SSFP (bSSFP) acquisition with segmented 3D stack-of-spirals for functional brain studies at ultrahigh field.

METHODS

Functional experiments were performed with an accelerated 3D stack-of-spirals sequence with water excitation for fat suppression. The resulting data were reconstructed using an iterative algorithm with corrections for system imperfections such as trajectory deviations and B0 inhomogeneity. In the first set of experiments, we evaluated the signal change and stability with respect to echo and TR for a full-field checkerboard stimulus. To demonstrate the high spatio-temporal resolution of the developed method, the results of three optimized protocols at submillimeter resolution (0.6-mm isotropic and 0.8-mm isotropic) and at 1.2 mm isotropic resolution for whole-brain coverage were shown.

RESULTS

Water excitation and the model-based iterative reconstruction improved image quality. The BOLD-related signal changes increased with longer TE and longer TR. We observed an increase in thermal noise performance at lower TE and higher TR. However, signal stability deteriorates at higher TE and TR. Therefore, optimized protocols used shorter TE and moderately long TR to maximize the sensitivity and speed. Reproducible activations were detected along the gray-matter gyri in the submillimeter protocols with a median signal change of approximately 4% across subjects.

CONCLUSIONS

Three-dimensional stack-of-spirals enables passband balanced SSFP functional imaging at a much higher spatial and temporal scale, compared with conventional spoiled gradient-echo train sequences.

Optimization of flip angle and radiofrequency pulse phase to maximize steady-state magnetization in three-dimensional missing pulse steady-state free precession

Missing pulse (MP) steady-state free precession (SSFP) is a magnetic resonance imaging (MRI) pulse sequence that is highly tolerant to the magnetic field inhomogeneity. In this study, optimal flip angle and radiofrequency (RF) phase scheduling in three-dimensional (3D) MP-SSFP is introduced to maximize the steady-state magnetization while keeping broadband excitation to cover widely distributed frequencies generated by inhomogeneous magnetic fields. Numerical optimization based on extended phase graph (EPG) simulation was performed to maximize the MP-SSFP steady-state magnetization. To limit the specific absorption rate (SAR) associated with the broadband excitation in 3D MP-SSFP, SAR constraint was introduced in the numerical optimization. Optimized flip angle and RF phase settings were experimentally tested by introducing a linear inhomogeneous magnetic field in a range of 10-20 mT/m and using a phantom with known T1/T2 relaxation and diffusion parameters at 3 T. The experimental results were validated through comparisons with EPG simulation. Image contrasts and molecular diffusion effects were investigated in in vivo human brain imaging with 3D MP-SSFP with the optimal flip angle and RF phase settings. In the phantom measurements, the optimal flip angle and RF phase settings improved the MP-SSFP steady-state magnetization/signal-to-noise ratio by up to 41% under the fixed SAR conditions, which matched well with EPG simulation results. In vivo brain imaging with the optimal RF pulse settings provided T2-like image contrasts. Diffusion effects were relatively minor with the linear inhomogeneous field of 10-20 mT/m for white and gray matter, but cerebrospinal fluid showed conspicuous signal intensity attenuation as the linear inhomogeneous field increased. Numerical optimization achieved significant improvement in the steady-state magnetization in MP-SSFP compared with the RF pulse settings used in previous studies. The proposed flip angle and RF phase optimization is promising to improve 3D MP-SSFP image quality for MRI in inhomogeneous magnetic fields.

High-resolution, respiratory-resolved coronary MRA using a Phyllotaxis-reordered variable-density 3D cones trajectory

PURPOSE

To develop a respiratory-resolved motion-compensation method for free-breathing, high-resolution coronary magnetic resonance angiography (CMRA) using a 3D cones trajectory.

METHODS

To achieve respiratory-resolved 0.98 mm resolution images in a clinically relevant scan time, we undersample the imaging data with a variable-density 3D cones trajectory. For retrospective motion compensation, translational estimates from 3D image-based navigators (3D iNAVs) are used to bin the imaging data into four phases from end-expiration to end-inspiration. To ensure pseudo-random undersampling within each respiratory phase, we devise a phyllotaxis readout ordering scheme mindful of eddy current artifacts in steady state free precession imaging. Following binning, residual 3D translational motion within each phase is computed using the 3D iNAVs and corrected for in the imaging data. The noise-like aliasing characteristic of the combined phyllotaxis and cones sampling pattern is leveraged in a compressed sensing reconstruction with spatial and temporal regularization to reduce aliasing in each of the respiratory phases.

RESULTS

In initial studies of six subjects, respiratory motion compensation using the proposed method yields improved image quality compared to non-respiratory-resolved approaches with no motion correction and with 3D translational correction. Qualitative assessment by two cardiologists and quantitative evaluation with the image edge profile acutance metric indicate the superior sharpness of coronary segments reconstructed with the proposed method (P < 0.01).

CONCLUSION

We have demonstrated a new method for free-breathing, high-resolution CMRA based on a variable-density 3D cones trajectory with modified phyllotaxis ordering and respiratory-resolved motion compensation with 3D iNAVs.

Magnetic resonance imaging with gradient sound respiration guide

Respiratory motion management is crucial for high-resolution MRI of the heart, lung, liver and kidney. In this article, respiration guide using acoustic sound generated by pulsed gradient waveforms was introduced in the pulmonary ultrashort echo time (UTE) sequence and validated by comparing with retrospective respiratory gating techniques. The validated sound-guided respiration was implemented in non-contrast enhanced renal angiography. In the sound-guided respiration, breathe−in and–out instruction sounds were generated with sinusoidal gradient waveforms with two different frequencies (602 and 321 Hz). Performance of the sound-guided respiration was evaluated by measuring sharpness of the lung-liver interface with a 10–90% rise distance, w_10-90, and compared with three respiratory motion managements in a free-breathing UTE scan: without respiratory gating (w/o gating), 0-dimensional k-space navigator (k-point navigator), and image-based self-gating (Img-SG). The sound-guided respiration was implemented in stack-of-stars balanced steady-state free precession with inversion recovery preparation for renal angiography. No subjects reported any discomfort or inconvenience with the sound-guided respiration in pulmonary or renal MRI scans. The lung-liver interface of the UTE images for sound-guided respiration (w_10-90 = 6.99 ± 2.90 mm), k-point navigator (8.51 ± 2.71 mm), and Img-SG (7.01 ± 2.06 mm) was significantly sharper than that for w/o gating (17.13 ± 2.91 mm; p < 0.0001 for all of sound-guided respiration, k-point navigator and Img-SG). Sharpness of the lung-liver interface was comparable between sound-guided respiration and Img-SG (p = 0.99), but sound-guided respiration achieved better visualization of pulmonary vasculature. Renal angiography with the sound-guided respiration clearly delineated renal, segmental and interlobar arteries. In conclusion, the gradient sound guided respiration can facilitate a consistent diaphragm position in every breath and achieve performance of respiratory motion management comparable to image-based self-gating.

Balanced steady-state free precession MRCP is a robust alternative to respiration-navigated 3D turbo-spin-echo MRCP

Background

Despite synchronization to respiration, respiration-navigated (RN) 3D turbo-spin-echo MRCP is limited by susceptibility to motion artifacts. The aim of this study was to assess the quality of pancreaticobiliary duct visualization of a non-RN MRCP alternative based on balanced steady-state free precession imaging (BSSFP) with overlapping slices compared with RN-MRCP.

Methods

This is a retrospective study on 50 patients without pancreaticobiliary duct disease receiving MRCP at 1.5 T. We performed an intraindividual comparison of coronal RN-MRCP with combined coronal and transverse BSSFP-MRCP. Image quality was scored by 3 readers for 6 pancreaticobiliary duct segments (3 pancreatic, 3 biliary) using a 6-point scale. A segment score of 3 or lower as assessed by at least 2 of 3 readers was defined as insufficient segment visualization. Nonparametric tests and interrater reliability testing were used for statistical analysis.

Results

Overall duct visualization averaged over all readers was scored with 4.5 ± 1.1 for RN-MRCP (pancreatic, 4.1 ± 0.5; biliary, 5.0 ± 0.4) and 4.9 ± 0.9 for combined coronal and transverse BSSFP-MRCP (pancreatic, 4.6 ± 0.6; biliary, 5.1 ± 0.6), respectively (p < 0.001). The number of segments visualized insufficiently was 81/300 for RN-MRCP and 43/300 for BSSFP-MRCP (p < 0.001). Segments visualized insufficiently only in RN-MRCP had a mean score of 4.4 ± 0.8 in BSSFP-MRCP. Overall interrater agreement on superiority of BSSFP-MRCP segment scores over corresponding RN-MRCP was 0.70. Mean acquisition time was 98% longer for RN-MRCP (198.0 ± 98.7 s) than for combined coronal and transverse BSSFP-MRCP (100.2 ± 0.4 s).

Conclusions

Non-RN BSSFP-MRCP with overlapping slices is a fast alternative to RN-MRCP, frequently providing sufficient duct visualization when RN-MRCP fails.

Rapid 3D whole-heart cine imaging using golden ratio stack of spirals

Three-dimensional cine imaging provides a wealth of information about cardiac anatomy and function, but its use in the clinical environment is limited because data acquisition is very time consuming. In this work, a free-breathing 3D whole-heart cine imaging framework was developed using a time-efficient stack of spirals trajectory and accelerated reconstruction. Two suitable view ordering methods are considered with different spacing between k-space readouts in the partition dimension: uniform and tiny golden ratio based. A simulation study suggested the latter did not present any benefits in terms of similarity to the true image. The proposed method was subsequently tested on 10 prospective subjects and compared with conventional multi-slice breath-hold imaging. Image quality was evaluated using objective and subjective scores and ventricular measurements were compared to assess clinical accuracy. Image quality was lower in the proposed technique than in breath-hold images but good agreement was found in clinically relevant ventricular measurements. In addition, the proposed method was fast to acquire, required minimal planning and provided full anatomical coverage with isotropic resolution.

Combined T2 -preparation and multidimensional outer volume suppression for coronary artery imaging with 3D cones trajectories

PURPOSE

To develop a modular magnetization preparation sequence for combined T2 -preparation and multidimensional outer volume suppression (OVS) for coronary artery imaging.

METHODS

A combined T2 -prepared 1D OVS sequence with fat saturation was defined to contain a 90°-60 180°60 composite nonselective tip-down pulse, two 180°Y hard pulses for refocusing, and a -90° spectral-spatial sinc tip-up pulse. For 2D OVS, 2 modules were concatenated, selective in X and then Y. Bloch simulations predicted robustness of the sequence to B0 and B1 inhomogeneities. The proposed sequence was compared with a T2 -prepared 2D OVS sequence proposed by Luo et al, which uses a spatially selective 2D spiral tip-up. The 2 sequences were compared in phantom studies and in vivo coronary artery imaging studies with a 3D cones trajectory.

RESULTS

Phantom results demonstrated superior OVS for the proposed sequence compared with the Luo sequence. In studies on 15 healthy volunteers, the proposed sequence had superior image edge profile acutance values compared with the Luo sequence for the right (P < .05) and left (P < .05) coronary arteries, suggesting superior vessel sharpness. The proposed sequence also had superior signal-to-noise ratio (P < .05) and passband-to-stopband ratio (P < .05). Reader scores and reader preference indicated superior coronary image quality of the proposed sequence for both the right (P < .05) and left (P < .05) coronary arteries.

CONCLUSION

The proposed sequence with concatenated 1D spatially selective tip-ups and integrated fat saturation has superior image quality and suppression compared with the Luo sequence with 2D spatially selective tip-up.

Reconstruction of undersampled 3D non-Cartesian image-based navigators for coronary MRA using an unrolled deep learning model

PURPOSE

To rapidly reconstruct undersampled 3D non-Cartesian image-based navigators (iNAVs) using an unrolled deep learning (DL) model, enabling nonrigid motion correction in coronary magnetic resonance angiography (CMRA).

METHODS

An end-to-end unrolled network is trained to reconstruct beat-to-beat 3D iNAVs acquired during a CMRA sequence. The unrolled model incorporates a nonuniform FFT operator in TensorFlow to perform the data-consistency operation, and the regularization term is learned by a convolutional neural network (CNN) based on the proximal gradient descent algorithm. The training set includes 6,000 3D iNAVs acquired from 7 different subjects and 11 scans using a variable-density (VD) cones trajectory. For testing, 3D iNAVs from 4 additional subjects are reconstructed using the unrolled model. To validate reconstruction accuracy, global and localized motion estimates from DL model-based 3D iNAVs are compared with those extracted from 3D iNAVs reconstructed with l 1 -ESPIRiT. Then, the high-resolution coronary MRA images motion corrected with autofocusing using the l 1 -ESPIRiT and DL model-based 3D iNAVs are assessed for differences.

RESULTS

3D iNAVs reconstructed using the DL model-based approach and conventional l 1 -ESPIRiT generate similar global and localized motion estimates and provide equivalent coronary image quality. Reconstruction with the unrolled network completes in a fraction of the time compared to CPU and GPU implementations of l 1 -ESPIRiT (20× and 3× speed increases, respectively).

CONCLUSIONS

We have developed a deep neural network architecture to reconstruct undersampled 3D non-Cartesian VD cones iNAVs. Our approach decreases reconstruction time for 3D iNAVs, while preserving the accuracy of nonrigid motion information offered by them for correction.

3D Cartesian fast interrupted steady-state (FISS) imaging

PURPOSE

To enable intrinsic and efficient fat suppression in 3D Cartesian fast interrupted steady-state (FISS) acquisitions.

METHODS

A periodic interruption of the balanced steady-state free precession (bSSFP) readout train (FISS) has been previously proposed for 2D radial imaging. FISS modulates the bSSFP frequency response pattern in terms of shape, width and location of stop band (attenuated transverse magnetization). Depending on the FISS interruption rate, the stop band characteristic can be exploited to suppress the fat spectrum at 3.5 ppm, thus yielding intrinsic fat suppression. For conventional 2D Cartesian sampling, ghosting/aliasing artifacts along phase-encoding direction have been reported. In this work, we propose to extend FISS to 3D Cartesian imaging and report countermeasures for the previously observed ghosting/aliasing artifacts. Key parameters (dummy prepulses, spatial resolution, and interruption rate) are investigated to optimize fat suppression and image quality. FISS behavior is examined using extended phase graph simulations to recommend parametrizations which are validated in phantom and in vivo measurements on a 1.5T MRI scanner for 3 applications: upper thigh angiography, abdominal imaging, and free-running 5D CINE.

RESULTS

Using optimized parameters, 3D Cartesian FISS provides homogeneous and consistent fat suppression for all 3 applications. In upper thigh angiography, vessel structures can be recovered in FISS that are obscured in bSSFP. Fat suppression in free-running cardiac CINE resulted in less fat-related motion aliasing and yielded better image quality.

CONCLUSION

3D Cartesian FISS is feasible and offers homogeneous intrinsic fat suppression for selected imaging parameters without the need for dedicated preparation pulses, making it a promising candidate for free-running fat-suppressed imaging.

Whole-heart coronary MR angiography using a 3D cones phyllotaxis trajectory

PURPOSE

To develop a 3D cones steady-state free precession sequence with improved robustness to respiratory motion while mitigating eddy current artifacts for free-breathing whole-heart coronary magnetic resonance angiography.

METHOD

The proposed sequence collects cone interleaves using a phyllotaxis pattern, which allows for more distributed k-space sampling for each heartbeat compared to a typical sequential collection pattern. A Fibonacci number of segments is chosen to minimize eddy current effects with the trade-off of an increased number of acquisition heartbeats. For verification, phyllotaxis-cones is compared to sequential-cones through simulations, phantom studies, and in vivo coronary scans with 8 subjects using 2D image-based navigators for retrospective motion correction.

RESULTS

Simulated point spread functions and moving phantom results show less coherent motion artifacts for phyllotaxis-cones compared to sequential-cones. Assessment of the right and left coronary arteries using reader scores and the image edge profile acutance vessel sharpness metric indicate superior image quality and sharpness for phyllotaxis-cones.

CONCLUSION

Phyllotaxis 3D cones results in improved qualitative image scores and coronary vessel sharpness for free-breathing whole-heart coronary magnetic resonance angiography compared to standard sequential ordering when using a steady-state free precession sequence.

Frequency-modulated bSSFP for phase-sensitive separation of water and fat

Our study proposes the use of a frequency-modulated acquisition which suppresses banding artefacts in combination with a phase-sensitive water-fat separation algorithm. The performance of the phase-sensitive separation for standard bSSFP, complex sum combination thereof, and frequency-modulated bSSFP were compared in in vivo measurements of the upper and lower legs at 1.5 and 3 T. It is shown, that the standard acquisition suffered from banding artefacts and major swaps between tissues. The dual-acquisition bSSFP could alleviate banding artefacts and only minor swaps occurred, but it comes at the expense of a doubled acquisition. In the frequency-modulated acquisitions all banding artefacts and the associated phase jumps were eliminated and no swaps between tissues occurred. It therefore provides a means to robustly separate water and fat, in one single radial bSSFP scan, using the phase-sensitive approach, even in the presence of high field inhomogeneities.

MRI-based Assessment of 3D Intrafractional Motion of Head and Neck Cancer for Radiation Therapy

PURPOSE

To determine the 3-dimensional (3D) intrafractional motion of head and neck squamous cell carcinoma (HNSCC).

METHODS AND MATERIALS

Dynamic contrast-enhanced magnetic resonance images from 56 patients with HNSCC in the treatment position were analyzed. Dynamic contrast-enhanced magnetic resonance imaging consisted of 3D images acquired every 2.9 seconds for 4 minutes 50 seconds. Intrafractional tumor motion was studied in the 3 minutes 43 seconds of images obtained after initial contrast enhancement. To assess tumor motion, rigid registration (translations only) was performed using a region of interest (ROI) mask around the tumor. The results were compared with bulk body motion from registration to all voxels. Motion was split into systematic motion and random motion. Correlations between the tumor site and random motion were tested. The within-subject coefficient of variation was determined from 8 patients with repeated baseline measures. Random motion was also assessed at the end of the first week (38 patients) and second week (25 patients) of radiation therapy to investigate trends of motion.

RESULTS

Tumors showed irregular occasional rapid motion (eg, swallowing or coughing), periodic intermediate motion (respiration), and slower systematic drifts throughout treatment. For 95% of the patients, displacements due to systematic and random motion were <1.4 mm and <2.1 mm, respectively, 95% of the time. The motion without an ROI mask was significantly (P<.0001, Wilcoxon signed rank test) less than the motion with an ROI mask, indicating that tumors can move independently from the bony anatomy. Tumor motion was significantly (P=.005, Mann-Whitney U test) larger in the hypopharynx and larynx than in the oropharynx. The within-subject coefficient of variation for random motion was 0.33. The average random tumor motion did not increase notably during the first 2 weeks of treatment.

CONCLUSIONS

The 3D intrafractional tumor motion of HNSCC is small, with systematic motion <1.4 mm and random motion <2.1 mm 95% of the time.

Spectrally selective three-dimensional dynamic balanced steady-state free precession for hyperpolarized C-13 metabolic imaging with spectrally selective radiofrequency pulses

PURPOSE

Balanced steady-state free precession (bSSFP) sequences can provide superior signal-to-noise ratio efficiency for hyperpolarized (HP) carbon-13 (¹³ C) magnetic resonance imaging by efficiently utilizing the nonrecoverable magnetization, but managing their spectral response is challenging in the context of metabolic imaging. A new spectrally selective bSSFP sequence was developed for fast imaging of multiple HP ¹³ C metabolites with high spatiotemporal resolution.

THEORY AND METHODS

This novel approach for bSSFP spectral selectivity incorporates optimized short-duration spectrally selective radiofrequency pulses within a bSSFP pulse train and a carefully chosen repetition time to avoid banding artifacts.

RESULTS

The sequence enabled subsecond 3D dynamic spectrally selective imaging of ¹³ C metabolites of copolarized [1-¹³ C]pyruvate and [¹³ C]urea at 2-mm isotropic resolution, with excellent spectral selectivity (∼100:1). The sequence was successfully tested in phantom studies and in vivo studies with normal mice.

CONCLUSION

This sequence is expected to benefit applications requiring dynamic volumetric imaging of metabolically active ¹³ C compounds at high spatiotemporal resolution, including preclinical studies at high field and, potentially, clinical studies. Magn Reson Med 78:963-975, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Nonrigid Motion Correction With 3D Image-Based Navigators for Coronary MR Angiography

PURPOSE

To develop a retrospective nonrigid motion-correction method based on 3D image-based navigators (iNAVs) for free-breathing whole-heart coronary magnetic resonance angiography (MRA).

METHODS

The proposed method detects global rigid-body motion and localized nonrigid motion from 3D iNAVs and compensates them with an autofocusing algorithm. To model the global motion, 3D rotation and translation are estimated from the 3D iNAVs. Two sets of localized nonrigid motions are obtained from deformation fields between 3D iNAVs and reconstructed binned images, respectively. A bank of motion-corrected images is generated and the final image is assembled pixel-by-pixel by selecting the best focused pixel from this bank. In vivo studies with six healthy volunteers were conducted to compare the performance of the proposed method with 3D translational motion correction and no correction.

RESULTS

In vivo studies showed that compared to no correction, 3D translational motion correction and the proposed method increased the vessel sharpness by 13% ± 13% and 19% ± 16%, respectively. Out of 90 vessel segments, 75 segments showed improvement with the proposed method compared to 3D translational correction.

CONCLUSION

We have developed a nonrigid motion-correction method based on 3D iNAVs and an autofocusing algorithm that improves the vessel sharpness of free-breathing whole-heart coronary MRA. Magn Reson Med 77:1884-1893, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

3D image-based navigators for coronary MR angiography

PURPOSE

To develop a method for acquiring whole-heart 3D image-based navigators (iNAVs) with isotropic resolution for tracking and correction of localized motion in coronary magnetic resonance angiography (CMRA).

METHODS

To monitor motion in all regions of the heart during a free-breathing scan, a variable-density cones trajectory was designed to collect a 3D iNAV every heartbeat in 176 ms with 4.4 mm isotropic spatial resolution. The undersampled 3D iNAV data were reconstructed with efficient self-consistent parallel imaging reconstruction (ESPIRiT). 3D translational and nonrigid motion-correction methods using 3D iNAVs were compared to previous translational and nonrigid methods using 2D iNAVs.

RESULTS

Five subjects were scanned with a 3D cones CMRA sequence, accompanied by both 2D and 3D iNAVs. The quality of the right and left anterior descending coronary arteries was assessed on 2D and 3D iNAV-based motion-corrected images using a vessel sharpness metric and qualitative reader scoring. This assessment showed that nonrigid motion correction based on 3D iNAVs produced results that were noninferior to correction based on 2D iNAVs.

CONCLUSION

The ability to acquire isotropic-resolution 3D iNAVs every heartbeat during a CMRA scan was demonstrated. Such iNAVs enabled direct measurement of localized motion for nonrigid motion correction in free-breathing whole-heart CMRA. Magn Reson Med 77:1874-1883, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Enhanced phase-sensitive SSFP reconstruction for fat-water separation in phased-array acquisitions

PURPOSE

To propose and assess a method to improve the reliability of phase-sensitive fat-water separation for phased-array balanced steady-state free precession (bSSFP) acquisitions. Phase-sensitive steady-state free precession (PS-SSFP) is an efficient fat-water separation technique that detects the phase difference between neighboring bands in the bSSFP magnetization profile. However, large spatial variations in the sensitivity profiles of phased-array coils can lead to noisy phase estimates away from the coil centers, compromising tissue classification.

MATERIALS AND METHODS

We first perform region-growing phase correction in individual coil images via unsupervised selection of a fat-voxel seed near the peak of each coil's sensitivity profile. We then use an optimal linear combination of phase-corrected images to segregate fat and water signals. The proposed method was demonstrated on noncontrast-enhanced SSFP angiograms of the thigh, lower leg, and foot acquired at 1.5T using an 8-channel coil. Individual coil PS-SSFP with a common seed selection for all coils, individual coil PS-SSFP with coil-wise seed selection, PS-SSFP after coil combination, and IDEAL reconstructions were also performed. Water images reconstructed via PS-SSFP methods were compared in terms of the level of fat suppression and the similarity to reference IDEAL images (signed-rank test).

RESULTS

While tissue misclassification was broadly evident across regular PS-SSFP images, the proposed method achieved significantly higher levels of fat suppression (P < 0.005) and increased similarity to reference IDEAL images (P < 0.005).

CONCLUSION

The proposed method enhances fat-water separation in phased-array acquisitions by producing improved phase estimates across the imaging volume. J. Magn. Reson. Imaging 2016;44:148-157.

Combined outer volume suppression and T2 preparation sequence for coronary angiography

PURPOSE

To develop a magnetization preparation sequence for simultaneous outer volume suppression (OVS) and T2 weighting in whole-heart coronary magnetic resonance angiography.

METHODS

A combined OVS and T2 preparation sequence (OVS-T2 Prep) was designed with a nonselective adiabatic 90° tipdown pulse, two adiabatic 180° refocusing pulses, and a 2D spiral -90° tipup pulse. The OVS-T2 Prep preserves the magnetization inside an elliptic cylinder with T2 weighting, while saturating the magnetization outside the cylinder. Its performance was tested on phantoms and on 13 normal subjects with coronary magnetic resonance angiography using 3D cones trajectories.

RESULTS

Phantom studies showed expected T2 -dependent signal amplitude in the spatial passband and suppressed signal in the spatial stopband. In vivo studies with full-field-of-view cones yielded a passband-to-stopband signal ratio of 3.18 ± 0.77 and blood-myocardium contrast-to-noise ratio enhancement by a factor of 1.43 ± 0.20 (P < 0.001). In vivo studies with reduced-field-of-view cones showed that OVS-T2 Prep well suppressed the aliasing artifacts, as supported by significantly reduced signal in the regions with no tissues compared to the images acquired without preparation (P < 0.0001).

CONCLUSION

OVS-T2 Prep is a compact sequence that can accelerate coronary magnetic resonance angiography by suppressing signals from tissues surrounding the heart while simultaneously enhancing the blood-myocardium contrast.

Water Selective Imaging and bSSFP Banding Artifact Correction in Humans and Small Animals at 3T and 7T, Respectively

INTRODUCTION

The purpose of this paper is to develop an easy method to generate both fat signal and banding artifact free 3D balanced Steady State Free Precession (bSSFP) images at high magnetic field.

METHODS

In order to suppress fat signal and bSSFP banding artifacts, two or four images were acquired with the excitation frequency of the water-selective binomial radiofrequency pulse set On Resonance or shifted by a maximum of 3/4TR. Mice and human volunteers were imaged at 7 T and 3 T, respectively to perform whole-body and musculoskeletal imaging. "Sum-Of-Square" reconstruction was performed and combined or not with parallel imaging.

RESULTS

The frequency selectivity of 1-2-3-2-1 or 1-3-3-1 binomial pulses was preserved after (3/4TR) frequency shifting. Consequently, whole body small animal 3D imaging was performed at 7 T and enabled visualization of small structures within adipose tissue like lymph nodes. In parallel, this method allowed 3D musculoskeletal imaging in humans with high spatial resolution at 3 T. The combination with parallel imaging allowed the acquisition of knee images with ~500 μm resolution images in less than 2 min. In addition, ankles, full head coverage and legs of volunteers were imaged, demonstrating the possible application of the method also for large FOV.

CONCLUSION

In conclusion, this robust method can be applied in small animals and humans at high magnetic fields. The high SNR and tissue contrast obtained in short acquisition times allows to prescribe bSSFP sequence for several preclinical and clinical applications.

Magnetic Resonance Sequences and Rapid Acquisition for MR-Guided Interventions

Interventional MR uses rapid imaging to guide diagnostic and therapeutic procedures. One of the attractions of MR-guidance is the abundance of inherent contrast mechanisms available. Dynamic procedural guidance with real-time imaging has pushed the limits of MR technology, demanding rapid acquisition and reconstruction paired with interactive control and device visualization. This article reviews the technical aspects of real-time MR sequences that enable MR-guided interventions.

Nuclear-Overhauser-enhanced MR imaging of (31)P-containing metabolites: multipoint-Dixon vs. frequency-selective excitation

The purpose of this study is to develop nuclear-Overhauser-enhanced (NOE) [(1)H]-(31)P magnetic resonance imaging (MRI) based on 3D fully-balanced steady-state free precession (fbSSFP). Therefore, two implementations of a 3D fbSSFP sequence are compared using frequency-selective excitation (FreqSel) and multipoint-Dixon (MP-Dixon). (31)P-containing model solutions and four healthy volunteers were examined at field strengths of B0=3T and 7T. Maps of the distribution of phosphocreatine (PCr), inorganic phosphate (Pi), and adenosine 5´-triphosphate (ATP) in the human calf were obtained with an isotropic resolution of 1.5cm (1.0cm) in an acquisition time of 5min (10min). NOE-pulses had the highest impact on the PCr acquisitions enhancing the signal up to (82 ± 13) % at 3T and up to (37 ± 9) % at 7T. An estimation of the level of PCr in muscle tissue from [(1)H]-(31)P MRI data yielded a mean value of (33 ± 8) mM. In conclusion, direct [(1)H]-(31)P imaging using FreqSel as well as MP-Dixon is possible in clinically feasible acquisition times. FreqSel should be preferred for measurements where only a single metabolite resonance is considered. MP-Dixon performs better in terms of SNR if a larger spectral width is of interest.

Utilization of a balanced steady state free precession signal model for improved fat/water decomposition

PURPOSE

Chemical shift based fat/water decomposition methods such as IDEAL are frequently used in challenging imaging environments with large B0 inhomogeneity. However, they do not account for the signal modulations introduced by a balanced steady state free precession (bSSFP) acquisition. Here we demonstrate improved performance when the bSSFP frequency response is properly incorporated into the multipeak spectral fat model used in the decomposition process.

THEORY AND METHODS

Balanced SSFP allows for rapid imaging but also introduces a characteristic frequency response featuring periodic nulls and pass bands. Fat spectral components in adjacent pass bands will experience bulk phase offsets and magnitude modulations that change the expected constructive and destructive interference between the fat spectral components. A bSSFP signal model was incorporated into the fat/water decomposition process and used to generate images of a fat phantom, and bilateral breast and knee images in four normal volunteers at 1.5 Tesla.

RESULTS

Incorporation of the bSSFP signal model into the decomposition process improved the performance of the fat/water decomposition.

CONCLUSION

Incorporation of this model allows rapid bSSFP imaging sequences to use robust fat/water decomposition methods such as IDEAL. While only one set of imaging parameters were presented, the method is compatible with any field strength or repetition time.

Spectrally selective imaging with wideband balanced steady-state free precession MRI

PURPOSE

Unwanted, bright fat signals in balanced steady-state free precession sequences are commonly suppressed using spectral shaping. Here, a new spectral-shaping method is proposed to significantly improve the uniformity of stopband suppression without compromising the level of passband signals.

METHODS

The proposed method combines binomial-pattern excitation pulses with a wideband balanced steady-state free precession sequence kernel. It thereby increases the frequency separation between the centers of pass and stopbands by π radians, enabling improved water-fat contrast. Simulations were performed to find the optimal flip angles and subpulse spacing for the binomial pulses that maximize contrast and signal efficiency.

RESULTS

Comparisons with a conventional binomial balanced steady-state free precession sequence were performed in simulations as well as phantom and in vivo experiments at 1.5 T and 3 T. Enhanced fat suppression is demonstrated in vivo with an average improvement of 58% in blood-fat and 68% in muscle-fat contrast (P < 0.001, Wilcoxon signed-rank test).

CONCLUSION

The proposed binomial wideband balanced steady-state free precession method is a promising candidate for spectrally selective imaging with enhanced reliability against field inhomogeneities.

Investigation of the saturation pulse artifact in non-enhanced MR angiography of the lower extremity arteries at 7 Tesla

When performing non-enhanced time-of-flight MR angiography of the lower extremity arteries at 7 T with cardiac triggering, the acquisition time is a crucial consideration. Therefore, in previous studies, saturation RF pulses were applied only every second TR. In the axial source images a slight artifact with an appearance similar to aliasing could be observed. The purpose of this study was to investigate the origin of this artifact. The reason for the artifact is supposed to be related to the two effective TRs during acquisition caused by the sparsely applied saturation RF pulse. Several sequence variants were simulated and implemented within the sequence source code to examine this hypothesis. An adaptation of the excitation flip angles for each TR as well as a correction factor for the k-space data was calculated. Additionally, a different ordering of the k-space data during acquisition was implemented as well as the combination of the latter with the k-space correction factor. The observations from the simulations were verified using both a static and a flow phantom and, finally, in a healthy volunteer using the same measurement setup as in previous volunteer and patient studies. Of all implemented techniques, only the reordering of the k-space was capable of suppressing the artifact almost completely at the cost of creating a ringing artifact. The phantom measurements showed the same results as the simulations and could thus confirm the hypothesis regarding the origin of the artifact. This was additionally verified in the healthy volunteer. The origin of the artifact could be confirmed to be the periodic signal variation caused by two effective TRs during acquisition.

Whole heart coronary imaging with flexible acquisition window and trigger delay

Coronary magnetic resonance imaging (MRI) requires a correctly timed trigger delay derived from a scout cine scan to synchronize k-space acquisition with the quiescent period of the cardiac cycle. However, heart rate changes between breath-held cine and free-breathing coronary imaging may result in inaccurate timing errors. Additionally, the determined trigger delay may not reflect the period of minimal motion for both left and right coronary arteries or different segments. In this work, we present a whole-heart coronary imaging approach that allows flexible selection of the trigger delay timings by performing k-space sampling over an enlarged acquisition window. Our approach addresses coronary motion in an interactive manner by allowing the operator to determine the temporal window with minimal cardiac motion for each artery region. An electrocardiogram-gated, k-space segmented 3D radial stack-of-stars sequence that employs a custom rotation angle is developed. An interactive reconstruction and visualization platform is then employed to determine the subset of the enlarged acquisition window for minimal coronary motion. Coronary MRI was acquired on eight healthy subjects (5 male, mean age = 37 ± 18 years), where an enlarged acquisition window of 166–220 ms was set 50 ms prior to the scout-derived trigger delay. Coronary visualization and sharpness scores were compared between the standard 120 ms window set at the trigger delay, and those reconstructed using a manually adjusted window. The proposed method using manual adjustment was able to recover delineation of five mid and distal right coronary artery regions that were otherwise not visible from the standard window, and the sharpness scores improved in all coronary regions using the proposed method. This paper demonstrates the feasibility of a whole-heart coronary imaging approach that allows interactive selection of any subset of the enlarged acquisition window for a tailored reconstruction for each branch region.

High-field MR imaging in pediatric congenital heart disease: initial results

BACKGROUND

Comprehensive assessment of pediatric congenital heart disease (CHD) at any field strength mandates evaluation of both vascular and dynamic cardiac anatomy for which diagnostic quality contrast-enhanced magnetic resonance angiography (CEMRA) and cardiac cine are crucial.

OBJECTIVE

To determine whether high-resolution (HR) CEMRA and steady-state free precession (SSFP) cine can be performed reliably at 3.0 T in children with CHD and to compare the image quality to similar techniques performed at 1.5 T.

MATERIALS AND METHODS

Twenty-eight patients with a median age of 5 months and average weight 9.0 ± 7.8 kg with suspected or known CHD were evaluated at 3.0 T. SSFP cine (n = 86 series) and HR-CEMRA (n = 414 named vascular segments) were performed and images were scored for image quality and artifacts. The findings were compared to those of 28 patients with CHD of similar weight who were evaluated at 1.5 T.

RESULTS

Overall image quality on HR-CEMRA was rated as excellent or good in 96% (397/414) of vascular segments at 3.0 T (k = 0.49) and in 94% (349/371) of vascular segments at 1.5 T (k = 0.36). Overall image quality of SSFP was rated excellent or good in 91% (78/86) of cine series at 3.0 T (k = 0.55) and in 81% (87/108) at 1.5 T (k = 0.47). Off-resonance artifact was common at both field strengths, varied over the cardiac cycle and was more prevalent at 3.0 T. At 3.0 T, off-resonance dark band artifact on SSFP cine was absent in 3% (3/86), mild in 69% (59/86), moderate in 27% (23/86) and severe in 1% (1/86) of images; at 1.5 T, dark band artifact was absent in 16% (17/108), mild in 69% (75/108), moderate in 12% (13/108) and severe in 3% (3/108) of cine images. The signal-to-noise ratio and contrast-to-noise ratio of both SSFP cine and HR-CEMRA images were significantly higher at 3.0 T than at 1.5 T (P < 0.001).

CONCLUSION

Signal-to-noise ratio and contrast-to-noise ratio of high-resolution contrast-enhanced magnetic resonance angiography and SSFP cine were higher at 3.0 T than at 1.5 T. Artifacts on SSFP cine were cardiac phase specific and more prevalent at 3.0 T such that frequency-tuning was required in one-third of exams. In neonates, high spatial resolution CEMRA was highly reliable in defining extracardiac vascular anatomy.

Ungated radial quiescent-inflow single-shot (UnQISS) magnetic resonance angiography using optimized azimuthal equidistant projections

PURPOSE

We hypothesized that non-contrast-enhanced MR angiography (NEMRA) could be performed without cardiac gating by using a variant of the quiescent-inflow single-shot (QISS) technique.

METHODS

Ungated QISS (UnQISS) MRA was evaluated in eight patients with peripheral arterial disease at 1.5T. The radial acquisition used optimized azimuthal equidistant projections, a long quiescent inflow time (1200 ms) to ensure replenishment of saturated in-plane spins irrespective of the cardiac phase, and a lengthy readout (1200 ms) so that a complete cardiac cycle was sampled for each slice. Venous and background tissue suppression was obtained using frequency-offset-corrected inversion radiofrequency pulses.

RESULTS

Scan time for UnQISS was 15.4 min for an eight-station whole-leg acquisition. The appearance of UnQISS MRA acquired using the body coil was comparable to electrocardiographic-gated QISS MRA using phased array coils. A small radial view angle increment minimized eddy current-related artifacts, whereas image quality was inferior with a golden view angle radial increment or Cartesian trajectory. In patient studies, ≥50% stenoses were consistently detected.

CONCLUSION

Using UnQISS, peripheral NEMRA can be performed without the need for cardiac gating. The use of fixed imaging parameters and body coil for signal reception further simplifies the scan procedure.

Balanced SSFP Dixon imaging with banding-artifact reduction at 3 Tesla

PURPOSE

To develop a three-dimensional (3D) balanced steady-state free-precession (bSSFP) two-point Dixon method with banding-artifact suppression to offer robust high-resolution 3D bright-fluid imaging.

METHODS

A complex sum reconstruction that combines phase-cycled bSSFP images acquired at specific echo times for robust fat/water separation without banding was investigated and compared with a magnitude-based method. Bloch simulations using both single-peak and multiple-peak fat models were performed to predict the performance of these methods for a wide range of echo times and repetition times. The quality and degree of fat/water separation was evaluated in both simulations and using in vivo imaging.

RESULTS

Simulations predicted that both effective banding-artifact suppression and substantial improvements in fat/water separation are possible at echo times that are different from conventional echo times, enabling improved spatial resolution. Comparisons between various echo times and repetition times in vivo validated the improved fat/water separation and effective banding-artifact removal predicted by the simulations.

CONCLUSION

The proposed complex sum Dixon 3D bSSFP method is able to effectively separate fat and water at different sets of echo times, while removing banding-artifacts, providing a fast, high-resolution, T2 -like sequence without blurring.

Strategies for improved 3D small-tip fast recovery imaging

PURPOSE

Small-tip fast recovery (STFR) imaging is a recently proposed steady-state sequence that has similar image contrast as balanced steady-state free precession but has the potential to simultaneously remove banding artifacts and transient fluctuation. STFR relies on a "tip-up" radiofrequency (RF) pulse tailored to the accumulated phase during the free precession (data acquisition) interval, designed to bring spins back to the longitudinal axis, thereby preserving transverse magnetization as longitudinal magnetization for the next pulse repetition time. We recently proposed an RF-spoiled STFR sequence suitable for thin slab imaging, however, in many applications, e.g., functional magnetic resonance imaging or isotropic-resolution structural imaging, three-dimensional (3D) steady-state imaging is desirable. Unfortunately, 3D STFR imaging is challenging due to the need for 3D tailored RF pulses. Here, we propose new strategies for improved 3D STFR imaging, based on (i) unspoiled imaging, and (ii) joint design of nonslice-selective tip-down/tip-up RF pulses.

THEORY AND METHODS

We derive an analytic signal model for the proposed unspoiled STFR sequence, and propose two strategies for designing the 3D tailored tip-down/tip-up RF pulses. We validate the analytic results using phantom and in vivo imaging experiments.

RESULTS

Our analytic model and imaging experiments demonstrate that the proposed unspoiled STFR sequence is less sensitive to tip-up excitation error compared to the corresponding spoiled sequence, and may, therefore, be an attractive candidate for 3D imaging. The proposed "joint" RF pulse design method, in which we formulate the tip-down/tip-up RF pulse design task as a magnitude least squares problem, produces modest improvement over a simpler "Separate" design approach. Using the proposed unspoiled sequence and joint RF pulse design, we demonstrate proof-of-principle 3D STFR brain images with balanced steady-state free precession-like signal properties but with reduced banding.

CONCLUSION

Using the proposed unspoiled sequence and joint RF pulse design, STFR brain images in a 3D region of interest with balanced steady-state free precession-like signal properties but with reduced banding can be obtained.

Self-gated fat-suppressed cardiac cine MRI

PURPOSE

To develop a self-gated alternating repetition time balanced steady-state free precession (ATR-SSFP) pulse sequence for fat-suppressed cardiac cine imaging.

METHODS

Cardiac gating is computed retrospectively using acquired magnetic resonance self-gating data, enabling cine imaging without the need for electrocardiogram (ECG) gating. Modification of the slice-select rephasing gradients of an ATR-SSFP sequence enables the acquisition of a one-dimensional self-gating readout during the unused short repetition time (TR). Self-gating readouts are acquired during every TR of segmented, breath-held cardiac scans. A template-matching algorithm is designed to compute cardiac trigger points from the self-gating signals, and these trigger points are used for retrospective cine reconstruction. The proposed approach is compared with ECG-gated ATR-SSFP and balanced steady-state free precession in 10 volunteers and five patients.

RESULTS

The difference of ECG and self-gating trigger times has a variability of 13 ± 11 ms (mean ± SD). Qualitative reviewer scoring and ranking indicate no statistically significant differences (P > 0.05) between self-gated and ECG-gated ATR-SSFP images. Quantitative blood-myocardial border sharpness is not significantly different among self-gated ATR-SSFP ( 0.61±0.15 mm -1), ECG-gated ATR-SSFP ( 0.61±0.15 mm -1), or conventional ECG-gated balanced steady-state free precession cine MRI ( 0.59±0.15 mm -1).

CONCLUSION

The proposed self-gated ATR-SSFP sequence enables fat-suppressed cardiac cine imaging at 1.5 T without the need for ECG gating and without decreasing the imaging efficiency of ATR-SSFP.

Rapid isotropic resolution cartilage assessment using radial alternating repetition time balanced steady-state free-precession imaging

PURPOSE

To compare a balanced steady-state free-precession sequence with a radial k-space trajectory and alternating repetition time fat suppression (Radial-ATR) with other currently used fat-suppressed 3D sequences for evaluating the articular cartilage of the knee joint at 3.0T.

MATERIALS AND METHODS

Radial-ATR, fast spin-echo (FSE-Cube), gradient recall-echo acquired in the steady-state (GRASS), and spoiled gradient recall-echo (SPGR) sequences with similar voxel volumes and identical scan times were performed at 3.0T on both knee joints of five volunteers. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) measurements were performed for all sequences using a double acquisition method and compared using Mann-Whitney Wilcoxon tests. Radial-ATR sequences with 0.3 mm and 0.4 mm isotropic resolution were also performed on the knee joints of seven volunteers and three patients with osteoarthritis.

RESULTS

Average SNR values for cartilage, synovial fluid, and bone marrow were 54.7, 153.3, and 12.9, respectively, for Radial ATR, 30.8, 44.1, and 1.9, respectively, for FSE-Cube, 13.3, 46.9, and 3.3, respectively, for GRASS, and 19.1, 8.1, and 2.1, respectively, for SPGR. Average CNR values between cartilage and synovial fluid and between cartilage and bone marrow were 98.6 and 41.8, respectively, for VIPR-ATR, 13.4 and 28.8, respectively, for FSE-Cube, 33.6 and 10.0, respectively, for GRASS, and 11.0 and 16.9, respectively, for SPGR. Radial-ATR had significantly higher (P < 0.001) cartilage, synovial fluid, and bone marrow SNR and significantly higher (P < 0.01) CNR between cartilage and synovial fluid and between cartilage and bone marrow than FSE-Cube, GRASS, and SPGR. Radial-ATR provided excellent visualization of articular cartilage at high isotropic resolution with no image degradation due to off-resonance banding artifacts.

CONCLUSION

Radial-ATR had superior SNR efficiency to other fat-suppressed 3D cartilage imaging sequences and produced high isotropic resolution images of the knee joint which could be used for evaluating articular cartilage at 3.0T.

In vivo visualization of polymer-based mesh implants using conventional magnetic resonance imaging and positive-contrast susceptibility imaging

PURPOSE

Polymer-based textile meshes for abdominal hernia treatment are invisible by conventional imaging methods, including magnetic resonance imaging (MRI). Integration of iron particles in the mesh base material allows MRI visualization of meshes. Positive-contrast susceptibility imaging (PCSI) was implemented to separate susceptibility-induced voids from proton-deficient voids. The purpose of this study was to compare PCSI with conventional gradient echo and turbo spin echo (TSE) sequences for the in vivo assessment of superparamagnetic iron oxide particle-loaded surgical meshes in an animal model.

METHODS AND MATERIALS

Iron-loaded polymer meshes were implanted into the abdominal wall of 10 rabbits. At days 1, 30, and 90 after surgery, conventional gradient echo, TSE, and PCSI were performed at 1.5 T in the sagittal and axial planes. Images were scored by 2 radiologists with respect to mesh visibility, delineation of the surrounding tissue, differentiation from other structures, and overall diagnostic use, on a 4-point scale ranging from 1 (insufficient) to 4 (excellent). The results were compared using Wilcoxon signed-rank tests. The mesh shape, possible deformation or fracture, and possible mesh migration were evaluated on the different pulse sequences and compared with the results at surgery and autopsy.

RESULTS

The iron-loaded meshes appeared as hypointense signal voids on gradient echo sequences, as a hyperintense line on PCSI, and as a very thin dark line on TSE images. In all animals, a precise depiction of the mesh location and its spatial configuration and integrity was possible by MRI and confirmed by surgical and autopsy results. In all 4 categories and at all 3 time points of imaging, image quality scores were significantly higher for gradient echo imaging (range, 3.60-3.80) compared with PCSI (range, 3.12-3.42) and TSE (range, 1.64-1.89). At day 90, the image quality ratings of gradient echo and PCSI were comparable. In 2 cases, the complete delineation of mesh borders was impossible because of signal voids of adjacent anatomical structures, whereas PCSI helped achieve this differentiation.

CONCLUSION

In this rabbit model of iron-loaded implanted abdominal meshes, standard gradient echo imaging was best suitable to assess implant location, integrity, and configuration. In 2 of 10 animals, PCSI helped achieve a complete delineation of mesh borders.

High-resolution 3D radial bSSFP with IDEAL

Radial trajectories facilitate high-resolution balanced steady state free precession (bSSFP) because the efficient gradients provide more time to extend the trajectory in k-space. A number of radial bSSFP methods that support fat-water separation have been developed; however, most of these methods require an environment with limited B0 inhomogeneity. In this work, high-resolution bSSFP with fat-water separation is achieved in more challenging B0 environments by combining a 3D radial trajectory with the IDEAL chemical species separation method. A method to maintain very high resolution within the timing constraints of bSSFP and IDEAL is described using a dual-pass pulse sequence. The sampling of a unique set of radial lines at each echo time is investigated as a means to circumvent the longer scan time that IDEAL incurs as a multiecho acquisition. The manifestation of undersampling artifacts in this trajectory and their effect on chemical species separation are investigated in comparison to the case in which each echo samples the same set of radial lines. This new bSSFP method achieves 0.63 mm isotropic resolution in a 5-min scan and is demonstrated in difficult in vivo imaging environments, including the breast and a knee with ACL reconstruction hardware at 1.5 T.

Small-tip fast recovery imaging using non-slice-selective tailored tip-up pulses and radiofrequency-spoiling

Small-tip fast recovery (STFR) imaging is a new steady-state imaging sequence that is a potential alternative to balanced steady-state free precession. Under ideal imaging conditions, STFR may provide comparable signal-to-noise ratio and image contrast as balanced steady-state free precession, but without signal variations due to resonance offset. STFR relies on a tailored "tip-up," or "fast recovery," radiofrequency pulse to align the spins with the longitudinal axis after each data readout segment. The design of the tip-up pulse is based on the acquisition of a separate off-resonance (B0) map. Unfortunately, the design of fast (a few ms) slice- or slab-selective radiofrequency pulses that accurately tailor the excitation pattern to the local B0 inhomogeneity over the entire imaging volume remains a challenging and unsolved problem. We introduce a novel implementation of STFR imaging based on "non-slice-selective" tip-up pulses, which simplifies the radiofrequency pulse design problem significantly. Out-of-slice magnetization pathways are suppressed using radiofrequency-spoiling. Brain images obtained with this technique show excellent gray/white matter contrast, and point to the possibility of rapid steady-state T(2)/T(1) -weighted imaging with intrinsic suppression of cerebrospinal fluid, through-plane vessel signal, and off-resonance artifacts. In the future, we expect STFR imaging to benefit significantly from parallel excitation hardware and high-order gradient shim systems.

Stages: sub-Fourier dynamic shim updating using nonlinear magnetic field phase preparation

Heterogeneity of the static magnetic field in magnetic resonance imaging may cause image artifacts and degradation in image quality. The field heterogeneity can be reduced by dynamically adjusting shim fields or dynamic shim updating, in which magnetic field homogeneity is optimized for each tomographic slice to improve image quality. A limitation of this approach is that a new magnetic field can be applied only once for each slice, otherwise image quality would improve somewhere to its detriment elsewhere in the slice. The motivation of this work is to overcome this limitation and develop a technique using nonlinear magnetic fields to dynamically shim the static magnetic field within a single Fourier-encoded volume or slice, called sub-Fourier dynamic shim updating. However, the nonlinear magnetic fields are not used as shim fields; instead, they impart a strong spatial dependence to the acquired MR signal by nonlinear phase preparation, which may be exploited to locally improve magnetic field homogeneity during acquisition. A theoretical description of the method is detailed, simulations and a proof-of-principle experiment are performed using a magnet coil with a known field geometry. The method is shown to remove artifacts associated with magnetic field homogeneity in balanced steady-state free-precession pulse sequences. We anticipate that this method will be useful to improve the quality of magnetic resonance images by removing deleterious artifacts associated with a heterogeneous static magnetic field.

Noncontrast-enhanced renal angiography using multiple inversion recovery and alternating TR balanced steady-state free precession

Noncontrast-enhanced renal angiography techniques based on balanced steady-state free precession avoid external contrast agents, take advantage of high inherent blood signal from the T 2 / T 1 contrast mechanism, and have short steady-state free precession acquisition times. However, background suppression is limited; inflow times are inflexible; labeling region is difficult to define when tagging arterial flow; and scan times are long. To overcome these limitations, we propose the use of multiple inversion recovery preparatory pulses combined with alternating pulse repetition time balanced steady-state free precession to produce renal angiograms. Multiple inversion recovery uses selective spatial saturation followed by four nonselective inversion recovery pulses to concurrently null a wide range of background T 1 species while allowing for adjustable inflow times; alternating pulse repetition time steady-state free precession maintains vessel contrast and provides added fat suppression. The high level of suppression enables imaging in three-dimensional as well as projective two-dimensional formats, the latter of which has a scan time as short as one heartbeat. In vivo studies at 1.5 T demonstrate the superior vessel contrast of this technique.

Rapid gradient-echo imaging

Gradient-echo sequences are widely used in magnetic resonance imaging (MRI) for numerous applications ranging from angiography to perfusion to functional MRI. Compared with spin-echo techniques, the very short repetition times of gradient-echo methods enable very rapid 2D and 3D imaging, but also lead to complicated "steady states." Signal and contrast behavior can be described graphically and mathematically, and depends strongly on the type of spoiling: fully balanced (no spoiling), gradient spoiling, or radiofrequency (RF)-spoiling. These spoiling options trade off between high signal and pure T(1) contrast, while the flip angle also affects image contrast in all cases, both of which can be demonstrated theoretically and in image examples. As with spin-echo sequences, magnetization preparation can be added to gradient-echo sequences to alter image contrast. Gradient-echo sequences are widely used for numerous applications such as 3D perfusion imaging, functional MRI, cardiac imaging, and MR angiography.

Frequency-specific SSFP for hyperpolarized ¹³C metabolic imaging at 14.1 T

Metabolic imaging of hyperpolarized [1-(13)C] pyruvate co-polarized with [(13)C]urea by dynamic nuclear polarization with rapid dissolution is a promising new method for assessing tumor metabolism and perfusion simultaneously in vivo. Novel pulse sequences are required to enable dynamic imaging of multiple (13)C spectral lines with high spatiotemporal resolution. The goal of this study was to investigate a new frequency-specific approach for rapid metabolic imaging of multiple (13)C resonances using the spectral selectivity of steady-state free precession pulse (SSFP) trains. Methods developed in simulations were implemented in a dynamic frequency-cycled balanced SSFP pulse sequence on a 14.1-T animal magnetic resonance imaging scanner. This acquisition was tested in thermal and hyperpolarized phantom imaging studies and in a transgenic mouse with prostate cancer.

ISMRM workshop on fat-water separation: insights, applications and progress in MRI

Approximately 130 attendees convened on February 19-22, 2012 for the first ISMRM-sponsored workshop on water-fat imaging. The motivation to host this meeting was driven by the increasing number of research publications on this topic over the past decade. The scientific program included an historical perspective and a discussion of the clinical relevance of water-fat MRI, a technical description of multiecho pulse sequences, a review of data acquisition and reconstruction algorithms, a summary of the confounding factors that influence quantitative fat measurements and the importance of MRI-based biomarkers, a description of applications in the heart, liver, pancreas, abdomen, spine, pelvis, and muscles, an overview of the implications of fat in diabetes and obesity, a discussion on MR spectroscopy, a review of childhood obesity, the efficacy of lifestyle interventional studies, and the role of brown adipose tissue, and an outlook on federal funding opportunities from the National Institutes of Health.

Free-breathing multiphase whole-heart coronary MR angiography using image-based navigators and three-dimensional cones imaging

Noninvasive visualization of the coronary arteries in vivo is one of the most important goals in cardiovascular imaging. Compared to other paradigms for coronary MR angiography, a free-breathing three-dimensional whole-heart iso-resolution approach simplifies prescription effort, requires less patient cooperation, reduces overall exam time, and supports retrospective reformats at arbitrary planes. However, this approach requires a long continuous acquisition and must account for respiratory and cardiac motion throughout the scan. In this work, a new free-breathing coronary MR angiography technique that reduces scan time and improves robustness to motion is developed. Data acquisition is accomplished using a three-dimensional cones non-Cartesian trajectory, which can reduce the number of readouts 3-fold or more compared to conventional three-dimensional Cartesian encoding and provides greater robustness to motion/flow effects. To further enhance robustness to motion, two-dimensional navigator images are acquired to directly track respiration-induced displacement of the heart and enable retrospective compensation of all acquired data (none discarded) for image reconstruction. In addition, multiple cardiac phases are imaged to support retrospective selection of the best phase(s) for visualizing each coronary segment. Experimental results demonstrate that whole-heart coronary angiograms can be obtained rapidly and robustly with this proposed technique.

Use of simulated annealing for the design of multiple repetition time balanced steady-state free precession imaging

Balanced steady-state free precession is an ultrafast sequence with high signal-to-noise efficiency, but it also generates a strong fat signal which can mask important features. One method of fat suppression is to modify the balanced steady-state free precession spectrum using multiple repetition times to create a wide stopband over the fat frequency. However, with three or more pulse repetition times, the number of parameters creates a vast search space with many local minima of a cost function. We report on the initial results of using simulated annealing to find optimal sequences for two applications of multiple-pulse repetition time balanced steady-state free precession: positive contrast imaging and fat suppression.

Nonuniform and multidimensional Shinnar-Le Roux RF pulse design method

The Shinnar-Le Roux (SLR) radiofrequency (RF) pulse design algorithm is widely used for designing slice-selective RF pulses due to its intuitiveness, optimality, and speed. SLR is limited, however, in that it is only capable of designing one-dimensional pulses played along constant gradients. We present a nonuniform SLR RF pulse design framework that extends most of the capabilities of classical SLR to nonuniform gradient trajectories and multiple dimensions. Specifically, like classical SLR, the new method is a hard pulse approximation-based technique that uses filter design relationships to produce the lowest power RF pulse that satisfies target magnetization ripple levels. The new method is validated and compared with methods conventionally used for nonuniform and multidimensional large-tip-angle RF pulse design.

The central signal singularity phenomenon in balanced SSFP and its application to positive-contrast imaging

Small perturbations of steady-state sequence parameters can induce very large spectral profile deviations that are localized to specific off-resonant frequencies, denoted critical frequencies. Although, a small number of studies have previously considered the use of these highly specific modulations for MR angiography and elastography, many potential applications still remain to be explored. An analysis of this phenomenon using a linear systems technique and a geometric magnetization trajectory technique shows that the critical frequencies correspond to singularities in the steady-state signal equation. An interleaved acquisition combined with a complex difference technique yields a spectral profile containing sharp peaks interleaved with wide stopbands, while a complex sum technique yields a spectral profile similar to that of balanced steady-state free precession. Simulations and phantom experiments are used to demonstrate a novel application of this technique for positive-contrast imaging of superparamagnetic iron-oxide nanoparticles. The technique is shown to yield images with high levels of positive contrast and good water and fat background suppression. The technique can also simultaneously yield images with contrast similar to balanced steady-state free precession.

Linear Response Equilibrium versus echo-planar encoding for fast high-spatial resolution 3D chemical shift imaging

In this work Linear Response Equilibrium (LRE) and Echo-planar spectroscopic imaging (EPSI) are compared in terms of sensitivity per unit time and power deposition. In addition an extended dual repetition time scheme to generate broad stopbands for improved inherent water suppression in LRE is presented. The feasibility of LRE and EPSI for assessing cholesterol esters in human carotid plaques with high spatial resolution of 1.95×1.15×1.15 mm(3) on a clinical 3T MR system is demonstrated. In simulations and phantom experiments it is shown that LRE has comparable but lower sensitivity per unit time relative to EPSI despite stronger signal generated. This relates to the lower sampling efficiency in LRE relative to EPSI as a result of limited gradient performance on clinical MR systems. At the same time, power deposition of LRE is significantly reduced compared to EPSI making it an interesting niche application for in vivo high field spectroscopic imaging of metabolites within a limited bandwidth.

Signal compensation and compressed sensing for magnetization-prepared MR angiography

Magnetization-prepared acquisitions offer a trade-off between image contrast and scan efficiency for magnetic resonance imaging. Because the prepared signals gradually decay, the contrast can be improved by frequently repeating the preparation, which in turn significantly increases the scan time. A common solution is to perform the data collection progressing from low- to high-spatial-frequency samples following each preparation. Unfortunately, this leads to loss of spatial resolution, and thereby image blurring. In this work, a new technique is proposed that first corrects the signal decay in high-frequency data to mitigate the resolution loss and improve the image contrast without reducing the scan efficiency. The proposed technique then employs a sparsity-based nonlinear reconstruction to further improve the image quality. In addition to reducing the amplified high-frequency noise, this reconstruction extrapolates missing k-space samples in the case of undersampled compressed-sensing acquisitions. The technique is successfully demonstrated for noncontrast-enhanced flow-independent angiography of the lower extremities, an application that substantially benefits from both the signal compensation and the nonlinear reconstruction.