Single-voxel MRS with prospective motion correction and retrospective frequency correction

Model-based frequency-and-phase correction of 1H MRS data with 2D linear-combination modeling

PURPOSE

Retrospective frequency-and-phase correction (FPC) methods attempt to remove frequency-and-phase variations between transients to improve the quality of the averaged MR spectrum. However, traditional FPC methods like spectral registration struggle at low SNR. Here, we propose a method that directly integrates FPC into a 2D linear-combination model (2D-LCM) of individual transients ("model-based FPC"). We investigated how model-based FPC performs compared to the traditional approach, i.e., spectral registration followed by 1D-LCM in estimating frequency-and-phase drifts and, consequentially, metabolite level estimates.

METHODS

We created synthetic in-vivo-like 64-transient short-TE sLASER datasets with 100 noise realizations at 5 SNR levels and added randomly sampled frequency and phase variations. We then used this synthetic dataset to compare the performance of 2D-LCM with the traditional approach (spectral registration, averaging, then 1D-LCM). Outcome measures were the frequency/phase/amplitude errors, the SD of those ground-truth errors, and amplitude Cramér Rao lower bounds (CRLBs). We further tested the proposed method on publicly available in-vivo short-TE PRESS data.

RESULTS

2D-LCM estimates (and accounts for) frequency-and-phase variations directly from uncorrected data with equivalent or better fidelity than the conventional approach. Furthermore, 2D-LCM metabolite amplitude estimates were at least as accurate, precise, and certain as the conventionally derived estimates. 2D-LCM estimation of FPC and amplitudes performed substantially better at low-to-very-low SNR.

CONCLUSION

Model-based FPC with 2D linear-combination modeling is feasible and has great potential to improve metabolite level estimation for conventional and dynamic MRS data, especially for low-SNR conditions, for example, long TEs or strong diffusion weighting.

Sequence building block for magnetic resonance spectroscopy on Siemens VE-series scanners

We present a sequence building block (SBB) that embeds magnetic resonance spectroscopy (MRS) into another sequence on the Siemens VE platform without any custom hardware. This enables dynamic studies such as functional MRS (fMRS), dynamic shimming and frequency correction, and acquisition of navigator images for motion correction. The SBB supports nonlocalised spectroscopy (free induction decay), STimulated Echo Acquisition Mode single voxel spectroscopy, and 1D, 2D and 3D phase-encoded chemical shift imaging. It can embed 1H or X-nuclear MRS into a 1H sequence; and 1H-MRS into an X-nuclear sequence. We demonstrate integration into the vendor's gradient-recalled echo sequence. We acquire test data in phantoms with three coils (31P/1H, 13C/1H and 2H/1H) and in two volunteers on a 7-T Terra MRI scanner. Fifteen lines of code are required to insert the SBB into a sequence. Spectra and images are acquired successfully in all cases in phantoms, and in human abdomen and calf muscle. Phantom comparison of signal-to-noise ratio and linewidth showed that the SBB has negligible effects on image and spectral quality, except that it sometimes produces a nuclear Overhauser effect (NOE) signal enhancement for multinuclear applications in line with conventional 1H NOE pulses. Our new SBB embeds MRS into a host imaging or spectroscopy sequence in 15 lines of code. It allows homonuclear and heteronuclear interleaving. The package is available through the standard C2P procedure. We hope this will lower the barrier for entry to studies applying dynamic fMRS and for online motion correction and B0-shim updating.

Fast high-resolution prospective motion correction for single-voxel spectroscopy

PURPOSE

To develop a fast high-resolution image-based motion correction method using spiral navigators with multislice-to-volume registration.

METHODS

A semi-LASER sequence was modified to include a multislice spiral navigator for prospective motion correction (∼305 ms including acquisition, processing, and feedback) as well as shim and frequency navigators for prospective shim and frequency correction (∼100 ms for each). MR spectra were obtained in the prefrontal cortex in five healthy subjects at 3 T with and without prospective motion and shim correction. The effect of key navigator parameters (number of slices, image resolution, and excitation flip angle) on registration accuracy was assessed using simulations.

RESULTS

Without prospective motion and shim correction, spectral quality degraded significantly in the presence of voluntary motion. In contrast, with prospective motion and shim correction, spectral quality was improved (metabolite linewidth = 6.7 ± 0.6 Hz, SNR= 67 ± 9) and in good agreement with baseline data without motion (metabolite linewidth = 6.9 ± 0.9 Hz, SNR = 73 ± 9). In addition, there was no significant difference in metabolites concentrations measured without motion and with prospective motion and shim correction in the presence of motion. Simulations showed that the registration precision was comparable when using three navigator slices with 3 mm resolution and when using the entire volume (all slices) with 8 mm resolution.

CONCLUSION

The proposed motion correction scheme allows fast and precise prospective motion and shim correction for single-voxel spectroscopy at 3 T. With 3 mm resolution, only a few navigator slices are necessary to achieve excellent motion correction performance.

Model-based frequency-and-phase correction of 1H MRS data with 2D linear-combination modeling

Purpose

Retrospective frequency-and-phase correction (FPC) methods attempt to remove frequency-and-phase variations between transients to improve the quality of the averaged MR spectrum. However, traditional FPC methods like spectral registration struggle at low SNR. Here, we propose a method that directly integrates FPC into a two-dimensional linear-combination model (2D-LCM) of individual transients ('model-based FPC'). We investigated how model-based FPC performs compared to the traditional approach, i.e., spectral registration followed by 1D-LCM in estimating frequency-and-phase drifts and, consequentially, metabolite level estimates.

Methods

We created synthetic in-vivo-like 64-transient short-TE sLASER datasets with 100 noise realizations at 5 SNR levels and added randomly sampled frequency and phase variations. We then used this synthetic dataset to compare the performance of 2D-LCM with the traditional approach (spectral registration, averaging, then 1D-LCM). Outcome measures were the frequency/phase/amplitude errors, the standard deviation of those ground-truth errors, and amplitude Cramér Rao Lower Bounds (CRLBs). We further tested the proposed method on publicly available in-vivo short-TE PRESS data.

Results

2D-LCM estimates (and accounts for) frequency-and-phase variations directly from uncorrected data with equivalent or better fidelity than the conventional approach. Furthermore, 2D-LCM metabolite amplitude estimates were at least as accurate, precise, and certain as the conventionally derived estimates. 2D-LCM estimation of frequency and phase correction and amplitudes performed substantially better at low-to-very-low SNR.

Conclusion

Model-based FPC with 2D linear-combination modeling is feasible and has great potential to improve metabolite level estimation for conventional and dynamic MRS data, especially for low-SNR conditions, e.g., long TEs or strong diffusion weighting.

Concurrent intrinsic optical imaging and fMRI at ultra-high field using magnetic field proof optical components

Intrinsic optical imaging (IOI) is a well established technique to quantify activation-related hemodynamical changes at the surface of the brain, which can be used to investigate the underlying processes of BOLD signal formation. To directly and quantitatively relate IOI and fMRI, simultaneous measurements with the two modalities are necessary. Here, a novel technical solution for a completely in-bore setup is presented, which uses only magnetic field proof components and thus allows concurrent recordings with a quality similar to that obtained in separate experiments. Measurements of the somatosensory cortex of rats with electrical forepaw stimulation were used to verify this approach. The high spatial and temporal resolution of the fMRI data, which is possible due to the high magnetic field of 14.1 T, the use of a point-spread function-based distortion correction and optimized additional anatomical images, allowed accurate colocalization of the images of the two modalities. Accordingly, detailed investigations of the temporal and spatial relationships between the hemodynamic parameters and the fMRI signal, which demonstrate the linear dependence of the BOLD effect on changes in the concentrations of oxygenated and deoxygenated hemoglobin, are possible. Comparisons between the signals emerging from arterial, venous and parenchymal areas are possible and show clearly distinct characteristics. The presented setup allows combining MRI measurements and optical recordings without serious losses in the data quality of either modality. While the proposed combination of fMRI and IOI can help to gain valuable insight into the generation of the BOLD effect, the setup can be easily modified to include different types of optical or MRI measurements.

Interleaved and simultaneous multi-nuclear magnetic resonance in vivo. Review of principles, applications and potential

Magnetic resonance signals from different nuclei can be excited or received at the same time,rendering simultaneous or rapidly interleaved multi-nuclear acquisitions feasible. The advan-tages are a reduction of total scan time compared to sequential multi-nuclear acquisitions or that additional information from heteronuclear data is obtained at thesame time and anatomical position. Information content can be qualitatively increased by delivering a more comprehensive MR-based picture of a transient state (such as an exercise bout). Also, combiningnon-proton MR acquisitions with ¹ Hinformation (e.g., dynamic shim updates and motion correction) can be used to improve data quality during long scans and benefits image coregistration. This work reviews the literature on interleaved and simultaneous multi-nuclear MRI and MRS in vivo. Prominent use cases for this methodology in clinical and research applications are brain and muscle, but studies have also been carried out in other targets, including the lung, knee, breast and heart. Simultaneous multi-nuclear measurements in the liver and kidney have also been performed, but exclusively in rodents. In this review, a consistent nomenclature is proposed, to help clarify the terminology used for this principle throughout the literature on in-vivo MR. An overview covers the basic principles, the technical requirements on the MR scanner and the implementations realised either by MR system vendors or research groups, from the early days until today. Considerations regarding the multi-tuned RF coils required and heteronuclear polarisation interactions are briefly discussed, and fields for future in-vivo applications for interleaved multi-nuclear MR pulse sequences are identified.

In-phase simultaneous spectral editing of lactate and alanine with suppression of J-coupled lipids by the modified selective multiple quantum coherence sequences

¹H magnetic resonance spectroscopy (MRS) with the multiple quantum coherence (MQC) technique allows for the detection of lactate, an end product of glycolysis, in the environment of lipids. The method can also be used to detect alanine, a byproduct of glutaminolysis. An issue is that when both lactate and alanine are detected together by the MQC technique, a phase mismatch arises between lactate and alanine signals due to off-resonance rotations and the difference in double quantum coherence frequencies between the two molecules. Such phase mismatch can cause errors in spectral fitting and metabolite quantification. In this study, we designed two pulse sequences that eliminate such phase differences of lactate and alanine while suppressing lipid signals by modifications of the Selective Multiple Quantum Coherence (Sel-MQC) sequence, a well-known MQC technique. Using the product operator formalism and the off-resonance rotation matrices, the phase evolutions of lactate and alanine during the spectrally selective pulses and the free precession times of the sequence at the single quantum, double quantum and zero quantum coherence states of these molecules were calculated. The multiple quantum (MQ) evolution time t₁ that can remove the phase difference of lactate and alanine at the echo was calculated and fine-tuned with experiments. The lactate and alanine signal intensities and the editing efficiencies from the two modified Sel-MQC sequences were theoretically predicted by using the product operator evolutions and compared with the experimental data. The J-coupled lipid signals were successfully suppressed by both sequences. One of the two developed sequences was applied to a human body with a phantom of lactate and alanine, which resulted in successful in-phase editing of lactate and alanine and suppression of the lipid signals from the body. The study sets an important foundation for the noninvasive detection of lactate and alanine from tumors of cancer patients.

Coherence pathway analysis of J-coupled lipids and lactate and effective suppression of lipids upon the selective multiple quantum coherence lactate editing sequence

Objective:The selective multiple quantum coherence (Sel-MQC) sequence is a magnetic resonance spectroscopy (MRS) technique used to detect lactate and suppress co-resonant lipid signalsin vivo. The coherence pathways of J-coupled lipids upon the sequence, however, have not been studied, hindering a logical design of the sequence to fully attenuate lipid signals. The objective of this study is to elucidate the coherence pathways of J-coupled lipids upon the Sel-MQC sequence and find a strategy to effectively suppress lipid signals from these pathways while keeping the lactate signal.Approach:The product operator formalism was used to express the evolutions of the J-coupled spins of lipids and lactate. The transformations of the product operators by the spectrally selective pulses of the sequence were calculated by using the off-resonance rotation matrices. The coherence pathways and the conversion rates of the individual pathways were derived from them. Experiments were performed on phantoms and two human subjects at 3T.Main results:The coherence pathways contributing to the various lipid resonance signals by the Sel-MQC sequence depending on the gradient ratios and RF pulse lengths were identified. Theoretical calculations of the signals from the determined coherence pathways and signal attenuations by gradients matched the experimental data very well. Lipid signals from fatty tissues of the subjects were successfully suppressed to the noise level by using the gradient ratio -0.8:-1:2 or 1:0.8:2. The new gradient ratios kept the lactate signal the same as with the previously used gradient ratio 0:-1:2.Significance:The study has elucidated the coherence pathways of J-coupled lipids upon the Sel-MQC sequence and demonstrated how lipid signals can be effectively suppressed while keeping lactate signals by using information from the coherence pathway analysis. The findings enable applying the Sel-MQC sequence to lactate detection in an environment of high concentrations of lipids.

Preprocessing, analysis and quantification in single-voxel magnetic resonance spectroscopy: experts' consensus recommendations

Once an MRS dataset has been acquired, several important steps must be taken to obtain the desired metabolite concentration measures. First, the data must be preprocessed to prepare them for analysis. Next, the intensity of the metabolite signal(s) of interest must be estimated. Finally, the measured metabolite signal intensities must be converted into scaled concentration units employing a quantitative reference signal to allow meaningful interpretation. In this paper, we review these three main steps in the post-acquisition workflow of a single-voxel MRS experiment (preprocessing, analysis and quantification) and provide recommendations for best practices at each step.

Motion correction methods for MRS: experts' consensus recommendations

Long acquisition times due to intrinsically low signal-to-noise ratio and the need for highly homogeneous B₀ field make MRS particularly susceptible to motion or scanner instability compared with MRI. Motion-induced changes in both localization and shimming (ie B₀ homogeneity) degrade MRS data quality. To mitigate the effects of motion three approaches can be employed: (1) subject immobilization, (2) retrospective correction, and (3) prospective real-time correction using internal and/or external tracking methods. Prospective real-time correction methods can simultaneously update localization and the B₀ field to improve MRS data quality. While localization errors can be corrected with both internal (navigators) and external (optical camera, NMR probes) tracking methods, the B₀ field correction requires internal navigator methods to measure the B₀ field inside the imaged volume and the possibility to update the scanner shim hardware in real time. Internal and external tracking can rapidly update the MRS localization with submillimeter and subdegree precision, while scanner frequency and first-order shims of scanner hardware can be updated by internal methods every sequence repetition. These approaches are most well developed for neuroimaging, for which rigid transformation is primarily applicable. Real-time correction greatly improves the stability of MRS acquisition and quantification, as shown in clinical studies on subjects prone to motion, including children and patients with movement disorders, enabling robust measurement of metabolite signals including those with low concentrations, such as gamma-aminobutyric acid and glutathione. Thus, motion correction is recommended for MRS users and calls for tighter integration and wider availability of such methods by MR scanner manufacturers.

Prospective frequency and motion correction for edited 1H magnetic resonance spectroscopy

The major inhibitory neurotransmitter gamma-aminobutyric acid (GABA) and the dominant antioxidant glutathione (GSH) both play a crucial role in brain functioning and are involved in several neurodegenerative and psychiatric diseases. Magnetic resonance spectroscopy (MRS) is a unique way to measure these neurometabolites non-invasively, but the measurement is highly sensitive to head movements, and especially in specific patient groups, motion stabilization in MRS could be valuable. Conventional MRS is acquired at relatively short echo times (TE), however, for unambiguous detection of GABA and GSH, spectral editing techniques are typically used. These depend on longer TEs and use frequency selective spectral editing pulses to separate the low-intensity peaks of GABA and GSH from overlapping resonances, but results in further increased motion sensitivity. Low-intensity metabolite peaks are usually edited one-by-one, however, simultaneous editing of multiple metabolites can be achieved using a Hadamard scheme, resulting in a substantial reduction in scan time. To investigate and correct for motion sensitivity in both conventional short-TE MRS (PRESS) and edited MRS (HERMES), we implemented a navigator-based prospective motion correction strategy including reacquisition of corrupted data. PRESS and HERMES spectra were acquired without motion, with motion with correction (repeated twice), and with motion without correction. Results indicate that when sufficient retrospective outlier removal is used, no significant differences in concentration and spectral quality were observed between motion conditions, even without prospective correction. HERMES spectral editing data showed to be more sensitive to motion, as significant differences in metabolite estimates and variability of spectral quality measures were observed for tCr, GABA+ and GSH when only retrospective outlier removal was applied. When using both prospective and retrospective correction, spectral quality was improved to almost the level of the no-motion acquisition. No differences in metabolite ratios for GABA and GSH could be observed when using motion correction. In conclusion, edited MRS showed to be more prone to motion artifacts, and prospective motion correction can restore most of the spectral quality in both conventional and edited MRS.

Motion correction in magnetic resonance spectroscopy

In vivo proton magnetic resonance spectroscopy and spectroscopic imaging (MRS/MRSI) are valuable tools to study normal and abnormal human brain physiology. However, they are sensitive to motion, due to strong crusher gradients, long acquisition times, reliance on high magnetic field homogeneity, and particular acquisition methods such as spectral editing. The effects of motion include incorrect spatial localization, phase fluctuations, incoherent averaging, line broadening, and ultimately quantitation errors. Several retrospective methods have been proposed to correct motion-related artifacts. Recent advances in hardware also allow prospective (real-time) correction of the effects of motion, including adjusting voxel location, center frequency, and magnetic field homogeneity. This article reviews prospective and retrospective methods available in the literature and their implications for clinical MRS/MRSI. In combination, these methods can attenuate or eliminate most motion-related artifacts and facilitate the acquisition of high-quality data in the clinical research setting.

Correcting frequency and phase offsets in MRS data using robust spectral registration

An algorithm for retrospective correction of frequency and phase offsets in MRS data is presented. The algorithm, termed robust spectral registration (rSR), contains a set of subroutines designed to robustly align individual transients in a given dataset even in cases of significant frequency and phase offsets or unstable lipid contamination and residual water signals. Data acquired by complex multiplexed editing approaches with distinct subspectral profiles are also accurately aligned. Automated removal of unstable lipid contamination and residual water signals is applied first, when needed. Frequency and phase offsets are corrected in the time domain by aligning each transient to a weighted average reference in a statistically optimal order using nonlinear least-squares optimization. The alignment of subspectra in edited datasets is performed using an approach that specifically targets subtraction artifacts in the frequency domain. Weighted averaging is then used for signal averaging to down-weight poorer-quality transients. Algorithm performance was assessed on one simulated and 67 in vivo pediatric GABA-/GSH-edited HERMES datasets and compared with the performance of a multistep correction method previously developed for aligning HERMES data. The performance of the novel approach was quantitatively assessed by comparing the estimated frequency/phase offsets against the known values for the simulated dataset or by examining the presence of subtraction artifacts in the in vivo data. Spectral quality was improved following robust alignment, especially in cases of significant spectral distortion. rSR reduced more subtraction artifacts than the multistep method in 64% of the GABA difference spectra and 75% of the GSH difference spectra. rSR overcomes the major challenges of frequency and phase correction.

Advanced single voxel 1 H magnetic resonance spectroscopy techniques in humans: Experts' consensus recommendations

Conventional proton MRS has been successfully utilized to noninvasively assess tissue biochemistry in conditions that result in large changes in metabolite levels. For more challenging applications, namely, in conditions which result in subtle metabolite changes, the limitations of vendor-provided MRS protocols are increasingly recognized, especially when used at high fields (≥3 T) where chemical shift displacement errors, B0 and B1 inhomogeneities and limitations in the transmit B1 field become prominent. To overcome the limitations of conventional MRS protocols at 3 and 7 T, the use of advanced MRS methodology, including pulse sequences and adjustment procedures, is recommended. Specifically, the semiadiabatic LASER sequence is recommended when TE values of 25-30 ms are acceptable, and the semiadiabatic SPECIAL sequence is suggested as an alternative when shorter TE values are critical. The magnetic field B0 homogeneity should be optimized and RF pulses should be calibrated for each voxel. Unsuppressed water signal should be acquired for eddy current correction and preferably also for metabolite quantification. Metabolite and water data should be saved in single shots to facilitate phase and frequency alignment and to exclude motion-corrupted shots. Final averaged spectra should be evaluated for SNR, linewidth, water suppression efficiency and the presence of unwanted coherences. Spectra that do not fit predefined quality criteria should be excluded from further analysis. Commercially available tools to acquire all data in consistent anatomical locations are recommended for voxel prescriptions, in particular in longitudinal studies. To enable the larger MRS community to take advantage of these advanced methods, a list of resources for these advanced protocols on the major clinical platforms is provided. Finally, a set of recommendations are provided for vendors to enable development of advanced MRS on standard platforms, including implementation of advanced localization sequences, tools for quality assurance on the scanner, and tools for prospective volume tracking and dynamic linear shim corrections.

Prospective motion and B0 shim correction for MR spectroscopy in human brain at 7T

PURPOSE

To demonstrate feasibility and performance of prospective motion and B₀ shim correction for MRS in human brain at 7T.

METHODS

Prospective motion correction using an optical camera and linear B₀ shim correction using FASTMAP-like navigators were implemented into a semi-LASER sequence. The effect of motion on spectral quality was assessed without and with prospective correction in prefrontal cortex in 11 subjects.

RESULTS

Without prospective motion and shim correction, motion resulted in considerable degradation of MR spectra (broader linewidth, lower signal-to-noise ratio, degraded water suppression). With prospective motion and shim correction, spectral quality remained excellent despite motion. Prospective motion correction alone was not sufficient to prevent degradation of spectral quality.

CONCLUSION

Prospective motion and B₀ shim correction is feasible at 7T and should help improve the robustness of MRS, particularly in motion-prone populations.

Effects of carrier frequency mismatch on frequency-selective spectral editing

OBJECTIVES

This study sought to investigate the effects of carrier frequency mismatch on spectral editing and its correction by frequency matching of basis functions.

MATERIALS AND METHODS

Full density matrix computations and Monte-Carlo simulations based on magnetic resonance spectroscopy (MRS) data collected from five healthy volunteers at 7 T were used to analyze the effects of carrier frequency mismatch on spectral editing. Relative errors in metabolite quantification were calculated with and without frequency matching of basis functions. The algorithm for numerical computation of basis functions was also improved for higher computational efficiency.

RESULTS

We found significant errors without frequency matching of basis functions when carrier frequency mismatch was generally considered negligible. By matching basis functions with the history of frequency deviation, the mean errors in glutamate, glutamine, γ-aminobutyric acid, and glutathione concentrations were reduced from 3.90%, 1.85%, 11.53%, and 3.43% to 0.18%, 0.34%, 0.40%, and 0.51%, respectively.

CONCLUSION

Matching basis functions to frequency deviation history was necessary even when frequency deviations during frequency-selective spectral editing were fairly small. Basis set frequency matching significantly improved accuracy in the quantification of glutamate, glutamine, γ-aminobutyric acid, and glutathione concentrations.

Prospective frequency correction using outer volume suppression-localized navigator for MR spectroscopy and spectroscopic imaging

PURPOSE

New frequency correction methods are required to achieve the accurate measurement of frequency drifts in MRS and MRSI. We present a prospective frequency correction method with outer volume suppression (OVS)-based localization and selective water excitation for effective frequency correction with better SNR improvement compared to other techniques.

METHODS

An OVS-localized navigator was developed to prospectively correct frequency drifts during MRS and MRSI measurements. The performance of the navigator was tested on the human brain and a solution phantom for frequency drifts induced by head motion or gradient heating by a preceding DWI experiment at 3T.

RESULTS

The OVS-localized navigator could accurately track motion-induced frequency drifts with an RMS error of 0.5 Hz. The SNR of MRS signals was not affected by use of the OVS-localized navigator when compared with and without the navigator (P > 0.05). The frequency drifts induced by DWI experiments were 5.1 ± 0.3 Hz/min during MRSI measurements on humans, resulting in increased spectral linewidth, significant bias in metabolite concentrations, and significantly increased Cramér-Rao lower bounds (P < 0.05). After prospective frequency corrections, the quality of MRSI was recovered to the level of those without any DWI-induced frequency drifts, judged by the spectral linewidth, metabolite concentrations, and Cramér-Rao lower bounds.

CONCLUSION

The OVS-localized navigator demonstrated effective prospective frequency corrections for large frequency drifts (5 Hz/min) without introducing any saturation-induced SNR loss. These benefits can be particularly beneficial for the acquisition of MRS signals with long T₁ and/or short TR, and spectral editing.

Real-time Correction of Motion and Imager Instability Artifacts during 3D γ-Aminobutyric Acid-edited MR Spectroscopic Imaging

Purpose To compare the involuntary head motion, frequency and B₀ shim changes, and effects on data quality during real-time-corrected three-dimensional γ-aminobutyric acid-edited magnetic resonance (MR) spectroscopic imaging in subjects with mild cognitive impairment (MCI), patients with Parkinson disease (PD), and young and older healthy volunteers. Materials and Methods In this prospective study, MR spectroscopic imaging datasets were acquired at 3 T after written informed consent was obtained. Translational and rotational head movement, frequency, and B₀ shim were determined with an integrated volumetric navigator. Head motion patterns and imager instability were investigated in 33 young healthy control subjects (mean age ± standard deviation, 31 years ± 5), 34 older healthy control subjects (mean age, 67 years ± 8), 34 subjects with MCI (mean age, 72 years ± 5), and 44 patients with PD (mean age, 64 years ± 8). Spectral quality was assessed by means of region-of-interest analysis. Group differences were evaluated with Bonferroni-corrected Mann-Whitney tests. Results Three patients with PD and four subjects with MCI were excluded because of excessive head motion (ie, > 0.8 mm translation per repetition time of 1.6 seconds throughout >10 minutes). Older control subjects, patients with PD, and subjects with MCI demonstrated 1.5, 2, and 2.5 times stronger head movement, respectively, than did young control subjects (1.79 mm ± 0.77) (P < .001). Of young control subjects, older control subjects, patients with PD, and subjects with MCI, 6%, 35%, 38%, and 51%, respectively, moved more than 3 mm during the MR spectroscopic imaging acquisition of approximately 20 minutes. The predominant movements were head nodding and "sliding out" of the imager. Frequency changes were 1.1- and 1.4-fold higher in patients with PD (P = .007) and subjects with MCI (P < .001), respectively, and B₀ shim changes were 1.3-, 1.5-, and 1.9-fold higher in older control subjects (P = .005), patients with PD (P < .001), and patients with MCI (P < .001), respectively, compared with those of young control subjects (12.59 Hz ± 2.49, 3.61 Hz · cm^-1 ± 1.25). Real-time correction provided high spectral quality in all four groups (signal-to-noise ratio >15, Cramér-Rao lower bounds < 20%). Conclusion Real-time motion and B₀ monitoring provides valuable information about motion patterns and B₀ field variations in subjects with different predispositions for head movement. Immediate correction improves data quality, particularly in patients who have difficulty avoiding movement. ^© RSNA, 2017 Online supplemental material is available for this article.

Retrospective correction of frequency drift in spectral editing: The GABA editing example

GABA levels can be measured using proton MRS with a two-step editing sequence. However due to the low concentration of GABA, long acquisition time is usually needed to achieve sufficient SNR to detect small differences in many psychiatric disorders. During this long scan time the frequency offset of the measured voxel can change because of magnetic field drift and patient movement. This drift will change the frequency of the editing pulse relative to that of metabolites, leading to errors in quantification. In this article we describe a retrospective method to correct for frequency drift in spectral editing. A series of reference signals for each metabolite was generated for a range of frequency offsets and then averaged together based on the history of frequency changes over the scan. These customized basis sets were used to fit the in vivo data. Our results demonstrate the effectiveness of the correction method and the remarkable robustness of a GABA editing technique with a top hat editing profile in the presence of frequency drift.

Measurement of Hypothalamic Glucose Under Euglycemia and Hyperglycemia by MRI at 3T

PURPOSE

To evaluate the feasibility of using a clinical magnetic resonance (MR) system and MR spectroscopy (MRS) to measure glucose concentration changes in the human hypothalamus, a structure central to whole-body glucose regulation.

SUBJECTS AND METHODS

A time series of MR spectra (semi-LASER, TE = 28 msec), localized to the bilateral hypothalamus (∼1.6 ml) were obtained at 3T in six healthy subjects at baseline (euglycemia) and during a ∼65-70-minute-long hyperglycemic clamp in 11-minute blocks with interleaved T1 FLASH images to retrospectively assess head motion, and track changes in cerebrospinal fluid (CSF) partial volume. The LCModel was used to quantify the sum of glucose and taurine concentrations, [Glc+Tau], along with their associated Cramér-Rao lower bounds (CRLB).

RESULTS

Spectral quality allowed quantification of [Glc+Tau] (sum reported due to high negative correlation between these metabolites) with CRLB <25% in 35/36 timepoints during hyperglycemia. Increased [Glc+Tau] was observed with hyperglycemia in all subjects, but most reliably in those with plasma glucose targets ≥300 mg/dl. For these subjects, [Glc+Tau]baseline (n = 4) was 1.5 (±0.3, SD) mM, and increased to 4.5 (±1.1) mM (n = 16) for timepoints acquired ≥25 minutes after onset of the clamp, with 15/16 timepoints having no overlap of 95% confidence intervals (CIs) between baseline and hyperglycemia. Preliminary analysis revealed a linear (1:5) relationship between hypothalamus-blood glucose concentrations.

CONCLUSION

It is feasible to measure glucose concentration changes in the human hypothalamus using a standard 3T scanner and a short-echo semi-LASER sequence by utilizing retrospective motion tracking, CSF correction, predetermined quality acceptance criteria, and hyperglycemic blood glucose levels ≥300 mg/dl.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:681-691.

Prospective frequency correction for macromolecule-suppressed GABA editing at 3T

PURPOSE

To investigate the effects of B₀ field offsets and drift on macromolecule (MM)-suppressed GABA-editing experiments, and to implement and test a prospective correction scheme. "Symmetric" editing schemes are proposed to suppress unwanted coedited MM signals in GABA editing.

MATERIALS AND METHODS

Full density-matrix simulations of both conventional (nonsymmetric) and symmetric MM-suppressed editing schemes were performed for the GABA spin system to evaluate their offset-dependence. Phantom and in vivo (15 subjects at 3T) GABA-edited experiments with symmetrical suppression of MM signals were performed to quantify the effects of field offsets on the total GABA+MM signal (designated GABA+). A prospective frequency correction method based on interleaved water referencing (IWR) acquisitions was implemented and its experimental performance evaluated during positive and negative drift.

RESULTS

Simulations show that the signal from MM-suppressed symmetrical editing schemes is an order of magnitude more susceptible to field offsets than the signal from nonsymmetric editing schemes. The MM-suppressed GABA signal changes by 8.6% per Hz for small field offsets. IWR significantly reduces variance in the field offset and measured GABA levels (both P < 0.001 by F-tests), maintaining symmetric suppression of MM signal.

CONCLUSION

Symmetrical editing schemes substantially increase the dependence of measurements on B₀ field offsets, which can arise due to patient movement and/or scanner instability. It is recommended that symmetrical editing should be used in combination with effective B₀ stabilization, such as that provided by IWR. J. Magn. Reson. Imaging 2016;44:1474-1482.

Basic Principles and Clinical Applications of Magnetic Resonance Spectroscopy in Neuroradiology

Magnetic resonance spectroscopy is a powerful tool to assist daily clinical diagnostics. This review is intended to give an overview on basic principles of the technology, discuss some of its technical aspects, and present typical applications in daily clinical routine in neuroradiology.

Volumetric navigated MEGA-SPECIAL for real-time motion and shim corrected GABA editing

Mescher-Garwood (MEGA) editing with spin echo full intensity acquired localization (MEGA-SPECIAL, MSpc) is a technique to acquire γ-aminobutyric acid (GABA) without macromolecule (MM) contamination at a TE of 68 ms. However, due to the requirement of multiple shot-localization, it is often susceptible to subject motion and B0 inhomogeneity. A method is presented for real-time shim and motion correction (ShMoCo) using volumetric navigators to correct for motion and motion-related B0 inhomogeneity during MSpc acquisition. A phantom experiment demonstrates that ShMoCo restores the GABA peak and improves spectral quality in the presence of motion and zero- and first-order shim changes. The ShMoCo scans were validated in three subjects who performed up-down and left-right head rotations. Qualitative assessment of these scans indicates effective reduction of subtraction artefacts and well edited GABA peaks, while quantitative analysis indicates superior fitting and spectral quality relative to scans with no correction. Copyright © 2015 John Wiley & Sons, Ltd.

Knee cartilage MRI with in situ mechanical loading using prospective motion correction

PURPOSE

To assess the feasibility of high resolution knee cartilage MRI with in situ mechanical loading using optical tracking to compensate for motion.

METHODS

In vivo cartilage MRI with in situ mechanical loading is demonstrated on a clinical 3T system for the patellofemoral as well as for the tibiofemoral knee joint using a T1-weighted spoiled three-dimensional gradient-echo sequence. Prospective motion correction is performed with a moiré phase tracking system consisting of an in-bore camera and a single tracking marker attached to the skin.

RESULTS

Rigid-body approximation required for prospective correction with optical motion tracking is fulfilled well enough for the patellofemoral as well as for the tibiofemoral joint when the tracking marker is attached to the knee cap and the shin, respectively. Presaturation proves to be efficient in suppressing pulsation artifacts from the popliteal artery and residual motion artifacts primarily arising from nonrigid motion of the posterior knee compartment.

CONCLUSION

The proposed technique enables knee cartilage imaging under in situ mechanical loading with submillimeter spatial resolution devoid of significant motion artifacts and thus appropriate for cartilage volumetry. It has the potential to provide new insight into the biomechanics of the knee and might complement the panoply of diagnostic MR methods for osteoarthritis.

Reverse retrospective motion correction

PURPOSE

One potential barrier for using prospective motion correction (PMC) in the clinic is the unpredictable nature of a scan because of the direct interference with the imaging sequence. We demonstrate that a second set of "de-corrected" images can be reconstructed from a scan with PMC that show how images would have appeared without PMC enabled.

THEORY AND METHODS

For three-dimensional scans, the effects of PMC can be undone by performing a retrospective reconstruction based on the inverse of the transformation matrix used for real time gradient feedback. Retrospective reconstruction is performed using a generalized SENSE approach with continuous head motion monitored using a single-marker optical camera system.

RESULTS

Reverse retrospective reconstruction is demonstrated for phantom and in vivo scans using an magnetization-prepared rapid gradient echo (MPRAGE) sequence including parallel and Partial Fourier acceleration.

CONCLUSION

Reverse retrospective reconstruction can almost perfectly undo the effects of prospective feedback, and thereby provide a second image data set with the effects of motion correction removed. In case of correct feedback, one can directly compare the quality of the corrected with that of the uncorrected scan. Additionally, because erroneous feedback during PMC may introduce artifacts, it is possible to eliminate artifacts in a corrupted scan by reversing the false gradient updates. Magn Reson Med 75:2341-2349, 2016. © 2015 Wiley Periodicals, Inc.

Motion artifacts in MRI: A complex problem with many partial solutions

Subject motion during magnetic resonance imaging (MRI) has been problematic since its introduction as a clinical imaging modality. While sensitivity to particle motion or blood flow can be used to provide useful image contrast, bulk motion presents a considerable problem in the majority of clinical applications. It is one of the most frequent sources of artifacts. Over 30 years of research have produced numerous methods to mitigate or correct for motion artifacts, but no single method can be applied in all imaging situations. Instead, a "toolbox" of methods exists, where each tool is suitable for some tasks, but not for others. This article reviews the origins of motion artifacts and presents current mitigation and correction methods. In some imaging situations, the currently available motion correction tools are highly effective; in other cases, appropriate tools still need to be developed. It seems likely that this multifaceted approach will be what eventually solves the motion sensitivity problem in MRI, rather than a single solution that is effective in all situations. This review places a strong emphasis on explaining the physics behind the occurrence of such artifacts, with the aim of aiding artifact detection and mitigation in particular clinical situations.

3D GABA imaging with real-time motion correction, shim update and reacquisition of adiabatic spiral MRSI

Gamma-aminobutyric acid (GABA) and glutamate (Glu) are the major neurotransmitters in the brain. They are crucial for the functioning of healthy brain and their alteration is a major mechanism in the pathophysiology of many neuro-psychiatric disorders. Magnetic resonance spectroscopy (MRS) is the only way to measure GABA and Glu non-invasively in vivo. GABA detection is particularly challenging and requires special MRS techniques. The most popular is MEscher-GArwood (MEGA) difference editing with single-voxel Point RESolved Spectroscopy (PRESS) localization. This technique has three major limitations: a) MEGA editing is a subtraction technique, hence is very sensitive to scanner instabilities and motion artifacts. b) PRESS is prone to localization errors at high fields (≥3T) that compromise accurate quantification. c) Single-voxel spectroscopy can (similar to a biopsy) only probe steady GABA and Glu levels in a single location at a time. To mitigate these problems, we implemented a 3D MEGA-editing MRS imaging sequence with the following three features: a) Real-time motion correction, dynamic shim updates, and selective reacquisition to eliminate subtraction artifacts due to scanner instabilities and subject motion. b) Localization by Adiabatic SElective Refocusing (LASER) to improve the localization accuracy and signal-to-noise ratio. c) K-space encoding via a weighted stack of spirals provides 3D metabolic mapping with flexible scan times. Simulations, phantom and in vivo experiments prove that our MEGA-LASER sequence enables 3D mapping of GABA+ and Glx (Glutamate+Gluatmine), by providing 1.66 times larger signal for the 3.02ppm multiplet of GABA+ compared to MEGA-PRESS, leading to clinically feasible scan times for 3D brain imaging. Hence, our sequence allows accurate and robust 3D-mapping of brain GABA+ and Glx levels to be performed at clinical 3T MR scanners for use in neuroscience and clinical applications.

Automatic voxel positioning for MRS at 7 T

OBJECT

The purpose of this study was to test, for the first time, whether spectroscopy voxels could be positioned automatically with high accuracy and reproducibility in ultrahigh-field longitudinal magnetic resonance spectroscopy (MRS) studies.

MATERIALS AND METHODS

MRS voxels were automatically positioned in two cingulate subregions of 12 healthy subjects using a vendor-provided automatic voxel positioning (AutoAlign) technique, and were manually placed in the same regions of 10 healthy subjects by an experienced technician in three 7 T MRS scan sessions. Different coils were used for manual (24-channel coil) and automatic (32-channel coil) voxel placement, and the effects of signal-to-noise-ratio differences on the spectra were considered.

RESULTS

Over three scan sessions and two regions scanned for each subject, a mean voxel geometric overlap ratio of 0.91 for automatic positioning reflected accurate voxel alignment, while the geometric overlap ratio was only 0.70 for voxels placed manually. Comparable voxel positions among the three scan sessions (p > 0.05) indicated high reproducibility of automatic voxel alignment. In comparison, significant voxel displacement among scan sessions (p < 0.05) was found using manual voxel positioning.

CONCLUSIONS

In view of the highly accurate and reproducible voxel alignment with automatic voxel positioning, we propose the application of automatic rather than manual voxel positioning in future ultrahigh-field longitudinal MRS studies.

The use of magnetic resonance spectroscopy as a tool for the measurement of bi-hemispheric transcranial electric stimulation effects on primary motor cortex metabolism

Transcranial direct current stimulation (tDCS) is a neuromodulation technique that has been increasingly used over the past decade in the treatment of neurological and psychiatric disorders such as stroke and depression. Yet, the mechanisms underlying its ability to modulate brain excitability to improve clinical symptoms remains poorly understood. To help improve this understanding, proton magnetic resonance spectroscopy ((1)H-MRS) can be used as it allows the in vivo quantification of brain metabolites such as γ-aminobutyric acid (GABA) and glutamate in a region-specific manner. In fact, a recent study demonstrated that (1)H-MRS is indeed a powerful means to better understand the effects of tDCS on neurotransmitter concentration. This article aims to describe the complete protocol for combining tDCS (NeuroConn MR compatible stimulator) with (1)H-MRS at 3 T using a MEGA-PRESS sequence. We will describe the impact of a protocol that has shown great promise for the treatment of motor dysfunctions after stroke, which consists of bilateral stimulation of primary motor cortices. Methodological factors to consider and possible modifications to the protocol are also discussed.

A comparison of spectral quality in magnetic resonance spectroscopy data acquired with and without a novel EPI-navigated PRESS sequence in school-aged children with fetal alcohol spectrum disorders

Single voxel spectroscopy (SVS) can generate useful information regarding metabolite concentrations provided that the MR signal can be averaged over several minutes during which the subject remains stationary. This requirement can be particularly challenging for children who cannot otherwise be scanned without sedation. To address this problem we developed an EPI volume navigated (vNav) SVS PRESS sequence, which applies real-time head pose (location and orientation), frequency, and first-order B0 shim adjustments. A water-independent preprocessing algorithm removes residual frequency and phase shifts resulting from within-TR movements. We compare results and performance of the standard and vNav PRESS sequences in a sample of 9- to 10-year-olds from a South African cohort of children with fetal alcohol spectrum disorders (FASD) and healthy controls. Magnetic resonance spectroscopy (MRS) data in the deep cerebellar nuclei were initially acquired with the standard PRESS sequence. The children were re-scanned 1 year later with the vNav PRESS sequence. Good quality data were acquired in 73% using the vNav PRESS sequence, compared to only 50% for the standard PRESS sequence. Additionally, tighter linewidths and smaller variances in the measured concentrations were observed. These findings confirm previous reports demonstrating the efficacy of our innovative vNav sequence with healthy volunteers and young children with HIV and expand its application to a school-aged population with FASD-disorders often associated with attention problems and hyperactivity. This study provides the most direct evidence to date regarding degree to which these new methods can improve data quality in research studies employing MRS.

Frequency and phase drift correction of magnetic resonance spectroscopy data by spectral registration in the time domain

PURPOSE

Frequency and phase drifts are a common problem in the acquisition of in vivo magnetic resonance spectroscopy (MRS) data. If not accounted for, frequency and phase drifts will result in artifactual broadening of spectral peaks, distortion of spectral lineshapes, and a reduction in signal-to-noise ratio (SNR). We present herein a new method for estimating and correcting frequency and phase drifts in in vivo MRS data.

METHODS

We used a simple method of fitting each spectral average to a reference scan (often the first average in the series) in the time domain through adjustment of frequency and phase terms. Due to the similarity with image registration, this method is referred to as "spectral registration." Using simulated data with known frequency and phase drifts, the performance of spectral registration was compared with two existing methods at various SNR levels.

RESULTS

Spectral registration performed well in comparison with the other methods tested in terms of both frequency and phase drift estimation.

CONCLUSIONS

Spectral registration provides an effective method for frequency and phase drift correction. It does not involve the collection of navigator echoes, and does not rely on any specific resonances, such as residual water or creatine, making it highly versatile.

Real-time motion- and B0-correction for LASER-localized spiral-accelerated 3D-MRSI of the brain at 3T

The full potential of magnetic resonance spectroscopic imaging (MRSI) is often limited by localization artifacts, motion-related artifacts, scanner instabilities, and long measurement times. Localized adiabatic selective refocusing (LASER) provides accurate B1-insensitive spatial excitation even at high magnetic fields. Spiral encoding accelerates MRSI acquisition, and thus, enables 3D-coverage without compromising spatial resolution. Real-time position- and shim/frequency-tracking using MR navigators correct motion- and scanner instability-related artifacts. Each of these three advanced MRI techniques provides superior MRSI data compared to commonly used methods. In this work, we integrated in a single pulse sequence these three promising approaches. Real-time correction of motion, shim, and frequency-drifts using volumetric dual-contrast echo planar imaging-based navigators were implemented in an MRSI sequence that uses low-power gradient modulated short-echo time LASER localization and time efficient spiral readouts, in order to provide fast and robust 3D-MRSI in the human brain at 3T. The proposed sequence was demonstrated to be insensitive to motion- and scanner drift-related degradations of MRSI data in both phantoms and volunteers. Motion and scanner drift artifacts were eliminated and excellent spectral quality was recovered in the presence of strong movement. Our results confirm the expected benefits of combining a spiral 3D-LASER-MRSI sequence with real-time correction. The new sequence provides accurate, fast, and robust 3D metabolic imaging of the human brain at 3T. This will further facilitate the use of 3D-MRSI for neuroscience and clinical applications.

Quality of 186 child brain spectra using motion and B0 shim navigated single voxel spectroscopy

PURPOSE

To evaluate B0 shim and motion navigated single voxel spectroscopy in children. Assess the repeatability of metabolite concentrations in three regions: medial frontal grey matter, peritrigonal white matter, and basal ganglia. Determine the extent of intra- and interacquisition movement in this population.

METHODS

Linewidth and signal to noise ratio were calculated to assess spectral quality of 186 spectra at 3 Tesla. Repeatability was assessed on 31 repeat scans. Navigator images were used to assess localization errors, while navigator motion and shim logs were used to demonstrate the efficacy of correction needed during the scans.

RESULTS

Average linewidths ± standard deviations of N-acetyl aspartate are 3.8 ± 0.6 Hz, 4.4 ± 0.5 Hz, and 4.7 ± 0.8 Hz in each region, respectively. Scan-to-scan measurement variance in metabolite concentrations closely resembled the expected variance. A total of 73% and 32% of children moved before and during the acquisition, causing a voxel shift of more than 10% of the voxel volume, 1.5 mm. The predominant movement directions were sliding out of the coil and nodding (up-down rotation). First-order B0 corrections were significant (>10 μT/m) in 18 % of acquisitions.

CONCLUSION

Prospective motion and B0 correction provides high quality repeatable spectra. The study found that most children moved between acquisitions and a substantial number moved during acquisitions.

Real-time automated spectral assessment of the BOLD response for neurofeedback at 3 and 7T

Echo-planar imaging is the dominant functional MRI data acquisition scheme for evaluating the BOLD signal. To date, it remains the only approach providing neurofeedback from spatially localized brain activity. Real-time functional single-voxel proton spectroscopy (fSVPS) may be an alternative for spatially specific BOLD neurofeedback at 7T because it allows for a precise estimation of the local T2* signal, EPI-specific artifacts may be avoided, and the signal contrast may increase. In order to explore and optimize this alternative neurofeedback approach, we tested fully automated real-time fSVPS spectral estimation procedures to approximate T2* BOLD signal changes from the unsuppressed water peak, i.e. lorentzian non-linear complex spectral fit (LNLCSF) in frequency and frequency-time domain. The proposed approaches do not require additional spectroscopic localizers in contrast to conventional T2* approximation based on linear regression of the free induction decay (FID). For methods comparison, we evaluated quality measures for signals from the motor and the visual cortex as well as a real-time feedback condition at high (3T) and at ultra-high (7T) magnetic field strengths. Using these methods, we achieved reliable and fast water peak spectral parameter estimations. At 7T, we observed an absolute increase of spectra line narrowing due to the BOLD effect, but quality measures did not improve due to artifactual line broadening. Overall, the automated fSVPS approach can be used to assess dynamic spectral changes in real-time, and to provide localized T2* neurofeedback at 3 and 7T.

Localization errors in MR spectroscopic imaging due to the drift of the main magnetic field and their correction

PURPOSE

To analyze the effect of B0 field drift on multivoxel MR spectroscopic imaging and to propose an approach for its correction.

THEORY AND METHODS

It is shown, both theoretically and in a phantom, that for ∼30 min acquisitions a linear B0 drift (∼0.1 ppm/h) will cause localization errors that can reach several voxels (centimeters) in the slower varying phase encoding directions. An efficient and unbiased estimator is proposed for tracking the drift by interleaving short (∼ T2*), nonlocalized acquisitions on the nonsuppressed water each pulse repetition time, as shown in 10 volunteers at 1.5 and 3 T.

RESULTS

The drift is shown to be predominantly linear in both the phantom and volunteers at both fields. The localization errors are observed and quantified in both phantom and volunteers. The unbiased estimator is shown to reliably track the instantaneous frequency in vivo despite only using a small portion of the FID.

CONCLUSION

Contrary to single-voxel MR spectroscopy, where it leads to line broadening, field drift can lead to localization errors in the longer chemical shift imaging experiments. Fortunately, this drift can be obtained at a negligible cost to sequence timing, and corrected for in post processing.

Prospective motion correction in brain imaging: a review

Motion correction in magnetic resonance imaging by real-time adjustment of the imaging pulse sequence was first proposed more than 20 years ago. Recent advances have resulted from combining real-time correction with new navigator and external tracking mechanisms capable of quantifying rigid-body motion in all 6 degrees of freedom. The technique is now often referred to as "prospective motion correction." This article describes the fundamentals of prospective motion correction and reviews the latest developments in its application to brain imaging and spectroscopy. Although emphasis is placed on the brain as the organ of interest, the same principles apply whenever the imaged object can be approximated as a rigid body. Prospective motion correction can be used with most MR sequences, so it has potential to make a large impact in clinical routine. To maximize the benefits obtained from the technique, there are, however, several challenges still to be met. These include practical implementation issues, such as obtaining tracking data with minimal delay, and more fundamental problems, such as the magnetic field distortions caused by a moving object. This review discusses these challenges and summarizes the state of the art. We hope that this work will motivate further developments in prospective motion correction and help the technique to reach its full potential.

Parallel and partial Fourier imaging with prospective motion correction

Subject motion during scan is a major source of artifacts in MR examinations. Prospective motion correction is a promising technique that tracks subject motion and adjusts the imaging volume in real time; however, additional retrospective correction may be necessary to achieve robust image quality and compatibility with other imaging options. Real-time realignment of the imaging volume by prospective motion correction changes the coil sensitivity weighting and the field inhomogeneity relative to the imaging volume. This can pose image reconstruction problems with parallel imaging and partial Fourier imaging, which rely on coil sensitivity and image phase information, respectively. This work presents a practical method for reconstructing images acquired using prospective motion correction with parallel imaging and/or partial Fourier imaging. Our proposed approach is data driven and noniterative; data are binned into several position bins based on motion measurements made during the prospective motion correction acquisition and the data in each bin are processed through intrabin operations such as parallel imaging reconstruction (in case of undersampling), phase correction, and coil combination before combination of the position bins. We demonstrate the effectiveness of our technique through simulation studies and in vivo experiments using a prospectively motion-corrected three-dimensional fast spin echo sequence.

Real-time motion and B0 correction for localized adiabatic selective refocusing (LASER) MRSI using echo planar imaging volumetric navigators

A method is presented to correct the effects of motion and motion-related B(0) perturbations on spectroscopic imaging in real time through the use of a volumetric navigator. It is demonstrated that, for an axial slice, lifting the chin significantly disrupts the B(0) homogeneity in the zero-order (frequency), first-order Y (coronal) axis and second-order ZY term. This volumetric navigator is able to measure and correct in real time both head pose and zero- to first-order B(0) inhomogeneities. The volumetric navigator was validated in six volunteers who deliberately lifted and then dropped their chin during the scan. These scans show that motion correction alone is not sufficient to recover the spectral quality. By applying real-time shim adjustments, spectral quality was fully recovered to linewidths below 0.08 ppm and the signal-to-noise ratio to within acceptable limits in five of six subjects. In the sixth subject, 83% of the spectra within the volume of interest were recovered, compared with the worst case nonshim-corrected scan, where none of the voxels fell within these quality bounds. It is shown that the use of a volumetric navigator comes at no additional cost to the scan time or spectral signal-to-noise ratio.

Real-time dynamic frequency and shim correction for single-voxel magnetic resonance spectroscopy

Subject motion during brain magnetic resonance spectroscopy acquisitions generally reduces the magnetic field (B₀) homogeneity across the volume of interest or voxel. This is the case even if prospective motion correction ensures that the voxel follows the head. We introduce a novel method for rapidly mapping linear variations in B₀ across a small volume using two-dimensional excitations. The new field mapping technique was integrated into a prospectively motion-corrected single-voxel ¹H magnetic resonance spectroscopy sequence. Interference with the magnetic resonance spectroscopy measurement was negligible, and there was no penalty in scan time. Frequency shifts were also measured continuously, and both frequency and first-order shim corrections were applied in real time. Phantom experiments and in vivo studies demonstrated that the resulting motion- and shim-corrected sequence is able to mitigate line broadening and maintain spectral quality even in the presence of large-amplitude subject motion.

Single Voxel Proton Spectroscopy for Neurofeedback at 7 Tesla

Echo-planar imaging (EPI) in fMRI is regularly used to reveal BOLD activation in presubscribed regions of interest (ROI). The response is mediated by relative changes in T2* which appear as changes in the image pixel intensities. We have proposed an application of functional single-voxel proton spectroscopy (fSVPS) for real-time studies at ultra-high MR field which can be comparable to the EPI BOLD fMRI technique. A spin-echo SVPS protocol without water suppression was acquired with 310 repetitions on a 7T Siemens MR scanner (TE/TR = 20/1000 ms, flip angle α = 90°, voxel size 10 × 10 × 10 mm³). Transmitter reference voltage was optimized for the voxel location. Spectral processing of the water signal free induction decay (FID) using log-linear regression was used to estimate the T2* change between rest and activation of a functional task. The FID spectrum was filtered with a Gaussian window around the water peak, and log-linear regression was optimized for the particular ROI by adoption of the linearization length. The spectroscopic voxel was positioned on an ROI defined from a real-time fMRI EPI BOLD localizer. Additional online signal processing algorithms performed signal drift removal (exponential moving average), despiking and low-pass filtering (modified Kalman filter) and, finally, the dynamic feedback signal normalization. Two functional tasks were used to estimate the sensitivity of the SVPS method compared to BOLD signal changes, namely the primary motor cortex (PMC, left hand finger tapping) and visual cortex (VC, blinking checkerboard). Four healthy volunteers performed these tasks and an additional session using real-time signal feedback modulating their activation level of the PMC. Results show that single voxel spectroscopy is able to provide a good and reliable estimation of the ΒΟLD signal changes. Small data size and FID signal processing instead of processing entire brain volumes as well as more information revealed from the acquired total water spectrum, i.e., direct estimation of the T2* values and B0 changes, make SVPS proton spectroscopy suitable and advantageous for real-time neurofeedback studies. Particular challenges of ultra-high field spectroscopy due to the non-linearity in the spectral information, e.g., poor main magnetic field homogeneity and the absence of motion correction for the SVPS sequence may lead to the special artifacts in the control signal which still need to be addressed. The contrast to noise ratio (CNR), experimental statistic (t-values) and percent signal change were used as quality parameters to estimate the method performance. The potential and challenges of the spectroscopic approach for fMRI studies needs to be further investigated.

Spectroscopic imaging with prospective motion correction and retrospective phase correction

Motion-induced artifacts are much harder to recognize in magnetic resonance spectroscopic imaging than in imaging experiments and can therefore lead to erroneous interpretation. A method for prospective motion correction based on an optical tracking system has recently been proposed and has already been successfully applied to single voxel spectroscopy. In this work, the utility of prospective motion correction in combination with retrospective phase correction is evaluated for spectroscopic imaging in the human brain. Retrospective phase correction, based on the interleaved reference scan method, is used to correct for motion-induced frequency shifts and ensure correct phasing of the spectra across the whole spectroscopic imaging slice. It is demonstrated that the presented correction methodology can reduce motion-induced degradation of spectroscopic imaging data.

Correction of frequency drifts induced by gradient heating in 1H spectra using interleaved reference spectroscopy

PURPOSE

To evaluate the impact of heating-induced frequency drifts on single-voxel spectroscopy and to demonstrate correction strategies based on the interleaved reference scan technique.

MATERIALS AND METHODS

Frequency drifts induced by gradient heating are assessed for two clinical 3 Tesla (T) whole body MR systems. The interleaved reference scan (IRS) method is used for correcting these frequency drifts in 1H spectra in vitro and in vivo. For severely drift-affected spectroscopy experiments, a feedback-based version of the IRS sequence is proposed, which adds the functionality of a frequency lock to prevent a degradation of the water suppression.

RESULTS

It is shown that the line widths of the spectral resonances can be largely reduced with the interleaved reference scan method, resulting in considerably improved peak resolution and signal-to-noise ratio. The feedback-based IRS method additionally allows for stable water suppression, even in the presence of very strong frequency drifts.

CONCLUSION

If spectroscopy scans are combined with imaging scans with a high gradient duty cycle such as diffusion-weighted imaging or functional MRI, a drift correction with IRS can considerably improve the validity of data analysis in research studies.

Prospective motion correction for magnetic resonance spectroscopy using single camera Retro-Grate reflector optical tracking

PURPOSE

To introduce and evaluate a method of prospective motion correction for localized proton magnetic resonance spectroscopy (1H-MRS) using a single-camera optical tracking system.

MATERIALS AND METHODS

Five healthy participants were scanned at 3T using a point-resolved spectroscopic sequence (PRESS) with a motion-tracking module and phase navigator. Head motion in six degrees was tracked with a Retro-Grate Reflector (RGR) tracking system and target via a mirror mounted inside the bore. Participants performed a series of three predetermined motion patterns during scanning.

RESULTS

Left-right rotation (Rz) (average 12°) resulted in an increase in the total choline to total creatine ratio (Cho/Cr) of +14.6 ± 1.5% (P = 0.0009) for scans without correction, but no change for scans with correction (+1.1 ± 1.5%; P = 0.76). Spectra with uncorrected Z-translations showed large lipid peaks (skull) with changes in Cho/Cr of -13.2 ± 1.6% (P = 0.02, no motion correction) and -2.2 ± 2.4% (P = 0.51) with correction enabled. There were no significant changes in the ratios of N-acetylaspartate, glutamate+glutamine, or myo-inositol to creatine compared to baseline scans for all experiments.

CONCLUSION

Prospective motion correction for 1H-MRS, using single-camera RGR tracking, can reduce spectral artifacts and quantitation errors in Cho/Cr ratios due to head motion and promises improved spectral quality and reproducibility.

Real-time motion and B0 corrected single voxel spectroscopy using volumetric navigators

In population groups where head pose cannot be assumed to be constant during a magnetic resonance spectroscopy examination or in difficult-to-shim regions of the brain, real-time volume of interest, frequency, and shim optimization may be necessary. We investigate the effect of pose change on the B0 homogeneity of a (2 cm)3 volume and observe typical first-order shim changes of 1 μT/m per 1° rotation (chin down to up) in four different volumes of interest in a single volunteer. An echo planar imaging volume navigator was constructed to measure and apply in real-time within each pulse repetition time: volume of interest positioning, frequency adjustment, and first-order shim adjustment. This volume navigator is demonstrated in six healthy volunteers and achieved a mean linewidth of 4.4 Hz, similar to that obtained by manual shim adjustment of 4.9 Hz. Furthermore, this linewidth is maintained by the volume navigator at 4.9 Hz in the presence of pose change. By comparison, a mean linewidth of 7.5 Hz was observed, when no correction was applied.

A novel phase and frequency navigator for proton magnetic resonance spectroscopy using water-suppression cycling

Magnetic resonance spectroscopy is sensitive to movements, in part, because of motion-induced phase and frequency variations that lead to incoherent averaging. For in vivo proton magnetic resonance spectroscopy, the unsuppressed or under-suppressed water signal can be used to restore coherent averaging; however, this approach results in baseline distortions due to the large water peak. Therefore, a novel water-suppression cycling scheme was developed that alternates between positive and negative residual water signal. Using the residual water signal, the method allows for shot-to-shot phase and frequency correction of individual free induction decays and restoration of signal losses due to incoherent averaging, yet near-complete elimination of residual water. It is demonstrated that the residual water signal can be used to restore metabolite peaks in a brain spectrum from a subject who performed intentional head movements. The ability to correct phase and frequency fluctuations during subject motion is vital for use with adaptive motion correction approaches that ensure proper voxel positioning during head movements.

Rapid six-degree-of-freedom motion detection using prerotated baseline spherical navigator echoes

A new spherical navigator echo (SNAV) registration technique is presented. This technique starts by collecting a set of SNAV templates at a reference position. These templates are acquired by rotating the gradient system to result in rotation angles that uniformly cover a predefined range of rotation. The rotation angles between an unknown physically transformed position and the reference position are subsequently determined by finding the template with the lowest sum of squared differences with SNAV at the transformed position. Translations are calculated from the phase differences between the best-match SNAV template and the SNAV acquired at the transformed position. In comparison with the conventional SNAV registration technique, the proposed technique is noniterative, robust, and can detect 3-dimensional rigid body motion in less than 50 msec. The technique was verified with phantom and in vivo experiments, which demonstrated subdegree rotational and submillimeter translational accuracy over a range of simultaneous ±20° and ±10° mm of motion.