Low power method for estimating the magnetization transfer bound-pool macromolecular fraction

Collagen degradation assessment with an in vitro rotator cuff tendinopathy model using multiparametric ultrashort-TE magnetization transfer (UTE-MT) imaging

PURPOSE

This study aims to assess ultrashort-TE magnetization transfer (UTE-MT) imaging of collagen degradation using an in vitro model of rotator cuff tendinopathy.

METHODS

Thirty-six supraspinatus tendon specimens were divided into three groups and treated with 600 U collagenase (Group 1), 150 U collagenase (Group 2), and phosphate buffer saline (Group 3). UTE-MT imaging was performed to assess changes in macromolecular fraction (MMF), macromolecule transverse relaxation time (T2m), water longitudinal relaxation rate constant (R1m), the magnetization exchange rate from the macromolecular to water pool (Rm0 w) and from water to the macromolecular pool (Rm0 m), and magnetization transfer ratio (MTR) at baseline and following digestion and their differences between groups. Biochemical and histological studies were conducted to determine the extent of collagen degradation. Correlation analyses were performed with MMF, T2m, R1m, Rm0 w, Rm0 m, and MTR, respectively. Univariate and multivariate linear regression analyses were performed to evaluate combinations of UTE-MT parameters to predict collagen degradation.

RESULTS

MMF, T2m, R1m, Rm0 m, and MTR decreased after digestion. MMF (r = -0.842, p < 0.001), MTR (r = -0.78, p < 0.001), and Rm0 m (r = -0.662, p < 0.001) were strongly negatively correlated with collagen degradation. The linear regression model of differences in MMF and Rm0 m before and after digestion explained 68.9% of collagen degradation variation in the tendon. The model of postdigestion in MMF and T2m and the model of MTR explained 54.2% and 52.3% of collagen degradation variation, respectively.

CONCLUSION

This study highlighted the potential of UTE-MT parameters for evaluation of supraspinatus tendinopathy.

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.

Quantitative magnetization transfer MRI of desmoplasia in pancreatic ductal adenocarcinoma xenografts

Quantitative assessment of desmoplasia in pancreatic ductal adenocarcinoma (PDAC) may be critical for staging or prediction of response to therapy. We performed quantitative magnetization transfer (qMT) MRI measurements in 18 mouse xenograft tumors generated from three PDAC cell lines. The qMT parameter bound proton fraction (BPF) was found to be significantly higher in tumors grown using the BxPC-3 cell line (5.31 ± 0.87, mean ± standard deviation) compared with the BPF measured for tumors grown from Panc-1 (3.65 ± 0.60) and Capan-1 (1.50 ± 0.58) cell lines (P < 0.05 for each comparison). Histologic measurements demonstrated a similar trend; BxPC-3 tumors had significantly higher fibrosis levels (percentage of fibrotic tissue area, 6.21 ± 2.10) compared with Panc-1 (2.88 ± 1.13) and Capan-1 (1.69 ± 1.01) tumors. BPF was well correlated with quantitative fibrosis levels (r = 0.77, P < 0.01). Our results indicate that qMT measurements offer the potential to noninvasively quantify fibrosis levels in PDAC mouse xenograft models and thus serve as a valuable in vivo biomarker of desmoplasia in PDAC.

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%.

Quantitative MR imaging: physical principles and sequence design in abdominal imaging

Quantitative magnetic resonance (MR) imaging seeks to quantify fundamental biologic and MR-inducible tissue properties, in contrast to the routine application of MR imaging in the clinic, in which differences in MR parameters are used to generate contrast for subsequent subjective image analysis. Fundamental parameters that are commonly quantified by using MR imaging include proton density, diffusion, T1 relaxation, T2 and T2* relaxation, and magnetization transfer. Applications of these MR imaging-quantifiable parameters to abdominal imaging include oncologic imaging, evaluation of diffuse liver disease, and assessment of splenic, renal, and pancreatic disease. An understanding of the inherent physical principles underlying the basic quantitative parameters as well as the commonly used pulse sequences requisite to their derivation is critical, as this field is rapidly growing and its use will likely continue to expand in the clinic. The full potential of quantitative MR imaging applied to abdominal imaging has yet to be realized, but the myriad applications reported to date will undoubtedly continue to grow.

Fast bound pool fraction mapping using stimulated echoes

Magnetization transfer imaging advanced to an indispensible tool for investigating white matter changes. Quantitative magnetization transfer imaging methods allow the determination of the bound pool fraction (BPF), which is thought to be directly linked to myelin integrity. Long acquisition times and high specific absorption rates are still inhibiting broad in vivo utilization of currently available BPF mapping techniques. Herewith, a stimulated echoes amplitude modulation-based, single-shot echo planar imaging technique for BPF and T(1) quantification is presented at 3T. It allows whole brain mapping in 10-15 min and is low in specific absorption rates. The method was validated with different concentrations of bovine serum albumin (BSA) phantoms. Intra- and inter-subject variability was assessed in vivo. Phantom measurements verified linearity between bovine serum albumin concentrations and measured BPF, which was independent of T(1) variations. T(1) values in the phantoms correlated well with values provided by standard T(1) mapping methods. Intrasubject variability was minimal and mean regional BPFs of 10 volunteers (e.g., left frontal white matter=0.135 ± 0.003, right frontal white matter=0.129 ± 0.006) were in line with previously published data. Assessment of interhemispheric BPF differences revealed significantly higher BPF for the left brain hemisphere. To sum up, these results suggest the proposed method useful for cross-sectional and longitudinal studies of white matter changes in the human brain.

Fast bound pool fraction imaging of the in vivo rat brain: association with myelin content and validation in the C6 glioma model

Cross-relaxation imaging (CRI) is a quantitative magnetic resonance technique that measures the kinetic parameters of magnetization transfer between protons bound to water and protons bound to macromolecules. In this study, in vivo, four-parameter CRI of normal rat brains (N=5) at 3.0 T was first directly compared to histology. The bound pool fraction, f, was strongly associated with myelin density (Pearson's r=0.99, p<0.001). The correlation persisted in separate analyses of gray matter (GM; r=0.89, p=0.046) and white matter (WM; r=0.97, p=0.029). Subsequently, a new time-efficient approach for solely capturing the whole-brain parametric map of f was proposed, validated with histology, and used to estimate myelin density. Since the described approach for the rapid acquisition of f applied constraints to other CRI parameters, a theoretical analysis of error was performed. Estimates of f in normal and pathologic tissue were expected to have <10% error. A comparison of values for f obtained from the traditional four-parameter fit of CRI data versus the proposed rapid acquisition of f was within this expected margin for in vivo rat brain gliomas (N=4; mean±SE; 3.9±0.2% vs. 4.0±0.2%, respectively). In both whole-brain f maps and myelin density maps, replacement of normal GM and WM by proliferating and invading tumor cells could be readily identified. The rapid, whole-brain acquisition of the bound pool fraction may provide a reliable method for detection of glioma invasion in both GM and WM during animal and human imaging.

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.

A simple correction for B1 field errors in magnetization transfer ratio measurements

B1 errors are a problem in magnetization transfer ratio (MTR) measurements because the MTR value is dependent on the amplitude of the magnetization transfer (MT) pulse. B1 errors can arise from radiofrequency (RF) nonuniformity (caused by the RF coil, or skin effect and dielectric resonance in the subject's head) and also from inaccurate setting of the transmitter output when compensating for varying amounts of loading of the RF coil. B1 errors, and hence MTR errors, may be up to 5-10%, a large source of error in quantitative MR measurements. Radiofrequency nonuniformity may cause MTR histograms to be broadened. The dependence of MTR on B1 was modeled using binary spin bath theory, with a continuous wave (CW) approximation. For B1 reductions of up to 20%, normalized plots for different brain tissue types could be approximated by a single line, indicating that a systematic correction could be applied to MTR measurements with a known B1 error, regardless of tissue type. On a 1.5-T scanner with a birdcage coil, MTR was measured in 18 tissue types in five controls. The MT pulse amplitude was reduced in steps from its nominal value by up to 20%. Averaging data over all controls and tissue types resulted in a line fitting mtr(normalized)=0.812b(1normalized)+0.193, where mtr(normalized) is the normalized value of MTR (relative to its value at the nominal B1) and b(1normalized) is the normalized value of B1 (relative to its nominal value). For a 20% reduction in MT pulse amplitude (i.e., b(1normalized)=0.80), the mean MTR value for the 18 tissue types was 7.0 percent units (pu) below the correct value. After correction using the single equation above for all tissue types, all MTR values were within 1.5 pu of their correct value [root mean square (rms) error=0.7 pu]. Magnetization transfer ratio values tended to be slightly overcorrected because the simple linear correction scheme is only an approximation to the true MTR dependence on B1. A B1 field mapping technique was implemented, based on the double angle method (DAM), with fast spin-echo (FSE) readout, and TR=15 s; this took a total of 6 min of imaging time. This was used to quantify B(1) errors and correct MTR maps and histograms. However, the cerebrospinal fluid (CSF) T1 is very long (approximately 4.2 s); thus, to achieve complete longitudinal relaxation (a requirement of the DAM B1 mapping method), an increase in TR and, hence, acquisition time would be required. In general, however, we are not interested in calculating the B1 in the CSF, although it is important that the B1 is determined in partial volume voxels around the CSF. Using our birdcage head coil, whole-brain B1 histograms were found to have full-width half maximums (FWHMs) ranging from just 6.8% to 11.5% of the nominal B1 value. The FSE DAM B1 field mapping technique was shown to be robust, although a longer TR time may be desirable to ensure complete elimination of CSF partial volume errors. The procedure can be applied on any scanner where the Euro-MT sequence is available, or alternatively, where the amplitude of B1 or of the MT pulse can be manually reduced in order to perform this type of "calibration" experiment for the particular MTR sequence used. The MTR is known to be highly dependent on the parameters of the sequence used, in particular, the MT pulse shape, flip angle, duration, and offset frequency, and the repetition time TR' between successive MT pulses. Therefore, correction schemes will differ for different MTR sequences, and new data sets would be required to calculate these different correction schemes.

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.

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.

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.

Quantitative T(1rho) and magnetization transfer magnetic resonance imaging of acute cerebral ischemia in the rat

It has been previously shown that T1 in the rotating frame (T(1rho)) is a very sensitive and early marker of cerebral ischemia and that, interestingly, it can provide prognostic information about the degree of subsequent neuronal damage. In the present study the authors have quantified T(1rho) together with the rate and other variables of magnetization transfer (MT) associated with spin interactions between the bulk and semisolid macromolecular pools by means of Z spectroscopy, to examine the possible overlap of mechanisms affecting these magnetic resonance imaging contrasts. Substantial prolongation of cerebral T(1rho) was observed minutes after induction of ischemia, this change progressing in a time-dependent manner. Difference Z spectra (contralateral nonischemic minus ischemic brain tissue) showed a significant positive reminder in the time points from 0.5 to 3 hours after induction of ischemia, the polarity of this change reversing by 24 hours. Detailed analysis of the MT variables showed that the initial Z spectral changes were due to concerted increase in the maximal MT (+3%) and amount of MT (+4%). Interestingly, the MT rates derived either from the entire frequency range of Z spectra or the time constant for the first-order forward exchange (k(sat)) were unchanged at this time, these variables reducing only one day after induction of ischemia. The authors conclude that T(1rho) changes in the acute phase of ischemia coincide with both elevated maximal MT and amount of MT. These changes occur independent of the overall MT rate and in the absence of net water gain to the tissue, whereas in the consolidating infarction the decrease in the rate and amount of MT, as well as the extensive prolongation of T(1rho), are associated with water accumulation.

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.

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.