Advances in human cardiac 31P-MR spectroscopy: SLOOP and clinical applications

Method for fast lipid reconstruction and removal processing in 1 H MRSI of the brain

PURPOSE

To develop a new rapid spatial filtering method for lipid removal, fast lipid reconstruction and removal processing (FLIP), which selectively isolates and removes interfering lipid signals from outside the brain in a full-FOV 2D MRSI and whole-brain 3D echo planar spectroscopic imaging (EPSI).

THEORY AND METHODS

FLIP uses regularized least-squares regression based on spatial prior information from MRI to selectively remove lipid signals originating from the scalp and measure the brain metabolite signals with minimum cross contamination. FLIP is a noniterative approach, thus allowing a rapid processing speed, and uses only spatial information without any spectral priors. The performance of FLIP was compared with the Papoulis-Gerchberg algorithm (PGA), Hankel singular value decomposition (HSVD), and fast image reconstruction with L2 regularization (L2).

RESULTS

FLIP in both 2D and 3D MRSI resulted in consistent metabolite quantification in a greater number of voxels with less concentration variation than other algorithms, demonstrating effective and robust lipid-removal performance. The percentage of voxels that met quality criteria with FLIP, PGA, HSVD, and L2 processing were 90%, 57%, 29%, and 42% in 2D MRSI, and 80%, 75%, 76%, and 74% in 3D EPSI, respectively. The quantification results of full-FOV MRSI using FLIP were comparable to those of volume-localized MRSI, while allowing significantly increased spatial coverage. FLIP performed the fastest in 2D MRSI.

CONCLUSION

FLIP is a new lipid-removal algorithm that promises fast and effective lipid removal with improved volume coverage in MRSI.

Deuterium Metabolic Imaging-Rediscovery of a Spectroscopic Tool

The growing demand for metabolism-specific imaging techniques has rekindled interest in Deuterium (2H) Metabolic Imaging (DMI), a robust method based on administration of a substrate (glucose, acetate, fumarate, etc.) labeled with the stable isotope of hydrogen and the observation of its metabolic fate in three-dimensions. This technique allows the investigation of multiple metabolic processes in both healthy and diseased states. Despite its low natural abundance, the short relaxation time of deuterium allows for rapid radiofrequency (RF) pulses without saturation and efficient image acquisition. In this review, we provide a comprehensive picture of the evolution of DMI over the course of recent decades, with a special focus on its potential clinical applications.

Fast, regional three-dimensional hybrid (1D-Hadamard 2D-rosette) proton MR spectroscopic imaging in the human temporal lobes

1 H-MRSI is commonly performed with gradient phase encoding, due to its simplicity and minimal radio frequency (RF) heating (specific absorption rate). Its two well-known main problems-(i) "voxel bleed" due to the intrinsic point-spread function, and (ii) chemical shift displacement error (CSDE) when slice-selective RF pulses are used, which worsens with increasing volume of interest (VOI) size-have long become accepted as unavoidable. Both problems can be mitigated with Hadamard multislice RF encoding. This is demonstrated and quantified with numerical simulations, in a multislice phantom and in five healthy young adult volunteers at 3 T, targeting a 2-cm thick temporal lobe VOI through the bilateral hippocampus. This frequently targeted region (e.g. in epilepsy and Alzheimer's disease) is subject to strong, 1-2 ppm.cm-1 regional B0, susceptibility gradients that can dramatically reduce the signal-to-noise ratio (SNR) and water suppression effectiveness. The chemical shift imaging (CSI) sequence used a 3-ms Shinnar-Le Roux (SLR) 90° RF pulse, acquiring eight steps in the slice direction. The Hadamard sequence acquired two overlapping slices using the same SLR 90° pulses, under twofold stronger gradients that proportionally halved the CSDE. Both sequences used 2D 20 × 20 rosette spectroscopic imaging (RSI) for in-plane spatial localization and both used RF and gradient performance characteristics that are easily met by all modern MRI instruments. The results show that Hadamard spectroscopic imaging (HSI) suffered dramatically less signal bleed within the VOI compared with CSI (<1% vs. approximately 26% in simulations; and 5%-8% vs. >50%) in a phantom specifically designed to test these effects. The voxels' SNR per unit volume per unit time was also 40% higher for HSI. In a group of five healthy volunteers, we show that HSI with in-plane 2D-RSI facilitates fast, 3D multivoxel encoding at submilliliter spatial resolution, over the bilateral human hippocampus, in under 10 min, with negligible CSDE, spectral and spatial contamination and more than 6% improved SNR per unit time per unit volume.

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.

Investigating the effect of trigger delay on cardiac 31P MRS signals

The heart’s geometry and its metabolic activity vary over the cardiac cycle. The effect of these fluctuations on phosphorus (³¹P) magnetic resonance spectroscopy (MRS) data quality and metabolite ratios was investigated. 12 healthy volunteers were measured using a 7 T MR scanner and a cardiac ³¹P-¹H loop coil. ³¹P chemical shift imaging data were acquired untriggered and at four different times during the cardiac cycle using acoustic triggering. Signals of adenosine-triphosphate (ATP), phosphocreatine (PCr), inorganic phosphate (Pi) and 2,3-diphosphoglycerate (2,3-DPG) and their fit quality as Cramér-Rao lower bounds (CRLB) were quantified including corrections for contamination by ³¹P signals from blood, flip angle, saturation and total acquisition time. The myocardial filling factor was estimated from cine short axis views. The corrected signals of PCr and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ-ATP were higher during end-systole and lower during diastasis than in untriggered acquisitions (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.05$$\end{document}P<0.05). Signal intensities of untriggered scans were between those with triggering to end-systole and diastasis. Fit quality of PCr and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ-ATP peaks was best during end-systole when blood contamination of ATP and Pi signals was lowest. While metabolite ratios and pH remained stable over the cardiac cycle, signal amplitudes correlated strongly with myocardial voxel filling. Triggering of cardiac ³¹P MRS acquisitions improves signal amplitudes and fit quality if the trigger delay is set to end-systole. We conclude that triggering to end-systole is superior to triggering to diastasis.

Use of cardiac magnetic resonance to detect changes in metabolism in heart failure

The heart has a massive adenosine triphosphate (ATP) requirement, produced from the oxidation of metabolic substrates such as fat and glucose. Magnetic resonance spectroscopy offers a unique opportunity to probe this biochemistry: ³¹Phosphorus spectroscopy can demonstrate the production of ATP and quantify levels of the transport molecule phosphocreatine while 13Carbon spectroscopy can demonstrate the metabolic fates of glucose in real time. These techniques allow the metabolic deficits in heart failure to be interrogated and can be a potential future clinical tool.

Imaging based magnetic resonance spectroscopy (MRS) localization for quantitative neurochemical analysis and cerebral metabolism studies

Accurate quantitative metabolic imaging of the brain presents significant challenges due to the complexity and heterogeneity of its structures and compositions with distinct compartmentations of brain tissue types (e.g., gray and white matter). The brain is compartmentalized into various regions based on their unique functions and locations. In vivo magnetic resonance spectroscopy (MRS) techniques allow non-invasive measurements of neurochemicals in either single voxel or multiple voxels, yet the spatial resolution and detection sensitivity of MRS are significantly lower compared with MRI. A fundamentally different approach, namely spectral localization by imaging (SLIM) provides a new framework that overcomes major limitations of conventional MRS techniques. Conventional MRS allows only rectangular voxel shapes that do not conform to the shapes of brain structures or lesions, while SLIM allows compartments with arbitrary shapes. However, the restrictive assumption proposed in the original concept of SLIM, i.e., compartmental homogeneity, led to spectral localization errors, which have limited its broad applications. This review focuses on the recent technical frontiers of image-based MRS localization techniques that overcome the limitations of SLIM through the development and implementation of various new strategies, including incorporation of magnetic field inhomogeneity corrections, the use of multiple receiver coils, and prospective optimization of data acquisition.

B0-adjusted and sensitivity-encoded spectral localization by imaging (BASE-SLIM) in the human brain in vivo

Spectral localization by imaging (SLIM) based magnetic resonance spectroscopy (MRS) provides a framework that overcomes major limitations of current MRS techniques, which allow only rectangular voxel shapes that do not conform to the shapes of brain structures or lesions. However, the restrictive assumption of compartmental homogeneity in SLIM can lead to localization errors, thus its applications have been very limited to date. SLIM-based localization is subject to errors due to inhomogeneous B0 and B1 fields, particularly in organs with complex compartmental geometry including the human brain. The limitations of SLIM were overcome through the development and implementation of B0-adjusted and sensitivity-encoded SLIM (BASE-SLIM) that includes corrections for inhomogeneities of both B0 and B1 fields throughout the volume of interest. In this study, we demonstrate significantly improved localization accuracy in compartments with arbitrary shapes and reliable quantification of metabolite concentrations in gray and white matter of the human brain using the BASE-SLIM technique.

Highly-accelerated quantitative 2D and 3D localized spectroscopy with linear algebraic modeling (SLAM) and sensitivity encoding

Noninvasive magnetic resonance spectroscopy (MRS) with chemical shift imaging (CSI) provides valuable metabolic information for research and clinical studies, but is often limited by long scan times. Recently, spectroscopy with linear algebraic modeling (SLAM) was shown to provide compartment-averaged spectra resolved in one spatial dimension with many-fold reductions in scan-time. This was achieved using a small subset of the CSI phase-encoding steps from central image k-space that maximized the signal-to-noise ratio. Here, SLAM is extended to two- and three-dimensions (2D, 3D). In addition, SLAM is combined with sensitivity-encoded (SENSE) parallel imaging techniques, enabling the replacement of even more CSI phase-encoding steps to further accelerate scan-speed. A modified SLAM reconstruction algorithm is introduced that significantly reduces the effects of signal nonuniformity within compartments. Finally, main-field inhomogeneity corrections are provided, analogous to CSI. These methods are all tested on brain proton MRS data from a total of 24 patients with brain tumors, and in a human cardiac phosphorus 3D SLAM study at 3T. Acceleration factors of up to 120-fold versus CSI are demonstrated, including speed-up factors of 5-fold relative to already-accelerated SENSE CSI. Brain metabolites are quantified in SLAM and SENSE SLAM spectra and found to be indistinguishable from CSI measures from the same compartments. The modified reconstruction algorithm demonstrated immunity to maladjusted segmentation and errors from signal heterogeneity in brain data. In conclusion, SLAM demonstrates the potential to supplant CSI in studies requiring compartment-average spectra or large volume coverage, by dramatically reducing scan-time while providing essentially the same quantitative results.

Magnetic resonance Spectroscopy with Linear Algebraic Modeling (SLAM) for higher speed and sensitivity

Speed and signal-to-noise ratio (SNR) are critical for localized magnetic resonance spectroscopy (MRS) of low-concentration metabolites. Matching voxels to anatomical compartments a priori yields better SNR than the spectra created by summing signals from constituent chemical-shift-imaging (CSI) voxels post-acquisition. Here, a new method of localized Spectroscopy using Linear Algebraic Modeling (SLAM) is presented, that can realize this additional SNR gain. Unlike prior methods, SLAM generates spectra from C signal-generating anatomic compartments utilizing a CSI sequence wherein essentially only the C central k-space phase-encoding gradient steps with highest SNR are retained. After MRI-based compartment segmentation, the spectra are reconstructed by solving a sub-set of linear simultaneous equations from the standard CSI algorithm. SLAM is demonstrated with one-dimensional CSI surface coil phosphorus MRS in phantoms, the human leg and the heart on a 3T clinical scanner. Its SNR performance, accuracy, sensitivity to registration errors and inhomogeneity, are evaluated. Compared to one-dimensional CSI, SLAM yielded quantitatively the same results 4-times faster in 24 cardiac patients and healthy subjects. SLAM is further extended with fractional phase-encoding gradients that optimize SNR and/or minimize both inter- and intra-compartmental contamination. In proactive cardiac phosphorus MRS of six healthy subjects, both SLAM and fractional-SLAM (fSLAM) produced results indistinguishable from CSI while preserving SNR gains of 36-45% in the same scan-time. Both SLAM and fSLAM are simple to implement and reduce the minimum scan-time for CSI, which otherwise limits the translation of higher SNR achievable at higher field strengths to faster scanning.

Spatial localization in nuclear magnetic resonance spectroscopy

The ability to select a discrete region within the body for signal acquisition is a fundamental requirement of in vivo NMR spectroscopy. Ideally, it should be possible to tailor the selected volume to coincide exactly with the lesion or tissue of interest, without loss of signal from within this volume or contamination with extraneous signals. Many techniques have been developed over the past 25 years employing a combination of RF coil properties, static magnetic field gradients and pulse sequence design in an attempt to meet these goals. This review presents a comprehensive survey of these techniques, their various advantages and disadvantages, and implications for clinical applications. Particular emphasis is placed on the reliability of the techniques in terms of signal loss, contamination and the effect of nuclear relaxation and J-coupling. The survey includes techniques based on RF coil and pulse design alone, those using static magnetic field gradients, and magnetic resonance spectroscopic imaging. Although there is an emphasis on techniques currently in widespread use (PRESS, STEAM, ISIS and MRSI), the review also includes earlier techniques, in order to provide historical context, and techniques that are promising for future use in clinical and biomedical applications.

The relationship between plasma BNP level and the myocardial phosphocreatine/adenosine triphosphate ratio determined by phosphorus-31 magnetic resonance spectroscopy in patients with dilated cardiomyopathy

The purpose of this study was to evaluate the correlation between the plasma B-type natriuretic peptide (BNP) level and the myocardial phosphocreatine/adenosine triphosphate ratio determined using rapid phosphorus-31 magnetic resonance spectroscopy (31P-MRS) in patients with dilated cardiomyopathy (DCM). Thirteen DCM patients, who had slight or moderate heart failure, were examined. The plasma BNP was measured on a day close to the rapid 31P-MRS study. 31P-MRS measurements were conducted with a 1.5-T MR instrument. The plasma BNP levels tended to be correlated negatively with the myocardial phosphocreatine/adenosine triphosphate, although the correlation did not reach statistical significance (r = -0.54, p = 0.06). By contrast, the log of the plasma BNP levels was correlated negatively with the myocardial phosphocreatine/adenosine triphosphate (r = -0.73, p < 0.01). Our results indicate that the myocardial energy metabolism evaluated using 31P-MRS tends to be correlated with the severity of heart failure and left ventricular dysfunction estimated using the plasma BNP levels in DCM patients. This paper provides additional information regarding the relationship between the BNP and myocardial energy metabolism in DCM patients.

Feasibility of rapid-sequence 31P magnetic resonance spectroscopy in cardiac patients

PURPOSE

To determine the clinical feasibility of rapid-sequence phosphorus-31 magnetic resonance spectroscopy (31P-MRS) of the heart with cardiac patients using a 1.5T clinical MR system.

MATERIAL AND METHODS

Twenty cardiac patients, i.e. dilated cardiomyopathy (DCM) 13 cases, hypertrophic cardiomyopathy (HCM) 3 cases, hypertensive heart diseases (HHD) 3 cases, and aortic regurgitation (AR) 1 case were examined using rapid cardiac 31P-MRS. Complete three-dimensional localization was performed using a two-dimensional phosphorus chemical-shift imaging sequence in combination with 30-mm axial slice-selective excitation. The rapid-sequence 31P-MRS procedure was phase encoded in arrays of 8 x 8 steps with an average of 4 acquisitions. The total examination time, including proton imaging and shimming, for the rapid cardiac 31P-MRS procedure, ranged from 10 to 15 min, depending on the heart rate. Student's t test was used to compare creatine phosphate (PCr)/adenosine triphosphate (ATP) ratios from the cardiac patients with those of the control subjects (n = 13).

RESULTS

The myocardial PCr/ATP ratio obtained by rapid 31P-MRS was significantly lower (P < 0.001) in DCM patients (1.82 +/- 0.33, mean +/- SD), and in patients with global myocardial dysfunction (combined data for 20 patients: 1.89 +/- 0.32) than in normal volunteers (2.96 +/- 0.59). These results are similar to previous studies.

CONCLUSION

Rapid-sequence 31P-MRS may be a valid diagnostic tool for patients with cardiac disease.

Detection of myocardial infarctions by acquisition-weighted 31P-MR spectroscopy in humans

PURPOSE

To determine whether the recently applied technique of acquisition-weighted 31P-MR spectroscopy (AW-MRS) allows for the detection of depressed energy metabolism in patients with inferior wall myocardial infarctions.

MATERIALS AND METHODS

Eight patients with subacute myocardial infarction and wall motion abnormalities restricted to the inferior wall were examined with a 1.5-T MR scanner. Global and regional left ventricular (LV) function was assessed by cine MRI, and the size and extent of myocardial infarction was assessed by late enhancement (LE). MRS was performed with an AW three-dimensional chemical shift imaging sequence. Phosphocreatine/ATP ratios were determined with the postprocessing model AMARES for four voxels positioned in the anterior, lateral, inferior, and septal parts of the LV.

RESULTS

The LV ejection fraction (EF) was reduced to 37.5%+/-9.0%. Seven of eight patients had transmural LE in the inferior wall, and one patient showed subendocardial enhancement in the inferior-lateral parts. Phosphocreatine/ATP ratios of the inferior wall were significantly reduced (P <0.05) compared to all other parts of the LV (1.03 +/- 0.39 (inferior), 1.67 +/- 0.81 (lateral), 1.73 +/- 0.29 (anterior), and 1.49 +/- 0.31 (septal)). The ratios in five of seven patients with transmural enhancement were <1.00 in the inferior wall.

CONCLUSION

Acquisition weighting allows for the detection of inferior wall infarctions in patients. Transmural signal enhancement is associated with significant depression of phosphocreatine/ATP ratios.

Assessing gene and cell therapies applied in striated skeletal and cardiac muscle: is there a role for nuclear magnetic resonance?

Gene and cell therapies convey high hopes for treatment of skeletal and heart muscle diseases. In the experimental protocols under development as well as in the first clinical trials, longitudinal control by an atraumatic procedure is needed. Nuclear magnetic resonance (NMR), via its two modalities, imaging or spectroscopy, should play a major role both for in vivo animal and human studies, because of the great number of parameters that can be measured, sequentially or simultaneously, and because of its aptitude to monitor several steps of protocols, in particular to detect physiological modifications induced by therapies. We review here the many possible applications of nuclear magnetic resonance in gene/cell therapies where muscle is the target organ, with emphasis on the application of nuclear magnetic resonance to functional studies.

Volume tracking cardiac 31P spectroscopy

The limited reliability and accuracy of cardiac spectroscopy have been partly attributed to effects from respiratory motion. In this work, we developed a prospective volume tracking method for respiratory motion compensation based on multiple navigator echoes and demonstrated its application in cardiac (31)P spectroscopy. The sequence consists of two 2D selective excitation pulses preceding the spectroscopic experiment to sample respiratory motion components. The navigator information is evaluated in real-time to calculate the shift of the heart from respiration. Based on the displacement information, the spectroscopic volume and/or grid position is prospectively corrected to track the volume of interest. The method was validated with a moving compartment phantom simulating in vivo respiratory motion. With volume tracking, no signal contamination was apparent. Spectra obtained in 14 healthy volunteers were evaluated using time-domain fitting procedures. The fitting accuracy improved consistently with volume tracking compared to data from non-navigated reference acquisitions. Compared to other gating approaches available for spectroscopy, the current technique does not degrade the scan efficiency, thus allowing effective use of scan time.