Rapid quantitation of magnetization transfer using pulsed off-resonance irradiation and echo planar imaging

B1 inhomogeneity-corrected T1 mapping and quantitative magnetization transfer imaging via simultaneously estimating Bloch-Siegert shift and magnetization transfer effects

PURPOSE

To introduce a method of inducing Bloch-Siegert shift and magnetization Transfer Simultaneously (BTS) and demonstrate its utilization for measuring binary spin-bath model parameters free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular fraction (f $$ f $$ ), magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) and local transmit field (B 1 + $$ {B}_1^{+} $$ ).

THEORY AND METHODS

Bloch-Siegert shift and magnetization transfer is simultaneously induced through the application of off-resonance irradiation in between excitation and acquisition of an RF-spoiled gradient-echo scheme. Applying the binary spin-bath model, an analytical signal equation is derived and verified through Bloch simulations. Monte Carlo simulations were performed to analyze the method's performance. The estimation of the binary spin-bath parameters withB 1 + $$ {B}_1^{+} $$ compensation was further investigated through experiments, both ex vivo and in vivo.

RESULTS

Comparing BTS with existing methods, simulations showed that existing methods can significantly biasT 1 $$ {T}_1 $$ estimation when not accounting for transmitB 1 $$ {B}_1 $$ heterogeneity and MT effects that are present. Phantom experiments further showed that the degree of this bias increases with increasing macromolecular proton fraction. Multi-parameter fit results from an in vivo brain study generated values in agreement with previous literature. Based on these studies, we confirmed that BTS is a robust method for estimating the binary spin-bath parameters in macromolecule-rich environments, even in the presence ofB 1 + $$ {B}_1^{+} $$ inhomogeneity.

CONCLUSION

A method of estimating Bloch-Siegert shift and magnetization transfer effect has been developed and validated. Both simulations and experiments confirmed that BTS can estimate spin-bath parameters (T 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ) that are free fromB 1 + $$ {B}_1^{+} $$ bias.

High-resolution magnetization-transfer imaging of post-mortem marmoset brain: Comparisons with relaxometry and histology

Cell membranes and macromolecules or paramagnetic compounds interact with water proton spins, which modulates magnetic resonance imaging (MRI) contrast providing information on tissue composition. For a further investigation, quantitative magnetization transfer (qMT) parameters (at 3T), including the ratio of the macromolecular and water proton pools, F, and the exchange-rate constant as well as the (observed) longitudinal and the effective transverse relaxation rates (at 3T and 7T), R₁^obs and R₂^*, respectively, were measured at high spatial resolution (200 µm) in a slice of fixed marmoset brain and compared to histology results obtained with Gallyas' myelin stain and Perls' iron stain. R₁^obs and R₂^* were linearly correlated with the iron content for the entire slice, whereas distinct differences were obtained between gray and white matter for correlations of relaxometry and qMT parameters with myelin content. The combined results suggest that the macromolecular pool interacting with water consists of myelin and (less efficient) non-myelin contributions. Despite strong correlation of F and R₁^obs, none of these parameters was uniquely specific to myelination. Due to additional sensitivity to iron stores, R₁^obs and R₂^* were more sensitive for depicting microstructural differences between cortical layers than F.

Age-related differences in myeloarchitecture measured at 7 T

We have used the magnetisation transfer (MT) MRI measure as a primary measure of myelination in both the gray matter (GM) of the 78 cortical automated anatomical labeling (AAL) regions of the brain, and the underlying white matter in each region, in a cohort of healthy adults (aged 19-62 year old). The results revealed a significant quadratic trend in myelination with age, with average global myelination peaking at 42.9 year old in gray matter, and at 41.7 year old in white matter. We also explored the possibility of using the Nuclear Overhauser Enhancement (NOE) effect, which is acquired in a similar method to MT, as an additional measure of myelination. We found that the MT and NOE signals were strongly correlated in the brain and that the NOE effects displayed similar (albeit weaker) parabolic trends with age. We also investigated differences in cortical thickness with age, and confirmed a previous result of a linear decline of 4.5 ± 1.2 μm/y.

In vivo measurement of a new source of contrast, the dipolar relaxation time, T1D , using a modified inhomogeneous magnetization transfer (ihMT) sequence

PURPOSE

This paper describes a technique that can be used in vivo to measure the dipolar relaxation time, T_1D , of macromolecular protons contributing to magnetization transfer (MT) in tissues and to produce quantitative T_1D maps.

THEORY AND METHODS

The technique builds upon the inhomogeneous MT (ihMT) technique that is particularly sensitive to tissue components with long T_1D . A standard ihMT experiment was altered to introduce a variable time for switching between positive and negative offset frequencies for RF saturation. A model for the dependence of ihMT was developed and used to fit data acquired in vivo.

RESULTS

Application of the method to images from brains of healthy volunteers produced values of T_1D  = (5.9 ± 1.2) ms in gray matter and T_1D  = (6.2 ± 0.4) ms in white matter regions and provided maps of the T_1D parameter.

CONCLUSION

The model and experiments described provide access to a new relaxation characteristic of tissue with potentially unique diagnostic information. Magn Reson Med 78:1362-1372, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

High-resolution imaging of magnetisation transfer and nuclear Overhauser effect in the human visual cortex at 7 T

The aim of this study was to optimise a pulse sequence for high-resolution imaging sensitive to the effects of conventional macromolecular magnetisation transfer (MT(m)) and nuclear Overhauser enhancement (NOE), and to use it to investigate variations in these parameters across the cerebral cortex. A high-spatial-resolution magnetisation transfer-prepared turbo field echo (MT-TFE) sequence was designed to have high sensitivity to MT(m) and NOE effects, whilst being robust to B0 and B1 inhomogeneities, and producing a good point spread function across the cortex. This was achieved by optimising the saturation and imaging components of the sequence using simulations based on the Bloch equations, including exchange and an image simulator. This was used to study variations in these parameters across the cortex. Using the sequence designed to be sensitive to NOE and MT(m), a variation in signals corresponding to a variation in MT(m) and NOE across the cortex, consistent with a reduction in myelination from the white matter surface to the pial surface of the cortex, was observed. In regions in which the stria was visible on T2*-weighted images, it could also be detected in signals sensitive to MT(m) and NOE. There was greater variation in signals sensitive to NOE, suggesting that the NOE signal is more sensitive to myelination. A sequence has been designed to image variations in MT(m) and NOE at high spatial resolution and has been used to investigate variations in contrast in these parameters across the cortex.

In vivo imaging of glucose uptake and metabolism in tumors

Tumors have a greater reliance on anaerobic glycolysis for energy production than normal tissues. We developed a noninvasive method for imaging glucose uptake in vivo that is based on magnetic resonance imaging and allows the uptake of unlabeled glucose to be measured through the chemical exchange of protons between hydroxyl groups and water. This method differs from existing molecular imaging methods because it permits detection of the delivery and uptake of a metabolically active compound in physiological quantities. We show that our technique, named glucose chemical exchange saturation transfer (glucoCEST), is sensitive to tumor glucose accumulation in colorectal tumor models and can distinguish tumor types with differing metabolic characteristics and pathophysiologies. The results of this study suggest that glucoCEST has potential as a useful and cost-effective method for characterizing disease and assessing response to therapy in the clinic.

Matrix-algebra-based calculations of the time evolution of the binary spin-bath model for magnetization transfer

Quantification of magnetization-transfer (MT) experiments are typically based on the assumption of the binary spin-bath model. This model allows for the extraction of up to six parameters (relative pool sizes, relaxation times, and exchange rate constants) for the characterization of macromolecules, which are coupled via exchange processes to the water in tissues. Here, an approach is presented for estimating MT parameters acquired with arbitrary saturation schemes and imaging pulse sequences. It uses matrix algebra to solve the Bloch-McConnell equations without unwarranted simplifications, such as assuming steady-state conditions for pulsed saturation schemes or neglecting imaging pulses. The algorithm achieves sufficient efficiency for voxel-by-voxel MT parameter estimations by using a polynomial interpolation technique. Simulations, as well as experiments in agar gels with continuous-wave and pulsed MT preparation, were performed for validation and for assessing approximations in previous modeling approaches. In vivo experiments in the normal human brain yielded results that were consistent with published data.

Higher order reconstruction for MRI in the presence of spatiotemporal field perturbations

Despite continuous hardware advances, MRI is frequently subject to field perturbations that are of higher than first order in space and thus violate the traditional k-space picture of spatial encoding. Sources of higher order perturbations include eddy currents, concomitant fields, thermal drifts, and imperfections of higher order shim systems. In conventional MRI with Fourier reconstruction, they give rise to geometric distortions, blurring, artifacts, and error in quantitative data. This work describes an alternative approach in which the entire field evolution, including higher order effects, is accounted for by viewing image reconstruction as a generic inverse problem. The relevant field evolutions are measured with a third-order NMR field camera. Algebraic reconstruction is then formulated such as to jointly minimize artifacts and noise in the resulting image. It is solved by an iterative conjugate-gradient algorithm that uses explicit matrix-vector multiplication to accommodate arbitrary net encoding. The feasibility and benefits of this approach are demonstrated by examples of diffusion imaging. In a phantom study, it is shown that higher order reconstruction largely overcomes variable image distortions that diffusion gradients induce in EPI data. In vivo experiments then demonstrate that the resulting geometric consistency permits straightforward tensor analysis without coregistration.

Magnetization transfer phenomenon in the human brain at 7 T

Magnetization transfer is an important source of contrast in magnetic resonance imaging which is sensitive to the concentration of macromolecules and other solutes present in the tissue. Magnetization transfer effects can be visualized in magnetization transfer ratio images or quantified via the z-spectrum. This paper presents methods of measuring the z-spectrum and of producing high-resolution MTR images and maps of z-spectrum asymmetry in vivo at 7 T, within SAR limits. It also uses a 3-compartment model to measure chemical exchange and magnetization transfer parameters from the z-spectrum data. The peak in the z-spectrum associated with chemical exchange between amide and water protons (amide proton transfer, APT, effects) is much more apparent at 7 T than at 3 T. Furthermore at 7 T quantitative APT results varied between the corpus callosum and other white matter structures, suggesting that quantitative APT imaging could be used as a method of measuring myelination. The results also suggest that chemical exchange is not responsible for the phase shift observed in susceptibility weighted images between grey matter and white matter.

MT effects and T1 quantification in single-slice spoiled gradient echo imaging

We investigated magnetization transfer (MT) effects on the steady-state MR signal for a sample subjected to a series of identical on-resonance RF pulses, such as would be experienced while imaging a single slice using a spoiled gradient echo sequence. The MT coupling terms for a two-pool system were added to the Bloch equations and we derived the resulting steady-state signal equation and compared this result to the conventional signal equation without MT effects. The steady-state signal is increased by a few percent of the equilibrium magnetization because of MT. One consequence of this MT effect is inaccuracy in T(1) values determined via conventional steady-state gradient echo methods. (Theory predicts greater than 10% errors in T(1) for white matter when using short TR.) A second consequence is the ability to quantify the relaxation and MT parameters by fitting the gradient echo steady state signal to the signal equation appropriately modified to include MT effects. The theory was tested in samples of MnCl(2), cross-linked bovine serum albumin (BSA), and cross-linked BSA + MnCl(2).

Quantitative magnetization transfer imaging via selective inversion recovery with short repetition times

Quantitative magnetization transfer imaging (qMTI) methods are able to estimate fundamental sample parameters, such as the relative size of the solid-like macromolecular proton pool and the spin exchange rate between this pool and the directly measured free water protons. One such method is selective inversion recovery (SIR), in which the free water protons are selectively inverted and the signal is fit to a biexponential function of the inversion time (TI). SIR uses only low-power pulses and requires no separate RF (B1) or static field (B0) field maps, and the analysis is largely independent of the macromolecular pool lineshape. These are all advantages over steady-state off-resonance saturation qMTI methods. However, up to now, SIR has been implemented only with repetition times TR>>T1. This paper describes a modification of SIR with smaller TR values and a greater signal-to-noise ratio (SNR) efficiency.

Quantitative magnetization transfer by trains of radio frequency pulses in human brain: extension of a free evolution model to continuous-wave-like conditions

A theoretical model of free evolution between repeated magnetic transfer (MT) pulses was extended to continuous-wave (CW)-like conditions showing that only the repetitive "direct" saturation of bulk water changes the transient and stationary behavior. The influence of the pulse repetition period (PR) on progressive saturation was studied in cortical gray matter (GM) and central white matter (WM) under conditions of short periods of free evolution and strong macromolecular saturation. Interpulse delays of 3 ms were achieved in vivo on a 1.5-T MR system with bell-shaped MT pulses of 12-ms duration and nominal flip angles of up to 1440 degrees and single-shot readout by a stimulated echo acquisition mode localization sequence. The frequency offset was chosen between 1 and 3 kHz to avoid excessive direct saturation. The stationary MT ratio (MTR) followed an inverse linear PR dependence, showing a consistent partial saturation of about 90% at zero PR for both WM and GM. Comparison to a relaxation-matched liquid indicated the presence of MT, but not necessarily of direct saturation. The transient behavior indicated considerable direct saturation, but this could also be explained by MT. These inconsistencies showed that the intervals of time evolution in our experiments were too long to be modeled by CW-like conditions. Free evolution takes place during the whole PR rather than during the interpulse delay only. Quantification using the rates of free evolution theory yielded the saturations and rate constants necessary to explain the observed behavior. The theory of rapid CW-like pulsing provides an upper limit for the rate of progressive saturation. This limit is approached at PR below an estimated value of 5 ms. The phenomenological PR dependence of the steady-state MTR may indicate that MT exceeded the direct saturation. Unlike to an idealized CW experiment, the extrapolated value at zero PR is subject to direct effects and not a physically meaningful constant.