High-resolution hyperpolarized metabolic imaging of the rat heart using k-t PCA and k-t SPARSE

Model-constrained reconstruction accelerated with Fourier-based undersampling for hyperpolarized [1-13 C] pyruvate imaging

PURPOSE

Model-constrained reconstruction with Fourier-based undersampling (MoReFUn) is introduced to accelerate the acquisition of dynamic MRI using hyperpolarized [1-13 C]-pyruvate.

METHODS

The MoReFUn method resolves spatial aliasing using constraints introduced by a pharmacokinetic model that describes the signal evolution of both pyruvate and lactate. Acceleration was evaluated on three single-channel data sets: a numerical digital phantom that is used to validate the accuracy of reconstruction and model parameter restoration under various SNR and undersampling ratios, prospectively and retrospectively sampled data of an in vitro dynamic multispectral phantom, and retrospectively undersampled imaging data from a prostate cancer patient to test the fidelity of reconstructed metabolite time series.

RESULTS

All three data sets showed successful reconstruction using MoReFUn. In simulation and retrospective phantom data, the restored time series of pyruvate and lactate maintained the image details, and the mean square residual error of the accelerated reconstruction increased only slightly (< 10%) at a reduction factor up to 8. In prostate data, the quantitative estimation of the conversion-rate constant of pyruvate to lactate was achieved with high accuracy of less than 10% error at a reduction factor of 2 compared with the conversion rate derived from unaccelerated data.

CONCLUSION

The MoReFUn technique can be used as an effective and reliable imaging acceleration method for metabolic imaging using hyperpolarized [1-13 C]-pyruvate.

Hyperpolarized Metabolic MRI-Acquisition, Reconstruction, and Analysis Methods

Hyperpolarized metabolic MRI with 13C-labeled agents has emerged as a powerful technique for in vivo assessments of real-time metabolism that can be used across scales of cells, tissue slices, animal models, and human subjects. Hyperpolarized contrast agents have unique properties compared to conventional MRI scanning and MRI contrast agents that require specialized imaging methods. Hyperpolarized contrast agents have a limited amount of available signal, irreversible decay back to thermal equilibrium, bolus injection and perfusion kinetics, cellular uptake and metabolic conversion kinetics, and frequency shifts between metabolites. This article describes state-of-the-art methods for hyperpolarized metabolic MRI, summarizing data acquisition, reconstruction, and analysis methods in order to guide the design and execution of studies.

Accelerated MR spectroscopic imaging-a review of current and emerging techniques

Over more than 30 years in vivo MR spectroscopic imaging (MRSI) has undergone an enormous evolution from theoretical concepts in the early 1980s to the robust imaging technique that it is today. The development of both fast and efficient sampling and reconstruction techniques has played a fundamental role in this process. State-of-the-art MRSI has grown from a slow purely phase-encoded acquisition technique to a method that today combines the benefits of different acceleration techniques. These include shortening of repetition times, spatial-spectral encoding, undersampling of k-space and time domain, and use of spatial-spectral prior knowledge in the reconstruction. In this way in vivo MRSI has considerably advanced in terms of spatial coverage, spatial resolution, acquisition speed, artifact suppression, number of detectable metabolites and quantification precision. Acceleration not only has been the enabling factor in high-resolution whole-brain 1 H-MRSI, but today is also common in non-proton MRSI (31 P, 2 H and 13 C) and applied in many different organs. In this process, MRSI techniques had to constantly adapt, but have also benefitted from the significant increase of magnetic field strength boosting the signal-to-noise ratio along with high gradient fidelity and high-density receive arrays. In combination with recent trends in image reconstruction and much improved computation power, these advances led to a number of novel developments with respect to MRSI acceleration. Today MRSI allows for non-invasive and non-ionizing mapping of the spatial distribution of various metabolites' tissue concentrations in animals or humans, is applied for clinical diagnostics and has been established as an important tool for neuro-scientific and metabolism research. This review highlights the developments of the last five years and puts them into the context of earlier MRSI acceleration techniques. In addition to 1 H-MRSI it also includes other relevant nuclei and is not limited to certain body regions or specific applications.

Cardiac T 2 ∗ measurement of hyperpolarized 13 C metabolites using metabolite-selective multi-echo spiral imaging

PURPOSE

Noninvasive imaging with hyperpolarized (HP) pyruvate can capture in vivo cardiac metabolism. For proper quantification of the metabolites and optimization of imaging parameters, understanding MR characteristics such as T 2 ∗ s of the HP signals is critical. This study is to measure in vivo cardiac T 2 ∗ s of HP [1-¹³ C]pyruvate and the products in rodents and humans.

METHODS

A dynamic ¹³ C multi-echo spiral imaging sequence that acquires [¹³ C]bicarbonate, [1-¹³ C]lactate, and [1-¹³ C]pyruvate images in an interleaved manner was implemented for a clinical 3 Tesla system. T 2 ∗ of each metabolite was calculated from the multi-echo images by fitting the signal decay of each region of interest mono-exponentially. The performance of measuring T 2 ∗ using the sequence was first validated using a ¹³ C phantom and then with rodents following a bolus injection of HP [1-¹³ C]pyruvate. In humans, T 2 ∗ of each metabolite was calculated for left ventricle, right ventricle, and myocardium.

RESULTS

Cardiac T 2 ∗ s of HP [1-¹³ C]pyruvate, [1-¹³ C]lactate, and [¹³ C]bicarbonate in rodents were measured as 24.9 ± 5.0, 16.4 ± 4.7, and 16.9 ± 3.4 ms, respectively. In humans, T 2 ∗ of [1-¹³ C]pyruvate was 108.7 ± 22.6 ms in left ventricle and 129.4 ± 8.9 ms in right ventricle. T 2 ∗ of [1-¹³ C]lactate was 40.9 ± 8.3, 44.2 ± 5.5, and 43.7 ± 9.0 ms in left ventricle, right ventricle, and myocardium, respectively. T 2 ∗ of [¹³ C]bicarbonate in myocardium was 64.4 ± 2.5 ms. The measurements were reproducible and consistent over time after the pyruvate injection.

CONCLUSION

The proposed metabolite-selective multi-echo spiral imaging sequence reliably measures in vivo cardiac T 2 ∗ s of HP [1-¹³ C]pyruvate and products.

Correlation between lactate dehydrogenase/pyruvate dehydrogenase activities ratio and tissue pH in the perfused mouse heart: A potential noninvasive indicator of cardiac pH provided by hyperpolarized magnetic resonance

Cardiovascular diseases account for more than 30% of all deaths worldwide and many could be ameliorated with early diagnosis. Current cardiac imaging modalities can assess blood flow, heart anatomy and mechanical function. However, for early diagnosis and improved treatment, further functional biomarkers are needed. One such functional biomarker could be the myocardium pH. Although tissue pH is already determinable via MR techniques, and has been since the early 1990s, it remains elusive to use practically. The objective of this study was to explore the possibility to evaluate cardiac pH noninvasively, using in-cell enzymatic rates of hyperpolarized [1-¹³ C]pyruvate metabolism (ie, moles of product produced per unit time) determined directly in real time using magnetic resonance spectroscopy in a perfused mouse heart model. As a gold standard for tissue pH we used ³¹ P spectroscopy and the chemical shift of the inorganic phosphate (Pi) signal. The nonhomogenous pH distribution of the perfused heart was analyzed using a multi-parametric analysis of this signal, thus taking into account the heterogeneous nature of this characteristic. As opposed to the signal ratio of hyperpolarized [¹³ C]bicarbonate to [¹³ CO₂ ], which has shown correlation to pH in other studies, we investigated here the ratio of two intracellular enzymatic rates: lactate dehydrogenase (LDH) and pyruvate dehydrogenase (PDH), by way of determining the production rates of [1-¹³ C]lactate and [¹³ C]bicarbonate, respectively. The enzyme activities determined here are intracellular, while the pH determined using the Pi signal may contain an extracellular component, which could not be ruled out. Nevertheless, we report a strong correlation between the tissue pH and the LDH/PDH activities ratio. This work may pave the way for using the LDH/PDH activities ratio as an indicator of cardiac intracellular pH in vivo, in an MRI examination.

Accelerating CEST MRI in the mouse brain at 9.4 T by exploiting sparsity in the Z-spectrum domain

PURPOSE

Chemical exchange saturation transfer (CEST) is an MR contrast modality offering an enhanced sensitivity for the detection of dilute metabolites with exchangeable protons. Quantitative analysis requires the acquisition of a number of images (usually between 20 and 50 RF offsets) per Z-spectrum, leading to long acquisition times of the order of 5-40 min in practice. In this work, we explore the possibility of employing sparsity in the Z-spectrum domain (irradiation offset dimension) to provide an accelerated acquisition scheme without compromising the quality of reconstructed CEST spectra.

METHOD AND THEORY

Ex vivo and in vivo data were acquired on an experimental, small animal 9.4 T system. Three different reconstruction methods were tested: k-Z SPARSE, k-Z SLR and k-Z principal component analysis (PCA) using retrospective undersampling with net acceleration factors R = 2, 3, 5. The quality of the reconstructed data was compared with respect to CEST spectra and full magnetization transfer ratio (MTR) asymmetry maps.

RESULTS

In both phantom and in vivo data, CEST spectra and the resulting MTR asymmetry maps were reconstructed without significant deterioration in data quality. For a low acceleration factor (R = 2, 3) all applied methods resulted in similar data quality, while for high acceleration factor (R = 5) only k-Z PCA and k-Z SLR could be used. Loss in spatial resolution was observed in reconstruction with k-Z PCA for all acceleration factors. An example of prospective undersampling with acceleration factor R = 3 and k-Z PCA reconstruction demonstrates improved CEST maps when compared with fully sampled data acquisition with either three times longer scan duration or threefold prolonged acquisition window per frequency offset.

CONCLUSION

The acquisition time of CEST spectra can be significantly accelerated by exploiting the sparsity of the Z-domain. For prospective and retrospective analysis using k-Z PCA, an acceleration factor of up to R = 3 can be used without significant loss in data quality.

Acquisition strategies for spatially resolved magnetic resonance detection of hyperpolarized nuclei

Hyperpolarization is an emerging method in magnetic resonance imaging that allows nuclear spin polarization of gases or liquids to be temporarily enhanced by up to five or six orders of magnitude at clinically relevant field strengths and administered at high concentration to a subject at the time of measurement. This transient gain in signal has enabled the non-invasive detection and imaging of gas ventilation and diffusion in the lungs, perfusion in blood vessels and tissues, and metabolic conversion in cells, animals, and patients. The rapid development of this method is based on advances in polarizer technology, the availability of suitable probe isotopes and molecules, improved MRI hardware and pulse sequence development. Acquisition strategies for hyperpolarized nuclei are not yet standardized and are set up individually at most sites depending on the specific requirements of the probe, the object of interest, and the MRI hardware. This review provides a detailed introduction to spatially resolved detection of hyperpolarized nuclei and summarizes novel and previously established acquisition strategies for different key areas of application.

Fast Imaging for Hyperpolarized MR Metabolic Imaging

MRI with hyperpolarized carbon-13 agents has created a new type of noninvasive, in vivo metabolic imaging that can be applied in cell, animal, and human studies. The use of ¹³ C-labeled agents, primarily [1-¹³ C]pyruvate, enables monitoring of key metabolic pathways with the ability to image substrate and products based on their chemical shift. Over 10 sites worldwide are now performing human studies with this new approach for studies of cancer, heart disease, liver disease, and kidney disease. Hyperpolarized metabolic imaging studies must be performed within several minutes following creation of the hyperpolarized agent due to irreversible decay of the net magnetization back to equilibrium, so fast imaging methods are critical. The imaging methods must include multiple metabolites, separated based on their chemical shift, which are also undergoing rapid metabolic conversion (via label exchange), further exacerbating the challenges of fast imaging. This review describes the state-of-the-art in fast imaging methods for hyperpolarized metabolic imaging. This includes the approach and tradeoffs between three major categories of fast imaging methods-fast spectroscopic imaging, model-based strategies, and metabolite specific imaging-as well additional options of parallel imaging, compressed sensing, tailored RF flip angles, refocused imaging methods, and calibration methods that can improve the scan coverage, speed, signal-to-noise ratio (SNR), resolution, and/or robustness of these studies. To date, these approaches have produced extremely promising initial human imaging results. Improvements to fast hyperpolarized metabolic imaging methods will provide better coverage, SNR, resolution, and reproducibility for future human imaging studies. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

Extended quantitative dynamic contrast-enhanced cardiac perfusion imaging in mice using accelerated data acquisition and spatially distributed, two-compartment exchange modeling

The objective of the present work was to improve data acquisition and quantification of dynamic contrast-enhanced perfusion imaging in the in vivo murine heart. Four-fold undersampled data were acquired in 14 mice and reconstructed using k-t SPARSE. A two-compartment exchange model was employed to provide additional characterization of myocardial tissue based on compartment volumes and the permeability surface area product. The feasibility of the proposed method was tested using compartment-based analysis of contrast-enhanced perfusion data acquired with intravascular and extracellular contrast agents. A significantly different permeability surface area product was measured for the intravascular versus extracellular contrast agent (0.13-0.15 ml/g/min vs 0.86-0.88 ml/g/min). The reduced extravasation also resulted in significantly smaller interstitial volumes of the intravascular versus extracellular agent (9.8-11% vs 45-47%). No difference was found for myocardial blood flow (6.5-7.2 ml/g/min vs 6.0-7.0 ml/g/min). The results presented here show that two-compartment exchange modeling in the in vivo murine heart is feasible and gives access to tissue parameters beyond myocardial blood flow.

A spin-thermodynamic approach to characterize spin dynamics in TEMPO-based samples for dissolution DNP at 7 T field

The spin dynamics of dissolution DNP samples consisting of 4.5 M [¹³C]urea in a mixture of (1/1)_Vol glycerol/water using 4-Oxo-TEMPO as a radical was investigated. We analyzed the DNP dynamics as function of radical concentration at 7 T and 3.4 T static magnetic field as well as function of deuteration of the solvent matrix at the high field. The spin dynamics could be reproduced in all cases, at least qualitatively, by a thermodynamic model based on spin temperatures of the nuclear Zeeman baths and an electron non-Zeeman (dipolar) bath. We find, however, that at high field (7 T) and low radical concentrations (25 mM) the nuclear spins do not reach the same spin temperature indicating a weak coupling of the two baths. At higher radical concentrations, as well as for all radical concentrations at low field (3.4 T), the two nuclear Zeeman baths reach the same spin temperature within experimental errors. Additionally, the spin system was prepared with different initial conditions. For these cases, the thermodynamic model was able to predict the time evolution of the system well. While the DNP profiles do not give clear indications to a specific polarization transfer mechanism, at high field (7 T) increased coupling is seen. The EPR line shapes cannot clarify this in absence of ELDOR type experiments, nevertheless DNP profiles and dynamics under frequency-modulated microwave irradiation illustrate the expected increase in coupling between electrons with increasing radical concentration.

Production of highly polarized [1-13 C]acetate by rapid decarboxylation of [2-13 C]pyruvate - application to hyperpolarized cardiac spectroscopy and imaging

PURPOSE

The objective of the present work was to develop and implement an efficient approach to hyperpolarize [1-¹³ C]acetate and apply it to in vivo cardiac spectroscopy and imaging.

METHODS

Rapid hydrogen peroxide induced decarboxylation was used to convert hyperpolarized [2-¹³ C]pyruvate into highly polarized [1-¹³ C]acetate employing an additional step following rapid dissolution of [2-¹³ C]pyruvate in a home-built multi-sample dissolution dynamic nuclear polarization system. Phantom dissolution experiments were conducted to determine optimal parameters of the decarboxylation reaction, retaining polarization and T₁ of [1-¹³ C]acetate. In vivo feasibility of detecting [1-¹³ C]acetate metabolism is demonstrated using slice-selective spectroscopy and multi-echo imaging of [1-¹³ C]acetate and [1-¹³ C]acetylcarnitine in the healthy rat heart.

RESULTS

The first in vivo signal was observed ~23 s after dissolution. At the corresponding time point in the phantom experiments, 97.9 ± 0.4% of [2-¹³ C]pyruvate were converted into [1-¹³ C]acetate by the decarboxylation reaction. T₁ and polarization of [1-¹³ C]acetate was determined to be 29.7 ± 1.9% and a 47.7 ± 0.5 s. Polarization levels of [2-¹³ C]pyruvate and [1-¹³ C]acetate were not significantly different after transfer to the scanner. In vivo, [1-¹³ C]acetate and [1-¹³ C]acetylcarnitine could be detected using spectroscopy and imaging.

CONCLUSION

Decarboxylation of hyperpolarized [2-¹³ C]pyruvate enables the efficient production of highly polarized [1-¹³ C]acetate that is applicable to study short-chain fatty acid metabolism in the in vivo heart.