Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Magn Reson Med. Author manuscript; available in PMC 2010 June 1.
Published in final edited form as:
PMCID: PMC2860191

WAter Saturation Shift Referencing (WASSR) for chemical exchange saturation transfer experiments


Chemical exchange saturation transfer (CEST) is a contrast mechanism exploiting exchange-based magnetization transfer (MT) between solute and water protons. CEST effects compete with direct water saturation and conventional MT processes and generally can only be quantified through an asymmetry analysis of the water saturation spectrum (Z-spectrum) with respect to the water frequency, a process that is exquisitely sensitive to magnetic field inhomogeneities. Here, it is shown that direct water saturation imaging allows measurement of the absolute water frequency in each voxel, allowing proper centering of Z-spectra on a voxel-by-voxel basis independent of spatial B0 field variations. Optimal acquisition parameters for this “water saturation shift referencing” or “WASSR” approach were estimated using Monte Carlo simulations and later confirmed experimentally. The optimal ratio of the WASSR sweep width to the linewidth of the direct saturation curve was found to be 3.3–4.0, requiring a sampling of 16–32 points. The frequency error was smaller than 1 Hz at signal to noise ratios of 40 or higher. The WASSR method was applied to study glycogen, where the chemical shift difference between the hydroxyl (OH) protons and bulk water protons at 3T is so small (0.75–1.25 ppm) that the CEST spectrum is inconclusive without proper referencing.

Keywords: CEST, absolute water frequency referencing, direct water saturation, glycogen, muscle


Chemical exchange saturation transfer (CEST) has recently emerged as an alternative contrast mechanism for MRI (19). In this approach, selective radiofrequency (RF) irradiation of exchangeable solute protons, such as amide (NH) and hydroxyl (OH), groups (1013), is detected through progressive saturation of the water signal consequential to chemical exchange. Chemical exchange is a well-known magnetization transfer (MT) mechanism with a single transfer pathway. As such CEST should be distinguished from conventional magnetization transfer contrast (MTC) imaging based on semi-solid protons, where saturation transfer occurs in multiple steps, both through dipolar coupling and chemical exchange. Such indirect detection through magnetization transfer (MT) processes allows detection of species normally invisible in MRI with an effect size of a few percent on the water signal, i.e. sensitivity enhancement to the molar concentration range. For instance, CEST observation of solutes and particles in the millimolar to nanomolar range has been demonstrated both in vitro (3,6,1417) and in vivo (10,18,19). Similar to MTC studies, CEST effects can be detected through irradiation of the proton spectrum as a function of frequency offset and measurement of the ratio of the water signal with saturation and without. A simulated example of such a water saturation spectrum, more commonly called a Z-spectrum (20) or CEST spectrum (2), is given in Fig. 1A, where a solute (amide) proton at 3.5 ppm from water shows a clear CEST effect. The large saturation effect around 0 ppm is due to direct water saturation, with the water frequency assigned to 0 ppm for reference.

Figure 1
Simulated Z-spectra (dashed lines) and MTRasym plots (solid lines) for an amide proton (A, C, E) and a glycogen OH proton (B, D, F)

CEST effects can be quite small and, if the resonance frequency of the solute protons is sufficiently close to the water frequency, the CEST effects have to compete with direct saturation (Fig. 1B). In vivo the situation is even worse, due to the occurrence of MTC effects. In order to minimize these interferences, CEST imaging generally requires an asymmetry analysis of the MT effects with respect to the water frequency, providing an MT asymmetry ratio (MTRasym) spectrum. This parameter is defined as


in which Δω is the shift difference between the irradiation frequency and the the water center frequency and Sω) and S0 are the water intensities after a long presaturation pulse at the offset frequency and without a presaturation pulse, respectively. Examples of such MTRasym spectra for typical examples of the NH group in proteins or the OH-group in glycogen are demonstrated by the solid curves in Figs. 1A–D. It can be seen that when the direct saturation curve is quite narrow (Figs. 1A and B), an asymmetry analysis works quite well. When the line broadens, the distinction of the CEST effect from direct saturation becomes problematic if the chemical shift difference between the solute proton and water frequency (Δω) is relatively small, as is the case for the hydroxyl (OH) protons of glycogen (Δω = 0.75–1.25 ppm, Fig. 1D). However, if the water frequency is known exactly, an asymmetry analysis still provides the correct CEST effect. Unfortunately, this seemingly simple procedure renders CEST imaging extremely sensitive to magnetic field inhomogeneities, as illustrated in Figs. 1E and F. Due to the steep slope of the direct saturation curve, even a small B0 field difference and a concomitant shift in the Z-spectrum may cause a large change in MTRasym. In vivo, where local magnetic susceptibility differences are commonplace, this effect will lead to the occurrence of artifactual signal spikes in CEST images (18,21). Depending on the solute studied, this may lead to simple errors in quantification of the magnitude of the effect, as for the amide proton in Fig. 1E, or sometimes even to total removal of the CEST effect, such as for glycogen in Fig. 1F. In addition to improving shimming, which is not likely to be sufficient in vivo, another approach to reduce this problem is to acquire images over a range of saturation frequency offsets (Z-images, typically with high saturation power), followed by polynomial or cubic-spline fitting and centering of the Z-spectrum in each voxel (10,12). This may work in vitro or even in vivo if the CEST and water saturation curves are sufficiently separated, because the direct saturation curve is symmetric and the center can be easily assigned, even if the bottom of the water saturation curve is flat and wide. However, if the CEST peak overlaps with the water resonance, as in the case for OH groups in sugars (Figs. 1D and F) (19) this approach to centering is not possible. In addition, if studies are done in vivo where the conventional MT effect has been shown to be asymmetric (10,22,23), such an approach may be less accurate. At first glance it seems that traditional field mapping (2428) should be able to solve this issue, but the need for at least one voxel where the exact water frequency is known prohibits application of this approach.

We here propose that absolute water frequency referencing can be accomplished by acquiring a pure direct water saturation image. This can be done by using RF irradiation that is of sufficiently small power and duration to have negligible interference of magnetization transfer effects (both MTC and CEST). To test this water saturation shift referencing (WASSR) approach, we first performed simulations to evaluate the experimental parameters needed to accurately measure the absolute water frequency in each image voxel, namely the WASSR sweep width, the saturation line width, and the number of points needed to sufficiently sample the frequency at different signal to noise ratios (SNRs). To test its feasibility, we applied WASSR to the CEST measurement of glycogen (glycoCEST, (19)), where the OH proton CEST peak cannot be visually separated from direct saturation at a field strength of 3T. This was done both in vitro and in vivo.


Theory: WASSR with maximum–symmetry algorithm

We consider the WASSR Z-spectrum, which isolates the effect of direct water saturation. Note that the shape of the direct water saturation Z-spectrum is not affected by field inhomogeneities making it symmetric with respect to its center frequency (point of minimum intensity). Thus, the center of the direct saturation curve can be found by equalizing the frequency differences between its rising and falling slopes using a symmetry analysis. It is important to realize that the position of this symmetric line is affected by the local magnetic field, which varies between voxels due to differences in magnetic susceptibility. To determine the offset in each voxel, we used the procedure outlined in Fig. 2B. The observed signal, f(x), is reflected (mirror image) with respect to the experimentally found water frequency. We then minimize the mean squared error between the measured (original) intensities and the cubic-spline interpolated intensities of the reflected curve with the following discrete energy function [2],

Figure 2
Data processing procedures


in which xi are the sampled WASSR frequency locations among N number of points (i = 1…N) and f(2Cxi) is the reflected and interpolated signal of f(xi) when C is the estimated center frequency of the original curve. The maximum-symmetry center frequency (MSCF) estimator executed argument of minimum (argmin) over a sum (< >) of mean squared error and numerical optimization was performed with Nelder-Mead simplex direct search method (29). Search initialization was achieved by selecting the median frequency of points that are at half the median intensity height. Other initializations are possible; these heuristics were chosen for convenience. Note that the mean estimated power is taken only over the points that are sampled and that, therefore, the estimate of C is restricted to the lie within the range of frequencies that are sampled. Symmetry analysis is not possible without acquiring at least some samples around the point of symmetry, which therefore requires acquisition of a minimum number of points depending on the width of the saturation curve and the chosen spectral width of the Z-spectrum. The latter depends on the range of frequency shifts over the sample.


Monte Carlo (MC) simulations were performed to estimate the optimal parameters for WASSR imaging, including the WASSR Sweep Width (WSW) and the number of WASSR points (WN) needed to properly cover the WSW. To determine the WSW and the minimum WN, it is necessary to know the width and shape of the direct saturation curve, which was simulated using the exact steady state solution for water magnetization during saturation as provided by Mulkern and Williams (30). In terms of signal intensity, this gives:


in which R1= 1/T1 (longitudinal relaxation), R2=1/T2 (transverse relaxation), ω1=γB1 and Δω=ωω0. Δω denotes the offset of the signal point (x) with respect to the center frequency. We investigated the optimal WSW in terms of a ratio with respect to the WASSR Line Width (WLW) as defined by the Full-Width at Half Minimum of the saturation lineshape. Calculation of the line width of a lineshape such as Eq. [3] is normally done by determining the frequency offset (Δω) at half height (S/S0 = 0.5 for full saturation). However, for the direct saturation lineshapes used here, this is not straightforward because the saturation is generally not 100% at Δω = 0. We therefore first calculated the actual signal intensity at S(Δω =0)/S0, subtracted this from the reference intensity (S/S0 = 1) and determined the frequency at half height from that. The result is:


The WASSR points were applied equally spaced within the WSW as this is most practical when the field varies over a sample and shifts the saturation curve over multiple frequencies. In WASSR, the center frequency of the saturation curve needs to be found using the maximum-symmetry algorithm described above. We simulated the absolute error in this frequency (offset error) for different ratios WSW/WLW as a function of WN and SNR. The noise used was Rician and SNR was taken for water signal without saturation, defined as SNR(S0). To estimate the frequency shift, we modeled the center frequencies as a zero-mean normal distribution (N) with standard deviation of WLW at each spectrum as follows:


For each WSW/WLW ratio and WN, the expected mean absolute error was estimated with 1275 MC iterations at a simulated SNR(S0) of 100:1. The set of pairs of best WSW/WLW ratios (i.e. minimum offset error) at a certain WN was fitted with a logarithmic form to generate an expression for the optimal ratio when using a certain number of points:


where A and B are the fitting parameters. Using this expression, the expected WASSR offset error was computed as a function of WN for SNR(S0) values of 10:1, 20:1, 40:1, 80:1, 100:1, and 1000:1. Again, each expected WASSR error computation employed 1275 MC iterations. To determine the variability of parameter fit, these MC experiments were repeated 25 times.


A 200mM solution of bovine liver glycogen (Type 1X G0885-25G; Sigma, St. Louis, MO) with an average molecular mass of about 50 kD was prepared by dissolving 4.98g of glycogen in 151 ml of phosphate buffered saline (pH 7.3). The glycogen concentration is expressed in glycosyl units, with each glucosyl unit contributing 168 g/mol [(4980 mg/168 mg/mmol)/151 ml) = 0.2 M].

Human subjects

Human studies were approved by the Johns Hopkins Medicine Institutional Review Board (IRB) and the Kennedy Krieger Institute, and all subjects provided signed consent. Four healthy subjects (two females and two males, 33 ± 10 years and 68 ± 9 kg, Mean ± STD) were kept in the resting state about 0.5 to 1 hour to maintain the baseline before they were scanned to measure the glycogen CEST effect in calf muscle.

MRI Acquisition

All images were acquired using a whole-body Philips 3T Achieva scanner (Philips Medical System, Best, The Netherlands) equipped with 80 mT/m gradients. RF was transmitted using the body coil and SENSE reception (31) was employed. A series of consecutive direct saturation and CEST scans were performed using the 8-element knee coil for both the glycogen phantom and in vivo human calf muscle. To minimize leg motion, foam padding was placed between the subject’s lower leg and the knee coil. In all cases, second order shims over the entire muscle on the imaging slice were optimized to minimize B0 field inhomogeneity. Notice that the width of Z-spectra depends on T2 and that the WASSR procedure provides an absolute field frequency map so that there is no need for higher order shimming for the CEST acquisition. Clinical imagers generally employ a prescan to center the bulk water signal of the object/subject, optimize the flip angle and shim the field. Note that no such “prescan” should be made between direct saturation and CEST scans to maintain the same field reference conditions. For both scans, saturation was accomplished using a rectangular RF pulse before the turbo spin echo (TSE) image acquisition, as previously described by Jones et al. (21).

The power level needed for each saturation experiment depended on the load and was optimized by measuring sets of Z-spectra under these different conditions. For WASSR, the power and pulse lengths were chosen as small as possible to have sufficient direct saturation, while minimizing any MT effects. For CEST, the maximum pulse length allowed for the body coil within the protected clinical software (500 ms) was used and the power was optimized for maximum effect at the phantom and muscle loads. WASSR was obtained at higher frequency resolution than CEST, but over a smaller frequency range as only the direct saturation region needs to be covered. The WASSR range was chosen sufficiently large to validate the simulated results, consequently leading to a larger number of frequencies needed in vivo than for the phantom.

Single-slice glycogen phantom imaging was performed using SENSE factor = 2, TSE factor [number of refocusing pulses] = 34 (two-shots TSE), TR = 3000 ms, TE = 11 ms, matrix = 128 × 122, FOV = 100 × 100 mm2, slice thickness = 5 mm, NSA = 1. Imaging parameters for human calf muscle experiments were identical to those in phantom experiments except for the following: FOV = 160 × 160 mm2. The saturation spectral parameters for WASSR and CEST are indicated in Table 1.

Table 1
WASSR and CEST imaging parameters utilized in this study.

Data analysis

Data were processed using a custom-written program in Matlab (Mathworks, Natick, MA, USA). For efficiency of the processing, the background around the object was masked by excluding data points below the 50th percentile of voxel intensity. The data processing can be divided into reconstruction of a WASSR shift map and a CEST map, as shown in Fig. 2A. To create a WASSR shift map, the WASSR spectrum from each voxel was cubic-spline interpolated and the center frequency was optimized by the maximum–symmetry algorithm followed by search initialization of the center frequency (see details in Theory). The frequency difference between optimized center frequency and 0 ppm (experimental water frequency) for each voxel was stored to the WASSR shift map. To create a CEST map, the CEST spectrum in each voxel was interpolated using fast fourier transformation (FFT). First, 33 original raw CEST data points in the range of −4 to 4 ppm (frequency-domain) were transformed to the time-domain. Second, the transformed CEST data in time-domain were processed by 2 Hz line-broadening and zero-padding. Third, the processed CEST data in time-domain were inverse-transformed to frequency-domain. Overall, 33 points were interpolated to 2048 points. Subsequently, the frequency difference from the WASSR shift map was applied to correct the CEST Z-spectrum on a voxel-by-voxel basis. Finally, the MTRasym intensity computed by Eq. [1] was integrated over the 0–2 ppm range on a voxel-by-voxel basis to produce the CEST map. CEST data in different regions of interest (ROIs) were quantified for 4 subjects. Four different regions of 35 ± 0.8 voxels (mean ± STD) were chosen for calf muscle in each subject.



In Fig. 3A, simulations of the error in the determined center frequency of the direct saturation curve (WASSR shift map frequency) are shown as a function of the sweepwidth-to-linewidth ratio (WSW/WLW) under different sampling conditions at SNR(S0) = 100:1. It can be observed that 12 WASSR points within 3.25 times WSW/LW result already in sub-Hz error. The relationship between the predicted optimal ratio (WSW/WLW at minimized mean absolute error) and number of points (BOLD line in Fig. 3A) was fitted to be:

Figure 3
Simulations of offset error (logarithmic scale) in the WASSR center frequency

When plotting this optimal ratio at different WASSR resolution as a function of SNR (Fig. 3B), it is found that the estimated absolute offset error already falls into sub-Hz range when using 16 points at a SNR(S0) of 40:1 and smaller at higher SNR. Considering that the SNR(S0) for typical WASSR data is around 100:1 for humans and over 1000:1 for phantoms, these results imply that using a WSW/WLW ratio of 3.3–4 with sampling of 16–32 points is a good general guideline for WASSR imaging if the requirement for the offset error is to be less than about 0.3 Hz.

Glycogen Phantom Studies

Fig. 4A shows a typical series of CEST maps obtained for a glycogen phantom after optimized shimming, placing the transmitter offset on the average water resonance frequency, and voxel-by-voxel determination of MTRasym. The CEST images are shown for the offset frequency where glycogen CEST (glycoCEST) signals are expected to be optimal (19), namely between 0.75 and 1.5 ppm from water. A clear spatial intensity fluctuation is visible even for this homogeneous solution. When looking at the corresponding CEST Z-spectrum acquired in a single voxel placed as indicated at the accompanying image, it is not clear whether the data are centered correctly, because only a single, rather broad saturation curve is visible representing the combined effects of CEST and direct saturation. However, the corresponding narrow WASSR Z-spectrum, which has negligible MTC and CEST effects due to the use of brief low power saturation, clearly shows that the spectrum was not centered properly, resulting in a non-zero WASSR MTRasym. The frequency shift of the center of the direct saturation curve in this particular voxel was only 0.08 ppm (or 10.6 Hz) off center, illustrating the exquisite sensitivity of the WASSR approach due to the steep slope of the saturation curve. Even though this shift difference is very small, the CEST MTRasym spectrum (middle of Fig. 4A, solid line) shows a clear deviation from zero. The CEST maps after correction (left of Fig. 4B) show remarkable improvement in image homogeneity, compared to maps without (left of Fig. 4A). The spatial distribution of signal intensity is more uniform over the phantom and CEST maps show maximum signal in the range of about 0.75–1.25 ppm frequency offset, known to correspond to the glycogen OH proton frequency range (19). Note that the erroneous MTRasym peak (right of Fig. 4A, solid line) disappears in the corrected WASSR spectrum (right of Fig. 4B, solid line), confirming that the symmetric direct saturation spectrum was centered properly.

Figure 4
CEST maps, single voxel-based Z-spectra and MTRasym plots from CEST and WASSR obtained for the 200 mM glycogen phantom before (A) and after (B) WASSR correction

Fig. 5 shows WASSR shift maps and CEST MTRasym signal processed as a function of WN. While visual inspection of the shift maps in Fig. 5A does not show any obvious difference, the effect of WN choice on quantification is revealed clearly when looking at the actual MTRasym numbers (Fig. 5B and Table 2). When using WN = 16, a glycoCEST effect that is consistent around 11.2–12.2% can be found for ROIs reflecting regions with different local fields as indicated in the accompanying reference image. The average glycoCEST effect for this phantom was 11.5%.

Figure 5
WASSR shift maps and quantification of glycoCEST MTRasym as a function of number of WASSR points for the glycogen phantom
Table 2
Dependence of WASSR-corrected glycoCEST MTRasym in the 200 mM glycogen phantom on the number of WASSR points.

Human Studies

To evaluate whether the WASSR method can be used to elucidate the hidden glycoCEST effects in tissue at 3T, we performed WASSR experiment on human calf muscle. Fig. 6A shows the initial CEST maps (left) for subject S1 obtained without correction. Similar to the phantom, large signal fluctuations are found, which may be erroneously interpreted as variations in glycogen content in different muscle regions. When looking at the corresponding CEST Z-spectrum in a single voxel placed in a gastrocnemius calf muscle, it appears to be centered properly. This is due to the convolution of CEST, MTC and direct saturation effects, which, at the power level needed to see a CEST effect, merge together into a single broad saturation curve at this low field strength. However, the water offset can be easily found in the WASSR Z-spectrum, similar to the phantom results. The frequency shift of this particular voxel is about 0.08 ppm (or 10.29 Hz), while shift variation over the entire cross-section of calf muscle was found to be around 0.47 ppm (or 60 Hz) in WASSR shift maps for the four subjects (Fig. 6D). When applying the WASSR-based shift correction in each voxel, the glycoCEST maps show improved spatial homogeneity over the calf muscle (right of Fig. 6A). The corresponding CEST MTRasym plot clearly shows positive signal intensity of about 6% within the expected frequency range for glycogen and close-to-zero intensity outside that range. In Fig. 6C, the integrated glycoCEST signal (range 0–2 ppm) before (left) and after (right) WASSR correction is compared for four subjects. While uncorrected glycoCEST numbers show large variations, the WASSR corrected values are more consistent and of lower standard error. The results clearly represent that glycoCEST values in different muscle regions can be quantified if proper absolute water frequency correction is applied, while uncorrected CEST values show random variation, including erroneous negative numbers.

Figure 6
glycoCEST maps in human calf muscle before and after WASSR correction


The results show that it is possible to use direct water saturation to precisely (sub-Hz accuracy) map the absolute water frequency on a voxel-by-voxel basis. This so-called WASSR acquisition employs a brief low-power RF saturation pulse, for which magnetization transfer effects (both CEST and MTC) are expected to be negligible and for which the direct water saturation spectrum is much narrower than during the actual CEST acquisition, where higher power broadens the curve. For this approach to be successful as reference for CEST studies, it is required to have an unchanged magnetic field between the WASSR and CEST acquisitions, which is easily feasible in practice by not performing a separate prescan (shimming and water offset determination) before the CEST acquisition that follows the WASSR mapping. As an ultimate demonstration of this approach, it was shown that glycoCEST effects could be detected at 3T in vivo in muscle even though, at this field strength, the glycoCEST resonance overlaps so much with the direct saturation curve that a separate resonance cannot be distinguished.

Technical considerations

As any other MRI approach, the accuracy of the WASSR method depends highly on the MRI parameters used to acquire the data and the experimental conditions. The parameters of interest are the line width of the direct saturation curve and the spread of the magnetic field over the region of interest, which together determine the choice of sweep width and the number of points needed to sample the frequency spread. Based on our MC simulations, the optimum range of WSW was found to be 3.3–4 times the WLW with a minimum sampling requirement of 16 points in order to achieve a sub-Hertz (< 0.3Hz) offset error for SNR(S0) values of 80 or more. The use of only a small WN is possible due to the use of a simple maximum-symmetry algorithm that takes into account the complete WASSR Z-spectrum to estimate the center frequency (Fig. 2B). The above sampling guideline was found to be well applicable in the glycogen phantom and in vivo in muscle, but it is possible that the WSW/WLW may have to be increased for organs with large field variations and motion effects, such as the liver. When respiratory motion is an issue, data acquisition probably will need to be triggered to the respiration to secure similar field patterns between WASSR and the CEST acquisitions. In addition, the time between the two acquisitions should be kept as short as possible. For the current study, we used comparable resolution for CEST and WASSR, but a reduction of the number of phase-encodes in WASSR by at least a factor of two should be well possible, reducing the reference scan time to only 1–2 min.

The observant reader will notice that the WASSR shift map clearly resembles a typical field map (32,33) as acquired routinely on MR scanners and may wonder whether it is possible to perform the shift centering using standard gradient-echo based field mapping. The answer is that it is not, because there is no reference voxel where the exact center frequency is known. Even though the frequency differences between all voxels are correct, the center frequency in the voxels may still be off. With respect to this, it is important to realize that, when optimizing the scanner frequency offset, the RF is placed on the water signal maximum, which is summed over the complete field of view. So there is no indication as to where the exact center is for a particular voxel within a few Hz, the precision needed for these CEST experiments. Recently Sun et al. (34) reported that B0 field map can be used in combination with other information, such as T1, T2 and B1 field maps. Although this algorithm should work for studies with small ranges of B0 shifts, it may not be applicable for large frequency shifts that are commonplace for in vivo human data. In comparison, the algorithm described in this study requires only a WASSR map (additional scan with length of around 1–2 min), while providing reasonable correction for any shift range in B0. An additional advantage of using a WASSR shift map instead of a field map is that WASSR and CEST experiments are obtained using identical acquisition schemes, avoiding a mismatch in field-based image distortions between the two data sets. Another suitable approach to obtain an image of absolute water frequency would be water-based spectroscopic imaging (35). This would require 2D phase encoding followed by spectral acquisition. A typical matrix of 64 × 64 at a TR of 1 s with a SENSE factor of 2 in each direction would take over 15min. As there is sufficient signal, TR can probably be reduced and a combination of such fast acquisition and time-domain analyses may be possible in the future.

Previous frequency correction approaches used cubic-spline and polynomial fitting applied to a z-spectrum acquired using the CEST acquisition parameters. WASSR, on the other hand, is an MRI acquisition method employing acquisition parameters that differ from CEST and provides an absolute water frequency map. This new principle that can be used in combination with many fitting and interpolation approaches and, as such, completely different from these previous z-spectral fitting approaches. It can be used independent of any interfering effects, such as MT asymmetry and destructive interference of CEST and direct saturation effects. One question is whether it will always be necessary to use WASSR and the answer is no. When the CEST and direct saturation effects are sufficiently separated, it may be possible to fit the exact minimum without the WASSR scans. This will be the case under conditions of sufficient SNR, spectral resolution, and sufficiently large CEST effects. However, depending on the situation studied, this may require a lengthy acquisition of the whole z-spectrum to achieve sufficient SNR or spectral resolution, which may not be feasible in vivo, especially in humans. The WASSR approach on the contrary, always provides a correct map and can be acquired rapidly (~1min), if necessary even before and after the high-SNR CEST acquisition. This will allow reduction of the number of points in the z-spectra to a few in a small frequency range around the CEST offset of interest and the mirrored reference range, so that multiple acquisitions can be done to increase SNR. Thus, WASSR is definitely necessary for glycoCEST at lower field and may be the method of choice for APT and PARACEST in cases of limited SNR and available scan time or with complex MT asymmetry patterns. Also noteworthy is that WASSR allows accuracy of a few Hz or less, i.e. a few hundredths of a ppm while rapidly acquiring only 16–32 spectral points with single acquisition.

Application to glycoCEST MRI

As a first demonstration of the WASSR method, we studied a phantom with bovine liver glycogen (200 mM) and tried to detect glycogen signals in the human calf muscle in vivo. The glycogen hydroxyl groups resonate around 0.75 and 1.25 ppm from water (19) and, at the RF power levels needed for proper CEST magnetization transfer, cannot be distinguished as a separate resonance in the Z-spectrum at the field of 3T. In the phantom, the glycoCEST effect thus appears as an asymmetric broadening of the direct water saturation curve, thereby complicating its detection (Fig. 4) due to uncertainty about the exact water frequency offset. Without the WASSR referencing the CEST maps showed positive and negative MTR asymmetries indicating that it is not possible to reproducibly discern glycogen. When correcting for the absolute water frequency using WASSR (Fig. 5), a homogeneous glycogen distribution could be detected when using 16 or more points for the frequency mapping, in line with expectations based on our simulations. The averaged MTRasym for 200 mM glycogen was about 11.5 ± 0.6%.

As a first in vivo application, we studied the potential of WASSR to assist in detection glycogen in the calf muscle. The results for four subjects showed that WASSR correction was essential to find reproducibly positive glycoCEST effects (Fig. 6C). This should open up the way for in vivo glycogen MRI studies, which may add diagnostic capability for many common diseases, including obesity, insulin resistance and Type-2 diabetes. The resulting quantified CEST effects in Table 3 give a range of MTRasym values of 4.8–5.5% for four different muscles, the soleus, the medial and lateral heads of the gastrocnemius, and the tibialis anterior. Based on the phantom result, this would correspond to a glycogen concentration of around 100 mM or less in the leg, which is not unreasonable in view of the literature value of around 50–100 mM (36). As the CEST effect is complex and its quantification depends on many parameters, including T1, T2, pH, scan repetition time, and even B1, a more detailed study comparing glycoCEST with the gold-standard 13C studies (37,38) is needed to better assess absolute quantification. Such a study is currently being designed.

Table 3
WASSR-corrected glycoCEST MTRasym for in vivo calf muscle


A new absolute water frequency mapping approach was demonstrated that employs the direct water saturation obtained with a brief low power RF pulse. Guidelines for choice of sweep-width and sampling rate were provided. As a first demonstration, this so-called WASSR approach was applied to center the Z-spectrum of glycogen in phantoms, where the OH groups resonate so close to the water protons that the CEST and direct water saturation resonances are convolved even at 3T. The resulting data show that WASSR allowed identification of glycoCEST effects in both phantoms and the human calf muscle. The WASSR correction method is expected to be beneficial for accurate quantification of CEST effects with wide range of application both in vivo tissue and in vitro with MRI and MRS.


The authors thank Ms. T. Brawner and Dr. C. Jones for technical assistance and Drs. P. Barker, M. McMahon and Q. Qin for helpful discussions. Equipment used in this study was manufactured by Philips. Dr. Kim is supported by a grant from Philips Medical Systems to Kennedy Krieger Research Institute. Dr. van Zijl is a paid lecturer for Philips Medical Systems and is the inventor of technology that is licensed to Philips. This arrangement has been approved by Johns Hopkins University in accordance with its conflict of interest policies. This publication was made possible in part by a grant from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH). Its contents are solely the responsibility of the authors and do not necessarily represent the official view of NCRR or NIH.

Grant support: NIH-NCRR P41-RR15241, NIH-NIBIB R01-EB02634, R21-EB02666, Philips Medical System


1. Ward KM, Aletras AH, Balaban RS. A new class of contrast agents for MRI based on proton chemical exchange dependent saturation transfer (CEST) J Magn Reson. 2000;143(1):79–87. [PubMed]
2. Ward KM, Balaban RS. Determination of pH using water protons and chemical exchange dependent saturation transfer (CEST) Magn Reson Med. 2000;44(5):799–802. [PubMed]
3. Goffeney N, Bulte JW, Duyn J, Bryant LH, Jr, van Zijl PC. Sensitive NMR detection of cationic-polymer-based gene delivery systems using saturation transfer via proton exchange. J Am Chem Soc. 2001;123(35):8628–8629. [PubMed]
4. Zhou J, van Zijl PC. Chemical exchange saturation transfer imaging and spectroscopy. Progr NMR Spectr. 2006;48(2–3):109–136.
5. Zhang S, Merritt M, Woessner DE, Lenkinski RE, Sherry AD. PARACEST agents: modulating MRI contrast via water proton exchange. Acc Chem Res. 2003;36(10):783–790. [PubMed]
6. Aime S, Barge A, Delli Castelli D, Fedeli F, Mortillaro A, Nielsen FU, Terreno E. Paramagnetic lanthanide(III) complexes as pH-sensitive chemical exchange saturation transfer (CEST) contrast agents for MRI applications. Magn Reson Med. 2002;47(4):639–648. [PubMed]
7. Aime S, Delli Castelli D, Fedeli F, Terreno E. A paramagnetic MRI-CEST agent responsive to lactate concentration. J Am Chem Soc. 2002;124(32):9364–9365. [PubMed]
8. Yoo B, Pagel MD. A PARACEST MRI contrast agent to detect enzyme activity. J Am Chem Soc. 2006;128(43):14032–14033. [PubMed]
9. Liu G, Li Y, Pagel MD. Design and characterization of a new irreversible responsive PARACEST MRI contrast agent that detects nitric oxide. Magn Reson Med. 2007;58(6):1249–1256. [PubMed]
10. Zhou J, Payen JF, Wilson DA, Traystman RJ, van Zijl PC. Using the amide proton signals of intracellular proteins and peptides to detect pH effects in MRI. Nat Med. 2003;9(8):1085–1090. [PubMed]
11. van Zijl PC, Zhou J, Mori N, Payen JF, Wilson D, Mori S. Mechanism of magnetization transfer during on-resonance water saturation. A new approach to detect mobile proteins, peptides, and lipids. Magn Reson Med. 2003;49(3):440–449. [PubMed]
12. Zhou J, Lal B, Wilson DA, Laterra J, van Zijl PC. Amide proton transfer (APT) contrast for imaging of brain tumors. Magn Reson Med. 2003;50(6):1120–1126. [PubMed]
13. Ling W, Regatte RR, Navon G, Jerschow A. Assessment of glycosaminoglycan concentration in vivo by chemical exchange-dependent saturation transfer (gagCEST) Proc Natl Acad Sci U S A. 2008;105(7):2266–2270. [PubMed]
14. Snoussi K, Bulte JW, Gueron M, van Zijl PC. Sensitive CEST agents based on nucleic acid imino proton exchange: detection of poly(rU) and of a dendrimer-poly(rU) model for nucleic acid delivery and pharmacology. Magn Reson Med. 2003;49(6):998–1005. [PubMed]
15. Aime S, Delli Castelli D, Lawson D, Terreno E. Gd-loaded liposomes as T1, susceptibility, and CEST agents, all in one. J Am Chem Soc. 2007;129(9):2430–2431. [PubMed]
16. Zhang S, Winter P, Wu K, Sherry AD. A novel europium(III)-based MRI contrast agent. J Am Chem Soc. 2001;123(7):1517–1518. [PubMed]
17. Woods M, Woessner DE, Zhao P, Pasha A, Yang MY, Huang CH, Vasalitiy O, Morrow JR, Sherry AD. Europium(III) macrocyclic complexes with alcohol pendant groups as chemical exchange saturation transfer agents. J Am Chem Soc. 2006;128(31):10155–10162. [PMC free article] [PubMed]
18. Gilad AA, McMahon MT, Walczak P, Winnard PT, Jr, Raman V, van Laarhoven HW, Skoglund CM, Bulte JW, van Zijl PC. Artificial reporter gene providing MRI contrast based on proton exchange. Nat Biotechnol. 2007;25(2):217–219. [PubMed]
19. van Zijl PC, Jones CK, Ren J, Malloy CR, Sherry AD. MRI detection of glycogen in vivo by using chemical exchange saturation transfer imaging (glycoCEST) Proc Natl Acad Sci U S A. 2007;104(11):4359–4364. [PubMed]
20. Bryant RG. The dynamics of water-protein interactions. Annu Rev Biophys Biomol Struct. 1996;25:29–53. [PubMed]
21. Jones CK, Schlosser MJ, van Zijl PC, Pomper MG, Golay X, Zhou J. Amide proton transfer imaging of human brain tumors at 3T. Magn Reson Med. 2006;56(3):585–592. [PubMed]
22. Pekar J, Jezzard P, Roberts DA, Leigh JS, Jr, Frank JA, McLaughlin AC. Perfusion imaging with compensation for asymmetric magnetization transfer effects. Magn Reson Med. 1996;35(1):70–79. [PubMed]
23. Swanson SD. Protein mediated magnetic coupling between lactate and water protons. J Magn Reson. 1998;135(1):248–255. [PubMed]
24. Prammer MG, Haselgrove JC, Shinnar M, Leigh JS., Jr A new approach to automatic shimming. J Magn Reson. 1988;77(11):40–52.
25. Schneider E, Glover G. Rapid in vivo proton shimming. Magn Reson Med. 1991;18(2):335–347. [PubMed]
26. Webb P, Macovski A. Rapid, fully automatic, arbitrary-volume in vivo shimming. Magn Reson Med. 1991;20(1):113–122. [PubMed]
27. Maudsley AA, Simon HE, Hilal SK. Magnetic field measurement by NMR imaging. J Phys E: Sci Instrum. 1984;17(3):216–220.
28. Tropp JS, Derby KA, Hawryszko C, Sugiura S, Yamagata H. Automated shimming of B0 for spectroscopic imaging. J Magn Reson. 1989;85:244–254.
29. Lagarias JC, Reeds JA, Wright MH, Wright PE. Convergence Properties of the Nelder--Mead Simplex Method in Low Dimensions. SIAM J on Optimization. 1998;9(1):112–147.
30. Mulkern RV, Williams ML. The general solution to the Bloch equation with constant rf and relaxation terms: application to saturation and slice selection. Med Phys. 1993;20(1):5–13. [PubMed]
31. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: sensitivity encoding for fast MRI. Magn Reson Med. 1999;42(5):952–962. [PubMed]
32. Webb P, Spielman D, Macovski A. Inhomogeneity correction for in vivo spectroscopy by high-resolution water referencing. Magn Reson Med. 1992;23(1):1–11. [PubMed]
33. Dagher AP, Balaban RS, Aletras A. BEST-CEST: B0-Error Shift Tracking Chemical Exchange Saturation Transfer Imaging. Proceedings of the 10th Annual Meeting of ISMRM; Honolulu, Hawaii, USA. 2002. p. 611.
34. Sun PZ, Farrar CT, Sorensen AG. Correction for artifacts induced by B(0) and B(1) field inhomogeneities in pH-sensitive chemical exchange saturation transfer (CEST) imaging. Magn Reson Med. 2007;58(6):1207–1215. [PubMed]
35. Spielman D, Webb P, Macovski A. Water referencing for spectroscopic imaging. Magn Reson Med. 1989;12(1):38–49. [PubMed]
36. Price TB, Rothman DL, Avison MJ, Buonamico P, Shulman RG. 13C-NMR measurements of muscle glycogen during low-intensity exercise. J Appl Physiol. 1991;70(4):1836–1844. [PubMed]
37. Cline GW, Petersen KF, Krssak M, Shen J, Hundal RS, Trajanoski Z, Inzucchi S, Dresner A, Rothman DL, Shulman GI. Impaired glucose transport as a cause of decreased insulin-stimulated muscle glycogen synthesis in type 2 diabetes. N Engl J Med. 1999;341(4):240–246. [PubMed]
38. Magnusson I, Rothman DL, Katz LD, Shulman RG, Shulman GI. Increased rate of gluconeogenesis in type II diabetes mellitus. A 13C nuclear magnetic resonance study. J Clin Invest. 1992;90(4):1323–1327. [PMC free article] [PubMed]