Mapping of brain macromolecules and their use for spectral processing of 1 H-MRSI data with an ultra-short acquisition delay at 7 T

GABA-edited MEGA-PRESS at 3 T: Does a measured macromolecule background improve linear combination modeling?

PURPOSE

The J-difference edited γ-aminobutyric acid (GABA) signal is contaminated by other co-edited signals-the largest of which originates from co-edited macromolecules (MMs)-and is consequently often reported as "GABA+." MM signals are broader and less well-characterized than the metabolites, and are commonly approximated using a Gaussian model parameterization. Experimentally measured MM signals are a consensus-recommended alternative to parameterized modeling; however, they are relatively under-studied in the context of edited MRS.

METHODS

To address this limitation in the literature, we have acquired GABA-edited MEGA-PRESS data with pre-inversion to null metabolite signals in 13 healthy controls. An experimental MM basis function was derived from the mean across subjects. We further derived a new parameterization of the MM signals from the experimental data, using multiple Gaussians to accurately represent their observed asymmetry. The previous single-Gaussian parameterization, mean experimental MM spectrum and new multi-Gaussian parameterization were compared in a three-way analysis of a public MEGA-PRESS dataset of 61 healthy participants.

RESULTS

Both the experimental MMs and the multi-Gaussian parameterization exhibited reduced fit residuals compared to the single-Gaussian approach (p = 0.034 and p = 0.031, respectively), suggesting they better represent the underlying data than the single-Gaussian parameterization. Furthermore, both experimentally derived models estimated larger MM fractional contribution to the GABA+ signal for the experimental MMs (58%) and multi-Gaussian parameterization (58%), compared to the single-Gaussian approach (50%).

CONCLUSIONS

Our results indicate that single-Gaussian parameterization of edited MM signals is insufficient and that both experimentally derived GABA+ spectra and their parameterized replicas improve the modeling of GABA+ spectra.

7 Tesla magnetic resonance spectroscopic imaging predicting IDH status and glioma grading

INTRODUCTION

With the application of high-resolution 3D 7 Tesla Magnetic Resonance Spectroscopy Imaging (MRSI) in high-grade gliomas, we previously identified intratumoral metabolic heterogeneities. In this study, we evaluated the potential of 3D 7 T-MRSI for the preoperative noninvasive classification of glioma grade and isocitrate dehydrogenase (IDH) status. We demonstrated that IDH mutation and glioma grade are detectable by ultra-high field (UHF) MRI. This technique might potentially optimize the perioperative management of glioma patients.

METHODS

We prospectively included 36 patients with WHO 2021 grade 2-4 gliomas (20 IDH mutated, 16 IDH wildtype). Our 7 T 3D MRSI sequence provided high-resolution metabolic maps (e.g., choline, creatine, glutamine, and glycine) of these patients' brains. We employed multivariate random forest and support vector machine models to voxels within a tumor segmentation, for classification of glioma grade and IDH mutation status.

RESULTS

Random forest analysis yielded an area under the curve (AUC) of 0.86 for multivariate IDH classification based on metabolic ratios. We distinguished high- and low-grade tumors by total choline (tCho) / total N-acetyl-aspartate (tNAA) ratio difference, yielding an AUC of 0.99. Tumor categorization based on other measured metabolic ratios provided comparable accuracy.

CONCLUSIONS

We successfully classified IDH mutation status and high- versus low-grade gliomas preoperatively based on 7 T MRSI and clinical tumor segmentation. With this approach, we demonstrated imaging based tumor marker predictions at least as accurate as comparable studies, highlighting the potential application of MRSI for pre-operative tumor classifications.

Proton Free Induction Decay MRSI at 7T in the Human Brain Using an Egg-Shaped Modified Rosette K-Space Trajectory

Purpose

1.1 Proton ( 1 H)-MRSI via spatial-spectral encoding poses high demands on gradient hardware at ultra-high fields and high-resolutions. Rosette trajectories help alleviate these problems, but at reduced SNR-efficiency due to their k-space densities not matching any desired k-space filter. We propose modified rosette trajectories, which more closely match a Hamming filter, and thereby improve SNR performance while still staying within gradient hardware limitations and without prolonging scan time.

Methods

1.2Analytical and synthetic simulations were validated with phantom and in vivo measurements at 7 T. The rosette and modified rosette trajectories were measured in five healthy volunteers in six minutes in a 2D slice in the brain. A 3D sequence was measured in one volunteer within 19 minutes. The SNR, linewidth, CRLBs, lipid contamination and data quality of the proposed modified rosette trajectory were compared to the rosette trajectory.

Results

1.3Using the modified rosette trajectories, an improved k-space weighting function was achieved resulting in an increase of up to 12% in SNR compared to rosette's dependent on the two additional trajectory parameters. Similar results were achieved for the theoretical SNR calculation based on k-space densities, as well as when using the pseudo-replica method for simulated, in-vivo and phantom data. The CRLBs improved slightly, but non-significantly for the modified rosette trajectories, while the linewidths and lipid contamination remained similar.

Conclusion

1.4By improving the rosette trajectory's shape, modified rosette trajectories achieved higher SNR at the same scan time and data quality.

A Comparison of 7 Tesla MR Spectroscopic Imaging and 3 Tesla MR Fingerprinting for Tumor Localization in Glioma Patients

This paper investigated the correlation between magnetic resonance spectroscopic imaging (MRSI) and magnetic resonance fingerprinting (MRF) in glioma patients by comparing neuro-oncological markers obtained from MRSI to T1/T2 maps from MRF. Data from 12 consenting patients with gliomas were analyzed by defining hotspots for T1, T2, and various metabolic ratios, and comparing them using Sørensen-Dice similarity coefficients (DSCs) and the distances between their centers of intensity (COIDs). The median DSCs between MRF and the tumor segmentation were 0.73 (T1) and 0.79 (T2). The DSCs between MRSI and MRF were the highest for Gln/tNAA (T1: 0.75, T2: 0.80, tumor: 0.78), followed by Gly/tNAA (T1: 0.57, T2: 0.62, tumor: 0.54) and tCho/tNAA (T1: 0.61, T2: 0.58, tumor: 0.45). The median values in the tumor hotspot were T1 = 1724 ms, T2 = 86 ms, Gln/tNAA = 0.61, Gly/tNAA = 0.28, Ins/tNAA = 1.15, and tCho/tNAA = 0.48, and, in the peritumoral region, were T1 = 1756 ms, T2 = 102 ms, Gln/tNAA = 0.38, Gly/tNAA = 0.20, Ins/tNAA = 1.06, and tCho/tNAA = 0.38, and, in the NAWM, were T1 = 950 ms, T2 = 43 ms, Gln/tNAA = 0.16, Gly/tNAA = 0.07, Ins/tNAA = 0.54, and tCho/tNAA = 0.20. The results of this study constitute the first comparison of 7T MRSI and 3T MRF, showing a good correspondence between these methods.

High-resolution proton metabolic mapping of the human brain at 7 T using free induction decay rosette spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) provides information about the spatial distribution of metabolites in the brain. These metabolite maps can be valuable in diagnosing central nervous system pathology. However, MRSI generally suffers from a long acquisition time, poor spatial resolution, and a low metabolite signal-to-noise ratio (SNR). Ultrahigh field strengths (≥ 7 T) can benefit MRSI with an improved SNR and allow high-resolution metabolic mapping. Non-Cartesian spatial-spectral encoding techniques, such as rosette spectroscopic imaging, can efficiently sample spatial and temporal domains, which significantly reduces the imaging time and enables high-resolution metabolic mapping in a clinically relevant scan time. In the current study, high-resolution (in-plane resolution of 2 × 2 mm2 ) mapping of proton (1 H) metabolites in the human brain at 7 T, is demonstrated. Five healthy subjects participated in the study. Using a time-efficient rosette trajectory and short TR/TE free induction decay MRSI, high-resolution maps of 1 H metabolites were obtained in a clinically relevant imaging time (6 min). Suppression of the water signal was achieved with an optimized water suppression enhanced through T1 effects approach and lipid removal was performed using L2 -regularization in the postprocessing. Spatial distributions of N-acetyl-aspartate, total choline, creatine, N-acetyl-aspartyl glutamate, myo-inositol, and glutamate were generated with Cramer-Rao lower bounds of less than 20%.

Macromolecule modelling for improved metabolite quantification using short echo time brain 1 H MRS at 3 T and 7 T: The PRaMM Model

Purpose

To improve reliability of metabolite quantification at both, 3 T and 7 T, we propose a novel parametrized macromolecules quantification model (PRaMM) for brain 1 H MRS, in which the ratios of macromolecule peak intensities are used as soft constraints.

Methods

Full- and metabolite-nulled spectra were acquired in three different brain regions with different ratios of grey and white matter from six healthy volunteers, at both 3 T and 7 T. Metabolite-nulled spectra were used to identify highly correlated macromolecular signal contributions and estimate the ratios of their intensities. These ratios were then used as soft constraints in the proposed PRaMM model for quantification of full spectra. The PRaMM model was validated by comparison with a single component macromolecule model and a macromolecule subtraction technique. Moreover, the influence of the PRaMM model on the repeatability and reproducibility compared to those other methods was investigated.

Results

The developed PRaMM model performed better than the two other approaches in all three investigated brain regions. Several estimates of metabolite concentration and their Cramér-Rao lower bounds were affected by the PRaMM model reproducibility, and repeatability of the achieved concentrations were tested by evaluating the method on a second repeated acquisitions dataset. While the observed effects on both metrics were not significant, the fit quality metrics were improved for the PRaMM method (p≤0.0001). Minimally detectable changes are in the range 0.5 - 1.9 mM and percent coefficients of variations are lower than 10% for almost all the clinically relevant metabolites. Furthermore, potential overparameterization was ruled out.

Conclusion

Here, the PRaMM model, a method for an improved quantification of metabolites was developed, and a method to investigate the role of the MM background and its individual components from a clinical perspective is proposed.

Quantification of short echo time MRS signals with improved version of QUantitation based on quantum ESTimation algorithm

Magnetic resonance spectroscopy offers information about metabolite changes in the organism, which can be used in diagnosis. While short echo time proton spectra exhibit more distinguishable metabolites compared with proton spectra acquired with long echo times, their quantification (and providing estimates of metabolite concentrations) is more challenging. They are hampered by a background signal, which originates mainly from macromolecules (MM) and mobile lipids. An improved version of the quantification algorithm QUantitation based on quantum ESTimation (QUEST), with MM prior knowledge (QUEST-MM), dedicated to proton signals and invoking appropriate prior knowledge on MM, is proposed and tested. From a single acquisition, it enables better metabolite quantification, automatic estimation of the background, and additional automatic quantification of MM components, thus improving its applicability in the clinic. The proposed algorithm may facilitate studies that involve patients with pathological MM in the brain. QUEST-MM and three QUEST-based strategies for quantifying short echo time signals are compared in terms of bias-variance trade-off and Cramér-Rao lower bound estimates. The performances of the methods are evaluated through extensive Monte Carlo studies. In particular, the histograms of the metabolite and MM amplitude distributions demonstrate the performances of the estimators. They showed that QUEST-MM works better than QUEST (Subtract approach) and is a good alternative to QUEST when measured MM signal is unavailable or unsuitable. Quantification with QUEST-MM is shown for 1 H in vivo rat brain signals obtained with the SPECIAL pulse sequence at 9.4 T, and human brain signals obtained, respectively, with STEAM at 4 T and PRESS at 3 T. QUEST-MM is implemented in jMRUI and will be available for public use from version 7.1.

Metabolic Insights into Iron Deposition in Relapsing-Remitting Multiple Sclerosis via 7 T Magnetic Resonance Spectroscopic Imaging

OBJECTIVE

To investigate the metabolic pattern of different types of iron accumulation in multiple sclerosis (MS) lesions, and compare metabolic alterations within and at the periphery of lesions and newly emerging lesions in vivo according to iron deposition.

METHODS

7 T MR spectroscopic imaging and susceptibility-weighted imaging was performed in 31 patients with relapsing-remitting MS (16 female/15 male; mean age, 36.9 ± 10.3 years). Mean metabolic ratios of four neuro-metabolites were calculated for regions of interest (ROI) of normal appearing white matter (NAWM), "non-iron" (lesion without iron accumulation on SWI), and three distinct types of iron-laden lesions ("rim": distinct rim-shaped iron accumulation; "area": iron deposition across the entire lesions; "transition": transition between "area" and "rim" accumulation shape), and for lesion layers of "non-iron" and "rim" lesions. Furthermore, newly emerging "non-iron" and "iron" lesions were compared longitudinally, as measured before their appearance and one year later.

RESULTS

Thirty-nine of 75 iron-containing lesions showed no distinct paramagnetic rim. Of these, "area" lesions exhibited a 65% higher mIns/tNAA (p = 0.035) than "rim" lesions. Comparing lesion layers of both "non-iron" and "rim" lesions, a steeper metabolic gradient of mIns/tNAA ("non-iron" +15%, "rim" +40%) and tNAA/tCr ("non-iron" -15%, "rim" -35%) was found in "iron" lesions, with the lesion core showing +22% higher mIns/tNAA (p = 0.005) and -23% lower tNAA/tCr (p = 0.048) in "iron" compared to "non-iron" lesions. In newly emerging lesions, 18 of 39 showed iron accumulation, with the drop in tNAA/tCr after lesion formation remaining significantly lower compared to pre-lesional tissue over time in "iron" lesions (year 0: p = 0.013, year 1: p = 0.041) as opposed to "non-iron" lesions (year 0: p = 0.022, year 1: p = 0.231).

CONCLUSION

7 T MRSI allows in vivo characterization of different iron accumulation types each presenting with a distinct metabolic profile. Furthermore, the larger extent of neuronal damage in lesions with a distinct iron rim was reconfirmed via reduced tNAA/tCr concentrations, but with metabolic differences in lesion development between (non)-iron-containing lesions. This highlights the ability of MRSI to further investigate different types of iron accumulation and suggests possible implications for disease monitoring.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct (2H) and indirect (1H) detection of deuterium labeled compounds at 7T and clinical 3T

INTRODUCTION

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'-2H2]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using 2H MRSI (DMI) and 1H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

METHODS

Five volunteers (4 m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8 g/kg oral [6,6'-2H2]-glucose administration using time-resolved 3D 2H FID-MRSI with elliptical phase encoding at 7T and 3D 1H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

RESULTS

One hour after oral tracer administration regionally averaged deuterium labeled Glx4 concentrations and the dynamics were not significantly different over all participants between 7T 2H DMI and 3T 1H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc6 data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual 2H and 1H data points a weak to moderate negative correlation was observed for Glx4 concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc6 data GM (r=-0.61, p<0.001) and WM (r=-0.70, p<0.001).

CONCLUSION

This study demonstrates that indirect detection of deuterium labeled compounds using 1H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to 2H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

1H magnetic resonance spectroscopic imaging of deuterated glucose and of neurotransmitter metabolism at 7 T in the human brain

Impaired glucose metabolism in the brain has been linked to several neurological disorders. Positron emission tomography and carbon-13 magnetic resonance spectroscopic imaging (MRSI) can be used to quantify the metabolism of glucose, but these methods involve exposure to radiation, cannot quantify downstream metabolism, or have poor spatial resolution. Deuterium MRSI (2H-MRSI) is a non-invasive and safe alternative for the quantification of the metabolism of 2H-labelled substrates such as glucose and their downstream metabolic products, yet it can only measure a limited number of deuterated compounds and requires specialized hardware. Here we show that proton MRSI (1H-MRSI) at 7 T has higher sensitivity, chemical specificity and spatiotemporal resolution than 2H-MRSI. We used 1H-MRSI in five volunteers to differentiate glutamate, glutamine, γ-aminobutyric acid and glucose deuterated at specific molecular positions, and to simultaneously map deuterated and non-deuterated metabolites. 1H-MRSI, which is amenable to clinically available magnetic-resonance hardware, may facilitate the study of glucose metabolism in the brain and its potential roles in neurological disorders.

Noninvasive 3-Dimensional 1 H-Magnetic Resonance Spectroscopic Imaging of Human Brain Glucose and Neurotransmitter Metabolism Using Deuterium Labeling at 3T : Feasibility and Interscanner Reproducibility

OBJECTIVES

Noninvasive, affordable, and reliable mapping of brain glucose metabolism is of critical interest for clinical research and routine application as metabolic impairment is linked to numerous pathologies, for example, cancer, dementia, and depression. A novel approach to map glucose metabolism noninvasively in the human brain has been presented recently on ultrahigh-field magnetic resonance (MR) scanners (≥7T) using indirect detection of deuterium-labeled glucose and downstream metabolites such as glutamate, glutamine, and lactate. The aim of this study was to demonstrate the feasibility to noninvasively detect deuterium-labeled downstream glucose metabolites indirectly in the human brain via 3-dimensional (3D) proton ( 1 H) MR spectroscopic imaging on a clinical 3T MR scanner without additional hardware.

MATERIALS AND METHODS

This prospective, institutional review board-approved study was performed in 7 healthy volunteers (mean age, 31 ± 4 years, 5 men/2 women) after obtaining written informed consent. After overnight fasting and oral deuterium-labeled glucose administration, 3D metabolic maps were acquired every ∼4 minutes with ∼0.24 mL isotropic spatial resolution using real-time motion-, shim-, and frequency-corrected echo-less 3D 1 H-MR spectroscopic Imaging on a clinical routine 3T MR system. To test the interscanner reproducibility of the method, subjects were remeasured on a similar 3T MR system. Time courses were analyzed using linear regression and nonparametric statistical tests. Deuterium-labeled glucose and downstream metabolites were detected indirectly via their respective signal decrease in dynamic 1 H MR spectra due to exchange of labeled and unlabeled molecules.

RESULTS

Sixty-five minutes after deuterium-labeled glucose administration, glutamate + glutamine (Glx) signal intensities decreased in gray/white matter (GM/WM) by -1.63 ± 0.3/-1.0 ± 0.3 mM (-13% ± 3%, P = 0.02/-11% ± 3%, P = 0.02), respectively. A moderate to strong negative correlation between Glx and time was observed in GM/WM ( r = -0.64, P < 0.001/ r = -0.54, P < 0.001), with 60% ± 18% ( P = 0.02) steeper slopes in GM versus WM, indicating faster metabolic activity. Other nonlabeled metabolites showed no significant changes. Excellent intrasubject repeatability was observed across scanners for static results at the beginning of the measurement (coefficient of variation 4% ± 4%), whereas differences were observed in individual Glx dynamics, presumably owing to physiological variation of glucose metabolism.

CONCLUSION

Our approach translates deuterium metabolic imaging to widely available clinical routine MR scanners without specialized hardware, offering a safe, affordable, and versatile (other substances than glucose can be labeled) approach for noninvasive imaging of glucose and neurotransmitter metabolism in the human brain.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct ( 2 H) and indirect ( 1 H) detection of deuterium labeled compounds at 7T and clinical 3T

Introduction

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'- ² H ₂ ]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using ² H MRSI (DMI) and ¹ H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

Methods

Five volunteers (4m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8g/kg oral [6,6'- ² H ₂ ]-glucose administration using time-resolved 3D ² H FID-MRSI with elliptical phase encoding at 7T and 3D ¹ H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

Results

One hour after oral tracer administration regionally averaged deuterium labeled Glx ₄ concentrations and the dynamics were not significantly different over all participants between 7T ² H DMI and 3T ¹ H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc ₆ data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual ² H and ¹ H data points a weak to moderate negative correlation was observed for Glx ₄ concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc ₆ data GM (r=- 0.61, p<0.001) and WM (r=-0.70, p<0.001).

Conclusion

This study demonstrates that indirect detection of deuterium labeled compounds using ¹ H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to ² H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

The influence of season on glutamate and GABA levels in the healthy human brain investigated by magnetic resonance spectroscopy imaging

Seasonal changes in neurotransmitter systems have been demonstrated in imaging studies and are especially noticeable in diseased states such as seasonal affective disorder (SAD). These modulatory neurotransmitters, such as serotonin, are influencing glutamatergic and GABAergic neurotransmission. Furthermore, central components of the circadian pacemaker are regulated by GABA (the suprachiasmatic nucleus) or glutamate (e.g., the retinohypothalamic tract). Therefore, we explored seasonal differences in the GABAergic and glutamatergic system in 159 healthy individuals using magnetic resonance spectroscopy imaging with a GABA-edited 3D-MEGA-LASER sequence at 3T. We quantified GABA+/tCr, GABA+/Glx, and Glx/tCr ratios (GABA+, GABA+ macromolecules; Glx, glutamate + glutamine; tCr, total creatine) in five different subcortical brain regions. Differences between time periods throughout the year, seasonal patterns, and stationarity were tested using ANCOVA models, curve fitting approaches, and unit root and stationarity tests, respectively. Finally, Spearman correlation analyses between neurotransmitter ratios within each brain region and cumulated daylight and global radiation were performed. No seasonal or monthly differences, seasonal patterns, nor significant correlations could be shown in any region or ratio. Unit root and stationarity tests showed stable patterns of GABA+/tCr, GABA+/Glx, and Glx/tCr levels throughout the year, except for hippocampal Glx/tCr. Our results indicate that neurotransmitter levels of glutamate and GABA in healthy individuals are stable throughout the year. Hence, despite the important correction for age and gender in the analyses of MRS derived GABA and glutamate, a correction for seasonality in future studies does not seem necessary. Future investigations in SAD and other psychiatric patients will be of high interest.

Lesion-Specific Metabolic Alterations in Relapsing-Remitting Multiple Sclerosis Via 7 T Magnetic Resonance Spectroscopic Imaging

BACKGROUND

Magnetic resonance spectroscopic imaging (MRSI) of the brain enables in vivo assessment of metabolic alterations in multiple sclerosis (MS). This provides complementary insights into lesion pathology that cannot be obtained via T1- and T2-weighted conventional magnetic resonance imaging (cMRI).

PURPOSE

The aims of this study were to assess focal metabolic alterations inside and at the periphery of lesions that are visible or invisible on cMRI, and to correlate their metabolic changes with T1 hypointensity and the distance of lesions to cortical gray matter (GM).

METHODS

A 7 T MRSI was performed on 51 patients with relapsing-remitting MS (30 female/21 male; mean age, 35.4 ± 9.9 years). Mean metabolic ratios were calculated for segmented regions of interest (ROIs) of normal-appearing white matter, white matter lesions, and focal regions of increased mIns/tNAA invisible on cMRI. A subgroup analysis was performed after subdividing based on T1 relaxation and distance to cortical GM. Metabolite ratios were correlated with T1 and compared between different layers around cMRI-visible lesions.

RESULTS

Focal regions of, on average, 2.8-fold higher mIns/tNAA than surrounding normal-appearing white matter and with an appearance similar to that of MS lesions were found, which were not visible on cMRI (ie, ~4% of metabolic hotspots). T1 relaxation was positively correlated with mIns/tNAA ( P ≤ 0.01), and negatively with tNAA/tCr ( P ≤ 0.01) and tCho/tCr ( P ≤ 0.01). mIns/tCr was increased outside lesions, whereas tNAA/tCr distributions resembled macroscopic tissue damage inside the lesions. mIns/tCr was -21% lower for lesions closer to cortical GM ( P ≤ 0.05).

CONCLUSIONS

7 T MRSI allows in vivo visualization of focal MS pathology not visible on cMRI and the assessment of metabolite levels in the lesion center, in the active lesion periphery and in cortical lesions. This demonstrated the potential of MRSI to image mIns as an early biomarker in lesion development.

Artifact suppression in readout-segmented consistent K-t space EPSI (RS-COKE) for fast 1 H spectroscopic imaging at 7 T

PURPOSE

Fast proton (¹ H) MRSI is an important diagnostic tool for clinical investigations, providing metabolic and spatial information. MRSI at 7 T benefits from increased SNR and improved separation of peaks but requires larger spectral widths. RS-COKE (Readout-Segmented Consistent K-t space Epsi) is an echo planar spectroscopic imaging (Epsi) variant capable to support the spectral width required for human brain metabolites spectra at 7 T. However, mismatches between readout segments lead to artifacts, particularly when subcutaneous lipid signals are not suppressed. In this study, these mismatches and their effects are analyzed and reduced.

METHODS

The following corrections to the data were performed: i) frequency-dependent phase corrections; ii) k-space trajectory corrections, derived from short reference scans; and iii) smoothing of data at segment transitions to mitigate the effect of residual mismatches. The improvement was evaluated by performing single-slice RS-COKE on a head-shaped phantom with a "lipid" layer and healthy subjects, using varying resolutions and durations ranging from 4.1 × 4.7 × 15 mm³ in 5:46 min to 3.1 × 3.3 × 15 mm³ in 13:07 min.

RESULTS

Artifacts arising from the readout-segmented acquisition were substantially reduced, thus providing high-quality spectroscopic imaging in phantom and human scans. LCModel fitting of the human data resulted in a relative Cramer-Rao lower bounds within 6% for NAA, Cr, and Cho images in the majority of the voxels.

CONCLUSION

Using the new reference scans and reconstruction steps, RS-COKE was able to deliver fast ¹ H MRSI at 7 T, overcoming the spectral width limitation of standard EPSI at this field strength.

Test-retest reproducibility of human brain multi-slice 1 H FID-MRSI data at 9.4T after optimization of lipid regularization, macromolecular model, and spline baseline stiffness

PURPOSE

This study analyzes the effects of retrospective lipid suppression, a simulated macromolecular prior knowledge and different spline baseline stiffness values on 9.4T multi-slice proton FID-MRSI data spanning the whole cerebrum of human brain and the reproducibility of respective metabolite ratio to total creatine (/tCr) maps for 10 brain metabolites.

METHODS

Measurements were performed twice on 5 volunteers using a short TR and TE FID MRSI 2D sequence at 9.4T. The effects of retrospective lipid L2-regularization, macromolecular spectrum and different LCModel baseline flexibilities on SNR, FWHM, fitting residual, Cramér-Rao lower bound, and metabolite ratio maps were investigated. Intra-subject, inter-session coefficient of variation and the test-retest reproducibility of the mean metabolite ratios (/tCr) of each slice was calculated.

RESULTS

Transversal, sagittal, and coronal slices of many metabolite ratio maps correspond to the anatomically expected concentration relations in gray and white matter for the majority of the cerebrum when using a flexible baseline in LCModel fit. Results from the second measurements of the same subjects show that slice positioning and data quality correlate significantly to the first measurement. L2-regularization provided effective suppression of lipid-artifacts, but should be avoided if no artifacts are detected.

CONCLUSION

Reproducible concentration ratio maps (/tCr) for 4 metabolites (total choline, N-acetylaspartate, glutamate, and myoinositol) spanning the majority of the cerebrum and 6 metabolites (N-acetylaspartylglutamate, γ-aminobutyric acid, glutathione, taurine, glutamine, and aspartate) covering 32 mm in the upper part of the brain were acquired at 9.4T using multi-slice FID MRSI with retrospective lipid suppression, a macromolecular spectrum and a flexible LCModel baseline.

Quantitative T1-relaxation corrected metabolite mapping of 12 metabolites in the human brain at 9.4 T

Magnetic resonance spectroscopic imaging (MRSI) is a non-invasive imaging modality that enables observation of metabolites. Applications of MRSI for neuroimaging have shown promise for monitoring and detecting various diseases. This study builds off previously developed techniques of short TR, ¹H FID MRSI by correcting for T₁-weighting of the metabolites and utilizing an internal water reference to produce quantitative (mmol kg^-1) metabolite maps. This work reports and shows quantitative metabolite maps for 12 metabolites for a single slice. Voxel-specific T₁-corrections for water are common in MRSI studies; however, most studies use either averaged T₁-relaxation times to correct for T₁-weighting of metabolites or omit this correction step entirely. This work employs the use of voxel-specific T₁-corrections for metabolites in addition to water. Utilizing averaged T₁-relaxation times for metabolites can bias metabolite maps for metabolites that have strong differences between T₁-relaxation for GM and WM (i.e. Glu). This work systematically compares quantitative metabolite maps to single voxel quantitative results and qualitatively compares metabolite maps to previous works.

Short symmetric and highly selective asymmetric first and second order gradient modulated offset independent adiabaticity (GOIA) pulses for applications in clinical MRS and MRSI

Gradient modulated RF pulses, especially gradient offset independent adiabaticity (GOIA) pulses, are increasingly gaining attention for high field clinical magnetic resonance spectroscopy and spectroscopic imaging (MRS/MRSI) due to the lower peak B₁ amplitude and associated power demands achievable relative to its non-modulated adiabatic full passage counterparts. In this work we describe the development of two GOIA RF pulses: 1) A power efficient, 3.0 ms wideband uniform rate with smooth truncation (WURST) modulated RF pulse with 15 kHz bandwidth compatible with a clinically feasible peak B₁ amplitude of 0.87 kHz (or 20 µT), and 2) A highly selective asymmetric 6.66 ms RF pulse with 20 kHz bandwidth designed to achieve a single-sided, fractional transition width of only 1.7%. Effects of potential asynchrony between RF and gradient-modulated (GM) waveforms for 3 ms GOIA-WURST RF pulses was evaluated by simulation and experimentally. Results demonstrate that a 20+ µs asynchrony between RF and GM functions substantially degrades inversion performance when using large RF offsets to achieve translation. A projection-based method is presented that allows a quick calibration of RF and GM asynchrony on pre-clinical/clinical MR systems. The asymmetric GOIA pulse was implemented within a multi-pulse OVS sequence to achieve power efficient, highly-selective, and B₁ and T₁-independent signal suppression for extracranial lipid suppression. The developed GOIA pulses were utilized with linear gradient modulation (X, Y, Z gradient fields), and with second-order-field modulations (Z2, X2Y2 gradient fields) to provide elliptically-shaped regions-of-interest for MRS and MRSI acquisitions. Both described GOIA-RF pulses have substantial clinical value; specifically, the 3.0 ms GOIA-WURST pulse is beneficial to realize short TE sLASER localized proton MRS/MRSI sequences, and the asymmetric GOIA RF pulse has applications in highly selective outer volume signal suppression to allow interrogation of tissue proximal to extracranial lipids with full-intensity.

The Impact of Theta-Burst Stimulation on Cortical GABA and Glutamate in Treatment-Resistant Depression: A Surface-Based MRSI Analysis Approach

Background: Theta burst stimulation (TBS) belongs to one of the biological antidepressant treatment options. When applied bilaterally, excitatory intermittent TBS (iTBS) is commonly targeted to the left and inhibitory continuous TBS (cTBS) to the right dorsolateral prefrontal cortex. TBS was shown to influence neurotransmitter systems, while iTBS is thought to interfere with glutamatergic circuits and cTBS to mediate GABAergic neurotransmission.

Objectives: We aimed to expand insights into the therapeutic effects of TBS on the GABAergic and glutamatergic system utilizing 3D-multivoxel magnetic resonance spectroscopy imaging (MRSI) in combination with a novel surface-based MRSI analysis approach to investigate changes of cortical neurotransmitter levels in patients with treatment-resistant depression (TRD).

Methods: Twelve TRD patients (five females, mean age ± SD = 35 ± 11 years) completed paired MRSI measurements, using a GABA-edited 3D-multivoxel MEGA-LASER sequence, before and after 3 weeks of bilateral TBS treatment. Changes in cortical distributions of GABA+/tNAA (GABA+macromolecules relative to total N-acetylaspartate) and Glx/tNAA (Glx = mixed signal of glutamate and glutamine), were investigated in a surface-based region-of-interest (ROI) analysis approach.

Results: ANCOVAs revealed a significant increase in Glx/tNAA ratios in the left caudal middle frontal area (p_corr. = 0.046, F = 13.292), an area targeted by iTBS treatment. Whereas, contralateral treatment with cTBS evoked no alterations in glutamate or GABA concentrations.

Conclusion: This study demonstrates surface-based adaptions in the stimulation area to the glutamate metabolism after excitatory iTBS but not after cTBS, using a novel surface-based analysis of 3D-MRSI data. The reported impact of facilitatory iTBS on glutamatergic neurotransmission provides further insight into the neurobiological effects of TBS in TRD.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

7T HR FID-MRSI Compared to Amino Acid PET: Glutamine and Glycine as Promising Biomarkers in Brain Tumors

Simple Summary

Magnetic resonance spectroscopic imaging is an imaging method that can map the distribution of multiple biochemicals in the human brain in one scan. Using stronger magnetic fields, such as 7 Tesla, allows for higher resolution images and more biochemical maps. To test these results, we compared it to positron emission tomography, the established clinical standard for metabolic imaging. This comparison mainly looked at the overlap between regions with increased signal between both methods. We found that the molecules glutamine and glycine, only mappable at 7 Tesla, corresponded better to positron emission tomography than the commonly used choline.

Abstract

(1) Background: Recent developments in 7T magnetic resonance spectroscopic imaging (MRSI) made the acquisition of high-resolution metabolic images in clinically feasible measurement times possible. The amino acids glutamine (Gln) and glycine (Gly) were identified as potential neuro-oncological markers of importance. For the first time, we compared 7T MRSI to amino acid PET in a cohort of glioma patients. (2) Methods: In 24 patients, we co-registered 7T MRSI and routine PET and compared hotspot volumes of interest (VOI). We evaluated dice similarity coefficients (DSC), volume, center of intensity distance (CoI), median and threshold values for VOIs of PET and ratios of total choline (tCho), Gln, Gly, myo-inositol (Ins) to total N-acetylaspartate (tNAA) or total creatine (tCr). (3) Results: We found that Gln and Gly ratios generally resulted in a higher correspondence to PET than tCho. Using cutoffs of 1.6-times median values of a control region, DSCs to PET were 0.53 ± 0.36 for tCho/tNAA, 0.66 ± 0.40 for Gln/tNAA, 0.57 ± 0.36 for Gly/tNAA, and 0.38 ± 0.31 for Ins/tNAA. (4) Conclusions: Our 7T MRSI data corresponded better to PET than previous studies at lower fields. Our results for Gln and Gly highlight the importance of future research (e.g., using Gln PET tracers) into the role of both amino acids.

Inter-subject stability and regional concentration estimates of 3D-FID-MRSI in the human brain at 7 T

PURPOSE

Recently, a 3D-concentric ring trajectory (CRT)-based free induction decay (FID)-MRSI sequence was introduced for fast high-resolution metabolic imaging at 7 T. This technique provides metabolic ratio maps of almost the entire brain within clinically feasible scan times, but its robustness has not yet been thoroughly investigated. Therefore, we have assessed quantitative concentration estimates and their variability in healthy volunteers using this approach.

METHODS

We acquired whole-brain 3D-CRT-FID-MRSI at 7 T in 15 min with 3.4 mm nominal isometric resolution in 24 volunteers (12 male, 12 female, mean age 27 ± 6 years). Concentration estimate maps were calculated for 15 metabolites using internal water referencing and evaluated in 55 different regions of interest (ROIs) in the brain. Data quality, mean metabolite concentrations, and their inter-subject coefficients of variation (CVs) were compared for all ROIs.

RESULTS

Of 24 datasets, one was excluded due to motion artifacts. The concentrations of total choline, total creatine, glutamate, myo-inositol, and N-acetylaspartate in 44 regions were estimated within quality thresholds. Inter-subject CVs (mean over 44 ROIs/minimum/maximum) were 9%/5%/19% for total choline, 10%/6%/20% for total creatine, 11%/7%/24% for glutamate, 10%/6%/19% for myo-inositol, and 9%/6%/19% for N-acetylaspartate.

DISCUSSION

We defined the performance of 3D-CRT-based FID-MRSI for metabolite concentration estimate mapping, showing which metabolites could be robustly quantified in which ROIs with which inter-subject CVs expected. However, the basal brain regions and lesser-signal metabolites in particular remain as a challenge due susceptibility effects from the proximity to nasal and auditory cavities. Further improvement in quantification and the mitigation of B₀ /B₁ -field inhomogeneities will be necessary to achieve reliable whole-brain coverage.

Achieving high-resolution 1H-MRSI of the human brain with compressed-sensing and low-rank reconstruction at 7 Tesla

Low sensitivity MR techniques such as magnetic resonance spectroscopic imaging (MRSI) greatly benefit from the gain in signal-to-noise provided by ultra-high field MR. High-resolution and whole-slab brain MRSI remains however very challenging due to lengthy acquisition, low signal, lipid contamination and field inhomogeneity. In this study, we propose an acquisition-reconstruction scheme that combines 1H free-induction-decay (FID)-MRSI sequence, short TR acquisition, compressed sensing acceleration and low-rank modeling with total-generalized-variation constraint to achieve metabolite imaging in two and three dimensions at 7 Tesla. The resulting images and volumes reveal highly detailed distributions that are specific to each metabolite and follow the underlying brain anatomy. The MRSI method was validated in a high-resolution phantom containing fine metabolite structures, and in five healthy volunteers. This new application of compressed sensing acceleration paves the way for high-resolution MRSI in clinical setting with acquisition times of 5 min for 2D MRSI at 2.5 mm and of 20 min for 3D MRSI at 3.3 mm isotropic.

The cross-sectional interplay between neurochemical profile and brain connectivity

Neurochemical profile and brain connectivity are both critical aspects of brain function. However, our knowledge of their interplay across development is currently poor. We combined single-voxel magnetic resonance spectroscopy and resting functional magnetic resonance imaging in a cross-sectional sample spanning from childhood to adulthood which was reassessed in ~1.5 years (N = 293). We revealed the developmental trajectories of 20 neurochemicals in two key developmental brain regions (the intraparietal sulcus, IPS, and the middle frontal gyrus, MFG). We found that certain neurochemicals exhibited similar developmental trajectories across the two regions, while other trajectories were region-specific. Crucially, we mapped the connectivity of the brain regions IPS and MFG to the rest of the brain across development as a function of regional glutamate and GABA concentration. We demonstrated that glutamate concentration within the IPS is modulated by age in explaining IPS connectivity with frontal, temporal and parietal regions. In mature participants, higher glutamate within the IPS was related to more negative connectivity while the opposite pattern was found for younger participants. Our findings offer specific developmental insights on the interplay between the brain's resting activity and the glutamatergic system both of which are crucial for regulating normal functioning and are dysregulated in several clinical conditions.

Contribution of macromolecules to brain 1 H MR spectra: Experts' consensus recommendations

Proton MR spectra of the brain, especially those measured at short and intermediate echo times, contain signals from mobile macromolecules (MM). A description of the main MM is provided in this consensus paper. These broad peaks of MM underlie the narrower peaks of metabolites and often complicate their quantification but they also may have potential importance as biomarkers in specific diseases. Thus, separation of broad MM signals from low molecular weight metabolites enables accurate determination of metabolite concentrations and is of primary interest in many studies. Other studies attempt to understand the origin of the MM spectrum, to decompose it into individual spectral regions or peaks and to use the components of the MM spectrum as markers of various physiological or pathological conditions in biomedical research or clinical practice. The aim of this consensus paper is to provide an overview and some recommendations on how to handle the MM signals in different types of studies together with a list of open issues in the field, which are all summarized at the end of the paper.

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.

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.

Separation of Metabolites and Macromolecules for Short-TE 1H-MRSI Using Learned Component-Specific Representations

Short-echo-time (TE) proton magnetic resonance spectroscopic imaging (MRSI) allows for simultaneously mapping a number of molecules in the brain, and has been recognized as an important tool for studying in vivo biochemistry in various neuroscience and disease applications. However, separation of the metabolite and macromolecule (MM) signals present in the short-TE data with significant spectral overlaps remains a major technical challenge. This work introduces a new approach to solve this problem by integrating imaging physics and representation learning. Specifically, a mixed unsupervised and supervised learning-based strategy was developed to learn the metabolite and MM-specific low-dimensional representations using deep autoencoders. A constrained reconstruction formulation is proposed to integrate the MRSI spatiospectral encoding model and the learned representations as effective constraints for signal separation. An efficient algorithm was developed to solve the resulting optimization problem with provable convergence. Simulation and experimental results have been obtained to demonstrate the component-specific representation power of the learned models and the capability of the proposed method in separating metabolite and MM signals for practical short-TE [Formula: see text]-MRSI data.

Lipid-suppressed and tissue-fraction corrected metabolic distributions in human central brain structures using 2D 1 H magnetic resonance spectroscopic imaging at 7 T

INTRODUCTION

Magnetic resonance spectroscopic imaging (MRSI) has the potential to add a layer of understanding of the neurobiological mechanisms underlying brain diseases, disease progression, and treatment efficacy. Limitations related to metabolite fitting of low signal-to-noise ratios data, signal variations due to partial-volume effects, acquisition and extracranial lipid artifacts, along with clinically relevant aspects such as scan time constraints, are among the challenges associated with in vivo MRSI.

METHODS

The aim of this work was to address some of these factors and to develop an acquisition, reconstruction, and postprocessing pipeline to derive lipid-suppressed metabolite values of central brain structures based on free-induction decay measurements made using a 7 T MR scanner. Anatomical images were used to perform high-resolution (1 mm³ ) partial-volume correction to account for gray matter, white matter (WM), and cerebral-spinal fluid signal contributions. Implementation of automatic quality control thresholds and normalization of metabolic maps from 23 subjects to the Montreal Neurological Institute (MNI) standard atlas facilitated the creation of high-resolution average metabolite maps of several clinically relevant metabolites in central brain regions, while accounting for macromolecular distributions. Partial-volume correction improved the delineation of deep brain nuclei. We report average metabolite values including glutamate + glutamine (Glx), glycerophosphocholine, choline and phosphocholine (tCho), (phospo)creatine, myo-inositol and glycine (mI-Gly), glutathione, N-acetyl-aspartyl glutamate(and glutamine), and N-acetyl-aspartate in the basal ganglia, central WM (thalamic radiation, corpus callosum) as well as insular cortex and intracalcarine sulcus.

CONCLUSION

MNI-registered average metabolite maps facilitate group-based analysis, thus offering the possibility to mitigate uncertainty in variable MRSI data.

Flexible MEGA editing scheme with asymmetric adiabatic pulses applied for T 2 measurement of lactate in human brain

PURPOSE

A flexible MEGA editing scheme which decouples the editing efficiency from TE is proposed and the utility of asymmetric adiabatic pulses for this new technique is explored. It is demonstrated that the method enables robust T 2 measurement of lactate in healthy human brain.

METHODS

The proposed variation of the MEGA scheme applies editing pulses in both acquired spectra, ensuring that the difference in J-evolution of the target resonance leads to maximal signal yield in the difference spectrum for arbitrary TE. A MEGA-sLASER sequence is augmented with asymmetric adiabatic editing pulses for enhanced flexibility and immunity to B 1 + miscalibration and inhomogeneities. The technique is validated and optimized for flexible lactate editing via a simple analytical model, numerical simulations and in vitro experiments. The T 2 relaxation constant of lactate is determined in vivo via multiple-TE measurements with the proposed method and a dedicated postprocessing and quantification approach.

RESULTS

Asymmetric adiabatic editing pulses improve robustness and facilitate efficient J-editing in sequences or protocols with strong timing constraints. Single voxel measurements using the proposed MEGA scheme in the occipital cortex of six healthy subjects yield a relaxation constant of T 2 = 171 ± 19  ms for the methyl resonance of lactate at a field strength of 3T.

CONCLUSIONS

The proposed MEGA editing scheme allows for novel kinds of J-editing experiments and promises to be an asset to robust T 2 measurement of lactate and potentially other J-coupled metabolites in vivo.

High-resolution metabolic imaging of high-grade gliomas using 7T-CRT-FID-MRSI

OBJECTIVES

Successful neurosurgical intervention in gliomas depends on the precision of the preoperative definition of the tumor and its margins since a safe maximum resection translates into a better patient outcome. Metabolic high-resolution imaging might result in improved presurgical tumor characterization, and thus optimized glioma resection. To this end, we validated the performance of a fast high-resolution whole-brain 3D-magnetic resonance spectroscopic imaging (MRSI) method at 7T in a patient cohort of 23 high-grade gliomas (HGG).

MATERIALS AND METHODS

We preoperatively measured 23 patients with histologically verified HGGs (17 male, 8 female, age 53 ± 15) with an MRSI sequence based on concentric ring trajectories with a 64 × 64 × 39 measurement matrix, and a 3.4 × 3.4 × 3.4 mm3 nominal voxel volume in 15 min. Quantification used a basis-set of 17 components including N-acetyl-aspartate (NAA), total choline (tCho), total creatine (tCr), glutamate (Glu), glutamine (Gln), glycine (Gly) and 2-hydroxyglutarate (2HG). The resultant metabolic images were evaluated for their reliability as well as their quality and compared to spatially segmented tumor regions-of-interest (necrosis, contrast-enhanced, non-contrast enhanced + edema, peritumoral) based on clinical data and also compared to histopathology (e.g., grade, IDH-status).

RESULTS

Eighteen of the patient measurements were considered usable. In these patients, ten metabolites were quantified with acceptable quality. Gln, Gly, and tCho were increased and NAA and tCr decreased in nearly all tumor regions, with other metabolites such as serine, showing mixed trends. Overall, there was a reliable characterization of metabolic tumor areas. We also found heterogeneity in the metabolic images often continued into the peritumoral region. While 2HG could not be satisfyingly quantified, we found an increase of Glu in the contrast-enhancing region of IDH-wildtype HGGs and a decrease of Glu in IDH1-mutant HGGs.

CONCLUSIONS

We successfully demonstrated high-resolution 7T 3D-MRSI in HGG patients, showing metabolic differences between tumor regions and peritumoral tissue for multiple metabolites. Increases of tCho, Gln (related to tumor metabolism), Gly (related to tumor proliferation), as well as decreases in NAA, tCr, and others, corresponded very well to clinical tumor segmentation, but were more heterogeneous and often extended into the peritumoral region.

Clinical High-Resolution 3D-MR Spectroscopic Imaging of the Human Brain at 7 T

OBJECTIVES

Available clinical magnetic resonance spectroscopic imaging (MRSI) sequences are hampered by long scan times, low spatial resolution, strong field inhomogeneities, limited volume coverage, and low signal-to-noise ratio. High-resolution, whole-brain mapping of more metabolites than just N-acetylaspartate, choline, and creatine within clinically attractive scan times is urgently needed for clinical applications. The aim is therefore to develop a free induction decay (FID) MRSI sequence with rapid concentric ring trajectory (CRT) encoding for 7 T and demonstrate its clinical feasibility for mapping the whole cerebrum of healthy volunteers and patients.

MATERIALS AND METHODS

Institutional review board approval and written informed consent were obtained. Time-efficient, 3-dimensional encoding of an ellipsoidal k-space by in-plane CRT and through-plane phase encoding was integrated into an FID-MRSI sequence. To reduce scan times further, repetition times were shortened, and variable temporal interleaves were applied. Measurements with different matrix sizes were performed to validate the CRT encoding in a resolution phantom. One multiple sclerosis patient, 1 glioma patient, and 6 healthy volunteers were prospectively measured. For the healthy volunteers, brain segmentation was performed to quantify median metabolic ratios, Cramér-Rao lower bounds (CRLBs), signal-to-noise ratios, linewidths, and brain coverage among all measured matrix sizes ranging from a 32 × 32 × 31 matrix with 6.9 × 6.9 × 4.2 mm nominal voxel size acquired in ~3 minutes to an 80 × 80 × 47 matrix with 2.7 × 2.7 × 2.7 mm nominal voxel size in ~15 minutes for different brain regions.

RESULTS

Phantom structures with diameters down to 3 to 4 mm were visible. In vivo MRSI provided high spectral quality (median signal-to-noise ratios, >6.3 and linewidths, <0.082 ppm) and fitting quality. Cramér-Rao lower bounds were ranging from less than 22% for glutamine (highest CRLB in subcortical gray matter) to less than 9.5% for N-acetylaspartate for the 80 × 80 × 47 matrix (highest CRLB in the temporal lobe). This enabled reliable mapping of up to 8 metabolites (N-acetylaspartate, N-acetylaspartyl glutamate, total creatine, glutamine, glutamate, total choline, myo-inositol, glycine) and macromolecules for all resolutions. Coverage of the whole cerebrum allowed visualization of the full extent of diffuse and local multiple sclerosis-related neurochemical changes (eg, up to 100% increased myo-inositol). Three-dimensional brain tumor metabolic maps provided valuable information beyond that of single-slice MRSI, with up to 200% higher choline, up to 100% increased glutamine, and increased glycine in tumor tissue.

CONCLUSIONS

Seven Tesla FID-MRSI with time-efficient CRT readouts offers clinically attractive acquisition protocols tailored either for speed or for the investigation of small pathologic details and low-abundant metabolites. This can complement clinical MR studies of various brain disorders. Significant metabolic anomalies were demonstrated in a multiple sclerosis and a glioma patient for myo-inositol, glutamine, total choline, glycine, and N-acetylaspartate concentrations.

Effects of different macromolecular models on reproducibility of FID-MRSI at 7T

PURPOSE

A properly characterized macromolecular (MM) contribution is essential for accurate metabolite quantification in FID-MRSI. MM information can be included into the fitting model as a single component or parameterized and included over several individual MM resonances, which adds flexibility when pathologic changes are present but is prone to potential overfitting. This study investigates the effects of different MM models on MRSI reproducibility.

METHODS

Clinically feasible, high-resolution FID-MRSI data were collected in ~5 min at 7 Tesla from 10 healthy volunteers and quantified via LCModel (version 6.3) with 3 basis sets, each with a different approach for how the MM signal was handled: averaged measured whole spectrum (full MM), 9 parameterized components (param MM) with soft constraints to avoid overparameterization, or without any MM information included in the fitting prior knowledge. The test-retest reproducibility of MRSI scans was assessed voxel-wise using metabolite coefficients of variation and intraclass correlation coefficients and compared between the basis sets. Correlations of concentration estimates were investigated for the param MM fitting model.

RESULTS

The full MM model provided the most reproducible quantification of total NAA, total Cho, myo-inositol, and glutamate + glutamine ratios to total Cr (coefficients of variations ≤ 8%, intraclass correlation coefficients ≥ 0.76). Using the param MM model resulted in slightly lower reproducibility (up to +3% higher coefficients of variations, up to -0.1 decreased intraclass correlation coefficients). The quantification of the parameterized macromolecules did not affect quantification of the overlapping metabolites.

CONCLUSION

Clinically feasible FID-MRSI with an experimentally acquired MM spectrum included in prior knowledge provides highly reproducible quantification for the most common neurometabolites in healthy volunteers. Parameterization of the MM spectrum may be preferred as a compromise between quantification accuracy and reproducibility when the MM content is expected to be pathologically altered.

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.

7 T Magnetic Resonance Spectroscopic Imaging in Multiple Sclerosis: How Does Spatial Resolution Affect the Detectability of Metabolic Changes in Brain Lesions?

OBJECTIVES

The aim of this study was to assess the utility of increased spatial resolution of magnetic resonance spectroscopic imaging (MRSI) at 7 T for the detection of neurochemical changes in multiple sclerosis (MS)-related brain lesions.

MATERIALS AND METHODS

This prospective, institutional review board-approved study was performed in 20 relapsing-remitting MS patients (9 women/11 men; mean age ± standard deviation, 30.8 ± 7.7 years) after receiving written informed consent. Metabolic patterns in MS lesions were compared at 3 different spatial resolutions of free induction decay MRSI with implemented parallel imaging acceleration: 2.2 × 2.2 × 8 mm; 3.4 × 3.4 × 8 mm; and 6.8 × 6.8 × 8 mm voxel volumes, that is, matrix sizes of 100 × 100, 64 × 64, and 32 × 32, respectively. The quality of data was assessed by signal-to-noise ratio and Cramér-Rao lower bounds. Statistical analysis was performed using Wilcoxon signed-rank tests with correction for multiple testing.

RESULTS

Seventy-seven T2-hyperintense MS lesions were investigated (median volume, 155.7 mm; range, 10.8-747.0 mm). The mean metabolic ratios in lesions differed significantly between the 3 MRSI resolutions (ie, 100 × 100 vs 64 × 64, 100 × 100 vs 32 × 32, and 64 × 64 vs 32 × 32; P < 0.001). With the ultra-high resolution (100 × 100), we obtained 40% to 80% higher mean metabolic ratios and 100% to 150% increase in maximum metabolic ratios in the MS lesions compared with the lowest resolution (32 × 32), while maintaining good spectral quality (signal-to-noise ratio >12, Cramér-Rao lower bounds <20%) and measurement time of 6 minutes. There were 83% of MS lesions that showed increased myo-inositol/N-acetylaspartate with the 100 × 100 resolution, but only 66% were distinguishable with the 64 × 64 resolution and 35% with the 32 × 32 resolution.

CONCLUSIONS

Ultra-high-resolution MRSI (~2 × 2 × 8 mm voxel volume) can detect metabolic alterations in MS, which cannot be recognized by conventional MRSI resolutions, within clinically acceptable time.

Fast high-resolution brain metabolite mapping on a clinical 3T MRI by accelerated 1 H-FID-MRSI and low-rank constrained reconstruction

PURPOSE

Epitomizing the advantages of ultra short echo time and no chemical shift displacement error, high-resolution-free induction decay magnetic resonance spectroscopic imaging (FID-MRSI) sequences have proven to be highly effective in providing unbiased characterizations of metabolite distributions. However, its merits are often overshadowed in high-resolution settings by reduced signal-to-noise ratios resulting from the smaller voxel volumes procured by extensive phase encoding and the related acquisition times.

METHODS

To address these limitations, we here propose an acquisition and reconstruction scheme that offers both implicit dataset denoising and acquisition acceleration. Specifically, a slice selective high-resolution FID-MRSI sequence was implemented. Spectroscopic datasets were processed to remove fat contamination, and then reconstructed using a total generalized variation (TGV) regularized low-rank model. We further measured reconstruction performance for random undersampled data to assess feasibility of a compressed-sensing SENSE acceleration scheme. Performance of the lipid suppression was assessed using an ad hoc phantom, while that of the low-rank TGV reconstruction model was benchmarked using simulated MRSI data. To assess real-world performance, 2D FID-MRSI acquisitions of the brain in healthy volunteers were reconstructed using the proposed framework.

RESULTS

Results from the phantom and simulated data demonstrate that skull lipid contamination is effectively removed and that data reconstruction quality is improved with the low-rank TGV model. Also, we demonstrated that the presented acquisition and reconstruction methods are compatible with a compressed-sensing SENSE acceleration scheme.

CONCLUSIONS

An original reconstruction pipeline for 2D 1 H-FID-MRSI datasets was presented that places high-resolution metabolite mapping on 3T MR scanners within clinically feasible limits.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Macromolecule-suppressed GABA measurements correlate more strongly with behavior than macromolecule-contaminated GABA+ measurements

The inhibitory neurotransmitter γ-aminobutyric acid (GABA) is known to be fundamental to the neuronal processes underlying visual orientation and vibrotactile frequency and amplitude discrimination. Previous studies have demonstrated that performance on visual and vibrotactile psychophysics tasks is associated with in vivo measurements of "GABA+" levels - a measure of GABA substantially contaminated by a macromolecular (MM) signal. Here, we establish that these prior findings are indeed driven by the GABA fraction of that signal. Edited magnetic resonance spectroscopy (MRS) was used to measure GABA with and without MM suppression in the sensorimotor (SM1) and occipital cortices in 14 healthy male adults. Volunteers also underwent psychophysical experiments to assess their performance on visual orientation discrimination and vibrotactile amplitude and frequency discrimination. We show that MM-suppressed GABA levels correlate more strongly with individual differences in vibrotactile (in the case of SM1 GABA; amplitude: r = -0.63, p = 0.03; frequency: r = -0.62, p = 0.02) and visual orientation (in the case of occipital GABA; r = -0.59, p = 0.05) discrimination thresholds than GABA levels contaminated by MM (vibrotactile amplitude: r = -0.36, p = 0.30; vibrotactile frequency: r = -0.53, p = 0.09; visual orientation: r = 0.21, p = 0.55). These findings further support the view that measurements of endogenous GABA acquired with edited MRS can usefully probe neurochemical-behavioral relationships in humans. Moreover, the more specific measurement of GABA used in this study provides increased statistical power to observe these regionally specific relationships.

Density-weighted concentric circle trajectories for high resolution brain magnetic resonance spectroscopic imaging at 7T

PURPOSE

Full-slice magnetic resonance spectroscopic imaging at ≥7 T is especially vulnerable to lipid contaminations arising from regions close to the skull. This contamination can be mitigated by improving the point spread function via higher spatial resolution sampling and k-space filtering, but this prolongs scan times and reduces the signal-to-noise ratio (SNR) efficiency. Currently applied parallel imaging methods accelerate magnetic resonance spectroscopic imaging scans at 7T, but increase lipid artifacts and lower SNR-efficiency further. In this study, we propose an SNR-efficient spatial-spectral sampling scheme using concentric circle echo planar trajectories (CONCEPT), which was adapted to intrinsically acquire a Hamming-weighted k-space, thus termed density-weighted-CONCEPT. This minimizes voxel bleeding, while preserving an optimal SNR.

THEORY AND METHODS

Trajectories were theoretically derived and verified in phantoms as well as in the human brain via measurements of five volunteers (single-slice, field-of-view 220 × 220 mm2 , matrix 64 × 64, scan time 6 min) with free induction decay magnetic resonance spectroscopic imaging. Density-weighted-CONCEPT was compared to (a) the originally proposed CONCEPT with equidistant circles (here termed e-CONCEPT), (b) elliptical phase-encoding, and (c) 5-fold Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration accelerated elliptical phase-encoding.

RESULTS

By intrinsically sampling a Hamming-weighted k-space, density-weighted-CONCEPT removed Gibbs-ringing artifacts and had in vivo +9.5%, +24.4%, and +39.7% higher SNR than e-CONCEPT, elliptical phase-encoding, and the Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration accelerated elliptical phase-encoding (all P < 0.05), respectively, which lead to improved metabolic maps.

CONCLUSION

Density-weighted-CONCEPT provides clinically attractive full-slice high-resolution magnetic resonance spectroscopic imaging with optimal SNR at 7T. Magn Reson Med 79:2874-2885, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Key clinical benefits of neuroimaging at 7T

The growing interest in ultra-high field MRI, with more than 35.000 MR examinations already performed at 7T, is related to improved clinical results with regard to morphological as well as functional and metabolic capabilities. Since the signal-to-noise ratio increases with the field strength of the MR scanner, the most evident application at 7T is to gain higher spatial resolution in the brain compared to 3T. Of specific clinical interest for neuro applications is the cerebral cortex at 7T, for the detection of changes in cortical structure, like the visualization of cortical microinfarcts and cortical plaques in Multiple Sclerosis. In imaging of the hippocampus, even subfields of the internal hippocampal anatomy and pathology may be visualized with excellent spatial resolution. Using Susceptibility Weighted Imaging, the plaque-vessel relationship and iron accumulations in Multiple Sclerosis can be visualized, which may provide a prognostic factor of disease. Vascular imaging is a highly promising field for 7T which is dealt with in a separate dedicated article in this special issue. The static and dynamic blood oxygenation level-dependent contrast also increases with the field strength, which significantly improves the accuracy of pre-surgical evaluation of vital brain areas before tumor removal. Improvement in acquisition and hardware technology have also resulted in an increasing number of MR spectroscopic imaging studies in patients at 7T. More recent parallel imaging and short-TR acquisition approaches have overcome the limitations of scan time and spatial resolution, thereby allowing imaging matrix sizes of up to 128×128. The benefits of these acquisition approaches for investigation of brain tumors and Multiple Sclerosis have been shown recently. Together, these possibilities demonstrate the feasibility and advantages of conducting routine diagnostic imaging and clinical research at 7T.

Altered macromolecular pattern and content in the aging human brain

The resonances originating from proteins underlie those of metabolites in brain ¹ H nuclear magnetic resonance (NMR) spectra. These resonances have different physical properties from those of metabolites, such as shorter T₁ and T₂ relaxation time constants. The age dependence of the macromolecular pattern and content in the human brain was investigated with a focus on adults over 66 years of age using ultrahigh-field in vivo magnetic resonance spectroscopy. Eighteen young and 23 cognitively normal older adults were studied at 7 T. Metabolite spectra were acquired in the occipital cortex and the posterior cingulate cortex with single-voxel stimulated echo acquisition mode (STEAM) spectroscopy in 14 young and 20 older adults. Macromolecular spectra were acquired in the occipital cortex using an inversion recovery STEAM sequence in four young and three older adults. The macromolecular pattern was apparent over the 0.5-4.5-ppm range in the inversion recovery spectra and the 0.5-2-ppm range in the metabolite spectra. Macromolecular content was quantified from metabolite spectra using LCModel and from inversion recovery spectra using integration. Age-associated differences in the macromolecular pattern were apparent via both types of spectra, with the largest difference observed for the 1.7- and 2-ppm macromolecular resonances. A higher macromolecular content was observed in the older adults for both brain regions. Age-specific macromolecular spectra are needed when comparing metabolite spectra from subjects of differing ages because of age-associated differences in macromolecular pattern. Age-associated pattern and content differences may provide information about the aging process.

Metabolic and Structural Imaging at 7 Tesla After Repetitive Mild Traumatic Brain Injury in Immature Rats

Mild traumatic brain injury (mTBI) in children is a common and serious public health problem. Traditional neuroimaging findings in children who sustain mTBI are often normal, putting them at risk for repeated mTBI (rmTBI). There is a need for more sensitive imaging techniques capable of detecting subtle neurophysiological alterations after injury. We examined neurochemical and white matter changes using diffusion tensor imaging of the whole brain and proton magnetic resonance spectroscopy of the hippocampi at 7 Tesla in 18-day-old male rats at 7 days after mTBI and rmTBI. Traumatic axonal injury was assessed by beta-amyloid precursor protein accumulation using immunohistochemistry. A significant decrease in fractional anisotropy and increase in axial and radial diffusivity were observed in several brain regions, especially in white matter regions, after a single mTBI versus sham and more prominently after rmTBI. In addition, we observed accumulation of beta-amyloid precursor protein in the external capsule after mTBI and rmTBI. mTBI and rmTBI reduced the N-acetylaspartate/creatine ratio (NAA/Cr) and increased the myoinositol/creatine ratio (Ins/Cr) versus sham. rmTBI exacerbated the reduction in NAA/Cr versus mTBI. The choline/creatine (Cho/Cr) and (lipid/Macro Molecule 1)/creatine (Lip/Cr) ratios were also decreased after rmTBI versus sham. Diffusion tensor imaging findings along with the decrease in Cho and Lip after rmTBI may reflect damage to axonal membrane. NAA and Ins are altered at 7 days after mTBI and rmTBI likely reflecting neuro-axonal damage and glial response, respectively. These findings may be relevant to understanding the extent of disability following mTBI and rmTBI in the immature brain and may identify possible therapeutic targets.