Quantitative imaging of magnetization transfer using multiple selective pulses

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

SIMPLEX: Multiple phase-cycled bSSFP quantitative magnetization transfer imaging with physic-guided simulation learning of neural network

Most quantitative magnetization transfer (qMT) imaging methods require acquiring additional quantitative maps (such as T1) for data fitting. A method based on multiple phase-cycled bSSFP was recently proposed to enable high-resolution 3D qMT imaging based on least square fitting without any extra acquisition, and thus has high potential for simplifying the qMT procedure. However, the quantification of qMT parameters with this method was suboptimal, limiting its potential for clinical application despite its simpler protocol and higher spatial resolution. To improve the fitting of qMT data obtained with multiple phase-cycled bSSFP, we propose SIMulation-based Physics-guided Learning of neural network for qMT parameters EXtraction, or SIMPLEX. In contrast to previous deep learning supervised approaches for quantitative MR that require the acquisition of input data and corresponding ground truth for training, we leveraged the MR signal model to generate training samples without expensive data curation. The network was trained exclusively with simulation data by predicting the simulation parameters. The same network was applied directly to in-vivo data without additional training. The approach was verified with both simulation and in-vivo data. SIMPLEX showed a decrease in fitting mean squared error for all simulation data compared to the existing least-square fitting method. The in-vivo experiment revealed that the network performed well with the real in vivo data unseen during training. For all experiments, we observed that SIMPLEX consistently improved the quantification quality of the qMT parameters whilst being more robust to noise compared to the prior technique. The proposed SIMPLEX will expedite the routine clinical application of qMT by providing qMT parameters (exchange rate, pool fraction) as well as T1, T2, and ΔB0 maps simultaneously with high spatial resolution, better reliability, and reduced processing time.

Macromolecular proton fraction mapping based on spin-lock magnetic resonance imaging

PURPOSE

In MRI, the macromolecular proton fraction (MPF) is a key parameter of magnetization transfer (MT). It represents the relative amount of immobile protons associated with semi-solid macromolecules involved in MT with free water protons. We aim to quantify MPF based on spin-lock MRI and explore its advantages over the existing MPF-mapping methods.

METHODS

In the proposed method, termed MPF quantification based on spin-lock (MPF-SL), off-resonance spin-lock is used to sensitively measure the MT effect. MPF-SL is designed to measure a relaxation rate (R_mpfsl ) that is specific to the MT effect by removing the R_1ρ relaxation due to the mobile water and chemical exchange pools. A theory is derived to quantify MPF from the measured R_mpfsl . No prior knowledge of tissue relaxation parameters, including T₁ or T₂ , is needed to quantify MPF using MPF-SL. The proposed approach is validated with Bloch-McConnell simulations, phantom, and in vivo liver studies at 3.0T.

RESULTS

Both Bloch-McConnell simulations and phantom experiments show that MPF-SL is insensitive to variations of the mobile water pool and the chemical exchange pool. MPF-SL is specific to the MT effect and can measure MPF reliably. In vivo liver studies show that MPF-SL can be used to detect collagen deposition in patients with liver fibrosis.

CONCLUSION

A novel MPF imaging method based on spin-lock MRI is proposed. The confounding factors are removed, and the measurement is specific to the MT effect. It holds promise for MPF-sensitive diagnostic imaging in clinical settings.

Towards an analytic solution for pulsed CEST

Chemical exchange saturation transfer (CEST) is an imaging method based on magnetization exchange between solutes and water. This exchange generates changes in the measured signal after off-resonance radiofrequency irradiation. Although the analytic solution for CEST with continuous wave (CW) irradiation has been determined, most studies are performed using pulsed irradiation. In this work, we derive an analytic solution for the CEST signal after pulsed irradiation that includes both short-time rotation effects and long-time saturation effects in a two-pool system corresponding to water and a low-concentration exchanging solute pool. Several approximations are made to balance the accuracy and simplicity of the resulting analytic form, which is tested against numerical solutions of the coupled Bloch equations and is found to be largely accurate for amides at high fields, but less accurate at the higher exchange rates, lower offsets and typically higher irradiation powers of amines.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Histological validation of fast macromolecular proton fraction mapping as a quantitative myelin imaging method in the cuprizone demyelination model

Cuprizone-induced demyelination in mice is a frequently used model in preclinical multiple sclerosis research. A recent quantitative clinically-targeted MRI method, fast macromolecular proton fraction (MPF) mapping demonstrated a promise as a myelin biomarker in human and animal studies with a particular advantage of sensitivity to both white matter (WM) and gray matter (GM) demyelination. This study aimed to histologically validate the capability of MPF mapping to quantify myelin loss in brain tissues using the cuprizone demyelination model. Whole-brain MPF maps were obtained in vivo on an 11.7T animal MRI scanner from 7 cuprizone-treated and 7 control С57BL/6 mice using the fast single-point synthetic-reference method. Brain sections were histologically stained with Luxol Fast Blue (LFB) for myelin quantification. Significant (p < 0.05) demyelination in cuprizone-treated animals was found according to both LFB staining and MPF in all anatomical structures (corpus callosum, anterior commissure, internal capsule, thalamus, caudoputamen, and cortex). MPF strongly correlated with quantitative histology in all animals (r = 0.95, p < 0.001) as well as in treatment and control groups taken separately (r = 0.96, p = 0.002 and r = 0.93, p = 0.007, respectively). Close agreement between histological myelin staining and MPF suggests that fast MPF mapping enables robust and accurate quantitative assessment of demyelination in both WM and GM.

Incorporating dixon multi-echo fat water separation for novel quantitative magnetization transfer of the human optic nerve in vivo

PURPOSE

The optic nerve (ON) represents the sole pathway between the eyes and brain; consequently, diseases of the ON can have dramatic effects on vision. However, quantitative magnetization transfer (qMT) applications in the ON have been limited to ex vivo studies, in part because of the fatty connective tissue that surrounds the ON, confounding the magnetization transfer (MT) experiment. Therefore, the aim of this study was to implement a multi-echo Dixon fat-water separation approach to remove the fat component from MT images.

METHODS

MT measurements were taken in a single slice of the ON and frontal lobe using a three-echo Dixon readout, and the water and out-of-phase images were applied to a two-pool model in ON tissue and brain white matter to evaluate the effectiveness of using Dixon fat-water separation to remove fatty tissue from MT images.

RESULTS

White matter data showed no significant differences between image types; however, there was a significant increase (p < 0.05) in variation in the out-of-phase images in the ON relative to the water images.

CONCLUSIONS

The results of this study demonstrate that Dixon fat-water separation can be robustly used for accurate MT quantification of anatomies susceptible to partial volume effects resulting from fat. Magn Reson Med 77:707-716, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Effects of magnetization transfer on T1 contrast in human brain white matter

MRI based on T1 relaxation contrast is increasingly being used to study brain morphology and myelination. Although it provides for excellent distinction between the major tissue types of gray matter, white matter, and CSF, reproducible quantification of T1 relaxation rates is difficult due to the complexity of the contrast mechanism and dependence on experimental details. In this work, we perform simulations and inversion-recovery MRI measurements at 3T and 7T to show that substantial measurement variability results from unintended and uncontrolled perturbation of the magnetization of MRI-invisible (1)H protons of lipids and macromolecules. This results in bi-exponential relaxation, with a fast component whose relative contribution under practical conditions can reach 20%. This phenomenon can strongly affect apparent relaxation rates, affect contrast between tissue types, and result in contrast variations over the brain. Based on this novel understanding, ways are proposed to minimize this experimental variability and its effect on T1 contrast, quantification accuracy and reproducibility.

A review of optimization and quantification techniques for chemical exchange saturation transfer MRI toward sensitive in vivo imaging

Chemical exchange saturation transfer (CEST) MRI is a versatile imaging method that probes the chemical exchange between bulk water and exchangeable protons. CEST imaging indirectly detects dilute labile protons via bulk water signal changes following selective saturation of exchangeable protons, which offers substantial sensitivity enhancement and has sparked numerous biomedical applications. Over the past decade, CEST imaging techniques have rapidly evolved owing to contributions from multiple domains, including the development of CEST mathematical models, innovative contrast agent designs, sensitive data acquisition schemes, efficient field inhomogeneity correction algorithms, and quantitative CEST (qCEST) analysis. The CEST system that underlies the apparent CEST-weighted effect, however, is complex. The experimentally measurable CEST effect depends not only on parameters such as CEST agent concentration, pH and temperature, but also on relaxation rate, magnetic field strength and more importantly, experimental parameters including repetition time, RF irradiation amplitude and scheme, and image readout. Thorough understanding of the underlying CEST system using qCEST analysis may augment the diagnostic capability of conventional imaging. In this review, we provide a concise explanation of CEST acquisition methods and processing algorithms, including their advantages and limitations, for optimization and quantification of CEST MRI experiments.

Quantitative MRI and ultrastructural examination of the cuprizone mouse model of demyelination

The cuprizone mouse model of demyelination was used to investigate the influence that white matter changes have on different magnetic resonance imaging results. In vivo T2 -weighted and magnetization transfer images (MTIs) were acquired weekly in control (n = 5) and cuprizone-fed (n = 5) mice, with significant increases in signal intensity in T2 -weighted images (p < 0.001) and lower magnetization transfer ratio (p < 0.001) in the corpus callosum of the cuprizone-fed mice starting at 3 weeks and peaking at 4 and 5 weeks, respectively. Diffusion tensor imaging (DTI), quantitative MTI (qMTI), and T1/T2 measurements were used to analyze freshly excised tissue after 6 weeks of cuprizone administration. In multicomponent T2 analysis with 10 ms echo spacing, there was no visible myelin water component associated with the short T2 value. Quantitative MTI metrics showed significant differences in the corpus callosum and external capsule of the cuprizone-fed mice, similar to previous studies of multiple sclerosis in humans and animal models of demyelination. Fractional anisotropy was significantly lower and mean, axial, and radial diffusivity were significantly higher in the cuprizone-fed mice. Cellular distributions measured in electron micrographs of the corpus callosum correlated strongly to several different quantitative MRI metrics. The largest Spearman correlation coefficient varied depending on cellular type: T1 versus the myelinated axon fraction (ρ = -0.90), the bound pool fraction (ƒ) versus the myelin sheath fraction (ρ = 0.93), and axial diffusivity versus the non-myelinated cell fraction (ρ = 0.92). Using Pearson's correlation coefficient, ƒ was strongly correlated to the myelin sheath fraction (r = 0.98) with a linear equation predicting myelin content (5.37ƒ - 0.25). Of the calculated MRI metrics, ƒ was the strongest indicator of myelin content, while longitudinal relaxation rates and diffusivity measurements were the strongest indicators of changes in tissue structure.

Quantitative MR imaging of two-pool magnetization transfer model parameters in myelin mutant shaking pup

Magnetization transfer (MT) imaging quantitatively assesses cerebral white matter disease through its sensitivity to macromolecule-bound protons including those associated with myelin proteins and lipid bilayers. However, traditional MT contrast measured by the MT ratio (MTR) lacks pathologic specificity as demyelination, axon loss, inflammation and edema all impact MTR, directly and/or indirectly through multiple covariances among imaging parameters (particularly MTR with T(1)) and tissue features (e.g. axon loss with demyelination). In this study, more complex modeling of MT phenomena ("quantitative" MT or qMT) was applied to a less complex disease model (the myelin mutant shaking [sh] pup, featuring hypomyelination but neither inflammation nor axon loss) in order to eliminate the covariances on both sides of the MR-pathology "equation" and characterize these important relationships free from the usual confounds. qMT measurements were acquired longitudinally in 6 sh pups and 4 age-matched controls ranging from 3 to 21 months of age and compared with histology. The qMT parameter, bound pool fraction (f), was the most distinctive between diseased and control animals; both f and longitudinal relaxation rate R(1) tracked myelination with normal aging, whereas MTR did not--presumably owing to counterbalancing MT and R(1) effects. qMT imaging provides a more accurate and potentially more specific non-invasive tissue characterization.

Characterization of cerebral white matter properties using quantitative magnetic resonance imaging stains

The image contrast in magnetic resonance imaging (MRI) is highly sensitive to several mechanisms that are modulated by the properties of the tissue environment. The degree and type of contrast weighting may be viewed as image filters that accentuate specific tissue properties. Maps of quantitative measures of these mechanisms, akin to microstructural/environmental-specific tissue stains, may be generated to characterize the MRI and physiological properties of biological tissues. In this article, three quantitative MRI (qMRI) methods for characterizing white matter (WM) microstructural properties are reviewed. All of these measures measure complementary aspects of how water interacts with the tissue environment. Diffusion MRI, including diffusion tensor imaging, characterizes the diffusion of water in the tissues and is sensitive to the microstructural density, spacing, and orientational organization of tissue membranes, including myelin. Magnetization transfer imaging characterizes the amount and degree of magnetization exchange between free water and macromolecules like proteins found in the myelin bilayers. Relaxometry measures the MRI relaxation constants T1 and T2, which in WM have a component associated with the water trapped in the myelin bilayers. The conduction of signals between distant brain regions occurs primarily through myelinated WM tracts; thus, these methods are potential indicators of pathology and structural connectivity in the brain. This article provides an overview of the qMRI stain mechanisms, acquisition and analysis strategies, and applications for these qMRI stains.

Fast macromolecular proton fraction mapping from a single off-resonance magnetization transfer measurement

A new method was developed for fast quantitative mapping of the macromolecular proton fraction defined within the two-pool model of magnetization transfer. The method utilizes a single image with off-resonance saturation, a reference image for data normalization, and T(1), B(0), and B(1) maps with the total acquisition time ~10 min for whole-brain imaging. Macromolecular proton fraction maps are reconstructed by iterative solution of the matrix pulsed magnetization transfer equation with constrained values of other model parameters. Theoretical error model describing the variance due to noise and the bias due to deviations of constrained parameters from their actual values was formulated based on error propagation rules. The method was validated by comparison with the conventional multiparameter multipoint fit of the pulsed magnetization transfer model based on data from two healthy subjects and two multiple sclerosis patients. It was demonstrated theoretically and experimentally that accuracy of the method depends on the offset frequency and flip angle of the saturation pulse, and optimal ranges of these parameters are 4-7 kHz and 600°-900°, respectively. At optimal sampling conditions, the single-point method enables <10% relative macromolecular proton fraction errors. Comparison with the multiparameter fitting method revealed very good agreement with no significant bias and limits of agreement around ± 0.7%.

Simulation and optimization of pulsed radio frequency irradiation scheme for chemical exchange saturation transfer (CEST) MRI-demonstration of pH-weighted pulsed-amide proton CEST MRI in an animal model of acute cerebral ischemia

Chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) is capable of measuring dilute labile protons and microenvironmental properties. However, the CEST contrast is dependent upon experimental conditions-particularly, the radiofrequency (RF) irradiation scheme. Although continuous-wave RF irradiation has been used conventionally, the limited RF pulse duration or duty cycle of most clinical systems requires the use of pulsed RF irradiation. Here, the conventional numerical simulation is extended to describe pulsed-CEST MRI contrast as a function of RF pulse parameters (i.e., RF pulse duration and flip angle) and labile proton properties (i.e., exchange rate and chemical shift). For diamagnetic CEST agents undergoing slow or intermediate chemical exchange, simulation shows a linear regression relationship between the optimal mean RF power of pulsed-CEST MRI and continuous-wave-CEST MRI. The optimized pulsed-CEST contrast is approximately equal to that of continuous-wave-CEST MRI for exchange rates less than 50 s(-1), as confirmed experimentally using a multicompartment pH phantom. In the acute stroke animals, we showed that pulsed- and continuous-wave-amide proton CEST MRI demonstrated similar contrast. In summary, our study elucidated the RF irradiation dependence of pulsed-CEST MRI contrast, providing useful insights to guide its experimental optimization and quantification.

The MT pool size ratio and the DTI radial diffusivity may reflect the myelination in shiverer and control mice

A quantitative magnetization transfer (qMT) technique was employed to quantify the ratio of the sizes of the bound and free water proton pools in ex vivo mouse brains. The goal was to determine the pool size ratio sensitivity to myelin. Fixed brains from both shiverer mice and control littermates were imaged. The pool size ratio in the corpus callosum of shiverer mice was substantially lower than that in the control mice, while there was no distinguishable difference in the pool size ratio in the gray matter. These results correlate with diffusion tensor imaging (DTI) derived radial diffusivity which previously was shown to reflect myelin integrity in this animal model. Histological study reveals the presence of myelin in control mice white matter and the absence of myelin in shiverer mice white matter, supporting the qMT and DTI results. Our findings support the view that qMT may be used for estimating myelin integrity.

Quantitative magnetization transfer measured pool-size ratio reflects optic nerve myelin content in ex vivo mice

Optic nerves from mice that have undergone retinal ischemia were examined using a newly implemented quantitative magnetization transfer (qMT) technique. Previously published results indicate that the optic nerve from retinal ischemia mice suffered significant axon degeneration without detectable myelin injury at 3 days after reperfusion. At this time point, we acquired ex vivo qMT parameters from both shiverer mice (which have nearly no myelin) and control mice that have undergone retinal ischemia, and these qMT measures were compared with diffusion tensor imaging (DTI) results. Our findings suggests that the qMT estimated ratio of the pool sizes of the macromolecular and free water protons reflected the different myelin contents in the optic nerves between the shiverer and control mice. This pool size ratio was specific to myelin content only and was not significantly affected by the presence of axon injury in mouse optic nerve 3 days after retinal ischemia.

Magnetization transfer MR imaging in multiple sclerosis

We introduce the fundamental aspects of MT, of MT MR imaging, and the respective analysis techniques. We then review the applications of MT MR imaging to multiple sclerosis. Finally we review the technique's contribution to our understanding of this disease.

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.

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

A technique for producing a quantitative measure of magnetization transfer parameters in a clinically feasible time scale is proposed. The combination of pulsed off-resonance irradiation and echo planar imaging has produced an imaging sequence that negates the need for continuous wave irradiation and allows the approach to steady-state conditions to be studied. Data analysis involves the step-by-step numerical solution of the modified Bloch equations to generate a quantitative model of the measured signal intensity based on the relative size of the bound proton pool and the bound proton pool transverse relaxation time. The sequence and model are applied to the study of a series of agar gels of varying concentrations and the results are compared to those from the literature.

Simultaneous measurement of saturation and relaxation in human brain by repetitive magnetization transfer pulses

Magnetization transfer (MT) by equidistant pulse trains can be described as being analogous to progressive partial saturation, where 'direct' saturation of water is amplified by MT contributions that are dependent on macromolecular content and differential saturation. This concept was applied to study the transition to steady state in the human brain using similar MT-pulses as in imaging. Up to 41 bell-shaped MT-pulses of 12 ms duration were applied at frequency offsets between 0.5 and 15 kHz with flip angles between 1080 and 1440 degrees . Central white and parietal gray matter was studied in human subjects using STEAM for localized read-out (TE = 30 ms, TM = 13.7 ms). The apparent degree of saturation, delta(app), and the longitudinal relaxation of the water pool during the pulse repetition period (PR) were fitted to the transient behavior after signal correction for cerebro-spinal fluid. PR was varied between 15 and 100 ms to assess the PR-dependence of the fitted parameters. The MT-term in delta(app) exceeded the direct saturation and attained its maximum at PR > or = 100 ms. The macromolecular pool was only partially saturated by a single MT-pulse. The offset may be increased to 2.5 kHz to reduce direct saturation without sacrificing MT in white matter. The estimated relaxation rates (1.04 +/- 0.11 s(-1) in WM; 0.76 +/- 0.13 s(-1) in GM) were faster than are commonly observed at 1.5 T. The apparent saturation is a measure for MT that is not confounded by relaxation. To maximize MT in brain tissue, MT-pulses should be applied at PR = 100 ms or longer. At shorter PR, a larger steady state saturation is obtained at the cost of increased contributions from direct saturation. Since this accelerates the convergence, PR should be decreased to reach the steady state within a specified time. A faster transition can always be achieved at a reduced frequency offset via increased direct saturation.

Cross-relaxation imaging reveals detailed anatomy of white matter fiber tracts in the human brain

Cross-relaxation imaging is a new quantitative MRI modality, which allows mapping of fundamental parameters determining the magnetization transfer (MT) effect in tissues, cross-relaxation rate constant (k) and bound pool fraction (f). This study introduces a new time-efficient technique for cross-relaxation imaging, which obtains three-dimensional (3D) whole-brain k and f maps with scan time of <30 min and isotropic spatial resolution of 1.4 mm. The technical principle of the method is based on four-point fit of a matrix model of pulsed MT to imaging data obtained with variable offset frequency saturation while using a complimentary R1 (=1 / T1) map. Anatomical correlations of in vivo cross-relaxation parametric maps were evaluated in three healthy subjects. The f maps revealed correspondence of areas with highly elevated f = 12-15% to major fiber tracts such as corpus callosum, anterior commissure, optic radiations, and major brain fasciculi. The rest of white matter (WM) demonstrated lower f = 9-11%, resulting in clear visual contrast of fiber tracts. Even lower f = 6.5-8.5% were found in gray matter (GM) with the highest f = 8.5% in the anterior thalamus. Distribution of k was relatively uniform in WM and produced sharp contrast between GM and WM (k = 1.6 and 3.3 s(-1), respectively). The most marked feature of k maps was their ability to visualize the corticospinal tract, which had elevated k = 3.4-3.8 s(-1) but appeared invisible on f maps. The observed patterns on f maps can be explained by variations in the density of myelinated fibers, while the trends of k may reflect regional differences in axonal organization. Cross-relaxation imaging can be used in various clinical studies focused on brain development and white matter diseases.

Assessment and correction ofB1-induced errors in magnetization transfer ratio measurements

The magnetization transfer ratio (MTR) is strongly related to the field strength (B(1)) of the saturation pulse. B(1) variations therefore can result in significant MTR variations and can affect histogram analysis, particularly if data from a large volume of interest are included. A multicenter study was performed to determine the typical range of B(1) errors and the corresponding MTR variations in brain tissue of healthy volunteers. Seven subjects were included at each center resulting in a total cohort of 28 subjects. Additionally, numerical simulations were done to study this relationship more generally for pulsed saturation. It could be demonstrated, both theoretically and empirically, that for typical B(1) errors there is a linear relationship between B(1) error and the corresponding MTR change. In addition, for proton density-weighted sequences, this relationship seems to be largely independent of the underlying relaxation properties. Mean B(1) errors in the entire brain were typically in the range between -3% and -7%. Due to different coil characteristics, significant MTR differences between different scanners and sites were observed. Using a simple correction scheme that is based on a linear regression analysis between MTR and B(1) data it was possible to reduce the intersubject variation by approximately 50%. Furthermore, interscanner variation could be reduced such that no significant differences between scanners could be detected. The correction scheme may be useful when investigating MTR as an outcome measure in single or multicenter studies.

Method for quantitative imaging of the macromolecular 1H fraction in tissues

A new method was developed for mapping the relative density of the macromolecular protons involved in magnetization transfer (MT). This method employs a stimulated echo preparation scheme in order to modulate the phase distribution within a spin ensemble. This labeled spin ensemble is then used as an intrinsic indicator, which is diluted due to magnetization exchange with macromolecular protons. A pulse sequence is presented which compensates for longitudinal relaxation, allows observation of the dilution effect only, and provides for calculation of parameter maps using indicator dilution theory. Compared to other quantitative MT techniques, neither additional relaxation time measurements nor knowledge regarding the lineshape of the macromolecular proton pool are required. Moreover, the inherent low specific absorption rate and the low sensitivity for B(1) errors make this method favorable in a clinical setting. This sequence was used to measure the macromolecular proton density in cross-linked bovine serum albumin. Using a navigated echo planar readout, the sequence was also employed to visualize the macromolecular content of human brain in vivo.

Quantitative imaging of magnetization transfer using an inversion recovery sequence

A new imaging method has been developed for quantitatively measuring magnetization transfer (MT). It uses a simple inversion recovery sequence, although one with very short (milliseconds) inversion times, and thus can be implemented on clinical imaging systems with little modification to existing pulse sequences. The sequence requires an inversion pulse with a length much longer than T(2m) (typically 10 micros) and much shorter than T(2f) (typically tens of ms) and 1/k(mf) (typically tens of ms), where T(2m) and T(2f) are the transverse relaxation times of the immobile macromolecular and free water protons, respectively, and k(mf) is the rate of MT between these populations. The resultant NMR signal is sensitive to MT when this inversion pulse affects the mobile and immobile proton pools to different degrees and by appropriate analysis of the signals obtained for different inversion times, quantitative information can be derived on the macromolecular content and exchange rates within the sample. The method has been used in conjunction with echo planar imaging to produce maps of the spatial distribution of the macromolecular content and MT rate in cross-linked bovine serum albumin. Comparisons between this method and other quantitative MT techniques are discussed.

Precise estimate of fundamental in-vivo MT parameters in human brain in clinically feasible times

A methodology is presented for extracting precise quantitative MT parameters using a magnetisation-prepared spoiled gradient echo sequence. This method, based on a new mathematical model, provides relaxation parameters for human brain in-vitro and in-vivo. The in-vivo parameters have been obtained from three different regions of normal white matter: occipital white matter, frontal white matter and centrum semiovale; two regions of normal grey matter: cerebral cortex and cerebellum, and from five regions with MS lesions. All this has been achieved using MT images collected within a timeframe that is clinically feasible. We hope that this new technique will shed light on the properties and dynamics of water compartments within the brain.

Pulsed Z-spectroscopic imaging of cross-relaxation parameters in tissues for human MRI: theory and clinical applications

A new method of pulsed Z-spectroscopic imaging is proposed for in vivo visualization and quantification of the parameters describing cross-relaxation between protons with liquid-like and solid-like relaxation properties in tissues. The method is based on analysis of the magnetization transfer (MT) effect as a function of the offset frequency and amplitude of a pulsed off- resonance saturation incorporated in a spoiled gradient-echo MRI pulse sequence. The theoretical concept of the method relies on an approximated analytical model of pulsed MT that provides a simple three-parameter equation for a pulsed steady-state Z-spectrum taken far from resonance. Using this model, the parametric images of cross-relaxation rate constant, content, and T(2) of the semisolid proton fraction can be reconstructed from a series of MT-weighted images and a coregistered T(1) map. The method was implemented on a 0.5 T clinical MRI scanner, and it provided high-quality 3D parametric maps within an acceptable scanning time. The estimates of cross-relaxation parameters in brain tissues were shown to be quantitatively consistent with the literature data. Clinical examples of the parametric images of human brain pathologies (multiple sclerosis and glioma) demonstrated high tissue contrast and clear visualization of the lesions.

Quantitative imaging of magnetization transfer exchange and relaxation properties in vivo using MRI

We describe a novel imaging technique that yields all of the observable properties of the binary spin-bath model for magnetization transfer (MT) and demonstrate this method for in vivo studies of the human head. Based on a new model of the steady-state behavior of the magnetization during a pulsed MT-weighted imaging sequence, this approach yields parametric images of the fractional size of the restricted pool, the magnetization exchange rate, the T(2) of the restricted pool, as well as the relaxation times in the free pool. Validated experimentally on agar gels and samples of uncooked beef, we demonstrate the method's application on two normal subjects and a patient with multiple sclerosis.

Dose resolution optimization of polymer gel dosimeters using different monomers

Polymer gel dosimeters of different formulations were manufactured from different monomers of acrylamide, acrylic acid, methacrylic acid, 1-vinyl-2-pyrrolidinone, 2-hydroxyethyl methacrylate and 2-hydroxyethyl acrylate. Gelatin and agarose were used as the gelling agents and N,N'-methylene-bis-acrylamide was used as a co-monomer in each polymer gel dosimeter. The T2 dependence of each dosimeter was analysed using a model of fast exchange of magnetization. The influence of the half-dose and the apparent T2 of the polymer-proton pool on the dose resolution (Dpdelta) were examined. Comparisons are made with the commonly employed R2-dose sensitivity. Differences exist suggesting that experiments reported in the literature using what were thought to be more optimal dosimeters may not actually be so. Based on Dpdelta of each formulation, conclusions are drawn on the optimal formulation required for a specific range of absorbed doses. In addition, information about the extent of polymerization of the monomers used along with some characteristics of the polymer network formed are reported. The influence of the concentration of monomers and gelling agent was subsequently evaluated using a model of fast exchange of magnetization. Based on these calculations, further improvement in Dpdelta can be expected.

The relationship of problems in biomedical MRI to the study of porous media

The NMR methods that are used to characterize inanimate porous media measure relaxation times and related phenomena and material transport, fluid displacement and flow. Biological tissues are comprised of multiple small, fluid-filled compartments, such as cells, that restrict the movement of the bulk solvent water and whose constituents influence water proton relaxation times via numerous interactions with macromolecular surfaces. Several of the methods and concepts that have been developed in one field of application are also of great value in the other, and it may be expected that technical developments that have been spurred by biomedical applications of MR imaging will be used in the continuing study of porous media. Some recent specific studies from our laboratory include the development of multiple quantum coherence methods for studies of ordered water in anisotropic macromolecular assemblies, studies of the degree of restriction of water diffusion in cellular systems, multiple selective inversion imaging to depict the ratios of proton pool sizes and rates of magnetization transfer between proton populations, and diffusion tensor imaging to depict tissue anisotropies. These illustrate how approaches to obtain structural information from biological media are also relevant to porous media. For example, the recent development of oscillating gradient spin echo techniques (OGSE), an approach that extends our ability to resolve apparent diffusion changes over different time scales in tissues, has also been used to compute surface to volume measurements in assemblies of pores. Each of the new methods can be adapted to provide spatially resolved quantitative measurements of properties of interest, and these can be efficiently acquired with good accuracy using fast imaging methods such as echo planar imaging. The community of NMR scientists focused on applications to porous media should remain in close communication with those who use MRI to study problems in biomedicine, to their mutual benefits.

Studies of magnetization transfer and relaxation in irradiated polymer gels--interpretation of MRI-based dosimetry

Magnetization transfer and NMR relaxation rates were measured for water protons in two types of polymer gels developed for radiation dosimetry with MRI in order to quantify the contributions of different relaxation processes to the radiation response in such gels. Measurements included the rate of magnetization transfer between proton pools and the ratio of the sizes of exchanging pools, R1 and R2. A model of relaxation in irradiated gels is presented to explain their properties. The model incorporates three proton pools: free water, macromolecular and interfacial. Two pools are insufficient to model the data. In these systems, radiation-induced polymerization appears to increase the size of a solid-like macromolecular proton pool but does not affect the rate constant of magnetization transfer per proton from macromolecular protons to the free water protons. The relation between R1 and the pool size ratio is consistent with free water exchanging with a macromolecular pool with an R1 of approximately 8 Hz. In addition, the rate of magnetization transfer is not limited by the rate of chemical exchange between the free water and the interfacial protons, and magnetization transfer most probably occurs via labile proton exchange rather than via bound water molecules.

Estimation of magnetization transfer rates from PACE experiments with pulsed RF saturation

A new imaging method has been developed for estimating the magnetization transfer rate (MTR) in a biologic two-pool system such as the brain tissue. The transfer rate is calculated from the ratio of the MTR to T(1sat), where T(1sat) is the apparent longitudinal relaxation time under complete saturation of the macromolecular pool. MTR and T(1sat) maps were obtained with a phase acquisition of composite echo (PACE) technique combined with pulsed radiofrequency (RF) saturation. The influences of RF saturation power and frequency offset on quantitative results were investigated with phantom and in vivo measurements. In white matter of seven healthy volunteers we found a mean transfer rate of 1.5 sec(-1), where the highest transfer rate was found in the genu of the corpus callosum (k(f) = 1. 9 sec(-1)). It could be shown that conditions near to complete saturation can also be reached under common restrictions by the specific absorption rate.

Quantitative interpretation of magnetization transfer in spoiled gradient echo MRI sequences

A method for analyzing general pulsed magnetization transfer (MT) experiments in which off-resonance saturation pulses are interleaved with on-resonance excitation pulses is presented. We apply this method to develop a steady-state signal equation for MT-weighted spoiled gradient echo sequences and consider approximations that facilitate its rapid computation. Using this equation, we assess various experimental designs for quantitatively imaging the fractional size of the restricted pool, cross-relaxation rate, and T(1) and T(2) relaxation times of the two pools in a binary spin bath system. From experiments on agar gel, this method is shown to reliably and accurately estimate the exchange and relaxation properties of a material in an imaging context, suggesting the feasibility of using this technique in vivo.