Accelerated proton echo planar spectroscopic imaging (PEPSI) using GRAPPA with a 32-channel phased-array coil

Downfield proton MRSI with whole-brain coverage at 3T

PURPOSE

To develop a 3D downfield (DF) MRSI protocol with whole brain coverage and post-processing pipeline for creation of metabolite maps.

METHODS

A 3D, circularly phase-encoded version of the previously developed 2D DF MRSI sequence with1 3 ‾ 3 1 ‾ $$ 1\overline{3}3\overline{1} $$ spectral-spatial excitation and frequency selective refocusing was implemented and tested in five healthy volunteers at 3T. The DF metabolite maps with a nominal spatial resolution of 0.7 cm3 were recorded in eight slices at 3T in a scan time of 22 m 40 s. An MRSI post-processing pipeline was developed to create DF metabolite maps. Metabolite concentrations and uncertainty estimates were compared between region differences for nine DF peaks.

RESULTS

LCModel analysis showed Cramer Rao lower bounds average values of 3%-4% for protein amide resonances in the three selected regions (anterior cingulate, dorsolateral prefrontal cortex, and centrum semiovale); Cramer Rao lower bounds were somewhat higher for individual peaks but for the most part were less than 20%. While DF concentration maps were visually quite homogeneous throughout the brain, general linear regression analysis corrected for multiple comparisons found significant differences between centrum semiovale and dorsolateral prefrontal cortex for peaks at 7.09 ppm (p = 0.014), 7.90 ppm (p = 0.009), 8.18 ppm (p = 0.009), combined amides (p = 0.009), and between anterior cingulate and dorsolateral prefrontal cortex for the 7.30 ppm peak (p = 0.020). Cramer Rao lower bounds values were not significantly different between brain regions for any of the DF peaks.

CONCLUSION

The 3D DF MRSI of the human brain at 3T with wide spatial coverage for the mapping of exchangeable amide and other resonances is feasible at a nominal spatial resolution of 0.7 cm3 .

Downfield Proton MRSI with whole-brain coverage at 3T

Purpose

To develop a 3D downfield magnetic resonance spectroscopic imaging (DF-MRSI) protocol with whole brain coverage and post-processing pipeline for creation of metabolite maps.

Methods

A 3D, circularly phase-encoded version of the previously developed 2D DF-MRSI sequence with spectral-spatial excitation and frequency selective refocusing was implemented and tested in 5 healthy volunteers at 3T. Downfield metabolite maps with a nominal spatial resolution of 0.7 cm 3 were recorded in 8 slices at 3T in a scan time of 22m 40s. An MRSI post-processing pipeline was developed to create DF metabolite maps. Metabolite concentrations and uncertainty estimates were compared between region differences for nine downfield peaks.

Results

LCModel analysis showed CRLB average values of 3-4% for protein amide resonances in the three selected regions (anterior cingulate (ACC), dorsolateral prefrontal cortex (DLPFC), and centrum semiovale (CSO)); CRLBs were somewhat higher for individual peaks but for the most part were less than 20%. While DF concentration maps were visually quite homogeneous throughout the brain, general linear regression analysis corrected for multiple comparisons found significant differences between CSO and DLPFC for peaks at 7.09 ppm (p= 0.014), 7.90 ppm (p=0.009), 8.18 ppm (p=0.009), combined amides (p=0.009), and between ACC and DLPFC for the 7.30 ppm peak (p=0.020). CRLB values were not significantly different between brain regions for any of the DF peaks.

Conclusion

3D DF-MRSI of the human brain at 3T with wide spatial coverage for the mapping of exchangeable amide and other resonances is feasible at a nominal spatial resolution of 0.7 cm 3 .

Improved signal-to-noise performance of MultiNet GRAPPA 1 H FID MRSI reconstruction with semi-synthetic calibration data

PURPOSE

To further develop MultiNet GRAPPA, a neural-network-based reconstruction, for lower SNR proton MRSI (¹ H MRSI) data using adapted undersampling schemes and improved training sets.

METHODS

¹ H FID-MRSI data and an anatomical image for GRAPPA reconstruction were acquired in two slices in the human brain (n = 6) at 7T. MRSI data were retrospectively undersampled for a 4×, 6×, and 7× acceleration rate. Signal-to-noise, relative error (RE) between accelerated and fully sampled metabolic maps, RMS of the lipid artifacts, and fitting reliability were compared across acceleration rates, to the fully sampled data, and with different kinds and amounts of training images.

RESULTS

Training with semi-synthetic images resulted in higher SNR and lower lipid RMS relative to training with acquired images from one or several subjects. SNR increased with the number of semi-synthetic training images and the 4× accelerated data retains ∼30% more SNR than other accelerated data. Spectra reconstructed with 20 semi-synthetic averages retained ∼100% more SNR and had ∼5% lower lipid RMS than those reconstructed with the center k-space points of one image as was originally proposed for very high SNR MRSI data and had higher fitting reliability. The metabolite RE was lowest when training with 20-semi-synthetic training images and highest when training with the center k-space points of one image.

CONCLUSION

MultiNet GRAPPA is feasible with lower SNR ¹ H MRSI data if 20-semi-synthetic training images are used at a 4× acceleration rate. This acceleration rate provided the best trade-off between scan time and spectral SNR.

High-resolution, 3D multi-TE 1 H MRSI using fast spatiospectral encoding and subspace imaging

PURPOSE

To develop a novel method to achieve fast, high-resolution, 3D multi-TE ¹ H-MRSI of the brain.

METHODS

A new multi-TE MRSI acquisition strategy was developed that integrates slab selective excitation with adiabatic refocusing for better volume coverage, rapid spatiospectral encoding, sparse multi-TE sampling, and interleaved water navigators for field mapping and calibration. Special data processing strategies were developed to interpolate the sparsely sampled data, remove nuisance signals, and reconstruct multi-TE spatiospectral distributions with high SNR. Phantom and in vivo experiments have been carried out to demonstrate the capability of the proposed method.

RESULTS

The proposed acquisition can produce multi-TE ¹ H-MRSI data with three TEs at a nominal spatial resolution of 3.4 × 3.4 × 5.3 mm³ in around 20 min. High-SNR brain metabolite spatiospectral reconstructions can be obtained from both a metabolite phantom and in vivo experiments by the proposed method.

CONCLUSION

High-resolution, 3D multi-TE ¹ H-MRSI of the brain can be achieved within clinically feasible time. This capability, with further optimizations, could be translated to clinical applications and neuroscience studies where simultaneously mapping metabolites and neurotransmitters and TE-dependent molecular spectral changes are of interest.

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

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

Clinical Applications of Magnetic Resonance Spectroscopy in Brain Tumors: From Diagnosis to Treatment

Magnetic resonance spectroscopy (MRS) is a valuable tool for imaging brain tumors, primarily as an adjunct to conventional imaging and clinical presentation. MRS is useful in initial diagnosis of brain tumors, helping differentiate tumors from possible mimics such as metastatic disease, lymphoma, demyelination, and infection, as well as in the subsequent follow-up of patients after resection and chemoradiation. Unfortunately, the spectroscopic appearance of many pathologies can overlap, and ultimately follow-up or biopsy may be required to make a definitive diagnosis. Future developments may continue to increase the value of MRS for initial diagnosis, treatment planning, and early detection of recurrence.

Advanced magnetic resonance spectroscopic neuroimaging: Experts' consensus recommendations

Magnetic resonance spectroscopic imaging (MRSI) offers considerable promise for monitoring metabolic alterations associated with disease or injury; however, to date, these methods have not had a significant impact on clinical care, and their use remains largely confined to the research community and a limited number of clinical sites. The MRSI methods currently implemented on clinical MRI instruments have remained essentially unchanged for two decades, with only incremental improvements in sequence implementation. During this time, a number of technological developments have taken place that have already greatly benefited the quality of MRSI measurements within the research community and which promise to bring advanced MRSI studies to the point where the technique becomes a true imaging modality, while making the traditional review of individual spectra a secondary requirement. Furthermore, the increasing use of biomedical MR spectroscopy studies has indicated clinical areas where advanced MRSI methods can provide valuable information for clinical care. In light of this rapidly changing technological environment and growing understanding of the value of MRSI studies for biomedical studies, this article presents a consensus from a group of experts in the field that reviews the state-of-the-art for clinical proton MRSI studies of the human brain, recommends minimal standards for further development of vendor-provided MRSI implementations, and identifies areas which need further technical development.

Accelerated J-resolved 1 H-MRSI with limited and sparse sampling of ( k , t 1 , t 2 -space

PURPOSE

To accelerate the acquisition of J-resolved proton magnetic resonance spectroscopic imaging (¹ H-MRSI) data for high-resolution mapping of brain metabolites and neurotransmitters.

METHODS

The proposed method used a subspace model to represent multidimensional spatiospectral functions, which significantly reduced the number of parameters to be determined from J-resolved ¹ H-MRSI data. A semi-LASER-based (Localization by Adiabatic SElective Refocusing) echo-planar spectroscopic imaging (EPSI) sequence was used for data acquisition. The proposed data acquisition scheme sampled k , t 1 , t 2 -space in variable density, where t₁ and t₂ specify the J-coupling and chemical-shift encoding times, respectively. Selection of the J-coupling encoding times (or, echo time values) was based on a Cramer-Rao lower bound analysis, which were optimized for gamma-aminobutyric acid (GABA) detection. In image reconstruction, parameters of the subspace-based spatiospectral model were determined by solving a constrained optimization problem.

RESULTS

Feasibility of the proposed method was evaluated using both simulated and experimental data from a spectroscopic phantom. The phantom experimental results showed that the proposed method, with a factor of 12 acceleration in data acquisition, could determine the distribution of J-coupled molecules with expected accuracy. In vivo study with healthy human subjects also showed that 3D maps of brain metabolites and neurotransmitters can be obtained with a nominal spatial resolution of 3.0 × 3.0 × 4.8 mm³ from J-resolved ¹ H-MRSI data acquired in 19.4 min.

CONCLUSIONS

This work demonstrated the feasibility of highly accelerated J-resolved ¹ H-MRSI using limited and sparse sampling of k , t 1 , t 2 -space and subspace modeling. With further development, the proposed method may enable high-resolution mapping of brain metabolites and neurotransmitters in clinical applications.

Brain metabolites in chronic migraine patients with medication overuse headache

BACKGROUND

Medication overuse headache may be associated with widespread alterations along the thalamocortical pathway, a pathway involved in pain perception and disease progression. This study addressed whether brain metabolites in key regions of the thalamocortical pathway differed between chronic migraine patients with medication overuse headache and without medication overuse headache.

METHODS

Magnetic resonance spectroscopic imaging was used to map metabolites in the bilateral anterior cingulate cortices, mid cingulate cortices, posterior cingulate cortices, and the thalami. Sixteen patients with medication overuse headache were compared with 16 matched patients without medication overuse headache and 16 matched healthy controls.

RESULTS

Glutamate and glutamine in the right mid cingulate cortex and myo-inositol in the left anterior cingulate cortex were significantly higher in patients with medication overuse headache than patients without medication overuse headache, but similar to healthy controls. Both patient groups exhibited reduced N-acetyl-aspartate and creatine in the thalamus, reduced myo-inositol in the right anterior cingulate cortex, and elevated choline in the right mid cingulate cortex. Finally, a negative association between myo-inositol laterality index in the anterior cingulate cortices and number of days per month with acute medication use was found across all patients.

CONCLUSIONS

Patients with medication overuse headache were characterized by a distinct concentration profile of myo-inositol, a glial marker, in the anterior cingulate cortices that may have arisen from medication overuse and could contribute to the development of medication overuse headache.

Ultrafast magnetic resonance spectroscopic imaging using SPICE with learned subspaces

PURPOSE

To develop a subspace learning method for the recently proposed subspace-based MRSI approach known as SPICE, and achieve ultrafast 1 H-MRSI of the brain.

THEORY AND METHODS

A novel strategy is formulated to learn a low-dimensional subspace representation of MR spectra from specially acquired training data and use the learned subspace for general MRSI experiments. Specifically, the subspace learning problem is formulated as learning "empirical" distributions of molecule-specific spectral parameters (e.g., concentrations, lineshapes, and frequency shifts) by integrating physics-based model and the training data. The learned spectral parameters and quantum mechanical simulation basis can then be combined to construct acquisition-specific subspace for spatiospectral encoding and processing. High-resolution MRSI acquisitions combining ultrashort-TE/short-TR excitation, sparse sampling, and the elimination of water suppression have been performed to evaluate the feasibility of the proposed method.

RESULTS

The accuracy of the learned subspace and the capability of the proposed method in producing high-resolution 3D 1 H metabolite maps and high-quality spatially resolved spectra (with a nominal resolution of ∼2.4 × 2.4 × 3 mm3 in 5 minutes) were demonstrated using phantom and in vivo studies. By eliminating water suppression, we are also able to extract valuable information from the water signals for data processing ( B 0 map, frequency drift, and coil sensitivity) as well as for mapping tissue susceptibility and relaxation parameters.

CONCLUSIONS

The proposed method enables ultrafast 1 H-MRSI of the brain using a learned subspace, eliminating the need of acquiring subject-dependent navigator data (known as D 1 ) in the original SPICE technique. It represents a new way to perform MRSI experiments and an important step toward practical applications of high-resolution MRSI.

Non-Cartesian GRAPPA and coil combination using interleaved calibration data - application to concentric-ring MRSI of the human brain at 7T

PURPOSE

Proton MR spectroscopic imaging (MRSI) benefits from B₀ ≥ 7T and multichannel receive coils, promising substantial resolution improvements. However, MRSI acquisition with high spatial resolution requires efficient acceleration and coil combination. To speed up the already-fast sampling via concentric rings, we implemented additional, non-Cartesian, hybrid through-time/through-k-space (tt/tk)-generalized autocalibrating partially parallel acquisition (GRAPPA). A new multipurpose interleaved calibration scan (interleaved MUSICAL) acquires reference data for both coil combination and PI. This renders the reconstruction process (especially PI) less sensitive to instabilities.

METHODS

Six healthy volunteers were scanned at 7T. Three calibration datasets for coil combination and PI were recorded: a) iMUSICAL, b) static MUSICAL as prescan, c) moved MUSICAL as prescan with misaligned head position. The coil combination performance, including motion sensitivity, of iMUSICAL was compared to MUSICAL for single-slice free induction decay (FID)-MRSI. Through-time/through-k-space-GRAPPA with constant/variable-density undersampling was evaluated on the same data, comparing the three calibration datasets. Additionally, the proposed method was successfully applied to 3D whole-brain FID-MRSI.

RESULTS

Using iMUSICAL for coil combination yielded the highest signal-to-noise ratio (SNR) (+9%) and lowest Cramer-Rao lower bounds (CRLBs) (-6%) compared to both MUSICAL approaches, with similar metabolic map quality. Also, excellent mean g-factors of 1.07 and low residual lipid aliasing were obtained when using iMUSICAL as calibration data for two-fold, variable-density undersampling, while significantly degraded metabolic maps were obtained using the misaligned MUSICAL calibration data.

CONCLUSION

Through-time/through-k-space-GRAPPA can accelerate already time-efficient non-Cartesian spatial-spectral 2D/3D-MRSI encoding even further. Particularly promising results have been achieved using iMUSICAL as a robust, interleaved multipurpose calibration for MRSI reconstruction, without extra calibration prescan.

Methodological consensus on clinical proton MRS of the brain: Review and recommendations

Proton MRS (1 H MRS) provides noninvasive, quantitative metabolite profiles of tissue and has been shown to aid the clinical management of several brain diseases. Although most modern clinical MR scanners support MRS capabilities, routine use is largely restricted to specialized centers with good access to MR research support. Widespread adoption has been slow for several reasons, and technical challenges toward obtaining reliable good-quality results have been identified as a contributing factor. Considerable progress has been made by the research community to address many of these challenges, and in this paper a consensus is presented on deficiencies in widely available MRS methodology and validated improvements that are currently in routine use at several clinical research institutions. In particular, the localization error for the PRESS localization sequence was found to be unacceptably high at 3 T, and use of the semi-adiabatic localization by adiabatic selective refocusing sequence is a recommended solution. Incorporation of simulated metabolite basis sets into analysis routines is recommended for reliably capturing the full spectral detail available from short TE acquisitions. In addition, the importance of achieving a highly homogenous static magnetic field (B0 ) in the acquisition region is emphasized, and the limitations of current methods and hardware are discussed. Most recommendations require only software improvements, greatly enhancing the capabilities of clinical MRS on existing hardware. Implementation of these recommendations should strengthen current clinical applications and advance progress toward developing and validating new MRS biomarkers for clinical use.

Fast data acquisition techniques in magnetic resonance spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) is an important technique for assessing the spatial variation of metabolites in vivo. The long scan times in MRSI limit clinical applicability due to patient discomfort, increased costs, motion artifacts, and limited protocol flexibility. Faster acquisition strategies can address these limitations and could potentially facilitate increased adoption of MRSI into routine clinical protocols with minimal addition to the current anatomical and functional acquisition protocols in terms of imaging time. Not surprisingly, a lot of effort has been devoted to the development of faster MRSI techniques that aim to capture the same underlying metabolic information (relative metabolite peak areas and spatial distribution) as obtained by conventional MRSI, in greatly reduced time. The gain in imaging time results, in some cases, in a loss of signal-to-noise ratio and/or in spatial and spectral blurring. This review examines the current techniques and advances in fast MRSI in two and three spatial dimensions and their applications. This review categorizes the acceleration techniques according to their strategy for acquisition of the k-space. Techniques such as fast/turbo-spin echo MRSI, echo-planar spectroscopic imaging, and non-Cartesian MRSI effectively cover the full k-space in a more efficient manner per T_R . On the other hand, techniques such as parallel imaging and compressed sensing acquire fewer k-space points and employ advanced reconstruction algorithms to recreate the spatial-spectral information, which maintains statistical fidelity in test conditions (ie no statistically significant differences on voxel-wise comparisions) with the fully sampled data. The advantages and limitations of each state-of-the-art technique are reviewed in detail, concluding with a note on future directions and challenges in the field of fast spectroscopic imaging.

Reduction of lipid contamination in MR spectroscopy imaging using signal space projection

PURPOSE

Lipid contamination can complicate the metabolite quantification in MR spectroscopic imaging (MRSI). In addition to various experimental methods demonstrated to be feasible for lipid suppression, the postprocessing method is beneficial in the flexibility of applications. In this study, the signal space projection (SSP) algorithm is proposed to suppress the lipid signal in the MRSI.

METHODS

The performance of lipid suppression using SSP and SSP combined with the Papoulis-Gerchberg (PG) algorithm (PG+SSP) is examined in 2D MRSI data and the results were compared with outer volume saturation (OVS) methods. Up to 10 lipid spatial components were extracted by SSP from lipid signals in the range of 0.8~1.5 ppm.

RESULTS

Our results show that most lipid signals were found in the first 4 to 5 components and that lipid signals on the spectra can be suppressed using 4 to 5 components. Metabolites concentrations were quantified using LCModel. Two regions of interest (ROIs) were manually selected on the peripheral and inner brain regions. The quantification of metabolites in terms of fitting reliability (CRLB) and spatial variations within ROIs (SpaVar) is improved using SSP. When 5 to 6 components were used in SSP and PG+SSP, the metabolite concentrations and the associated SpaVar and CRLB are at the same level as those from the OVS.

CONCLUSION

We have demonstrated that the SSP method can be used to suppress the lipid signals of MRSI and SSP with 5 to 6 components is suggested to have a similar suppression performance as the OVS method.

Multiparametric (mp) MRI of prostate cancer

Prostate cancer (PCa) is one of the most prevalent cancers in men. A large number of men are detected with PCa; however, the clinical behavior ranges from low-grade indolent tumors that never develop into a clinically significant disease to aggressive, invasive tumors that may rapidly progress to metastatic disease. The challenges in clinical management of PCa are at levels of screening, diagnosis, treatment, and follow-up after treatment. Magnetic resonance imaging (MRI) methods have shown a potential role in detection, localization, staging, assessment of aggressiveness, targeting biopsies, etc. in PCa patients. Multiparametric MRI (mpMRI) is emerging as a better option compared to the individual imaging methods used in the evaluation of PCa. There are attempts to improve the reproducibility and reliability of mpMRI by using an objective scoring system proposed in the prostate imaging reporting and data system (PIRADS) for standardized reporting. Prebiopsy mpMRI may be used to detect PCa in men with elevated prostate-specific antigen or abnormal digital rectal examination and to enable targeted biopsies. mpMRI can also be used to decide on clinical management of patients, for example active surveillance, and may help in detecting only the pathology that requires detection. It can potentially not only guide patient selection for initial and repeat biopsy but also reduce false-negative biopsies. This review presents a description of the MR methods most commonly applied for investigations of prostate. The anatomical, functional and metabolic parameters obtained from these MR methods are discussed with regard to their physical basis and their contribution to mpMRI investigations of PCa.

Recent advances in parallel imaging for MRI

Magnetic Resonance Imaging (MRI) is an essential technology in modern medicine. However, one of its main drawbacks is the long scan time needed to localize the MR signal in space to generate an image. This review article summarizes some basic principles and recent developments in parallel imaging, a class of image reconstruction techniques for shortening scan time. First, the fundamentals of MRI data acquisition are covered, including the concepts of k-space, undersampling, and aliasing. It is demonstrated that scan time can be reduced by sampling a smaller number of phase encoding lines in k-space; however, without further processing, the resulting images will be degraded by aliasing artifacts. Nearly all modern clinical scanners acquire data from multiple independent receiver coil arrays. Parallel imaging methods exploit properties of these coil arrays to separate aliased pixels in the image domain or to estimate missing k-space data using knowledge of nearby acquired k-space points. Three parallel imaging methods-SENSE, GRAPPA, and SPIRiT-are described in detail, since they are employed clinically and form the foundation for more advanced methods. These techniques can be extended to non-Cartesian sampling patterns, where the collected k-space points do not fall on a rectangular grid. Non-Cartesian acquisitions have several beneficial properties, the most important being the appearance of incoherent aliasing artifacts. Recent advances in simultaneous multi-slice imaging are presented next, which use parallel imaging to disentangle images of several slices that have been acquired at once. Parallel imaging can also be employed to accelerate 3D MRI, in which a contiguous volume is scanned rather than sequential slices. Another class of phase-constrained parallel imaging methods takes advantage of both image magnitude and phase to achieve better reconstruction performance. Finally, some applications are presented of parallel imaging being used to accelerate MR Spectroscopic Imaging.

Fast magnetic resonance spectroscopic imaging techniques in human brain- applications in multiple sclerosis

Multi voxel magnetic resonance spectroscopic imaging (MRSI) is an important imaging tool that combines imaging and spectroscopic techniques. MRSI of the human brain has been beneficially applied to different clinical applications in neurology, particularly in neurooncology but also in multiple sclerosis, stroke and epilepsy. However, a major challenge in conventional MRSI is the longer acquisition time required for adequate signal to be collected. Fast MRSI of the brain in vivo is an alternative approach to reduce scanning time and make MRSI more clinically suitable.Fast MRSI can be categorised into spiral, echo-planar, parallel and turbo imaging techniques, each with its own strengths. After a brief introduction on the basics of non-invasive examination (¹H-MRS) and localization techniques principles, different fast MRSI techniques will be discussed from their initial development to the recent innovations with particular emphasis on their capacity to record neurochemical changes in the brain in a variety of pathologies.The clinical applications of whole brain fast spectroscopic techniques, can assist in the assessment of neurochemical changes in the human brain and help in understanding the roles they play in disease. To give a good example of the utilities of these techniques in clinical context, MRSI application in multiple sclerosis was chosen. The available up to date and relevant literature is discussed and an outline of future research is presented.

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

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

Fast and efficient free induction decay MR spectroscopic imaging of the human brain at 9.4 Tesla

PURPOSE

The purpose of this work was to develop a fast and efficient MRSI-FID acquisition scheme and test its performance in vivo. The aim was to find a trade-off between the minimal total acquisition time and signal-to-noise ratio of the acquired spectra.

METHODS

Measurements were performed on a 9.4 Tesla system. Sequence optimization included redesign of water suppression, optimization of the sequence gradients, and improvement of the sampling efficiency by minimizing the read-out time. This resulted in an acquisition time of 2:47 and 22:13 minutes for 2D (TR = 57 ms; 3-mm in-plane resolution) and 3D MRSI (TR = 57 ms; 16 slices; 3-mm isotropic resolution), respectively.

RESULTS

Despite strong T₁ weighting and first-order phase problems, it was possible to obtain spectra of an acceptable quality. The average line width calculated for the tCr peak across the entire field of view was 26.9 ± 9.6 Hz for 2D and 30.0 ± 11.3 Hz for 3D MRSI. In 3D measurements, the percent fraction of voxels fitted with Cramer-Rao lower bounds below 10% was 53.3 ± 4.1%, 63.4 ± 8.4%, and 81.0 ± 2.9% for Glu, tCr, and tNAA, respectively.

CONCLUSION

Considering the typically long duration of high-resolution MRSI, the proposed technique may be of interest for clinical applications and/or studies that focus on following the biochemistry of dynamic processes. Magn Reson Med 78:1281-1295, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Sparse Parallel MRI Based on Accelerated Operator Splitting Schemes

Recently, the sparsity which is implicit in MR images has been successfully exploited for fast MR imaging with incomplete acquisitions. In this paper, two novel algorithms are proposed to solve the sparse parallel MR imaging problem, which consists of l₁ regularization and fidelity terms. The two algorithms combine forward-backward operator splitting and Barzilai-Borwein schemes. Theoretically, the presented algorithms overcome the nondifferentiable property in l₁ regularization term. Meanwhile, they are able to treat a general matrix operator that may not be diagonalized by fast Fourier transform and to ensure that a well-conditioned optimization system of equations is simply solved. In addition, we build connections between the proposed algorithms and the state-of-the-art existing methods and prove their convergence with a constant stepsize in Appendix. Numerical results and comparisons with the advanced methods demonstrate the efficiency of proposed algorithms.

(2 + 1)D-CAIPIRINHA accelerated MR spectroscopic imaging of the brain at 7T

Purpose

To compare a new parallel imaging (PI) method for multislice proton magnetic resonance spectroscopic imaging (¹H‐MRSI), termed (2 + 1)D‐CAIPIRINHA, with two standard PI methods: 2D‐GRAPPA and 2D‐CAIPIRINHA at 7 Tesla (T).

Methods

(2 + 1)D‐CAIPIRINHA is a combination of 2D‐CAIPIRINHA and slice‐CAIPIRINHA. Eight healthy volunteers were measured on a 7T MR scanner using a 32‐channel head coil. The best undersampling patterns were estimated for all three PI methods. The artifact powers, g‐factors, Cramér–Rao lower bounds (CRLB), and root mean square errors (RMSE) were compared quantitatively among the three PI methods. Metabolic maps and spectra were compared qualitatively.

Results

(2 + 1)D‐CAIPIRINHA allows acceleration in three spatial dimensions in contrast to 2D‐GRAPPA and 2D‐CAIPIRINHA. Thus, this sequence significantly decreased the RMSE of the metabolic maps by 12.1 and 6.9%, on average, for 4 < R < 11, compared with 2D‐GRAPPA and 2D‐CAIPIRINHA, respectively. The artifact power was 22.6 and 8.4% lower, and the CRLB were 3.4 and 0.6% lower, respectively.

Conclusion

(2 + 1)‐CAIPIRINHA can be implemented for multislice MRSI in the brain, enabling higher accelerations than possible with two‐dimensional (2D) parallel imaging methods. An eight‐fold acceleration was still feasible in vivo with negligible PI artifacts with lipid decontamination, thus decreasing the measurement time from 120 to 15 min for a 64 × 64 × 4 matrix. Magn Reson Med 78:429–440, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Regional distributions of brain glutamate and glutamine in normal subjects

Glutamate (Glu) and glutamine (Gln) play an important role in neuronal regulation and are of value as MRS-observable diagnostic biomarkers. In this study the relative concentrations of these metabolites have been measured in multiple regions in the normal brain using a short-TE whole-brain MRSI measurement at 3 T combined with a modified data analysis approach that used spatial averaging to obtain high-SNR spectra from atlas-registered anatomic regions or interest. By spectral fitting of high-SNR spectra this approach yielded reliable measurements across a wide volume of the brain. Spectral averaging also demonstrated increased SNR and improved fitting accuracy for the sum of Glu and Gln (Glx) compared with individual voxel fitting. Results in 26 healthy controls showed relatively constant Glu/Cr and Gln/Cr throughout the cerebrum, although with increased values in the anterior cingulum and paracentral lobule, and increased Gln/Cr in the superior motor area. The deep gray-matter regions of thalamus, putamen, and pallidum show lower Glu/Cr compared with cortical white-matter regions. Lobar measurements demonstrated reduced Glu/Cr and Gln/Cr in the cerebellum as compared with the cerebrum, where white-matter regions show significantly lower Glu/Cr and Gln/Cr as compared with gray-matter regions across multiple brain lobes. Regression analysis showed no significant effect of gender on Glu/Cr or Gln/Cr measurement; however, Glx/Cr ratio was found to be significantly negatively correlated with age in some lobar brain regions. In summary, this methodology provides the spectral quality necessary for reliable separation of Glu and Gln at 3 T from a single MRSI acquisition enabling generation of regional distributions of metabolites over a large volume of the brain, including cortical regions. Copyright © 2016 John Wiley & Sons, Ltd.

Sparse spectral deconvolution algorithm for noncartesian MR spectroscopic imaging

PURPOSE

To minimize line shape distortions and spectral leakage artifacts in MR spectroscopic imaging (MRSI).

METHODS

A spatially and spectrally regularized non-Cartesian MRSI algorithm that uses the line shape distortion priors, estimated from water reference data, to deconvolve the spectra is introduced. Sparse spectral regularization is used to minimize noise amplification associated with deconvolution. A spiral MRSI sequence that heavily oversamples the central k-space regions is used to acquire the MRSI data. The spatial regularization term uses the spatial supports of brain and extracranial fat regions to recover the metabolite spectra and nuisance signals at two different resolutions. Specifically, the nuisance signals are recovered at the maximum resolution to minimize spectral leakage, while the point spread functions of metabolites are controlled to obtain acceptable signal-to-noise ratio.

RESULTS

The comparisons of the algorithm against Tikhonov regularized reconstructions demonstrates considerably reduced line-shape distortions and improved metabolite maps.

CONCLUSION

The proposed sparsity constrained spectral deconvolution scheme is effective in minimizing the line-shape distortions. The dual resolution reconstruction scheme is capable of minimizing spectral leakage artifacts.

Combining parallel detection of proton echo planar spectroscopic imaging (PEPSI) measurements with a data-consistency constraint improves SNR

One major challenge of MRSI is the poor signal-to-noise ratio (SNR), which can be improved by using a surface coil array. Here we propose to exploit the spatial sensitivity of different channels of a coil array to enforce the k-space data consistency (DC) in order to suppress noise and consequently to improve MRSI SNR. MRSI data were collected using a proton echo planar spectroscopic imaging (PEPSI) sequence at 3 T using a 32-channel coil array and were averaged with one, two and eight measurements (avg-1, avg-2 and avg-8). The DC constraint was applied using a regularization parameter λ of 1, 2, 3, 5 or 10. Metabolite concentrations were quantified using LCModel. Our results show that the suppression of noise by applying the DC constraint to PEPSI reconstruction yields up to 32% and 27% SNR gain for avg-1 and avg-2 data with λ = 5, respectively. According to the reported Cramer-Rao lower bounds, the improvement in metabolic fitting was significant (p < 0.01) when the DC constraint was applied with λ ≥ 2. Using the DC constraint with λ = 3 or 5 can minimize both root-mean-square errors and spatial variation for all subjects using the avg-8 data set as reference values. Our results suggest that MRSI reconstructed with a DC constraint can save around 70% of scanning time to obtain images and spectra with similar SNRs using λ = 5.

Lipid suppression via double inversion recovery with symmetric frequency sweep for robust 2D-GRAPPA-accelerated MRSI of the brain at 7 T

This work presents a new approach for high-resolution MRSI of the brain at 7 T in clinically feasible measurement times. Two major problems of MRSI are the long scan times for large matrix sizes and the possible spectral contamination by the transcranial lipid signal. We propose a combination of free induction decay (FID)-MRSI with a short acquisition delay and acceleration via in-plane two-dimensional generalised autocalibrating partially parallel acquisition (2D-GRAPPA) with adiabatic double inversion recovery (IR)-based lipid suppression to allow robust high-resolution MRSI. We performed Bloch simulations to evaluate the magnetisation pathways of lipids and metabolites, and compared the results with phantom measurements. Acceleration factors in the range 2-25 were tested in a phantom. Five volunteers were scanned to verify the value of our MRSI method in vivo. GRAPPA artefacts that cause fold-in of transcranial lipids were suppressed via double IR, with a non-selective symmetric frequency sweep. The use of long, low-power inversion pulses (100 ms) reduced specific absorption rate requirements. The symmetric frequency sweep over both pulses provided good lipid suppression (>90%), in addition to a reduced loss in metabolite signal-to-noise ratio (SNR), compared with conventional IR suppression (52-70%). The metabolic mapping over the whole brain slice was not limited to a rectangular region of interest. 2D-GRAPPA provided acceleration up to a factor of nine for in vivo FID-MRSI without a substantial increase in g-factors (<1.1). A 64 × 64 matrix can be acquired with a common repetition time of ~1.3 s in only 8 min without lipid artefacts caused by acceleration. Overall, we present a fast and robust MRSI method, using combined double IR fat suppression and 2D-GRAPPA acceleration, which may be used in (pre)clinical studies of the brain at 7 T.

High-resolution (1) H-MRSI of the brain using SPICE: Data acquisition and image reconstruction

PURPOSE

To develop data acquisition and image reconstruction methods to enable high-resolution (1) H MR spectroscopic imaging (MRSI) of the brain, using the recently proposed subspace-based spectroscopic imaging framework called SPICE (SPectroscopic Imaging by exploiting spatiospectral CorrElation).

THEORY AND METHODS

SPICE is characterized by the use of a subspace model for both data acquisition and image reconstruction. For data acquisition, we propose a novel spatiospectral encoding scheme that provides hybrid data sets for determining the subspace structure and for image reconstruction using the subspace model. More specifically, we use a hybrid chemical shift imaging /echo-planar spectroscopic imaging sequence for two-dimensional (2D) MRSI and a dual-density, dual-speed echo-planar spectroscopic imaging sequence for three-dimensional (3D) MRSI. For image reconstruction, we propose a method that can determine the subspace structure and the high-resolution spatiospectral reconstruction from the hybrid data sets generated by the proposed sequences, incorporating field inhomogeneity correction and edge-preserving regularization.

RESULTS

Phantom and in vivo brain experiments were performed to evaluate the performance of the proposed method. For 2D MRSI experiments, SPICE is able to produce high-SNR spatiospectral distributions with an approximately 3 mm nominal in-plane resolution from a 10-min acquisition. For 3D MRSI experiments, SPICE is able to achieve an approximately 3 mm in-plane and 4 mm through-plane resolution in about 25 min.

CONCLUSION

Special data acquisition and reconstruction methods have been developed for high-resolution (1) H-MRSI of the brain using SPICE. Using these methods, SPICE is able to produce spatiospectral distributions of (1) H metabolites in the brain with high spatial resolution, while maintaining a good SNR. These capabilities should prove useful for practical applications of SPICE. Magn Reson Med 76:1059-1070, 2016. © 2015 Wiley Periodicals, Inc.

Proton MRS and MRSI of the brain without water suppression

Water suppression (WS) techniques have played a vital role in the commencement and development of in vivo proton magnetic resonance spectroscopy (MRS, including spectroscopic imaging - MRSI). WS not only made in vivo proton MRS functionally available but also made its applications conveniently accessible, and it has become an indispensable tool in most of the routine applications of in vivo proton MR spectroscopy. On the other hand, WS brought forth some challenges. Therefore, various techniques of proton MRS without WS have been developed since the pioneering work in the late 1990s. After more than one and a half decades of advances in both hardware and software, non-water-suppressed proton MRS is coming to the stage of maturity and seeing increasing application in biomedical research and clinical diagnosis. In this article, we will review progress in the technical development and applications of proton MRS without WS.

Feasibility and reproducibility of echo planar spectroscopic imaging on the quantification of hepatic fat

OBJECTIVES

1H-MRS is widely regarded as the most accurate noninvasive method to quantify hepatic fat content (HFC). When practical period of breath holding, and acquisition of HFC over multiple liver areas is considered, a fast MR spectroscopic imaging technique is desired. The aim of this study is to examine the feasibility and reproducibility of echo planar spectroscopic imaging (EPSI) on the quantification of HFC in subject with various HFCs.

METHODS

Twenty two volunteers were examined in a 3T MR system. The acquisition time of proposed EPSI protocol was 18 seconds. The EPSI scans were repeated 8 times for each subject to test reproducibility. The peak of water and individual peaks of fat including methyl, methylene, and allylic peaks at 0.9, 1.3, and 2.0 ppm were fitted. Calculated amount of water and fat content were corrected for T2 relaxation. The total HFC was defined as the combination of individual peaks. Standard deviation (SD), coefficient of variance (COV) and fitting reliability of HFC quantified by LCModel were calculated.

RESULTS

Our results show that the SDs of total HFC for all subjects are less than 2.5%. Fitting reliability is mostly under 10% and positively correlates with COV. Subjects separated into three subgroups according to quantified total HFC show that improved fitting reliability and reproducibility can be achieved on subjects with higher total HFC.

CONCLUSIONS

We have demonstrated feasibility of the proposed EPSI protocols on the quantification of HFC over a whole slice of liver with scan time in a single breath hold.

Statistical mapping of metabolites in the medial wall of the brain: a proton echo planar spectroscopic imaging study

With magnetic resonance spectroscopic imaging (MRSI), it is possible to simultaneously map distributions of several brain metabolites with relatively good spatial resolution in a short time. Although other functional imaging modalities have taken advantage of population-based inferences using spatially extended statistics, this approach remains little utilized for MRSI. In this study, statistical nonparametric mapping (SnPM) was applied to two-dimensional MRSI data from the medial walls of the human brain to assess the effect of normal aging on metabolite concentrations. The effects of different preprocessing steps on these results were then explored. Short echo time MRSI of left and right medial walls was acquired in conjunction with absolute quantification of total choline, total creatine (tCr), glutamate and glutamine, myo-inositol, and N-acetyl-aspartate. Individual images were spatially warped to a common anatomical frame of reference. Age effects were assessed within SnPM as were the effects of voxel subsampling, variance smoothing, and spatial smoothing. The main findings were: (1) regions in the bilateral dorsal anterior cingulate and in the left posterior cingulate exhibited higher tCr concentrations with age; (2) voxel subsampling but not spatial smoothing enhanced the cluster-level statistical sensitivity; and (3) variance smoothing was of little benefit in this study. Our study shows that spatially extended statistics can yield information about regional-specific changes in metabolite concentrations obtained by short echo time MRSI. This opens up the possibility for systematic comparisons of metabolites in the medial wall of the brain.

A pilot validation of multi-echo based echo-planar correlated spectroscopic imaging in human calf muscles

A current limitation of MR spectroscopic imaging of multiple skeletal muscles is prolonged scan duration. A significant reduction in the total scan duration using the echo-planar correlated spectroscopic imaging (EP-COSI) sequence was accomplished using two bipolar readout trains with different phase-encoded echoes for one of two spatial dimensions within a single repetition time (TR). The second bipolar readout was used for spatially encoding the outer k-space, whereas the first readout was used for the central k-space only. The performance of this novel sequence, called multi-echo based echo-planar correlated spectroscopic imaging (ME-EPCOSI), was demonstrated by localizing specific key features in calf muscles and bone marrow of 11 healthy volunteers and five subjects with type 2 diabetes (T2D). A 3 T MRI-MRS scanner equipped with a transmit-receive extremity coil was used. Localization of the ME-EPCOSI sequence was in good agreement with the earlier single-readout based EP-COSI sequence and the required scan time was reduced by a factor of two. In agreement with an earlier report using single-voxel based 2D MRS, significantly increased unsaturated pools of intramyocellular lipid (IMCL) and extramyocellular lipid (EMCL) and decreased IMCL and EMCL unsaturation indices (UIs) were observed in the soleus and tibialis anterior muscle regions of subjects with T2D compared with healthy controls. In addition, significantly decreased choline content was observed in the soleus of T2D subjects compared with healthy controls. Multi-voxel characterization of IMCL and EMCL ratios and UI in the calf muscle may be useful for the non-invasive assessment of altered lipid metabolism in the pathophysiology of T2D.

Accelerated 1 H MRSI using randomly undersampled spiral-based k-space trajectories

PURPOSE

To develop and evaluate the performance of an acquisition and reconstruction method for accelerated MR spectroscopic imaging (MRSI) through undersampling of spiral trajectories.

THEORY AND METHODS

A randomly undersampled spiral acquisition and sensitivity encoding (SENSE) with total variation (TV) regularization, random SENSE+TV, is developed and evaluated on single-slice numerical phantom, in vivo single-slice MRSI, and in vivo three-dimensional (3D)-MRSI at 3 Tesla. Random SENSE+TV was compared with five alternative methods for accelerated MRSI.

RESULTS

For the in vivo single-slice MRSI, random SENSE+TV yields up to 2.7 and 2 times reduction in root-mean-square error (RMSE) of reconstructed N-acetyl aspartate (NAA), creatine, and choline maps, compared with the denoised fully sampled and uniformly undersampled SENSE+TV methods with the same acquisition time, respectively. For the in vivo 3D-MRSI, random SENSE+TV yields up to 1.6 times reduction in RMSE, compared with uniform SENSE+TV. Furthermore, by using random SENSE+TV, we have demonstrated on the in vivo single-slice and 3D-MRSI that acceleration factors of 4.5 and 4 are achievable with the same quality as the fully sampled data, as measured by RMSE of reconstructed NAA map, respectively.

CONCLUSION

With the same scan time, random SENSE+TV yields lower RMSEs of metabolite maps than other methods evaluated. Random SENSE+TV achieves up to 4.5-fold acceleration with comparable data quality as the fully sampled acquisition. Magn Reson Med 74:13-24, 2015. © 2014 Wiley Periodicals, Inc.

Error decomposition for parallel imaging reconstruction using modulation-domain representation of undersampled data

This paper presents a quantitative approach to evaluating and optimizing parallel imaging reconstruction for a clinical requirement. By introducing a "modulation domain representation" for undersampled data, the presented approach decomposes parallel imaging reconstruction error into multiple error components that can be grouped into three categories: image fidelity error, residue aliasing artifacts, and amplified noise. It is experimentally found that these error components have different image-space patterns that compromise imaging quality in different fashions. An error function may be defined as the weighted summation of these error components. By choosing a set of weighting coefficients that can quantify desirable image quality, parallel imaging may be optimized for a clinical requirement. It is found that error decomposition model may improve clinical utility of parallel imaging, providing an application-oriented approach to clinical parallel imaging.

Reproducibility and reliability of short-TE whole-brain MR spectroscopic imaging of human brain at 3T

PURPOSE

A feasibility study of an echo-planar spectroscopic imaging (EPSI) using a short echo time (TE) that trades off sensitivity, compared with other short-TE methods, to achieve whole brain coverage using inversion recovery and spatial oversampling to control lipid bleeding.

METHODS

Twenty subjects were scanned to examine intersubject variance. One subject was scanned five times to examine intrasubject reproducibility. Data were analyzed to determine coefficients of variance (COV) and intraclass correlation coefficient (ICC) for N-acetylaspartate (NAA), total creatine (tCr), total choline (tCho), glutamine/glutamate (Glx), and myo-inositol (mI). Regional metabolite concentrations were derived by using multi-voxel analysis based on lobar-level anatomic regions.

RESULTS

For whole-brain mean values, the intrasubject COVs were 14%, 15%, and 20% for NAA, tCr, and tCho, respectively, and 31% for Glx and mI. The intersubject COVs were up to 6% higher. For regional distributions, the intrasubject COVs were ≤ 5% for NAA, tCr, and tCho; ≤ 9% for Glx; and ≤15% for mI, with about 6% higher intersubject COVs. The ICCs of 5 metabolites were ≥ 0.7, indicating the reliability of the measurements.

CONCLUSION

The present EPSI method enables estimation of the whole-brain metabolite distributions, including Glx and mI with small voxel size, and a reasonable scan time and reproducibility.

Comparison of sagittal and transverse echo planar spectroscopic imaging on the quantification of brain metabolites

PURPOSE

We quantitatively compared sagittal and transverse echo planar spectroscopic imaging (EPSI) on the quantification of metabolite concentrations with consideration of tissue variation. A quantification strategy is proposed to collect the necessary information for quantification of concentrations in a minimized acquisition time.

METHODS

Six transverse and six sagittal EPSI data were collected on healthy volunteers. Metabolite concentrations of N-acetyl-aspartate (NAA), total creatine (tCr), total choline (tCho), myo-inositol (mI), and glutamate and glutamine complex (Glx) were quantified using water scaling with partial volume and relaxation correction. Linear regression analysis was performed to extract concentrations in gray matter (GM) and white matter (WM). The inter- and intrasubject coefficients of variance (CV) were estimated.

RESULTS

Concentrations and fitting errors of sagittal and transverse EPSI were at same level. GM to WM contrast of concentrations was found in NAA, tCr, and tCho. The intersubject CVs revealed greater variability in the sagittal EPSI than in the transverse EPSI. The intrasubject CVs of the transverse EPSI were below 5% for NAA, tCr, and tCho.

CONCLUSION

We showed that quantified concentrations of sagittal and transverse EPSI after partial volume correction are comparable and reproducible. The proposed quantification strategy can be conveniently adapted into various MRI protocols.

A subspace approach to high-resolution spectroscopic imaging

PURPOSE

To accelerate spectroscopic imaging using sparse sampling of (k,t)-space and subspace (or low-rank) modeling to enable high-resolution metabolic imaging with good signal-to-noise ratio.

METHODS

The proposed method, called SPectroscopic Imaging by exploiting spatiospectral CorrElation, exploits a unique property known as partial separability of spectroscopic signals. This property indicates that high-dimensional spectroscopic signals reside in a very low-dimensional subspace and enables special data acquisition and image reconstruction strategies to be used to obtain high-resolution spatiospectral distributions with good signal-to-noise ratio. More specifically, a hybrid chemical shift imaging/echo-planar spectroscopic imaging pulse sequence is proposed for sparse sampling of (k,t)-space, and a low-rank model-based algorithm is proposed for subspace estimation and image reconstruction from sparse data with the capability to incorporate prior information and field inhomogeneity correction.

RESULTS

The performance of the proposed method has been evaluated using both computer simulations and phantom studies, which produced very encouraging results. For two-dimensional spectroscopic imaging experiments on a metabolite phantom, a factor of 10 acceleration was achieved with a minimal loss in signal-to-noise ratio compared to the long chemical shift imaging experiments and with a significant gain in signal-to-noise ratio compared to the accelerated echo-planar spectroscopic imaging experiments.

CONCLUSION

The proposed method, SPectroscopic Imaging by exploiting spatiospectral CorrElation, is able to significantly accelerate spectroscopic imaging experiments, making high-resolution metabolic imaging possible.

Glutamate and glutamine: a review of in vivo MRS in the human brain

Our understanding of the roles that the amino acids glutamate (Glu) and glutamine (Gln) play in the mammalian central nervous system has increased rapidly in recent times. Many conditions are known to exhibit a disturbance in Glu-Gln equilibrium, and the exact relationships between these changed conditions and these amino acids are not fully understood. This has led to increased interest in Glu/Gln quantitation in the human brain in an array of conditions (e.g. mental illness, tumor, neuro-degeneration) as well as in normal brain function. Accordingly, this review has been undertaken to describe the increasing number of in vivo techniques available to study Glu and Gln separately, or pooled as 'Glx'. The present MRS methods used to assess Glu and Gln vary in approach, complexity, and outcome, thus the focus of this review is on a description of MRS acquisition approaches, and an indication of relative utility of each technique rather than brain pathologies associated with Glu and/or Gln perturbation. Consequently, this review focuses particularly on (1) one-dimensional (1)H MRS, (2) two-dimensional (1)H MRS, and (3) one-dimensional (13)C MRS techniques.

Impact of reduced k-space acquisition on pathologic detectability for volumetric MR spectroscopic imaging

PURPOSE

To assess the impact of accelerated acquisitions on the spectral quality of volumetric magnetic resonance spectroscopic imaging (MRSI) and to evaluate their ability in detecting metabolic changes with mild injury.

MATERIALS AND METHODS

The implementation of a generalized autocalibrating partially parallel acquisition (GRAPPA) method for a high-resolution whole-brain echo planar SI (3D-EPSI) sequence is first described and the spectral accuracy of the GRAPPA-EPSI method is investigated using lobar and voxel-based analyses for normal subjects and patients with mild traumatic brain injuries (mTBI). The performance of GRAPPA was compared with that of fully encoded EPSI for five datasets collected from normal subjects at the same scanning session, as well as on 45 scans (20 normal subjects and 25 mTBI patients) for which the reduced k-space sampling was simulated. For comparison, a central k-space lower-resolution 3D-EPSI acquisition was also simulated. Differences in individual metabolites and metabolite ratio distributions of the mTBI group relative to those of age-matched control subjects were statistically evaluated using analyses divided into hemispheric brain lobes and tissue types.

RESULTS

GRAPPA-EPSI with 16-minute scan time yielded robust and similar results in terms of MRSI quantitation, spectral fitting, and accuracy with that of fully sampled 3D-EPSI acquisitions and was more accurate than central k-space acquisition. Primary findings included high correlations (accuracy of 92.6%) between the GRAPPA and fully sampled results.

CONCLUSION

Although the reduced encoding method is associated with lower signal-to-noise ratio (SNR) that impacts the quality of spectral analysis, the use of the parallel imaging method can lead to the same diagnostic outcomes as the fully sampled data when using the sensitivity-limited volumetric MRSI.

Strategies for rapid in vivo 1H and hyperpolarized 13C MR spectroscopic imaging

In vivo MRSI is an important imaging modality that has been shown in numerous research studies to give biologically relevant information for assessing the underlying mechanisms of disease and for monitoring response to therapy. The increasing availability of high field scanners and multichannel radiofrequency coils has provided the opportunity to acquire in vivo data with significant improvements in sensitivity and signal to noise ratio. These capabilities may be used to shorten acquisition time and provide increased coverage. The ability to acquire rapid, volumetric MRSI data is critical for examining heterogeneity in metabolic profiles and for relating serial changes in metabolism within the same individual during the course of the disease. In this review we discuss the implementation of strategies that use alternative k-space sampling trajectories and parallel imaging methods in order to speed up data acquisition. The impact of such methods is demonstrated using three recent examples of how these methods have been applied. These are to the acquisition of robust 3D (1)H MRSI data within 5-10 min at a field strength of 3 T, to obtaining higher sensitivity for (1)H MRSI at 7 T and to using ultrafast volumetric and dynamic (13)C MRSI for monitoring the changes in signals that occur following the injection of hyperpolarized (13)C agents.

Short- and long-term quantitation reproducibility of brain metabolites in the medial wall using proton echo planar spectroscopic imaging

Proton echo planar spectroscopic imaging (PEPSI) is a fast magnetic resonance spectroscopic imaging (MRSI) technique that allows mapping spatial metabolite distributions in the brain. Although the medial wall of the cortex is involved in a wide range of pathological conditions, previous MRSI studies have not focused on this region. To decide the magnitude of metabolic changes to be considered significant in this region, the reproducibility of the method needs to be established. The study aims were to establish the short- and long-term reproducibility of metabolites in the right medial wall and to compare regional differences using a constant short-echo time (TE30) and TE averaging (TEavg) optimized to yield glutamatergic information. 2D sagittal PEPSI was implemented at 3T using a 32 channel head coil. Acquisitions were repeated immediately and after approximately 2 weeks to assess the coefficients of variation (COV). COVs were obtained from eight regions-of-interest (ROIs) of varying size and location. TE30 resulted in better spectral quality and similar or lower quantitation uncertainty for all metabolites except glutamate (Glu). When Glu and glutamine (Gln) were quantified together (Glx) reduced quantitation uncertainty and increased reproducibility was observed for TE30. TEavg resulted in lowered quantitation uncertainty for Glu but in less reliable quantification of several other metabolites. TEavg did not result in a systematically improved short- or long-term reproducibility for Glu. The ROI volume was a major factor influencing reproducibility. For both short- and long-term repetitions, the Glu COVs obtained with TEavg were 5-8% for the large ROIs, 12-17% for the medium sized ROIs and 16-26% for the smaller cingulate ROIs. COVs obtained with TE30 for the less specific Glx were 3-5%, 8-10% and 10-15%. COVs for N-acetyl aspartate, creatine and choline using TE30 with long-term repetition were between 2-10%. Our results show that the cost of more specific glutamatergic information (Glu versus Glx) is the requirement of an increased effect size especially with increasing anatomical specificity. This comes in addition to the loss of sensitivity for other metabolites. Encouraging results were obtained with TE30 compared to other previously reported MRSI studies. The protocols implemented here are reliable and may be used to study disease progression and intervention mechanisms.

High-resolution mapping of human brain metabolites by free induction decay (1)H MRSI at 7 T

This work describes a new approach for high-spatial-resolution (1)H MRSI of the human brain at 7 T. (1)H MRSI at 7 T using conventional approaches, such as point-resolved spectroscopy and stimulated echo acquisition mode with volume head coils, is limited by technical difficulties, including chemical shift displacement errors, B(0)/B(1) inhomogeneities, a high specific absorption rate and decreased T(2) relaxation times. The method presented here is based on free induction decay acquisition with an ultrashort acquisition delay (TE*) of 1.3 ms. This allows full signal detection with negligible T(2) decay or J-modulation. Chemical shift displacement errors were reduced to below 5% per part per million in the in-slice direction and were eliminated in-plane. The B(1) sensitivity was reduced significantly and further corrected using flip angle maps. Specific absorption rate requirements were well below the limit (~20 % = 0.7 W/kg). The suppression of subcutaneous lipid signals was achieved by substantially improving the point-spread function. High-quality metabolic mapping of five important brain metabolites was achieved with high in-plane resolution (64 × 64 matrix with a 3.4 × 3.4 × 12 mm(3) nominal voxel size) in four healthy subjects. The ultrashort TE* increased the signal-to-noise ratio of J-coupled resonances, such as glutamate and myo-inositol, several-fold to enable the mapping of even these metabolites with high resolution. Four measurement repetitions in one healthy volunteer provided proof of the good reproducibility of this method. The high spatial resolution allowed the visualization of several anatomical structures on metabolic maps. Free induction decay MRSI is insensitive to T(2) decay, J-modulation, B(1) inhomogeneities and chemical shift displacement errors, and overcomes specific absorption rate restrictions at ultrahigh magnetic fields. This makes it a promising method for high-resolution (1)H MRSI at 7 T and above.

Neurologic 3D MR spectroscopic imaging with low-power adiabatic pulses and fast spiral acquisition

PURPOSE

To improve clinical three-dimensional (3D) MR spectroscopic imaging with more accurate localization and faster acquisition schemes.

MATERIALS AND METHODS

Institutional review board approval and patient informed consent were obtained. Data were acquired with a 3-T MR imager and a 32-channel head coil in phantoms, five healthy volunteers, and five patients with glioblastoma. Excitation was performed with localized adiabatic spin-echo refocusing (LASER) by using adiabatic gradient-offset independent adiabaticity wideband uniform rate and smooth truncation (GOIA-W[16,4]) pulses with 3.5-msec duration, 20-kHz bandwidth, 0.81-kHz amplitude, and 45-msec echo time. Interleaved constant-density spirals simultaneously encoded one frequency and two spatial dimensions. Conventional phase encoding (PE) (1-cm3 voxels) was performed after LASER excitation and was the reference standard. Spectra acquired with spiral encoding at similar and higher spatial resolution and with shorter imaging time were compared with those acquired with PE. Metabolite levels were fitted with software, and Bland-Altman analysis was performed.

RESULTS

Clinical 3D MR spectroscopic images were acquired four times faster with spiral protocols than with the elliptical PE protocol at low spatial resolution (1 cm3). Higher-spatial-resolution images (0.39 cm3) were acquired twice as fast with spiral protocols compared with the low-spatial-resolution elliptical PE protocol. A minimum signal-to-noise ratio (SNR) of 5 was obtained with spiral protocols under these conditions and was considered clinically adequate to reliably distinguish metabolites from noise. The apparent SNR loss was not linear with decreasing voxel sizes because of longer local T2* times. Improvement of spectral line width from 4.8 Hz to 3.5 Hz was observed at high spatial resolution. The Bland-Altman agreement between spiral and PE data is characterized by narrow 95% confidence intervals for their differences (0.12, 0.18 of their means). GOIA-W(16,4) pulses minimize chemical-shift displacement error to 2.1%, reduce nonuniformity of excitation to 5%, and eliminate the need for outer volume suppression.

CONCLUSION

The proposed adiabatic spiral 3D MR spectroscopic imaging sequence can be performed in a standard clinical MR environment. Improvements in image quality and imaging time could enable more routine acquisition of spectroscopic data than is possible with current pulse sequences.

Derivative encoding for parallel magnetic resonance imaging

PURPOSE

To introduce a linear shift-invariant relationship between the partial derivatives of k space signals acquired using multichannel receive coils and to demonstrate that k space derivatives can be used for image unwrapping.

METHODS

Fourier transform of k space derivatives contains information on the spatial origins of aliased pixels; therefore, images can be reconstructed by k space derivatives. Fully sampled phantom and brain images acquired at 3 T using a standard eight channel receive coil were used to validate the k space derivatives theorem by unwrapping aliased images.

RESULTS

Derivative encoding leads to new methods for parallel imaging reconstruction in both k space and image domains. Noise amplification in sensitivity encoding image reconstruction, which is considered to produce the optimal SNR, can be further reduced using k space derivative encoding without making any assumptions on the characteristics of the images to be reconstructed.

CONCLUSIONS

This work demonstrated that the partial derivative of the k space signal acquired from one coil with respect to one direction can be expressed as a sum of partial derivatives of signals from multiple coils with respect to the perpendicular k space direction(s). This relationship between the partial derivatives of k space signals is linear and shift-invariant in the Cartesian coordinate system.

MR spectroscopy and spectroscopic imaging of the brain

Magnetic resonance spectroscopy (MRS) and the related technique of magnetic resonance spectroscopic imaging (MRSI) are widely used in both clinical and preclinical research for the non-invasive evaluation of brain metabolism. They are also used in medical practice, although their ultimate clinical value continues to be a source of discussion. This chapter reviews the general information content of brain spectra and commonly used protocols for both MRS and MRSI and also touches on data analysis methods and quantitation. The main focus is on proton MRS for application in humans, but many of the methods are also applicable to other nuclei and studies of animal models as well.

Potential impact of a 32-channel receiving head coil technology on the results of a functional MRI paradigm

PURPOSE

The authors investigated the potential of a 32-channel (32ch) receiving head coil for functional magnetic resonance imaging (fMRI) compared to a standard eight-channel (8ch) coil using a motor task.

MATERIAL AND METHODS

Brain activation was analyzed in 14 healthy right-handed subjects performing finger tapping with the right index finger (block design) during two experimental sessions, one with the 8ch and one with the 32ch coil (applied in a pseudorandomized order). Additionally, a phantom study was performed to compare signal-to-noise ratios (SNRs) of both coils.

RESULTS

During both fMRI sessions, analysis of motor conditions resulted in an activation of the left "hand knob" (precentral gyrus). Application of the 32ch coil obtained additional activation clusters in the right cerebellum, left superior frontal gyrus (SMA), left supramarginal gyrus, and left postcentral gyrus. The phantom study revealed a significantly higher SNR for the 32ch coil compared to the 8ch coil in superficial cortical areas located near the surface of the brain.

CONCLUSION

The 32ch technology has a potential impact on fMRI studies, especially in paradigms that result in activation of cortical areas located near the surface of the brain.