Spin Echo Dynamics as Codes for Quantitative Imaging
Gigi Galiana1

1Diagnostic Radiology, Yale University, New Haven, CT, United States

Synopsis

Recent work has highlighted the complicated and distinctive dynamics that shape the signal observed during an irregularly-spaced train of spin echoes. Here, we use those signals as codes that allows us to identify mixtures of molecules that are too similar to distinguish by standard spectroscopy. Extensions to water based systems, where relaxation mechanisms are the primary driver of signal dynamics, will also be presented.

Purpose

Small amounts of diglycerides (DAGs) in lipid depots are highly significant to metabolic health.1 But diglycerides are difficult to distinguish from triglycerides (TAGs) with MR, since the small differences in chemical shift and J-coupling lead to nearly identical spectra. However, other studies of lipids have shown that even small differences in chemical shift, J-coupling and relaxation dynamics lead to large differences in the signal observed under spin echo trains with varied spacings.2,3

This work shows that pseudo-random pulse sequences, like irregularly spaced spin echo trains, can rapidly generate signal codes that better distinguish spectrally similar molecules. The resulting codes (determined either experimentally or with product operator simulations) can be used to analyze data from a mixture of compounds with either dictionary based methods or linear decomposition approaches.4,5 Though initial work is focused on demonstrating the feasibility of distinguishing DAG and TAG, the general approach may be useful for other metabolites or compounds that are spectrally too similar to distinguish through free evolution alone. It may also provide a new source of contrast for waters in structured media like tissue.

Methods

Spinach, a product operator based spin simulator, was used to model spin networks based on the glycerol backbones of DAG and TAG.6 These spin systems were used to generate free induction decays as well as signals resulting from trains of 16 spin echoes generated with random spacings of the 180 degree refocusing pulses. A number of pulse sequences were studied with random echo spacings chosen from different distributions. Each sequence was simulated on both compounds and one sequence was then used to analyze mixtures of various concentrations where the T2 of the measured signal differed from the modeled T2.

As preliminary experimental verification, signals were collected from diacetin, triacetin, and mixtures of these two. Though commercial technical grade diacetin is a mixture of several molecules, we regarded it as a pure compound, and bottles were prepared with diacetin/triacetin volume % of 100/0, 0/100, 50/50, 10/90 and 5/95. These were studied in a 3T Siemens Trio using both single voxel spectroscopy and multiecho acquisition with 1D spatial phase encoding.

Results

Simulations suggest that families of pulse sequences with echo spacings chosen randomly from the same distribution can have very different ability to distinguish the glyceride backbones. A histogram of the correlation between diglyceride and triglyceride signals across all experiments is shown in Figure 1 and compares favorably against the correlation between simulated spectra (i.e. free evolution), which was found to be 0.79. A sequence with median performance (correlation = .29) was used to study the effect of erroneous T2 for a range of mixtures. As shown in Figure 2, errors of up to 30% in the presumed T2 have little effect on the quantification of diglycerides, though some bias can be detected.

Experimentally measured spectra and multiecho signals for “pure” triacetin and diacetin are shown in Figure 3. These codes for the “pure” liquids were then used to analyze mixtures of diacetin/triacetin using standard linear decomposition. The codes from this first experiment were then used to analyze the data acquired in a second study. Those results are shown in Figure 4, using volume % as the x-axis. (Molar fractions are difficult to calculate without a more thorough chemical analysis of the diacetin.) While neither method gives quantitative agreement with the ground truth, the multiecho approach shows considerably greater accuracy.

Discussion and Conclusion

Modeling spin behavior under arbitrary pulse sequences may better distinguish species than free evolution, and it may be a promising new approach to MRS in vivo. Furthermore, extensions to water spins in structured samples, where relaxation dynamics are the primary driver of echo dynamics, will also be explored. Together these methods could provide a powerful new source of quantitative contrast.

Acknowledgements

We gratefully acknowledge Hemant Tagare for providing expertise on the decomposition methods studied in this work and Todd Constable for useful commentary on the experiments.

References

[1] Shulman G, Ectopic Fat in Insulin Resistance, Dyslipidemia, and Cardiometabolic Disease, N Engl J Med 2014;371:1131-41;

[2] Jenista ER, Stokes AM, Branca RT, Warren WS, Optimized, unequal pulse spacing in multiple echo sequences improves refocusing in magnetic resonance, J Chemical Physics 2009;131:2045101-7.

[3] Stokes, AM, Feng Y, Mitropoulos T, Warren WS, Enhanced Refocusing of Fat Signals Using Optimized Multipulse Echo Sequences, MRM 2013;69:1044-4055.

[4] Ma D, Gulani V, Seiberlich N, Liu K, Sunshine JL, Duerk JL, Griswold MA, Nature 2013;495:187-192.

[5] Galiana G, Constable RT, PLOS One 2014; DOI: 10.1371/journal.pone.0086008.

[6] Hatzakis E, Agiomyrgianaki A, Kostidis S, Dais P, High Resolution NMR Spectroscopy: An Alternative Fast Tool for Qualitative and Quantitative Analysis of Diacylglycerol (DAG) Oil, J Am Oil Chem Soc 2011;88:1695-1708.

Figures

Figure 1: Histogram of correlation between simulated DAG/TAG signals for various 8 echo sequences (multiple random spacings chosen from distributions of various means and variances) shows a wide range.

Figure 2: T2 errors of up to 10% had negligible effect on the estimated concentration, while errors up to 30% only slightly bias the quantification of DAG.

Figure 3: Experimental spectra from triacetin (red) and diacetin (blue) are less distinguishable than echo trains.

Figure 4: Multiecho signals (blue) yield more accurate experimental quantification of triacetin/diacetin than single voxel spectroscopy (red). Here, the data acquired from pure triacetin and diacetin in experiment 1 (either multiecho codes or spectra) were used to analyze data from mixtures of triacetin and diacetin. That same data was used to analyze data from all five vials in a second experiment. The plot above summarizes results from both methods acquired in the two separate experiments.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3000