GRASP-Pro: imProving GRASP DCE-MRI through self-calibrating subspace-modeling and contrast phase automation

Technical Note: Swing golden angle - A navigator-interleaved golden angle trajectory with eddy current suppression - Application in free-running cardiac MRI

BACKGROUND

Golden angle (GA) radial trajectory is advantageous for dynamic magnetic resonance imaging (MRI). Recently, several advanced algorithms have been developed based on navigator-interleaved GA trajectory to realize free-running cardiac MRI. However, navigator-interleaved GA trajectory suffers from the eddy-current effect, which reduces the image quality.

PURPOSE

This work aims to integrate the navigator-interleaved GA trajectory with clinical cardiac MRI acquisition, with the minimum eddy-current artifacts. The ultimate goal is to realize a high-quality free-running cardiac imaging technique.

METHODS

In this paper, we propose a new "swing golden angle" (swingGA) radial profile order. SwingGA samples the k-space by rotating back and forth at the generalized golden ratio interval, with smoothly interleaved navigator readouts. The sampling efficiency and angle increment distributions were investigated by numerical simulations. Static phantom imaging experiments were conducted to evaluate the eddy current effect, compared with cartesian, golden angle radial (GA), and tiny golden angle (tGA) trajectories. Furthermore, 12 heart-healthy subjects (aged 21-25 years) were recruited for free-running cardiac imaging with different sampling trajectories. Dynamic images were reconstructed by a low-rank subspace-constrained algorithm. The image quality was evaluated by signal-to-noise-ratio and spectrum analysis in the heart region, and compared with traditional clinical cardiac MRI images.

RESULTS

SwingGA pattern achieves the highest sampling efficiency (mSE > 0.925) and the minimum azimuthal angle increment (mAD < 1.05). SwingGA can effectively suppress eddy currents in static phantom images, with the lowest normalized root mean square error (nRMSE) values among radial trajectories. For the in-vivo cardiac images, swingGA enjoys the highest SNR both in the blood pool and myocardium, and contains the minimum level of high-frequency artifacts. The free-running cardiac images have good consistency with traditional clinical cardiac MRI, and the swingGA sampling pattern achieves the best image quality among all sampling patterns.

CONCLUSIONS

The proposed swingGA sampling pattern can effectively improve the sampling efficiency and reduce the eddy currents for the navigator-interleaved GA sequence. SwingGA is a promising sampling pattern for free-running cardiac MRI.

Blind deconvolution decreases requirements on temporal resolution of DCE-MRI: Application to 2nd generation pharmacokinetic modeling

PURPOSE

Dynamic Contrast-Enhanced (DCE) MRI with 2nd generation pharmacokinetic models provides estimates of plasma flow and permeability surface-area product in contrast to the broadly used 1st generation models (e.g. the Tofts models). However, the use of 2nd generation models requires higher frequency with which the dynamic images are acquired (around 1.5 s per image). Blind deconvolution can decrease the demands on temporal resolution as shown previously for one of the 1st generation models. Here, the temporal-resolution requirements achievable for blind deconvolution with a 2nd generation model are studied.

METHODS

The 2nd generation model is formulated as the distributed-capillary adiabatic-tissue-homogeneity (DCATH) model. Blind deconvolution is based on Parker's model of the arterial input function. The accuracy and precision of the estimated arterial input functions and the perfusion parameters is evaluated on synthetic and real clinical datasets with different levels of the temporal resolution.

RESULTS

The estimated arterial input functions remained unchanged from their reference high-temporal-resolution estimates (obtained with the sampling interval around 1 s) when increasing the sampling interval up to about 5 s for synthetic data and up to 3.6-4.8 s for real data. Further increasing of the sampling intervals led to systematic distortions, such as lowering and broadening of the 1st pass peak. The resulting perfusion-parameter estimation error was below 10% for the sampling intervals up to 3 s (synthetic data), in line with the real data perfusion-parameter boxplots which remained unchanged up to the sampling interval 3.6 s.

CONCLUSION

We show that use of blind deconvolution decreases the demands on temporal resolution in DCE-MRI from about 1.5 s (in case of measured arterial input functions) to 3-4 s. This can be exploited in increased spatial resolution or larger organ coverage.

High-resolution free-breathing hepatobiliary phase MRI of the liver using XD-GRASP

PURPOSE

To evaluate the performance of high-resolution free-breathing (FB) hepatobiliary phase imaging of the liver using the eXtra-Dimension Golden-angle RAdial Sparse Parallel (XD-GRASP) MRI technique.

METHODS

Fifty-eight clinical patients (41 males, mean age = 52.9 ± 12.9) with liver lesions who underwent dynamic contrast-enhanced MRI with a liver-specific contrast agent were prospectively recruited for this study. Both breath-hold volumetric interpolated examination (BH-VIBE) imaging and FB imaging were performed during the hepatobiliary phase. FB images were acquired using a stack-of-stars golden-angle radial sequence and were reconstructed using the XD-GRASP method. Two experienced radiologists blinded to acquisition schemes independently scored the overall image quality, liver edge sharpness, hepatic vessel clarity, conspicuity of lesion, and overall artifact level of each image. The non-parametric paired two-tailed Wilcoxon signed-rank test was used for statistical analysis.

RESULTS

Compared to BH-VIBE images, XD-GRASP images received significantly higher scores (P < 0.05) for the liver edge sharpness (4.83 ± 0.45 vs 4.29 ± 0.46), the hepatic vessel clarity (4.64 ± 0.67 vs 4.15 ± 0.56) and the conspicuity of lesion (4.75 ± 0.53 vs 4.31 ± 0.50). There were no significant differences (P > 0.05) between BH-VIBE and XD-GRASP images for the overall image quality (4.61 ± 0.50 vs 4.74 ± 0.47) and the overall artifact level (4.13 ± 0.44 vs 4.05 ± 0.61).

CONCLUSION

Compared to conventional BH-VIBE MRI, FB radial acquisition combined with XD-GRASP reconstruction facilitates higher spatial resolution imaging of the liver during the hepatobiliary phase. This enhancement can significantly improve the visualization and evaluation of the liver.

Dynamic Radial MR Imaging for Endoleak Surveillance after Endovascular Repair of Abdominal Aortic Aneurysms with Inconclusive CT Angiography: A Prospective Study

Background/Objectives: To assess free-breathing, dynamic radial magnetic resonance angiography (MRA) for detecting endoleaks post-endovascular aortic repair (EVAR) in cases with inconclusive computed tomography angiography (CTA). Methods: This prospective single-center study included 17 participants (mean age, 70 ± 9 years; 13 males) who underwent dynamic radial MRI (Golden-angle RAdial Sparse Parallel-Volumetric Interpolated BrEath-hold, GRASP-VIBE) after inconclusive multiphasic CT for the presence of endoleaks during the follow-up of EVAR-treated abdominal aortic aneurysms. CT and MRI datasets were independently assessed by two radiologists for image quality, diagnostic confidence, and the presence/type of endoleak. Statistical analyses included interrater and intermethod agreement, and diagnostic performance (sensitivity, specificity, area under the curve (AUC)). Results: Subjective image analysis demonstrated good image quality and interrater agreement (k ≥ 0.6) for both modalities, while diagnostic confidence was significantly higher in MRA (p = 0.03). There was significantly improved accuracy for detecting type II endoleaks on MRA (AUC 0.97 [95% CI: 0.87, 1.0]) compared to CTA (AUC 0.66 [95% CI: 0.41, 0.91]; p = 0.03). Although MRA demonstrated higher values for sensitivity, specificity, AUC, and interrater agreement, none of the other types nor the overall detection rate for endoleaks showed differences in the diagnostic performance over CT (p ≥ 0.12). CTA and MRA revealed slight to moderate intermethod concordance in endoleak detection (k = 0.3-0.64). Conclusions: The GRASP-VIBE MRA characterized by high spatial and temporal resolution demonstrates clinical feasibility with good image quality and superior diagnostic confidence. It notably enhances diagnostic performance in detecting and classifying endoleaks, particularly type II, compared to traditional multiphase CTA with inconclusive findings.

Feasibility of measuring blood-brain barrier permeability using ultra-short echo time radial magnetic resonance imaging

BACKGROUND AND PURPOSE

The purpose of this study is to evaluate the feasibility of using 3-dimensional (3D) ultra-short echo time (UTE) radial imaging method for measurement of the permeability of the blood-brain barrier (BBB) to gadolinium-based contrast agent. In this study, we propose to use the golden-angle radial sparse parallel (GRASP) method with 3D center-out trajectories for UTE, hence named as 3D UTE-GRASP. We first examined the feasibility of using 3D UTE-GRASP dynamic contrast-enhanced (DCE)-magnetic resonance imaging (MRI) for differentiating subtle BBB disruptions induced by focused ultrasound (FUS). Then, we examined the BBB permeability changes in Alzheimer's disease (AD) pathology using Alzheimer's disease transgenic mice (5xFAD) at different ages.

METHODS

For FUS experiments, we used four Sprague Dawley rats at similar ages where we compared BBB permeability of each rat receiving the FUS sonication with different acoustic power (0.4-1.0 MPa). For AD transgenic mice experiments, we included three 5xFAD mice (6, 12, and 16 months old) and three wild-type mice (4, 8, and 12 months old).

RESULTS

The result from FUS experiments showed a progressive increase in BBB permeability with increase of acoustic power (p < .05), demonstrating the sensitivity of DCE-MRI method for detecting subtle changes in BBB disruption. Our AD transgenic mice experiments suggest an early BBB disruption in 5xFAD mice, which is further impaired with aging.

CONCLUSION

The results in this study substantiate the feasibility of using the proposed 3D UTE-GRASP method for detecting subtle BBB permeability changes expected in neurodegenerative diseases, such as AD.

Incorporation of view sharing and KWIC filtering into GRASP-Pro improves spatial resolution of single-shot, multi-TI, late gadolinium enhancement MRI

While single-shot late gadolinium enhancement (LGE) is useful for imaging patients with arrhythmia and/or dyspnea, it produces low spatial resolution. One approach to improve spatial resolution is to accelerate data acquisition using compressed sensing (CS). Our previous work described a single-shot, multi-inversion time (TI) LGE pulse sequence using radial k-space sampling and CS, but over-regularization resulted in significant image blurring that muted the benefits of data acceleration. The purpose of the present study was to improve the spatial resolution of the single-shot, multi-TI LGE pulse sequence by incorporating view sharing (VS) and k-space weighted contrast (KWIC) filtering into a GRASP-Pro reconstruction. In 24 patients (mean age = 61 ± 16 years; 9/15 females/males), we compared the performance of our improved multi-TI LGE and standard multi-TI LGE, where clinical standard LGE was used as a reference. Two clinical raters independently graded multi-TI images and clinical LGE images visually on a five-point Likert scale (1, nondiagnostic; 3, clinically acceptable; 5, best) for three categories: the conspicuity of myocardium or scar, artifact, and noise. The summed visual score (SVS) was defined as the sum of the three scores. Myocardial scar volume was quantified using the full-width at half-maximum method. The SVS was not significantly different between clinical breath-holding LGE (median 13.5, IQR 1.3) and multi-TI LGE (median 12.5, IQR 1.6) (P = 0.068). The myocardial scar volumes measured from clinical standard LGE and multi-TI LGE were strongly correlated (coefficient of determination, R2  = 0.99) and in good agreement (mean difference = 0.11%, lower limit of the agreement = -2.13%, upper limit of the agreement = 2.34%). The inter-rater agreement in myocardial scar volume quantification was strong (intraclass correlation coefficient = 0.79). The incorporation of VS and KWIC into GRASP-Pro improved spatial resolution. Our improved 25-fold accelerated, single-shot LGE sequence produces clinically acceptable image quality, multi-TI reconstruction, and accurate myocardial scar volume quantification.

Intravoxel incoherent motion magnetic resonance imaging reconstruction from highly under-sampled diffusion-weighted PROPELLER acquisition data via physics-informed residual feedback unrolled network

Objective. The acquisition of diffusion-weighted images for intravoxel incoherent motion (IVIM) imaging is time consuming. This work aims to accelerate the scan through a highly under-sampling diffusion-weighted turbo spin echo PROPELLER (DW-TSE-PROPELLER) scheme and to develop a reconstruction method for accurate IVIM parameter mapping from the under-sampled data.Approach.The proposed under-sampling DW-TSE-PROPELLER scheme for IVIM imaging is that a few blades perb-value are acquired and rotated along theb-value dimension to cover high-frequency information. A physics-informed residual feedback unrolled network (PIRFU-Net) is proposed to directly estimate distortion-free and artifact-free IVIM parametric maps (i.e., the perfusion-free diffusion coefficientDand the perfusion fractionf) from highly under-sampled DW-TSE-PROPELLER data. PIRFU-Net used an unrolled convolution network to explore data redundancy in the k-q space to remove under-sampling artifacts. An empirical IVIM physical constraint was incorporated into the network to ensure that the signal evolution curves along theb-value follow a bi-exponential decay. The residual between the realistic and estimated measurements was fed into the network to refine the parametric maps. Meanwhile, the use of synthetic training data eliminated the need for genuine DW-TSE-PROPELLER data.Main results.The experimental results show that the DW-TSE-PROPELLER acquisition was six times faster than full k-space coverage PROPELLER acquisition and within a clinically acceptable time. Compared with the state-of-the-art methods, the distortion-freeDandfmaps estimated by PIRFU-Net were more accurate and had better-preserved tissue boundaries on a simulated human brain and realistic phantom/rat brain/human brain data.Significance.Our proposed method greatly accelerates IVIM imaging. It is capable of directly and simultaneously reconstructing distortion-free, artifact-free, and accurateDandfmaps from six-fold under-sampled DW-TSE-PROPELLER data.

Rapid 3D T1 mapping using deep learning-assisted Look-Locker inversion recovery MRI

PURPOSE

Conventional 3D Look-Locker inversion recovery (LLIR) T1 mapping requires multi-repetition data acquisition to reconstruct images at different inversion times for T1 fitting. To ensure B1 robustness, sufficient time of delay (TD) is needed between repetitions, which prolongs scan time. This work proposes a novel deep learning-assisted LLIR MRI approach for rapid 3D T1 mapping without TD.

THEORY AND METHODS

The proposed approach is based on the fact thatT 1 * $$ {\mathrm{T}}_1^{\ast } $$ , the effective T1 in LLIR imaging, is independent of TD and can be estimated from both LLIR imaging with and without TD, while accurate conversion ofT 1 * $$ {\mathrm{T}}_1^{\ast } $$ to T1 requires TD. Therefore, deep learning can be used to learn the conversion ofT 1 * $$ {\mathrm{T}}_1^{\ast } $$ to T1 , which eliminates the need for TD. This idea was implemented for inversion-recovery-prepared Golden-angel RAdial Sparse Parallel T1 mapping (GraspT1 ). 39 GraspT1 datasets with a TD of 6 s (GraspT1 -TD6) were used for training, which also incorporates additional anatomical images. The trained network was applied for T1 estimation in 14 GraspT1 datasets without TD (GraspT1 -TD0). The robustness of the trained network was also tested.

RESULTS

Deep learning-based T1 estimation from GraspT1 -TD0 is accurate compared to the reference. Incorporation of additional anatomical images improves the accuracy of T1 estimation. The technique is also robust against slight variation in spatial resolution, imaging orientation and scanner platform.

CONCLUSION

Our approach eliminates the need for TD in 3D LLIR imaging without affecting the T1 estimation accuracy. It represents a novel use of deep learning towards more efficient and robust 3D LLIR T1 mapping.

Latent signal models: Learning compact representations of signal evolution for improved time-resolved, multi-contrast MRI

PURPOSE

To improve time-resolved reconstructions by training auto-encoders to learn compact representations of Bloch-simulated signal evolution and inserting the decoder into the forward model.

METHODS

Building on model-based nonlinear and linear subspace techniques, we train auto-encoders on dictionaries of simulated signal evolution to learn compact, nonlinear, latent representations. The proposed latent signal model framework inserts the decoder portion of the auto-encoder into the forward model and directly reconstructs the latent representation. Latent signal models essentially serve as a proxy for fast and feasible differentiation through the Bloch equations used to simulate signal. This work performs experiments in the context of T₂ -shuffling, gradient echo EPTI, and MPRAGE-shuffling. We compare how efficiently auto-encoders represent signal evolution in comparison to linear subspaces. Simulation and in vivo experiments then evaluate if reducing degrees of freedom by incorporating our proxy for the Bloch equations, the decoder portion of the auto-encoder, into the forward model improves reconstructions in comparison to subspace constraints.

RESULTS

An auto-encoder with 1 real latent variable represents single-tissue fast spin echo, EPTI, and MPRAGE signal evolution to within 0.15% normalized RMS error, enabling reconstruction problems with 3 degrees of freedom per voxel (real latent variable + complex scaling) in comparison to linear models with 4-8 degrees of freedom per voxel. In simulated/in vivo T₂ -shuffling and in vivo EPTI experiments, the proposed framework achieves consistent quantitative normalized RMS error improvement over linear approaches. From qualitative evaluation, the proposed approach yields images with reduced blurring and noise amplification in MPRAGE-shuffling experiments.

CONCLUSION

Directly solving for nonlinear latent representations of signal evolution improves time-resolved MRI reconstructions.

Kz-accelerated variable-density stack-of-stars MRI

This work aimed to develop a modified stack-of-stars golden-angle radial sampling scheme with variable-density acceleration along the slice (kz) dimension (referred to as VD-stack-of-stars) and to test this new sampling trajectory with multi-coil compressed sensing reconstruction for rapid motion-robust 3D liver MRI. VD-stack-of-stars sampling implements additional variable-density undersampling along the kz dimension, so that slice resolution (or volumetric coverage) can be increased without prolonging scan time. The new sampling trajectory (with increased slice resolution) was compared with standard stack-of-stars sampling with fully sampled kz (with standard slice resolution) in both non-contrast-enhanced free-breathing liver MRI and dynamic contrast-enhanced MRI (DCE-MRI) of the liver in volunteers. For both sampling trajectories, respiratory motion was extracted from the acquired radial data, and images were reconstructed using motion-compensated (respiratory-resolved or respiratory-weighted) dynamic radial compressed sensing reconstruction techniques. Qualitative image quality assessment (visual assessment by experienced radiologists) and quantitative analysis (as a metric of image sharpness) were performed to compare images acquired using the new and standard stack-of-stars sampling trajectories. Compared to standard stack-of-stars sampling, both non-contrast-enhanced and DCE liver MR images acquired with VD-stack-of-stars sampling presented improved overall image quality, sharper liver edges and increased hepatic vessel clarity in all image planes. The results have suggested that the proposed VD-stack-of-stars sampling scheme can achieve improved performance (increased slice resolution or volumetric coverage with better image quality) over standard stack-of-stars sampling in free-breathing DCE-MRI without increasing scan time. The reformatted coronal and sagittal images with better slice resolution may provide added clinical value.

GRASPNET: Fast spatiotemporal deep learning reconstruction of golden-angle radial data for free-breathing dynamic contrast-enhanced magnetic resonance imaging

The purpose of the current study was to develop a deep learning technique called Golden-angle RAdial Sparse Parallel Network (GRASPnet) for fast reconstruction of dynamic contrast-enhanced 4D MRI acquired with golden-angle radial k-space trajectories. GRASPnet operates in the image-time space and does not use explicit data consistency to minimize the reconstruction time. Three different network architectures were developed: (1) GRASPnet-2D: 2D convolutional kernels (x,y) and coil and contrast dimensions collapsed into a single combined dimension; (2) GRASPnet-3D: 3D kernels (x,y,t); and (3) GRASPnet-2D + time: two 3D kernels to first exploit spatial correlations (x,y,1) followed by temporal correlations (1,1,t). The networks were trained using iterative GRASP reconstruction as the reference. Free-breathing 3D abdominal imaging with contrast injection was performed on 33 patients with liver lesions using a T1-weighted golden-angle stack-of-stars pulse sequence. Ten datasets were used for testing. The three GRASPnet architectures were compared with iterative GRASP results using quantitative and qualitative analysis, including impressions from two body radiologists. The three GRASPnet techniques reduced the reconstruction time to about 13 s with similar results with respect to iterative GRASP. Among the GRASPnet techniques, GRASPnet-2D + time compared favorably in the quantitative analysis. Spatiotemporal deep learning enables reconstruction of dynamic 4D contrast-enhanced images in a few seconds, which would facilitate translation to clinical practice of compressed sensing methods that are currently limited by long reconstruction times.

Deep learning-based super-resolution gradient echo imaging of the pancreas: Improvement of image quality and reduction of acquisition time

PURPOSE

The purpose of this study was to evaluate the impact of a deep learning-based super-resolution technique on T1-weighted gradient-echo acquisitions (volumetric interpolated breath-hold examination; VIBE) on the assessment of pancreatic MRI at 1.5 T compared to standard VIBE imaging (VIBESTD).

MATERIALS AND METHODS

This retrospective single-center study was conducted between April 2021 and October 2021. Fifty patients with a total of 50 detectable pancreatic lesion entities were included in this study. There were 27 men and 23 women, with a mean age of 69 ± 13 (standard deviation [SD]) years (age range: 33-89 years). VIBESTD (precontrast, dynamic, postcontrast) was retrospectively processed with a deep learning-based super-resolution algorithm including a more aggressive partial Fourier setting leading to a simulated acquisition time reduction (VIBESR). Image analysis was performed by two radiologists regarding lesion detectability, noise levels, sharpness and contrast of pancreatic edges, as well as regarding diagnostic confidence using a 5-point Likert-scale with 5 being the best.

RESULTS

VIBESR was rated better than VIBESTD by both readers regarding lesion detectability (5 [IQR: 5, 5] vs. 5 [IQR: 4, 5], for reader 1; 5 [IQR: 5, 5] vs. 4 [IQR: 4, 5]) for reader 2; both P <0.001), noise levels (5 [IQR: 5, 5] vs. 5 [IQR: 4, 5] for reader 1; 5 [IQR: 5, 5] vs. 4 [IQR: 4, 5] for reader 2; both P <0.001), sharpness and contrast of pancreatic edges (5 [IQR: 5, 5] vs. 5 [IQR: 4, 5] for reader 1; 5 [IQR: 5, 5] vs. 4 [IQR: 4, 5] for reader 2; both P <0.001), as well as regarding diagnostic confidence (5 [IQR: 5, 5] vs. 5 [IQR: 4, 5] for reader 1; 5 [IQR: 5, 5] vs. 4 [IQR: 4, 5] for reader 2; both P <0.001). There were no significant differences between lesion sizes as measured by the two readers on VIBESR and VIBESTD images (P > 0.05). The mean acquisition time for VIBESTD (15 ± 1 [SD] s; range: 11-16 s) was longer than that for VIBESR (13 ± 1 [SD] s; range: 11-14 s) (P < 0.001).

CONCLUSION

Our results indicate that the newly developed deep learning-based super-resolution algorithm adapted to partial Fourier acquisitions has a positive influence not only on shortening the examination time but also on improvement of image quality in pancreatic MRI.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Intraframe motion blurring, as a major challenge in free-breathing dynamic MRI, can be reduced if high temporal resolution can be achieved. To address this challenge, this work proposes a highly accelerated 4D (3D + time) dynamic MRI framework with subsecond temporal resolution that does not require explicit motion compensation. The method combines standard stack-of-stars golden-angle radial sampling and tailored GRASP-Pro (Golden-angle RAdial Sparse Parallel imaging with imProved performance) reconstruction. Specifically, 4D dynamic MRI acquisition is performed continuously without motion gating or sorting. The k-space centers in stack-of-stars radial data are organized to guide estimation of a temporal basis, with which GRASP-Pro reconstruction is employed to enforce joint low-rank subspace and sparsity constraints. This new basis estimation strategy is the new feature proposed for subspace-based reconstruction in this work to achieve high temporal resolution (e.g., subsecond/3D volume). It does not require sequence modification to acquire additional navigation data, it is compatible with commercially available stack-of-stars sequences, and it does not need an intermediate reconstruction step. The proposed 4D dynamic MRI approach was tested in abdominal motion phantom, free-breathing abdominal MRI, and dynamic contrast-enhanced MRI (DCE-MRI). Our results have shown that GRASP-Pro reconstruction with the new basis estimation strategy enables highly-accelerated 4D dynamic imaging at subsecond temporal resolution (with five spokes or less for each dynamic frame per image slice) for both free-breathing non-DCE-MRI and DCE-MRI. In the abdominal phantom, better image quality with lower root mean square error and higher structural similarity index was achieved using GRASP-Pro compared with standard GRASP. With the ability to acquire each 3D image in less than 1 s, intraframe respiratory blurring can be intrinsically reduced for body applications with our approach, which eliminates the need for explicit motion detection and motion compensation.

A Slice-Low-Rank Plus Sparse (slice-L + S) Reconstruction Method for k-t Undersampled Multiband First-Pass Myocardial Perfusion MRI

Purpose

The synergistic use of k‐t undersampling and multiband (MB) imaging has the potential to provide extended slice coverage and high spatial resolution for first‐pass perfusion MRI. The low‐rank plus sparse (L + S) model has shown excellent performance for accelerating single‐band (SB) perfusion MRI.

Methods

A MB data consistency method employing ESPIRiT maps and through‐plane coil information was developed. This data consistency method was combined with the temporal L + S constraint to form the slice‐L + S method. Slice‐L + S was compared to SB L + S and the sequential operations of split slice‐GRAPPA and SB L + S (seq‐SG‐L + S) using synthetic data formed from multislice SB images. Prospectively k‐t undersampled MB data were also acquired and reconstructed using seq‐SG‐L + S and slice‐L + S.

Results

Using synthetic data with total acceleration rates of 6–12, slice‐L + S outperformed SB L + S and seq‐SG‐L + S (N = 7 subjects) with respect to normalized RMSE and the structural similarity index (P < 0.05 for both). For the specific case with MB factor = 3 and rate 3 undersampling, or for SB imaging with rate 9 undersampling (N = 7 subjects), the normalized RMSE values were 0.037 ± 0.007, 0.042 ± 0.005, and 0.031 ± 0.004; and the structural similarity index values were 0.88 ± 0.03, 0.85 ± 0.03, and 0.89 ± 0.02 for SB L + S, seq‐SG‐L + S, and slice‐L + S, respectively (P < 0.05 for both). For prospectively undersampled MB data, slice‐L + S provided better image quality than seq‐SG‐L + S for rate 6 (N = 7) and rate 9 acceleration (N = 7) as scored by blinded experts.

Conclusion

Slice‐L + S outperformed SB‐L + S and seq‐SG‐L + S and provides 9 slice coverage of the left ventricle with a spatial resolution of 1.5 mm × 1.5 mm with good image quality.

Repeatability and robustness of MP-GRASP T1 mapping

PURPOSE

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain.

THEORY AND METHODS

Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC).

RESULTS

MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996).

CONCLUSION

MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

Estimation of the capillary level input function for dynamic contrast-enhanced MRI of the breast using a deep learning approach

PURPOSE

To develop a deep learning approach to estimate the local capillary-level input function (CIF) for pharmacokinetic model analysis of DCE-MRI.

METHODS

A deep convolutional network was trained with numerically simulated data to estimate the CIF. The trained network was tested using simulated lesion data and used to estimate voxel-wise CIF for pharmacokinetic model analysis of breast DCE-MRI data using an abbreviated protocol from women with malignant (n = 25) and benign (n = 28) lesions. The estimated parameters were used to build a logistic regression model to detect the malignancy.

RESULT

The pharmacokinetic parameters estimated using the network-predicted CIF from our breast DCE data showed significant differences between the malignant and benign groups for all parameters. Testing the diagnostic performance with the estimated parameters, the conventional approach with arterial input function (AIF) showed an area under the curve (AUC) between 0.76 and 0.87, and the proposed approach with CIF demonstrated similar performance with an AUC between 0.79 and 0.81.

CONCLUSION

This study shows the feasibility of estimating voxel-wise CIF using a deep neural network. The proposed approach could eliminate the need to measure AIF manually without compromising the diagnostic performance to detect the malignancy in the clinical setting.

Golden-Angle Radial MRI: Basics, Advances, and Applications

In recent years, golden‐angle radial sampling has received substantial attention and interest in the magnetic resonance imaging (MRI) community, and it has become a popular sampling trajectory for both research and clinical use. However, although the number of relevant techniques and publications has grown rapidly, there is still a lack of a review paper that provides a comprehensive overview and summary of the basics of golden‐angle rotation, the advantages and challenges/limitations of golden‐angle radial sampling, and recommendations in using different types of golden‐angle radial trajectories for MRI applications. Such a review paper is expected to be helpful both for clinicians who are interested in learning the potential benefits of golden‐angle radial sampling and for MRI physicists who are interested in exploring this research direction. The main purpose of this review paper is thus to present an overview and summary about golden‐angle radial MRI sampling. The review consists of three sections. The first section aims to answer basic questions such as: what is a golden angle; how is the golden angle calculated; why is golden‐angle radial sampling useful, and what are its limitations. The second section aims to review more advanced trajectories of golden‐angle radial sampling, including tiny golden‐angle rotation, stack‐of‐stars golden‐angle radial sampling, and three‐dimensional (3D) kooshball golden‐angle radial sampling. Their respective advantages and limitations and potential solutions to address these limitations are also discussed. Finally, the third section reviews MRI applications that can benefit from golden‐angle radial sampling and provides recommendations to readers who are interested in implementing golden‐angle radial trajectories in their MRI studies.

Dynamic pulmonary MRI using motion-state weighted motion-compensation (MostMoCo) reconstruction with ultrashort TE: A structural and functional study

PURPOSE

To improve the quality of structural images and the quantification of ventilation in free-breathing dynamic pulmonary MRI.

METHODS

A 3D radial ultrashort TE (UTE) sequence with superior-inferior navigators was used to acquire pulmonary data during free breathing. All acquired data were binned into different motion states according to the respiratory signal extracted from superior-inferior navigators. Motion-resolved images were reconstructed using eXtra-Dimensional (XD) UTE reconstruction. The initial motion fields were generated by registering images at each motion state to other motion states in motion-resolved images. A motion-state weighted motion-compensation (MostMoCo) reconstruction algorithm was proposed to reconstruct the dynamic UTE images. This technique, termed as MostMoCo-UTE, was compared with XD-UTE and iterative motion-compensation (iMoCo) on a porcine lung and 10 subjects.

RESULTS

MostMoCo reconstruction provides higher peak SNR (37.0 vs. 35.4 and 34.2) and structural similarity (0.964 vs. 0.931 and 0.947) compared to XD-UTE and iMoCo in the porcine lung experiment. Higher apparent SNR and contrast-to-noise ratio are achieved using MostMoCo in the human experiment. MostMoCo reconstruction better preserves the temporal variations of signal intensity of parenchyma compared to iMoCo, shows reduced random noise and improved sharpness of anatomical structures compared to XD-UTE. In the porcine lung experiment, the quantification of ventilation using MostMoCo images is more accurate than that using XD-UTE and iMoCo images.

CONCLUSION

The proposed MostMoCo-UTE provides improved quality of structural images and quantification of ventilation for free-breathing pulmonary MRI. It has the potential for the detection of structural and functional disorders of the lung in clinical settings.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T₁ ), the spin-spin relaxation time (T₂ ), and the spin-lattice relaxation in the rotating frame (T_1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T₁ -, T₂ -, or T_1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

Data-Consistent non-Cartesian deep subspace learning for efficient dynamic MR image reconstruction

Non-Cartesian sampling with subspace-constrained image reconstruction is a popular approach to dynamic MRI, but slow iterative reconstruction limits its clinical application. Data-consistent (DC) deep learning can accelerate reconstruction with good image quality, but has not been formulated for non-Cartesian subspace imaging. In this study, we propose a DC non-Cartesian deep subspace learning framework for fast, accurate dynamic MR image reconstruction. Four novel DC formulations are developed and evaluated: two gradient decent approaches, a directly solved approach, and a conjugate gradient approach. We applied a U-Net model with and without DC layers to reconstruct T1-weighted images for cardiac MR Multitasking (an advanced multidimensional imaging method), comparing our results to the iteratively reconstructed reference. Experimental results show that the proposed framework significantly improves reconstruction accuracy over the U-Net model without DC, while significantly accelerating the reconstruction over conventional iterative reconstruction.

Bridging cell-scale simulations and radiologic images to explain short-time intratumoral oxygen fluctuations

Radiologic images provide a way to monitor tumor development and its response to therapies in a longitudinal and minimally invasive fashion. However, they operate on a macroscopic scale (average value per voxel) and are not able to capture microscopic scale (cell-level) phenomena. Nevertheless, to examine the causes of frequent fast fluctuations in tissue oxygenation, models simulating individual cells’ behavior are needed. Here, we provide a link between the average data values recorded for radiologic images and the cellular and vascular architecture of the corresponding tissues. Using hybrid agent-based modeling, we generate a set of tissue morphologies capable of reproducing oxygenation levels observed in radiologic images. We then use these in silico tissues to investigate whether oxygen fluctuations can be explained by changes in vascular oxygen supply or by modulations in cellular oxygen absorption. Our studies show that intravascular changes in oxygen supply reproduce the observed fluctuations in tissue oxygenation in all considered regions of interest. However, larger-magnitude fluctuations cannot be recreated by modifications in cellular absorption of oxygen in a biologically feasible manner. Additionally, we develop a procedure to identify plausible tissue morphologies for a given temporal series of average data from radiology images. In future applications, this approach can be used to generate a set of tissues comparable with radiology images and to simulate tumor responses to various anti-cancer treatments at the tissue-scale level.

Low levels of oxygen, called hypoxia, are observable in many solid tumors. They are associated with more aggressive malignant cells that are resistant to chemo-, radio-, and immunotherapies. Recently developed imaging techniques provide a way to measure the magnitude of frequent short-term oxygen fluctuations, but they operate on a macro-scale voxel level. To examine the possible causes of rapid oxygen fluctuations at the cell level, we developed a hybrid agent-based mathematical model. We tested two different mechanisms that may be responsible for these cyclic effects on tissue oxygenation: temporal variations in vascular influx of oxygen and modulations in cellular oxygen absorption. Additionally, we developed a procedure to identify plausible tissue morphologies from data collected from radiological images. This can provide a bridge between the micro-scale simulations with individual cells and the longitudinal medical images containing average values. In future applications, this approach can be used to generate a set of tissues compatible with radiology images and to simulate tumor responses to various anticancer treatments at the cell-scale level.

Magnetization-prepared GRASP MRI for rapid 3D T1 mapping and fat/water-separated T1 mapping

PURPOSE

This study aimed to (i) develop Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI using a stack-of-stars trajectory for rapid free-breathing T1 mapping and (ii) extend MP-GRASP to multi-echo acquisition (MP-Dixon-GRASP) for fat/water-separated (water-specific) T1 mapping.

METHODS

An adiabatic non-selective 180° inversion-recovery pulse was added to a gradient-echo-based golden-angle stack-of-stars sequence for magnetization-prepared 3D single-echo or 3D multi-echo acquisition. In combination with subspace-based GRASP-Pro reconstruction, the sequence allows for standard T1 mapping (MP-GRASP) or fat/water-separated T1 mapping (MP-Dixon-GRASP), respectively. The accuracy of T1 mapping using MP-GRASP was evaluated in a phantom and volunteers (brain and liver) against clinically accepted reference methods. The repeatability of T1 estimation was also assessed in the phantom and volunteers. The performance of MP-Dixon-GRASP for water-specific T1 mapping was evaluated in a fat/water phantom and volunteers (brain and liver).

RESULTS

ROI-based mean T1 values are correlated between the references and MP-GRASP in the phantom (R2 = 1.0), brain (R2 = 0.96), and liver (R2 = 0.73). MP-GRASP achieved good repeatability of T1 estimation in the phantom (R2 = 1.0), brain (R2 = 0.99), and liver (R2 = 0.82). Water-specific T1 is different from in-phase and out-of-phase composite T1 (composite T1 when fat and water signal are mixed in phase or out of phase) both in the phantom and volunteers.

CONCLUSION

This work demonstrated the initial performance of MP-GRASP and MP-Dixon-GRASP MRI for rapid 3D T1 mapping and 3D fat/water-separated T1 mapping in the brain (without motion) and in the liver (during free breathing). With fat/water-separated T1 estimation, MP-Dixon-GRASP could be potentially useful for imaging patients with fatty-liver diseases.

Simultaneous evaluation of perfusion and morphology using GRASP MRI in hepatic fibrosis

OBJECTIVES

To determine if golden-angle radial sparse parallel (GRASP) dynamic contrast-enhanced (DCE)-MRI allows simultaneous evaluation of perfusion and morphology in liver fibrosis.

METHODS

Participants who were scheduled for liver biopsy or resection were enrolled (NCT02480972). Images were reconstructed at 12-s temporal resolution for morphologic assessment and at 3.3-s temporal resolution for quantitative evaluation. The image quality of the morphologic images was assessed on a four-point scale, and the Liver Imaging Reporting and Data System score was recorded for hepatic observations. Comparisons were made between quantitative parameters of DCE-MRI for the different fibrosis stages, and for hepatocellular carcinoma (HCCs) with different LR features.

RESULTS

DCE-MRI of 64 participants (male = 48) were analyzed. The overall image quality consistently stood at 3.5 ± 0.4 to 3.7 ± 0.4 throughout the exam. Portal blood flow significantly decreased in participants with F2-F3 (n = 18, 175 ± 110 mL/100 mL/min) and F4 (n = 12, 98 ± 47 mL/100 mL/min) compared with those in participants with F0-F1 (n = 34, 283 ± 178 mL/100 mL/min, p < 0.05 for all). In participants with F4, the arterial fraction and extracellular volume were significantly higher than those in participants with F0-F1 and F2-F3 (p < 0.05). Compared with HCCs showing non-LR-M features (n = 16), HCCs with LR-M (n = 5) had a significantly prolonged mean transit time and lower arterial blood flow (p < 0.05).

CONCLUSIONS

Liver MRI using GRASP obtains both sufficient spatial resolution for confident diagnosis and high temporal resolution for pharmacokinetic modeling. Significant differences were found between the MRI-derived portal blood flow at different hepatic fibrosis stages.

KEY POINTS

A single MRI examination is able to provide both images with sufficient spatial resolution for anatomic evaluation and those with high temporal resolution for pharmacokinetic modeling. Portal blood flow was significantly lower in clinically significant hepatic fibrosis and mean transit time and extracellular volume increased in cirrhosis, compared with those in no or mild hepatic fibrosis. HCCs with different LR features showed different quantitative parameters of DCE-MRI: longer mean transit time and lower arterial flow were observed in HCCs with LR-M features.

Magnetic resonance parameter mapping using model-guided self-supervised deep learning

PURPOSE

To develop a model-guided self-supervised deep learning MRI reconstruction framework called reference-free latent map extraction (RELAX) for rapid quantitative MR parameter mapping.

METHODS

Two physical models are incorporated for network training in RELAX, including the inherent MR imaging model and a quantitative model that is used to fit parameters in quantitative MRI. By enforcing these physical model constraints, RELAX eliminates the need for full sampled reference data sets that are required in standard supervised learning. Meanwhile, RELAX also enables direct reconstruction of corresponding MR parameter maps from undersampled k-space. Generic sparsity constraints used in conventional iterative reconstruction, such as the total variation constraint, can be additionally included in the RELAX framework to improve reconstruction quality. The performance of RELAX was tested for accelerated T1 and T2 mapping in both simulated and actually acquired MRI data sets and was compared with supervised learning and conventional constrained reconstruction for suppressing noise and/or undersampling-induced artifacts.

RESULTS

In the simulated data sets, RELAX generated good T1 /T2 maps in the presence of noise and/or undersampling artifacts, comparable to artifact/noise-free ground truth. The inclusion of a spatial total variation constraint helps improve image quality. For the in vivo T1 /T2 mapping data sets, RELAX achieved superior reconstruction quality compared with conventional iterative reconstruction, and similar reconstruction performance to supervised deep learning reconstruction.

CONCLUSION

This work has demonstrated the initial feasibility of rapid quantitative MR parameter mapping based on self-supervised deep learning. The RELAX framework may also be further extended to other quantitative MRI applications by incorporating corresponding quantitative imaging models.

On the use of low-dimensional temporal subspace constraints to reduce reconstruction time and improve image quality of accelerated 4D-MRI

BACKGROUND AND PURPOSE

The purpose of this work is to investigate the use of low-dimensional temporal subspace constraints for 4D-MRI reconstruction from accelerated data in the context of MR-guided online adaptive radiation therapy (MRgOART).

MATERIALS AND METHODS

Subspace basis functions are derived directly from the accelerated golden angle radial stack-of-stars 4D-MRI data. The reconstruction times, image quality, and motion estimates are investigated as a function of the number of subspace coefficients and compared with a conventional frame-by-frame reconstruction. These experiments were performed in five patients with four 4D-MRI scans per patient on a 1.5T MR-Linac.

RESULTS

If two or three subspace coefficients are used, the iterative reconstruction time is reduced by 32% and 18%, respectively, compared to conventional parallel imaging with compressed sensing reconstructions. No significant difference was found between motion estimates made with the subspace-constrained reconstructions (p > 0.08). Qualitative improvements in image quality included reduction in apparent noise and reductions in streaking artifacts from the radial k-space coverage.

CONCLUSION

Incorporating subspace constraints for accelerated 4D-MRI reconstruction reduces noise and residual undersampling artifacts in the images while reducing computation time, making it a strong candidate for use in clinical MRgOART workflows.

Rapid golden-angle diffusion-weighted propeller MRI for simultaneous assessment of ADC and IVIM

Purpose:

Golden-angle single-shot PROPLLER (GA-SS-PROP) is proposed to accelerate the PROPELLER acquisition for distortion-free diffusion-weighted (DW) imaging. Acceleration is achieved by acquiring one-shot per b-value and several b-values can be acquired along a diffusion direction, where the DW signal follows a bi-exponential decay (i.e. IVIM). Sparse reconstruction is used to reconstruct full resolution DW images. Consequently, apparent diffusion coefficient (ADC) map and IVIM maps (i.e., perfusion fraction (f) and the perfusion-free diffusion coefficient (D)) are obtained simultaneously. The performance of GA-SS-PROP was demonstrated with simulation and human experiments.

Methods:

A realistic numerical phantom of high-quality diffusion images of the brain was developed. The error of the reconstructed DW images and quantitative maps were compared to the ground truth. The pulse sequence was developed to acquire human brain data. For comparison, fully sampled PROPELLER and conventional single-shot echo planar imaging (SS-EPI) acquisitions were performed.

Results:

GA-SS-PROP was 5 times faster than conventional PROPELLER acquisition with comparable image quality. The simulation demonstrated that sparse reconstruction is effective in restoring contrast and resolution. The human experiments demonstrated that GA-SS-PROP achieved superior image fidelity compared to SS-EPI for the same acquisition time and same in-plane resolution (1 × 1 mm²).

Conclusion:

GA-SS-PROP offers fast, high-resolution and distortion-free DW images. The generated quantitative maps (f, D and ADC) can provide valuable information on tissue perfusion and diffusion properties simultaneously, which are desirable in many applications, especially in oncology. As a turbo spin-echo based technique, it can be applied in most challenging regions where SS-EPI is problematic.

Extreme MRI: Large-scale volumetric dynamic imaging from continuous non-gated acquisitions

PURPOSE

To develop a framework to reconstruct large-scale volumetric dynamic MRI from rapid continuous and non-gated acquisitions, with applications to pulmonary and dynamic contrast-enhanced (DCE) imaging.

THEORY AND METHODS

The problem considered here requires recovering 100 gigabytes of dynamic volumetric image data from a few gigabytes of k-space data, acquired continuously over several minutes. This reconstruction is vastly under-determined, heavily stressing computing resources as well as memory management and storage. To overcome these challenges, we leverage intrinsic three-dimensional (3D) trajectories, such as 3D radial and 3D cones, with ordering that incoherently cover time and k-space over the entire acquisition. We then propose two innovations: (a) A compressed representation using multiscale low-rank matrix factorization that constrains the reconstruction problem, and reduces its memory footprint. (b) Stochastic optimization to reduce computation, improve memory locality, and minimize communications between threads and processors. We demonstrate the feasibility of the proposed method on DCE imaging acquired with a golden-angle ordered 3D cones trajectory and pulmonary imaging acquired with a bit-reversed ordered 3D radial trajectory. We compare it with "soft-gated" dynamic reconstruction for DCE and respiratory-resolved reconstruction for pulmonary imaging.

RESULTS

The proposed technique shows transient dynamics that are not seen in gating-based methods. When applied to datasets with irregular, or non-repetitive motions, the proposed method displays sharper image features.

CONCLUSIONS

We demonstrated a method that can reconstruct massive 3D dynamic image series in the extreme undersampling and extreme computation setting.

MRSIGMA: Magnetic Resonance SIGnature MAtching for real-time volumetric imaging

PURPOSE

To propose a real-time 3D MRI technique called MR SIGnature MAtching (MRSIGMA) for high-resolution volumetric imaging and motion tracking with very low imaging latency.

METHODS

MRSIGMA consists of two steps: (1) offline learning of a database of possible 3D motion states and corresponding motion signature ranges and (2) online matching of new motion signatures acquired in real time with prelearned motion states. Specifically, the offline learning step (non-real-time) reconstructs motion-resolved 4D images representing different motion states and assigns a unique motion range to each state. The online matching step (real-time) acquires motion signatures only and selects one of the prelearned 3D motion states for each newly acquired signature, which generates 3D images efficiently in real time. The MRSIGMA technique was evaluated on 15 golden-angle stack-of-stars liver data sets, and the performance of respiratory motion tracking with the online-generated real-time 3D MRI was compared with the corresponding 2D projections acquired in real time.

RESULTS

The total latency of generating each 3D image during online matching was about 300 ms, including acquisition of the motion signature data (~138 ms) and corresponding matching process (~150 ms). Linear correlation assessment suggested excellent correlation (R²  = 0.948) between motion displacement measured from the online-generated real-time 3D images and the 2D real-time projections.

CONCLUSION

This proof-of-concept study demonstrates the feasibility of MRSIGMA for high-resolution real-time volumetric imaging, which shifts the acquisition and reconstruction burden to an offline learning step and leaves fast online matching for online imaging with very low imaging latency. The MRSIGMA technique can potentially be used for real-time motion tracking in MRI-guided radiation therapy.

Improving dynamic contrast-enhanced MRI of the lung using motion-weighted sparse reconstruction: Initial experiences in patients

PURPOSE

The purpose of this study was to evaluate the performance of motion-weighted Golden-angle RAdial Sparse Parallel MRI (motion-weighted GRASP) for free-breathing dynamic contrast-enhanced MRI (DCE-MRI) of the lung.

METHODS

Motion-weighted GRASP incorporates a soft-gating motion compensation algorithm into standard GRASP reconstruction, so that motion-corrupted motion k-space (e.g., k-space acquired in inspiratory phases) contributes less to the final reconstructed images. Lung MR data from 20 patients (mean age = 57.9 ± 13.5) with known pulmonary lesions were retrospectively collected for this study. Each subject underwent a free-breathing DCE-MR scan using a fat-statured T₁-weighted stack-of-stars golden-angle radial sequence and a post-contrast breath-hold MR scan using a Cartesian volumetric-interpolated imaging sequence (BH-VIBE). Each radial dataset was reconstructed using GRASP without motion compensation and motion-weighted GRASP. All MR images were visually evaluated by two experienced radiologists blinded to reconstruction and acquisition schemes independently. In addition, the influence of motion-weighted reconstruction on dynamic contrast-enhancement patterns was also investigated.

RESULTS

For image quality assessment, motion-weighted GRASP received significantly higher visual scores than GRASP (P < 0.05) for overall image quality (3.68 vs. 3.39), lesion conspicuity (3.54 vs. 3.18) and overall artifact level (3.53 vs. 3.15). There was no significant difference (P > 0.05) between the breath-hold BH-VIBE and motion-weighted GRASP images. For assessment of temporal fidelity, motion-weighted GRASP maintained a good agreement with respect to GRASP.

CONCLUSION

Motion-weighted GRASP achieved better reconstruction performance in free-breathing DCE-MRI of the lung compared to standard GRASP, and it may enable improved assessment of pulmonary lesions.

Six-dimensional quantitative DCE MR Multitasking of the entire abdomen: Method and application to pancreatic ductal adenocarcinoma

PURPOSE

To develop a quantitative DCE MRI technique enabling entire-abdomen coverage, free-breathing acquisition, 1-second temporal resolution, and T₁ -based quantification of contrast agent concentration and kinetic modeling for the characterization of pancreatic ductal adenocarcinoma (PDAC).

METHODS

Segmented FLASH readouts following saturation-recovery preparation with randomized 3D Cartesian undersampling was used for incoherent data acquisition. MR Multitasking was used to reconstruct 6-dimensional images with 3 spatial dimensions, 1 T₁ recovery dimension for dynamic T₁ quantification, 1 respiratory dimension to resolve respiratory motion, and 1 DCE time dimension to capture the contrast kinetics. Sixteen healthy subjects and 14 patients with pathologically confirmed PDAC were recruited for the in vivo studies, and kinetic parameters v_p , K^trans , v_e , and K_ep were evaluated for each subject. Intersession repeatability of Multitasking DCE was assessed in 8 repeat healthy subjects. One-way unbalanced analysis of variance was performed between control and patient groups.

RESULTS

In vivo studies demonstrated that v_p , K^trans , and K_ep of PDAC were significantly lower compared with nontumoral regions in the patient group (P = .002, .003, .004, respectively) and normal pancreas in the control group (P = .011, <.001, <.001, respectively), while v_e was significantly higher than nontumoral regions (P < .001) and healthy pancreas (P < .001). The kinetic parameters showed good in vivo repeatability (interclass correlation coefficient: v_p , 0.95; K^trans , 0.98; v_e , 0.96; K_ep , 0.99).

CONCLUSION

The proposed Multitasking DCE is promising for the quantification of vascular properties of PDAC. Quantitative DCE parameters were repeatable in vivo and showed significant differences between normal pancreas and both tumor and nontumoral regions in patients with PDAC.