Hyperpolarized 13C spectroscopic imaging using single‐shot 3D sequences with unpaired adiabatic refocusing pulses

Hyperpolarized MRI with 13C‐labeled metabolites has enabled metabolic imaging of tumors in vivo. The heterogeneous nature of tumors and the limited lifetime of the hyperpolarization require high resolution, both temporally and spatially. We describe two sequences that make more efficient use of the 13C polarization than previously described single‐shot 3D sequences. With these sequences, the target metabolite resonances were excited using spectral‐spatial pulses and the data acquired using spiral readouts from a series of echoes created using a fast‐spin‐echo sequence employing adiabatic 180° pulses. The third dimension was encoded with blipped gradients applied in an interleaved order to the echo train. Adiabatic inversion pulses applied in the absence of slice selection gradients allowed acquisition of signal from odd echoes, formed by unpaired adiabatic pulses, as well as from even echoes. The sequences were tested on tumor‐bearing mice following intravenous injection of hyperpolarized [1‐13C]pyruvate. [1‐13C] pyruvate and [1‐13C] lactate images were acquired in vivo with a 4 × 4 × 2 cm3 field of view and a 32 × 32 × 16 matrix, leading to a nominal resolution of 1.25 × 1.25 × 1.25 mm3 and an effective resolution of 1.25 × 1.25 × 4.5 mm3 when the z‐direction point spread function was taken into account. The acquisition of signal from more echoes also allowed for an improvement in the signal‐to‐noise ratio for resonances with longer T 2 relaxation times. The pulse sequences described here produced hyperpolarized 13C images with improved resolution and signal‐to‐noise ratio when compared with similar sequences described previously.

lactate can be used to assess treatment response, 7 detect occult disease, 4 and monitor disease progression. 8,9 There are four main challenges in pulse sequence design for hyperpolarized 13 C MRI 3 : (i) the number of pulses should be kept to a minimum, as each excitation depletes polarization; (ii) spatial resolution should be as high as possible; (iii) image acquisition should be fast to cope with the rapidly decaying polarization; and (iv) spectral resolution should be sufficient to resolve signals from different metabolites in the resulting images. Many imaging pulse sequences have been developed to meet these demands. Cunningham et al. proposed an echo-planar-based spectroscopic imaging (EPSI) technique for mapping multiple metabolites 10

and
Larson and colleagues combined this with multiband spectral-spatial excitations in order to optimize use of the polarization in individual metabolites. 11 Compressed sensing has been used to enhance imaging efficiency, 12 exploiting the sparse nature of the hyperpolarized signal, both spatially and spectrally. This led subsequently to the development of a highly accelerated EPSI technique that exploited both multiband spectral-spatial excitation (SpSp) pulses and compressed sensing. 13 Others have used single-band SpSp pulses to image a single metabolite at a time, for example with 2D and 3D echo planar imaging (EPI), 14,15 which can be accelerated with compressed sensing. 16 The EPI trajectory can also be designed to enable separation of metabolite resonances in the resulting images. 17 Spiral chemical shift imaging (CSI) has been studied extensively, [18][19][20] while traditional CSI with compressed sensing has also shown promise for dynamic measurements. 21 2D spatial and 1D spectral information can also be acquired after a single excitation by exploiting spatiotemporal encoding. 22 We recently proposed a single-shot 3D sequence that makes efficient use of the transverse magnetization after each excitation. 23 The sequence exploits the B 1 insensitivity of adiabatic pulses, following a dual-spin-echo (DSE) scheme that is widely used in other hyperpolarized MRI applications. [10][11][12]17,24,25 The adiabatic pulses in the DSE scheme result in better refocusing and hence less signal loss due to dephasing, but the pulse sequence is inefficient in its use of spin hyperpolarization, since signal is acquired from only half of the spin echoes. Here we propose an enhanced version of this single-shot 3D sequence that uses unpaired adiabatic inversion pulses for refocusing and doubles the number of encoding steps in the z direction with the same number of pulses. We also describe an alternative design that can significantly enhance the signal-tonoise ratio (SNR), at the expense of a shorter imaging time window.

| Tumor imaging
Tumors were obtained by subcutaneous implantation of EL4 murine lymphoma cells, as described in Reference 23. Mice were fasted for 6-8 h before imaging. 26 Images were acquired at 7 T and the [1-13 C] pyruvate was hyperpolarized as described previously. 23

| MRI scanner
Experiments were performed at 7 T (Agilent, Palo Alto, CA, USA), with a 42 mm diameter bird-cage volume coil for 1 H transmission and reception, and a similar volume coil for 13 C transmission and a 20 mm diameter surface coil for 13 C detection (Rapid Biomedical, Rimpar, Germany).

| Pulse sequence
Target metabolite resonances were excited using a 10.016 ms SpSp pulse 23 with a bandwidth of 350 Hz and 1645 Hz (center to center) between replicate excitation bands. A 90°pulse required a B 1 field of 131.11 μT and used a fly-back trajectory to reduce vulnerability to eddy currents, with maximum gradient strengths of 0.19 and 0.2993 T/m at positive and negative polarities, giving a 1.2 cm minimum slab thickness. A fast-spin-echo (FSE) train, with four refocusing pulses, followed the excitation pulse, with a stack of four spiral acquisitions during each spin echo ( Figure 1A).
Removal of slice selection gradients from the refocusing pulses used previously in a similar sequence 23 allowed signal acquisition from both even and odd echoes. Previously spirals were not acquired from odd echoes because a fully refocused spin echo, in the presence of slice selection gradients, requires a pair of hyperbolic secant (HS) adiabatic pulses, 27 otherwise a quadratic phase is left across the swept frequency range, as illustrated in Figure 2C. This non-linear phase generated by an unpaired HS pulse results in signal loss as each k-space point is the summation of signals from a range of frequencies. When this range of frequencies is large, as is the case when a slice selection gradient is used, then signal loss is considerable. In the absence of a slice selection gradient the frequency range over which there is phase variation depends only on local B 0 field variation. In a reasonably well-shimmed B 0 field this frequency range is much smaller (of the order of ±50 Hz as compared with ±5 kHz when a slice selection gradient is present) and moreover shows a smooth variation, leading to only mild signal loss ( Figure 2D). Each spiral is encoded in the xy plane and each stack of spirals is phase encoded in the z direction with blipped gradients such that a 3D k-space is acquired. The z-dimension information is phase encoded in an interleaved manner and the k z index of each spiral ( Figure 1B)  The phase-encoding gradients on the z axis are self-refocusing in each acquisition interval in order to comply with the Carr-Purcell-Meiboom-Gill (CPMG) condition. In each acquisition interval the spiral dimensions in the k xy plane are 4 × 4, 16 × 16, 32 × 32, 8 × 8, in the order displayed in Figure 1A, and the spiral dimensions for the 3D k-space are thus ( Figure 1B, 3.120, and 8.500 ms respectively. Each spiral encodes a field of view (FOV) of 4 × 4 cm 2 in the xy plane. This sequence is referred to as FSE-I hereafter. An alternative sequence (hereafter referred to as FSE-II) can be constructed such that the same z-dimension phase encoding as in Figure 1A is repeated in additional stacks 5 to 8 and the sampled points are averaged with the datasets acquired through stacks 1 to 4 (see Figure S1). A limitation of this trajectory design is that the less densely sampled regions in the corners of 3D k-space lead to a widened point spread function (PSF) in the z direction, which was investigated with the PSF simulations described below. A spherical or cylindrical 3D k-space would improve the z-direction PSF and enhance resolution in the z direction. The pulse was designed to have a bandwidth of 10 kHz and to achieve more than 99.9% inversion over a frequency band of ±4 kHz at a B 1 field strength between 100 and 6400 μT. C, The quadratic phase induced in transverse magnetization when the pulse is applied to magnetization lying on the y axis. The same quadratic phase could be observed across the whole slice when a slice selection gradient was applied. D, The quadratic phase variation is small within a frequency band of ±50 Hz of the frequency offset of the pulse The proposed pulse sequence and its k-space trajectory in the k z direction. A, FSE-I pulse sequence, which starts with an SpSp pulse and a train of four adiabatic refocusing pulses. A stack of four spiral acquisitions is acquired after each refocusing pulse, and each spiral is phase encoded in the z direction by blipped gradients on the z axis. For the FSE-II sequence four extra stacks of spiral acquisitions are acquired with the same refocusing pulses and phase encodings in the z direction as in the first four stacks. The two signal blocks are then averaged in order to improve SNR. The FSE-II sequence hence shares the same k-space trajectories as the FSE-I pulse sequence. B, k-space trajectory in the k z direction. Each spiral index in the sequence in A corresponds to a k z position in B Imperfect refocusing pulses can dramatically deplete the polarization, especially within their transition bands. Therefore adiabatic inversion pulses with hyperbolic-secant modulation (HSn, n = 8) 28 were used for refocusing, where each had a duration of 10 ms ( Figure 2A) and a bandwidth of 10 000 Hz ( Figure 2B shows the effect of the pulse on the longitudinal magnetization at a B 1 field strength of 100 μT). The pulses were designed to achieve perfect inversion (>99.9% of the longitudinal magnetization inverted) when the B 1 field was between 100 and 6400 μT. The large time-bandwidth product also ensured insensitivity to B 0 inhomogeneity and that there was only a small and smooth phase variation in the frequency response so that signal loss was minimized. The phase response of the pulse, when applied along the y axis, to the y magnetization is shown in Figure 2C. The central part of this phase response is expanded in Figure 2D.
To preserve the CPMG condition, the delay τ in Figure 1A was 15.284 ms, which was determined by the durations of the excitation and inversion pulses and the spirals. Since there are two spirals (eighth and ninth) acquired at the center of the k z direction, the resulting contrast is contributed by two different TE times (61.136 and 122.272 ms), which, in the FSE-II sequence, are averaged with another two echoes at 183.408 and 244.544 ms.

| Simulations
The PSF of the acquisition trajectory was simulated for the proposed sequences, where a constant k-space was sampled by the proposed k-space trajectory and 3D Fourier transformed into the image domain. 29 This was then compared with the PSF obtained using the forerunner of these sequences, the DSE sequence. 23 Despite additional signal obtained from the later echoes, the FSE-II sequence has the same k-space trajectory and hence the same PSF as the FSE-I sequence when relaxation is ignored.

| Pulse calibration
The adiabatic pulse was calibrated on a cylindrical phantom (7 mm inner diameter) filled with 5 M [1-13 C]lactate at thermal equilibrium. A pulseacquire sequence without any gradients was used, where 2048 data points were acquired at 4 μs intervals immediately after the pulse. The sequence swept a range of frequencies (from 0 Hz to 4500 Hz offset, with 300 Hz step size) and B 1 field strengths (from 40 to 470 μT). The inversion pulses were also tested on a cylindrical phantom of the same size injected with hyperpolarized [1-13 C] pyruvate and imaged using both the FSE-I and FSE-II sequences. All the gradients were turned off except for the slice selection gradient for the SpSp excitation pulse. To mimic the RF exposure in experiments in vivo, the SpSp pulse was targeted alternately to the [1-13 C] pyruvate and [1-13 C] lactate resonance frequencies, with flip angles of 7°and 45°respectively, a slice thickness of 20 mm, and a TR of 1 s (hence 2 s TR at each frequency). The whole acquisition lasted for 90 s.

| Phantom imaging
The original pulse sequence (DSE sequence 23 ) and the FSE-I sequence were tested on the cylindrical phantom containing thermally polarized [1-13 C] lactate. One end of the phantom was positioned at the iso-center (also the center of the FOV in the z direction) to examine the definition of the resulting images. The FOVs for both sequences were 4 × 4 × 2 cm 3 in x, y, and z dimensions, with eight phase-encoding steps in the z direction in the original DSE sequence and 16 steps in the FSE-I sequence. Identical spiral designs were used for the DSE sequence, which were the same as those described in Reference 23, except for the optimized refocusing lobe at the end of each spiral. The flip angle was 90°for both sequences, exciting a slab of 12 mm. The excited slab was thinner than the z-direction FOV to avoid wrap-around artifacts from both the transition band of the excitation pulse and the side lobes of the PSF in the z direction. The phantom was then repositioned with the center of the phantom at the iso-center in order to compare the SNR between the three sequences. To investigate signal loss caused by the quadratic phase imparted by unpaired adiabatic pulses, signal was acquired using the DSE and FSE-I sequences with the encoding gradients turned off on all three axes. The PSFs of the FSE-I and FSE-II sequences were simulated again by adding a weighting factor for signal decay to each k z plane using the measured signal loss due to unpaired adiabatic pulses, the [1-13 C] lactate T 2 value measured in the phantom (1 s), and a T 2 value (300 ms) estimated in vivo at 7 T. 25 The simulation was then repeated using a T 2 * decay of 12 ms measured at 7 T. 30 FSE proton images were acquired as 16 slices to cover the same volume as the 13  45°respectively. This helps to preserve the polarization for [1-13 C] pyruvate and enhances the image SNR for [1-13 C] lactate. The eighth pair of acquisitions was used as a reference scan, where the z-axis gradients were turned off except for the slice selection gradient and the spoiler.
The FOVs in all three dimensions were 4 × 4 × 2 cm 3 . The excitation bandwidth was 12 mm at half maximum to avoid wrap-around artifacts in the z direction. The total acquisition window was 90 s. The water proton frequency was used to estimate the [1-13 C] pyruvate and [1-13 C] lactate resonance frequencies. T 2 -weighted proton images (FSE, 128 × 128 matrix, 16 slices, 1.25 mm slice thickness) were acquired for positional reference. 13 C MRI images were reconstructed in MATLAB (MathWorks, Natick, MA, USA), as described in Reference 23.

| RESULTS
The adiabatic pulses achieved almost complete inversion at ±3500 Hz when the B 1 field was greater than 120 μT, matching their expected performance (Supplementary Figure S2A). Signals acquired from a phantom injected with hyperpolarized [1- 13   was almost no signal outside the excited slab in the FSE image, whereas signal could be observed in the DSE image. However, the full width at half maximum (FWHM) of the PSF, even for the FSE sequence, is about 4 times the nominal z-direction resolution (4.5 mm versus 1.25 mm). This is because the corners of 3D k-space were under-sampled with this spiral trajectory.
Axial images acquired using the DSE and FSE-II sequences, when the center of the phantom was at magnet center ( Figure 4A), showed that the SNRs of the two central z slices acquired using the FSE-I sequence were 0.86 and 0.83 times the SNRs acquired using the DSE sequence, reflecting the improved resolution, while the SNRs in slices acquired using the FSE-II sequence were 1.12 and 1.07 times the DSE SNRs ( Figure 4B). Expressing the SNR per unit volume, the FSE-I SNRs were 1.72 and 1.66 times the DSE SNRs and the FSE-II SNRs were 2.25 and 2.13 times greater. The SNR curves in Figure 4B show the slab profile for each sequence, which was in agreement with convolution of the slice response of the excitation pulse ( Figure 1C of Reference 23) with the PSFs in the z direction for each sequence.
Phantom data were also acquired with encoding gradients on all three axes turned off. The two acquired echoes in the DSE sequence and the four acquired echoes in the FSE-I sequence are shown in Figure 5A. The T 2 decay curve was obtained by fitting to the amplitude of the even echoes in the FSE-I sequence, as these were formed with paired adiabatic pulses and hence free of quadratic phase twist and thus loss of signal. Signals from the odd echoes in the FSE sequence were approximately 10% lower than expected from signal decay due to T 2 relaxation. This decrease was greater than that expected from the non-linear phase shift imparted by the adiabatic pulses ( Figure 5A) and imperfections in local B 0 and B 1 fields, and could be due to mis-setting of the frequency of the RF pulse and small shifts in the timing of the echoes in the presence of the adiabatic pulses. As shown in Figure 5B, signal loss due to unpaired adiabatic pulses increases slightly the width of the main lobe of the PSF by 2.3% while reducing the amplitude and bandwidth of the side lobes by 12.3% and 5.5%, respectively. However, these effects disappear when T 2 decay is taken into account for both FSE-I ( Figure 5B) and FSE-II sequences (not shown). When T 2 * decay was considered, the main lobe of the PSF was increased by 2.9% and the side lobe by 3.2%, compared with the T 2 decay only case, for both FSE-I and FSE-II sequences. The echoes from the central k-space in the FSE sequence were smaller in amplitude than their counterparts in the DSE sequence because the timing of the pulse train in the FSE sequence is different from that in the DSE sequence. In the DSE sequence, the minimum interval between the excitation pulse and the first refocusing pulse is determined solely by the pulse durations, while in the FSE sequence this interval also depends on the durations of the first and second spirals.

| DISCUSSION
While the previous DSE sequence provided high temporal and spatial resolution, 23 the FSE sequences described here doubled the spatial resolution in the z direction with no sacrifice of temporal resolution. Although the FSE-II sequence doubled the sequence length from 125 ms to 247 ms, in practice the frame rate is rarely set to below 250 ms/frame, as this would lead to fast depletion of the polarization.
Although the PSF of the FSE sequences is much sharper in the z direction, it has a pair of side lobes toward the ends of the z-direction FOV, which may result in image artifacts. Such artifacts can be observed in Figure 6. The slices at both ends of the z-direction FOV (20 mm) should show no signal due to the limited thickness of the excitation slab (12 mm). These artifacts are less obvious in images acquired with the DSE sequence, as they appear as image blur in the z direction. A window function could be applied to the k-space in the k z direction to reduce these artifacts, but this will render a PSF more similar to that of the DSE sequence and compromise the enhancement in image definition.
Alternatively, the end slices could be discarded from the final set of images. Simulations of the PSF in the phantom experiments showed that modulation, caused by both T 2 decay and signal loss in odd echoes, was negligible ( Figure 5B). Faster T 2 decay in vivo may result in greater modulation of the PSF, which is a common concern for FSE sequences, while the effects of signal acquisition from odd echoes should be relatively minor as the B 0 field is usually well shimmed on pre-clinical systems. In further simulations, the widths of the main lobes of the PSFs of the FSE sequences at a T 2 of 50 ms were half way between the widths in the FSE and DSE sequences when relaxation was ignored, demonstrating a resolution benefit for the FSE over the DSE sequence, even for metabolites with very short T 2 times. T 2 * decay during each acquisition interval can cause signal variations between the four spirals in each interval, which in general suppresses signals from the edge spirals in the 3D k-space and hence blurring in the z direction. However, the short durations of each spiral reduce this effect considerably.
Moreover, simulations with T 2 * = 12 ms, which was measured for [1-13 C] pyruvate at 7 T, 30 showed that the effect of a short T 2 * on the PSF was negligible.
While the nominal z-direction resolution of the FSE sequence (1.25 mm) was double that of the DSE sequence, the true resolution was approximately four times the nominal resolution due to the PSF of the spiral acquisition design. This is because the sampling density was much lower in the corners of the 3D k-space when compared with a fully sampled k-space. This could be improved, as discussed in Reference 23, by FIGURE 5 A, Signal loss in echoes formed with unpaired adiabatic refocusing pulses. DSE and FSE-I sequences were used to acquire signal from the lactate phantom shown in Figure 4, with all encoding gradients turned off. The second and fourth echoes acquired with FSE-I were fitted to give a curve that shows T 2 relaxation. The fitted curve matched well with the echoes acquired using the DSE pulse sequence. The first and third echoes in the FSE-I sequence, which were formed by unpaired adiabatic pulses, were 9.8% smaller than the echoes that would have been formed had there been no adiabatic phase twist (as indicated by the curve showing signal loss due to T 2 relaxation). B, Simulated PSFs, showing the theoretical case (solid line), the effects of signal loss at odd echoes (dashed line), and the effects of T 2 relaxation, as measured in the phantom (dotted line) and in vivo (dash-dotted line)  The key to implementation of these FSE sequences was the removal of the slice-selection gradients from the adiabatic refocusing pulses. The phantom experiments showed that signal loss caused by residual quadratic phase variation was less than 10% for both the first and third echoes ( Figure 5). However, this requires that the central frequencies of the adiabatic pulses closely match the resonance frequencies of the target metabolites in order to ensure that the majority of the resonance frequencies fall into the bottom of the phase variation well ( Figure 2D). Dynamic measurements of the metabolite frequencies would be desirable, and these could be embedded between the excitation pulse and the first refocusing pulse.
A concern for these FSE sequences is that off-resonance polarizations, for example, from untargeted metabolites or metabolites located in regions where their frequencies are off resonance, may be depleted by the transition bands of the inversion pulses. This could be avoided by using adiabatic pulses with larger bandwidths. Initial experiments using adiabatic pulses with a 1.5 kHz bandwidth were much less successful than those using pulses with a 10 kHz bandwidth. The lifetime of the observable signal with a 1.5 kHz bandwidth pulse was only half of that with 10 kHz pulses. Another potential problem is the fringe-field effect, 18,31 where the B 1 field drops rapidly at the ends of the transmit coil and results in imperfect adiabatic pulses that quickly deplete spin polarization in these areas. This phenomenon was observed in the experiments performed in vivo (Figure 8), where signals acquired using the FSE-I sequence persisted for much longer (60 s as compared with 30 s) than those acquired using the FSE-II sequence, which used double the number of inversion pulses. In contrast, in the phantom experiment, where the hyperpolarized [1-13 C] pyruvate remained at the iso-center, both sequences resulted in signals that decayed at the same rate (supplementary Figure S2B). The fringe-field effect could again be relieved by using pulses with a wider bandwidth, which would result in greater tolerance to B 1 field inhomogeneity. A representative B 1 map acquired using a Bloch-Siegert sequence, 32 with the same coil setup as used for the phantom and in vivo experiments, is shown in Supplementary Figure S4. A larger bandwidth requires increased RF power, and while this was not a problem in our preclinical system, where the B 1 field can go beyond 500 μT, on a clinical system pulse bandwidth may be limited by the available B 1 field. Spiral-based trajectories are relatively insensitive to motion artifacts because of rapid traversal of k-space and dense sampling of the center, but are prone to off-resonance effects. These can arise from B 0 field inhomogeneity and eddy currents and usually lead to blurred and rotated images as a result of phase error accumulation during a single spiral. 33 Moreover, with blipped z gradients, phase error accumulations within each stack of spirals can lead to gross image shifts and distortions. Reference scans, 34 where the z-encoding gradients were turned off, were used to correct phase errors between spirals in the acquisition train, while the phase error accumulation within each spiral was ignored as all the spirals were very short in duration. Calculations, similar to those used previously, 23 showed that the sequences described here are immune to B 0 field inhomogeneity within ± π=2 3:12 ms ¼ ±80:12 Hz (3.12 ms is the duration of Spirals 5 and 12 in Figure 1A, which contribute most to the phase error accumulated at the center of k-space). We can usually achieve a water proton linewidth in vivo of 40-70 Hz when a spectrum is acquired from the imaging volume; therefore, 13 C off-resonance frequencies should be much less than ±80 Hz. In addition, eddy current effects could be measured separately and removed retrospectively in the reconstruction process. 35 The lower SNR for the FSE-I pulse sequence, when compared with the DSE sequence, reflects its better resolution ( Figure 4B). Theoretically the DSE sequence should give an SNR that is ffiffiffi 2 p times that of the FSE-I SNR, if signal loss due toT 2and T 2 *-dependent relaxation is neglected and rectangular slice profiles are assumed. However, the DSE SNR for the central slices was only about 1.2 times that of FSE-I sequence. This is because the signals are also weighted by the PSFs of the two sequences. The DSE/FSE-I signal ratio for the central slice, calculated from the areas under the simulated PSFs, was about 1.7, which is less than the value of 2 expected from the lower nominal spatial resolution. Since twice the number of echoes were acquired in the FSE-I sequence and ffiffiffi 2 p times the noise, the DSE/FSE-I SNR ratio should be 1:7 2= ffiffiffi 2 p ≈1:2, in good agreement with what we observed experimentally. Other factors can also contribute to the SNR differences, including the non-linear phase in the FSE-I sequence and imperfect slice selection in the DSE sequence. The FSE-II sequence can be used to further improve the SNR of the FSE-I sequence, although there may be faster polarization decay as the FSE-II sequence uses more pulses. This need not be a problem for a kinetic analysis if the extra polarization loss due to imperfect refocusing pulses is incorporated into a general RF depletion term in the modified Bloch equations 36 (e.g., take all the RF depletion caused by an excitation pulse and the following train of inversion pulses as a single term). The FSE-II sequence improved the SNR substantially in phantom experiments and in two out of the three experiments performed in vivo with hyperpolarized [1-13 C] pyruvate. The reason that we failed to see an SNR improvement in the third mouse was not clear, but we infer that in this animal the T 2 of the [1-13 C] lactate resonance was shorter. Faster-decaying signals during the echo train may eliminate the benefit of acquiring signal from later echoes in the FSE-II sequence, where the first echo in the repeated echo train appears 152.84 ms after the excitation pulse. Assuming an echo train with perfectly refocused spin echoes, the T 2 will need to be at least 140 ms for there to be a gain in SNR by averaging signal from the later echoes. We have measured a T 2 for [1-13 C] lactate in this tumor model of 170 ms. 37 These pulse sequences have the potential for clinical translation, where they would benefit from the longer T 2 times at lower fields. Specific absorption rate (SAR) is a potential problem; however, the pulse bandwidth could be reduced and the adiabatic refocusing pulses could be replaced with non-selective composite pulses, 38 although obtaining phase uniformity and refocusing performance in inhomogeneous B 0 fields would be challenging. 39 The number of pulses in these sequences is already much smaller when compared with clinical FSE sequences, such as single-shot FSE and CUBE. 40,41 The repetition rate could also be reduced; for example, in a clinical study, where 2D imaging was used, the TR was 5 s. 4 Clinical scanners are usually limited to much lower gradient strengths and slew rates than those used here, which may lead to long spiral durations and hence unacceptable image distortion in the z direction. However, translation could still be possible if the FOV were kept small to cover a restricted region of interest so that short spirals could be used. Fringe-field effects, where the subject is not fully covered by the transmission coil, could be more of a problem. Spin polarization in a labeled molecule could be depleted rapidly if moving through a region with a rapidly varying B 1 field. Potential solutions to this problem include better-designed 13 C transmission coils that achieve a well-defined profile, such that the transition region, where the B 1 field drops off rapidly, is limited, and multichannel transmission to achieve more uniform in-plane polarization inversion.

| CONCLUSIONS
We have described a single-shot 3D FSE pulse sequence, which made more efficient use of the spin polarization and doubled the nominal zdirection resolution when compared with a similar sequence employing a DSE design. An alternative sequence that regained the SNR lost due to the improved spatial resolution has been described. Both sequences depend on the use of echoes formed by unpaired adiabatic refocusing pulses, where the quadratic phase twist is minimized by removing slice-selection gradients from the refocusing pulses. Although the primary applications of these sequences will be in pre-clinical studies, unpaired adiabatic pulses have the potential to be used for 13 C MRSI measurements in the clinic.