Triple echo steady-state (TESS) relaxometry

Improved gradient echo magnitude- and phase-based mapping of T 2 $$ {\mathrm{T}}_2 $$ using multiple RF spoiling increments at 3T and 7T

PURPOSE

The transverse relaxation time T  2 $$ {}_2 $$ holds significant relevance in clinical applications and research studies. Conventional T  2 $$ {}_2 $$ mapping approaches rely on spin-echo sequences, which require lengthy acquisition times and involve high radiofrequency (RF) power deposition. An alternative gradient echo (GRE) phase-based T  2 $$ {}_2 $$ mapping method, utilizing steady-state acquisitions at one small RF spoil phase increment, was recently demonstrated. Here, a modified magnitude- and phase-based T  2 $$ {}_2 $$ mapping approach is proposed, which improvesT 2 $$ {\mathrm{T}}_2 $$ estimations by simultaneous fitting ofT 1 $$ {\mathrm{T}}_1 $$ and signal amplitude (A ∝ P D $$ A\propto PD $$ ) at three or more RF spoiling phase increments, instead of assuming a fixedT 1 $$ {\mathrm{T}}_1 $$ value.

METHODS

The feasibility of the magnitude-phase-based method was assessed by simulations, in phantom and in vivo experiments using skipped-CAIPI three-dimensional-echo-planar imaging (3D-EPI) for rapid GRE imaging.T 2 $$ {\mathrm{T}}_2 $$ ,T 1 $$ {\mathrm{T}}_1 $$ and PD estimations obtained by our method were compared toT 2 $$ {\mathrm{T}}_2 $$ of the phase-based method andT 1 $$ {\mathrm{T}}_1 $$ and PD of spoiled GRE-based multi-parameter mapping using a multi-echo version of the same sequence.

RESULTS

The agreement of the proposedT 2 $$ {\mathrm{T}}_2 $$ with ground truth and referenceT 2 $$ {\mathrm{T}}_2 $$ values was higher than that of phase-basedT 2 $$ {\mathrm{T}}_2 $$ in simulations and in phantom data. While phase-basedT 2 $$ {\mathrm{T}}_2 $$ overestimation increases with actualT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ , the proposed method is accurate over a large range of physiologically meaningfulT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ values. At the same time, precision is improved. In vivo results were in line with these observations.

CONCLUSION

Accurate magnitude-phase-based T  2 $$ {}_2 $$ mapping is feasible in less than 5 min scan time for 1 mm nominal isotropic whole-head coverage at 3T and 7T.

Cartilage compositional MRI-a narrative review of technical development and clinical applications over the past three decades

Articular cartilage damage and degeneration are among hallmark manifestations of joint injuries and arthritis, classically osteoarthritis. Cartilage compositional MRI (Cart-C MRI), a quantitative technique, which aims to detect early-stage cartilage matrix changes that precede macroscopic alterations, began development in the 1990s. However, despite the significant advancements over the past three decades, Cart-C MRI remains predominantly a research tool, hindered by various technical and clinical hurdles. This paper will review the technical evolution of Cart-C MRI, delve into its clinical applications, and conclude by identifying the existing gaps and challenges that need to be addressed to enable even broader clinical application of Cart-C MRI.

Intra-scan RF power amplifier drift correction

PURPOSE

The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated.

METHODS

RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan.

RESULTS

RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability.

CONCLUSION

Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Fast and accessible T2 mapping using off-resonance corrected DESPOT2 with application to 3D prostate

PURPOSE

Most T1 and T2 mapping take long acquisitions or needs specialized sequences not widely accessible on clinical scanners. An available solution is DESPOT1/T2 (Driven equilibrium single pulse observation of T1/T2). DESPOT1/T2 uses Spoiled gradient-echo (SPGR) and balanced Steady-State Free Precession (bSSFP) sequences, offering an accessible and reliable way for 3D accelerated T1/T2 mapping. However, bSSFP is prone to off-resonance artifacts, limiting the application of DESPOT2 in regions with high susceptibility contrasts, like the prostate. Our proposal, DESPO+, employs the full bSSFP and SPGR models with a dictionary-based method to reconstruct 3D T1/T2 maps in the prostate region without off-resonance banding.

METHODS

DESPO+ modifies the bSSFP acquisition of the original variable flip angle DESPOT2. DESPO+ uses variable repetition and echo times, employing a dictionary-based method of the full bSSFP and SPGR models to reconstruct T1, T2, and Proton Density (PD) simultaneously. The proposed DESPO+ method underwent testing through simulations, T1/T2 phantoms, and on fourteen healthy subjects.

RESULTS

The results reveal a significant reduction in T2 map banding artifacts compared to the original DESPOT2 method. DESPO+ approach reduced T2 errors by up to seven times compared to DESPOT2 in simulations and phantom experiments. We also synthesized in-vivo T1-weighted/T2-weighted images from the acquired maps using a spin-echo model to verify the map's quality when lacking a reference. For in-vivo imaging, the synthesized images closely resemble those from the clinical MRI protocol, reducing scan time by around 50% compared to traditional spin-echo T1-weighted/T2-weighted acquisitions.

CONCLUSION

DESPO+ provides an off-resonance insensitive and clinically available solution, enabling high-resolution 3D T1/T2 mapping and synthesized T1-weighted/T2-weighted images for the entire prostate, all achieved within a short scan time of 3.6 min, similar to DESPOT1/T2.

Robust T2 estimation with balanced steady state free precession

PURPOSE

To develop a novel signal representation for balanced steady state free precession (bSSFP) displaying its T2 independence on B1 and on magnetization transfer (MT) effects.

METHODS

A signal model for bSSFP is developed that shows only an explicit dependence (up to a scaling factor) on E2 (and, therefore, T2) and a novel parameter c (with implicit dependence on the flip angle and E1). Moreover, it is shown that MT effects, entering the bSSFP signal via a binary spin bath model, can be captured by a redefinition of T1 and, therefore, leading to modification of E1, resulting in the same signal model. Various sets of phase-cycled bSSFP brain scans (different flip angles, different TR, different RF pulse durations, and different number of phase cycles) were recorded at 3 T. The parameters T2 (E2) and c were estimated using a variable projection (VARPRO) method and Monte-Carlo simulations were performed to assess T2 estimation precision.

RESULTS

Initial experiments confirmed the expected independence of T2 on various protocol settings, such as TR, the flip angle, B1 field inhomogeneity, and the RF pulse duration. Any variation (within the explored range) appears to directly affect the estimation of the parameter c only-in agreement with theory.

CONCLUSION

BSSFP theory predicts an extraordinary feature that all MT and B1-related variational aspects do not enter T2 estimation, making it a potentially robust methodology for T2 quantification, pending validation against existing standards.

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.

Magnetic resonance coherence pathway unraveling

Efficiently acquiring multi-contrast magnetic resonance imaging data is crucial for patient comfort and clinical throughput. Developing scan acceleration methods tailored for specific applications drastically improves the value of an MRI examination. Here, we propose a novel method to control the aliasing of simultaneously acquired images of multiple spin echo coherence pathways with the goal of producing high quality multi-contrast images from a single acquisition. Modulating the radiofrequency phase of several pulses applied in brief succession also uniquely modulates the phase of spin echo coherence pathways. A method, termed magnetic resonance coherence pathway unraveling (MR-CPU), to control the aliasing of simultaneously acquired coherence pathway images is developed here along with parallel imaging-based reconstruction methods to separate them. MR-CPU was validated in phantom experiments and tested in vivo. High levels of correlation between reference pathway images and MR-CPU-derived coherence pathway images were found from the phantom experiments. Minimal artifacts arising from the separation of the overlapped coherence pathway images were observed in vivo. MR-CPU provides a novel mechanism through which to acquire and separate multiple overlapped coherence pathway images, thus adding to the diagnostic potential of an MRI exam without the penalty of additional scan time.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

Approximate B 1 + scaling of the SSFP steady state

PURPOSE

It is shown that the steady state of rapid, TR-periodic steady-state free precession (SSFP) sequences at small to moderate flip angles exhibits a universal, approximate scaling law with respect to variations of B 1 + $$ {B}_1^{+} $$ . Implications for the accuracy and precision of relaxometry experiments are discussed.

METHODS

The approximate scaling law is derived from and numerically tested against known analytical solutions. To assess the attainable estimator precision in a typical relaxometry experiment, we calculate the Cramér-Rao bound (CRB) and perform Monte Carlo (MC) simulations.

RESULTS

The approximate universal scaling holds well up to moderate flip angles. For pure steady state relaxometry, we observe a significant precision penalty for simultaneous estimation of R 1 $$ {R}_1 $$ and B 1 + $$ {B}_1^{+} $$ , whereas good R 2 $$ {R}_2 $$ estimates can be obtained without even knowing the correct actual flip angle.

CONCLUSION

Simultaneous estimation of R 1 $$ {R}_1 $$ and B 1 + $$ {B}_1^{+} $$ from a set of SSFP steady states alone is not advisable. Apart from separate B 1 + $$ {B}_1^{+} $$ measurements, the problem can be addressed by adding transient state information, but, depending on the situation, residual effects due to the scaling may still require some attention.

Quantitative MRI for Evaluation of Musculoskeletal Disease: Cartilage and Muscle Composition, Joint Inflammation, and Biomechanics in Osteoarthritis

ABSTRACT

Magnetic resonance imaging (MRI) is a valuable tool for evaluating musculoskeletal disease as it offers a range of image contrasts that are sensitive to underlying tissue biochemical composition and microstructure. Although MRI has the ability to provide high-resolution, information-rich images suitable for musculoskeletal applications, most MRI utilization remains in qualitative evaluation. Quantitative MRI (qMRI) provides additional value beyond qualitative assessment via objective metrics that can support disease characterization, disease progression monitoring, or therapy response. In this review, musculoskeletal qMRI techniques are summarized with a focus on techniques developed for osteoarthritis evaluation. Cartilage compositional MRI methods are described with a detailed discussion on relaxometric mapping (T 2 , T 2 *, T 1ρ ) without contrast agents. Methods to assess inflammation are described, including perfusion imaging, volume and signal changes, contrast-enhanced T 1 mapping, and semiquantitative scoring systems. Quantitative characterization of structure and function by bone shape modeling and joint kinematics are described. Muscle evaluation by qMRI is discussed, including size (area, volume), relaxometric mapping (T 1 , T 2 , T 1ρ ), fat fraction quantification, diffusion imaging, and metabolic assessment by 31 P-MR and creatine chemical exchange saturation transfer. Other notable technologies to support qMRI in musculoskeletal evaluation are described, including magnetic resonance fingerprinting, ultrashort echo time imaging, ultrahigh-field MRI, and hybrid MRI-positron emission tomography. Challenges for adopting and using qMRI in musculoskeletal evaluation are discussed, including the need for metal artifact suppression and qMRI standardization.

Primary Multiparametric Quantitative Brain MRI: State-of-the-Art Relaxometric and Proton Density Mapping Techniques

This review on brain multiparametric quantitative MRI (MP-qMRI) focuses on the primary subset of quantitative MRI (qMRI) parameters that represent the mobile ("free") and bound ("motion-restricted") proton pools. Such primary parameters are the proton densities, relaxation times, and magnetization transfer parameters. Diffusion qMRI is also included because of its wide implementation in complete clinical MP-qMRI application. MP-qMRI advances were reviewed over the past 2 decades, with substantial progress observed toward accelerating image acquisition and increasing mapping accuracy. Areas that need further investigation and refinement are identified as follows: (a) the biologic underpinnings of qMRI parameter values and their changes with age and/or disease and (b) the theoretical limitations implicitly built into most qMRI mapping algorithms that do not distinguish between the different spatial scales of voxels versus spin packets, the central physical object of the Bloch theory. With rapidly improving image processing techniques and continuous advances in computer hardware, MP-qMRI has the potential for implementation in a wide range of clinical applications. Currently, three emerging MP-qMRI applications are synthetic MRI, macrostructural qMRI, and microstructural tissue modeling.

Phase-based fast 3D high-resolution quantitative T2 MRI in 7 T human brain imaging

Magnetic resonance imaging (MRI) is a powerful and versatile technique that offers a range of physiological, diagnostic, structural, and functional measurements. One of the most widely used basic contrasts in MRI diagnostics is transverse relaxation time (T₂)-weighted imaging, but it provides only qualitative information. Realizing quantitative high-resolution T₂ mapping is imperative for the development of personalized medicine, as it can enable the characterization of diseases progression. While ultra-high-field (≥ 7 T) MRI offers the means to gain new insights by increasing the spatial resolution, implementing fast quantitative T₂ mapping cannot be achieved without overcoming the increased power deposition and radio frequency (RF) field inhomogeneity at ultra-high-fields. A recent study has demonstrated a new phase-based T₂ mapping approach based on fast steady-state acquisitions. We extend this new approach to ultra-high field MRI, achieving quantitative high-resolution 3D T₂ mapping at 7 T while addressing RF field inhomogeneity and utilizing low flip angle pulses; overcoming two main ultra-high field challenges. The method is based on controlling the coherent transverse magnetization in a steady-state gradient echo acquisition; achieved by utilizing low flip angles, a specific phase increment for the RF pulses, and short repetition times. This approach simultaneously extracts both T₂ and RF field maps from the phase of the signal. Prior to in vivo experiments, the method was assessed using a 3D head-shaped phantom that was designed to model the RF field distribution in the brain. Our approach delivers fast 3D whole brain images with submillimeter resolution without requiring special hardware, such as multi-channel transmit coil, thus promoting high usability of the ultra-high field MRI in clinical practice.

A relaxometry method that emphasizes practicality and availability

PURPOSE

Although many methods have been proposed to quantitatively map the main MRI parameters (e.g., T₁ , T₂ , C × M₀ ), these methods often involve special sequences not readily available on clinical scanners and/or may require long scan times. In contrast, the proposed method can readily run on most scanners, offer flexible tradeoffs between scan time and image quality, and map MRI parameters jointly to ensure spatial alignment.

METHODS

The approach is based on the multi-shot spin-echo (SE) EPI sequence. The corresponding signal equation was derived and strategies for solving it were developed. As usual with multi-shot EPI, scan time can readily be traded-off against image quality by adjusting the echo train length. Validation was performed against reference relaxometry methods, in gel phantoms with varying concentrations of gadobutrol and gadoterate meglumine contrast agents. In vivo examples are further presented, from 3 neuroradiology patients.

RESULTS

Bland-Altman analysis was performed: for T₂ , as compared to 2D SE, bias was 0.29 ms and the 95% limits of agreement ranged from -1.15 to +1.73 ms. For T₁ , compared to inversion-recovery SE (and MOLLI), bias was -20.2 ms (and -14.5 ms) and the limits of agreement ranged from -62.4 to +22.0 ms (and -53.8 to +24.9 ms). The mean relative T₁ error between the proposed method and each of the 2 reference methods was similar to that of the reference methods among themselves.

CONCLUSION

In the constellation of existing relaxometry methods, the proposed method is meant to stand out in terms of its practicality and availability.

BUDA-MESMERISE: Rapid acquisition and unsupervised parameter estimation for T1 , T2 , M0 , B0 , and B1 maps

PURPOSE

Rapid acquisition scheme and parameter estimation method are proposed to acquire distortion-free spin- and stimulated-echo signals and combine the signals with a physics-driven unsupervised network to estimate T₁ , T₂ , and proton density (M₀ ) parameter maps, along with B₀ and B₁ information from the acquired signals.

THEORY AND METHODS

An imaging sequence with three 90° RF pulses is utilized to acquire spin- and stimulated-echo signals. We utilize blip-up/-down acquisition to eliminate geometric distortion incurred by the effects of B₀ inhomogeneity on rapid EPI acquisitions. For multislice imaging, echo-shifting is applied to utilize dead time between the second and third RF pulses to encode information from additional slice positions. To estimate parameter maps from the spin- and stimulated-echo signals with high fidelity, 2 estimation methods, analytic fitting and a novel unsupervised deep neural network method, are developed.

RESULTS

The proposed acquisition provided distortion-free T₁ , T₂ , relative proton density (M0), B₀ , and B₁ maps with high fidelity both in phantom and in vivo brain experiments. From the rapidly acquired spin- and stimulated-echo signals, analytic fitting and the network-based method were able to estimate T₁ , T₂ , M₀ , B₀ , and B₁ maps with high accuracy. Network estimates demonstrated noise robustness owing to the fact that the convolutional layers take information into account from spatially adjacent voxels.

CONCLUSION

The proposed acquisition/reconstruction technique enabled whole-brain acquisition of coregistered, distortion-free, T₁ , T₂ , M₀ , B₀ , and B₁ maps at 1 × 1 × 5 mm³ resolution in 50 s. The proposed unsupervised neural network provided noise-robust parameter estimates from this rapid acquisition.

Reproducibility of an Automated Quantitative MRI Assessment of Low-Grade Knee Articular Cartilage Lesions

OBJECTIVE

The goal of this study was to assess the reproducibility of an automated knee cartilage segmentation of 21 cartilage regions with a model-based algorithm and to compare the results with manual segmentation.

DESIGN

Thirteen patients with low-grade femoral cartilage defects were included in the study and were scanned twice on a 7-T magnetic resonance imaging (MRI) scanner 8 days apart. A 3-dimensional double-echo steady-state (3D-DESS) sequence was used to acquire MR images for automated cartilage segmentation, and T2-mapping was performed using a 3D triple-echo steady-state (3D-TESS) sequence. Cartilage volume, thickness, and T2 and texture features were automatically extracted from each knee for each of the 21 subregions. DESS was used for manual cartilage segmentation and compared with automated segmentation using the Dice coefficient. The reproducibility of each variable was expressed using standard error of measurement (SEM) and smallest detectable change (SDC).

RESULTS

The Dice coefficient for the similarity between manual and automated segmentation ranged from 0.83 to 0.88 in different cartilage regions. Test-retest analysis of automated cartilage segmentation and automated quantitative parameter extraction revealed excellent reproducibility for volume measurement (mean SDC for all subregions of 85.6 mm³), for thickness detection (SDC = 0.16 mm) and also for T2 values (SDC = 2.38 ms) and most gray-level co-occurrence matrix features (SDC = 0.1 a.u.).

CONCLUSIONS

The proposed technique of automated knee cartilage evaluation based on the segmentation of 3D MR images and correlation with T2 mapping provides highly reproducible results and significantly reduces the segmentation effort required for the analysis of knee articular cartilage in longitudinal studies.

A unifying view on extended phase graphs and Bloch simulations for quantitative MRI

Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.

Efficiency analysis for quantitative MRI of T1 and T2 relaxometry methods

This study presents a comparison of quantitative MRI methods based on an efficiency metric that quantifies their intrinsic ability to extract information about tissue parameters. Under a regime of unbiased parameter estimates, an intrinsic efficiency metricηwas derived for fully-sampled experiments which can be used to both optimize and compare sequences. Here we optimize and compare several steady-state and transient gradient-echo based qMRI methods, such as magnetic resonance fingerprinting (MRF), for jointT1andT2mapping. The impact of undersampling was also evaluated, assuming incoherent aliasing that is treated as noise by parameter estimation.In vivovalidation of the efficiency metric was also performed. Transient methods such as MRF can be up to 3.5 times more efficient than steady-state methods, when spatial undersampling is ignored. If incoherent aliasing is treated as noise during least-squares parameter estimation, the efficiency is reduced in proportion to the SNR of the data, with reduction factors of 5 often seen for practical SNR levels.In vivovalidation showed a very good agreement between the theoretical and experimentally predicted efficiency. This work presents and validates an efficiency metric to optimize and compare the performance of qMRI methods. Transient methods were found to be intrinsically more efficient than steady-state methods, however the effect of spatial undersampling can significantly erode this advantage.

Multiple-echo steady-state (MESS): Extending DESS for joint T2 mapping and chemical-shift corrected water-fat separation

Purpose

To extend the double echo steady‐state (DESS) sequence to enable chemical‐shift corrected water‐fat separation.

Methods

This study proposes multiple‐echo steady‐state (MESS), a sequence that modifies the readouts of the DESS sequence to acquire two echoes each with bipolar readout gradients with higher readout bandwidth. This enables water‐fat separation and eliminates the need for water‐selective excitation that is often used in combination with DESS, without increasing scan time. An iterative fitting approach was used to perform joint chemical‐shift corrected water‐fat separation and T₂ estimation on all four MESS echoes simultaneously. MESS and water‐selective DESS images were acquired for five volunteers, and were compared qualitatively as well as quantitatively on cartilage T₂ and thickness measurements. Signal‐to‐noise ratio (SNR) and T₂ quantification were evaluated numerically using pseudo‐replications of the acquisition.

Results

The water‐fat separation provided by MESS was robust and with quality comparable to water‐selective DESS. MESS T₂ estimation was similar to DESS, albeit with slightly higher variability. Noise analysis showed that SNR in MESS was comparable to DESS on average, but did exhibit local variations caused by uncertainty in the water‐fat separation.

Conclusion

In the same acquisition time as DESS, MESS provides water‐fat separation with comparable SNR in the reconstructed water and fat images. By providing additional image contrasts in addition to the water‐selective DESS images, MESS provides a promising alternative to DESS.

Streamlined magnetic resonance fingerprinting: Fast whole-brain coverage with deep-learning based parameter estimation

Magnetic resonance fingerprinting (MRF) is a quantitative MRI (qMRI) framework that provides simultaneous estimates of multiple relaxation parameters as well as metrics of field inhomogeneity in a single acquisition. However, current challenges exist in the forms of (1) scan time; (2) need for custom image reconstruction; (3) large dictionary sizes; (4) long dictionary-matching time. This study aims to introduce a novel streamlined magnetic-resonance fingerprinting (sMRF) framework based on a single-shot echo-planar imaging (EPI) sequence to simultaneously estimate tissue T1, T2, and T2* with integrated B1⁺ correction. Encouraged by recent work on EPI-based MRF, we developed a method that combines spin-echo EPI with gradient-echo EPI to achieve T2 in addition to T1 and T2* quantification. To this design, we add simultaneous multi-slice (SMS) acceleration to enable full-brain coverage in a few minutes. Moreover, in the parameter-estimation step, we use deep learning to train a deep neural network (DNN) to accelerate the estimation process by orders of magnitude. Notably, due to the high image quality of the EPI scans, the training process can rely simply on Bloch-simulated data. The DNN also removes the need for storing large dictionaries. Phantom scans along with in-vivo multi-slice scans from seven healthy volunteers were acquired with resolutions of 1.1×1.1×3 mm³ and 1.7×1.7×3 mm³, and the results were validated against ground truth measurements. Excellent correspondence was found between our T1, T2, and T2* estimates and results obtained from standard approaches. In the phantom scan, a strong linear relationship (R = 1-1.04, R²>0.96) was found for all parameter estimates, with a particularly high agreement for T2 estimation (R²>0.99). Similar findings are reported for the in-vivo human data for all of our parameter estimates. Incorporation of DNN results in a reduction of parameter estimation time on the order of 1000 x and a reduction in storage requirements on the order of 2500 x while achieving highly similar results as conventional dictionary matching (%differences of 7.4 ± 0.4%, 3.6 ± 0.3% and 6.0 ± 0.4% error in T1, T2, and T2* estimation). Thus, sMRF has the potential to be the method of choice for future MRF studies by providing ease of implementation, fast whole-brain coverage, and ultra-fast T1/T2/T2* estimation.

Clinical implementation of accelerated T2 mapping: Quantitative magnetic resonance imaging as a biomarker for annular tear and lumbar disc herniation

OBJECTIVES

This study evaluates GRAPPATINI, an accelerated T2 mapping sequence combining undersampling and model-based reconstruction to facilitate the clinical implementation of T2 mapping of the lumbar intervertebral disc.

METHODS

Fifty-eight individuals (26 females, 32 males, age 23.3 ± 8.0 years) were prospectively examined at 3 T. This cohort study consisted of 19 patients, 20 rowers, and 19 volunteers. GRAPPATINI was conducted with the same parameters as a conventional 2D multi-echo spin-echo (MESE) sequence in 02:27 min instead of 13:18 min. Additional T2 maps were calculated after discarding the first echo (T2-WO1ST) and only using even echoes (T2-EVEN). Segmentation was done on the four most central slices. The resulting T2 values were compared for all four measurements.

RESULTS

T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of bulging discs or herniated discs (all p < 0.001). For the posterior annular region, only T2-GRAPPATINI showed a significant difference (p = 0.011) between normal and herniated discs. There was a significant difference between T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of discs with and without an annular tear for the nucleus pulposus (all p < 0.001). The nucleus pulposus' T2 at different degeneration states showed significant differences between all group comparisons of Pfirrmann grades for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

CONCLUSIONS

GRAPPATINI facilitates the use of T2 values as quantitative imaging biomarkers to detect disc pathologies such as degeneration, lumbar disc herniation, and annular tears while simultaneously shortening the acquisition time from 13:18 to 2:27 min.

KEY POINTS

• T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of discs with bulging or herniation (all p < 0.001). • The investigated T2 mapping techniques differed significantly in discs with and without annular tearing (all p < 0.001). • The nucleus pulposus' T2 showed significant differences between different stages of degeneration in all group comparisons for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

Alternating unbalanced SSFP for 3D R 2 ' mapping of the human brain

PURPOSE

Measuring the transverse-relaxation rate R 2 ' provides valuable information in quantitative evaluation of tissue microstructure, for example, in terms of oxygenation levels. Here, we propose an alternating unbalanced SSFP pulse sequence for rapid whole-brain 3D R 2 ' mapping.

METHODS

Unlike currently practiced, spin echo-based R 2 ' measurement techniques, the proposed method alternates between SSFP-FID and SSFP-ECHO modes for rapid 3D encoding of transverse relaxation rates expressed as R₂ + R 2 ' and R₂ - R 2 ' . Z-shimming gradients embedded into multi-echo trains of each SSFP module are designed to achieve relative immunity to large-scale magnetic-field variations (ΔB₀ ). Appropriate models for the temporal evolution of the two groups of SSFP signals were constructed with ΔB₀ -induced modulations accounted for, leading to ΔB₀ -corrected estimation of R₂ , R 2 ' , and R 2 ∗ (= R₂ + R 2 ' ). Additionally, relative magnetic susceptibility (Δχ) maps were obtained by quantitative susceptibility mapping of the phase data. Numerical simulations were performed to optimize scan parameters, followed by in vivo studies at 3 T in 7 healthy subjects. Measured parameters were evaluated in six brain regions, and subjected to interparameter correlation analysis.

RESULTS

The resultant maps of R 2 ' and additionally derived R₂ , R 2 ∗ , and Δχ all demonstrated the expected contrast across brain territories (eg, deep brain structures versus cortex), with the measured values in good agreement with previous reports. Furthermore, regression analyses yielded strong linear relationships for the transverse relaxation parameters ( R 2 ' , R₂ , and R 2 ∗ ) against Δχ.

CONCLUSION

Results suggest feasibility of the proposed method as a practical and reliable means for measuring R 2 ' , R₂ , R 2 ∗ , and Δχ across the entire brain.

Musculoskeletal MR Imaging Applications at Ultra-High (7T) Field Strength

Regulatory approval of ultrahigh field (UHF) MR imaging scanners for clinical use has opened new opportunities for musculoskeletal imaging applications. UHF MR imaging has unique advantages in terms of signal-to-noise ratio, contrast-to-noise ratio, spectral resolution, and multinuclear applications, thus providing unique information not available at lower field strengths. But UHF also comes with a set of technical challenges that are yet to be resolved and may not be suitable for all imaging applications. This review focuses on the latest research in musculoskeletal MR imaging applications at UHF including morphologic imaging, T₂, T₂∗, and T_1ρ mapping, chemical exchange saturation transfer, sodium imaging, and phosphorus spectroscopy imaging applications.

Configuration-based electrical properties tomography

PURPOSE

To introduce phase-based conductivity mapping from a configuration space analysis.

METHODS

The frequency response function of balanced SSFP (bSSFP) is used to perform a configuration space analysis. It is shown that the transceive phase for conductivity mapping can be directly obtained by a simple fast Fourier transform of a series of phase-cycled bSSFP scans. For validation, transceive phase and off-resonance mapping with fast Fourier transform is compared with phase estimation using a recently proposed method, termed PLANET. Experiments were performed in phantoms and for in vivo brain imaging at 3 T using a quadrature head coil.

RESULTS

For fast Fourier transform, aliasing can lead to systematic phase errors. This bias, however, decreases rapidly with increasing sampling points. Interestingly, Monte Carlo simulations revealed a lower uncertainty for the transceive phase and the off-resonance using fast Fourier transform as compared with PLANET. Both methods, however, essentially retrieve the same phase information from a set of phase-cycled bSSFP scans. As a result, configuration-based conductivity mapping was successfully performed using eight phase-cycled bSSFP scans in the phantoms and for brain tissues. Overall, the retrieved values were in good agreement with expectations. Conductivity estimation and mapping of the field inhomogeneities can therefore be performed in conjunction with the estimation of other quantitative parameters, such as relaxation, using configuration theory.

CONCLUSIONS

Phase-based conductivity mapping can be estimated directly from a simple Fourier analysis, such as in conjunction with relaxometry, using a series of phase-cycled bSSFP scans.

Musculoskeletal MRI at 7 T: do we need more or is it more than enough?

Ultra-high field magnetic resonance imaging (UHF-MRI) provides important diagnostic improvements in musculoskeletal imaging. The higher signal-to-noise ratio leads to higher spatial and temporal resolution which results in improved anatomic detail and higher diagnostic confidence. Several methods, such as T2, T2*, T1rho mapping, delayed gadolinium-enhanced, diffusion, chemical exchange saturation transfer, and magnetisation transfer techniques, permit a better tissue characterisation. Furthermore, UHF-MRI enables in vivo measurements by low-γ nuclei (²³Na, ³¹P, ¹³C, and ³⁹K) and the evaluation of different tissue metabolic pathways. European Union and Food and Drug Administration approvals for clinical imaging at UHF have been the first step towards a more routinely use of this technology, but some drawbacks are still present limiting its widespread clinical application. This review aims to provide a clinically oriented overview about the application of UHF-MRI in the different anatomical districts and tissues of musculoskeletal system and its pros and cons. Further studies are needed to consolidate the added value of the use of UHF-MRI in the routine clinical practice and promising efforts in technology development are already in progress.

Modeling an equivalent b-value in diffusion-weighted steady-state free precession

PURPOSE

Diffusion-weighted steady-state free precession (DW-SSFP) is shown to provide a means to probe non-Gaussian diffusion through manipulation of the flip angle. A framework is presented to define an effective b-value in DW-SSFP.

THEORY

The DW-SSFP signal is a summation of coherence pathways with different b-values. The relative contribution of each pathway is dictated by the flip angle. This leads to an apparent diffusion coefficient (ADC) estimate that depends on the flip angle in non-Gaussian diffusion regimes. By acquiring DW-SSFP data at multiple flip angles and modeling the variation in ADC for a given form of non-Gaussianity, the ADC can be estimated at a well-defined effective b-value.

METHODS

A gamma distribution is used to model non-Gaussian diffusion, embedded in the Buxton signal model for DW-SSFP. Monte-Carlo simulations of non-Gaussian diffusion in DW-SSFP and diffusion-weighted spin-echo sequences are used to verify the proposed framework. Dependence of ADC on flip angle in DW-SSFP is verified with experimental measurements in a whole, human postmortem brain.

RESULTS

Monte-Carlo simulations reveal excellent agreement between ADCs estimated with diffusion-weighted spin-echo and the proposed framework. Experimental ADC estimates vary as a function of flip angle over the corpus callosum of the postmortem brain, estimating the mean and standard deviation of the gamma distribution as 1.50 · 10 - 4  mm2 /s and 2.10 · 10 - 4  mm2 /s.

CONCLUSION

DW-SSFP can be used to investigate non-Gaussian diffusion by varying the flip angle. By fitting a model of non-Gaussian diffusion, the ADC in DW-SSFP can be estimated at an effective b-value, comparable to more conventional diffusion sequences.

Accelerated T2 Mapping of the Lumbar Intervertebral Disc: Highly Undersampled K-Space Data for Robust T2 Relaxation Time Measurement in Clinically Feasible Acquisition Times

T2 mapping of the intervertebral disc (IVD) can depict quantitative changes reflecting biochemical change due to loss of glycosaminoglycan content. Conventional T2 mapping is usually performed using a 2-dimensional multi-echo-spin echo sequence (2D-MESE) with long acquisition times that are generally not compatible with clinical routine. This study investigates the applicability of GRAPPATINI, a T2 mapping sequence combining undersampling, model-based reconstruction, and parallel imaging, to offer clinically feasible acquisition times in T2 mapping of the lumbar IVD.

MATERIALS AND METHODS

Fifty-eight individuals (26 female; mean age, 23.3 ± 8.1 years) were prospectively studied at 3 T. GRAPPATINI was conducted with the same parameters as the 2D-MESE while shortening the acquisition time from 13:18 to 2:27 minutes. The setup was also validated in a phantom experiment using a 6.48-hour-long single echo-spin echo sequence as reference. The IVDs were manually segmented on 4 central slices.

RESULTS

The median nucleus pulposus showed a strong Pearson correlation coefficient between T2GRAPPATINI and T2MESE (rp = 0.919; P < 0.001). There was also a significant correlation for the ventral (rp = 0.241; P < 0.001) and posterior (rp = 0.418; P < 0.001) annular regions.In the single spin-echo phantom experiment, the most accurate T2 estimation was achieved using T2GRAPPATINI with a median absolute deviation of 15.3 milliseconds as compared with T2MESE with 26.5 milliseconds.

CONCLUSIONS

GRAPPATINI facilitates precise T2 mapping at 3 T in accordance with clinical standards and reference methods using the same parameters while shortening acquisition times from 13:18 to 2:27 minutes with the same parameters.

Multi-parametric artificial neural network fitting of phase-cycled balanced steady-state free precession data

PURPOSE

Standard relaxation time quantification using phase-cycled balanced steady-state free precession (bSSFP), eg, motion-insensitive rapid configuration relaxometry (MIRACLE), is subject to a considerable underestimation of tissue T₁ and T₂ due to asymmetric intra-voxel frequency distributions. In this work, an artificial neural network (ANN) fitting approach is proposed to simultaneously extract accurate reference relaxation times (T₁ , T₂ ) and robust field map estimates ( B 1 + , ΔB₀ ) from the bSSFP profile.

METHODS

Whole-brain bSSFP data acquired at 3T were used for the training of a feedforward ANN with N = 12, 6, and 4 phase-cycles. The magnitude and phase of the Fourier transformed complex bSSFP frequency response served as input and the multi-parametric reference set [T₁ , T₂ , B 1 + , ∆B₀ ] as target. The ANN predicted relaxation times were validated against the target and MIRACLE.

RESULTS

The ANN prediction of T₁ and T₂ for trained and untrained data agreed well with the reference, even for only four acquired phase-cycles. In contrast, relaxometry based on 4-point MIRACLE was prone to severe off-resonance-related artifacts. ANN predicted B 1 + and ∆B₀ maps showed the expected spatial inhomogeneity patterns in high agreement with the reference measurements for 12-point, 6-point, and 4-point bSSFP phase-cycling schemes.

CONCLUSION

ANNs show promise to provide accurate brain tissue T₁ and T₂ values as well as reliable field map estimates. Moreover, the bSSFP acquisition can be accelerated by reducing the number of phase-cycles while still delivering robust T₁ , T₂ , B 1 + , and ∆B₀ estimates.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Unipolar MR elastography: Theory, numerical analysis and implementation

In MR elastography (MRE), zeroth moment balanced motion-encoding gradients (MEGs) are incorporated into MRI sequences to induce a phase shift proportional to the local displacement caused by external actuation. To maximize the signal-to-noise ratio (SNR), fractional encoding is employed, i.e., the MEG duration is reduced below the wave period. Here, gradients encode primarily the velocity of the motion-reducing encoding efficiency. Thus, in GRE-MRE, T₂ * decay and motion sensitivity have to be balanced, imposing a lower limit on repetition times (TRs). We propose to use a single trapezoidal gradient, a "unipolar gradient", to directly encode spin displacement. Such gradients cannot be used in conventional sequences as they exhibit a large zeroth moment and dephase magnetization. By time-reversing a spoiled SSFP sequence, the spoiling gradient becomes an efficient unipolar MEG. The proposed "unipolar MRE" technique benefits from this approach in three ways: first, displacement encoding is split over multiple TRs increasing motion sensitivity; second, spoiler and MEG coincide, allowing a reduction in TR; third, motion sensitivity of a typical unipolar lobe is of an order of magnitude higher than a bipolar MEG of equal duration. In this work, motion encoding using unipolar MRE is analyzed using the extended phase graph (EPG) formalism with a periodic motion propagator. As an approximation, the two-transverse TR approximation for diffusion-weighted SSFP is extended to incorporate cyclic motion. A complex encoding efficiency metric is introduced to compare the displacement fields of unipolar and conventional GRE-MRE sequences in both magnitude and phase. The derived theoretical encoding equations are used to characterize the proposed sequence using an extensive parameter study. Unipolar MRE is validated against conventional GRE-MRE in a phantom study showing excellent agreement between measured displacement fields. In addition, unipolar MRE yields significantly increased octahedral shear strain-SNR relative to conventional GRE-MRE and allows for the recovery of high stiffness inclusions, where conventional GRE-MRE fails.

Phase-based T2 mapping with gradient echo imaging

PURPOSE

Transverse relaxation time (T₂ ) mapping with MRI has a plethora of clinical and research applications. Current T₂ mapping techniques are based primarily on spin-echo (SE) relaxometry strategies that rely on the signal magnitude, and often suffer from lengthy acquisition times. In this work, we propose a phase-based T₂ mapping technique where T₂ information is encoded into the signal phase of rapid gradient echo (GRE) acquisitions.

THEORY

Bloch equation simulations demonstrate that the phase of GRE acquisitions obtained with a very small inter-repetition RF phase increment has a strong monotonic dependence on T₂ , resulting from coherent transverse magnetization. This T₂ -dependent phase behavior forms the basis of the proposed T₂ mapping technique. To isolate T₂ -dependent phase from background phase, at least 2 data sets with different RF phase increments are acquired. The proposed method can also be combined with chemical shift encoded MRI to separate water and fat signals.

METHODS

The feasibility of the proposed technique was validated in a phantom experiment. In vivo feasibility was demonstrated in the brain, knee, abdomen, and pelvis. Comparisons were made with SE-based T₂ mapping, spectroscopy, and T₂ values from the literature.

RESULTS

The proposed method produced accurate T₂ maps compared with SE-based T₂ mapping in the phantom. Good qualitative agreement was observed in vivo between the proposed method and the reference. T₂ measured in various anatomies agreed well with values reported in the literature.

CONCLUSION

A phase-based T₂ mapping technique was developed and its feasibility demonstrated in phantoms and in vivo.

Multi-pathway multi-echo acquisition and neural contrast translation to generate a variety of quantitative and qualitative image contrasts

PURPOSE

Clinical exams typically involve acquiring many different image contrasts to help discriminate healthy from diseased states. Ideally, 3D quantitative maps of all of the main MR parameters would be obtained for improved tissue characterization. Using data from a 7-min whole-brain multi-pathway multi-echo (MPME) scan, we aimed to synthesize several 3D quantitative maps (T₁ and T₂ ) and qualitative contrasts (MPRAGE, FLAIR, T₁ -weighted, T₂ -weighted, and proton density [PD]-weighted). The ability of MPME acquisitions to capture large amounts of information in a relatively short amount of time suggests it may help reduce the duration of neuro MR exams.

METHODS

Eight healthy volunteers were imaged at 3.0T using a 3D isotropic (1.2 mm) MPME sequence. Spin-echo, MPRAGE, and FLAIR scans were performed for training and validation. MPME signals were interpreted through neural networks for predictions of different quantitative and qualitative contrasts. Predictions were compared to reference values at voxel and region-of-interest levels.

RESULTS

Mean absolute errors (MAEs) for T₁ and T₂ maps were 216 ms and 11 ms, respectively. In ROIs containing white matter (WM) and thalamus tissues, the mean T₁ /T₂ predicted values were 899/62 ms and 1139/58 ms, consistent with reference values of 850/66 ms and 1126/58 ms, respectively. For qualitative contrasts, signals were normalized to those of WM, and MAEs for MPRAGE, FLAIR, T₁ -weighted, T₂ -weighted, and PD-weighted contrasts were 0.14, 0.15, 0.13, 0.16, and 0.05, respectively.

CONCLUSIONS

Using an MPME sequence and neural-network contrast translation, whole-brain results were obtained with a variety of quantitative and qualitative contrast in ~6.8 min.

Magnetic Resonance Imaging of the Musculoskeletal System at 7T: Morphological Imaging and Beyond

In 2017, a whole-body 7T magnetic resonance imaging (MRI) device was given regulatory approval for clinical use in both the EU and United States for neuro and musculoskeletal applications. As 7 Tesla allows for higher signal-to-noise , which results in higher resolution images than those obtained on lower-field-strength scanners, it has attracted considerable attention from the musculoskeletal field, as evidenced by the increasing number of publications in the last decade. Besides morphological imaging, the quantitative MR methods, such as T₂, T₂∗, T_1ρ mapping, sodium imaging, chemical-exchange saturation transfer, and spectroscopy, substantially benefit from ultrahigh field scanning. In this review, we provide technical considerations for the individual techniques and an overview of (mostly) clinical applications for the assessment of cartilage, tendon, meniscus, and muscle. The first part of the review is dedicated to morphological applications at 7T, and the second part describes the most recent developments in quantitative MRI at 7T.

Investigation of the influence of B0 drift on the performance of the PLANET method and an algorithm for drift correction

Purpose

The PLANET method was designed to simultaneously reconstruct maps of T₁ and T₂, the off‐resonance, the RF phase, and the banding free signal magnitude. The method requires a stationary B₀ field over the course of a phase‐cycled balanced SSFP acquisition. In this work we investigated the influence of B₀ drift on the performance of the PLANET method for single‐component and two‐component signal models, and we propose a strategy for drift correction.

Methods

The complex phase‐cycled balanced SSFP signal was modeled with and without frequency drift. The behavior of the signal influenced by drift was mathematically interpreted as a sum of drift‐dependent displacement of the data points along an ellipse and drift‐dependent rotation around the origin. The influence of drift on parameter estimates was investigated experimentally on a phantom and on the brain of healthy volunteers and was verified by numerical simulations. A drift correction algorithm was proposed and tested on a phantom and in vivo.

Results

Drift can be assumed to be linear over the typical duration of a PLANET acquisition. In a phantom (a single‐component signal model), drift induced errors of 4% and 8% in the estimated T₁ and T₂ values. In the brain, where multiple components are present, drift only had a minor effect. For both single‐component and two‐component signal models, drift‐induced errors were successfully corrected by applying the proposed drift correction algorithm.

Conclusion

We have demonstrated theoretically and experimentally the sensitivity of the PLANET method to B₀ drift and have proposed a drift correction method.

Accelerated 3D T2 mapping with dictionary-based matching for prostate imaging

PURPOSE

To develop a fast and accurate method for 3D T2 mapping of prostate cancer using undersampled acquisition and dictionary-based fitting.

METHODS

3D high-resolution T2 -weighted images (0.9 × 0.9 × 3 mm3 ) were obtained with a multishot T2 -prepared balanced steady-state free precession (T2 -prep-bSSFP) acquisition sequence using a 3D variable density undersampled Cartesian trajectory. Each T2 -weighted image was reconstructed using total variation regularized sensitivity encoding. A flexible simulation framework based on extended phase graphs generated a dictionary of magnetization signals, which was customized to the proposed sequence. The dictionary was matched to the acquired T2 -weighted images to retrieve quantitative T2 values, which were then compared to gold-standard spin echo acquisition values using monoexponential fitting. The proposed approach was validated in simulations and a T1 /T2 phantom, and feasibility was tested in 8 healthy subjects.

RESULTS

The simulation analysis showed that the proposed T2 mapping approach is robust to noise and typically observed T1 variations. T2 values obtained in the phantom with T2 prep-bSSFP and the acquisition-specific, dictionary-based matching were highly correlated with the gold-standard spin echo method (r = 0.99). Furthermore, no differences were observed with the accelerated acquisition compared to the fully sampled acquisition (r = 0.99). T2 values obtained in prostate peripheral zone, central gland, and muscle in healthy subjects (age, 26 ± 6 years) were 97 ± 14, 76 ± 7, and 36 ± 3 ms, respectively.

CONCLUSION

3D quantitative T2 mapping of the whole prostate can be achieved in 3 minutes.

Frequency-modulated SSFP with radial sampling and subspace reconstruction: A time-efficient alternative to phase-cycled bSSFP

PURPOSE

A novel subspace-based reconstruction method for frequency-modulated balanced steady-state free precession (fmSSFP) MRI is presented. In this work, suitable data acquisition schemes, subspace sizes, and efficiencies for banding removal are investigated.

THEORY AND METHODS

By combining a fmSSFP MRI sequence with a 3D stack-of-stars trajectory, scan efficiency is maximized as spectral information is obtained without intermediate preparation phases. A memory-efficient reconstruction routine is implemented by introducing the low-frequency Fourier transform as a subspace which allows for the formulation of a convex reconstruction problem. The removal of banding artifacts is investigated by comparing the proposed acquisition and reconstruction technique to phase-cycled bSSFP MRI. Aliasing properties of different undersampling schemes are analyzed and water/fat separation is demonstrated by reweighting the reconstructed subspace coefficients to generate virtual spectral responses in a post-processing step.

RESULTS

A simple root-of-sum-of-squares combination of the reconstructed subspace coefficients yields high-SNR images with the characteristic bSSFP contrast but without banding artifacts. Compared to Golden-Angle trajectories, turn-based sampling schemes were superior in minimizing aliasing across reconstructed subspace coefficients. Water/fat separated images of the human knee were obtained by reweighting subspace coefficients.

CONCLUSIONS

The novel subspace-based fmSSFP MRI technique emerges as a time-efficient alternative to phase-cycled bSFFP. The method does not need intermediate preparation phases, offers high SNR and avoids banding artifacts. Reweighting of the reconstructed subspace coefficients allows for generating virtual spectral responses with applications to water/fat separation.

bSSFP phase correction and its use in magnetic resonance electrical properties tomography

PURPOSE

Balanced steady-state free precession (bSSFP) sequence is widely used because of its high SNR and high speed. However, bSSFP images suffer from "banding artifact" caused by B₀ inhomogeneity. In this article, we propose a method to remove this artifact in bSSFP phase images and investigate the usage of the corrected phase images in phase-based magnetic resonance electrical properties tomography (MREPT).

THEORY AND METHODS

Two bSSFP phase images, obtained with different excitation frequencies, are collaged to get rid of the regions containing banding artifacts. Phase of the collaged bSSFP image is the sum of the transceive phase of the RF system and an error term that depends on B₀ and T₂ . By using B₀ and T₂ maps, this error is eliminated from bSSFP phase images by using pixel-wise corrections. Conductivity maps are obtained from the uncorrected and the corrected phase images using the phase-based cr-MREPT method.

RESULTS

Phantom and human experiment results of the proposed method are illustrated for both phase images and conductivity maps. It is shown that uncorrected phase images yield unacceptable conductivity images. When only B₀ information is used for phase correction conductivity, reconstructions are substantially improved, and yet T₂ information is still needed to fully recover accurate and undistorted conductivity images.

CONCLUSIONS

With the proposed technique, B₀ sensitivity of the bSSFP phase images can be removed by using B₀ and T₂ maps. It is also shown that corrected bSSFP phase images are of sufficient quality to be used in conductivity imaging.

Multipathway multi-echo (MPME) imaging: all main MR parameters mapped based on a single 3D scan

PURPOSE

Quantitative parameter maps, as opposed to qualitative grayscale images, may represent the future of diagnostic MRI. A new quantitative MRI method is introduced here that requires a single 3D acquisition, allowing good spatial coverage to be achieved in relatively short scan times.

METHODS

A multipathway multi-echo sequence was developed, and at least 3 pathways with 2 TEs were needed to generate T₁ , T₂ , T₂ ^* , B₁ ⁺ , and B₀ maps. The method required the central k-space region to be sampled twice, with the same sequence but with 2 very different nominal flip angle settings. Consequently, scan time was only slightly longer than that of a single scan. The multipathway multi-echo data were reconstructed into parameter maps, for phantom as well as brain acquisitions, in 5 healthy volunteers at 3 T. Spatial resolution, matrix size, and FOV were 1.2 × 1.0 × 1.2 mm³ , 160 × 192 × 160, and 19.2 × 19.2 × 19.2 cm³ (whole brain), acquired in 11.5 minutes with minimal acceleration. Validation was performed against T₁ , T₂ , and T₂ ^* maps calculated from gradient-echo and spin-echo data.

RESULTS

In Bland-Altman plots, bias and limits of agreement for T₁ and T₂ results in vivo and in phantom were -2.9/±125.5 ms (T₁ in vivo), -4.8/±20.8 ms (T₂ in vivo), -1.5/±18.1 ms (T₁ in phantom), and -5.3/±7.4 ms (T₂ in phantom), for regions of interest including given brain structures or phantom compartments. Due to relatively high noise levels, the current implementation of the approach may prove more useful for region of interest-based as opposed to pixel-based interpretation.

CONCLUSIONS

We proposed a novel approach to quantitatively map MR parameters based on a multipathway multi-echo acquisition.

Simultaneous fat-free isotropic 3D anatomical imaging and T2 mapping of knee cartilage with lipid-insensitive binomial off-resonant RF excitation (LIBRE) pulses

BACKGROUND

Improved knee cartilage morphological delineation and T₂ mapping precision necessitates isotropic 3D high-resolution and efficient fat suppression.

PURPOSE

To develop and assess an isotropic 3D lipid-insensitive T₂ mapping technique of the knee for improved cartilage delineation and precise measurement of T₂ relaxation times.

STUDY TYPE

Prospective.

PHANTOM/SUBJECTS

Phantoms (n = 6) used in this study were designed to mimic the T₁ and T₂ relaxation times of cartilage and fat. The study cohort comprised healthy volunteers (n = 7) for morphometry and T₂ relaxation time measurements.

FIELD STRENGTH/SEQUENCE

A high-resolution isotropic 3D T₂ mapping technique that uses sequential T₂ -prepared segmented gradient-recalled echo (Iso3DGRE) images and lipid-insensitive binomial off-resonant radiofrequency (RF) excitation (LIBRE) at 3T.

ASSESSMENT

Numerical simulations and phantom experiments were performed to optimize the LIBRE pulse. Phantom studies were carried out to test the accuracy of the technique against reference standard spin-echo (SE) T₂ mapping. Subsequently, T₂ maps with and without LIBRE pulses were acquired in knees of healthy volunteers and the T₂ relaxation time values in different cartilage compartments were compared.

STATISTICAL TESTS

A two-tailed paired Student's t-test was used to compare the average T₂ values and the relative standard deviations (inverse measurement of the precision) obtained with and without LIBRE pulses.

RESULTS

A LIBRE pulse of 1 msec suppressed fat with an RF excitation frequency offset of 1560 Hz and optimal RF excitation angle of 35°. These results were corroborated by phantom and knee experiments. Robust and homogeneous fat suppression was obtained (a fat signal-to-noise ratio (SNR) decrease of 86.4 ± 2.4%). In phantoms, T₂ values were found in good agreement when comparing LIBRE-Iso3DGRE with SE (slope 0.93 ± 0.04, intercept 0.11 ± 1.6 msec, R² >0.99). In vivo, LIBRE excitation resulted in more precise T₂ estimation (23.7 ± 7.4%) than normal excitation (30.5 ± 9.9%, P < 0.0001).

DATA CONCLUSION

Homogeneous LIBRE fat signal suppression was achieved with a total RF pulse duration of 1 msec, allowing for the removal of chemical shift artifacts and resulting in improved cartilage delineation and precise T₂ values.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:1275-1284.

On the accuracy and precision of PLANET for multiparametric MRI using phase-cycled bSSFP imaging

Purpose

In this work we demonstrate how sequence parameter settings influence the accuracy and precision in T₁, T₂, and off‐resonance maps obtained with the PLANET method for a single‐component signal model. In addition, the performance of the method for the particular case of a two‐component relaxation model for white matter tissue was assessed.

Methods

Numerical simulations were performed to investigate the influence of sequence parameter settings on the accuracy and precision in the estimated parameters for a single‐component model, as well as for a two‐component white matter model. Phantom and in vivo experiments were performed for validation. In addition, the effects of Gibbs ringing were investigated.

Results

By making a proper choice for sequence parameter settings, accurate and precise parameter estimation can be achieved for a single‐component signal model over a wide range of relaxation times at realistic SNR levels. Due to the presence of a second myelin‐related signal component in white matter, an underestimation of approximately 30% in T₁ and T₂ was observed, predicted by simulations and confirmed by measurements. Gibbs ringing artifacts correction improved the precision and accuracy of the parameter estimates.

Conclusion

For a single‐component signal model there is a broad “sweet spot” of sequence parameter combinations for which a high accuracy and precision in the parameter estimates is achieved over a wide range of relaxation times. For a multicomponent signal model, the single‐component PLANET reconstruction results in systematic errors in the parameter estimates as expected.

The comparison of the performance of 3 T and 7 T T2 mapping for untreated low-grade cartilage lesions

OBJECTIVE

To investigate T₂ mapping as a possible marker for low-grade human articular cartilage lesions during a one-year follow-up, possible changes during the follow-up and compare the reliability and sensitivity of these measurements on high-field (3 T) and ultra-high-field (7 T) MRI scanners.

DESIGN

Twenty-one patients with femoral, tibial and patellar cartilage defect in the knee joint participated in the study. The MRI protocol consisted of morphological, as well as three-dimensional triple-echo steady-state (3D-TESS) T₂ mapping sequences with similar parameters at 3T and 7T. Patients were scanned at five time-points up to 12 months. T₂ values were evaluated in the lesion and healthy-appearing regions for superficial and deep cartilage zone. The repeated ANOVA was used to determine differences in T₂ values at various time points.

RESULTS

A significant decrease in T₂ values was observed between baseline and six months in the superficial layer of the lesion in patients at 3 T (decrease from 41.89 ± 9.3 ms to 31.21 ± 7.2 ms, which is a difference of -5.67 ± 2.2 ms (p = 0.031)), and at 12 months in the superficial layer of the lesion in patients at 3 T (decrease from 41.89 ± 9.3 ms to 35.28 ± 4.9 ms, which is a difference of -6.60 ± 4.4 ms (p = 0.044). No significant differences were recorded at 7 T.

CONCLUSION

The change in T₂ values acquired with 3 T 3D-TESS appears to be reflecting subtle changes of cartilage composition in the course of low-grade lesion development. 7 T T₂ mapping does not reflect these changes probably due to completely decayed short T₂ component.

Fast quantitative MRI using controlled saturation magnetization transfer

Purpose

This study demonstrates magnetization transfer (MT) effects directly affect relaxometry measurements and develops a framework that allows single‐pool models to be valid in 2‐pool MT systems.

Methods

A theoretical framework is developed in which a 2‐pool MT system effectively behaves as a single‐pool if the RMS RF magnetic field (B1rms{\text{B}}_{1}^{{{\text{rms}}}}) is kept fixed across all measurements. A practical method for achieving controlled saturation magnetization transfer (CSMT) using multiband RF pulses is proposed. Numerical, Phantom, and in vivo validations were performed directly comparing steady state (SS) estimation approaches that under correct single‐pool assumptions would be expected to vary in precision but not accuracy.

Results

Numerical simulations predict single‐pool estimates obtained from MT model generated data are not consistent for different SS estimation methods, and a systematic underestimation of T₂ is expected. Neither effect occurs under the proposed CSMT approach. Both phantom and in vivo experiments corroborate the numerical predictions. Experimental data highlights that even when using the same relaxometry method, different estimates are obtained depending on which combination of flip angles (FAs) and TRs are used if the CSMT approach is not used. Using CSMT, stable measurements of both T₁ and T₂ are obtained. The measured T₁(T1CSMT)) depends on B1rms{\text{B}}_{1}^{{{\text{rms}}}}, which is therefore an important parameter to specify.

Conclusion

This work demonstrates that conventional single pool relaxometry, which is highly efficient for human studies, results in unreliable parameter estimates in biological tissues because of MT effects. The proposed CSMT framework is shown to allow single‐pool assumptions to be valid, enabling reliable and efficient quantitative imaging to be performed.

Simultaneous multi-parametric mapping of total sodium concentration, T1, T2 and ADC at 7 T using a multi-contrast unbalanced SSFP

PURPOSE

Quantifying multiple NMR properties of sodium could be of benefit to assess changes in cellular viability in biological tissues. A proof of concept of Quantitative Imaging using Configuration States (QuICS) based on a SSFP sequence with multiple contrasts was implemented to extract simultaneously 3D maps of applied flip angle (FA), total sodium concentration, T₁, T₂, and Apparent Diffusion Coefficient (ADC).

METHODS

A 3D Cartesian Gradient Recalled Echo (GRE) sequence was used to acquire 11 non-balanced SSFP contrasts at a 6 × 6 × 6 mm³ isotropic resolution with carefully-chosen gradient spoiling area, RF amplitude and phase cycling, with TR/TE = 20/3.2 ms and 25 averages, leading to a total acquisition time of 1 h 18 min. A least-squares fit between the measured and the analytical complex signals was performed to extract quantitative maps from a mono-exponential model. Multiple sodium phantoms with different compositions were studied to validate the ability of the method to measure sodium NMR properties in various conditions.

RESULTS

Flip angle maps were retrieved. Relaxation times, ADC and sodium concentrations were estimated with controlled precision below 15%, and were in accordance with measurements from established methods and literature.

CONCLUSION

The results illustrate the ability to retrieve sodium NMR properties maps, which is a first step toward the estimation of FA, T₁, T₂, concentration and ADC of ²³Na for clinical research. With further optimization of the acquired QuICS contrasts, scan time could be reduced to be suitable with in vivo applications.

Diffusion effect on T2 relaxometry in triple-echo steady state free precession sequence

PURPOSE

The purpose of this study is to evaluate the effect of diffusion on SSFP (Steady-state Free Precession) signals in triple-echo steady state (TESS) sequence and ultimately on the accuracy of T2 relaxometry.

METHODS

The extended phase graph (EPG) algorithm was used to study the effect of diffusion on SSFP signals and T2 relaxometry. The simulation results were verified by a commercial phantom and in vivo studies. Based on the simulation results, a correction scheme was proposed to correct the estimated T2 values.

RESULTS

T2 underestimation in TESS was evident in case of small flip angle and large unbalanced gradient moment on objects with large T2 and D values. The T2 underestimation mainly originated from the diffusion sensitivity of SSFP-echo. It was also observed that SSFP-FID (Free Induction Decay) signals increased with increasing diffusion weighting under some specific conditions. The proposed correction scheme corrected the T2 underestimation, which verified that the underestimation was due to the neglect of diffusion effect. For clinical practice of TESS in tissues with short T2 such as cartilage and muscle, the diffusion effect of TESS is negligible.

CONCLUSION

The effect of diffusion cannot be neglected during TESS T2 quantification as it is the main source of T2 underestimation when small flip angle and large unbalanced gradient moment is used, especially for objects with large T2 and D values.

Skeletal Muscle Quantitative Nuclear Magnetic Resonance Imaging and Spectroscopy as an Outcome Measure for Clinical Trials

Recent years have seen tremendous progress towards therapy of many previously incurable neuromuscular diseases. This new context has acted as a driving force for the development of novel non-invasive outcome measures. These can be organized in three main categories: functional tools, fluid biomarkers and imagery. In the latest category, nuclear magnetic resonance imaging (NMRI) offers a considerable range of possibilities for the characterization of skeletal muscle composition, function and metabolism. Nowadays, three NMR outcome measures are frequently integrated in clinical research protocols. They are: 1/ the muscle cross sectional area or volume, 2/ the percentage of intramuscular fat and 3/ the muscle water T2, which quantity muscle trophicity, chronic fatty degenerative changes and oedema (or more broadly, “disease activity”), respectively. A fourth biomarker, the contractile tissue volume is easily derived from the first two ones. The fat fraction maps most often acquired with Dixon sequences have proven their capability to detect small changes in muscle composition and have repeatedly shown superior sensitivity over standard functional evaluation. This outcome measure will more than likely be the first of the series to be validated as an endpoint by regulatory agencies. The versatility of contrast generated by NMR has opened many additional possibilities for characterization of the skeletal muscle and will result in the proposal of more NMR biomarkers. Ultra-short TE (UTE) sequences, late gadolinium enhancement and NMR elastography are being investigated as candidates to evaluate skeletal muscle interstitial fibrosis. Many options exist to measure muscle perfusion and oxygenation by NMR. Diffusion NMR as well as texture analysis algorithms could generate complementary information on muscle organization at microscopic and mesoscopic scales, respectively. ³¹P NMR spectroscopy is the reference technique to assess muscle energetics non-invasively during and after exercise. In dystrophic muscle, ³¹P NMR spectrum at rest is profoundly perturbed, and several resonances inform on cell membrane integrity. Considerable efforts are being directed towards acceleration of image acquisitions using a variety of approaches, from the extraction of fat content and water T2 maps from one single acquisition to partial matrices acquisition schemes. Spectacular decreases in examination time are expected in the near future. They will reinforce the attractiveness of NMR outcome measures and will further facilitate their integration in clinical research trials.

[When is cartilage repair successful?]

Focal cartilage lesions are a cause of long-term disability and morbidity. After cartilage repair, it is crucial to evaluate long-term progression or failure in a reproducible, standardized manner. This article provides an overview of the different cartilage repair procedures and important characteristics to look for in cartilage repair imaging. Specifics and pitfalls are pointed out alongside general aspects. After successful cartilage repair, a complete, but not hypertrophic filling of the defect is the primary criterion of treatment success. The repair tissue should also be completely integrated to the surrounding native cartilage. After some months, the transplants signal should be isointense compared to native cartilage. Complications like osteophytes, subchondral defects, cysts, adhesion and chronic bone marrow edema or joint effusion are common and have to be observed via follow-up. Radiological evaluation and interpretation of postoperative changes should always take the repair method into account.

PLANET: An ellipse fitting approach for simultaneous T1 and T2 mapping using phase-cycled balanced steady-state free precession

Purpose

To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T₁ and T₂ from a single 3D phase‐cycled balanced steady‐state free precession (bSSFP) sequence.

Methods

A method is presented in which the elliptical signal model is used to describe the phase‐cycled bSSFP steady‐state signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T₁ and T₂, the banding free magnitude, and the off‐resonance are calculated from the fitting results.

Results

Maps of T₁ and T₂, as well as an off‐resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase‐cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T₁ was observed in the brain.

Conclusion

The presented method allows for accurate mapping of relaxation times and off‐resonance, and for the reconstruction of banding free magnitude images at realistic signal‐to‐noise ratios. Magn Reson Med 79:711–722, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.