Monitoring tumor cell death in murine tumor models using deuterium magnetic resonance spectroscopy and spectroscopic imaging.

2H magnetic resonance spectroscopic imaging has been shown recently to be a viable technique for metabolic imaging in the clinic. We show here that 2H MR spectroscopy and spectroscopic imaging measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate can be used to detect tumor cell death in vivo via the production of labeled malate. Production of [2,3-2H2]malate, following injection of [2,3-2H2]fumarate (1 g/kg) into tumor-bearing mice, was measured in a murine lymphoma (EL4) treated with etoposide, and in human breast (MDA-MB-231) and colorectal (Colo205) xenografts treated with a TRAILR2 agonist, using surface-coil localized 2H MR spectroscopy at 7 T. Malate production was also imaged in EL4 tumors using a fast 2H chemical shift imaging sequence. The malate/fumarate ratio increased from 0.016 ± 0.02 to 0.16 ± 0.14 in EL4 tumors 48 h after drug treatment (P = 0.0024, n = 3), and from 0.019 ± 0.03 to 0.25 ± 0.23 in MDA-MB-231 tumors (P = 0.0001, n = 5) and from 0.016 ± 0.04 to 0.28 ± 0.26 in Colo205 tumors (P = 0.0002, n = 5) 24 h after drug treatment. These increases were correlated with increased levels of cell death measured in excised tumor sections obtained immediately after imaging. 2H MR measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate provide a potentially less expensive and more sensitive method for detecting cell death in vivo than 13C MR measurements of hyperpolarized [1,4-13C2]fumarate metabolism, which have been used previously for this purpose.

Deuterium MRSI of tumor cell death in vivo following oral delivery of 2 H-labeled fumarate.

PURPOSE: There is an unmet clinical need for direct and sensitive methods to detect cell death in vivo, especially with regard to monitoring tumor treatment response. We have shown previously that tumor cell death can be detected in vivo from 2 H MRS and MRSI measurements of increased [2,3-2 H2 ]malate production following intravenous injection of [2,3-2 H2 ]fumarate. We show here that cell death can be detected with similar sensitivity following oral administration of the 2 H-labeled fumarate. METHODS: Mice with subcutaneously implanted EL4 tumors were fasted for 1 h before administration (200 µl) of [2,3-2 H2 ]fumarate (2 g/kg bodyweight) via oral gavage without anesthesia. The animals were then anesthetized, and after 30 min, tumor conversion of [2,3-2 H2 ]fumarate to [2,3-2 H2 ]malate was assessed from a series of 13 2 H spectra acquired over a period of 65 min. The 2 H spectra and 2 H spectroscopic images were acquired using a surface coil before and at 48 h after treatment with a chemotherapeutic drug (etoposide, 67 mg/kg). RESULTS: The malate/fumarate signal ratio increased from 0.022 ± 0.03 before drug treatment to 0.12 ± 0.04 following treatment (p = 0.023, n = 4). Labeled malate was undetectable in spectroscopic images acquired before treatment and increased in the tumor area following treatment. The increase in the malate/fumarate signal ratio was similar to that observed previously following intravenous administration of labeled fumarate. CONCLUSION: Orally administered [2,3-2 H2 ]fumarate can be used to detect tumor cell death noninvasively following treatment with a sensitivity that is similar to that obtained with intravenous administration.

Genetic algorithm-based optimization of pulse sequences.

PURPOSE: The performance of pulse sequences in vivo can be limited by fast relaxation rates, magnetic field inhomogeneity, and nonuniform spin excitation. We describe here a method for pulse sequence optimization that uses a stochastic numerical solver that in principle is capable of finding a global optimum. The method provides a simple framework for incorporating any constraint and implementing arbitrarily complex cost functions. Efficient methods for simulating spin dynamics and incorporating frequency selectivity are also described. METHODS: Optimized pulse sequences for polarization transfer between protons and X-nuclei and excitation pulses that eliminate J-coupling modulation were evaluated experimentally using a surface coil on phantoms, and also the detection of hyperpolarized [2-13 C]lactate in vivo in the case of J-coupling modulation-free excitation. RESULTS: The optimized polarization transfer pulses improved the SNR by ~50% with a more than twofold reduction in the B1 field, and J-coupling modulation-free excitation was achieved with a more than threefold reduction in pulse length. CONCLUSION: This process could be used to optimize any pulse when there is a need to improve the uniformity and frequency selectivity of excitation as well as to design new pulses to steer the spin system to any desired achievable state.

Comparison of 13 C MRI of hyperpolarized [1-13 C]pyruvate and lactate with the corresponding mass spectrometry images in a murine lymphoma model.

PURPOSE: To compare carbon-13 (13 C) MRSI of hyperpolarized [1-13 C]pyruvate metabolism in a murine tumor model with mass spectrometric (MS) imaging of the corresponding tumor sections in order to cross validate these metabolic imaging techniques and to investigate the effects of pyruvate delivery and tumor lactate concentration on lactate labeling. METHODS: [1-13 C]lactate images were obtained from tumor-bearing mice, following injection of hyperpolarized [1-13 C]pyruvate, using a single-shot 3D 13 C spectroscopic imaging sequence in vivo and using desorption electrospray ionization MS imaging of the corresponding rapidly frozen tumor sections ex vivo. The images were coregistered, and levels of association were determined by means of Spearman rank correlation and Cohen kappa coefficients as well as linear mixed models. The correlation between [1-13 C]pyruvate and [1-13 C]lactate in the MRS images and between [12 C] and [1-13 C]lactate in the MS images were determined by means of Pearson correlation coefficients. RESULTS: [1-13 C]lactate images generated by MS imaging were significantly correlated with the corresponding MRS images. The correlation coefficient between [1-13 C]lactate and [1-13 C]pyruvate in the MRS images was higher than between [1-13 C]lactate and [12 C]lactate in the MS images. CONCLUSION: The inhomogeneous distribution of labeled lactate observed in the MRS images was confirmed by MS imaging of the corresponding tumor sections. The images acquired using both techniques show that the rate of 13 C label exchange between the injected pyruvate and endogenous tumor lactate pool is more correlated with the rate of pyruvate delivery to the tumor cells and is less affected by the endogenous lactate concentration.

A multi spin echo pulse sequence with optimized excitation pulses and a 3D cone readout for hyperpolarized 13 C imaging.

PURPOSE: Imaging tumor metabolism in vivo using hyperpolarized [1-13 C]pyruvate is a promising technique for detecting disease, monitoring disease progression, and assessing treatment response. However, the transient nature of the hyperpolarization and its depletion following excitation limits the available time for imaging. We describe here a single-shot multi spin echo sequence, which improves on previously reported sequences, with a shorter readout time, isotropic point spread function (PSF), and better signal-to-noise ratio. METHODS: The sequence uses numerically optimized spectrally selective excitation pulses set to the resonant frequencies of pyruvate and lactate and a hyperbolic secant adiabatic refocusing pulse, all applied in the absence of slice selection gradients. The excitation pulses were designed to be resistant to the effects of B0 and B1 field inhomogeneity. The gradient readout uses a 3D cone trajectory composed of 13 cones, all fully refocused and distributed among 7 spin echoes. The maximal gradient amplitude and slew rate were set to 4 G/cm and 20 G/cm/ms, respectively, to demonstrate the feasibility of clinical translation. RESULTS: The pulse sequence gave an isotropic PSF of 2.8 mm. The excitation profiles of the optimized pulses closely matched simulations and a 46.10 ± 0.04% gain in image SNR was observed compared to a conventional Shinnar-Le Roux excitation pulse. The sequence was demonstrated with dynamic imaging of hyperpolarized [1-13 C]pyruvate and [1-13 C]lactate in vivo. CONCLUSION: The pulse sequence was capable of dynamic imaging of hyperpolarized 13 C labeled metabolites in vivo with relatively high spatial and temporal resolution and immunity to system imperfections.

Increasing the sensitivity of hyperpolarized [15 N2 ]urea detection by serial transfer of polarization to spin-coupled protons.

PURPOSE: Hyperpolarized 15 N-labeled molecules have been proposed as imaging agents for investigating tissue perfusion and pH. However, the sensitivity of direct 15 N detection is limited by the isotope's low gyromagnetic ratio. Sensitivity can be increased by transferring 15 N hyperpolarization to spin-coupled protons provided that there is not significant polarization loss during transfer. However, complete polarization transfer would limit the temporal window for imaging to the order of the proton T1 (2-3 s). To exploit the long T1 offered by storing polarization in 15 N and the higher sensitivity of 1 H detection, we have developed a pulse sequence for partial polarization transfer. METHODS: A polarization transfer pulse sequence was modified to allow partial polarization transfer, as is required for dynamic measurements, and that can be implemented with inhomogeneous B1 fields, as is often the case in vivo. The sequence was demonstrated with dynamic spectroscopy and imaging measurements with [15 N2 ]urea. RESULTS: When compared to direct 15 N detection, the sequence increased the signal-to-noise ratio (SNR) by a factor of 1.72 ± 0.25, where both experiments depleted ~20% of the hyperpolarization (>10-fold when 100% of the hyperpolarization is used). Simulations with measured cross relaxation rates showed that this sequence gave up to a 50-fold increase in urea proton polarization when compared to spontaneous polarization transfer via cross relaxation. CONCLUSION: The sequence gave an SNR increase that was close to the theoretical limit and can give a significant SNR benefit when compared to direct 13 C detection of hyperpolarized [13 C]urea.

Imaging Glioblastoma Response to Radiotherapy Using 2H Magnetic Resonance Spectroscopy Measurements of Fumarate Metabolism.

Early detection of tumor cell death in glioblastoma following treatment with chemoradiation has the potential to distinguish between true disease progression and pseudoprogression. Tumor cell death can be detected noninvasively in vivo by imaging the production of [2,3-2H2]malate from [2,3-2H2]fumarate using 2H magnetic resonance (MR) spectroscopic imaging. We show here that 2H MR spectroscopy and spectroscopic imaging measurements of [2,3-2H2]fumarate metabolism can detect tumor cell death in orthotopically implanted glioblastoma models within 48 hours following the completion of chemoradiation. Following the injection of [2,3-2H2]fumarate into tumor-bearing mice, production of [2,3-2H2]malate was measured in a human cell line-derived model and in radiosensitive and radioresistant patient-derived models of glioblastoma that were treated with temozolomide followed by targeted fractionated irradiation. The increase in the [2,3-2H2]malate/[2,3-2H2]fumarate signal ratio posttreatment, which correlated with histologic assessment of cell death, was a more sensitive indicator of treatment response than diffusion-weighted and contrast agent-enhanced 1H MRI measurements, which have been used clinically to detect responses of glioblastoma to chemoradiation. Overall, early detection of glioblastoma cell death using 2H MRI of malate production from fumarate could help improve the clinical evaluation of response to chemoradiation. SIGNIFICANCE: 2H magnetic resonance imaging of labeled fumarate metabolism can detect early evidence of tumor cell death following chemoradiation, meeting a clinical need to reliably detect treatment response in glioblastoma.

Singular value decomposition analysis of back projection operator of maximum likelihood expectation maximization PET image reconstruction.

Background In emission tomography maximum likelihood expectation maximization reconstruction technique has replaced the analytical approaches in several applications. The most important drawback of this iterative method is its linear rate of convergence and the corresponding computational burden. Therefore, simplifications are usually required in the Monte Carlo simulation of the back projection step. In order to overcome these problems, a reconstruction code has been developed with graphical processing unit based Monte Carlo engine which enabled full physical modelling in the back projection. Materials and methods Code performance was evaluated with simulations on two geometries. One is a sophisticated scanner geometry which consists of a dodecagon with inscribed circle radius of 8.7 cm, packed on each side with an array of 39 × 81 LYSO detector pixels of 1.17 mm sided squares, similar to a Mediso nanoScan PET/CT scanner. The other, simplified geometry contains a 38,4mm long interval as a voxel space, detector pixels are assigned in two parallel sections each containing 81 crystals of a size 1.17×1.17 mm. Results We have demonstrated that full Monte Carlo modelling in the back projection step leads to material dependent inhomogeneities in the reconstructed image. The reasons behind this apparently anomalous behaviour was analysed in the simplified system by means of singular value decomposition and explained by different speed of convergence. Conclusions To still take advantage of the higher noise stability of the full physical modelling, a new filtering technique is proposed for convergence acceleration. Some theoretical considerations for the practical implementation and for further development are also presented.