An optimal radial profile order based on the Golden Ratio for time-resolved MRI

Stability of conventional and machine learning-based tumor auto-segmentation techniques using undersampled dynamic radial bSSFP acquisitions on a 0.35 T hybrid MR-linac system

PURPOSE

Hybrid MRI-linear accelerator systems (MR-linacs) allow for the incorporation of MR images with high soft-tissue contrast into the radiation therapy procedure prior to, during, or post irradiation. This allows not only for the optimization of the treatment planning, but also for real-time monitoring of the tumor position using cine MRI, from which intrafractional motion can be compensated. Fast imaging and accurate tumor tracking are crucial for effective compensation. This study investigates the application of cine MRI with a radial acquisition scheme on a low-field MR-linac to accelerate the acquisition rate and evaluates the effect on tracking accuracy.

METHODS

An MR sequence using tiny golden-angle radial k-space sampling was developed and applied to cine imaging on patients with liver tumors on a 0.35 T MR-linac. Tumor tracking was assessed for accuracy and stability from the cine images with increasing k-space undersampling factors. Tracking was achieved using two different auto-segmentation algorithms: a deformable image registration B-spline similar to that implemented on the MR-linac and a convolutional neural network approach known as U-Net.

RESULTS

Radial imaging allows for increased temporal resolution with reliable tumor tracking, although tracking robustness decreases as temporal resolution increases. Additional acquisition-based artifacts can be avoided by reducing the angle increment using tiny golden-angles. The U-net algorithm was found to have superior auto-segmentation metrics compared to B-spline. U-net was able to track two well-defined tumors, imaged with just 30 spokes per image (10.6 frames per second), with an average Dice coefficient ≥ 83%, Hausdorff distance ≤ 1.4 pixel, and mean contour distance ≤ 0.5 pixel.

CONCLUSIONS

Radial acquisitions are commonplace in dynamic imaging; however, in MR-guided radiotherapy, robust tumor tracking is also required. This study demonstrates the in vivo feasibility of tumor tracking from radially acquired images on a low-field MR-linac. Radial imaging allows for decreased image acquisition times while maintaining robust tracking. The U-net algorithm can track a tumor with higher accuracy in images with undersampling artifacts than a conventional deformable B-spline algorithm and is a promising tool for tracking in MR-guided radiation therapy.

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.

Closed-loop control of k-space sampling via physiologic feedback for cine MRI

BACKGROUND

Segmented cine cardiac MRI combines data from multiple heartbeats to achieve high spatiotemporal resolution cardiac images, yet predefined k-space segmentation trajectories can lead to suboptimal k-space sampling. In this work, we developed and evaluated an autonomous and closed-loop control system for radial k-space sampling (ARKS) to increase sampling uniformity.

METHODS

The closed-loop system autonomously selects radial k-space sampling trajectory during live segmented cine MRI and attempts to optimize angular sampling uniformity by selecting views in regions of k-space that were not previously well-sampled. Sampling uniformity and the ability to detect cardiac phase in vivo was assessed using ECG data acquired from 10 normal subjects in an MRI scanner. The approach was then implemented with a fast gradient echo sequence on a whole-body clinical MRI scanner and imaging was performed in 4 healthy volunteers. The closed-loop k-space trajectory was compared to random, uniformly distributed and golden angle view trajectories via measurement of k-space uniformity and the point spread function. Lastly, an arrhythmic dataset was used to evaluate a potential application of the approach.

RESULTS

The autonomous trajectory increased k-space sampling uniformity by 15±7%, main lobe point spread function (PSF) signal intensity by 6±4%, and reduced ringing relative to golden angle sampling. When implemented, the autonomous pulse sequence prescribed radial view angles faster than the scan TR (0.98 ± 0.01 ms, maximum = 1.38 ms) and increased k-space sampling mean uniformity by 10±11%, decreased uniformity variability by 44±12%, and increased PSF signal ratio by 6±6% relative to golden angle sampling.

CONCLUSION

The closed-loop approach enables near-uniform radial sampling in a segmented acquisition approach which was higher than predetermined golden-angle radial sampling. This can be utilized to increase the sampling or decrease the temporal footprint of an acquisition and the closed-loop framework has the potential to be applied to patients with complex heart rhythms.

Evaluation of Golden-Angle-Sampled Dynamic Contrast-Enhanced MRI Reconstruction Using Objective Image Quality Measures: A Simulated Phantom Study

We aim to extend the use of image quality metrics (IQMs) from static magnetic resonance imaging (MRI) applications to dynamic MRI studies. We assessed the use of 2 IQMs, the root mean square error and structural similarity index, in evaluating the reconstruction of quantitative dynamic contrast-enhanced (DCE) MRI data acquired using golden-angle sampling and compressed sensing (CS). To address the difficulty of obtaining ground-truth knowledge of parameters describing dynamics in real patient data, we developed a Matlab simulation framework to assess quantitative CS-DCE-MRI. We began by validating the response of each IQM to the CS-MRI reconstruction process using static data and the performance of our simulation framework with simple dynamic data. We then extended the simulations to the more realistic extended Tofts model. When assessing the Tofts model, we tested 4 different methods of selecting a reference image for the IQMs. Results from the retrospective static CS-MRI reconstructions showed that each IQM is responsive to the CS-MRI reconstruction process. Simulations of a simple contrast evolution model validated the performance of our framework. Despite the complexity of the Tofts model, both IQM scores correlated well with the recovery accuracy of a central model parameter for all reference cases studied. This finding may form the basis of algorithms for automated selection of image reconstruction aspects, such as temporal resolution, in golden-angle-sampled CS-DCE-MRI. These further suggest that objective measures of image quality may find use in general dynamic MRI applications.

Adaptive sparsity level and dictionary size estimation for image reconstruction in accelerated 2D radial cine MRI

PURPOSE

In the past, dictionary learning (DL) and sparse coding (SC) have been proposed for the regularization of image reconstruction problems. The regularization is given by a sparse approximation of all image patches using a learned dictionary, that is, an overcomplete set of basis functions learned from data. Despite its competitiveness, DL and SC require the tuning of two essential hyperparameters: the sparsity level S - the number of basis functions of the dictionary, called atoms, which are used to approximate each patch, and K - the overall number of such atoms in the dictionary. These two hyperparameters usually have to be chosen a priori and are determined by repetitive and computationally expensive experiments. Furthermore, the final reported values vary depending on the specific situation. As a result, the clinical application of the method is limited, as standardized reconstruction protocols have to be used.

METHODS

In this work, we use adaptive DL and propose a novel adaptive sparse coding algorithm for two-dimensional (2D) radial cine MR image reconstruction. Using adaptive DL and adaptive SC, the optimal dictionary size K as well as the optimal sparsity level S are chosen dependent on the considered data.

RESULTS

Our three main results are the following: First, adaptive DL and adaptive SC deliver results which are comparable or better than the most widely used nonadaptive version of DL and SC. Second, the time needed for the regularization is accelerated due to the fact that the sparsity level S is never overestimated. Finally, the a priori choice of S and K is no longer needed but is optimally chosen dependent on the data under consideration.

CONCLUSIONS

Adaptive DL and adaptive SC can highly facilitate the application of DL- and SC-based regularization methods. While in this work we focused on 2D radial cine MR image reconstruction, we expect the method to be applicable to different imaging modalities as well.

Nonrigid 3D motion estimation at high temporal resolution from prospectively undersampled k-space data using low-rank MR-MOTUS

Purpose

With the recent introduction of the MR‐LINAC, an MR‐scanner combined with a radiotherapy LINAC, MR‐based motion estimation has become of increasing interest to (retrospectively) characterize tumor and organs‐at‐risk motion during radiotherapy. To this extent, we introduce low‐rank MR‐MOTUS, a framework to retrospectively reconstruct time‐resolved nonrigid 3D+t motion fields from a single low‐resolution reference image and prospectively undersampled k‐space data acquired during motion.

Theory

Low‐rank MR‐MOTUS exploits spatiotemporal correlations in internal body motion with a low‐rank motion model, and inverts a signal model that relates motion fields directly to a reference image and k‐space data. The low‐rank model reduces the degrees‐of‐freedom, memory consumption, and reconstruction times by assuming a factorization of space‐time motion fields in spatial and temporal components.

Methods

Low‐rank MR‐MOTUS was employed to estimate motion in 2D/3D abdominothoracic scans and 3D head scans. Data were acquired using golden‐ratio radial readouts. Reconstructed 2D and 3D respiratory motion fields were, respectively, validated against time‐resolved and respiratory‐resolved image reconstructions, and the head motion against static image reconstructions from fully sampled data acquired right before and right after the motion.

Results

Results show that 2D+t respiratory motion can be estimated retrospectively at 40.8 motion fields per second, 3D+t respiratory motion at 7.6 motion fields per second and 3D+t head‐neck motion at 9.3 motion fields per second. The validations show good consistency with image reconstructions.

Conclusions

The proposed framework can estimate time‐resolved nonrigid 3D motion fields, which allows to characterize drifts and intra and inter‐cycle patterns in breathing motion during radiotherapy, and could form the basis for real‐time MR‐guided radiotherapy.

Quantitative T2 mapping using accelerated 3D stack-of-spiral gradient echo readout

PURPOSE

To develop a rapid T2 mapping protocol using optimized spiral acquisition, accelerated reconstruction, and model fitting.

MATERIALS AND METHODS

A T2-prepared stack-of-spiral gradient echo (GRE) pulse sequence was applied. A model-based approach joined with compressed sensing was compared with the two methods applied separately for accelerated reconstruction and T2 mapping. A 2-parameter-weighted fitting method was compared with 2- or 3-parameter models for accurate T2 estimation under the influences of noise and B1 inhomogeneity. The performance was evaluated using both digital phantoms and healthy volunteers. Mitigating partial voluming with cerebrospinal fluid (CSF) was also tested.

RESULTS

Simulations demonstrates that the 2-parameter-weighted fitting approach was robust to a large range of B1 scales and SNR levels. With an in-plane acceleration factor of 5, the model-based compressed sensing-incorporated method yielded around 8% normalized errors compared to references. The T2 estimation with and without CSF nulling was consistent with literature values.

CONCLUSION

This work demonstrated the feasibility of a T2 quantification technique with 3D high-resolution and whole-brain coverage in 2-3 min. The proposed iterative reconstruction method, which utilized the model consistency, data consistency and spatial sparsity jointly, provided reasonable T2 estimation. The technique also allowed mitigation of CSF partial volume effect.

MR fingerprinting for rapid simultaneous T1 , T2 , and T1ρ relaxation mapping of the human articular cartilage at 3T

PURPOSE

To implement a novel technique for simultaneous, quantitative multiparametric mapping of the knee articular cartilage.

METHODS

A novel MR fingerprinting pulse sequence is proposed and implemented for simultaneous measurements of proton density, T1 , T2, and T1ρ relaxation times at 3T. The repeatability and reproducibility of the proposed technique were assessed in model phantoms. Institutional review board-approved MR fingerprinting imaging sequence was performed on healthy volunteers and patients with mild knee osteoarthritis. The Wilcoxon test was used to compare healthy controls and patients. The intra- and intersubject repeatability were assessed with coefficient of variation and the RMS coefficient of variation, respectively RESULTS: The Bland-Altman plots demonstrated an average difference of 4.67 ms, -0.09 ms, and 0.05 ms between 2 scans in the same scanner; and 9.68 ms, 0.29 ms, and -0.72 ms between the scans acquired on 2 different scanners for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with RMS coefficient of variation less than 3%, 6%, and 5% for T1 , T2 , and T1ρ , respectively. The Wilcoxon test showed a significant difference between control and mild osteoarthritis patients for T1 (P = .04), T2 (P = .01), and T1ρ (P = .02) relaxation time in medial tibial cartilage, as well as for T2 relaxation time (P = .02) in medial femoral cartilage.

CONCLUSION

The proposed MRF sequence is fast and can simultaneously measure the T1 , T2 , T1ρ , and B 1 + maps in a single scan. It is able to discriminate between mild osteoarthritis patients and healthy volunteers.

Deep learning reconstruction for cardiac magnetic resonance fingerprinting T1 and T2 mapping

PURPOSE

To develop a deep learning method for rapidly reconstructing T₁ and T₂ maps from undersampled electrocardiogram (ECG) triggered cardiac magnetic resonance fingerprinting (cMRF) images.

METHODS

A neural network was developed that outputs T₁ and T₂ values when given a measured cMRF signal time course and cardiac RR interval times recorded by an ECG. Over 8 million cMRF signals, corresponding to 4000 random cardiac rhythms, were simulated for training. The training signals were corrupted by simulated k-space undersampling artifacts and random phase shifts to promote robust learning. The deep learning reconstruction was evaluated in Monte Carlo simulations for a variety of cardiac rhythms and compared with dictionary-based pattern matching in 58 healthy subjects at 1.5T.

RESULTS

In simulations, the normalized root-mean-square error (nRMSE) for T₁ was below 1% in myocardium, blood, and liver for all tested heart rates. For T₂ , the nRMSE was below 4% for myocardium and liver and below 6% for blood for all heart rates. The difference in the mean myocardial T₁ or T₂ observed in vivo between dictionary matching and deep learning was 3.6 ms for T₁ and -0.2 ms for T₂ . Whereas dictionary generation and pattern matching required more than 4 min per slice, the deep learning reconstruction only required 336 ms.

CONCLUSION

A neural network is introduced for reconstructing cMRF T₁ and T₂ maps directly from undersampled spiral images in under 400 ms and is robust to arbitrary cardiac rhythms, which paves the way for rapid online display of cMRF maps.

Real-time exercise stress cardiac MRI with Fourier-series reconstruction from golden-angle radial data

Magnetic resonance imaging (MRI) can measure cardiac response to exercise stress for evaluating and managing heart patients in the practice of clinical cardiology. However, exercise stress cardiac MRI have been clinically limited by the ability of available MRI techniques to quantitatively measure fast and unstable cardiac dynamics during exercise. The presented work is to develop a new real-time MRI technique for improved quantitative performance of exercise stress cardiac MRI. This technique seeks to represent real-time cardiac images as a sparse Fourier-series along the time. With golden-angle radial acquisition, parallel imaging and compressed sensing can be integrated into a linear system of equations for resolving Fourier coefficients that are in turn used to generate real-time cardiac images from the Fourier-series representation. Fourier-series reconstruction from golden-angle radial data can effectively address data insufficiency due to MRI speed limitation, providing a real-time approach to exercise stress cardiac MRI. To demonstrate the feasibility, an exercise stress cardiac MRI experiment was run to investigate biventricular response to in-scanner biking exercise in a cohort of sixteen healthy volunteers. It was found that Fourier-series reconstruction from golden-angle radial data effectively detected exercise-induced increase in stroke volume and ejection fraction in a healthy heart. The presented work will improve the applications of exercise stress cardiac MRI in the practice of clinical cardiology.

Binding site exchange kinetics revealed through efficient spin-spin dephasing of hyperpolarized 129Xe

Spin exchange between different chemical environments is an important observable for characterizing chemical exchange kinetics in various contexts, including protein folding, chelation chemistry, and host-guest interactions. Such spins experience effective spin-spin relaxation rate, R _2,eff, that typically shows a dispersive behavior which requires detailed analysis. Here, we describe a class of highly simplified R _2,eff behavior by relying on hyperpolarized ¹²⁹Xe as a freely exchanging ligand reporter. It provides large chemical shift separations that yield reduced expressions of both the Swift-Connick and the Carver-Richards treatment of exchange-induced relaxation. Despite observing a diamagnetic system, R _2,eff is dominated by large Larmor frequency jumps and thus allows detection of otherwise inaccessible analyte concentrations with a single spin echo train (only 0.01% of the overall hyperpolarized spins need to be transiently bound to the molecule). The two Xe hosts cryptophane-A monoacid (CrA-ma) and cucurbit[6]uril (CB6) represent two exemplary families of container molecules (the latter one also serving as drug delivery vehicles) that act as highly efficient phase shifters for which we observed unprecedented exchange-induced relaxivity r ₂ (up to 866 s^-1 mM^-1). By including methods of spatial encoding, multiple data points can be collected simultaneously to isolate the exchange contribution and determine the effective exchange rate in partially occupied binding sites with a single delivery of hyperpolarized nuclei. The relaxivity is directly related to the guest turnover in these systems and temperature-dependent measurements yield an activation energy of E _A = 41 kJ mol^-1 for Xe@CrA-ma from simple relaxometry analysis. The concept is transferable to many applications where Xe is known to exhibit large chemical shifts.

Non-contrast coronary magnetic resonance angiography: current frontiers and future horizons

Coronary magnetic resonance angiography (coronary MRA) is advantageous in its ability to assess coronary artery morphology and function without ionizing radiation or contrast media. However, technical limitations including reduced spatial resolution, long acquisition times, and low signal-to-noise ratios prevent it from clinical routine utilization. Nonetheless, each of these limitations can be specifically addressed by a combination of novel technologies including super-resolution imaging, compressed sensing, and deep-learning reconstruction. In this paper, we first review the current clinical use and motivations for non-contrast coronary MRA, discuss currently available coronary MRA techniques, and highlight current technical developments that hold unique potential to optimize coronary MRA image acquisition and post-processing. In the final section, we examine the various research-based coronary MRA methods and metrics that can be leveraged to assess coronary stenosis severity, physiological function, and atherosclerotic plaque characterization. We specifically discuss how such technologies may contribute to the clinical translation of coronary MRA into a robust modality for routine clinical use.

Fully self-gated whole-heart 4D flow imaging from a 5-minute scan

PURPOSE

To develop and validate an acquisition and processing technique that enables fully self-gated 4D flow imaging with whole-heart coverage in a fixed 5-minute scan.

THEORY AND METHODS

The data are acquired continuously using Cartesian sampling and sorted into respiratory and cardiac bins using the self-gating signal. The reconstruction is performed using a recently proposed Bayesian method called ReVEAL4D. ReVEAL4D is validated using data from 8 healthy volunteers and 2 patients and compared with compressed sensing technique, L1-SENSE.

RESULTS

Healthy subjects-Compared with 2D phase-contrast MRI (2D-PC), flow quantification from ReVEAL4D shows no significant bias. In contrast, the peak velocity and peak flow rate for L1-SENSE are significantly underestimated. Compared with traditional parallel MRI-based 4D flow imaging, ReVEAL4D demonstrates small but significant biases in net flow and peak flow rate, with no significant bias in peak velocity. All 3 indices are significantly and more markedly underestimated by L1-SENSE. Patients-Flow quantification from ReVEAL4D agrees well with the 2D-PC reference. In contrast, L1-SENSE markedly underestimated peak velocity.

CONCLUSIONS

The combination of highly accelerated 5-minute Cartesian acquisition, self-gating, and ReVEAL4D enables whole-heart 4D flow imaging with accurate flow quantification.

Quantitative Skeletal Muscle Imaging Using 3D MR Fingerprinting With Water and Fat Separation

BACKGROUND

Quantitative muscle MRI is a robust tool to monitor intramuscular fatty replacement and disease activity in patients with neuromuscular disorders (NMDs).

PURPOSE

To implement a 3D sequence for quantifying simultaneously fat fraction (FF) and water T₁ (T_1,H2O ) in the skeletal muscle, evaluate regular undersampling in the partition-encoding direction, and compare it to a recently proposed 2D MR fingerprinting sequence with water and fat separation (MRF T₁ -FF).

STUDY TYPE

Prospective.

PHANTOM/SUBJECTS

Seventeen-vial phantom at different FF and T_1,H2O , 11 healthy volunteers, and 6 subjects with different NMDs.

FIELD STRENGTH/SEQUENCE

3T/3D MRF T₁ -FF, 2D MRF T₁ -FF, STEAM MRS ASSESSMENT: FF and T_1,H2O measured with the 2D and 3D sequences were compared in the phantom and in vivo at different undersampling factors (US). Data were acquired in healthy subjects before and after plantar dorsiflexions and at rest in patients.

STATISTICAL TESTS

Linear correlations, Bland-Altman analysis, two-way repeated measures analysis of variance (ANOVA), Student's t-test.

RESULTS

Up to a US factor of 3, the undersampled acquisitions were in good agreement with the fully sampled sequence (R²  ≥ 0.98, T_1,H2O bias ≤10 msec, FF bias ≤4 × 10^-4 ) both in phantom and in vivo. The 2D and 3D MRF T₁ -FF sequences provided comparable T_1,H2O and FF values (R²  ≥ 0.95, absolute T_1,H2O bias ≤35 msec, and absolute FF bias ≤0.003). The plantar dorsiflexion induced a significant increase of T_1,H2O in the tibialis anterior and extensor digitorum (relative increase of +10.8 ± 1.7% and + 7.7 ± 1.4%, respectively, P < 0.05), that was accompanied by a significant reduction of FF in the tibialis anterior (relative decrease of -16.3 ± 4.0%, P < 0.05). Some subjects with NMDs presented increased and heterogeneous T_1,H2O and FF values throughout the leg.

DATA CONCLUSION

Quantitative 3D T_1,H2O and FF maps covering the entire leg were obtained within acquisition times compatible with clinical research (4 minutes 20 seconds) and a 1 × 1 × 5 mm³ spatial resolution.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Achieving high spatial and temporal resolution with perfusion MRI in the head and neck region using golden-angle radial sampling

Objectives

Conventional perfusion-weighted MRI sequences often provide poor spatial or temporal resolution. We aimed to overcome this problem in head and neck protocols using a golden-angle radial sparse parallel (GRASP) sequence.

Methods

We prospectively included 58 patients for examination on a 3.0-T MRI using a study protocol. GRASP (A) was applied to a volumetric interpolated breath-hold examination (VIBE) with 135 reconstructed pictures and high temporal (2.5 s) and spatial resolution (0.94 × 0.94 × 3.00 mm). Additional sequences of matching temporal resolution (B: 2.5 s, 1.88 × 1.88 × 3.00 mm), with a compromise between temporal and spatial resolution (C: 7.0 s, 1.30 × 1.30 × 3.00 mm) and with matching spatial resolution (D: 145 s, 0.94 × 0.94 × 3.00 mm), were subsequently without GRASP. Instant inline-image reconstructions (E) provided one additional series of averaged contrast information throughout the entire acquisition duration of A. Overall diagnostic image quality, edge sharpness and contrast of soft tissues, vessels and lesions were subjectively rated using 5-point Likert scales. Objective image quality was measured as contrast-to-noise ratio in D and E.

Results

Overall, the anatomic and pathologic image quality was substantially better with the GRASP sequence for the temporally (A/B/C, all p < 0.001) and spatially resolved comparisons (D/E, all p < 0.002 except lesion edge sharpness with p = 0.291). Image artefacts were also less likely to occur with GRASP. Differences in motion, aliasing and truncation were mainly significant, but pulsation and fat suppression were comparable. In addition, the contrast-to-noise ratio of E was significantly better than that of D (p_D-E < 0.001).

Conclusions

High temporal and spatial resolution can be obtained synchronously using a GRASP-VIBE technique for perfusion evaluation in head and neck MRI.

Key Points

• Golden-angle radial sparse parallel (GRASP) sampling allows for temporally resolved dynamic acquisitions with a very high image quality.

• Very low-contrast structures in the head and neck region can benefit from using the GRASP sequence.

• Inline-image reconstruction of dynamic and static series from one single acquisition can replace the conventional combination of two acquisitions, thereby saving examination time.

Electronic supplementary material

The online version of this article (10.1007/s00330-020-07263-0) contains supplementary material, which is available to authorized users.

AZTEK: Adaptive zero TE k-space trajectories

PURPOSE

Because of short signal lifetimes and respiratory motion, 3D lung MRI is still challenging today. Zero-TE (ZTE) pulse sequences offer promising solutions as they overcome the issue of short T 2 ∗ . Nevertheless, as they rely on continuous readout gradients, the trajectories they follow in k-space are not adapted to retrospective gating and inferred motion correction.

THEORY AND METHODS

We propose AZTEK (adaptive ZTE k-space trajectories), a set of 3D radial trajectories featuring three tuning parameters, to adapt the acquisition to any moving organ while keeping seamless transitions between consecutive spokes. Standard ZTE and AZTEK trajectories were compared for static and moving phantom acquisitions as well as for human thoracic imaging performed on 3 volunteers (1 healthy and 2 patients with lung cancer).

RESULTS

For the static phantom, we observe comparable image qualities with standard and AZTEK trajectories. For the moving phantom, spatially coherent undersampling artifacts observed on gated images with the standard trajectory are alleviated with AZTEK. The same improvement in image quality is obtained in human, so details are more delineated in the lung with the use of the adaptive trajectory.

CONCLUSION

The AZTEK technique opens the possibility for 3D dynamic ZTE lung imaging with retrospective gating. It enables us to uniformly sample the k-space for any arbitrary respiratory motion gate, while preserving static image quality, improving dynamic image quality and guaranteeing continuous readout gradient transitions between spokes, which makes it appropriate to ZTE.

Silent 3D MR sequence for quantitative and multicontrast T1 and proton density imaging

This study aims to develop a silent, fast and 3D method for T1 and proton density (PD) mapping, while generating time series of T1-weighted (T1w) images with bias-field correction. Undersampled T1w images at different effective inversion times (TIs) were acquired using the inversion recovery prepared RUFIS sequence with an interleaved k-space trajectory. Unaliased images were reconstructed by constraining the signal evolution to a temporal subspace which was learned from the signal model. Parameter maps were obtained by fitting the data to the signal model, and bias-field correction was conducted on T1w images. Accuracy and repeatability of the method was accessed in repeated experiments with phantom and volunteers. For the phantom study, T1 values obtained by the proposed method were highly consistent with values from the gold standard method, R² = 0.9976. Coefficients of variation (CVs) ranged from 0.09% to 0.83%. For the volunteer study, T1 values from gray and white matter regions were consistent with literature values, and peaks of gray and white matter can be clearly delineated on whole-brain T1 histograms. CVs ranged from 0.01% to 2.30%. The acoustic noise measured at the scanner isocenter was 2.6 dBA higher compared to the in-bore background. Rapid and with low acoustic noise, the proposed method is shown to produce accurate T1 and PD maps with high repeatability by reconstructing sparsely sampled T1w images at different TIs using temporal subspace. Our approach can greatly enhance patient comfort during examination and therefore increase the acceptance of the procedure.

3D phase-resolved functional lung ventilation MR imaging in healthy volunteers and patients with chronic pulmonary disease

PURPOSE

To test the feasibility of 3D phase-resolved functional lung (PREFUL) MRI in healthy volunteers and patients with chronic pulmonary disease, to compare 3D to 2D PREFUL, and to investigate the required temporal resolution to obtain stable 3D PREFUL measurement.

METHODS

Sixteen participants underwent MRI using 2D and 3D PREFUL. Retrospectively, the spatial resolution of 3D PREFUL (4 × 4 × 4 mm³ ) was decreased to match the spatial resolution of 2D PREFUL (4 × 4 × 15 mm³ ), abbreviated as 3D_lowres . In addition to regional ventilation (RVent), flow-volume loops were computed and rated by a cross-correlation (CC). Ventilation defect percentage (VDP) maps were obtained. RVent, CC, VDP_RVent , and VDP_CC were compared for systematic differences between 2D, 3D_lowres , and 3D PREFUL. Dividing the 3D PREFUL data into 4- (≈ 20 phases), 8- (≈ 40 phases), and 12-min (≈ 60 phases) acquisition pieces, the ventilation parameter maps, including the heterogeneity of ventilation time to peak, were tested regarding the required temporal resolution.

RESULTS

RVent, CC, VDP_RVent , and VDP_CC  presented significant correlations between 2D and 3D PREFUL (r = 0.64-0.94). CC and VDP_CC  of 2D and 3D_lowres  PREFUL were significantly different (P < .0113). Comparing 3D_lowres  and 3D PREFUL, all parameters were found to be statistically different (P < .0045).

CONCLUSION

3D PREFUL MRI depicts the whole lung volume and breathing cycle with superior image resolution and with likely more precision compared to 2D PREFUL. Furthermore, 3D PREFUL is more sensitive to detect regions of hypoventilation and ventilation heterogeneity compared to 3D_lowres  PREFUL, which is important for early detection and improved monitoring of patients with chronic lung disease.

Rotated spiral RARE for high spatial and temporal resolution volumetric arterial spin labeling acquisition

Background:

Arterial Spin Labeling (ASL) MRI can provide quantitative images that are sensitive to both time averaged blood flow and its temporal fluctuations. 3D image acquisitions for ASL are desirable because they are more readily compatible with background suppression to reduce noise, can reduce signal loss and distortion, and provide uniform flow sensitivity across the brain. However, single-shot 3D acquisition for maximal temporal resolution typically involves degradation of image quality through blurring or noise amplification by parallel imaging. Here, we report a new approach to accelerate a common stack of spirals 3D image acquisition by pseudo golden-angle rotation and compressed sensing reconstruction without any degradation of time averaged blood flow images.

Methods:

28 healthy volunteers were imaged at 3T with background-suppressed unbalanced pseudo-continuous ASL combined with a pseudo golden-angle Stack-of-Spirals 3D RARE readout. A fully-sampled perfusion-weighted volume was reconstructed by 3D non-uniform Fast Fourier Transform (nuFFT) followed by sum-of-squares combination of the 32 individual channels. Coil sensitivities were estimated followed by reconstruction of the 39 single-shot volumes using an L₁-wavelet Compressed-Sensing reconstruction. Finally, brain connectivity analyses were performed in regions where BOLD signal suffers from low signal-to-noise ratio and susceptibility artifacts.

Results:

Image quality, assessed with a non-reference 3D blurring metric, of full time averaged blood flow was comparable to a conventional interleaved acquisition. The temporal resolution provided by the acceleration enabled identification and quantification of resting-state networks even in inferior regions such as the amygdala and inferior frontal lobes, where susceptibility artifacts can degrade conventional resting-state fMRI acquisitions.

Conclusion:

This approach can provide measures of blood flow modulations and resting-state networks for free within any research or clinical protocol employing ASL for resting blood flow.

A 3D hybrid-shot spiral sequence for hyperpolarized 13 C imaging

Purpose

Hyperpolarized imaging experiments have conflicting requirements of high spatial, temporal, and spectral resolution. Spectral-spatial RF excitation has been shown to form an attractive magnetization-efficient method for hyperpolarized imaging, but the optimum readout strategy is not yet known.

Methods

In this work, we propose a novel 3D hybrid-shot spiral sequence which features two constant density regions that permit the retrospective reconstruction of either high spatial or high temporal resolution images post hoc, (adaptive spatiotemporal imaging) allowing greater flexibility in acquisition and reconstruction.

Results

We have implemented this sequence, both via simulation and on a preclinical scanner, to demonstrate its feasibility, in both a ¹H phantom and with hyperpolarized ¹³C pyruvate in vivo.

Conclusions

This sequence forms an attractive method for acquiring hyperpolarized imaging datasets, providing adaptive spatiotemporal imaging to ameliorate the conflict of spatial and temporal resolution, with significant potential for clinical translation.

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.

3D pediatric cranial bone imaging using high-resolution MRI for visualizing cranial sutures: a pilot study

OBJECTIVE

There is an unmet need to perform imaging in young children and obtain CT-equivalent cranial bone images without subjecting the patients to radiation. In this study, the authors propose using a high-resolution fast low-angle shot golden-angle 3D stack-of-stars radial volumetric interpolated breath-hold examination (GA-VIBE) MRI sequence that is intrinsically robust to motion and has enhanced bone versus soft-tissue contrast.

METHODS

Patients younger than 11 years of age, who underwent clinical head CT scanning for craniosynostosis or other cranial malformations, were eligible for the study. 3D reconstructed images created from the GA-VIBE MRI sequence and the gold-standard CT scan were randomized and presented to 3 blinded reviewers. For all image sets, each reviewer noted the presence or absence of the 6 primary cranial sutures and recorded on 5-point Likert scales whether they recommended a second scan be performed.

RESULTS

Eleven patients (median age 1.8 years) underwent MRI after clinical head CT scanning was performed. Five of the 11 patients were sedated. Three clinicians reviewed the images, and there were no cases, either with CT scans or MR images, in which a reviewer agreed a repeat scan was required for diagnosis or surgical planning. The reviewers reported clear imaging of the regions of interest on 99% of the CT reviews and 96% of the MRI reviews. With CT as the standard, the sensitivity and specificity of the GA-VIBE MRI sequence to detect suture closure were 97% and 96%, respectively (n = 198 sutures read).

CONCLUSIONS

The 3D reconstructed images using the GA-VIBE sequence in comparison to the CT scans created clinically acceptable cranial images capable of detecting cranial sutures. Future directions include reducing the scan time, improving motion correction, and automating postprocessing for clinical utility.

Variable anisotropic FOV for 3D radial imaging with spiral phyllotaxis (VASP)

PURPOSE

To develop a new 3D radial trajectory based on the natural spiral phyllotaxis (SP), with variable anisotropic FOV.

THEORY & METHODS

A 3D radial trajectory based on the SP with favorable interleaving properties for cardiac imaging has been proposed by Piccini et al (Magn Reson Med. 2011;66:1049-1056), which supports a FOV with a fixed anisotropy. However, a fixed anisotropy can be inefficient when sampling objects with different anisotropic dimensions. We extend Larson's 3D radial method to provide variable anisotropic FOV for spiral phyllotaxis (VASP). Simulations were performed to measure distance between successive projections, analyze point spread functions, and compare aliasing artifacts for both VASP and conventional SP. VASP was fully implemented on a whole-body clinical MR scanner. Phantom and in vivo cardiac images were acquired at 1.5 tesla.

RESULTS

Simulations, phantom, and in vivo experiments confirmed that VASP can achieve variable anisotropic FOV while maintaining the favorable interleaving properties of SP. For an anisotropic FOV with 100:100:35 ratio, VASP required ~65% fewer radial projections than the conventional SP to satisfy Nyquist criteria. Alternatively, when the same number of radial projections were used as in conventional SP, VASP produced fewer aliasing artifacts for anisotropic objects within the excited imaging volumes.

CONCLUSION

We have developed a new method (VASP), which enables variable anisotropic FOV for 3D radial trajectory with SP. For anisotropic objects within the excited imaging volumes, VASP can reduce scan times and/or reduce aliasing artifacts.

Frontiers of Sodium MRI Revisited: From Cartilage to Brain Imaging

Sodium magnetic resonance imaging (²³ Na-MRI) is a highly promising imaging modality that offers the possibility to noninvasively quantify sodium content in the tissue, one of the most relevant parameters for biochemical investigations. Despite its great potential, due to the intrinsically low signal-to-noise ratio (SNR) of sodium imaging generated by low in vivo sodium concentrations, low gyromagnetic ratio, and substantially shorter relaxation times than for proton (¹ H) imaging, ²³ Na-MRI is extremely challenging. In this article, we aim to provide a comprehensive overview of the literature that has been published in the last 10-15 years and which has demonstrated different technical designs for a range of ²³ Na-MRI methods applicable for disease diagnoses and treatment efficacy evaluations. Currently, a wider use of 3.0T and 7.0T systems provide imaging with the expected increase in SNR and, consequently, an increased image resolution and a reduced scanning time. A great interest in translational research has enlarged the field of sodium MRI applications to almost all parts of the body: articular cartilage tendons, spine, heart, breast, muscle, kidney, and brain, etc., and several pathological conditions, such as tumors, neurological and degenerative diseases, and others. The quantitative parameter, tissue sodium concentration, which reflects changes in intracellular sodium concentration, extracellular sodium concentration, and intra-/extracellular volume fractions is becoming acknowledged as a reliable biomarker. Although the great potential of this technique is evident, there must be steady technical development for ²³ Na-MRI to become a standard imaging tool. The future role of sodium imaging is not to be considered as an alternative to ¹ H MRI, but to provide early, diagnostically valuable information about altered metabolism or tissue function associated with disease genesis and progression. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 1.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

Sliding window reduced FOV reconstruction for real-time cardiac imaging

BACKGROUND

Current functional cardiovascular imaging protocols mostly rely on electrocardiogram (ECG) gating and breathholding. The resulting image quality can substantially suffer from insufficient patient cooperation or severe arrhythmia. Real-time imaging can mitigate these effects but requires highly accelerated techniques, usually relying on non-cartesian trajectories and Compressed Sensing (CS).

METHODS

We investigate a sliding window reduced field of view (FOV) Echo Planar Imaging (EPI) technique for real-time cardiac MRI. Segmented EPI has been combined with a subtraction technique for reducing the FOV in cardiac applications to the region of the beating heart. Residual respiratory motion, potentially impairing the image quality, has been addressed by continuous update of the static image fraction, which is derived from a low-temporal resolution sliding window reconstruction. For further acceleration, the proposed technique was combined with parallel imaging.

RESULTS

The sliding window reduced FOV technique was proven feasible to reconstruct images of diagnostic image quality at a temporal resolution of 36.5ms per image. Semi-quantitative evaluation of image quality showed significant improvement over the existing rFOV method (p=0.039). Derived functional parameters show comparable results as with the BH-CINE reference. However, a trend to a slight underestimation of the largest and smallest in-plane volumes is observed.

CONCLUSION

The proposed technique is feasible of providing real-time cardiac MRI with a temporal resolution better than 40ms without the need of computably complex reconstruction techniques.

Rapid Imaging: Recent Advances in Abdominal MRI for Reducing Acquisition Time and Its Clinical Applications

Magnetic resonance imaging (MRI) plays an important role in abdominal imaging. The high contrast resolution offered by MRI provides better lesion detection and its capacity to provide multiparametric images facilitates lesion characterization more effectively than computed tomography. However, the relatively long acquisition time of MRI often detrimentally affects the image quality and limits its accessibility. Recent developments have addressed these drawbacks. Specifically, multiphasic acquisition of contrast-enhanced MRI, free-breathing dynamic MRI using compressed sensing technique, simultaneous multi-slice acquisition for diffusion-weighted imaging, and breath-hold three-dimensional magnetic resonance cholangiopancreatography are recent notable advances in this field. This review explores the aforementioned state-of-the-art techniques by focusing on their clinical applications and potential benefits, as well as their likely future direction.

Delayed contrast dynamics as marker of regional impairment in pulmonary fibrosis using 5D MRI - a pilot study

OBJECTIVE

To analyse delayed contrast dynamics of fibrotic lesions in interstitial lung disease (ILD) using five dimensional (5D) MRI and to correlate contrast dynamics with disease severity.

METHODS

20 patients (mean age: 71 years; M:F, 13:7), with chronic fibrosing ILD: n = 12 idiopathic pulmonary fibrosis (IPF) and n = 8 non-IPF, underwent thin-section multislice CT as part of the standard diagnostic workup and additionally MRI of the lung. 2 min after contrast injection, a radial gradient echo sequence with golden-angle spacing was acquired during 5 min of free-breathing, followed by 5D image reconstruction. Disease was categorized as severe or non-severe according to CT morphological regional severity. For each patient, 10 lesions were analysed.

RESULTS

IPF lesions showed later peak enhancement compared to non-IPF (severe: p = 0.01, non-severe: p = 0.003). Severe lesions showed later peak enhancement compared to non-severe lesions, in non-IPF (p = 0.04), but not in IPF (p = 0.35). There was a tendency towards higher accumulation and washout rates in IPF compared to non-IPF in non-severe disease. Severe lesions had lower washout rate than non-severe ones in both IPF (p = 0.003) and non-IPF (p = 0.005). Continuous contrast agent accumulation, without washout, was found only in IPF lesions.

CONCLUSIONS

Contrast agent dynamics are influenced by type and severity of pulmonary fibrosis, which might enable a more thorough characterisation of disease burden. The regional impairment is of particular interest in the context of antifibrotic treatments and was characterised using a non-invasive, non-irradiating, free-breathing method.

ADVANCES IN KNOWLEDGE

Delayed contrast enhancement patterns allow the assessment of regional lung impairment which could represent different disease stages or phenotypes in ILD.

Neural networks-based regularization for large-scale medical image reconstruction

In this paper we present a generalized Deep Learning-based approach for solving ill-posed large-scale inverse problems occuring in medical image reconstruction. Recently, Deep Learning methods using iterative neural networks (NNs) and cascaded NNs have been reported to achieve state-of-the-art results with respect to various quantitative quality measures as PSNR, NRMSE and SSIM across different imaging modalities. However, the fact that these approaches employ the application of the forward and adjoint operators repeatedly in the network architecture requires the network to process the whole images or volumes at once, which for some applications is computationally infeasible. In this work, we follow a different reconstruction strategy by strictly separating the application of the NN, the regularization of the solution and the consistency with the measured data. The regularization is given in the form of an image prior obtained by the output of a previously trained NN which is used in a Tikhonov regularization framework. By doing so, more complex and sophisticated network architectures can be used for the removal of the artefacts or noise than it is usually the case in iterative NNs. Due to the large scale of the considered problems and the resulting computational complexity of the employed networks, the priors are obtained by processing the images or volumes as patches or slices. We evaluated the method for the cases of 3D cone-beam low dose CT and undersampled 2D radial cine MRI and compared it to a total variation-minimization-based reconstruction algorithm as well as to a method with regularization based on learned overcomplete dictionaries. The proposed method outperformed all the reported methods with respect to all chosen quantitative measures and further accelerates the regularization step in the reconstruction by several orders of magnitude.

Echo Time-Dependence of Observed Lung T1 in Patients With Cystic Fibrosis and Correlation With Clinical Metrics

BACKGROUND

Noninvasive monitoring of early abnormalities and therapeutic intervention in cystic fibrosis (CF) lung disease using MRI is important. Lung T₁ mapping has shown potential for local functional imaging without contrast material. Recently, it was discovered that observed lung T₁ depends on the measurement echo time (TE).

PURPOSE

To examine TE-dependence of observed T₁ in patients with CF and its correlation with clinical metrics.

STUDY TYPE

Prospective.

POPULATION

In all, 75 pediatric patients with CF (8.6 ± 6.1 years, range 0.1-23 years), with 32 reexamined after 1 year.

FIELD STRENGTH/SEQUENCE

Patients were examined at 1.5T using an established MRI protocol and a multiecho inversion recovery 2D ultrashort echo time (UTE) sequence for T₁ (TE) mapping at five TEs including TE₁ = 70 μs.

ASSESSMENT

Morphological and perfusion MRI were assessed by a radiologist (M.W.) with 11 years of experience using an established CF-MRI scoring system. T₁ (TE) was quantified automatically. Clinical data including spirometry (FEV1pred%) and lung clearance index (LCI) were collected.

STATISTICAL TESTS

T₁ (TE) was correlated with the CF-MRI score, clinical data, and LCI.

RESULTS

T₁ (TE) showed a different curvature in CF than in healthy adults: T₁ at TE₁ was shorter in CF (1157 ms ± 73 ms vs. 1047 ms ± 70 ms, P < 0.001), but longer at TE₃ (1214 ms ± 72 ms vs. 1314 ms ± 68 ms, P < 0.001) and later TEs. The correlations of T₁ (TE) with patient age (ρ_TE1-TE5 = -0.55, -0.44, -0.24, -0.30, -0.22), and LCI (ρ_TE1-TE5 = -0.43, -0.42, -0.33, 0.27, -0.22) were moderate at ultra-short to short TE (P < 0.001) but decreased for longer TE. Moderate but similar correlations at all TE were found with MRI perfusion score (ρ_TE1-TE5 = -0.43, -0.51, -0.47, -0.46, -0.44) and FEV1pred% (ρ_TE1-TE5 = +0.44, +0.44, +0.43, +0.40, +0.39) (P < 0.05).

DATA CONCLUSION

TE should be considered when measuring lung T₁ , since observed differences between CF and healthy subjects strongly depend on TE. The different variation of correlation coefficients with TE for structural vs. functional metrics implies that TE-dependence holds additional information which may help to discern effects of tissue structural abnormalities and abnormal perfusion.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1 J. MAGN. RESON. IMAGING 2020;52:1645-1654.

Feasibility of real-time cardiac MRI in mice using tiny golden angle radial sparse

Cardiovascular magnetic resonance imaging has proven valuable for the assessment of structural and functional cardiac abnormalities. Even although it is an established imaging method in small animals, the long acquisition times of gated or self-gated techniques still limit its widespread application. In this study, the application of tiny golden angle radial sparse MRI (tyGRASP) for real-time cardiac imaging was tested in 12 constitutive nexilin (Nexn) knock-out (KO) mice, both heterozygous (Het, N = 6) and wild-type (WT, N = 6), and the resulting functional parameters were compared with a well-established self-gating approach. Real-time images were reconstructed for different temporal resolutions of between 16.8 and 79.8 ms per image. The suggested approach was additionally tested for dobutamine stress and qualitative first-pass perfusion imaging. Measurements were repeated twice within 2 weeks for reproducibility assessment. In direct comparison with the high-quality, self-gated technique, the real-time approach did not show any significant differences in global function parameters for acquisition times below 50 ms (rest) and 31.5 ms (stress). Compared with WT, the end-diastolic volume (EDV) and end-systolic volume (ESV) were markedly higher (P < 0.05) and the ejection fraction (EF) was significantly lower in the Het Nexn-KO mice at rest (P < 0.001). For the stress investigation, a clear decrease of EDV and ESV, and an increase in EF, but maintained stroke volume, could be observed in both groups. Combined with ECG-triggering, tyGRASP provided first-pass perfusion data with a temporal resolution of one image per heartbeat, allowing the quantitative assessment of upslope curves in the blood-pool and myocardium. Excellent inter-study reproducibility was achieved in all the functional parameters. The tyGRASP is a valuable real-time MRI technique for mice, which significantly reduces the scan time in preclinical cardiac functional imaging, providing sufficient image quality for deriving accurate functional parameters, and has the potential to investigate real-time and beat-to-beat changes.

Generalization of three-dimensional golden-angle radial acquisition to reduce eddy current artifacts in bSSFP CMR imaging

Purpose

We propose a novel generalization of the three-dimensional double-golden-angle profile ordering, which allows for whole-heart volumetric imaging with retrospective binning and reduced eddy current artifacts.

Methods

A novel theory bridging the gap between the three-dimensional double golden-angle trajectory, and the two-dimensional tiny-golden-angle trajectory is presented. This enables a class of double golden-angle profile orderings with a smaller angular distance between successive k-space readouts. The novel profile orderings were evaluated through simulations, phantom experiments, and in vivo imaging. Comparisons were made to the original double-golden-angle trajectory. Image uniformity and off-resonance sensitivity were evaluated using phantom measurements, and qualitative image quality was assessed using in vivo images acquired in a healthy volunteer.

Results

The proposed theory successfully reduced the angular step while maintaining image uniformity after binning. Simulations revealed a slow degradation with decreasing angular steps and an increasing number of physiological bins. The phantom images showed a definite improvement in image uniformity and increased robustness to off-resonance, and in vivo imaging corroborated those findings.

Conclusion

Reducing the angular step in cardio-respiratory-binned golden-angle imaging shows potential for overcoming eddy current-induced image artifacts associated with 3D golden-angle radial imaging.

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

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

Impact of respiratory motion correction on lesion visibility and quantification in thoracic PET/MR imaging

The impact of a method for MR-based respiratory motion correction of PET data on lesion visibility and quantification in patients with oncologic findings in the lung was evaluated. Twenty patients with one or more lesions in the lung were included. Hybrid imaging was performed on an integrated PET/MR system using 18F-FDG as radiotracer. The standard thoracic imaging protocol was extended by a free-breathing self-gated acquisition of MR data for motion modelling. PET data was acquired simultaneously in list-mode for 5-10 mins. One experienced radiologist and one experienced nuclear medicine specialist evaluated and compared the post-processed data in consensus regarding lesion visibility (scores 1-4, 4 being best), image noise levels (scores 1-3, 3 being lowest noise), SUVmean and SUVmax. Motion-corrected (MoCo) images were additionally compared with gated images. Non-motion-corrected free-breathing data served as standard of reference in this study. Motion correction generally improved lesion visibility (3.19 ± 0.63) and noise ratings (2.95 ± 0.22) compared to uncorrected (2.81 ± 0.66 and 2.95 ± 0.22, respectively) or gated PET data (2.47 ± 0.93 and 1.30 ± 0.47, respectively). Furthermore, SUVs (mean and max) were compared for all methods to estimate their respective impact on the quantification. Deviations of SUVmax were smallest between the uncorrected and the MoCo lesion data (average increase of 9.1% of MoCo SUVs), while SUVmean agreed best for gated and MoCo reconstructions (MoCo SUVs increased by 1.2%). The studied method for MR-based respiratory motion correction of PET data combines increased lesion sharpness and improved lesion activity quantification with high signal-to-noise ratio in a clinical setting. In particular, the detection of small lesions in moving organs such as the lung and liver may thus be facilitated. These advantages justify the extension of the PET/MR imaging protocol by 5-10 minutes for motion correction.

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.

MRI Reconstruction Using Markov Random Field and Total Variation as Composite Prior

Reconstruction of magnetic resonance images (MRI) benefits from incorporating a priori knowledge about statistical dependencies among the representation coefficients. Recent results demonstrate that modeling intraband dependencies with Markov Random Field (MRF) models enable superior reconstructions compared to inter-scale models. In this paper, we develop a novel reconstruction method, which includes a composite prior based on an MRF model and Total Variation (TV). We use an anisotropic MRF model and propose an original data-driven method for the adaptive estimation of its parameters. From a Bayesian perspective, we define a new position-dependent type of regularization and derive a compact reconstruction algorithm with a novel soft-thresholding rule. Experimental results show the effectiveness of this method compared to the state of the art in the field.

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.

CArtesian sampling with Variable density and Adjustable temporal resolution (CAVA)

PURPOSE

To develop a variable density Cartesian sampling method that allows retrospective adjustment of temporal resolution for dynamic MRI applications and to validate it in real-time phase contrast MRI (PC-MRI).

THEORY AND METHODS

The proposed method, called CArtesian sampling with Variable density and Adjustable temporal resolution (CAVA), begins by producing a sequence of phase encoding indices based on the golden ratio increment. Then, variable density is introduced by nonlinear stretching of the indices. Finally, the elements of the resulting sequence are rounded up to the nearest integer. The performance of CAVA is evaluated using PC-MRI data from a pulsatile flow phantom and real-time, free-breathing data from ten healthy volunteers.

RESULTS

CAVA enabled image recovery at various temporal resolutions that were selected retrospectively. For the pulsatile flow phantom, image quality and flow quantification accuracy from CAVA were comparable to that from another pseudo-random sampling pattern with fixed temporal resolution. In addition, flow quantification results based on CAVA were in good agreement with a breath-held segmented acquisition.

CONCLUSIONS

By allowing retrospective binning of the MRI data, CAVA provides an avenue to retrospectively adjust the temporal resolution of PC-MRI.

MRI evaluation of cerebrovascular reactivity in obstructive sleep apnea

Obstructive sleep apnea (OSA) is characterized by intermittent obstruction of the airways during sleep. Cerebrovascular reactivity (CVR) is an index of cerebral vessels' ability to respond to a vasoactive stimulus, such as increased CO₂. We hypothesized that OSA alters CVR, expressed as a breath-hold index (BHI) defined as the rate of change in CBF or BOLD signal during a controlled breath-hold stimulus mimicking spontaneous apneas by being both hypercapnic and hypoxic. In 37 OSA and 23 matched non sleep apnea (NSA) subjects, we obtained high temporal resolution CBF and BOLD MRI data before, during, and between five consecutive BH stimuli of 24 s, each averaged to yield a single BHI value. Greater BHI was observed in OSA relative to NSA as derived from whole-brain CBF (78.6 ± 29.6 vs. 60.0 ± 20.0 mL/min²/100 g, P = 0.010) as well as from flow velocity in the superior sagittal sinus (0.48 ± 0.18 vs. 0.36 ± 0.10 cm/s², P = 0.014). Similarly, BOLD-based BHI was greater in OSA in whole brain (0.19 ± 0.08 vs. 0.15 ± 0.03%/s, P = 0.009), gray matter (0.22 ± 0.09 vs. 0.17 ± 0.03%/s, P = 0.011), and white matter (0.14 ± 0.06 vs. 0.10 ± 0.02%/s, P = 0.010). The greater CVR is not currently understood but may represent a compensatory mechanism of the brain to maintain oxygen supply during intermittent apneas.

Simultaneous proton magnetic resonance fingerprinting and sodium MRI

PURPOSE

The goal of this work is to demonstrate a method for the simultaneous acquisition of proton multiparametric maps (T1 , T2 , and proton density) and sodium density images in 1 single scan. We hope that the development of such capabilities will help to ease the implementation of sodium MRI in clinical trials and provide more opportunities for researchers to investigate metabolism through sodium MRI.

METHODS

We developed a sequence based on magnetic resonance fingerprinting (MRF), which contains interleaved proton (1 H) and sodium (23 Na) excitations followed by a simultaneous center-out radial readout for both nuclei. The receive chain of a 7T scanner was modified to enable simultaneous acquisition of 1 H and 23 Na signal. The obtained signal-to-noise ratio (SNR) was evaluated, and the accuracy of both proton T1 , T2 , and B 1 + and sodium density maps were verified in phantoms. Finally, the method was demonstrated in 2 healthy subjects.

RESULTS

The SNR obtained using the simultaneous measurement was almost identical to single-nucleus measurements (<1% change). Similarly, the proton T1 and T2 maps remained stable (normalized root mean square error in T1 ≈ 2.2%, in T2 ≈ 1.4%, and B 1 + ≈ 5.4%), which indicates that the proposed sequence and hardware have no significant effects on the signal from either nucleus. In vivo measurements corroborated these results and demonstrated the feasibility of our method for in vivo application.

CONCLUSIONS

With the proposed approach, we were able to simultaneously acquire sodium density images in addition to proton T1 , T2 , and B 1 + maps as well as proton density images.

Rosette Trajectories Enable Ungated, Motion-Robust, Simultaneous Cardiac and Liver T2 * Iron Assessment

BACKGROUND

Quantitative T₂ * MRI is the standard of care for the assessment of iron overload. However, patient motion corrupts T₂ * estimates.

PURPOSE

To develop and evaluate a motion-robust, simultaneous cardiac and liver T₂ * imaging approach using non-Cartesian, rosette sampling and a model-based reconstruction as compared to clinical-standard Cartesian MRI.

STUDY TYPE

Prospective.

PHANTOM/POPULATION

Six ferumoxytol-containing phantoms (26-288 μg/mL). Eight healthy subjects and 18 patients referred for clinically indicated iron overload assessment.

FIELD STRENGTH/SEQUENCE

1.5T, 2D Cartesian and rosette gradient echo (GRE) ASSESSMENT: GRE T₂ * values were validated in ferumoxytol phantoms. In healthy subjects, test-retest and spatial coefficient of variation (CoV) analysis was performed during three breathing conditions. Cartesian and rosette T₂ * were compared using correlation and Bland-Altman analysis. Images were rated by three experienced radiologists on a 5-point scale.

STATISTICAL TESTS

Linear regression, analysis of variance (ANOVA), and paired Student's t-testing were used to compare reproducibility and variability metrics in Cartesian and rosette scans. The Wilcoxon rank test was used to assess reader score comparisons and reader reliability was measured using intraclass correlation analysis.

RESULTS

Rosette R2* (1/T₂ *) was linearly correlated with ferumoxytol concentration (r² = 1.00) and not significantly different than Cartesian values (P = 0.16). During breath-holding, ungated rosette liver and heart T₂ * had lower spatial CoV (liver: 18.4 ± 9.3% Cartesian, 8.8% ± 3.4% rosette, P = 0.02, heart: 37.7% ± 14.3% Cartesian, 13.4% ± 1.7% rosette, P = 0.001) and higher-quality scores (liver: 3.3 [3.0-3.6] Cartesian, 4.7 [4.1-4.9] rosette, P = 0.005, heart: 3.0 [2.3-3] Cartesian, 4.5 [3.8-5.0] rosette, P = 0.005) compared to Cartesian values. During free-breathing and failed breath-holding, Cartesian images had very poor to average image quality with significant artifacts, whereas rosette remained very good, with minimal artifacts (P = 0.001).

DATA CONCLUSION

Rosette k-sampling with a model-based reconstruction offers a clinically useful motion-robust T₂ * mapping approach for iron quantification. J. MAGN. RESON. IMAGING 2020;52:1688-1698.

Direct measurement of atrioventricular valve regurgitant jets using 4D flow cardiovascular magnetic resonance is accurate and reliable for children with congenital heart disease: a retrospective cohort study

BACKGROUND

3D-time resolved flow (4DF) cardiovascular magnetic resonance (CMR) with retrospective analysis of atrioventricular valve regurgitation (AVVR) allows for internal validation by multiple direct and indirect methods. Limited data exist on direct measurement of AVVR by 4DF CMR in pediatric congenital heart disease (CHD). We aimed to validate direct measurement of the AVVR jet as accurate and reliable compared to the volumetric method (clinical standard by 2D CMR) and as a superior method of internal validation than the annular inflow method.

METHODS

We identified 44 consecutive patients with diverse CHD referred for evaluation of AVVR by CMR. 1.5 T or 3 T scanners, intravenous contrast, and a combination of parallel imaging and compressed sensing were used. Four methods of measuring AVVR volume (RVol) were used: volumetric method (VOL; the clinical standard) = stroke volume by 2D balanced steady-state free precession - semilunar valve forward flow (SLFF); annular inflow method (AIM) = atrioventricular valve forward flow [AVFF] - semilunar valve net flow (SLNF); and direct measurement (JET). AVFF was measured using static and retrospective valve tracking planes. SLFF, SLNF, AVFF, and JET were measured by 4DF phase contrast. Regurgitant fraction was calculated as [RVol/(RVol+SLNF)]× 100. Statistical methods included Spearman, Wilcoxon rank sum test/Student paired t-test, Bland Altman analysis, and intra-class coefficient (ICC), where appropriate.

RESULTS

Regurgitant fraction by JET strongly correlated with the indirect methods (VOL and AIM) (ρ = 0.73-0.80, p < 0.001) and was similar to VOL with a median difference (interquartile range) of - 1.5% (- 8.3-7.2%; p = 0.624). VOL had weaker correlations with AIM and JET (ρ = 0.69-0.73, p < 0.001). AIM underestimated RF by 3.6-6.9% compared to VOL and JET, p < 0.03. Intra- and inter- observer reliability were excellent for all methods (ICC 0.94-0.99). The mean (±standard deviation) inter-observer difference for VOL was 2.4% (±5.1%), p < 0.05.

CONCLUSIONS

In a diverse cohort of pediatric CHD, measurement of AVVR using JET is accurate and reliable to VOL and is a superior method of internal validation compared to AIM. This study supports use of 4DF CMR for measurement of AVVR, obviating need for expert prospective prescription during image acquisition by 2D CMR.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain

PURPOSE

To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF).

METHODS

Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps.

RESULTS

A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T₁ and T₂ estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T₁ and 18 ms for T₂ values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.

CONCLUSION

We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

Highly accelerated, real-time phase-contrast MRI using radial k-space sampling and GROG-GRASP reconstruction: a feasibility study in pediatric patients with congenital heart disease

Retrospective electrocardiogram-gated, 2D phase-contrast (PC) flow MRI is routinely used in clinical evaluation of valvular/vascular disease in pediatric patients with congenital heart disease (CHD). In patients not requiring general anesthesia, clinical standard PC is conducted with free breathing for several minutes per slice with averaging. In younger patients under general anesthesia, clinical standard PC is conducted with breath-holding. One approach to overcome this limitation is using either navigator gating or self-navigation of respiratory motion, at the expense of lengthening scan times. An alternative approach is using highly accelerated, free-breathing, real-time PC (rt-PC) MRI, which to date has not been evaluated in CHD patients. The purpose of this study was to develop a 38.4-fold accelerated 2D rt-PC pulse sequence using radial k-space sampling and compressed sensing with 1.5 × 1.5 × 6.0 mm³ nominal spatial resolution and 40 ms nominal temporal resolution, and evaluate whether it is capable of accurately measuring flow in 17 pediatric patients (aortic valve, pulmonary valve, right and left pulmonary arteries) compared with clinical standard 2D PC (either breath-hold or free breathing). For clinical translation, we implemented an integrated reconstruction pipeline capable of producing DICOMs of the order of 2 min per time series (46 frames). In terms of association, forward volume, backward volume, regurgitant fraction, and peak velocity at peak systole measured with standard PC and rt-PC were strongly correlated (R² > 0.76; P < 0.001). Compared with clinical standard PC, in terms of agreement, forward volume (mean difference = 1.4% (3.0% of mean)) and regurgitant fraction (mean difference = -2.5%) were in good agreement, whereas backward volume (mean difference = -1.1 mL (28.2% of mean)) and peak-velocity at peak systole (mean difference = -21.3 cm/s (17.2% of mean)) were underestimated by rt-PC. This study demonstrates that the proposed rt-PC with the said spatial resolution and temporal resolution produces relatively accurate forward volumes and regurgitant fractions but underestimates backward volumes and peak velocities at peak systole in pediatric patients with CHD.

Free-breathing fat and R2 * quantification in the liver using a stack-of-stars multi-echo acquisition with respiratory-resolved model-based reconstruction

Purpose

To develop a free‐breathing hepatic fat and R2∗ quantification method by extending a previously described stack‐of‐stars model‐based fat‐water separation technique with additional modeling of the transverse relaxation rate R2∗.

Methods

The proposed technique combines motion‐robust radial sampling using a stack‐of‐stars bipolar multi‐echo 3D GRE acquisition with iterative model‐based fat‐water separation. Parallel‐Imaging and Compressed‐Sensing principles are incorporated through modeling of the coil‐sensitivity profiles and enforcement of total‐variation (TV) sparsity on estimated water, fat, and R2∗ parameter maps. Water and fat signals are used to estimate the confounder‐corrected proton‐density fat fraction (PDFF). Two strategies for handling respiratory motion are described: motion‐averaged and motion‐resolved reconstruction. Both techniques were evaluated in patients (n = 14) undergoing a hepatobiliary research protocol at 3T. PDFF and R2∗ parameter maps were compared to a breath‐holding Cartesian reference approach.

Results

Linear regression analyses demonstrated strong (r > 0.96) and significant (P ≪ .01) correlations between radial and Cartesian PDFF measurements for both the motion‐averaged reconstruction (slope: 0.90; intercept: 0.07%) and the motion‐resolved reconstruction (slope: 0.90; intercept: 0.11%). The motion‐averaged technique overestimated hepatic R2∗ values (slope: 0.35; intercept: 30.2 1/s) compared to the Cartesian reference. However, performing a respiratory‐resolved reconstruction led to better R2∗ value consistency (slope: 0.77; intercept: 7.5 1/s).

Conclusions

The proposed techniques are promising alternatives to conventional Cartesian imaging for fat and R2∗ quantification in patients with limited breath‐holding capabilities. For accurate R2∗ estimation, respiratory‐resolved reconstruction should be used.

A Perspective on MR Fingerprinting

This article reviews the basic concept of MR fingerprinting (MRF) with the goal of highlighting MRF's key contributions, putting them in the context of other quantitative MRI literature, and refining MRF's terminology. The article discusses the robustness and flexibility of MRF's signature dictionary-matching reconstruction along with more advanced MRF reconstructions. A key feature of MRF is the lack of assumptions about the signal evolution, which gives scientists the flexibility to tailor sequences for their needs. The article argues that the concept of unique fingerprints does not capture the requirements for successful parameter mapping and that an analysis of the signal's derivatives with respect to biophysical parameters, such as relaxation times, is more informative, as it allows one to evaluate the efficiency of a pulse sequence. The article points at the source of MRF's efficiency, namely, flip angle variations at the time scale of the relaxation times, and reveals that MRF's advantages are strongest at long scan times, as required for 3D imaging. Further, it outlines how MRF's flexibility can be used to design mutually tailored pulse sequences and biophysical models with the goal of improving the reproducibility of parameter mapping biological tissue. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.