Strategy for the spectral filtering of myo-inositol and other strongly coupled spins

Short-acquisition-time JPRESS and its application to paediatric brain tumours

OBJECTIVE

To develop and assess a short-duration JPRESS protocol for detection of overlapping metabolite biomarkers and its application to paediatric brain tumours at 3 Tesla.

MATERIALS AND METHODS

The short-duration protocol (6 min) was optimised and compared for spectral quality to a high-resolution (38 min) JPRESS protocol in a phantom and five healthy volunteers. The 6-min JPRESS was acquired from four paediatric brain tumours and compared with short-TE PRESS.

RESULTS

Metabolite identification between the 6- and 38-min protocols was comparable in phantom and volunteer data. For metabolites with Cramer-Rao lower bounds > 50%, interpretation of JPRESS increased confidence in assignment of lactate, myo-Inositol and scyllo-Inositol. JPRESS also showed promise for the detection of glycine and taurine in paediatric brain tumours when compared to short-TE MRS.

CONCLUSION

A 6-min JPRESS protocol is well tolerated in paediatric brain tumour patients. Visual inspection of a 6-min JPRESS spectrum enables identification of a range of metabolite biomarkers of clinical interest.

Distinguishing neuronal from astrocytic subcellular microstructures using in vivo Double Diffusion Encoded 1H MRS at 21.1 T

Measuring cellular microstructures non-invasively and achieving specificity towards a cell-type population within an interrogated in vivo tissue, remains an outstanding challenge in brain research. Magnetic Resonance Spectroscopy (MRS) provides an opportunity to achieve cellular specificity via the spectral resolution of metabolites such as N-Acetylaspartate (NAA) and myo-Inositol (mI), which are considered neuronal and astrocytic markers, respectively. Yet the information typically obtained with MRS describes metabolic concentrations, diffusion coefficients or relaxation rates rather than microstructures. Understanding how these metabolites are compartmentalized is a challenging but important goal, which so far has been mainly addressed using diffusion models. Here, we present direct in vivo evidence for the confinement of NAA and mI within sub-cellular components, namely, the randomly oriented process of neurons and astrocytes, respectively. Our approach applied Relaxation Enhanced MRS at ultrahigh (21.1 T) field, and used its high ¹H sensitivity to measure restricted diffusion correlations for NAA and mI using a Double Diffusion Encoding (DDE) filter. While very low macroscopic anisotropy was revealed by spatially localized Diffusion Tensor Spectroscopy, DDE displayed characteristic amplitude modulations reporting on confinements in otherwise randomly oriented anisotropic microstructures for both metabolites. This implies that for the chosen set of parameters, the DDE measurements had a biased sensitivity towards NAA and mI sited in the more confined environments of neurites and astrocytic branches, than in the cell somata. These measurements thus provide intrinsic diffusivities and compartment diameters, and revealed subcellular neuronal and astrocytic morphologies in normal in vivo rat brains. The relevance of these measurements towards human applications—which could in turn help understand CNS plasticity as well as diagnose brain diseases—is discussed.

On the Utility of Short Echo Time (TE) Single Voxel 1H-MRS in Non-Invasive Detection of 2-Hydroxyglutarate (2HG); Challenges and Potential Improvement Illustrated with Animal Models Using MRUI and LCModel

Mutations in isocitrate dehydrogenase 1 and 2 (IDH1/2) are frequently found in brain tumors, and the resulting onco–metabolite, 2–hydroxyglutarate (2HG), has been suggested to be a potential diagnostic and prognostic biomarker of the diseases. Indeed, recent studies have demonstrated the feasibility of non–invasively detecting 2HG by using proton magnetic resonance spectroscopy (1H–MRS). Due to severe spectral overlaps of 2HG with its background metabolites and spectral baselines, however, the majority of those previous studies employed spectral editing methods with long echo times (TEs) instead of the most commonly used short TE approach with spectral fitting. Consequently, the results obtained with spectral editing methods may potentially be prone to errors resulting from substantial signal loss due to relaxation. Given that the spectral region where the main signal of 2HG resides is particularly sensitive to spectral baseline in metabolite quantification, we have investigated the impact of incorporating voxel–specifically measured baselines into the spectral basis set on the performance of the conventional short TE approach in 2HG detection in rodent models (Fisher 344 rats; n = 19) of IDH1/2 mutant–overexpressing F98 glioma at 9.4T. Metabolite spectra were acquired (SPECIAL sequence) for a tumor region and the contralateral normal region of the brain for each animal. For the estimation of spectral baselines metabolite–nulled spectra were obtained (double–inversion–recovery SPECIAL sequence) for each individual voxels. Data were post–processed with and without the measured baselines using MRUI and LCModel—the two most widely used data post–processing packages. Our results demonstrate that in–vivo detection of 2HG using the conventional short TE approach is challenging even at 9.4T. However, incorporation of voxel–specifically measured spectral baselines may potentially improve its performance. Upon more thorough validation in a larger number of animals and more importantly in human patients, the potential utility of the proposed short TE acquisition with voxel–specific baseline measurement approach in 2HG detection may need to be considered in the study design.

Experiments for the editing of singlet peaks and simplification of ¹H NMR spectra of complex mixtures

A 1D NMR singlet filtered experiment (s-filter) is presented for the clean editing of singlets peaks in crowded ¹H NMR spectra. The experiment affords a clean editing of singlet peaks and weak coupled signals ((n)JHH < 3 Hz), while the remaining multiplets peaks are strongly attenuated to a factor over 90% of its original intensity. The experiment provided a useful simplification of the typically congested ¹H NMR spectra of foodstuffs, and the edited signals retained most of the sensitivity and quantitativity. The versatility of the s-filter scheme permitted its easy incorporation as excitation building block into other 1D and 2D experiments, providing a further editing of the singlet responses. The combination of the 1D s-filter experiment and difference spectroscopy provided a 1D reverse singlet filtered experiment (s-filter(rev)) for the clean editing of multiplet peaks and strong suppression of singlets. The experiments proposed could be valuable to facilitate the detection and/or quantification of target compounds in complex mixtures of metabolites.

Enhancement of spectral editing efficacy of multiple quantum filters in in vivo proton magnetic resonance spectroscopy

The performance of multiple quantum filters (MQFs) can be disappointing when the background signal also arises from coupled spins. Moreover, at 3.0 T and even higher fields the majority of the spin systems of key brain metabolites fall into the strong-coupling regime. In this manuscript we address comprehensively, the importance of the phase of the multiple quantum coherence-generating pulse (MQ-pulse) in the design of MQFs, using both product operator and numerical analysis, in both zero and double quantum filter designs. The theoretical analyses were experimentally validated with the examples of myo-inositol editing and the separation of glutamate from glutamine. The results demonstrate that the phase of the MQ-pulse per se provides an additional spectral discrimination mechanism based on the degree of coupling beyond the conventional level-of-coherence approach of MQFs. To obtain the best spectral discrimination of strongly-coupled spin systems, therefore, the phase of the MQ-pulse must be included in the portfolio of the sequence parameters to be optimized.

Role of proton magnetic resonance spectroscopy in differentiating oligodendrogliomas from astrocytomas

BACKGROUND AND PURPOSE

Preoperative differentiation of astrocytomas from oligodendrogliomas is clinically important, as oligodendrogliomas are more sensitive to chemotherapy. The purpose of this study was to assess the role of proton magnetic resonance spectroscopy in distinguishing astrocytomas from oligodendrogliomas.

METHODS

Forty-six patients [astrocytomas (n= 17) and oligodendrogliomas (n= 29)] underwent magnetic resonance imaging and multi voxel proton magnetic resonance spectroscopic imaging before treatment. Peak areas for N-acetylaspartate (NAA), creatine (Cr), choline (Cho), myo-inositol (mI), glutamate/glutamine (Glx), and lipids + lactate (Lip+Lac) were analyzed from voxels that exhibited hyperintensity on fluid-attenuated inversion recovery images and were normalized to Cr from each voxel. The average metabolite/Cr ratios from these voxels were then compared between astrocytomas and oligodendrogliomas. Receiver-operating curve analyses were used as measures of differentiation accuracy of metabolite ratios. A threshold value for a metabolite ratio was estimated by maximizing the sum of sensitivity and specificity.

RESULTS

A significant difference in mI/Cr was observed between astrocytomas and oligodendrogliomas (.50 +/- .18 vs. 0.66 +/- 0.20, P < .05). Using a threshold value of .56 for mI/Cr ratio, it was possible to differentiate oligodendrogliomas from astrocytomas with a sensitivity of 72.4% and specificity of 76.4%.

CONCLUSION

These results suggest that mI/Cr might aid in distinguishing oligodendrogliomas from astrocytomas.

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.

Difference spectroscopy using PRESS asymmetry: application to glutamate, glutamine, and myo-inositol

A simple, clinically viable technique utilizing PRESS and strong coupling properties is presented for discrimination of coupled brain metabolites. The method relies on signal variation due to alteration of inter-echo timings (PRESS asymmetry) while maintaining a constant total echo time. Spin response of singlets and weakly coupled spins is unchanged due to PRESS asymmetry, allowing difference spectroscopy to detect unobstructed strongly coupled resonances. No changes to the standard PRESS sequence are required except variation of inter-echo timings. The procedure is illustrated for the separate detection of glutamate from glutamine and the detection of myo-inositol in simulation, phantom, and in vivo experiments at 4.7 T. The subtraction yields calculated from the simulation were 53% for glutamate and 75% for myo-inositol, and a resultant contribution of 96% glutamate to the total glutamate/glutamine multiplet in the 2.04-2.14 ppm range. To extend the treatment to other field strengths and metabolites, an analytical approximation based on a strongly coupled AB system was used to model individual spin groups. Subtraction spectroscopy yields for different combinations of coupling parameters were calculated for the detection of various strongly coupled metabolites at common clinical field strengths. The approximation also predicts adequate glutamate/glutamine discrimination at 3.0 T using the difference spectroscopy method.

Contamination of single-voxel multiple quantum filters by external water signals arising from intermolecular multiple quantum coherences

Multiple-quantum filtered pulse sequences simplify overlapping metabolite spectra by the elimination of peaks from uncoupled spin species, most notably from methyl groups and water, and the minimization of unwanted coupled-spin peaks. However, it is shown in this study that a significant contaminant water signal can pass through this family of filters in the form of intermolecular multiple-quantum coherences. An imaging evaluation of a single-voxel multiple quantum filter experiment confirms that the water contamination is excited from outside of the voxel of interest, thus having an increased potential for broad spectral contamination. Phantom and in vivo experiments at 3.0 T are used to illustrate, first, significant water contamination of a single-voxel double quantum filter experiment optimized for the observation of glutamate, and second, the elimination of the unwanted water signal with conventional phase cycling and optimized filter gradient orientations.

Which pulse sequence is optimal for myo-inositol detection at 3T?

Optimized myo-inositol (mI) detection is important for diagnosing and monitoring a multitude of pathological conditions of the brain. Simulations are presented in this work, performed to decide which pulse sequence has the most significant advantage in terms of improving repeatability and accuracy of mI measurements at 3T over the pulse sequence used typically in the clinic, a TE = 35 ms PRESS sequence. Five classes of pulse sequences, four previously suggested for optimized mI detection (a short TE PRESS, a Carr-Purcell PRESS sequence, an optimized STEAM sequence, an optimized zero quantum filter), and one optimized for mI detection in this work (a single quantum filter) were compared to a standard, TE = 35 ms pulse sequence. While limiting the SNR of an acquisition to the equivalent SNR of a spectrum acquired in 5 min from an 8 cc voxel, it was found through simulations that the most repeatable mI measurements would be obtained with a Carr-Purcell sequence. This sequence was implemented in a clinical scanner, and improved mI measurements were demonstrated in vivo.

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.

1H MR spectroscopic imaging with short and long echo time to discriminate glycine in glial tumours

OBJECT

To investigate glycine (Gly) concentrations in low- and high-grade gliomas based on (1)H MR spectroscopic imaging (MRSI) with short and long echo time (TE). Myoinositol (MI) and Gly appear at the same resonance frequency of 3.56 ppm, but due to strong coupling the MI signal dephases more rapidly. Therefore, their contribution to the 3.56 ppm signal should be distinguishable comparing MRSI data acquired at short and long TE.

MATERIALS AND METHODS

(1)H MRSI (TE = 30 and 144 ms) was performed at 3 T in 29 patients with histopathological confirmed World Health Organization (WHO) grade II-IV gliomas and in FIVE healthy subjects. All spectra from the gliomas revealed increase of the 3.56 ppm resonance in the short TE spectra. Signal intensities of Gly and MI were differentiated either by analysing the short to long TE ratio of the resonance or by performing a weighted difference. Gly concentrations were compared between high-grade (WHO III-IV) and low-grade gliomas.

RESULTS

High-grade gliomas showed significantly higher Gly concentrations compared to low-grade gliomas.

CONCLUSION

Appropriate data processing of short and long TE (1)H MRSI provides a tool to distinguish and to quantify Gly and MI concentrations in gliomas. As Gly seems to be a marker of malignancy, more dedicated spectroscopic methods to differentiate these metabolites are justified.

Incorporating homonuclear polarization transfer into PRESS for proton spectral editing: illustration with lactate and glutathione

A proton spectral editing pulse sequence for the detection of metabolites with spin systems that involve weak coupling is presented. The sequence is based on homonuclear polarization transfer incorporated into the standard PRESS (Point RESolved Spectroscopy) sequence, which is a volume-selective double spin echo method, to enable spatial localization. All peaks in the region of interest are initially suppressed whether they are peaks from the target metabolite or from contaminating background. The target signal is then restored by polarization transfer from a proton that has a resonance outside the suppressed region and to which the target spins are weakly coupled. This is achieved by the application of a 90 degrees hard pulse with phase orthogonal to that of the PRESS excitation pulse at the location of the first echo in PRESS and by optimizing the two PRESS timings, TE(1) and TE(2), for most efficient yield. Background signal not coupled to any protons outside the initially saturated region remains suppressed. The advantage of this sequence compared to multiple quantum filters is that signal from singlet peaks outside the suppressed area are preserved and can thus be used as a reference. The efficacy of the sequence was verified experimentally on phantom solutions of lactate and glutathione at 3.0 T. For the AX(3) spin system of lactate, the sequence timings were optimized by product operator calculations whereas for the ABX spin system of the cysteinyl group of glutathione numerical calculations were performed for sequence timing optimization.

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.

Detection of the myo-inositol 4.06-ppm resonance by selective J rewinding: application to human prefrontal cortex in vivo

A new proton NMR single-voxel spectral editing strategy for the rapid measurement of myo-inositol in human brain is proposed. The spectral editing detects the 4.06-ppm, weakly coupled resonance by means of selective J rewinding. An 84.6-ms-long quadruple-resonance selective 180 degrees radiofrequency pulse, implemented within an adiabatic-refocused localization sequence, induces an in-phase triplet at 4.06 ppm, while eliminating the contribution from creatine, phosphorylethanolamine, lactate, and serine in this spectral region. The myo-inositol concentration in human prefrontal cortex is estimated to be 5.7 +/- 0.9 mumol/g (mean +/- SD, n = 7), with reference to NAA at 10 micromol/g.

Optimized diffusion-weighted spectroscopy for measuring brain glutamate apparent diffusion coefficient on a whole-body MR system

A diffusion-weighted stimulated echo acquisition mode sequence was implemented in order to measure the glutamate apparent diffusion coefficient (ADC) in the monkey brain on a whole-body 3 T system. TE and TM were adjusted for maximizing glutamate signal intensity. Glutamate ADC was measured in a 5.8 mL voxel made of gray and white matter in macaque monkeys. The effect of post-processing on the estimated ADC was carefully assessed and appeared to be critical. Individual scan phasing and macromolecule subtraction corrected for approximately 25% and approximately 15% biases in glutamate ADC, respectively. Proper data processing yielded ADC values of 0.21 +/- 0.03 microm(2)/ms for glutamate, 0.15 +/- 0.04 microm(2)/ms for N-acetylaspartate + N-acetylaspartylglutamate, 0.12 +/- 0.03 microm(2)/ms for creatine, 0.11 +/- 0.05 microm(2)/ms for choline and 0.18 +/- 0.04 microm(2)/ms for myo-inositol.

Decreased prefrontal Myo-inositol in major depressive disorder

BACKGROUND

Postmortem studies have shown robust prefrontal cortex glial losses and more subtle neuronal changes in major depressive disorder (MDD). Earlier proton magnetic resonance spectroscopy (1H-MRS) studies of the glial marker myo-inositol in MDD were subject to potential confounds. The primary hypothesis of this study was that MDD patients would show reduced prefrontal/anterior cingulate cortex levels of myo-inositol.

METHODS

Thirteen nonmedicated moderate-severe MDD patients and 13 matched control subjects were studied (six male, seven female per group). Proton magnetic resonance spectroscopy stimulated echo acquisition mode spectra (3.0 T; echo time=168 msec; mixing time=28 msec; repetition time=3000 msec) were obtained from prefrontal/anterior cingulate cortex. Metabolite data were adjusted for tissue composition.

RESULTS

Patients with MDD showed significantly lower myo-inositol/creatine ratios (.94+/-.23) than control subjects (1.32+/-.37) [F(1,23)=6.9; p=.016].

CONCLUSIONS

These data suggest a reduction of myo-inositol in prefrontal/anterior cingulate cortex in MDD, which could be a consequence of glial loss or altered glial metabolism. Additional in vivo studies of glial markers could add to the understanding of the pathophysiology of MDD.

Variability of metabolite yield using STEAM or PRESS sequences in vivo at 3.0 T, illustrated with myo-inositol

Using as an example the myo-inositol (mI) band at approximately 3.6 ppm in the proton spectrum from brain, an evaluation is presented that highlights the difficulties of quantifying metabolites with strongly coupled spins with either STEAM or PRESS and demonstrates some advantages of prospective sequence analysis when measuring their concentrations. The analysis emphasizes the variation in coupled-spin signal yield and lineshape, compared with that of uncoupled singlets such as N-acetylaspartate, a variation that differs from one metabolite spin system to another. This difference in variation between a target metabolite (e.g., mI) and its contaminating background metabolites (e.g., glutamate and taurine, etc.) is shown to provide in certain circumstances a substantial reduction in background contamination (both metabolite and macromolecule) while maintaining sufficient signal-to-noise ratio for precise quantification. For example, sequence times are demonstrated, both for STEAM and for PRESS, that, relative to the short echo-time sequences typical in the literature, enhance the signal to metabolite background of the 3.6-ppm band of mI by factors of 1.7 and 1.3, respectively, essentially eliminate the macromolecular baseline, and yet in vivo retain an S/N approximately 10 in both cases.

Effect of strong homonuclear proton coupling on localized13C detection using PRESS

The effect of strong homonuclear proton coupling on (13)C incorporation measurements by either indirect or direct means was investigated (and illustrated with glutamate) both numerically and experimentally at 3.0 T. In particular, two sequences were considered, each using a proton PRESS sequence for localization. The indirect (13)C detection method incorporated the POCE (proton observe carbon edited) technique onto PRESS, and for direct (13)C detection a DEPT (distortionless enhancement by polarization transfer) sequence was appended to the PRESS localization. Both analysis and experiment demonstrate that when strong homonuclear coupling of protons is additional to heteronuclear coupling with (13)C spins, the (13)C measures derived from either the indirect PRESS-POCE sequence or the direct-but-enhanced PRESS-DEPT sequence are significantly modified. Specifically, the MR lineshapes of both (13)C-bonded and nonbonded protons are changed during (13)C incorporation, giving rise, for example, to a potential cross-contamination of < or =30% between glutamate (13)C(3) and (13)C(4) measures from the PRESS-POCE indirect method. During direct-but-enhanced detection, the DEPT enhancement is reduced for glutamate (13)C(2), (13)C(3), and (13)C(4) but not equally, and the reduction is further exacerbated by proton PRESS localization, which gives rise to enhancements that are strong functions of PRESS TE(1) and TE(2).