Application of double spin echo spiral chemical shift imaging to rapid metabolic mapping of hyperpolarized [1-¹³C]-pyruvate

A metabolite specific 3D stack-of-spirals bSSFP sequence for improved bicarbonate imaging in hyperpolarized [1-13C]Pyruvate MRI

13C-bicarbonate is a crucial measure of pyruvate oxidation and TCA cycle flux, but is challenging to measure due to its relatively low concentration and thus will greatly benefit from improved signal-to-noise ratio (SNR). To address this, we developed and investigated the feasibility of a 3D stack-of-spirals metabolite-specific balanced steady-state free precession (MS-bSSFP) sequence for improving the SNR and spatial resolution of dynamic 13C-bicarbonate imaging in hyperpolarized [1-13C]pyruvate studies. The bicarbonate MS-bSSFP sequence was evaluated by simulations, phantoms studies, preclinical studies on five rats, brain studies on two healthy volunteers and renal study on one renal cell carcinoma patient. The simulations and phantom results showed that the bicarbonate-specific pulse had minimal perturbation of other metabolites (<1%). In the animal studies, the MS-bSSFP sequence provided an approximately 2.6-3 × improvement in 13C-bicarbonate SNR compared to a metabolite-specific gradient echo (MS-GRE) sequence without altering the bicarbonate or pyruvate kinetics, and the shorter spiral readout in the MS-bSSFP approach reduced blurring. Using the SNR ratio between MS-bSSFP and MS-GRE, the T2 values of bicarbonate and lactate in the rat kidneys were estimated as 0.5 s and 1.1 s, respectively. The in-vivo feasibility of bicarbonate MS-bSSFP sequence was demonstrated in two human brain studies and one renal study. These studies demonstrate the potential of the sequence for in-vivo applications, laying the foundation for future studies to observe this relatively low concentration metabolite with high-quality images and improve measurements of pyruvate oxidation.

Hyperpolarized Carbon (13C) MRI of the Kidney: Experimental Protocol

Alterations in renal metabolism are associated with both physiological and pathophysiologic events. The existing noninvasive analytic tools including medical imaging have limited capability for investigating these processes, which potentially limits current understanding of kidney disease and the precision of its clinical diagnosis. Hyperpolarized 13C MRI is a new medical imaging modality that can capture changes in the metabolic processing of certain rapidly metabolized substrates, as well as changes in kidney function. Here we describe experimental protocols for renal metabolic [1-13C]pyruvate and functional 13C-urea imaging step-by-step. These methods and protocols are useful for investigating renal blood flow and function as well as the renal metabolic status of rodents in vivo under various experimental (patho)physiological conditions.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This experimental protocol is complemented by two separate chapters describing the basic concept and data analysis.

Iterative motion-compensation reconstruction ultra-short TE (iMoCo UTE) for high-resolution free-breathing pulmonary MRI

PURPOSE

To develop a high-scanning efficiency, motion-corrected imaging strategy for free-breathing pulmonary MRI by combining an iterative motion-compensation reconstruction with a ultrashort echo time (UTE) acquisition called iMoCo UTE.

METHODS

An optimized golden-angle ordering radial UTE sequence was used to continuously acquire data for 5 minutes. All readouts were grouped to different respiratory motion states based on self-navigator signals, and then motion-resolved data was reconstructed by XD golden-angle radial sparse parallel reconstruction. One state from the motion-resolved images was selected as a reference, and then motion fields from the other states to the reference were derived via nonrigid registration. Finally, all motion-resolved data and motion fields were reconstructed by using an iterative motion-compensation (MoCo) reconstruction with a total generalized variation sparse constraint.

RESULTS

The iMoCo UTE strategy was evaluated in volunteers and nonsedated pediatric patient (4-6 years old) studies. Images reconstructed with iMoCo UTE provided sharper anatomical lung structures and higher apparent SNR and contrast-to-noise ratio compared to using other motion-correction strategies, such as soft-gating, motion-resolved reconstruction, and nonrigid MoCo. iMoCo UTE also showed promising results in an infant study.

CONCLUSION

The proposed iMoCo UTE combines self-navigation, motion modeling, and a compressed sensing reconstruction to increase scan efficiency and SNR and to reduce respiratory motion in lung MRI. This proposed strategy shows improvements in free-breathing lung MRI scans, especially in very challenging application situations such as pediatric MRI studies.

Pulse sequence considerations for quantification of pyruvate-to-lactate conversion kPL in hyperpolarized 13 C imaging

Hyperpolarized ¹³ C MRI takes advantage of the unprecedented 50 000-fold signal-to-noise ratio enhancement to interrogate cancer metabolism in patients and animals. It can measure the pyruvate-to-lactate conversion rate, k_PL , a metabolic biomarker of cancer aggressiveness and progression. Therefore, it is crucial to evaluate k_PL reliably. In this study, three sequence components and parameters that modulate k_PL estimation were identified and investigated in model simulations and through in vivo animal studies using several specifically designed pulse sequences. These factors included a magnetization spoiling effect due to RF pulses, a crusher gradient-induced flow suppression, and intrinsic image weightings due to relaxation. Simulation showed that the RF-induced magnetization spoiling can be substantially improved using an inputless k_PL fitting. In vivo studies found a significantly higher apparent k_PL with an additional gradient that leads to flow suppression (k_{PL,FID-Delay,Crush} /k_PL,FID-Delay  = 1.37 ± 0.33, P < 0.01, N = 6), which agrees with simulation outcomes (12.5% k_PL error with Δv = 40 cm/s), indicating that the gradients predominantly suppressed flowing pyruvate spins. Significantly lower k_PL was found using a delayed free induction decay (FID) acquisition versus a minimum-T_E version (k_PL,FID-Delay /k_PL,FID  = 0.67 ± 0.09, P < 0.01, N = 5), and the lactate peak had broader linewidth than pyruvate (Δω_lactate /Δω_pyruvate  = 1.32 ± 0.07, P < 0.000 01, N = 13). This illustrated that lactate's T₂ *, shorter than that of pyruvate, can affect calculated k_PL values. We also found that an FID sequence yielded significantly lower k_PL versus a double spin-echo sequence that includes spin-echo spoiling, flow suppression from crusher gradients, and more T₂ weighting (k_PL,DSE /k_PL,FID  = 2.40 ± 0.98, P < 0.0001, N = 7). In summary, the pulse sequence, as well as its interaction with pharmacokinetics and the tissue microenvironment, can impact and be optimized for the measurement of k_PL . The data acquisition and analysis pipelines can work synergistically to provide more robust and reproducible k_PL measures for future preclinical and clinical studies.

Fast data acquisition techniques in magnetic resonance spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) is an important technique for assessing the spatial variation of metabolites in vivo. The long scan times in MRSI limit clinical applicability due to patient discomfort, increased costs, motion artifacts, and limited protocol flexibility. Faster acquisition strategies can address these limitations and could potentially facilitate increased adoption of MRSI into routine clinical protocols with minimal addition to the current anatomical and functional acquisition protocols in terms of imaging time. Not surprisingly, a lot of effort has been devoted to the development of faster MRSI techniques that aim to capture the same underlying metabolic information (relative metabolite peak areas and spatial distribution) as obtained by conventional MRSI, in greatly reduced time. The gain in imaging time results, in some cases, in a loss of signal-to-noise ratio and/or in spatial and spectral blurring. This review examines the current techniques and advances in fast MRSI in two and three spatial dimensions and their applications. This review categorizes the acceleration techniques according to their strategy for acquisition of the k-space. Techniques such as fast/turbo-spin echo MRSI, echo-planar spectroscopic imaging, and non-Cartesian MRSI effectively cover the full k-space in a more efficient manner per T_R . On the other hand, techniques such as parallel imaging and compressed sensing acquire fewer k-space points and employ advanced reconstruction algorithms to recreate the spatial-spectral information, which maintains statistical fidelity in test conditions (ie no statistically significant differences on voxel-wise comparisions) with the fully sampled data. The advantages and limitations of each state-of-the-art technique are reviewed in detail, concluding with a note on future directions and challenges in the field of fast spectroscopic imaging.

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

Hyperpolarized MRI with 13 C-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 13 C 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-13 C]pyruvate. [1-13 C] pyruvate and [1-13 C] 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 T2 relaxation times. The pulse sequences described here produced hyperpolarized 13 C images with improved resolution and signal-to-noise ratio when compared with similar sequences described previously.

Dynamic diffusion-weighted hyperpolarized 13 C imaging based on a slice-selective double spin echo sequence for measurements of cellular transport

PURPOSE

To develop a pulse sequence to dynamically measure the ADC of hyperpolarized substrates during their perfusion, metabolic conversion, and transport.

METHODS

We proposed a slice-selective double spin echo sequence for dynamic hyperpolarized ¹³ C diffusion-weighted imaging. The proposed pulse sequence was optimized for a high field preclinical scanner through theoretical analysis and simulation. The performance of the method was compared to non-slice-selective double spin echo via in vivo studies. We also validated the sequence for dynamic ADC measurement in both phantom studies and transgenic mouse model of prostate cancer studies.

RESULTS

The optimized pulse sequence outperforms the traditional sequence with smaller saturation effects on the magnetization of hyperpolarized compounds that allowed more dynamic imaging frames covering a longer imaging time window. In pre-clinical studies (N = 8), the dynamic hyperpolarized lactate ADC maps of 6 studies in the prostate tumors showed an increase measured ADC over time, which might be related to lactate efflux from the tumor cells.

CONCLUSIONS

The proposed sequence was validated and shown to improve dynamic diffusion weighted imaging compared to the traditional double spin echo sequence, providing ADC maps of lactate through time.

Development of a Symmetric Echo-Planar Spectroscopy Imaging Framework for Hyperpolarized 13C Imaging in a Clinical PET/MR Scanner

Here, we developed a symmetric echo-planar spectroscopic imaging (EPSI) sequence for hyperpolarized 13C imaging on a clinical hybrid positron emission tomography/magnetic resonance imaging system. The pulse sequence uses parallel reconstruction pipelines to separately reconstruct data from odd-and-even gradient echoes to reduce artifacts from gradient imbalances. The ramp-sampled data in the spatiotemporal frequency space are regridded to compensate for the chemical-shift displacements. Unaliasing of nonoverlapping peaks outside of the sampled spectral width was performed to double the effective spectral width. The sequence was compared with conventional phase-encoded chemical-shift imaging (CSI) in phantoms, and it was evaluated in a canine cancer patient with ameloblastoma after injection of hyperpolarized [1-13C]pyruvate. The relative signal-to-noise ratio of EPSI with respect to CSI was 0.88, which is consistent with the decrease in sampling efficiency due to ramp sampling. Data regridding in the spatiotemporal frequency space significantly reduced spatial blurring compared with direct fast Fourier transform. EPSI captured the spatial distributions of both metabolites and their temporal dynamics in vivo with an in-plane spatial resolution of 5 × 9 mm2 and a temporal resolution of 3 seconds. Significantly higher spatial and temporal resolution for delineating anatomical structures in vivo was achieved for EPSI metabolic maps than for CSI maps, which suffered spatiotemporal blurring. The EPSI sequence showed promising results in terms of short acquisition time and sufficient spectral bandwidth of 500 Hz, allowing to adjust the trade-off between signal-to-noise ratio and encoding speed.

3D hyperpolarized C-13 EPI with calibrationless parallel imaging

With the translation of metabolic MRI with hyperpolarized 13C agents into the clinic, imaging approaches will require large volumetric FOVs to support clinical applications. Parallel imaging techniques will be crucial to increasing volumetric scan coverage while minimizing RF requirements and temporal resolution. Calibrationless parallel imaging approaches are well-suited for this application because they eliminate the need to acquire coil profile maps or auto-calibration data. In this work, we explored the utility of a calibrationless parallel imaging method (SAKE) and corresponding sampling strategies to accelerate and undersample hyperpolarized 13C data using 3D blipped EPI acquisitions and multichannel receive coils, and demonstrated its application in a human study of [1-13C]pyruvate metabolism.

A referenceless Nyquist ghost correction workflow for echo planar imaging of hyperpolarized [1-13 C]pyruvate and [1-13 C]lactate

Single-shot echo planar imaging (EPI), which allows an image to be acquired using a single excitation pulse, is used widely for imaging the metabolism of hyperpolarized 13 C-labelled metabolites in vivo as the technique is rapid and minimizes the depletion of the hyperpolarized signal. However, EPI suffers from Nyquist ghosting, which normally is corrected for by acquiring a reference scan. In a dynamic acquisition of a series of images, this results in the sacrifice of a time point if the reference scan involves a full readout train with no phase encoding. This time penalty is negligible if an integrated navigator echo is used, but at the cost of a lower signal-to-noise ratio (SNR) as a result of prolonged T2 * decay. We describe here a workflow for hyperpolarized 13 C EPI that requires no reference scan. This involves the selection of a ghost-containing background from a 13 C image of a single metabolite at a single time point, the identification of phase correction coefficients that minimize signal in the selected area, and the application of these coefficients to images acquired at all time points and from all metabolites. The workflow was compared in phantom experiments with phase correction using a 13 C reference scan, and yielded similar results in situations with a regular field of view (FOV), a restricted FOV and where there were multiple signal sources. When compared with alternative phase correction methods, the workflow showed an SNR benefit relative to integrated 13 C reference echoes (>15%) or better ghost removal relative to a 1 H reference scan. The residual ghosting in a slightly de-shimmed B0 field was 1.6% using the proposed workflow and 3.8% using a 1 H reference scan. The workflow was implemented with a series of dynamically acquired hyperpolarized [1-13 C]pyruvate and [1-13 C]lactate images in vivo, resulting in images with no observable ghosting and which were quantitatively similar to images corrected using a 13 C reference scan.

High resolution hyperpolarized 13 C MRSI using SPICE at 9.4T

PURPOSE

To test the feasibility of using the SPICE (SPectroscopic Imaging by exploiting spatiospectral CorrElation) technique, which uses the partial separability of spectroscopic data, for high resolution hyperpolarized (HP) ¹³ C spectroscopic imaging.

METHODS

Numerical simulations were performed to investigate the impact of transient HP signals on SPICE reconstruction. Furthermore, spectroscopic imaging exams from SPICE and conventional EPSI (echo-planar spectroscopic imaging) were simulated for comparison. For in vivo experiments, HP ¹³ C SPICE was performed in a mouse kidney by means of the injection of HP [1-¹³ C] pyruvate at 9.4T.

RESULTS

The variation of lactate/pyruvate from the simulated SPICE was less than 4% under various factors that affect the transient HP signal, suggesting that the impact is negligible. We found that while HP ¹³ C EPSI was limited to the low signal-to-noise ratio (SNR) of lactate, these limitations were mitigated through HP ¹³ C SPICE, facilitating the improved SNR of lactate and the distinction of tissues. Acquisition of a high resolution HP ¹³ C spectroscopic image was possible for the in vivo experiments. With the fine structural information, the acquired image showed higher signal of pyruvate and lactate in the renal cortices than in the medullas, which is known to be attributed to higher activity of lactate dehydrogenase.

CONCLUSION

The feasibility of HP ¹³ C SPICE was investigated. Simulation studies were conducted and in vivo experiments were performed in the mouse kidney at 9.4T. Results confirmed that a high resolution HP ¹³ C spectroscopic image with adequate spectral resolution can be obtained. Magn Reson Med 80:703-710, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

Flow-suppressed hyperpolarized 13 C chemical shift imaging using velocity-optimized bipolar gradient in mouse liver tumors at 9.4 T

PURPOSE

To optimize and investigate the influence of bipolar gradients for flow suppression in metabolic quantification of hyperpolarized ¹³ C chemical shift imaging (CSI) of mouse liver at 9.4 T.

METHODS

The trade-off between the amount of flow suppression using bipolar gradients and T2* effect from static spins was simulated. A free induction decay CSI sequence with alternations between the flow-suppressed and non-flow-suppressed acquisitions for each repetition time was developed and was applied to liver tumor-bearing mice via injection of hyperpolarized [1-¹³ C] pyruvate.

RESULTS

The in vivo results from flow suppression using the velocity-optimized bipolar gradient were comparable with the simulation results. The vascular signal was adequately suppressed and signal loss in stationary tissue was minimized. Application of the velocity-optimized bipolar gradient to tumor-bearing mice showed reduction in the vessel-derived pyruvate signal contamination, and the average lactate/pyruvate ratio increased by 0.095 (P < 0.05) in the tumor region after flow suppression.

CONCLUSION

Optimization of the bipolar gradient is essential because of the short ¹³ C T2* and high signal in venous flow in the mouse liver. The proposed velocity-optimized bipolar gradient can suppress the vascular signal, minimizing T2*-related signal loss in stationary tissues at 9.4 T. Magn Reson Med 78:1674-1682, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Imaging Renal Urea Handling in Rats at Millimeter Resolution using Hyperpolarized Magnetic Resonance Relaxometry

In vivo spin spin relaxation time (T2) heterogeneity of hyperpolarized [13C,15N2]urea in the rat kidney was investigated. Selective quenching of the vascular hyperpolarized 13C signal with a macromolecular relaxation agent revealed that a long-T2 component of the [13C,15N2]urea signal originated from the renal extravascular space, thus allowing the vascular and renal filtrate contrast agent pools of the [13C,15N2]urea to be distinguished via multi-exponential analysis. The T2 response to induced diuresis and antidiuresis was performed with two imaging agents: hyperpolarized [13C,15N2]urea and a control agent hyperpolarized bis-1,1-(hydroxymethyl)-1-13C-cyclopropane-2H8. Large T2 increases in the inner-medullar and papilla were observed with the former agent and not the latter during antidiuresis. Therefore, [13C,15N2]urea relaxometry is sensitive to two steps of the renal urea handling process: glomerular filtration and the inner-medullary urea transporter (UT)-A1 and UT-A3 mediated urea concentrating process. Simple motion correction and subspace denoising algorithms are presented to aid in the multi exponential data analysis. Furthermore, a T2-edited, ultra long echo time sequence was developed for sub-2 mm3 resolution 3D encoding of urea by exploiting relaxation differences in the vascular and filtrate pools.

Imaging Renal Urea Handling in Rats at Millimeter Resolution using Hyperpolarized Magnetic Resonance Relaxometry

In vivo spin spin relaxation time (T₂) heterogeneity of hyperpolarized [¹³C,¹⁵N₂]urea in the rat kidney was investigated. Selective quenching of the vascular hyperpolarized ¹³C signal with a macromolecular relaxation agent revealed that a long-T₂ component of the [¹³C,¹⁵N₂]urea signal originated from the renal extravascular space, thus allowing the vascular and renal filtrate contrast agent pools of the [¹³C,¹⁵N₂]urea to be distinguished via multi-exponential analysis. The T₂ response to induced diuresis and antidiuresis was performed with two imaging agents: hyperpolarized [¹³C,¹⁵N₂]urea and a control agent hyperpolarized bis-1,1-(hydroxymethyl)-1-¹³C-cyclopropane-²H₈. Large T₂ increases in the inner-medullar and papilla were observed with the former agent and not the latter during antidiuresis. Therefore, [¹³C,¹⁵N₂]urea relaxometry is sensitive to two steps of the renal urea handling process: glomerular filtration and the inner-medullary urea transporter (UT)-A1 and UT-A3 mediated urea concentrating process. Simple motion correction and subspace denoising algorithms are presented to aid in the multi exponential data analysis. Furthermore, a T₂-edited, ultra long echo time sequence was developed for sub-2 mm³ resolution 3D encoding of urea by exploiting relaxation differences in the vascular and filtrate pools.

Rapid in vivo apparent diffusion coefficient mapping of hyperpolarized (13) C metabolites

PURPOSE

Hyperpolarized (13) C magnetic resonance allows for the study of real-time metabolism in vivo, including significant hyperpolarized (13) C lactate production in many tumors. Other studies have shown that aggressive and highly metastatic tumors rapidly transport lactate out of cells. Thus, the ability to not only measure the production of hyperpolarized (13) C lactate but also understand its compartmentalization using diffusion-weighted MR will provide unique information for improved tumor characterization.

METHODS

We used a bipolar, pulsed-gradient, double spin echo imaging sequence to rapidly generate diffusion-weighted images of hyperpolarized (13) C metabolites. Our methodology included a simultaneously acquired B1 map to improve apparent diffusion coefficient (ADC) accuracy and a diffusion-compensated variable flip angle scheme to improve ADC precision.

RESULTS

We validated this sequence and methodology in hyperpolarized (13) C phantoms. Next, we generated ADC maps of several hyperpolarized (13) C metabolites in a normal rat, rat brain tumor, and prostate cancer mouse model using both preclinical and clinical trial-ready hardware.

CONCLUSION

ADC maps of hyperpolarized (13) C metabolites provide information about the localization of these molecules in the tissue microenvironment. The methodology presented here allows for further studies to investigate ADC changes due to disease state that may provide unique information about cancer aggressiveness and metastatic potential.

Single voxel localization for dynamic hyperpolarized (13)C MR spectroscopy

The PRESS technique has been widely used to achieve voxel localization for in vivo(1)H MRS acquisitions. However, for dynamic hyperpolarized (13)C MRS experiments, the transition bands of the refocusing pulses may saturate the pre-polarized substrate spins flowing into the voxel. This limitation may be overcome by designing refocusing pulses that do not perturb the resonance of the hyperpolarized substrate, but selectively refocuses the spins of the metabolic products. In this study, a PRESS pulse sequence incorporating spectral-spatial refocusing pulses that have a stop band ('notch') at the substrate resonance is tested in vivo using hyperpolarized [1-(13)C]pyruvate. Higher metabolite SNR was observed in experiments using the spectral-spatial refocusing pulses as compared to conventional refocusing pulses.

A 2DRF pulse sequence for bolus tracking in hyperpolarized 13C imaging

PURPOSE

A novel application of two-dimensional (2D) spatially selective radiofrequency (2DRF) excitation pulses in hyperpolarized 13C imaging is proposed for monitoring the bolus injection with highly efficient sampling of the initially polarized substrate, thus leaving more polarization available for detection of the subsequently generated metabolic products.

METHODS

A 2DRF pulse was designed with a spiral trajectory and conventional clinical gradient performance. To demonstrate the ability of our 2DRF bolus tracking pulse sequence, hyperpolarized [1-(13)ruvate in vivo imaging experiments were performed in normal rats, with a comparison to 1DRF excitation pulses.

RESULTS

Our designed 2DRF pulse was able to rapidly and efficiently monitor the injected bolus dynamics in vivo, with an 8-fold enhanced time resolution in comparison with 1DRF in our experimental settings. When applied at the pyruvate frequency for bolus tracking, our 2DRF pulse demonstrated reduced saturation of the hyperpolarization for the substrate and metabolic products compared to a 1DRF pulse, while being immune to ±0.5 ppm magnetic field inhomogeneity at 3T.

CONCLUSION

2DRF pulses in hyperpolarized 13C imaging can be used to efficiently monitor the bolus injection with reduced hyperpolarization saturation compared to 1DRF pulses. The parameters of our design are based on clinical scanner limits, which allows for rapid translation to human studies.

Comparison of acquisition schemes for hyperpolarised ¹³C imaging

The aim of this study was to characterise and compare widely used acquisition strategies for hyperpolarised (13)C imaging. Free induction decay chemical shift imaging (FIDCSI), echo-planar spectroscopic imaging (EPSI), IDEAL spiral chemical shift imaging (ISPCSI) and spiral chemical shift imaging (SPCSI) sequences were designed for two different regimes of spatial resolution. Their characteristics were studied in simulations and in tumour-bearing rats after injection of hyperpolarised [1-(13)C]pyruvate on a clinical 3-T scanner. Two or three different sequences were used on the same rat in random order for direct comparison. The experimentally obtained lactate signal-to-noise ratio (SNR) in the tumour matched the simulations. Differences between the sequences were mainly found in the encoding efficiency, gradient demand and artefact behaviour. Although ISPCSI and SPCSI offer high encoding efficiencies, these non-Cartesian trajectories are more prone than EPSI and FIDCSI to artefacts from various sources. If the encoding efficiency is sufficient for the desired application, EPSI has been proven to be a robust choice. Otherwise, faster spiral acquisition schemes are recommended. The conclusions found in this work can be applied directly to clinical applications.

Application of flow sensitive gradients for improved measures of metabolism using hyperpolarized (13) c MRI

PURPOSE

To develop the use of bipolar gradients to suppress partial-volume and flow-related artifacts from macrovascular, hyperpolarized spins.

THEORY AND METHODS

Digital simulations were performed over a range of spatial resolutions and gradient strengths to determine the optimal bipolar gradient strength and duration to suppress flowing spins while minimizing signal loss from static tissue. In vivo experiments were performed to determine the efficacy of this technique to suppress vascular signal in the study of hyperpolarized [1-(13)C]pyruvate renal metabolism.

RESULTS

Digital simulations showed that in the absence of bipolar gradients, partial-volume artifacts from the vasculature were still present, causing underestimation of the apparent reaction rate of pyruvate to lactate (kP). The addition of a bipolar gradient with b = 32 s/mm(2) sufficiently suppressed the vascular signal without a substantial decrease in signal from static tissue. In vivo results corroborate digital simulations, with similar peak lactate signal to noise ratio (SNR) but substantially different kP in the presence of bipolar gradients.

CONCLUSION

The proposed approach suppresses signal from flowing spins while minimizing signal loss from static tissue, removing contaminating signal from the vasculature and increasing kinetic modeling accuracy without substantially sacrificing SNR or temporal resolution.

In vivo measurement of apparent diffusion coefficients of hyperpolarized ¹³C-labeled metabolites

The combination of hyperpolarized MRS with diffusion weighting (dw) allows for determination of the apparent diffusion coefficient (ADC), which is indicative of the intra- or extracellular localization of the metabolite. Here, a slice-selective pulsed-gradient spin echo sequence was implemented to acquire a series of dw spectra from rat muscle in vivo to determine the ADCs of multiple metabolites after a single injection of hyperpolarized [1- ¹³C]pyruvate. An optimal control optimized universal-rotation pulse was used for refocusing to minimize signal loss caused by B1 imperfections. Non-dw spectra were acquired interleaved with the dw spectra and these were used to correct for signal decay during the acquisition as a result of T1 decay, pulse imperfections, flow etc. The data showed that the ADC values for [1- ¹³C]lactate (0.4-0.7 µm² /ms) and [1- ¹³C]alanine (0.4-0.9 µm² /ms) were about a factor of two lower than the ADC of [1- ¹³C]pyruvate (1.1-1.5 µm²/ms). This indicates a more restricted diffusion space for the former two metabolites consistent with lactate and alanine being intracellular. The higher ADC for pyruvate (similar to the proton ADC) reflected that the injected substance was not confined inside the muscle cells but also present extracellular.

Kinetic and perfusion modeling of hyperpolarized (13)C pyruvate and urea in cancer with arbitrary RF flip angles

The accurate detection and characterization of cancerous tissue is still a major problem for the clinical management of individual cancer patients and for monitoring their response to therapy. MRI with hyperpolarized agents is a promising technique for cancer characterization because it can non-invasively provide a local assessment of the tissue metabolic profile. In this work, we measured the kinetics of hyperpolarized [1-(13)C] pyruvate and (13)C-urea in prostate and liver tumor models using a compressed sensing dynamic MRSI method. A kinetic model fitting method was developed that incorporated arbitrary RF flip angle excitation and measured a pyruvate to lactate conversion rate, Kpl, of 0.050 and 0.052 (1/s) in prostate and liver tumors, respectively, which was significantly higher than Kpl in healthy tissues [Kpl =0.028 (1/s), P<0.001]. Kpl was highly correlated to the total lactate to total pyruvate signal ratio (correlation coefficient =0.95). We additionally characterized the total pyruvate and urea perfusion, as in cancerous tissue there is both existing vasculature and neovascularization as different kinds of lesions surpass the normal blood supply, including small circulation disturbance in some of the abnormal vessels. A significantly higher perfusion of pyruvate (accounting for conversion to lactate and alanine) relative to urea perfusion was seen in cancerous tissues (liver cancer and prostate cancer) compared to healthy tissues (P<0.001), presumably due to high pyruvate uptake in tumors.

Perfusion and diffusion sensitive 13C stimulated-echo MRSI for metabolic imaging of cancer

Metabolic imaging with hyperpolarized [1-(13)C]-pyruvate can rapidly probe tissue metabolic profiles in vivo and has been shown to provide cancer imaging biomarkers for tumor detection, progression, and response to therapy. This technique uses a bolus injection followed by imaging within 1-2 minutes. The observed metabolites include vascular components and their generation is also influenced by cellular transport. These factors complicate image interpretation, especially since [1-(13)C]lactate, a metabolic product that is a biomarker of cancer, is also produced by red blood cells. It would be valuable to understand the distribution of metabolites between the vasculature, interstitial space, and intracellular compartments. The purpose of this study was to better understand this compartmentalization by using a perfusion and diffusion-sensitive stimulated-echo acquisition mode (STEAM) MRSI acquisition method tailored to hyperpolarized substrates. Our results in mouse models showed that among metabolites, the injected substrate (13)C-pyruvate had the largest vascular fraction overall while (13)C-alanine had the smallest vascular fraction. We observed a larger vascular fraction of pyruvate and lactate in the kidneys and liver when compared to back muscle and prostate tumor tissue. Our data suggests that (13)C-lactate in prostate tumor tissue voxels was the most abundant labeled metabolite intracellularly. This was shown in STEAM images that highlighted abnormal cancer cell metabolism and suppressed vascular (13)C metabolite signals.

Spin echo measurements of the extravasation and tumor cell uptake of hyperpolarized [1-(13) C]lactate and [1-(13) C]pyruvate

PURPOSE

To assess the blood-tissue distribution of hyperpolarized (13) C-labeled molecules in vivo.

METHODS

Spin-echo experiments with simultaneous acquisition of the free induction decay (FID) signal following the excitation pulse and the spin-echo signal, were used to monitor hyperpolarized [1-(13) C]lactate, [1-(13) C]pyruvate, and the perfusion marker, [(13) C]HP001, following their intravenous injection into tumor-bearing mice. Apparent T2 relaxation times and diffusion coefficients were also measured.

RESULTS

An increasing tumor echo/FID ratio was observed for all three molecules, which could be explained by their extravasation into the tumor interstitial space, where T2 relaxation times were longer and diffusion coefficients smaller. Inhibition of the monocarboxylate transporter, which decreased by 40% the label exchange between pyruvate and lactate, reduced the increase in the echo/FID ratio for pyruvate and lactate, but not for HP001, demonstrating that some of the increase in the echo/FID ratio was due to cell uptake of lactate and pyruvate. The different relaxation and diffusion behavior of the intravascular and extravascular signals affected measurements of the apparent label exchange rate constants.

CONCLUSION

Simultaneous collection of both FID and echo signals can provide information on cell uptake thus giving further insight into the kinetics of hyperpolarized (13) C label exchange. Care is needed when comparing exchange rate constants determined in spin-echo-based studies.

How the signal-to-noise ratio influences hyperpolarized 13C dynamic MRS data fitting and parameter estimation

MRS of hyperpolarized (13) C-labeled compounds represents a promising technique for in vivo metabolic studies. However, robust quantification and metabolic modeling are still important areas of investigation. In particular, time and spatial resolution constraints may lead to the analysis of MRS signals with low signal-to-noise ratio (SNR). The relationship between SNR and the precision of quantitative analysis for the evaluation of the in vivo kinetic behavior of metabolites is unknown. In this article, this topic is addressed by Monte Carlo simulations, covering the problem of MRS signal model parameter estimation, with strong emphasis on the peak amplitude and kinetic model parameters. The results of Monte Carlo simulation were confirmed by in vivo experiments on medium-sized animals injected with hyperpolarized [1-(13) C]pyruvate. The results of this study may be useful for the establishment of experimental planning and for the optimization of kinetic model estimation as a function of the SNR value.