Myocardial perfusion MRI with an undersampled 3D stack-of-stars sequence

Arterial Input Function (AIF) Correction Using AIF Plus Tissue Inputs with a Bi-LSTM Network

Background: The arterial input function (AIF) is vital for myocardial blood flow quantification in cardiac MRI to indicate the input time-concentration curve of a contrast agent. Inaccurate AIFs can significantly affect perfusion quantification. Purpose: When only saturated and biased AIFs are measured, this work investigates multiple ways of leveraging tissue curve information, including using AIF + tissue curves as inputs and optimizing the loss function for deep neural network training. Methods: Simulated data were generated using a 12-parameter AIF mathematical model for the AIF. Tissue curves were created from true AIFs combined with compartment-model parameters from a random distribution. Using Bloch simulations, a dictionary was constructed for a saturation-recovery 3D radial stack-of-stars sequence, accounting for deviations such as flip angle, T2* effects, and residual longitudinal magnetization after the saturation. A preliminary simulation study established the optimal tissue curve number using a bidirectional long short-term memory (Bi-LSTM) network with just AIF loss. Further optimization of the loss function involves comparing just AIF loss, AIF with compartment-model-based parameter loss, and AIF with compartment-model tissue loss. The optimized network was examined with both simulation and hybrid data, which included in vivo 3D stack-of-star datasets for testing. The AIF peak value accuracy and ktrans results were assessed. Results: Increasing the number of tissue curves can be beneficial when added tissue curves can provide extra information. Using just the AIF loss outperforms the other two proposed losses, including adding either a compartment-model-based tissue loss or a compartment-model parameter loss to the AIF loss. With the simulated data, the Bi-LSTM network reduced the AIF peak error from -23.6 ± 24.4% of the AIF using the dictionary method to 0.2 ± 7.2% (AIF input only) and 0.3 ± 2.5% (AIF + ten tissue curve inputs) of the network AIF. The corresponding ktrans error was reduced from -13.5 ± 8.8% to -0.6 ± 6.6% and 0.3 ± 2.1%. With the hybrid data (simulated data for training; in vivo data for testing), the AIF peak error was 15.0 ± 5.3% and the corresponding ktrans error was 20.7 ± 11.6% for the AIF using the dictionary method. The hybrid data revealed that using the AIF + tissue inputs reduced errors, with peak error (1.3 ± 11.1%) and ktrans error (-2.4 ± 6.7%). Conclusions: Integrating tissue curves with AIF curves into network inputs improves the precision of AI-driven AIF corrections. This result was seen both with simulated data and with applying the network trained only on simulated data to a limited in vivo test dataset.

Quantitative myocardial perfusion with a hybrid 2D simultaneous multi-slice sequence

PURPOSE

To evaluate a novel 2D simultaneous multi-slice (SMS) myocardial perfusion acquisition and compare directly to a published quantitative 3D stack-of-stars (SoS) acquisition.

METHODS

A hybrid saturation recovery radial 2D SMS sequence following a single saturation was created for the quantification of myocardial blood flow (MBF). This sequence acquired three slices simultaneously and generated an arterial input function (AIF) using the first 24 rays. Validation was done in a novel way by alternating heartbeats between the hybrid 2D SMS and the 3D SoS acquisitions. Initial studies were done to study the effects of using only every other beat for the 2D SMS in two subjects, and for the 3D SoS in four subjects. The proposed alternating acquisitions were then performed in ten dog studies at rest, four dog studies at adenosine stress, and two human resting studies. Quantitative MBF analysis was performed for 2D SMS and 3D SoS separately, using a compartment model.

RESULTS

Acquiring every-other-beat data resulted in 6 ± 5% ("ideal") and 11 ± 8% ("practical") perfusion changes for both 2D SMS and 3D SoS methods. For alternating acquisitions, 2D SMS and 3D SoS quantitative perfusion values were comparable for both the twelve rest studies (2D SMS: 0.69 ± 0.16 vs 3D: 0.69 ± 0.15 ml/g/min, p = 0.55) and the four stress studies (2D SMS: 1.28 ± 0.22 vs 3D: 1.30 ± 0.24 ml/g/min, p = 0.61).

CONCLUSION

Every-other-beat acquisition changed estimated perfusion values relatively little for both sequences. The quantitative hybrid radial 2D SMS myocardial first-pass perfusion imaging sequence gave results similar to 3D perfusion when compared directly with an alternating beat acquisition.

Free-breathing motion-informed locally low-rank quantitative 3D myocardial perfusion imaging

Purpose

To propose respiratory motion‐informed locally low‐rank reconstruction (MI‐LLR) for robust free‐breathing single‐bolus quantitative 3D myocardial perfusion CMR imaging. Simulation and in‐vivo results are compared to locally low‐rank (LLR) and compressed sensing reconstructions (CS) for reference.

Methods

Data were acquired using a 3D Cartesian pseudo‐spiral in‐out k‐t undersampling scheme (R = 10) and reconstructed using MI‐LLR, which encompasses two stages. In the first stage, approximate displacement fields are derived from an initial LLR reconstruction to feed a motion‐compensated reference system to a second reconstruction stage, which reduces the rank of the inverse problem. For comparison, data were also reconstructed with LLR and frame‐by‐frame CS using wavelets as sparsifying transform (ℓ1‐wavelet).

Reconstruction accuracy relative to ground truth was assessed using synthetic data for realistic ranges of breathing motion, heart rates, and SNRs. In‐vivo experiments were conducted in healthy subjects at rest and during adenosine stress. Myocardial blood flow (MBF) maps were derived using a Fermi model.

Results

Improved uniformity of MBF maps with reduced local variations was achieved with MI‐LLR. For rest and stress, intra‐volunteer variation of absolute and relative MBF was lower in MI‐LLR (±0.17 mL/g/min [26%] and ±1.07 mL/g/min [33%]) versus LLR (±0.19 mL/g/min [28%] and ±1.22 mL/g/min [36%]) and versus ℓ1‐wavelet (±1.17 mL/g/min [113%] and ±6.87 mL/g/min [115%]). At rest, intra‐subject MBF variation was reduced significantly with MI‐LLR.

Conclusion

The combination of pseudo‐spiral Cartesian undersampling and dual‐stage MI‐LLR reconstruction improves free‐breathing quantitative 3D myocardial perfusion CMR imaging under rest and stress condition.

Click here for author‐reader discussions

High-Resolution Free-Breathing Quantitative First-Pass Perfusion Cardiac MR Using Dual-Echo Dixon With Spatio-Temporal Acceleration

Introduction

To develop and test the feasibility of free-breathing (FB), high-resolution quantitative first-pass perfusion cardiac MR (FPP-CMR) using dual-echo Dixon (FOSTERS; Fat-water separation for mOtion-corrected Spatio-TEmporally accelerated myocardial peRfuSion).

Materials and Methods

FOSTERS was performed in FB using a dual-saturation single-bolus acquisition with dual-echo Dixon and a dynamically variable Cartesian k-t undersampling (8-fold) approach, with low-rank and sparsity constrained reconstruction, to achieve high-resolution FPP-CMR images. FOSTERS also included automatic in-plane motion estimation and T2* correction to obtain quantitative myocardial blood flow (MBF) maps. High-resolution (1.6 x 1.6 mm²) FB FOSTERS was evaluated in eleven patients, during rest, against standard-resolution (2.6 x 2.6 mm²) 2-fold SENSE-accelerated breath-hold (BH) FPP-CMR. In addition, MBF was computed for FOSTERS and spatial wavelet-based compressed sensing (CS) reconstruction. Two cardiologists scored the image quality (IQ) of FOSTERS, CS, and standard BH FPP-CMR images using a 4-point scale (1–4, non-diagnostic – fully diagnostic).

Results

FOSTERS produced high-quality images without dark-rim and with reduced motion-related artifacts, using an 8x accelerated FB acquisition. FOSTERS and standard BH FPP-CMR exhibited excellent IQ with an average score of 3.5 ± 0.6 and 3.4 ± 0.6 (no statistical difference, p > 0.05), respectively. CS images exhibited severe artifacts and high levels of noise, resulting in an average IQ score of 2.9 ± 0.5. MBF values obtained with FOSTERS presented a lower variance than those obtained with CS.

Discussion

FOSTERS enabled high-resolution FB FPP-CMR with MBF quantification. Combining motion correction with a low-rank and sparsity-constrained reconstruction results in excellent image quality.

Direct imaging of white matter ultrashort T2∗ components at 7 Tesla

PURPOSE

To demonstrate direct imaging of the white matter ultrashort T₂^∗ components at 7 Tesla using inversion recovery (IR)-enhanced ultrashort echo time (UTE) MRI. To investigate its characteristics, potentials and limitations, and to establish a clinical protocol.

MATERIAL AND METHODS

The IR UTE technique suppresses long T₂^∗ signals within white matter by using adiabatic inversion in combination with dual-echo difference imaging. Artifacts arising at 7 T from long T₂^∗ scalp fat components were reduced by frequency shifting the IR pulse such that those frequencies were inverted likewise. For 8 healthy volunteers, the T₂^∗ relaxation times of white matter were then quantified. In 20 healthy volunteers, the UTE difference and fraction contrast were evaluated. Finally, in 6 patients with multiple sclerosis (MS), the performance of the technique was assessed.

RESULTS

A frequency shift of -1.2 ppm of the IR pulse (i.e. towards the fat frequency) provided a good suppression of artifacts. With this, an ultrashort compartment of (68 ± 6) % with a T₂^∗ time of (147 ± 58) μs was quantified with a chemical shift of (-3.6 ± 0.5) ppm from water. Within healthy volunteers' white matter, a stable ultrashort T₂^∗ fraction contrast was calculated. For the MS patients, a significant fraction reduction in the identified lesions as well as in the normal-appearing white matter was observed.

CONCLUSIONS

The quantification results indicate that the observed ultrashort components arise primarily from myelin tissue. Direct IR UTE imaging of the white matter ultrashort T₂^∗ components is thus feasible at 7 T with high quantitative inter-subject repeatability and good detection of signal loss in MS.

Initial investigation of free-breathing 3D whole-heart stress myocardial perfusion MRI

Objective: Myocardial first-pass perfusion imaging with MRI is well-established clinically. However, it is potentially weakened by limited myocardial coverage compared to nuclear medicine. Clinical evaluations of whole-heart MRI perfusion by 3D methods, while promising, have to date had the limit of breathhold requirements at stress. This work aims to develop a new free-breathing 3D myocardial perfusion method, and to test its performance in a small patient population.

Methods: This work required tolerance to respiratory motion for stress investigations, and therefore employed a “stack-of-stars” hybrid Cartesian-radial MRI acquisition method. The MRI sequence was highly optimised for rapid acquisition and combined with a compressed sensing reconstruction. Stress and rest datasets were acquired in four healthy volunteers, and in six patients with coronary artery disease (CAD), which were compared against clinical reference information.

Results: This free-breathing method produced datasets that appeared consistent with clinical reference data in detecting moderate-to-strong induced perfusion abnormalities. However, the majority of the mild defects identified clinically were not detected by the method, potentially due to the presence of transient myocardial artefacts present in the images.

Discussion: The feasibility of detecting CAD using this 3D first-pass perfusion sequence during free-breathing is demonstrated. Good agreement on typical moderate-to-strong CAD cases is promising, however, questions still remain on the sensitivity of the technique to milder cases.

Whole-heart, ungated, free-breathing, cardiac-phase-resolved myocardial perfusion MRI by using Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state (CRIMP)

PURPOSE

To develop a whole-heart, free-breathing, non-electrocardiograph (ECG)-gated, cardiac-phase-resolved myocardial perfusion MRI framework (CRIMP; Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state) and test its quantification feasibility.

METHODS

CRIMP used interleaved radial simultaneous multi-slice (SMS) slice groups to cover the whole heart in 9 or 12 short-axis slices. The sequence continuously acquired data without magnetization preparation, ECG gating or breath-holding, and captured multiple cardiac phases. Images were reconstructed by a motion-compensated patch-based locally low-rank reconstruction. Bloch simulations were performed to study the signal-to-noise ratio/contrast-to-noise ratio (SNR/CNR) for CRIMP and to study the steady-state signal under motion. Seven patients were scanned with CRIMP at stress and rest to develop the sequence. One human and two dogs were scanned at rest with a dual-bolus method to test the quantification feasibility of CRIMP. The dual-bolus scans were performed using both CRIMP and an ungated radial SMS saturation recovery (SMS-SR) sequence with injection dose = 0.075 mmol/kg to compare the sequences in terms of SNR, cardiac phase resolution and quantitative myocardial blood flow (MBF).

RESULTS

Perfusion images with multiple cardiac phases in all image slices with a temporal resolution of 72 ms/frame were obtained. Simulations and in-vivo acquisitions showed CRIMP kept the inner slices in steady-state regardless of motion. CRIMP outperformed SMS-SR in slice coverage (9 over 6), SNR (mean 20% improvement), and provided cardiac phase resolution. CRIMP and SMS-SR sequences provided comparable MBF values (rest systolic CRIMP = 0.58 ± 0.07, SMS-SR = 0.61 ± 0.16).

CONCLUSION

CRIMP allows for whole-heart, cardiac-phase-resolved myocardial perfusion images without ECG-gating or breath-holding. The sequence can provide MBF if an accurate arterial input function is obtained separately.

Quantitative 3D myocardial perfusion with an efficient arterial input function

PURPOSE

The purpose of this study was to further develop and combine several innovative sequence designs to achieve quantitative 3D myocardial perfusion. These developments include an optimized 3D stack-of-stars readout (150 ms per beat), efficient acquisition of a 2D arterial input function, tailored saturation pulse design, and potential whole heart coverage during quantitative stress perfusion.

THEORY AND METHODS

All studies were performed free-breathing on a Prisma 3T MRI scanner. Phantom validation was used to verify sequence accuracy. A total of 21 subjects (3 patients with known disease) were scanned, 12 with a rest only protocol and 9 with both stress (regadenoson) and rest protocols. First pass quantitative perfusion was performed with gadoteridol (0.075 mmol/kg).

RESULTS

Implementation and quantitative perfusion results are shown for healthy subjects and subjects with known coronary disease. Average rest perfusion for the 15 included healthy subjects was 0.79 ± 0.19 mL/g/min, the average stress perfusion for 6 healthy subject studies was 2.44 ± 0.61 mL/g/min, and the average global myocardial perfusion reserve ratio for 6 healthy subjects was 3.10 ± 0.24. Perfusion deficits for 3 patients with ischemia are shown. Average resting heart rate was 59 ± 7 bpm and the average stress heart rate was 81 ± 10 bpm.

CONCLUSION

This work demonstrates that a quantitative 3D myocardial perfusion sequence with the acquisition of a 2D arterial input function is feasible at high stress heart rates such as during stress. T1 values and gadolinium concentrations of the sequence match the reference standard well in a phantom, and myocardial rest and stress perfusion and myocardial perfusion reserve values are consistent with those published in literature.

Feasibility of multiple-view myocardial perfusion MRI using radial simultaneous multi-slice acquisitions

Purpose

Dynamic contrast enhanced MRI of the heart typically acquires 2–4 short-axis (SA) slices to detect and characterize coronary artery disease. This acquisition scheme is limited by incomplete coverage of the left ventricle. We studied the feasibility of using radial simultaneous multi-slice (SMS) technique to achieve SA, 2-chamber and/or 4-chamber long-axis (2CH LA and/or 4CH LA) coverage with and without electrocardiography (ECG) gating using a motion-robust reconstruction framework.

Methods

12 subjects were scanned at rest and/or stress, free breathing, with or without ECG gating. Multiple sets of radial SMS k-space were acquired within each cardiac cycle, and each SMS set sampled 3 parallel slices that were either SA, 2CH LA, or 4CH LA slices. The radial data was interpolated onto Cartesian space using an SMS GRAPPA operator gridding method. Self-gating and respiratory states binning of the data were done. The binning information as well as a pixel tracking spatiotemporal constrained reconstruction method were applied to obtain motion-robust image reconstructions. Reconstructions with and without the pixel tracking method were compared for signal-to-noise ratio and contrast-to-noise ratio.

Results

Full coverage of the heart (at least 3 SA and 3 LA slices) during the first pass of contrast at every heartbeat was achieved by using the radial SMS acquisition. The proposed pixel tracking reconstruction improves the average SNR and CNR by 21% and 30% respectively, and reduces temporal blurring for both gated and ungated acquisitions.

Conclusion

Acquiring simultaneous multi-slice SA, 2CH LA and/or 4CH LA myocardial perfusion images in every heartbeat is feasible in both gated and ungated acquisitions. This can add confidence when detecting and characterizing coronary artery disease by revealing ischemia in different views, and by providing apical coverage that is improved relative to SA slices alone. The proposed pixel tracking framework improves the reconstruction while adding little computational cost.

Sparsity adaptive reconstruction for highly accelerated cardiac MRI

PURPOSE

To enable parameter-free, accelerated cardiovascular magnetic resonance (CMR).

METHODS

Regularized reconstruction methods, such as compressed sensing (CS), can significantly accelerate MRI data acquisition but require tuning of regularization weights. In this work, a technique, called Sparsity adaptive Composite Recovery (SCoRe) that exploits sparsity in multiple, disparate sparsifying transforms is presented. A data-driven adjustment of the relative contributions of different transforms yields a parameter-free CS recovery process. SCoRe is validated in a dynamic digital phantom as well as in retrospectively and prospectively undersampled cine CMR data.

RESULTS

The results from simulation and 6 retrospectively undersampled datasets indicate that SCoRe with auto-tuned regularization weights yields lower root-mean-square error (RMSE) and higher structural similarity index (SSIM) compared to state-of-the-art CS methods. In 45 prospectively undersampled datasets acquired from 15 volunteers, the image quality was scored by 2 expert reviewers, with SCoRe receiving a higher average score (p  <  0.01) compared to other CS methods.

CONCLUSIONS

SCoRe enables accelerated cine CMR from highly undersampled data. In contrast to other acceleration techniques, SCoRe adapts regularization weights based on noise power and level of sparsity in each transform, yielding superior performance without admitting any free parameters.

Golden-ratio rotated stack-of-stars acquisition for improved volumetric MRI

PURPOSE

To develop and evaluate an improved stack-of-stars radial sampling strategy for reducing streaking artifacts.

METHODS

The conventional stack-of-stars sampling strategy collects the same radial angle for every partition (slice) encoding. In an undersampled acquisition, such an aligned acquisition generates coherent aliasing patterns and introduces strong streaking artifacts. We show that by rotating the radial spokes in a golden-angle manner along the partition-encoding direction, the aliasing pattern is modified, resulting in improved image quality for gridding and more advanced reconstruction methods. Computer simulations were performed and phantom as well as in vivo images for three different applications were acquired.

RESULTS

Simulation, phantom, and in vivo experiments confirmed that the proposed method was able to generate images with less streaking artifact and sharper structures based on undersampled acquisitions in comparison with the conventional aligned approach at the same acceleration factors. By combining parallel imaging and compressed sensing in the reconstruction, streaking artifacts were mostly removed with improved delineation of fine structures using the proposed strategy.

CONCLUSIONS

We present a simple method to reduce streaking artifacts and improve image quality in 3D stack-of-stars acquisitions by re-arranging the radial spoke angles in the 3D partition direction, which can be used for rapid volumetric imaging. Magn Reson Med 78:2290-2298, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Analysis of spatiotemporal fidelity in quantitative 3D first-pass perfusion cardiovascular magnetic resonance

BACKGROUND

Whole-heart first-pass perfusion cardiovascular magnetic resonance (CMR) relies on highly accelerated image acquisition. The influence of undersampling on myocardial blood flow (MBF) quantification has not been systematically investigated yet. In the present work, the effect of spatiotemporal scan acceleration on image reconstruction accuracy and MBF error was studied using a numerical phantom and validated in-vivo.

METHODS

Up to 10-fold scan acceleration using k-t PCA and k-t SPARSE-SENSE was simulated using the MRXCAT CMR numerical phantom framework. Image reconstruction results were compared to ground truth data in the k-f domain by means of modulation transfer function (MTF) analysis. In the x-t domain, errors pertaining to specific features of signal intensity-time curves and MBF values derived using Fermi model deconvolution were analysed. In-vivo first-pass CMR data were acquired in ten healthy volunteers using a dual-sequence approach assessing the arterial input function (AIF) and myocardial enhancement. 10x accelerated 3D k-t PCA and k-t SPARSE-SENSE were compared and related to non-accelerated 2D reference images.

RESULTS

MTF analysis revealed good recovery of data upon k-t PCA reconstruction at 10x undersampling with some attenuation of higher temporal frequencies. For 10x k-t SPARSE-SENSE the MTF was found to decrease to zero at high spatial frequencies for all temporal frequencies indicating a loss in spatial resolution. Signal intensity-time curve errors were most prominent in AIFs from 10x k-t PCA, thereby emphasizing the need for separate AIF acquisition using a dual-sequence approach. These findings were confirmed by MBF estimation based on AIFs from fully sampled and undersampled simulations. Average in-vivo MBF estimates were in good agreement between both accelerated and the fully sampled methods. Intra-volunteer MBF variation for fully sampled 2D scans was lower compared to 10x k-t PCA and k-t SPARSE-SENSE data.

CONCLUSION

Quantification of highly undersampled 3D first-pass perfusion CMR yields accurate MBF estimates provided the AIF is obtained using fully sampled or moderately undersampled scans as part of a dual-sequence approach. However, relative to fully sampled 2D perfusion imaging, intra-volunteer variation is increased using 3D approaches prompting for further developments.

First-pass myocardial perfusion MRI with reduced subendocardial dark-rim artifact using optimized Cartesian sampling

PURPOSE

The presence of subendocardial dark-rim artifact (DRA) remains an ongoing challenge in first-pass perfusion (FPP) cardiac magnetic resonance imaging (MRI). We propose a free-breathing FPP imaging scheme with Cartesian sampling that is optimized to minimize the DRA and readily enables near-instantaneous image reconstruction.

MATERIALS AND METHODS

The proposed FPP method suppresses Gibbs ringing effects-a major underlying factor for the DRA-by "shaping" the underlying point spread function through a two-step process: 1) an undersampled Cartesian sampling scheme that widens the k-space coverage compared to the conventional scheme; and 2) a modified parallel-imaging scheme that incorporates optimized apodization (k-space data filtering) to suppress Gibbs-ringing effects. Healthy volunteer studies (n = 10) were performed to compare the proposed method against the conventional Cartesian technique-both using a saturation-recovery gradient-echo sequence at 3T. Furthermore, FPP imaging studies using the proposed method were performed in infarcted canines (n = 3), and in two symptomatic patients with suspected coronary microvascular dysfunction for assessment of myocardial hypoperfusion.

RESULTS

Width of the DRA and the number of DRA-affected myocardial segments were significantly reduced in the proposed method compared to the conventional approach (width: 1.3 vs. 2.9 mm, P < 0.001; number of segments: 2.6 vs. 8.7; P < 0.0001). The number of slices with severe DRA was markedly lower for the proposed method (by 10-fold). The reader-assigned image quality scores were similar (P = 0.2), although the quantified myocardial signal-to-noise ratio was lower for the proposed method (P < 0.05). Animal studies showed that the proposed method can detect subendocardial perfusion defects and patient results were consistent with the gold-standard invasive test.

CONCLUSION

The proposed free-breathing Cartesian FPP imaging method significantly reduces the prevalence of severe DRAs compared to the conventional approach while maintaining similar resolution and image quality.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:542-555.

Radial simultaneous multi-slice CAIPI for ungated myocardial perfusion

OBJECTIVE

Simultaneous multi-slice (SMS) imaging is a slice acceleration technique that acquires multiple slices in the same time as a single slice. Radial controlled aliasing in parallel imaging results in higher acceleration (radial CAIPIRINHA or CAIPI) is a promising SMS method with less severe slice aliasing artifacts as compared to its Cartesian counterpart. Here we use radial CAIPI with data undersampling and constrained reconstruction to improve the utility of ungated cardiac perfusion acquisitions. We test the proposed framework with a traditional saturation recovery fast low-angle shot (turboFLASH) sequence and also without saturation recovery as a steady-state spoiled gradient echo (SPGR) sequence on animal and human studies.

METHODS

Simulations and phantom studies were performed for both the turboFLASH and the SPGR radial CAIPI methods. Ungated undersampled golden ratio radial CAIPI data with saturation recovery were acquired in 8 dogs and 2 human subjects. The CAIPI data without saturation pulses were acquired in 4 human subjects. For both methods, slice acceleration factors of two and three were used. A new spatio-temporal reconstruction using total variation and patch-based low rank constraints was used to jointly reconstruct the multi-slice multi-coil images.

RESULTS

Phantom scans and computer simulations showed that ungated SPGR generally provides better contrast to noise ratio (CNR) than the saturation recovery sequence if the saturation recovery time is less than 100ms. Both of the ungated radial CAIPI methods demonstrated promising image quality in terms of preserving dynamics of the contrast agent and maintaining anatomical structures, even with three slices acquired simultaneously.

CONCLUSION

Ungated simultaneous multi-slice acquisitions with either a saturation recovery turboFLASH sequence or a steady-state gradient echo SPGR sequence are feasible and provide increased slice coverage without loss of temporal resolution. Compared with a sensitivity encoding (SENSE) SMS reconstruction, the constrained reconstruction method provides better image quality for undersampled radial CAIPI data.

Quantitative three-dimensional myocardial perfusion cardiovascular magnetic resonance with accurate two-dimensional arterial input function assessment

BACKGROUND

Quantification of myocardial perfusion from first-pass cardiovascular magnetic resonance (CMR) images at high contrast agent (CA) dose requires separate acquisition of blood pool and myocardial tissue enhancement. In this study, a dual-sequence approach interleaving 2D imaging of the arterial input function with high-resolution 3D imaging for myocardial perfusion assessment is presented and validated for low and high CA dose.

METHODS

A dual-sequence approach interleaving 2D imaging of the aortic root and 3D imaging of the whole left ventricle using highly accelerated k-t PCA was implemented. Rest perfusion imaging was performed in ten healthy volunteers after administration of a Gadolinium-based CA at low (0.025 mmol/kg b.w.) and high dose (0.1 mmol/kg b.w.). Arterial input functions extracted from the 2D and 3D images were analysed for both doses. Myocardial contrast-to-noise ratios (CNR) were compared across volunteers and doses. Variations of myocardial perfusion estimates between volunteers and across myocardial territories were studied.

RESULTS

High CA dose imaging resulted in strong non-linearity of the arterial input function in the 3D images at peak CA concentration, which was avoided when the input function was derived from the 2D images. Myocardial CNR was significantly increased at high dose compared to low dose, with a 2.6-fold mean CNR gain. Most robust myocardial blood flow estimation was achieved using the arterial input function extracted from the 2D image at high CA dose. In this case, myocardial blood flow estimates varied by 24% between volunteers and by 20% between myocardial territories when analysed on a per-volunteer basis.

CONCLUSION

Interleaving 2D imaging for arterial input function assessment enables robust quantitative 3D myocardial perfusion imaging at high CA dose.

Evaluation of optimized breath-hold and free-breathing 3D ultrashort echo time contrast agent-free MRI of the human lung

PURPOSE

To evaluate an optimized stack of radials ultrashort echo time (UTE) 3D magnetic resonance imaging (MRI) sequence for breath-hold and free-breathing imaging of the human lung.

MATERIALS AND METHODS

A 3D stack of ultrashort echo time radials trajectory was optimized for coronal and axial lower-resolution breath-hold and higher-resolution free-breathing scans using Bloch simulations. The sequence was evaluated in 10 volunteers, without the use of contrast agents. Signal-to-noise ratio (SNR) mean and 95% confidence interval (CI) were determined from separate signal and noise images in a semiautomated fashion. The four scanning schemes were evaluated for significant differences in image quality using Student's t-test. Ten clinical patients were scanned with the sequence and findings were compared with concomitant computed tomography (CT) in nine patients. Breath-hold 3D spokes images were compared with 3D stack of radials in five volunteers. A Mann-Whitney U-test was performed to test significance in both cases.

RESULTS

Breath-hold imaging of the entire lung in volunteers was performed with SNR (mean = 42.5 [CI]: 35.5-49.5; mean = 34.3 [CI]: 28.6-40) in lung parenchyma for coronal and axial scans, respectively, which can be used as a quick scout scan. Longer respiratory triggered free-breathing scan enabled high-resolution UTE scanning with mean SNR of 14.2 ([CI]: 12.9-15.5) and 9.2 ([CI]: 8.2-10.2) for coronal and axial scans, respectively. Axial free-breathing scans showed significantly higher image quality (P = 0.008) than the three other scanning schemes. The mean score for comparison with CT was 1.67 (score 0: n = 0; 1: n = 3; 2: n = 6). There was no significant difference between CT and MRI (P = 0.25). 3D stack of radials images were significantly better than 3D spokes images (P < 0.001).

CONCLUSION

The optimized 3D stack of radials trajectory was shown to provide high-quality MR images of the lung parenchyma without the use of MRI contrast agents. The sequence may offer the possibility of breath-hold imaging and provides greater flexibility in trading off slice thickness and parallel imaging for scan time.

A review of 3D first-pass, whole-heart, myocardial perfusion cardiovascular magnetic resonance

A comprehensive review is undertaken of the methods available for 3D whole-heart first-pass perfusion (FPP) and their application to date, with particular focus on possible acceleration techniques. Following a summary of the parameters typically desired of 3D FPP methods, the review explains the mechanisms of key acceleration techniques and their potential use in FPP for attaining 3D acquisitions. The mechanisms include rapid sequences, non-Cartesian k-space trajectories, reduced k-space acquisitions, parallel imaging reconstructions and compressed sensing. An attempt is made to explain, rather than simply state, the varying methods with the hope that it will give an appreciation of the different components making up a 3D FPP protocol. Basic estimates demonstrating the required total acceleration factors in typical 3D FPP cases are included, providing context for the extent that each acceleration method can contribute to the required imaging speed, as well as potential limitations in present 3D FPP literature. Although many 3D FPP methods are too early in development for the type of clinical trials required to show any clear benefit over current 2D FPP methods, the review includes the small but growing quantity of clinical research work already using 3D FPP, alongside the more technical work. Broader challenges concerning FPP such as quantitative analysis are not covered, but challenges with particular impact on 3D FPP methods, particularly with regards to motion effects, are discussed along with anticipated future work in the field.

Comparison of centric and reverse-centric trajectories for highly accelerated three-dimensional saturation recovery cardiac perfusion imaging

PURPOSE

Highly undersampled three-dimensional (3D) saturation-recovery sequences are affected by k-space trajectory since the magnetization does not reach steady state during the acquisition and the slab excitation profile yields different flip angles in different slices. This study compares centric and reverse-centric 3D cardiac perfusion imaging.

METHODS

An undersampled (98 phase encodes) 3D ECG-gated saturation-recovery sequence that alternates centric and reverse-centric acquisitions each time frame was used to image phantoms and in vivo subjects. Flip angle variation across the slices was measured, and contrast with each trajectory was analyzed via Bloch simulation.

RESULTS

Significant variations in flip angle were observed across slices, leading to larger signal variation across slices for the centric acquisition. In simulation, severe transient artifacts were observed when using the centric trajectory with higher flip angles, placing practical limits on the maximum flip angle used. The reverse-centric trajectory provided less contrast, but was more robust to flip angle variations.

CONCLUSION

Both of the k-space trajectories can provide reasonable image quality. The centric trajectory can have higher CNR, but is more sensitive to flip angle variation. The reverse-centric trajectory is more robust to flip angle variation.

Accelerated cardiac MR stress perfusion with radial sampling after physical exercise with an MR-compatible supine bicycle ergometer

PURPOSE

To evaluate the feasibility of accelerated cardiac MR (CMR) perfusion with radial sampling using nonlinear image reconstruction after exercise on an MR-compatible supine bike ergometer.

METHODS

Eight healthy subjects were scanned on two separate days using radial and Cartesian CMR perfusion sequences in rest and exercise stress perfusion. Four different methods (standard gridding, conjugate gradient SENSE [CG-SENSE], nonlinear inversion with joint estimation of coil-sensitivity profiles [NLINV] and compressed sensing with a total variation constraint [TV]) were compared for the reconstruction of radial data. Cartesian data were reconstructed using SENSE. All images were assessed by two blinded readers in terms of image quality and diagnostic value.

RESULTS

CG-SENSE and NLINV were scored more favorably than TV (in both rest and stress perfusion cases, P < 0.05) and gridding (for rest perfusion cases, P < 0.05). TV images showed patchy artifacts, which negatively influenced image quality especially in the stress perfusion images acquired with a low number of radial spokes. Although CG-SENSE and NLINV received better scores than Cartesian sampling in both rest and exercise stress perfusion cases, these differences were not statistically significant (P > 0.05).

CONCLUSION

We have demonstrated the feasibility of accelerated CMR perfusion using radial sampling after physical exercise using a supine bicycle ergometer in healthy subjects. For reconstruction of undersampled radial perfusion, CG-SENSE and NLINV resulted in better image quality than standard gridding or TV reconstruction. Further technical improvements and clinical assessment are needed before using this approach in patients with suspected coronary artery disease.

Quantitative three-dimensional cardiovascular magnetic resonance myocardial perfusion imaging in systole and diastole

BACKGROUND

Two-dimensional (2D) perfusion cardiovascular magnetic resonance (CMR) remains limited by a lack of complete myocardial coverage. Three-dimensional (3D) perfusion CMR addresses this limitation and has recently been shown to be clinically feasible. However, the feasibility and potential clinical utility of quantitative 3D perfusion measurements, as already shown with 2D-perfusion CMR and positron emission tomography, has yet to be evaluated. The influence of systolic or diastolic acquisition on myocardial blood flow (MBF) estimates, diagnostic accuracy and image quality is also unknown for 3D-perfusion CMR. The purpose of this study was to establish the feasibility of quantitative 3D-perfusion CMR for the detection of coronary artery disease (CAD) and to compare systolic and diastolic estimates of MBF.

METHODS

Thirty-five patients underwent 3D-perfusion CMR with data acquired at both end-systole and mid-diastole. MBF and myocardial perfusion reserve (MPR) were estimated on a per patient and per territory basis by Fermi-constrained deconvolution. Significant CAD was defined as stenosis ≥70% on quantitative coronary angiography.

RESULTS

Twenty patients had significant CAD (involving 38 out of 105 territories). Stress MBF and MPR had a high diagnostic accuracy for the detection of CAD in both systole (area under curve [AUC]: 0.95 and 0.92, respectively) and diastole (AUC: 0.95 and 0.94). There were no significant differences in the AUCs between systole and diastole (p values >0.05). At stress, diastolic MBF estimates were significantly greater than systolic estimates (no CAD: 3.21 ± 0.50 vs. 2.75 ± 0.42 ml/g/min, p < 0.0001; CAD: 2.13 ± 0.45 vs. 1.98 ± 0.41 ml/g/min, p < 0.0001); but at rest, there were no significant differences (p values >0.05). Image quality was higher in systole than diastole (median score 3 vs. 2, p = 0.002).

CONCLUSIONS

Quantitative 3D-perfusion CMR is feasible. Estimates of MBF are significantly different for systole and diastole at stress but diagnostic accuracy to detect CAD is high for both cardiac phases. Better image quality suggests that systolic data acquisition may be preferable.

Three-dimensional dynamic contrast enhanced imaging of the carotid artery with direct arterial input function measurement

PURPOSE

Kinetic analysis using dynamic contrast enhanced MRI to assess neovascularization of carotid plaque requires images with high spatial and temporal resolution. This work demonstrates a new three-dimensional (3D) dynamic contrast enhanced imaging sequence, which directly measures the arterial input function with high temporal resolution yet maintains the high spatial resolution required to identify areas of increased adventitial neovascularity.

THEORY AND METHODS

The sequence consists of multiple rapid acquisitions of a saturation prepared dynamic 3D gradient recalled echo (GRE) sequence temporally interleaved with multiple acquisitions of a 2D slice. The saturation recovery time was adjusted to maintain signal linearity with the very different contrast agent concentrations in the 2D slice and 3D volume. The K(trans) maps were obtained from the 3D dynamic contrast measurements while the 2D slice was used to obtain the arterial input function. Calibration and dynamic studies are presented.

RESULTS

For contrast agent concentrations up to 5 mM, a saturation recovery time for the 2D slice of 20 ms resulted in less than a 10% deviation from the desired linear response of signal intensity with contrast agent concentration. The corresponding saturation recovery time of 83 ms for the 3D volume maintained less than a 10% deviation from the linear response up to contrast agent concentrations of 2 mM while a contrast agent concentration of 5 mM had almost a 30% deviation. There was a significant improvement in signal attenuation (9 ± 3% versus 23 ± 5% at 40 cm/s) when flow compensation was added to the slice select gradients. For patient studies, volume transfer and plasma fraction maps were calculated with data from the proposed sequence.

CONCLUSION

This work demonstrated a novel sequence for 3D dynamic contrast enhanced imaging with a simultaneously acquired 2D slice that directly measures the arterial input function with high temporal resolution. Acquisition parameters can be adjusted to accommodate the full range of contrast agent concentration values to be encountered and the kinetic parameters obtained were consistent with expected values.

Localized spatio-temporal constraints for accelerated CMR perfusion

PURPOSE

To develop and evaluate an image reconstruction technique for cardiac MRI (CMR) perfusion that uses localized spatio-temporal constraints.

METHODS

CMR perfusion plays an important role in detecting myocardial ischemia in patients with coronary artery disease. Breath-hold k-t-based image acceleration techniques are typically used in CMR perfusion for superior spatial/temporal resolution and improved coverage. In this study, we propose a novel compressed sensing-based image reconstruction technique for CMR perfusion, with applicability to free-breathing examinations. This technique uses local spatio-temporal constraints by regularizing image patches across a small number of dynamics. The technique was compared with conventional dynamic-by-dynamic reconstruction, and sparsity regularization using a temporal principal-component (pc) basis, as well as zero-filled data in multislice two-dimensional (2D) and three-dimensional (3D) CMR perfusion. Qualitative image scores were used (1 = poor, 4 = excellent) to evaluate the technique in 3D perfusion in 10 patients and five healthy subjects. On four healthy subjects, the proposed technique was also compared with a breath-hold multislice 2D acquisition with parallel imaging in terms of signal intensity curves.

RESULTS

The proposed technique produced images that were superior in terms of spatial and temporal blurring compared with the other techniques, even in free-breathing datasets. The image scores indicated a significant improvement compared with other techniques in 3D perfusion (x-pc regularization, 2.8 ± 0.5 versus 2.3 ± 0.5; dynamic-by-dynamic, 1.7 ± 0.5; zero-filled, 1.1 ± 0.2). Signal intensity curves indicate similar dynamics of uptake between the proposed method with 3D acquisition and the breath-hold multislice 2D acquisition with parallel imaging.

CONCLUSION

The proposed reconstruction uses sparsity regularization based on localized information in both spatial and temporal domains for highly accelerated CMR perfusion with potential use in free-breathing 3D acquisitions.

Steady-state first-pass perfusion (SSFPP): a new approach to 3D first-pass myocardial perfusion imaging

PURPOSE

To describe and characterize a new approach to first-pass myocardial perfusion utilizing balanced steady-state free precession acquisition without the use of saturation recovery or other magnetization preparation.

THEORY

The balanced steady-state free precession sequence is inherently sensitive to contrast agent enhancement of the myocardium. This sensitivity can be used to advantage in first-pass myocardial perfusion imaging by eliminating the need for magnetization preparation.

METHODS

Bloch equation simulations, phantom experiments, and in vivo 2D imaging studies were run comparing the proposed technique with three other methods: saturation recovery spoiled gradient echo, saturation recovery steady-state free precession, and steady-state spoiled gradient echo without magnetization preparation. Additionally, an acquisition-reconstruction strategy for 3D perfusion imaging is proposed and initial experience with this approach is demonstrated in healthy subjects and one patient.

RESULTS

Phantom experiments verified simulation results showing the sensitivity of the balanced steady-state free precession sequence to contrast agent enhancement in solid tissue is similar to that of magnetization-prepared acquisitions. Images acquired in normal volunteers showed the proposed technique provided superior signal and signal-to-noise ratio compared with all other sequences at baseline as well as postcontrast.

CONCLUSIONS

A new approach to first-pass myocardial perfusion is presented that obviates the need for magnetization preparation and provides high signal-to-noise ratio.