|
|
||||||||
Technical Developments |
1 From the Departments of Medical Biophysics (S.J.G., R.M.H.) and Medical Imaging (R.M.H.), University of Toronto, Ontario, Canada, and Imaging and Bioengineering Research, Sunnybrook and Women's College Health Sciences Centre, 2075 Bayview Ave, Toronto, Ontario, Canada M4N 3M5. Received September 10, 1998; revision requested November 3; revision received December 2; accepted March 1, 1999. Supported in part by the Terry Fox Foundation and GE Medical Systems of Canada. Address reprint requests to S.J.G. (e-mail: sgraham@sten.sunnybrook.utoronto).
| Abstract |
|---|
|
|
|---|
Index terms: Brain, abnormalities, 10.871 Brain, MR, 10.121417 Brain, white matter, 10.121417 Magnetic resonance (MR), magnetization transfer contrast, 10.121417 Sclerosis, multiple, 10.871
| Introduction |
|---|
|
|
|---|
Qualitatively, the MT effect is easily described. When radio-frequency (RF) energy is applied several kilohertz off resonance from the Larmor frequency of the water molecules in tissues, longitudinal magnetization associated with tissue proton macromolecules (the semisolid pool) is preferentially saturated in comparison with water protons (the liquid pool), which normally provide the dominant signal intensity in clinical MR imaging. Exchange of longitudinal magnetization between the pools, via macromolecular sites such as hydroxyl or amide groups (9,10), causes signal reduction in the liquid pool. For a variety of tissues, this MT effect is typically 10%50% of the signal intensity obtained at clinical imaging when MT pulses are not applied (11). Quantitatively, however, the MT effect is a complicated function of the relaxation times and absorption line shapes of the semisolid and liquid pools, their respective volume fractions, and the efficiency of MT exchange. Notwithstanding this complexity, quantitative understanding of MT provides important physical insight regarding how MT image contrast is generated for specific pulse sequences. From the fundamental MT and relaxation properties of tissues, the optimal contrast can be determined for specific clinical applications. Early in the development of MR imaging, similar analyses were performed to optimize T1- and T2-weighted contrast (12,13).
In the past decade, progress has been made in describing the MT properties of tissues mathematically. Initially, a two-pool model was developed to characterize the MT effect in an aqueous agar gel system under continuous RF irradiation off resonance (14). This model also applies for many clinically relevant tissues (gray matter, white matter, muscle, liver, and blood) when irradiation of the semisolid pool is correctly specified (11,15). More recently, the model was extended to include saturation pulses of short duration, which are used in clinical MT imaging to minimize patient power absorption and imaging time (16). In this study, we used this method to evaluate MT imaging protocols for two applications: static signal suppression in cerebral MR angiography and optimization of contrast in white matter disease. In particular, there is need for such an evaluation in white matter neuroimaging; the reported MT effect for normal white matter ranges from 16% (17) to 60% (18), which suggests that there are wide variations in the effectiveness of clinical MT protocols.
Simulations were performed for spoiled gradient-echo (GRE) imaging, a pulse sequence often used in conjunction with MT. The predicted signals, based on measurements of the MT and relaxation properties of excised tissues at body temperature, were compared with those in the literature.
| Materials and Methods |
|---|
|
|
|---|
Experiments were conducted at 1.5 T by using a superconducting magnet with a bore diameter of 20 cm (Nalorac Cryogenics, Martinez, Calif), controlled with a programmable console (SMIS, Surrey, England). Samples were maintained at body temperature with heated forced air. Experiments were performed without imaging to measure bulk sample properties. The rectangular pulse duration was approximately 10 µsec for a 90° pulse. Three MR experiments were performed: (a) a Carr-Purcell-Meiboom-Gill, or CPMG, sequence (21,22) to measure T2 (repetition time [TR] of 10,000 msec, 1-msec echo spacing, 5,000 even echoes sampled, four signals averaged); (b) an inversion-recovery sequence (21) to measure T1 relaxation (31 inversion time intervals spaced logarithmically from 1 to 32,000 msec, free induction decay amplitude acquired at each inversion time interval, two signals averaged, 10,000 msec between 90° and inversion pulses); and (c) continuous-wave RF irradiation off resonance (14) to measure the MT effect (5-second duration, 27 frequencies off resonance spaced logarithmically from 0.01 to 225 kHz, seven continuous-wave irradiation amplitudes increasing from pulse amplitude, or f1, of 83 Hz by factors of two). Because the duration of the MT experiments was approximately 84 minutes, the total set of MT data for blood was obtained by measuring seven well-mixed blood samples for 12 minutes each to eliminate potential error due to sedimentation of red blood cells. Measurements of the sample temperature indicated that negligible RF power absorption occurred during the course of these experiments (data not shown).
Analysis of these data first involved estimation of average relaxation times, or <T1> and <T2>, for each tissue by means of monoexponential fits to the relaxation data. For white matter, average T2 was estimated by means of a weighted sum of relaxation times obtained from a biexponential fit. MT parameters and relaxation times for the liquid pool A, T1A and T2A, were estimated from a least squares fit of the MT data to an established two-pool model describing the MT effect for continuous-wave RF irradiation (14) and knowledge of the average T1 and T2 values. An iterative procedure was included to estimate the absorption line shape of the semisolid pool (15), and it was performed at each irradiation frequency off resonance. (T1 for the semisolid pool B, T1B, was not estimated precisely from the data and was assumed to be 1,000 msec for all tissues. This assumption is justified because the T1B value has little influence on pulsed MT simulations [16]. Sensitivity analysis [data not shown] indicated that a 10% change in T1B produced a change of approximately 0.1% in MT effect.) All errors were reported either as the average SD obtained from the fitting procedure or as half the difference between the mean parameter estimates for the two samples of each tissue, whichever was largest.
Pulsed MT with Excitation, Absorption, Saturation, and Evolution, or "EASE"
Simulation of pulsed MT contrast requires an understanding of four key concepts: excitation, absorption, saturation, and evolution. These concepts are illustrated in Figure 1, and mathematic treatment is provided elsewhere (16). In response to an MT pulse, magnetization in the liquid pool, MA, nutates about the associated effective B1 field in a manner very similar to the well-known RF excitation described with the Bloch equations (Fig 1, A). The associated decrease in longitudinal magnetization can be mistaken with MT and is known as the direct effect. Absorption of RF energy by both pools, however, is characterized by the absorption line shapes, gA and gB, and the MT pulse amplitude, f1, applied at a given frequency offset from resonance,
(Fig 1, B). Because the semisolid pool has a very small T2 value (approximately 10 µsec), the line shape gB is much broader than that of gA, and absorption occurs preferentially in the semisolid pool and becomes less efficient (requires more power) as the offset frequency increases. Fast T2 relaxation in the semisolid pool also ensures that RF absorption manifests as saturation (attenuation) of longitudinal magnetization, MBz, with negligible generation of transverse magnetization (Fig 1, C). When the MT pulse is turned off, magnetization in both pools evolves owing to T1 recovery and MT exchange, characterized by a rate constant, R (Fig 1, D). The equilibrium magnetization in the liquid pool is taken as unity, whereas the relative size of the semisolid pool is indicated as MBo. The evolution occurs in two phases: a fast phase, dominated by the rate R, in which the liquid pool replenishes the longitudinal magnetization in the semisolid pool quite rapidly, and a slow phase, dominated by T1A, in which the coupled system returns to equilibrium.
|
Pulsed MT Imaging
Although the excitation, absorption, saturation, and evolution principles are applicable to arbitrary pulse sequences, MT pulses used with spoiled GRE imaging are of primary interest because this combination currently receives the most clinical application. The pulse sequence used in simulations in this study contained three distinct elements (Fig 2): MT pulse plus spoiler gradient (with total time duration, T); small flip angle excitation characterized by the angle,
, which produces transverse magnetization for the liquid pool but has negligible effect on the semisolid pool; and TR, which incorporates time for spatial encoding and evolution of image contrast. (The inclusion of MT pulses inevitably compromises the minimum TR achievable with spoiled GRE imaging). Transverse magnetization is assumed to be completely spoiled prior to application of the subsequent MT saturation pulse.
|
, and the conventional spoiled GRE parameters,
and TR. To investigate the relative importance of each of these parameters and the associations between them, simulations were performed for each clinical application by using a broad range of possible imaging protocols (f1, 0.21.0 kHz;
, 0.54.0 kHz;
, 5°40°; TR, 20125 msec). The spoiled GRE signal intensity without MT pulses, SI, was also calculated according to the appropriate equation (24):
Simulations of Cerebral MR Angiography
The pulse sequence shown in Figure 2 is easily used to simulate imaging contrast for cerebral MR angiography. Current state-of-the-art clinical imaging protocols combine pulsed MT, multiple offset thin-slab acquisition, or MOTSA (25); application of the three-dimensional fast spoiled GRE sequence; and tilted, optimized, nonsaturating, excitation, or TONE, pulses for ramped small flip angle excitation across each slab (26). A typical slab thickness of 20 mm was assumed, and the ratio of
values at the entrance and exit of the slab was 1:2.2 according to clinical prescriptions (27). Time-of-flight contrast between flowing blood and static tissues was assessed for vessels in the center of the slab. While magnetization associated with static tissues was allowed to reach steady state, blood was assumed to flow uniformly perpendicular through the slab and thus was subject to markedly fewer MT pulses and ramped RF excitation, which is dependent on blood velocity and TR (25). To assess image contrast for small 0.5-mm-diameter vessels, an empirical model relating vessel diameter and blood velocity was used (28) that predicted an average flow velocity of 37 mm/sec.
Because MR angiographic applications are subject to imaging time constraints, simulations were performed with TRs of 45 msec and 26 msec without MT pulses for comparison. These conditions represent the minimum TR values obtainable in a previous clinical study (27). Blood at the center of the image slab therefore received either six or 10 excitations with or without MT, respectively. In the absence of MT, over the middle half of the slab, the blood signal varied plus or minus 20% from the value at the center, assuming a flip angle of 25°.
To compare the efficiency of cerebral MR angiography with different MT pulse parameters and without MT, the signal difference between blood and static tissue was calculated and divided by the square root of TR. This figure of merit is proportional to the signal-differenceto-noise ratio independent of acquisition time. The calculation accounts for the increase in minimum imaging time that occurs when MT pulses are used, which could in principle be used to perform signal averaging when MT pulses are absent.
Simulations of White Matter Neuroimaging
In this example, simulations were performed with a TR value of 20 msec, which approximates the minimum TR that is currently achievable when MT pulses are included, and also for a TR value of 125 msec, a typical value used for MT imaging of multiple sclerosis (5). Results of imaging protocols were compared on the basis of the MT ratio difference between white and gray matters per square root of TR, in which the MT ratio was expressed as the percentage signal decrease due to MT relative to the signal achieved without MT at the same TR value.
| Results |
|---|
|
|
|---|
|
. The area under each line shape equals one and for clarity, error bars (obtained by varying
2 for the fit by ±1 SD) are indicated for white matter only. The absorption line shape for blood is much narrower than those for gray and white matters, displaying a fivefold larger amplitude on resonance (
= 0) and a half width at half maximum of approximately 500 Hz. The line shapes for white and gray matters are equivalent within experimental error and still have appreciable amplitude at 10 kHz off resonance (data not shown).
|
Mxy (blood signal minus static tissue signal) per square root of time is shown in Figure 4 for a range of MT imaging protocols at TRs of 45 msec (thick lines) and for protocols of 26 msec (thin lines) without MT pulses. For clarity, only MT protocols that were quite close to maximizing the signal difference per square root of time are shown. White matter was taken as the static tissue because it has a much shorter T1 than does gray matter and is more difficult to suppress. Figure 4, top, which illustrates the relationship with flip angle, indicates that use of an MT pulse of modest amplitude (f1 = 450 Hz,
= 1.2 kHz) resulted in an improvement of 11% in signal difference per square root of time compared with that for imaging without an MT pulse. Because the MT pulse increases the available minimum TR (in this case, 45 vs 26 msec), the imaging flip angle must be increased for maximum efficiency, from 25° without MT pulses to 30° with MT pulses. The gain in efficiency is only weakly dependent on MT pulse amplitude and off-resonance frequency (Fig 4, middle and bottom, thick lines), requiring f1 and
values above thresholds of only approximately 300 Hz and 1.5 kHz, respectively.
|
Gray MattertoWhite Matter Contrast
A similar family of curves is shown in Figure 5, which illustrates the MT ratio difference between white and gray matters per square root of time for different MT protocols with TRs of 20 and 125 msec. Again, only protocols producing MT ratio differences that were close to the most efficient at each TR value are shown. The maximum MT ratio difference was slightly less at longer TR (16% for TR of 20 msec, with respect to gray matter, 11% for TR of 125 msec [data not shown]). Expressed per square root of time, therefore, the MT ratio difference generated at TR of 20 msec was more efficient by a factor of approximately 3.6. The flip angle value for maximum MT was smaller at TR of 20 msec (
= 15°) and required adjustment within 5° for best efficiency, whereas
values for TR of 125 msec provided good sensitivity from 20° to 40° (Fig 5, top). A similar effect was exhibited with respect to MT pulse amplitude and frequency off resonance (Fig 5, middle and bottom). The maximum MT ratio difference at TR of 20 msec required f1 amplitudes of ap- proximately 300500 Hz at 1.52.5 kHz off resonance, whereas TR of 125 msec required f1 greater than approximately 500 Hz at 1.54.0 kHz off resonance.
|
| Discussion |
|---|
|
|
|---|
Cerebral MR Angiography
Previously, Goodrich et al (27) compared results with multiple offset thin-slab acquisition, or MOTSA, plus tilted, optimized, nonsaturating, excitation, or TONE, methods with MT pulses to those without MT pulses, with use of the same TR values investigated in our study and the same MT pulse envelope. For fixed MT pulse amplitude and frequency off resonance (f1 = 470 Hz,
= 1.2 kHz), various tilted, optimized, nonsaturating, excitation pulse flip angles and slab thicknesses were investigated. For 10 patients, the largest peak signal-differenceto-noise ratio with MT for 0.5- to 0.75-mm-diameter vessels occurred with a slab thickness of 16 mm and represented a gain of 29% compared to results with analogous imaging without MT pulses. When expressed as signal-differenceto-noise ratio per unit time, these two imaging protocols were in fact roughly equivalent in efficiency. In comparison, the simulations predicted 45% increased signal difference when MT pulses were used compared with when MT pulses were absent and predicted a modest gain in signal difference per unit time of 11%. Considering the inherent assumptions in the simulations, the simulation and clinical results were in quite good agreement. Furthermore, the simulation results indicate that vessel contrast is quite insensitive to the choice of MT protocol and suggest that appropriate technique is currently used in this clinical application (27). MT pulse amplitudes above a threshold of 400 Hz and offset frequencies above 1.5 kHz all prove satisfactory.
Imaging of White Matter Disease
McGowan et al (32) developed an MT imaging protocol that has been adopted in numerous neuroimaging studies (3,5,33). In initial experiments, sinc pulses were used, with nominal bandwidth of 100 Hz and duration of 19 msec; from investigations of several different f1 amplitudes, f1 of 160 Hz was found to produce the largest MT effect for frequency offsets ranging from 3 to 15 kHz. This MT pulse was coupled with GRE imaging (TR of 100 msec, minimum echo time,
of approximately 10°) to achieve primarily proton-density contrast in the absence of MT pulses. The average specific absorption rate of this pulse sequence was approximately 2.0 W/kg and generated an average MT ratio for normal white matter of 43% ± 1 (mean ± SD) (33). Alternatively, Finelli et al (18) proposed a technique with a 10-msec single-period sinusoid envelope and f1 of 700 Hz and
of 4 kHz coupled with a more rapid spoiled GRE sequence (TR of 25 msec, echo time of 4.4 msec, and flip angle of 6°). In this case, the average specific absorption rate was 4.0 W/kg, which produced an MT ratio of 59% ± 2 in normal white matter.
The large difference in MT ratios obtained with these two different imaging protocols can be easily explained theoretically. For a single ideal MT pulse that completely saturates the semisolid pool and leaves the liquid pool unaffected, the critical time, or tcrit, can be calculated when maximum signal decrease occurs in the liquid pool owing to MT exchange (16). From the Table, the critical time is approximately 100 msec for blood and approximately 70 msec for gray and white matters. The MT pulse repetition period should, therefore, be much less than the critical time to accumulate the MT effect efficiently over many pulses. Thus, imaging with shorter TR interval produces larger MT ratio values. If TR is much greater than the critical time, then MT effect is accumulated inefficiently; pulse power must be markedly increased to produce more direct effect. It is only because T1 recovery in the liquid pool occurs more rapidly than MT exchange under these conditions that appreciable MT ratio values are obtained at longer TR (16).
The simulation results suggest an alternative MT imaging protocol for neuroimaging to that previously advocated. Both McGowan et al (32) and Finelli et al (18) focused on maximizing the MT effect for white matter. However, if white matter MT is maximized, then it may also be maximized for white matter that is partially demyelinated or exhibits partial axonal loss. What is ultimately important is the MT ratio contrast between normal white matter and white matter lesions. With use of gray matter to represent axonal loss, we found that MT ratio contrast is maximized between white and gray matters with technique that does not provide the maximum MT effect in both tissues. Furthermore, this protocol, found by means of simulation (f1 = 400 Hz,
= 2 kHz,
= 15°, TR = 20 msec), exhibits an average specific absorption rate less than 4.0 W/kg when compared with that of Finelli et al (18).
General Recommendations
The results from these simulations were obtained on the basis of a number of simplifying assumptions. In particular, (a) MR properties of ex vivo tissues at body temperature were used to approximate those of living tissues in physiologic conditions, with bovine gray and white matters substituted for the human analogues; (b) T2* effects were neglected; (c) owing to the lack of detailed MR measurements of the MR properties of white matter lesions (which are variable in character), gray matter was used as a substitute; and (d) time-of-flight signal enhancement for blood accounted neither for vessel narrowing, branching, and tortuosity nor for the nonlocalized MT saturation in clinical application (27) (ie, where each MT pulse actually saturates the entire head, not just a single slab). The error introduced by these assumptions is difficult to assess. Assumption d, however, means that flowing blood is likely subject to more MT and imaging pulses than were simulated. This could account for the theoretically greater MT over that actually observed clinically (27).
Nevertheless, the simulation results were sufficiently compelling that they highlight the need to optimize MT imaging protocols by means of clinical studies and to investigate image contrast further in patient populations. For cerebral MR angiography, results in such studies have produced good MT imaging technique (27). There are several basic guidelines for conducting such investigations that are easily implemented in any MR imaging laboratory.
1. Choose the minimum TR that allows an acceptable specific absorption rate. Use of MT pulses alters the flip angle
that produces best image contrast per unit time, so image quality must be investigated as a function of
. This can be performed adequately with use of the standard MT pulse provided with the MR imager. Noise levels should be reduced to acceptable levels by means of signal averaging.
2. Choose an MT pulse: Both binomial on-resonance and shaped off-resonance pulses are acceptable. Findings in previous simulations and clinical studies indicate that appreciable MT effect can be obtained by means of both binomial on-resonance and shaped off-resonance pulses (16,18,34). The direct effect must be minimized in both cases, which is simpler with shaped pulses because they are less sensitive to static field inhomogeneity (34). For different shaped pulses, MT effect scales linearly with average pulse power (16,35); standard choices, such as sinc or Gaussian pulse envelopes should provide desired MT contrast within current specific absorption rate limits. To minimize the required power, the off-resonance frequency should be chosen just large enough to produce negligible direct effect, and results of simulations indicate that offset frequencies of 1.52.0 kHz are appropriate.
3. Investigate image contrast as a function of MT pulse amplitude. Weak MT pulses produce poor MT image contrast, whereas excess power produces undesirable direct effect. For neuroimaging, findings in this study suggest that the desirable MT pulse amplitude is insensitive to the imaging flip angle, so that
determined in guideline 1 should be sufficient. If MT pulse amplitude must be increased, then the TR value may also have to be increased slightly for specific absorption rate to remain acceptable, and guideline 1 should be revisited.
The simulation results also suggest that the MT protocol for cerebral MR angiography is flexible (a range of MT pulse parameters produce equivalent contrast). It may be possible to determine the best MT protocol for imaging white matter disease and then to use this protocol for cerebral MR angiography without much loss of efficiency. Such an approach would facilitate pulse sequence implementation but requires verification with further study.
Other Pulse Sequences
This study focused on use of MT pulses with spoiled GRE imaging, the predominant pulse sequence used in MT protocols. However, appreciable deliberate MT effect can also be generated with spin-echo imaging (34) as long as enough interleaved sections are acquired per TR interval to ensure that MT pulses are applied with a period much less than the critical time. This pulse sequence is more complicated to simulate because some saturation of the semisolid pool is caused by refocusing pulses at varying frequencies off resonance during multisection interleaving (36). An accurate comparison of pulsed MT with spoiled GRE versus spin-echo pulse sequences is beyond the scope of this article. Interestingly, use of MT pulses with fast spin-echo imaging produces little deliberate (17) MT effect, because the echo train length typically exceeds the critical time. However, unintended MT contrast is also produced when fast spin-echo imaging is performed without an MT saturation pulse, owing to the many refocusing pulses (36,37). This suggests that calculation of MT ratios is inadvisable with use of fast spin-echo imaging.
In conclusion, simulation results indicate that current use of MT in imaging protocols is close to optimal for cerebral MR angiography, whereas neuroimaging protocols can likely be made more effective. It is hoped that this report will encourage clinical neuroradiologists to investigate improvement of MT protocols for white matter lesions.
| Footnotes |
|---|
Author contributions: Guarantor of integrity of entire study, R.M.H.; study concepts and design, R.M.H., S.J.G.; definition of intellectual content, R.M.H., S.J.G.; literature research, S.J.G., R.M.H.; experimental studies, S.J.G.; data acquisition and analysis, S.J.G.; statistical analysis, S.J.G.; manuscript preparation and editing, S.J.G.; manuscript review, R.M.H.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |