Asymmetries of the balanced SSFP profile. Part I: theory and observation

Getting the phase consistent: The importance of phase description in balanced steady-state free precession MRI of multi-compartment systems

PURPOSE

Determine the correct mathematical phase description for balanced steady-state free precession (bSSFP) signals in multi-compartment systems.

THEORY AND METHODS

Based on published bSSFP signal models, different phase descriptions can be formulated: one predicting the presence and the other predicting the absence of destructive interference effects in multi-compartment systems. Numerical simulations of bSSFP signals of water and acetone were performed to evaluate the predictions of these different phase descriptions. For experimental validation, bSSFP profiles were measured at 3T using phase-cycled bSSFP acquisitions performed in a phantom containing mixtures of water and acetone, which replicates a system with two signal components. Localized single voxel MRS was performed at 7T to determine the relative chemical shift of the acetone-water mixtures.

RESULTS

Based on the choice of phase description, the simulated bSSFP profiles of water-acetone mixtures varied significantly, either displaying or lacking destructive interference effects, as predicted theoretically. In phantom experiments, destructive interference was consistently observed in the measured bSSFP profiles of water-acetone mixtures, supporting the theoretical description that predicts such interference effects. The connection between the choice of phase description and predicted observation enables unambiguous experimental identification of the correct phase description for multi-compartment bSSFP profiles, which is consistent with the Bloch equations.

CONCLUSION

The study emphasizes that consistent phase descriptions are crucial for accurately describing multi-compartment bSSFP signals, as incorrect phase descriptions result in erroneous predictions.

Myelin water imaging at 0.55 T using a multigradient-echo sequence

PURPOSE

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T.

METHODS

Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T.

RESULTS

The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate.

CONCLUSION

This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.

SPARCQ: A new approach for fat fraction mapping using asymmetries in the phase-cycled balanced SSFP signal profile

PURPOSE

To develop SPARCQ (Signal Profile Asymmetries for Rapid Compartment Quantification), a novel approach to quantify fat fraction (FF) using asymmetries in the phase-cycled balanced SSFP (bSSFP) profile.

METHODS

SPARCQ uses phase-cycling to obtain bSSFP frequency profiles, which display asymmetries in the presence of fat and water at certain TRs. For each voxel, the measured signal profile is decomposed into a weighted sum of simulated profiles via multi-compartment dictionary matching. Each dictionary entry represents a single-compartment bSSFP profile with a specific off-resonance frequency and relaxation time ratio. Using the results of dictionary matching, the fractions of the different off-resonance components are extracted for each voxel, generating quantitative maps of water and FF and banding-artifact-free images for the entire image volume. SPARCQ was validated using simulations, experiments in a water-fat phantom and in knees of healthy volunteers. Experimental results were compared with reference proton density FFs obtained with 1 H-MRS (phantoms) and with multiecho gradient-echo MRI (phantoms and volunteers). SPARCQ repeatability was evaluated in six scan-rescan experiments.

RESULTS

Simulations showed that FF quantification is accurate and robust for SNRs greater than 20. Phantom experiments demonstrated good agreement between SPARCQ and gold standard FFs. In volunteers, banding-artifact-free quantitative maps and water-fat-separated images obtained with SPARCQ and ME-GRE demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ FF was 0.0512.

CONCLUSION

SPARCQ demonstrates potential for fat quantification using asymmetries in bSSFP profiles and may be a promising alternative to conventional FF quantification techniques.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination ϴ) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination ϴ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.

Constrained Ellipse Fitting for Efficient Parameter Mapping With Phase-Cycled bSSFP MRI

Balanced steady-state free precession (bSSFP) imaging enables high scan efficiency in MRI, but differs from conventional sequences in terms of elevated sensitivity to main field inhomogeneity and nonstandard [Formula: see text]-weighted tissue contrast. To address these limitations, multiple bSSFP images of the same anatomy are commonly acquired with a set of different RF phase-cycling increments. Joint processing of phase-cycled acquisitions serves to mitigate sensitivity to field inhomogeneity. Recently phase-cycled bSSFP acquisitions were also leveraged to estimate relaxation parameters based on explicit signal models. While effective, these model-based methods often involve a large number of acquisitions (N ≈ 10-16), degrading scan efficiency. Here, we propose a new constrained ellipse fitting method (CELF) for parameter estimation with improved efficiency and accuracy in phase-cycled bSSFP MRI. CELF is based on the elliptical signal model framework for complex bSSFP signals; and it introduces geometrical constraints on ellipse properties to improve estimation efficiency, and dictionary-based identification to improve estimation accuracy. CELF generates maps of [Formula: see text], [Formula: see text], off-resonance and on-resonant bSSFP signal by employing a separate [Formula: see text] map to mitigate sensitivity to flip angle variations. Our results indicate that CELF can produce accurate off-resonance and banding-free bSSFP maps with as few as N = 4 acquisitions, while estimation accuracy for relaxation parameters is notably limited by biases from microstructural sensitivity of bSSFP imaging.

Pure balanced steady-state free precession imaging (pure bSSFP)

Purpose

To show that for tissues the conspicuous asymmetries in the frequency response function of bSSFP can be mitigated by using a short enough TR.

Theory and Methods

Configuration theory indicates that bSSFP becomes apparently “pure” (i.e., exhibiting a symmetric profile) in the limit of TR →0. To this end, the frequency profile of bSSFP was measured as a function of the TR using a manganese‐doped aqueous probe, as well as brain tissue that was shown to exhibit a pronounced asymmetry due to its microstructure. The frequency response function was sampled using N=72 (phantom) and N=36 (in vivo) equally distributed linear RF phase increments in the interval [0,2π). Imaging was performed with 2.0 mm isotropic resolution over a TR range of 1.5–8 ms at 3 and 1.5 T.

Results

As expected, pure substances showed a symmetric TR‐independent frequency profile, whereas brain tissue revealed a pronounced asymmetry. The observed asymmetry for the tissue, however, decreases with decreasing TR and gives strong evidence that the frequency response function of bSSFP becomes symmetric in the limit of TR →0, in agreement with theory. The limit of apparently pure bSSFP imaging can thus be achieved for a TR ∼ 1.5 ms at 1.5 T, whereas at 3 T, tissues still show some residual asymmetry.

Conclusion

In the limit of short enough TR, tissues become apparently pure for bSSFP. This limit can be reached for brain tissue at 1.5 T with TR ∼ 1–2 ms at clinically relevant resolutions.

From signal-based to comprehensive magnetic resonance imaging

We present and evaluate a new insight into magnetic resonance imaging (MRI). It is based on the algebraic description of the magnetization during the transient response—including intrinsic magnetic resonance parameters such as longitudinal and transverse relaxation times (T₁, T₂) and proton density (PD) and experimental conditions such as radiofrequency field (B₁) and constant/homogeneous magnetic field (B₀) from associated scanners. We exploit the correspondence among three different elements: the signal evolution as a result of a repetitive sequence of blocks of radiofrequency excitation pulses and encoding gradients, the continuous Bloch equations and the mathematical description of a sequence as a linear system. This approach simultaneously provides, in a single measurement, all quantitative parameters of interest as well as associated system imperfections. Finally, we demonstrate the in-vivo applicability of the new concept on a clinical MRI scanner.

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.

Sensitivity to myelin using model-free analysis of the water resonance line-shape in postmortem mouse brain

PURPOSE

Dysmyelinating diseases are characterized by abnormal myelin formation and function. Such microstructural abnormalities in myelin have been demonstrated to produce measurable effects on the MR signal. This work examines these effects on measurements of voxel-wise, high-resolution water spectra acquired using a 3D echo-planar spectroscopic imaging (EPSI) pulse sequence from both postmortem fixed control mouse brains and a dysmyelination mouse brain model.

METHODS

Perfusion fixed, resected control (n = 5) and shiverer (n = 4) mouse brains were imaged using 3D-EPSI with 100 µm isotropic resolution. The free induction decay (FID) was sampled every 2.74 ms over 192 echoes, for a total sampling duration of 526.08 ms. Voxel-wise FIDs were Fourier transformed to produce water spectra with 1.9 Hz resolution. Spectral asymmetry was computed and compared between the two tissue types.

RESULTS

The water resonance is more asymmetrically broadened in the white matter of control mouse brain compared with dysmyelinated white matter. In control brain, this is modulated by and consistent with previously reported orientationally dependent effects of white matter relative to B₀ . Similar sensitivity to orientation is observed in dysmyelinated white matter as well; however, the magnitude of the resonance asymmetry is much lower across all directions.

CONCLUSION

Results demonstrate that components of the spectra are specifically differentially affected by myelin concentration. This suggests that water proton spectra may be sensitive to the presence of myelin, and as such, could serve as a MRI-based biomarker of dysmyelinating disease, free of mathematical models.

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.

Structure or Exchange? On the Feasibility of Chemical Exchange Detection with Balanced Steady-State Free Precession in Tissue - An In Vitro Study

Balanced steady-state free precession imaging has recently been suggested for chemical exchange detection (bSSFPX). The objective of this work is to investigate the contributions of microstructural, chemical shift and chemical exchange effects to the asymmetry of the bSSFP profile at field strengths of 3 T and 9.4 T. To this end, in vitro bSSFPX experiments are performed for a range of repetition times and flip angles in glucose water solutions with different MnCl₂ concentrations and tissue homogenates obtained from the brainstem of pig brains. The experimental results are compared to multi-pool Bloch-McConnell simulations. Additionally, the influence of white matter tract geometry is analyzed ex vivo in pig brain hemispheres measured at two different angles with respect to B₀ . The detectable bSSFP profile asymmetry in glucose solutions with tissue-like relaxation times and white matter homogenates was consistent with Bloch-McConnell simulations but relatively small. In intact white matter tracts, the asymmetry was dominated by structural effects with a strong dependency on tract orientation relative to B₀ . In tracts perpendicular to B₀ , the asymmetry was ≈ 3-4 times higher than in the homogenates, thus barely affected by chemical exchange effects. In conclusion, chemical exchange-related bSSFP profile asymmetries are detectable in tissue homogenates, however, the observed asymmetry level is generally low and prone to confounding structural effects rendering in vivo chemical exchange detection with bSSFP challenging in the brain.

Magnetization transfer and frequency distribution effects in the SSFP ellipse

Purpose

To demonstrate that quantitative magnetization transfer (qMT) parameters can be extracted from steady‐state free‐precession (SSFP) data with no external T₁ map or banding artifacts.

Methods

SSFP images with multiple MT weightings were acquired and qMT parameters fitted with a two‐stage elliptical signal model.

Results

Monte Carlo simulations and data from a 3T scanner indicated that most qMT parameters could be recovered with reasonable accuracy. Systematic deviations from theory were observed in white matter, consistent with previous literature on frequency distribution effects.

Conclusions

qMT parameters can be extracted from SSFP data alone, in a manner robust to banding artifacts, despite several confounds.

Transceive phase mapping using the PLANET method and its application for conductivity mapping in the brain

Purpose

To demonstrate feasibility of transceive phase mapping with the PLANET method and its application for conductivity reconstruction in the brain.

Methods

Accuracy and precision of transceive phase (ϕ^±) estimation with PLANET, an ellipse fitting approach to phase‐cycled balanced steady state free precession (bSSFP) data, were assessed with simulations and measurements and compared to standard bSSFP. Measurements were conducted on a homogeneous phantom and in the brain of healthy volunteers at 3 tesla. Conductivity maps were reconstructed with Helmholtz‐based electrical properties tomography. In measurements, PLANET was also compared to a reference technique for transceive phase mapping, i.e., spin echo.

Results

Accuracy and precision of ϕ^± estimated with PLANET depended on the chosen flip angle and TR. PLANET‐based ϕ^± was less sensitive to perturbations induced by off‐resonance effects and partial volume (e.g., white matter + myelin) than bSSFP‐based ϕ^±. For flip angle = 25° and TR = 4.6 ms, PLANET showed an accuracy comparable to that of reference spin echo but a higher precision than bSSFP and spin echo (factor of 2 and 3, respectively). The acquisition time for PLANET was ~5 min; 2 min faster than spin echo and 8 times slower than bSSFP. However, PLANET simultaneously reconstructed T₁, T₂, B₀ maps besides mapping ϕ^±. In the phantom, PLANET‐based conductivity matched the true value and had the smallest spread of the three methods. In vivo, PLANET‐based conductivity was similar to spin echo‐based conductivity.

Conclusion

Provided that appropriate sequence parameters are used, PLANET delivers accurate and precise ϕ^± maps, which can be used to reconstruct brain tissue conductivity while simultaneously recovering T₁, T₂, and B₀ maps.

Multi-slice passband bSSFP for human and rodent fMRI at ultra-high field

Balanced steady-state free precession (bSSFP) can be used as an alternative to gradient-echo (GE) EPI for BOLD functional MRI when image distortions and signal drop-outs are severe such as at ultra-high field. However, 3D-bSSFP acquisitions have distinct drawbacks on either human or animal MR systems. On clinical scanners, 3D imaging is suboptimal for localized fMRI applications. It can also display distortions when acceleration methods such as spiral read-outs are used, and, compared to multi-slice acquisitions, suffers from increased sensitivity to motion or physiological noise which further results in blurring. On pre-clinical systems, 3D acquisitions have low temporal resolution due to limited acceleration options, while single slice often results in insufficient coverage. The aim of the present study was to implement a multi-slice bSSFP acquisition with Cartesian read-out to obtain non-distorted BOLD fMRI activation maps in the human and rat brain at ultra-high field. We show that, when using a new pseudo-steady-state, the bSSFP signal characteristics are preserved. In the human brain at 7 T, we demonstrate that both task- and resting-state fMRI can be performed with multi-slice bSSFP, with a temporal SNR that matches that of 3D-bSSFP, resulting in - at least - equal performance. In the rat brain at 14 T, we show that the multi-slice bSSFP protocol has similar sensitivity to gradient-echo EPI for task fMRI, while benefitting from much reduced distortions and drop-outs. The advantages of passband bSSFP at 14 T in comparison with GE-EPI are expected to be even more marked for mouse brain.

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.

Optimized quantification of spin relaxation times in the hybrid state

PURPOSE

The optimization and analysis of spin ensemble trajectories in the hybrid state-a state in which the direction of the magnetization adiabatically follows the steady state while the magnitude remains in a transient state.

METHODS

Numerical optimizations were performed to find spin ensemble trajectories that minimize the Cramér-Rao bound for T 1 -encoding, T 2 -encoding, and their weighted sum, respectively, followed by a comparison between the Cramér-Rao bounds obtained with our optimized spin-trajectories, Look-Locker sequences, and multi-spin-echo methods. Finally, we experimentally tested our optimized spin trajectories with in vivo scans of the human brain.

RESULTS

After a nonrecurring inversion segment on the southern half of the Bloch sphere, all optimized spin trajectories pursue repetitive loops on the northern hemisphere in which the beginning of the first and the end of the last loop deviate from the others. The numerical results obtained in this work align well with intuitive insights gleaned directly from the governing equation. Our results suggest that hybrid-state sequences outperform traditional methods. Moreover, hybrid-state sequences that balance T 1 - and T 2 -encoding still result in near optimal signal-to-noise efficiency for each relaxation time. Thus, the second parameter can be encoded at virtually no extra cost.

CONCLUSIONS

We provided new insights into the optimal encoding processes of spin relaxation times in order to guide the design of robust and efficient pulse sequences. We found that joint acquisitions of T 1 and T 2 in the hybrid state are substantially more efficient than sequential encoding techniques.

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.

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.

Balanced Steady-State Free Precession (bSSFP) from an effective field perspective: Application to the detection of chemical exchange (bSSFPX)

Chemical exchange saturation transfer (CEST) is a novel contrast mechanism and it is gaining increasing popularity as many promising applications have been proposed and investigated. Fast and quantitative CEST imaging techniques are further needed in order to increase the applicability of CEST for clinical use as well as to derive quantitative physiological and biological information. Steady-state methods for fast CEST imaging have been reported recently. Here, we observe that an extreme case of these methods is a balanced steady-state free precession (bSSFP) sequence. The bSSFP in itself is sensitive to the exchange processes; hence, no additional saturation or preparation is needed for CEST-like data acquisition. The bSSFP experiment can be regarded as observation during saturation, without separate saturation and acquisition modules as used in standard CEST and similar experiments. One of the differences from standard CEST methods is that the bSSFP spectrum is an XY-spectrum not a Z-spectrum. As the first proof-of-principle step, we have implemented the steady-state bSSFP sequence for chemical exchange detection (bSSFPX) and verified its feasibility in phantom studies. These studies have shown that bSSFPX can achieve exchange-mediated contrast comparable to the standard CEST experiment. Therefore, the bSSFPX method has a potential for fast and quantitative CEST data acquisition.

Characterization of T1 bias in skeletal muscle from fat in MOLLI and SASHA pulse sequences: Quantitative fat-fraction imaging with T1 mapping

PURPOSE

To characterize the effects of fat on commonly used T₁ mapping sequences and evaluate a new method of quantitative fat fraction (FF) imaging for low fractions based on the modulation of T₁ values by the fat pool.

METHODS

Bloch equation simulations and phantom and in vivo (skeletal muscle) experiments were used to characterize the response of the modified Look-Locker inversion recovery (MOLLI) and saturation recovery single-shot acquisition (SASHA) T₁ mapping sequences to fat-water systems with known FFs (0%-10%) at 1.5T. FFs were measured with single voxel spectroscopy and Dixon imaging methods. A new T₁ -based FF imaging method was evaluated using Monte Carlo simulations and phantom and in vivo experiments.

RESULTS

SASHA and MOLLI had similar T₁ dependence on FF, with characteristic under- or overestimation of T₁ values as a function of off-resonance frequency (30-70 ms variation in native T₁ per 1% FF). FF maps generated from the SASHA method yielded a low variability of ±0.25% for a signal-to-noise ratio of 150:1 in the nonsaturation image, with good agreement with spectroscopy and a performance that is superior to that of Dixon methods at low FFs.

CONCLUSION

Fat results in negative or positive shifts in native tissue T₁ measured with MOLLI and SASHA over a narrow range of off-resonance frequencies; T₁ shifts from fat can be used to accurately quantify FF. Magn Reson Med 77:237-249, 2017. © 2016 Wiley Periodicals, Inc.

Selective spectroscopic imaging of hyperpolarized pyruvate and its metabolites using a single-echo variable phase advance method in balanced SSFP

PURPOSE

In balanced steady state free precession (bSSFP), the signal intensity has a well-known dependence on the off-resonance frequency, or, equivalently, the phase advance between successive radiofrequency (RF) pulses. The signal profile can be used to resolve the contributions from the spectrally separated metabolites. This work describes a method based on use of a variable RF phase advance to acquire spatial and spectral data in a time-efficient manner for hyperpolarized 13C MRI.

THEORY AND METHODS

The technique relies on the frequency response from a bSSFP acquisition to acquire relatively rapid, high-resolution images that may be reconstructed to separate contributions from different metabolites. The ability to produce images from spectrally separated metabolites was demonstrated in vitro, as well as in vivo following administration of hyperpolarized 1-13C pyruvate in mice with xenograft tumors.

RESULTS

In vivo images of pyruvate, alanine, pyruvate hydrate, and lactate were reconstructed from four images acquired in 2 s with an in-plane resolution of 1.25 × 1.25 mm(2) and 5 mm slice thickness.

CONCLUSION

The phase advance method allowed acquisition of spectroscopically selective images with high spatial and temporal resolution. This method provides an alternative approach to hyperpolarized 13C spectroscopic MRI that can be combined with other techniques such as multiecho or fluctuating equilibrium bSSFP. Magn Reson Med 76:1102-1115, 2016. © 2015 Wiley Periodicals, Inc.

Dynamically phase-cycled radial balanced SSFP imaging for efficient banding removal

PURPOSE

Balanced steady-state free precession (bSSFP) imaging suffers from banding artifacts due to its inherent sensitivity to inhomogeneities in the main magnetic field. These artifacts can be removed by the acquisition of multiple images at different frequency offsets. However, conventional phase-cycling is hindered by a long scan time. The purpose of this work is to present a novel approach for efficient banding removal in bSSFP imaging.

THEORY AND METHODS

To this end, the phase-cycle during a single-shot radial acquisition of an image was dynamically changed. Thus, each projection is acquired with a different frequency offset. Using conventional radial gridding, an artifact-free image can be reconstructed out of this dataset.

RESULTS

The approach is validated at clinical field strength [3.0 Tesla (T)] as well as at ultrahigh field (9.4T). Robust elimination of banding artifacts was obtained for different imaging regions, including brain imaging at ultrahigh field with an in-plane resolution of 0.25 × 0.25 mm(2). Besides banding artifact-free imaging, the applicability of the proposed technique for fat-water separation is demonstrated.

CONCLUSION

Dynamically phase-cycled radial bSSFP has the potential for banding-free bSSFP imaging in a short scan time, in the presence of severe field inhomogeneities and at high resolution.

Challenges of High-resolution Diffusion Imaging of the Human Medial Temporal Lobe in Alzheimer Disease

The human medial temporal lobe performs an essential role in memory formation and retrieval. Diseases involving the hippocampus such as Alzheimer disease present a unique opportunity for advanced imaging techniques to detect abnormalities at an early stage. In particular, it is possible that diffusion imaging will measure abnormal microarchitecture beyond the realm of macroscopic imaging. However, this task is formidable because of the detailed anatomy of the medial temporal lobe, the difficulties in obtaining high-quality diffusion images of adequate resolution, and the challenges in diffusion data processing. Moreover, it is unclear if any differences will be significant for an individual patient or simply groups of patients. Successful endeavors will need to address each of these challenges in an integrated fashion. The rewards of such analysis may be detection of microscopic disease in vivo, which could represent a landmark accomplishment for the field of neuroradiology.

Magnetic Resonance Characterization of Porous Media Using Diffusion through Internal Magnetic Fields

When a porous material is inserted into a uniform magnetic field, spatially varying fields typically arise inside the pore space due to susceptibility contrast between the solid matrix and the surrounding fluid. As a result, direct measurement of the field variation may provide a unique opportunity to characterize the pore geometry. The sensitivity of nuclear magnetic resonance (NMR) to inhomogeneous field variations through their dephasing effects on diffusing spins is unique and powerful. Recent theoretical and experimental research sheds new light on how to utilize susceptibility-induced internal field gradients to quantitatively probe the microstructure of porous materials. This article reviews ongoing developments based on the stimulated echo-pulse sequence to extend the characterization of porous media using both spatially resolved and unresolved susceptibility-induced internal gradients that operate on a diffusing-spin ensemble.

Improving contrast to noise ratio of resonance frequency contrast images (phase images) using balanced steady-state free precession

Recent MRI studies have exploited subtle magnetic susceptibility differences between brain tissues to improve anatomical contrast and resolution. These susceptibility differences lead to resonance frequency shifts which can be visualized by reconstructing the signal phase in conventional gradient echo (GRE) acquisition techniques. In this work, a method is proposed to improve the contrast to noise ratio per unit time (CNR efficiency) of anatomical MRI based on resonance frequency contrast. The method, based on the balanced steady-state free precession (bSSFP) MRI acquisition technique, was evaluated in its ability to generate contrast between gray and white matter in human brain at 3T and 7T. The results show substantially improved CNR efficiency of bSSFP phase images (2.85±0.21 times at 3 T and 1.71±0.11 times at 7 T) compared to the GRE data in a limited spatial area. This limited spatial coverage is attributed to the sensitivity of bSSFP to macroscopic B(0) inhomogeneities. With this CNR improvement, high resolution bSSFP phase images (resolution=0.3×0.3×2 mm(3), acquisition time=10min) acquired at 3T had sufficient CNR to allow the visualization of cortical laminar structures in invivo human primary visual cortex. Practical application of the proposed method may require improvement of B(0) homogeneity and stability by additional preparatory scans and/or compensation schemes such as respiration and drift compensation. Without these additions, the CNR benefits of the method may be limited to studies at low field or limited regions of interest.