Numerical simulation of PRESS localized MR spectroscopy

Comparison of baseline correction algorithms for in vivo 1H-MRS

Proton MRS is used clinically to collect localized, quantitative metabolic data from living tissues. However, the presence of baselines in the spectra complicates accurate MRS data quantification. The occurrence of baselines is not specific to short-echo-time MRS data. In short-echo-time MRS, the baseline consists typically of a dominating macromolecular (MM) part, and can, depending on B0 shimming, poor voxel placement, and/or localization sequences, also contain broad water and lipid resonance components, indicated by broad components (BCs). In long-echo-time MRS, the MM part is usually much smaller, but BCs may still be present. The sum of MM and BCs is denoted by the baseline. Many algorithms have been proposed over the years to tackle these artefacts. A first approach is to identify the baseline itself in a preprocessing step, and a second approach is to model the baseline in the quantification of the MRS data themselves. This paper gives an overview of baseline handling algorithms and also proposes a new algorithm for baseline correction. A subset of suitable baseline removal algorithms were tested on in vivo MRSI data (semi-LASER at TE = 40 ms) and compared with the new algorithm. The baselines in all datasets were removed using the different methods and subsequently fitted using spectrIm-QMRS with a TDFDFit fitting model that contained only a metabolite basis set and lacked a baseline model. The same spectra were also fitted using a spectrIm-QMRS model that explicitly models the metabolites and the baseline of the spectrum. The quantification results of the latter quantification were regarded as ground truth. The fit quality number (FQN) was used to assess baseline removal effectiveness, and correlations between metabolite peak areas and ground truth models were also examined. The results show a competitive performance of our new proposed algorithm, underscoring its automatic approach and efficiency. Nevertheless, none of the tested baseline correction methods achieved FQNs as good as the ground truth model. All separately applied baseline correction methods introduce a bias in the observed metabolite peak areas. We conclude that all baseline correction methods tested, when applied as a separate preprocessing step, yield poorer FQNs and biased quantification results. While they may enhance visual display, they are not advisable for use before spectral fitting.

MEGA-PRESS of GABA+: Influences of acquisition parameters

γ-aminobutyric acid (GABA) was the first molecule that was edited with MEGA-PRESS. GABA edited spectroscopy is challenged by limited selectivity of editing pulses. Coediting of resonances from macromolecules (MM) is the greatest single limitation of GABA edited spectroscopy. In this contribution, relative signal contributions from GABA, MM and homocarnosine to the total MEGA-PRESS edited signal at ~3 ppm, i.e., GABA+, are simulated at 3 tesla using several acquisition schemes. The base scheme is modeled after those currently supplied by vendors: it uses typical pulse shapes and lengths, it minimizes the first echo time (TE), and the delay between the editing pulses is kept at TE/2. Edited spectra are simulated for imperfect acquisition parameters such as incorrect frequency, larger chemical shift displacement, incorrect transmit B₁ -field calibration for localization and editing pulses, and longer TE. An alternative timing scheme and longer editing pulses are also considered. Additional simulations are performed for symmetric editing around the MM frequency to suppress the MM signal. The relative influences of these acquisition parameters on the constituents of GABA+ are examined from the perspective of modern experimental designs for investigating brain GABA concentration differences in healthy and diseased humans. Other factors that influence signal contributions, such as T₁ and T₂ relaxation times are also considered.

Parameterization of metabolite and macromolecule contributions in interrelated MR spectra of human brain using multidimensional modeling

Macromolecular signals are crucial constituents of short echo-time ¹ H MR spectra with potential clinical implications in themselves as well as essential ramifications for the quantification of the usually targeted metabolites. Their parameterization, needed for general fitting models, is difficult because of their unknown composition. Here, a macromolecular signal parameterization together with metabolite signal quantification including relaxation properties is investigated by multidimensional modeling of interrelated 2DJ inversion-recovery (2DJ-IR) datasets. Simultaneous and iterative procedures for defining the macromolecular background (MMBG) as mono-exponentially or generally decaying signals over TE are evaluated. Varying prior knowledge and restrictions in the metabolite evaluation are tested to examine their impact on results and fitting stability for two sets of three-dimensional spectra acquired with metabolite-cycled PRESS from cerebral gray and white matter locations. One dataset was used for model optimization, and also examining the influence of prior knowledge on estimated parameters. The most promising model was applied to a second dataset. It turned out that the mono-exponential decay model appears to be inadequate to represent TE-dependent signal features of the MMBG. TE-adapted MMBG spectra were therefore determined. For a reliable overall quantification of implicated metabolite concentrations and relaxation times, a general fitting model had to be constrained in terms of the number of fitting variables and the allowed parameter space. With such a model in place, fitting precision for metabolite contents and relaxation times was excellent, while fitting accuracy is difficult to judge and bias was likely influenced by the type of fitting constraints enforced. In summary, the parameterization of metabolite and macromolecule contributions in interrelated MR spectra has been examined by using multidimensional modeling on complex 2DJ-IR datasets. A tightly restricted model allows fitting of individual subject data with high fitting precision documented in small Cramér-Rao lower bounds, good repeatability values and a relatively small spread of estimated concentration and relaxation values for a healthy subject cohort.

Effects of noise and linewidth on in vivo analysis of glutamate at 3 T

Magnetic resonance spectroscopy (MRS) can noninvasively detect metabolites in vivo, including glutamate (Glu). However, quantification is known to be affected by the overlaps among metabolite resonance lines and background macromolecule signals. We found that adding a moderate amount of noise or line broadening (2 Hz) caused large variations in concentration of Glu and other metabolites, when determined by LCModel analysis of in vivo short-echo time (TE) spectra. Theses variations were largely attributed to strong spectral baselines in short TE spectra, especially near 2.35 ppm, as well as overlapping metabolite resonance lines. To address this issue, we acquired in vivo data at 3 T using both short-TE and the multiple echo time J-resolved point-resolved spectroscopy (JPRESS) MRS techniques. We found that one-dimensional (1D) JPRESS, by simultaneously fitting the two cross-sections of JPRESS at J = 0 and J = 7.5 Hz, was highly resistant to variations in noise levels and spectral linewidths. Our results demonstrate that LCModel analysis of short-TE data is highly sensitive to variations in noise levels and spectral linewidths and this sensitivity is greatly reduced by 1D JPRESS given its substantially reduced baselines and enhanced spectral resolution.

Highly specific determination of IDH status using edited in vivo magnetic resonance spectroscopy

Background

Mutations in the isocitrate dehydrogenase (IDH) enzyme affect 40% of gliomas and represent a major diagnostic and prognostic marker. The goals of this study were to evaluate the performance of noninvasive magnetic resonance spectroscopy (MRS) methods to determine the IDH status of patients with brain gliomas through detection of the oncometabolite 2-hydroxyglutarate (2HG) and to compare performance of these methods with DNA sequencing and tissue 2HG analysis.

Methods

Twenty-four subjects with suspected diagnosis of low-grade glioma were included prospectively in the study. For all subjects, MRS data were acquired at 3T using 2 MRS methods, edited MRS using Mescher-Garwood point-resolved spectroscopy (MEGA-PRESS) sequence and a PRESS sequence optimized for 2HG detection, using a volume of interest larger than 6 mL. IDH mutational status was determined by a combination of automated immunohistochemical analysis and Sanger sequencing. Levels of 2HG in tissue samples measured by gas chromatography-mass spectrometry were compared with those estimated by MRS.

Results

Edited MRS provided 100% specificity and 100% sensitivity in the detection of 2HG. The 2HG levels estimated by this technique were in line with those derived from tissue samples. Optimized PRESS provided lower performance, in agreement with previous findings.

Conclusions

Our results suggest that edited MRS is one of the most reliable tools to predict IDH mutation noninvasively, showing high sensitivity and specificity for 2HG detection. Integrating edited MRS in clinical practice may be highly beneficial for noninvasive diagnosis of glioma, prognostic assessment, and treatment planning.

Apparent diffusion coefficients of the five major metabolites measured in the human brain in vivo at 3T

PURPOSE

To measure the apparent diffusion coefficients (ADC) of five main metabolites in the human brain at 3T with PRESS and STEAM, avoiding measurement biases because of cross-terms. Cross-terms arise from interactions between slice-selection and spoiler gradients in the localized spectroscopy sequence and the diffusion gradients.

METHODS

Diffusion-weighted spectra were acquired from the prefrontal cortex in five healthy subjects using STEAM (echo time [TE]/mixing time [TM]/pulse repetition time [TR] = 21.22/105/3000 ms, b-values = 0 and 3172 s/mm2 ) and PRESS (TE/TR = 54.2/3000 ms, b-values = 0 and 2204 s/mm2 ). Diffusion weighting was applied using bipolar gradients in three orthogonal directions. Post-processed spectra were analyzed with LCModel, and the trace/3 ADC values were calculated.

RESULTS

Comparable trace/3 ADC values (0.14-0.18 µm2 /ms) were obtained for five main metabolites with both methods. These metabolites were quantified with Cramér-Rao lower bounds below 15%.

CONCLUSION

The ADC values of the five main metabolites were successfully measured in the human brain at 3T with eliminated directional dependence. Both STEAM and PRESS can be used to probe the diffusivity of metabolites in normal brain and various pathologies on the clinical scanner with slightly higher precision achieved with STEAM for glutamate and myo-inositol. Magn Reson Med 79:2896-2901, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Improved localization, spectral quality, and repeatability with advanced MRS methodology in the clinical setting

PURPOSE

To investigate the utility of an advanced magnetic resonance spectroscopy (MRS) protocol in the clinical setting, and to compare the localization accuracy, spectral quality, and quantification repeatability between this advanced and the conventional vendor-provided MRS protocol on a clinical 3T platform.

METHODS

Proton spectra were measured from the posterior cingulate cortices in 30 healthy elderly subjects by clinical MR technologists using a vendor-provided (point resolved spectroscopy with advanced 3D gradient-echo B0 shimming) and an advanced (semi-LASER with FAST(EST)MAP shimming) protocol, in random order. Spectra were quantified with LCModel using standard pipelines for the clinical and research settings, respectively.

RESULTS

The advanced protocol outperformed the vendor-provided protocol in localization accuracy (chemical-shift-displacement error: 2.0%/ppm, semi-LASER versus 11.6%/ppm, point resolved spectroscopy), spectral quality (water linewidth: 6.1 ± 1.8 Hz, FAST(EST)MAP versus 10.5 ± 3.7 Hz, 3D gradient echo; P < 7e-6; residual water: 0.08 ± 0.12%, VAPOR versus 0.45 ± 0.50%, WET; P < 2e-5) and within-session repeatability of metabolite concentrations, particularly of low signal-to-noise ratio data with two to eight averages (test-retest coefficients of variance of metabolite concentrations, P < 0.01). Concentrations of J-coupled metabolites such as γ-aminobutyric acid and glutamate were biased when using the default pipeline with simulated macromolecules.

CONCLUSIONS

The quality of MRS data can be improved using advanced acquisition and analysis protocols on standard 3T hardware in the clinical setting, which can facilitate robust applications in central nervous system diseases. Magn Reson Med 79:1241-1250, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Simultaneous editing of GABA and glutathione at 7T using semi-LASER localization

PURPOSE

To demonstrate simultaneous editing of the two most commonly edited and overlapping signals, γ-aminobutyric acid (GABA), and glutathione (GSH), with Hadamard encoding and reconstruction of MEGA-edited spectroscopy (HERMES) using sLASER localization at 7T.

METHODS

Density matrix simulations of HERMES at 7T were carried out and compared with phantom experiments. Additional phantom experiments were performed to characterize the echo time (TE) -dependent modulation of GABA- and GSH-edited HERMES spectra at TE of 80-160 ms. In vivo experiments were performed in 10 healthy volunteers, comparing HERMES (11 min) to sequentially acquired MEGA-sLASER detection of GABA and GSH (2 × 11 min).

RESULTS

Simulations of HERMES show GABA- and GSH-edited spectra with negligible levels of crosstalk, and give modest agreement with phantom spectra. The TE series of GABA- and GSH-edited HERMES spectra modulate as a result of T₂ relaxation and coupling evolution, with GABA showing a stronger TE-dependence. In vivo HERMES experiments show well-edited GABA and GSH signals. Measured concentrations are not statistically different between HERMES and MEGA-sLASER for GABA (1. 051 ± 0.254 i.u. and 1.053 ± 0.248 i.u; P > 0.985) or GSH (0.300 ± 0.091 i.u. and 0.302 ± 0.093 i.u; P > 0.940).

CONCLUSION

Simulated, phantom and in vivo measurements of HERMES show excellent segregation of GABA- and GSH-edited signals, and excellent agreement with separately acquired MEGA-sLASER data. HERMES allows two-fold acceleration of editing while maintaining spectral quality compared with sequentially acquired MEGA-sLASER measurements. Magn Reson Med 80:474-479, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Transverse relaxation time constants of the five major metabolites in human brain measured in vivo using LASER and PRESS at 3 T

PURPOSE

The goal of this study was to measure and compare the apparent transverse relaxation time constants (T₂ ) of five intracellular metabolites using localization by adiabatic selective refocusing (LASER) and point-resolved spectroscopy (PRESS) sequences in the human brain at 3 T.

METHODS

Five healthy subjects were studied at 3 T. ¹ H spectra from the prefrontal cortex were acquired at six different echo times using LASER and PRESS sequences. Postprocessed data were analyzed with LCModel, and the resulting amplitudes were fitted using a mono-exponential decay function to determine the T₂ of metabolites.

RESULTS

Twenty-one percent higher apparent T₂ values for the singlet resonances of N-acetyl aspartate, total creatine, and total choline were measured with LASER as compared with PRESS, whereas comparable apparent T₂ values were measured for strongly coupled metabolites, glutamate, and myo-inositol, with both sequences.

CONCLUSIONS

Reliable T₂ measurements were obtained with both sequences for the five major intracellular metabolites. The LASER sequence appears to be more efficient in suppressing the diffusion component for singlets (having nonexchangeable protons) compared to J-coupled metabolites. Magn Reson Med 79:1260-1265, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Fast computation of full density matrix of multispin systems for spatially localized in vivo magnetic resonance spectroscopy

PURPOSE

Numerical simulations of three-dimensionally localized MRS spectra have been very time consuming for multispin systems because the current state-of-the-art method requires computation of a large ensemble of spins pixel-by-pixel in three dimensional space. This paper describes a highly accelerated technique for computing spatially localized MRS spectra using the full solution to the Liouville-von Neumann equation.

METHODS

The time evolution of spatially localized multispin density matrix as the full solution to the Liouville-von Neumann equation was analyzed. A new technique based on one dimensional spatial projection of the full density matrix was proposed. This method was implemented using a computer program written in Java language.

RESULTS

The MRS spectra calculated using the new method were found to be identical to conventional three-dimensional simulation for the same digitization of the voxel while the new method reduced computation time by orders of magnitude and led to not only improved speed but also accuracy. Applications of the new method to phantom studies of multispin systems and quantification of in vivo MRS spectra of brain were demonstrated.

CONCLUSION

The dramatically enhanced computational efficiency makes accurate simulation of localized MRS spectra highly accessible for calculating basis sets for spectral quantification and for optimizing pulse sequences.

Improved localisation for 2-hydroxyglutarate detection at 3T using long-TE semi-LASER

2-hydroxyglutarate (2-HG) has emerged as a biomarker of tumour cell IDH mutations that may enable the differential diagnosis of glioma patients. At 3 Tesla, detection of 2-HG with magnetic resonance spectroscopy is challenging because of metabolite signal overlap and a spectral pattern modulated by slice selection and chemical shift displacement. Using density matrix simulations and phantom experiments, an optimised semi-LASER scheme (TE = 110 ms) improves localisation of the 2-HG spin system considerably compared to an existing PRESS sequence. This results in a visible 2-HG peak in the in vivo spectra at 1.9 ppm in the majority of IDH mutated tumours. Detected concentrations of 2-HG were similar using both sequences, although the use of semi-LASER generated narrower confidence intervals. Signal overlap with glutamate and glutamine, as measured by pairwise fitting correlation was reduced. Lactate was readily detectable across glioma patients using the method presented here (mean CLRB: (10±2)%). Together with more robust 2-HG detection, long TE semi-LASER offers the potential to investigate tumour metabolism and stratify patients in vivo at 3T.

Localized 2D COSY sequences: Method and experimental evaluation for a whole metabolite quantification approach

Two-dimensional spectroscopy offers the possibility to unambiguously distinguish metabolites by spreading out the multiplet structure of J-coupled spin systems into a second dimension. Quantification methods that perform parametric fitting of the 2D MRS signal have recently been proposed for resolved PRESS (JPRESS) but not explicitly for Localized Correlation Spectroscopy (LCOSY). Here, through a whole metabolite quantification approach, correlation spectroscopy quantification performances are studied. The ability to quantify metabolite relaxation constant times is studied for three localized 2D MRS sequences (LCOSY, LCTCOSY and the JPRESS) in vitro on preclinical MR systems. The issues encountered during implementation and quantification strategies are discussed with the help of the Fisher matrix formalism. The described parameterized models enable the computation of the lower bound for error variance--generally known as the Cramér Rao bounds (CRBs), a standard of precision--on the parameters estimated from these 2D MRS signal fittings. LCOSY has a theoretical net signal loss of two per unit of acquisition time compared to JPRESS. A rapid analysis could point that the relative CRBs of LCOSY compared to JPRESS (expressed as a percentage of the concentration values) should be doubled but we show that this is not necessarily true. Finally, the LCOSY quantification procedure has been applied on data acquired in vivo on a mouse brain.

GABA concentration in posterior cingulate cortex predicts putamen response during resting state fMRI

The role of neurotransmitters in the activity of resting state networks has been gaining attention and has become a field of research with magnetic resonance spectroscopy (MRS) being one of the key techniques. MRS permits the measurement of γ-aminobutyric acid (GABA) and glutamate levels, the central biochemical constituents of the excitation-inhibition balance in vivo. The inhibitory effects of GABA in the brain have been largely investigated in relation to the activity of resting state networks in functional magnetic resonance imaging (fMRI). In this study GABA concentration in the posterior cingulate cortex (PCC) was measured using single voxel spectra acquired with standard point resolved spectroscopy (PRESS) from 20 healthy male volunteers at 3 T. Resting state fMRI was consecutively measured and the values of GABA/Creatine+Phosphocreatine ratio (GABA ratio) were included in a general linear model matrix as a step of dual regression analysis in order to identify voxels whose neuroimaging metrics during rest were related to individual levels of the GABA ratio. Our data show that the connection strength of putamen to the default-mode network during resting state has a negative linear relationship with the GABA ratio measured in the PCC. These findings highlight the role of PCC and GABA in segregation of the motor input, which is an inherent condition that characterises resting state.

Quantification of glutamate and glutamine using constant-time point-resolved spectroscopy at 3 T

Separate quantification of glutamate (Glu) and glutamine (Gln) using conventional MRS on clinical scanners is challenging. In previous work, constant-time point-resolved spectroscopy (CT-PRESS) was optimized at 3 T to detect Glu, but did not resolve Gln. To quantify Glu and Gln, a time-domain basis set was constructed taking into account metabolite T(2) relaxation times and dephasing from B(0) inhomogeneity. Metabolite concentrations were estimated by fitting the basis one-dimensional CT-PRESS diagonal magnitude spectra to the measured spectrum. This method was first validated using seven custom-built phantoms containing variable metabolite concentrations, and then applied to in vivo data acquired in rats exposed to vaporized ethanol and controls. Separate metabolite quantification revealed increased Gln after 16 weeks and increased Glu after 24 weeks of vaporized ethanol exposure in ethanol-treated compared with control rats. Without separate quantification, the signal from the combined resonances of Glu and Gln (Glx) showed an increase at both 16 and 24 weeks in ethanol-exposed rats, precluding the determination of the independent and differential contribution of each metabolite at each time.

Slice-selective broadband refocusing pulses for the robust generation of crushed spin-echoes

A major challenge for in vivo magnetic resonance spectroscopy with point-resolved spectroscopy (PRESS) is the low signal intensity for the measurement of weakly scalar coupled spins, for example lactate. The chemical-shift displacement error between the two coupling partners of the lactate molecule leads to a signal decrease. The chemical-shift displacement error is decreased and therefore the lactate signal is increased by using refocusing pulses with a broad bandwidth. Previously, slice-selective broadband universal rotation pulses (S-BURBOP) were designed and applied as refocusing pulses in the PRESS pulse sequence (Janich MA, et al., Journal of Magnetic Resonance, 2011, 213, 126-135). However, S-BURBOP pulses leave a phase error across the slice which is superimposed on the spectra when spatially resolving the PRESS voxel. In the present novel design of slice-selective broadband refocusing pulses (S-BREBOP) this phase error is avoided. S-BREBOP pulses obtain 2.5 times the bandwidth of conventional Shinnar-Le Roux pulses and are robust against ±20% miscalibration of the B(1) amplitude. S-BREBOP pulses were validated in phantoms and in a low-grade brain tumor of a patient. Compared to conventional Shinnar-Le Roux pulses they lead to a decrease of the chemical-shift displacement error and consequently a lactate signal increase.

Center-out echo-planar spectroscopic imaging with correction of gradient-echo phase and time shifts

A procedure to prevent the formation of image and spectral Nyquist ghosts in echo-planar spectroscopic imaging is introduced. It is based on a novel Cartesian center-out echo-planar spectroscopic imaging trajectory, referred to as EPSICO, and combined with a correction of the gradient-echo phase and time shifts. Processing of homogenous sets of forward and reflected echoes is no longer necessary, resulting in an optimized spectral width. The proposed center-out trajectory passively prevents the formation of Nyquist ghosts by privileging the acquisition of the center k-space line with forward echoes at the beginning of an echo-planar spectroscopic imaging dwell time and by ensuring that all k-space lines and their respective complex conjugates are acquired at equal time intervals. With the proposed procedure, concentrations of N-acetyl aspartate, creatine, choline, glutamate, and myo-inositol were reliably determined in human white matter at 3 T.

MR proton spectroscopy for myocardial lipid deposition quantification: a quantitative comparison between 1.5T and 3T

PURPOSE

To evaluate 3T magnetic resonance spectroscopy (MRS)-derived myocardial fat-signal fractions in comparison with those from 1.5T MRS.

MATERIALS AND METHODS

We conducted phantom, ex vivo and in vivo myocardial specimen evaluations at both 1.5T and 3T using (1)H-MRS. A phantom with nine fat-water emulsions was constructed to assess the accuracy of the spectroscopy measurements. Ex vivo spectroscopy data were acquired in 70 segments from 21 autopsy heart slices. In vivo spectroscopy data were acquired in the interventricular septum from 22 human volunteers.

RESULTS

Phantom experiments demonstrated that 1.5T and 3T measurements were highly correlated with the reference values (r = 0.78, P < 0.05). The ex vivo and in vivo experiments demonstrated an increase in signal-to-noise ratio (SNR) of 45 ± 73% and 76 ± 72% at 3T compared to 1.5T (P < 0.05). The mean fat-signal fraction was similar at 3T and 1.5T (1.11 ± 1.18 vs. 1.00 ± 1.09, respectively, P = NS) in ex vivo studies but were significantly different in the in vivo studies (2.47 ± 1.46 vs. 1.56 ± 1.34, P < 0.05). The fat-signal fractions from 3T and 1.5T correlated fairly well in all experiments.

CONCLUSION

3T MRS has significantly greater SNR and could potentially be more accurate as compared to 1.5T for quantification of myocardial fat fraction in in vivo studies.

Comparison of spectral fitting methods for overlapping J-coupled metabolite resonances

There is increasing interest in the use of two-dimensional J-resolved spectroscopic acquisition (multiecho) methods for in vivo proton magnetic resonance spectroscopy due to the improved discrimination of overlapping J-coupled multiplet resonances that is provided. Of particular interest is the potential for discrimination of the overlapping resonances of glutamate and glutamine. In this study, a new time-domain parametric spectral model that makes use of all available data is described for fitting the complete two-dimensional multiecho data, and the performance of this method was compared with fitting of one-dimensional spectra obtained following averaging multiecho data (echo time-averaged) and single-echo time PRESS (Point Resolved Spectroscopy) acquired spectra. These methods were compared using data obtained from a phantom containing typical brain metabolites and a human brain. Results indicate that improved performance and accuracy is obtained for the two-dimensional acquisition and spectral fitting model.

Field strength dependence of PRESS timings for simultaneous detection of glutamate and glutamine from 1.5 to 7T

An optimization of the PRESS sequence for magnetic resonance spectroscopy is presented to simultaneously detect the important brain metabolites of glutamate (Glu) and glutamine (Gln) at field strengths of 1.5, 3, 4.7, and 7T. Standard, clinical examinations typically use short echo times which in general are not ideal for the separation of Glu and Gln. The optimization procedure is based on numerical product operator simulations to produce yield and overlap measurements for all possible practical choices of PRESS inter-echo timings. The simulations illustrate the substantial modulations in Glu and Gln with field strength. At all field strengths, the optimized timings demonstrate a significant reduction in overlap compared to short echo PRESS, while maintaining a high metabolite signal, with Glu and Gln yields >90% when excluding T2 relaxation losses. Minimal overlap was attained at 7T (0.3% Gln contamination in the Glu signal), and 4.7T (1.2%). The optimized timings were applied in vivo on healthy volunteers at field strengths of 1.5 and 4.7T.

Optimized glutamate detection at 3T

PURPOSE

To identify the pulse sequence and acquisition parameters that result in the most accurate and repeatable measurements of glutamate (Glu) concentration in the brain at 3T.

MATERIALS AND METHODS

Simulations were performed to compare the accuracy and repeatability of 11 pulse sequences and acquisition parameters, within four general classes (PRESS, STEAM, Carr-Purcell PRESS [CPRESS] and TE averaged PRESS [JPRESS]), the majority of which were previously suggested as optimal for Glu detection. Three of the simulated acquisitions were implemented in a clinical scanner and measures of repeatability in vivo were compared to their simulated values.

RESULTS

Good agreement was demonstrated between simulated and experimentally determined measures of repeatability. Among the acquisitions considered, a CPRESS sequence with minimal echo time, together with, possibly, a short TE PRESS sequence, result in the most repeatable within session Glu measurements, while slightly overestimating the Glu concentration. Excellent accuracy is demonstrated by the simulations for a JPRESS sequence, at the expense of lower repeatability than optimal PRESS or CPRESS sequences.

CONCLUSION

Further proof of concept is presented toward validation of a simulation approach to understand pulse sequence performance in measuring the concentration of a given metabolite. Improved within session Glu measurement repeatability is predicted for CPRESS and demonstrated in vivo.

Spectral editing of weakly coupled spins using variable flip angles in PRESS constant echo time difference spectroscopy: application to GABA

A general in vivo magnetic resonance spectroscopy editing technique is presented to detect weakly coupled spin systems through subtraction, while preserving singlets through addition, and is applied to the specific brain metabolite gamma-aminobutyric acid (GABA) at 4.7 T. The new method uses double spin echo localization (PRESS) and is based on a constant echo time difference spectroscopy approach employing subtraction of two asymmetric echo timings, which is normally only applicable to strongly coupled spin systems. By utilizing flip angle reduction of one of the two refocusing pulses in the PRESS sequence, we demonstrate that this difference method may be extended to weakly coupled systems, thereby providing a very simple yet effective editing process. The difference method is first illustrated analytically using a simple two spin weakly coupled spin system. The technique was then demonstrated for the 3.01 ppm resonance of GABA, which is obscured by the strong singlet peak of creatine in vivo. Full numerical simulations, as well as phantom and in vivo experiments were performed. The difference method used two asymmetric PRESS timings with a constant total echo time of 131 ms and a reduced 120 degrees final pulse, providing 25% GABA yield upon subtraction compared to two short echo standard PRESS experiments. Phantom and in vivo results from human brain demonstrate efficacy of this method in agreement with numerical simulations.

Quantification of J-resolved proton spectra in two-dimensions with LCModel using GAMMA-simulated basis sets at 4 Tesla

A two-dimensional, J-resolved magnetic resonance spectroscopic extraction approach was developed employing GAMMA-simulated, LCModel basis-sets. In this approach, a two-dimensional J-resolved (2D-JPRESS) dataset was resolved into a series of one-dimensional spectra where each spectrum was modeled and fitted with its theoretically customized LCModel template. Metabolite levels were derived from the total integral across the J-series of spectra for each metabolite. Phantoms containing physiologic concentrations of the major brain chemicals were used for validation. Varying concentrations of glutamate and glutamine were evaluated at and around their accepted in vivo concentrations in order to compare the accuracy and precision of our method with 30 ms PRESS. We also assessed 2D-JPRESS and 30 ms PRESS in vivo, in a single voxel within the parieto-occipital cortex by scanning ten healthy volunteers once and a single healthy volunteer over nine repeated measures. Phantom studies demonstrated that serial fitting of 2D-JPRESS spectra with simulated LCModel basis sets provided accurate concentration estimates for common metabolites including glutamate and glutamine. Our in vivo results using 2D-JPRESS suggested superior reproducibility in measuring glutamine and glutamate relative to 30 ms PRESS. These novel methods have clear implications for clinical and research studies seeking to understand neurochemical dysfunction.

Effect of PRESS and STEAM sequences on magnetic resonance spectroscopic liver fat quantification

PURPOSE

To compare PRESS and STEAM MR spectroscopy for assessment of liver fat in human subjects.

MATERIALS AND METHODS

Single-voxel (20 x 20 x 20 mm) PRESS and STEAM spectra were obtained at 1.5T in 49 human subjects with known or suspected fatty liver disease. PRESS and STEAM sequences were obtained with fixed TR (1500 msec) and different TE (five PRESS spectra between TE 30-70 msec, five STEAM spectra between TE 20-60 msec). Spectra were quantified and T2 and T2-corrected peak area were calculated by different techniques. The values were compared for PRESS and STEAM.

RESULTS

Water T2 values from PRESS and STEAM were not significantly different (P = 0.33). Fat peak T2s were 25%-50% shorter on PRESS than on STEAM (P < 0.02 for all comparisons) and there was no correlation between T2s of individual peaks. PRESS systematically overestimated the relative fat peak areas (by 7%-263%) compared to STEAM (P < 0.005 for all comparisons). The peak area given by PRESS was more dependent on the T2-correction technique than STEAM.

CONCLUSION

Measured liver fat depends on the MRS sequence used. Compared to STEAM, PRESS underestimates T2 values of fat, overestimates fat fraction, and provides a less consistent fat fraction estimate, probably due to J coupling effects.

RF pulses for in vivo spectroscopy at high field designed under conditions of limited power using optimal control

Localized in vivo spectroscopy at high magnetic field strength (>3T) is susceptible to localization artifacts such as the chemical shift artifact and the spatial interference artifact for J-coupled spins. This latter artifact results in regions of anomalous phase for J-coupled spins. These artifacts are exacerbated at high magnetic field due to the increased frequency dispersion, coupled with the limited RF pulse bandwidths used for localization. Approaches to minimize these artifacts include increasing the bandwidth of the frequency selective excitation pulses, and the use of frequency selective saturation pulses to suppress the signals in the regions with anomalous phase. The goal of this article is to demonstrate the efficacy of optimal control methods to provide broader bandwidth frequency selective pulses for in vivo spectroscopy in the presence of limited RF power. It is demonstrated by examples that the use of optimal control methods enable the generation of (i) improved bandwidth selective excitation pulses, (ii) more efficient selective inversion pulses to be used for generation of spin echoes, and (iii) improved frequency selective saturation pulses. While optimal control also allows for the generation of frequency selective spin echo pulses, it is argued that it is more efficient to use dual inversion pulses for broadband generation of spin echoes. Finally, the optimal control routines and example RF pulses are made available for downloading.

Numerical simulations of localized high field 1H MR spectroscopy

The limited bandwidths of volume selective RF pulses in localized in vivo MRS experiments introduce spatial artifacts that complicate spectral quantification of J-coupled metabolites. These effects are commonly referred to as a spatial interference or "four compartment" artifacts and are more pronounced at higher field strengths. The main focus of this study is to develop a generalized approach to numerical simulations that combines full density matrix calculations with 3D localization to investigate the spatial artifacts and to provide accurate prior knowledge for spectral fitting. Full density matrix calculations with 3D localization using experimental pulses were carried out for PRESS (TE=20, 70 ms), STEAM (TE=20, 70 ms) and LASER (TE=70 ms) pulse sequences and compared to non-localized simulations and to phantom solution data at 4 T. Additional simulations at 1.5 and 7 T were carried out for STEAM and PRESS (TE=20 ms). Four brain metabolites that represented a range from weak to strong J-coupling networks were included in the simulations (lactate, N-acetylaspartate, glutamate and myo-inositol). For longer TE, full 3D localization was necessary to achieve agreement between the simulations and phantom solution spectra for the majority of cases in all pulse sequence simulations. For short echo time (TE=20 ms), ideal pulses without localizing gradients gave results that were in agreement with phantom results at 4 T for STEAM, but not for PRESS (TE=20). Numerical simulations that incorporate volume localization using experimental RF pulses are shown to be a powerful tool for generation of accurate metabolic basis sets for spectral fitting and for optimization of experimental parameters.

Spatial effects in the detection of gamma-aminobutyric acid: improved sensitivity at high fields using inner volume saturation

The MEGA-PRESS-IVS method has been developed, which combines MEGA (a frequency-selective editing technique) editing with the point-resolved spectroscopy sequence (PRESS) and inner volume saturation (IVS) localization, reducing the deleterious effects of spatial variation in coupling evolution. The IVS method has been previously described for improved efficiency of lactate detection. The current study demonstrates that the combination of MEGA-PRESS with IVS results in increased sensitivity for edited single-voxel measurements of glutamate and gamma-aminobutyric acid (GABA). A four-compartment model of coupling evolution is investigated through simple product operators and full spin-system simulations and the predicted pattern of signal evolution is demonstrated through MEGA-PRESS-MRSI. MEGA-PRESS-IVS is then compared to MEGA-PRESS in a phantom and an average signal increase of 24% is demonstrated in five healthy volunteers.

A detailed analysis of localized J-difference GABA editing: theoretical and experimental study at 4 T

The problem of low signal-to-noise ratio for gamma-aminobutyric acid (GABA) in vivo is exacerbated by inefficient detection schemes and non-optimal experimental parameters. To analyze the mechanisms for GABA signal loss of a MEGA-PRESS J-difference sequence at 4 T, numerical simulations were performed ranging from ideal to realistic experimental implementation, including volume selection and experimental radio frequency (RF) pulse shapes with a macromolecular minimization scheme. The simulations were found to be in good agreement with phantom and in vivo data from human brain. The overall GABA signal intensity for the simulations with realistic conditions for the MEGA-PRESS difference spectrum was calculated to be almost half of the signal simulated under ideal conditions (~43% signal loss). In contrast, creatine was reduced significantly less then GABA (~19% signal loss). The 'four-compartment' distribution due to J-coupling in the PRESS-based localization was one of the most significant sources of GABA signal loss, in addition to imperfect RF profiles for volume selection and editing. An alternative strategy that reduces signal loss due to the four-compartment distribution is suggested. In summary, a detailed analysis of J-difference editing is provided with estimates of the relative amounts of GABA signal losses due to various mechanisms. The numerical simulations presented in this study should facilitate both implementation of the more efficient acquisition and quantification process of J-coupled systems.

Proton MR spectroscopy of the prostate

PURPOSE

To summarize current technical and biochemical aspects and clinical applications of proton magnetic resonance spectroscopy (MRS) of the human prostate in vivo.

MATERIAL AND METHODS

Pertinent radiological and biochemical literature was searched and retrieved via electronic media (medline, pubmed. Basic concepts of MRS of the prostate and its clinical applications were extracted.

RESULTS

Clinical MRS is usually based on point resolved spectroscopy (PRESS) or spin echo (SE) sequences, along with outer volume suppression of signals from outside of the prostate. MRS of the prostate detects indicator lines of citrate, choline, and creatine. While healthy prostate tissue demonstrates high levels of citrate and low levels of choline that marks cell wall turnover, prostate cancer utilizes citrate for energy metabolism and shows high levels of choline. The ratio of (choline+creatine)/citrate distinguishes between healthy tissue and prostate cancer. Particularly when combined with magnetic resonance (MR) imaging, three-dimensional MRS imaging (3D-CSI, or 3D-MRSI) detects and localizes prostate cancer in the entire prostate with high sensitivity and specificity. Combined MR imaging and 3D-MRSI exceed the sensitivity and specificity of sextant biopsy of the prostate. When MRS and MR imaging agree on prostate cancer presence, the positive predictive value is about 80-90%. Distinction between healthy tissue and prostate cancer principally is maintained after various therapeutic treatments, including hormone ablation therapy, radiation therapy, and cryotherapy of the prostate.

CONCLUSIONS

Since it is non-invasive, reliable, radiation-free, and essentially repeatable, combined MR imaging and 3D-MRSI of the prostate lends itself to the planning of biopsy and therapy, and to post-therapeutic follow-up. For broad clinical acceptance, it will be necessary to facilitate MRS examinations and their evaluation and make MRS available to a wider range of institutions.

Simultaneous detection of resolved glutamate, glutamine, and gamma-aminobutyric acid at 4 T

A new approach is introduced to simultaneously detect resolved glutamate (Glu), glutamine (Gln), and gamma-aminobutyric acid (GABA) using a standard STEAM localization pulse sequence with the optimized sequence timing parameters. This approach exploits the dependence of the STEAM spectra of the strongly coupled spin systems of Glu, Gln, and GABA on the echo time TE and the mixing time TM at 4 T to find an optimized sequence parameter set, i.e., {TE, TM}, where the outer-wings of the Glu C4 multiplet resonances around 2.35 ppm, the Gln C4 multiplet resonances around 2.45 ppm, and the GABA C2 multiplet resonance around 2.28 ppm are significantly suppressed and the three resonances become virtual singlets simultaneously and thus resolved. Spectral simulation and optimization were conducted to find the optimized sequence parameters, and phantom and in vivo experiments (on normal human brains, one patient with traumatic brain injury, and one patient with brain tumor) were carried out for verification. The results have demonstrated that the Gln, Glu, and GABA signals at 2.2-2.5 ppm can be well resolved using a standard STEAM sequence with the optimized sequence timing parameters around {82 ms,48 ms} at 4 T, while the other main metabolites, such as N-acetyl aspartate (NAA), choline (tCho), and creatine (tCr), are still preserved in the same spectrum. The technique can be easily implemented and should prove to be a useful tool for the basic and clinical studies associated with metabolism of Glu, Gln, and/or GABA.

GAVA: spectral simulation for in vivo MRS applications

An application that provides a flexible and easy to use interface to the GAMMA spectral simulation package is described that is targeted at investigations using in vivo MR spectroscopic methods. The program makes available a number of widely used spatially localized MRS pulse sequences and NMR parameters for commonly observed tissue metabolites, enabling spectra to be simulated for any pulse sequence parameter and viewed in an integrated display. The application is interfaced with a database for storage of all simulation parameters and results of the simulations. This application provides a convenient method for generating a priori spectral information used in parametric spectral analyses and for visual examination of the effects of difference pulse sequences and parameter settings.

Elimination of spatial interference in PRESS-localized editing spectroscopy

Unambiguous detection of gamma-amino butyric acid (GABA) in the human brain is hindered by low concentration and spectral overlap with other metabolites. The popular MEGA-PRESS (PRESS: point-resolved spectroscopic sequence) method allows spectral separation of GABA from other metabolites, but suffers from a significant signal-to-noise ratio (SNR) reduction due to the 4-compartment artifact. An alternative PRESS localization technique (PRESS+4) was investigated and compared to MEGA-PRESS using numerical simulations, phantom, and in vivo experiments. It was shown that while the MEGA-PRESS method suffers significant signal loss ( approximately equal 20% for the difference spectrum), GABA signal intensity in PRESS+4 is reduced by only 2% compared to the nonlocalized condition at 4T. The improved method retains important features of the popular MEGA-PRESS such as additional water suppression and macromolecular elimination as demonstrated in human brain experiments. This method is not limited to GABA J-difference editing, but can be applied in any PRESS-based experiments. It should prove particularly useful at higher field, where the 4-compartment artifact is especially detrimental.

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.

ProFit: two-dimensional prior-knowledge fitting of J-resolved spectra

A two-dimensional fitting procedure is introduced, capable of extracting the full amount of information from 2D J-resolved magnetic resonance spectroscopic data. The fitting procedure uses a linear combination of 2D model spectra. For reducing the degrees of freedom and increasing robustness, it is divided into a non-linear outer loop and an inner linear least-squares fit for the concentrations. In vitro and in vivo experiments on a whole-body 3 T MR scanner show the detectability of a wide range of metabolites in the human brain, namely total creatine, N-acetylaspartate, N-acetylaspartylglutamate, choline-containing compounds, glutamate, myo-inositol, glutathione, scyllo-inositol, gamma-aminobutyric acid, alanine and ascorbic acid.