Metabolite‐cycled STEAM and semi‐LASER localization for MR spectroscopy of the human brain at 9.4T

High-field downfield MR spectroscopic imaging in the human brain

PURPOSE

To investigate the feasibility of downfield MR spectroscopic imaging (DF-MRSI) in the human brain at 7T.

METHODS

A 7T DF-MRSI pulse sequence was implemented based on the previously described methodology at 3T, with 3D phase-encoding,1 3 ‾ 3 1 ‾ $$ 1\overline{3}3\overline{1} $$ spectral-spatial excitation, and frequency selective refocusing. Data were pre-processed followed by analysis using the "LCModel" software package, and metabolite maps created from the LCModel results. Total scan time, including brain MRI and a water-reference MRSI, was 24 min. The sequence was tested in 10 normal volunteers. Estimated metabolite levels and uncertainty values (Cramer Rao lower bounds, CRLBs) for nine downfield peaks were compared between seven different brain regions, anterior cingulate cortex (ACC), centrum semiovale (CSO), corpus callosum (CC), cerebellar vermis (CV), dorsolateral prefrontal cortex (DLPFC), posterior cingulate cortex (PCC), and thalamus (Thal).

RESULTS

DF peaks were relatively uniformly distributed throughout the brain, with only a small number of peaks showing any significant regional variations. Most DF peaks had average CRLB<25% in most brain regions. Average SNR values were higher for the brain regions ACC and DLPFC (˜7 ± 0.95, mean ± SD) while in a range of 3.4-6.0 for other brain regions. Average linewidth (FWHM) values were greater than 35 Hz in the ACC, CV, and Thal, and 22 Hz in CC, CSO, DLPFC, and PCC.

CONCLUSION

High-field DF-MRSI is able to spatially map exchangeable protons in the human brain at high resolution and with near whole-brain coverage in acceptable scan times, and in the future may be used to study metabolism of brain tumors or other neuropathological disorders.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

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 .

Meta-analysis and open-source database for in vivo brain Magnetic Resonance spectroscopy in health and disease

Proton (1H) Magnetic Resonance Spectroscopy (MRS) is a non-invasive tool capable of quantifying brain metabolite concentrations in vivo. Prioritization of standardization and accessibility in the field has led to the development of universal pulse sequences, methodological consensus recommendations, and the development of open-source analysis software packages. One on-going challenge is methodological validation with ground-truth data. As ground-truths are rarely available for in vivo measurements, data simulations have become an important tool. The diverse literature of metabolite measurements has made it challenging to define ranges to be used within simulations. Especially for the development of deep learning and machine learning algorithms, simulations must be able to produce accurate spectra capturing all the nuances of in vivo data. Therefore, we sought to determine the physiological ranges and relaxation rates of brain metabolites which can be used both in data simulations and as reference estimates. Using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, we've identified relevant MRS research articles and created an open-source database containing methods, results, and other article information as a resource. Using this database, expectation values and ranges for metabolite concentrations and T2 relaxation times are established based upon a meta-analyses of healthy and diseased brains.

Meta-analysis and Open-source Database for In Vivo Brain Magnetic Resonance Spectroscopy in Health and Disease

Proton ( 1 H) Magnetic Resonance Spectroscopy (MRS) is a non-invasive tool capable of quantifying brain metabolite concentrations in vivo . Prioritization of standardization and accessibility in the field has led to the development of universal pulse sequences, methodological consensus recommendations, and the development of open-source analysis software packages. One on-going challenge is methodological validation with ground-truth data. As ground-truths are rarely available for in vivo measurements, data simulations have become an important tool. The diverse literature of metabolite measurements has made it challenging to define ranges to be used within simulations. Especially for the development of deep learning and machine learning algorithms, simulations must be able to produce accurate spectra capturing all the nuances of in vivo data. Therefore, we sought to determine the physiological ranges and relaxation rates of brain metabolites which can be used both in data simulations and as reference estimates. Using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, we've identified relevant MRS research articles and created an open-source database containing methods, results, and other article information as a resource. Using this database, expectation values and ranges for metabolite concentrations and T 2 relaxation times are established based upon a meta-analyses of healthy and diseased brains.

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 .

Measurement of glucose metabolism in the occipital lobe and frontal cortex after oral administration of [1-13C]glucose at 9.4 T

For the first time, labeling effects after oral intake of [1-13C]glucose are observed in the human brain with pure 1H detection at 9.4 T. Spectral time series were acquired using a short-TE 1H MRS MC-semiLASER (Metabolite Cycling semi Localization by Adiabatic SElective Refocusing) sequence in two voxels of 5.4 mL in the frontal cortex and the occipital lobe. High-quality time-courses of [4-13C]glutamate, [4-13C]glutamine, [3-13C]glutamate + glutamine, [2-13C] glutamate+glutamine and [3-13C]aspartate for individual volunteers and additionally, group-averaged time-courses of labeled and non-labeled brain glucose could be obtained. Using a one-compartment model, mean metabolic rates were calculated for each voxel position: The mean rate of the TCA-cycle (Vtca) value was determined to be 1.36 and 0.93 μmol min^-1 g^-1, the mean rate of glutamine synthesis (Vgln) was calculated to be 0.23 and 0.45 μmol min^-1 g^-1, the mean exchange rate between cytosolic amino acids and mitochondrial Krebs cycle intermediates (Vx) rate was found to be 0.57 and 1.21 μmol min^-1 g^-1 for the occipital lobe and the frontal cortex, respectively. These values were in agreement with previously reported data. Altogether, it can be shown that this most simple technique combining oral administration of [1-13C]Glc with pure 1H MRS acquisition is suitable to measure metabolic rates.

Simultaneous detection of metabolite concentration changes, water BOLD signal and pH changes during visual stimulation in the human brain at 9.4T

This study presents a method to directly link metabolite concentration changes and BOLD response in the human brain during visual stimulation by measuring the water and metabolite signals simultaneously. Therefore, the metabolite-cycling (MC) non-water suppressed semiLASER localization technique was optimized for functional ¹H MRS in the human brain at 9.4 T. Data of 13 volunteers were acquired during a 26:40 min visual stimulation block-design paradigm. Activation-induced BOLD signal was observed in the MC water signal as well as in the NAA-CH₃ and tCr-CH₃ singlets. During stimulation, glutamate concentration increased 2.3 ± 2.0% to a new steady-state, while a continuous increase over the whole stimulation period could be observed in lactate with a mean increase of 35.6 ± 23.1%. These increases of Lac and Glu during brain activation confirm previous findings reported in literature. A positive correlation of the MC water BOLD signal with glutamate and lactate concentration changes was found. In addition, a pH decrease calculated from a change in the ratio of PCr to Cr was observed during brain activation, particularly at the onset of the stimulation.

Quantification of NAD + in human brain with 1 H MR spectroscopy at 3 T: Comparison of three localization techniques with different handling of water magnetization

Purpose

The detection of nicotinamide‐adenine‐dinucleotide (NAD⁺) is challenging using standard ¹H MR spectroscopy, because it is of low concentration and affected by polarization‐exchange with water. Therefore, this study compares three techniques to access NAD⁺ quantification at 3 T–one with and two without water presaturation.

Methods

A large brain volume in 10 healthy subjects was investigated with three techniques: semi‐LASER with water‐saturation (WS) (TE = 35 ms), semi‐LASER with metabolite‐cycling (MC) (TE = 35 ms), and the non‐water‐excitation (nWE) technique 2D ISIS‐localization with chemical‐shift‐selective excitation (2D I‐CSE) (TE = 10.2 ms). Spectra were quantified with optimized modeling in FiTAID.

Results

NAD⁺ could be well quantified in cohort‐average spectra with all techniques. Obtained apparent NAD⁺ tissue contents are all lower than expected from literature confirming restricted visibility by ¹H MRS. The estimated value from WS‐MRS (58 μM) was considerably lower than those obtained with non‐WS techniques (146 μM for MC‐semi‐LASER and 125 μM for 2D I‐CSE). The nWE technique with shortest TE gave largest NAD⁺ signals but suffered from overlap with large amide signals. MC‐semi‐LASER yielded best estimation precision as reflected in relative Cramer‐Rao bounds (14%, 21 μM/146 μM) and also best robustness as judged by the coefficient‐of‐variance over the cohort (11%, 10 μM/146 μM). The MR‐visibility turned out as 16% with WS and 41% with MC.

Conclusion

Three methods to assess NAD⁺ in human brain at 3 T have been compared. NAD⁺ could be detected with a visibility of ∼41% for the MC method. This may open a new window for the observation of pathological changes in the clinical research setting.

Respiratory-triggered quantitative MR spectroscopy of the human cervical spinal cord at 7 T

PURPOSE

Ultra-high field 1 H MR spectroscopy (MRS) is of great interest to help characterizing human spinal cord pathologies. However, very few studies have been reported so far in this small size structure at these fields due to challenging experimental difficulties caused by static and radiofrequency field heterogeneities, as well as physiological motion. In this work, in line with the recent developments proposed to strengthen spinal cord MRS feasibility at 7 T, a respiratory-triggered acquisition approach was optimized to compensate for dynamic B 0 field heterogeneities and to provide robust cervical spinal cord MRS data.

METHODS

A semi-LASER sequence was purposely used, and a dedicated raw data processing algorithm was developed to enhance MR spectral quality by discarding corrupted scans. To legitimate the choices done during the optimization stage, additional tests were carried out to determine the impact of breathing, voluntary motion, body mass index, and fitting algorithm. An in-house quantification tool was concomitantly designed for accurate estimation of the metabolite concentration ratios for choline, N-acetyl-aspartate (NAA), myo-inositol and glutathione. The method was tested on a cohort of 14 healthy volunteers.

RESULTS

Average water linewidth and NAA signal-to-noise ratio reached 0.04 ppm and 11.01, respectively. The group-average metabolic ratios were in good agreement with previous studies and showed intersession reproducibility variations below 30%.

CONCLUSION

The developed approach allows a rise of the acquired MRS signal quality and of the quantification robustness as compared to previous studies hence offering strengthened possibilities to probe the metabolism of degenerative and traumatic spinal cord pathologies.

Magnetic resonance spectroscopic imaging of downfield proton resonances in the human brain at 3 T

PURPOSE

To develop an MRSI technique capable of mapping downfield proton resonances in the human brain.

METHODS

A spectral-spatial excitation and frequency-selective refocusing scheme, in combination with 2D phase encoding, was developed for mapping of downfield resonances without any perturbation of the water magnetization. An alternative scheme using spectral-spatial refocusing was also investigated for simultaneous detection of both downfield and upfield resonances. The method was tested in 5 healthy human volunteers.

RESULTS

Downfield metabolite maps with a nominal spatial resolution of 1.5 cm³ were recorded at 3 T in a scan time of 12 minutes. Cramer-Rao lower bounds for nine different downfield peaks were 20% or less over a single supraventricular slice. Downfield spectral profiles were similar to those in the literature recorded previously using single-voxel localization methods. The same approach was also used for upfield MRSI, and simultaneous upfield and downfield acquisitions.

CONCLUSION

The developed MRSI pulse sequence was shown to be an efficient way of rapidly mapping downfield resonances in the human brain at 3 T, maximizing sensitivity through the relaxation enhancement effect. Because the MRSI approach is efficient in terms of data collection and can be readily implemented at short TE, somewhat higher spatial resolution can be achieved than has been reported in previous single-voxel downfield MRS studies. With this approach, nine downfield resonances could be mapped in a single slice for the first time using MRSI at 3 T.

Absolute choline tissue concentration mapping for prostate cancer localization and characterization using 3D 1 H MRSI without water-signal suppression

PURPOSE

Until now, ¹ H MRSI of the prostate has been performed with suppression of the large water signal to avoid distortions of metabolite signals. However, this signal can be used for absolute quantification and spectral corrections. We investigated the feasibility of water-unsuppressed MRSI in patients with prostate cancer for water signal-mediated spectral quality improvement and determination of absolute tissue levels of choline.

METHODS

Eight prostate cancer patients scheduled for radical prostatectomy underwent multi-parametric MRI at 3 T, including 3D water-unsuppressed semi-LASER MRSI. A postprocessing algorithm was developed to remove the water signal and its artifacts and use the extracted water signal as intravoxel reference for phase and frequency correction of metabolite signals and for absolute metabolite quantification.

RESULTS

Water-unsuppressed MRSI with dedicated postprocessing produced water signal and artifact-free MR spectra throughout the prostate. In all patients, the absolute choline tissue concentration was significantly higher in tumorous than in benign tissue areas (mean ± SD: 7.2 ± 1.4 vs 3.8 ± 0.7 mM), facilitating tumor localization by choline mapping. Tumor tissue levels of choline correlated better with the commonly used (choline + spermine + creatine)/citrate ratio (r = 0.78 ± 0.1) than that of citrate (r = 0.21 ± 0.06). The highest maximum choline concentrations occurred in high-risk cancer foci.

CONCLUSION

This report presents the first successful water-unsuppressed MRSI of the whole prostate. The water signal enabled amelioration of spectral quality and absolute metabolite quantification. In this way, choline tissue levels were identified as tumor biomarker. Choline mapping may serve as a tool in prostate cancer localization and risk scoring in multi-parametric MRI for diagnosis and biopsy procedures.

Improved prospective frequency correction for macromolecule-suppressed GABA editing with metabolite cycling at 3T

PURPOSE

To combine metabolite cycling with J-difference editing (MC MEGA) to allow for prospective frequency correction at each transient without additional acquisitions and compare it to water-suppressed MEGA-PRESS (WS MEGA) editing with intermittent prospective frequency correction.

METHODS

Macromolecule-suppressed gamma aminobutyric acid (GABA)-edited experiments were performed in a phantom and in the occipital lobe (OCC) (n = 12) and medial prefrontal cortex (mPFC) (n = 8) of the human brain. Water frequency consistency and average offset over acquisition time were compared. GABA multiplet patterns, signal intensities, and choline subtraction artifacts were evaluated. In vivo GABA concentrations were compared and related to frequency offset in the OCC.

RESULTS

MC MEGA was more stable with 21% and 32% smaller water frequency SDs in the OCC and mPFC, respectively. MC MEGA also had 39% and 40% smaller average frequency offsets in the OCC and mPFC, respectively. Phantom GABA multiplet patterns and signal intensities were similar. In vivo GABA concentrations were smaller in MC MEGA than in WS MEGA, with median (interquartile range) of 2.52 (0.27) and 2.29 (0.19) institutional units (i.u.), respectively in the OCC scans without prior DTI, and 0.99 (0.3) and 1.72 (0.5), respectively in the mPFC. OCC WS MEGA GABA concentrations, but not MC MEGA GABA concentrations were moderately correlated with frequency offset. mPFC WS MEGA spectra contained significantly more subtraction artifacts than MC MEGA spectra.

CONCLUSION

MC MEGA is feasible and allows for prospective frequency correction at every transient. MC MEGA GABA concentrations were not biased by frequency offsets and contained less subtraction artifacts compared to WS MEGA.

ProFit-1D-A 1D fitting software and open-source validation data sets

PURPOSE

Accurate and precise MRS fitting is crucial for metabolite concentration quantification of ¹ H-MRS spectra. LCModel, a spectral fitting software, has shown to have certain limitations to perform advanced spectral fitting by previous literature. Herein, we propose an open-source spectral fitting algorithm with adaptive spectral baseline determination and more complex cost functions.

THEORY

The MRS spectra are characterized by several parameters, which reflect the environment of the contributing metabolites, properties of the acquisition sequence, or additional disturbances. Fitting parameters should accurately describe these parameters. Baselines are also a major contributor to MRS spectra, in which smoothness of the spline baselines used for fitting can be adjusted based on the properties of the spectra. Three different cost functions used for the minimization problem were also investigated.

METHODS

The newly developed ProFit-1D fitting algorithm is systematically evaluated for simulations of several types of possible in vivo parameter variations. Although accuracy and precision are tested with simulated spectra, spectra measured in vivo at 9.4 T are used for testing precision using subsets of averages. ProFit-1D fitting results are also compared with LCModel.

RESULTS

Both ProFit-1D and LCModel fitted the spectra well with induced parameter and baseline variations. ProFit-1D proved to be more accurate than LCModel for simulated spectra. However, LCModel showed a somewhat increased precision for some spectral simulations and for in vivo data.

CONCLUSION

The open-source ProFit-1D fitting algorithm demonstrated high accuracy while maintaining precise metabolite concentration quantification. Finally, through the newly proposed cost functions, new ways to improve fitting were shown.

Relaxation-corrected macromolecular model enables determination of 1 H longitudinal T1 -relaxation times and concentrations of human brain metabolites at 9.4T

PURPOSE

Ultrahigh field MRS has improved characterization of the neurochemical profile. To compare results obtained at 9.4T to those from lower field strengths, it is of interest to quantify the concentrations of metabolites measured. Thus, measuring T₁ -relaxation times is necessary to correct for T₁ -weighting that occurs in acquisitions for single-voxel spectroscopy and spectroscopic imaging. A macromolecule (MM) simulation model was developed to fit MM contributions to the short TE inversion series used to measure T₁ -relaxation times.

METHODS

An inversion series with seven time points was acquired with metabolite-cycled STEAM to estimate T₁ -relaxation times of metabolites. A short TE was employed in this study to retain signals from metabolites with short T₂ -relaxation times and J-couplings. The underlying macromolecule spectrum was corrected by developing a sequence-specific, relaxation-corrected simulated MM model. Quantification of metabolite peaks was performed using internal water referencing and relaxation corrections.

RESULTS

T₁ -relaxation times for metabolites range from approximately 750 to approximately 2000 ms and approximately 1000 to approximately 2400 ms in gray matter (GM)- and white matter (WM)- rich voxels, respectively. Quantification of metabolites was compared between GM and WM voxels, as well as between results that used a simulated MM spectrum against those that used an experimentally acquired MM spectrum. Metabolite concentrations are reported in mmol/kg quantities.

CONCLUSION

T₁ -relaxation times are reported for nonsinglet resonances for the first time at 9.4T by use of a MM simulation model to account for contributions from the MM spectrum. In addition to T₁ -relaxation times, quantification results of metabolites from GM- and WM-rich voxels are reported.

Water and lipid suppression techniques for advanced 1 H MRS and MRSI of the human brain: Experts' consensus recommendations

The neurochemical information provided by proton magnetic resonance spectroscopy (MRS) or MR spectroscopic imaging (MRSI) can be severely compromised if strong signals originating from brain water and extracranial lipids are not properly suppressed. The authors of this paper present an overview of advanced water/lipid-suppression techniques and describe their advantages and disadvantages. Moreover, they provide recommendations for choosing the most appropriate techniques for proper use. Methods of water signal handling are primarily focused on the VAPOR technique and on MRS without water suppression (metabolite cycling). The section on lipid-suppression methods in MRSI is divided into three parts. First, lipid-suppression techniques that can be implemented on most clinical MR scanners (volume preselection, outer-volume suppression, selective lipid suppression) are described. Second, lipid-suppression techniques utilizing the combination of k-space filtering, high spatial resolutions and lipid regularization are presented. Finally, three promising new lipid-suppression techniques, which require special hardware (a multi-channel transmit system for dynamic B₁⁺ shimming, a dedicated second-order gradient system or an outer volume crusher coil) are introduced.

Met carriers of the BDNF Val66Met polymorphism show reduced Glx/NAA in the pregenual ACC in two independent cohorts

The Met allele of the Val66Met SNP of the BDNF gene (rs6265) is associated with impaired activity-dependent release of brain-derived neurotrophic factor (BDNF), resulting in reduced synaptic plasticity, impaired glutamatergic neurotransmission, and morphological changes. While previous work has demonstrated Val66Met effects on magnetic resonance spectroscopy (MRS) markers of either glutamatergic metabolism (Glx) or neuronal integrity (NAA), no study has investigated Val66Met effects on these related processes simultaneously. As these metabolites share a metabolic pathway, the Glx/NAA ratio may be a more sensitive marker of changes associated with the Val66Met SNP. This ratio is increased in psychiatric disorders linked to decreased functioning in the anterior cingulate cortex (ACC). In this study, we investigated the correlation of the Val66Met polymorphism of the BDNF gene with Glx/NAA in the pregenual anterior cingulate cortex (pgACC) using MRS at 3 Tesla (T) (n = 30, all males) and 7 T (n = 98, 40 females). In both cohorts, Met carriers had lower Glx/NAA compared to Val homozygotes. Follow-up analyses using absolute quantification revealed that the Met carriers do not show decreased pgACC glutamate or glutamine levels, but instead show increased NAA compared to the Val homozygotes. This finding may in part explain conflicting evidence for Val66Met as a risk factor for developing psychiatric illnesses.

Comparison of four 31P single-voxel MRS sequences in the human brain at 9.4 T

PURPOSE

In this study, different single-voxel localization sequences were implemented and systematically compared for the first time for phosphorous MRS (³¹ P-MRS) in the human brain at 9.4 T.

METHODS

Two multishot sequences, image-selected in vivo spectroscopy (ISIS) and a conventional slice-selective excitation combined with localization by adiabatic selective refocusing (semiLASER) variant of the spin-echo full intensity-acquired localized spectroscopy (SPECIAL-semiLASER), and two single-shot sequences, semiLASER and stimulated echo acquisition mode (STEAM), were implemented and optimized for ³¹ P-MRS in the human brain at 9.4 T. Pulses and coil setup were optimized, localization accuracy was tested in phantom experiments, and absolute SNR of the sequences was compared in vivo. The SNR per unit time (SNR/t) was derived and compared for all four sequences and verified experimentally for ISIS in two different voxel sizes (3 × 3 × 3 cm³ , 5 × 5 × 5 cm³ , 10-minute measurement time). Metabolite signals obtained with ISIS were quantified. The possible spectral quality in vivo acquired in clinically feasible time (3:30 minutes, 3 × 3 × 3 cm³ ) was explored for two different coil setups.

RESULTS

All evaluated sequences performed with good localization accuracy in phantom experiments and provided well-resolved spectra in vivo. However, ISIS has the lowest chemical shift displacement error, the best localization accuracy, the highest SNR/t for most metabolites, provides metabolite concentrations comparable to literature values, and is the only one of the sequences that allows for the detection of the whole ³¹ P spectrum, including β-adenosine triphosphate, with the used setup. The SNR/t of STEAM is comparable to the SNR/t of ISIS. The semiLASER and SPECIAL-semiLASER sequences provide good results for metabolites with long T₂ .

CONCLUSION

At 9.4 T, high-quality single-voxel localized ³¹ P-MRS can be performed in the human brain with different localization methods, each with inherent characteristics suitable for different research issues.

Repeatability of proton magnetic resonance spectroscopy of the brain at 7 T: effect of scan time on semi-localized by adiabatic selective refocusing and short-echo time stimulated echo acquisition mode scans and their comparison

Background

Proton magnetic resonance spectroscopy (MRS) provides a unique opportunity for in vivo measurements of the brain's metabolic profile. Two methods of mainstream data acquisition are compared at 7 T, which provides certain advantages as well as challenges. The two representative methods have seldom been compared in terms of measured metabolite concentrations and different scan times. The current study investigated proton MRS of the posterior cingulate cortex using a semi-localized by adiabatic selective refocusing (sLASER) sequence and a short echo time (TE) stimulated echo acquisition mode (sSTEAM) sequence, and it compared their reliability and repeatability at 7 T using a 32-channel head coil.

Methods

Sixteen healthy subjects were prospectively enrolled and scanned twice with an off-bed interval between scans. The scan parameters for sLASER were a TR/TE of 6.5 s/32 ms and 32 and 48 averages (sLASER×32 and sLASER×48, respectively). The scan parameters for sSTEAM were a TR/TE of 4 s/5 ms and 32, 48, and 64 averages (sSTEAM4×32, sSTEAM4×48, and sSTEAM4×64, respectively) in addition to that with a TR/TE of 8 s/5 ms and 32 averages (sSTEAM8×32). Data were analyzed using LCModel. Metabolites quantified with Cramér-Rao lower bounds (CRLBs) >50% were classified as not detected, and metabolites quantified with mean or median CRLBs ≤20% were included for further analysis. The SNR, CRLBs, coefficient of variation (CV), and metabolite concentrations were statistically compared using the Shapiro-Wilk test, one-way ANOVA, or the Friedman test.

Results

The sLASER spectra for N-acetylaspartate + N-acetylaspartylglutamate (tNAA) and glutamate (Glu) had a comparable or higher SNR than sSTEAM spectra. Ten metabolites had lower CRLBs than prefixed thresholds: aspartate (Asp), γ-aminobutyric acid (GABA), glutamine (Gln), Glu, glutathione (GSH), myo-inositol (Ins), taurine (Tau), the total amount of phosphocholine + glycerophosphocholine (tCho), creatine + phosphocreatine (tCr), and tNAA. Performance of the two sequences was satisfactory except for GABA, for which sLASER yielded higher CRLBs (≥18%) than sSTEAM. Some significant differences in CRLBs were noted, but they were ≤2% except for GABA and Gln. Signal averaging significantly lowered CRLBs for some metabolites but only by a small amount. Measurement repeatability as indicated by median CVs was ≤10% for Gln, Glu, Ins, tCho, tCr, and tNAA in all scans, and that for Asp, GABA, GSH, and Tau was ≥10% under some scanning conditions. The CV for GABA according to sLASER was significantly higher than that according to sSTEAM, whereas the CV for Ins was higher according to sSTEAM. An increase in signal averaging contribute little to lower CVs except for Ins.

Conclusions

Both sequences quantified brain metabolites with a high degree of precision and repeatability. They are comparable except for GABA, for which sSTEAM would be a better choice.

Localized Prediction of Glutamate from Whole-Brain Functional Connectivity of the Pregenual Anterior Cingulate Cortex

Local measures of neurotransmitters provide crucial insights into neurobiological changes underlying altered functional connectivity in psychiatric disorders. However, noninvasive neuroimaging techniques such as magnetic resonance spectroscopy (MRS) may cover anatomically and functionally distinct areas, such as p32 and p24 of the pregenual anterior cingulate cortex (pgACC). Here, we aimed to overcome this low spatial specificity of MRS by predicting local glutamate and GABA based on functional characteristics and neuroanatomy in a sample of 88 human participants (35 females), using complementary machine learning approaches. Functional connectivity profiles of pgACC area p32 predicted pgACC glutamate better than chance (R2 = 0.324) and explained more variance compared with area p24 using both elastic net and partial least-squares regression. In contrast, GABA could not be robustly predicted. To summarize, machine learning helps exploit the high resolution of fMRI to improve the interpretation of local neurometabolism. Our augmented multimodal imaging analysis can deliver novel insights into neurobiology by using complementary information.SIGNIFICANCE STATEMENT Magnetic resonance spectroscopy (MRS) measures local glutamate and GABA noninvasively. However, conventional MRS requires large voxels compared with fMRI, because of its inherently low signal-to-noise ratio. Consequently, a single MRS voxel may cover areas with distinct cytoarchitecture. In the largest multimodal 7 tesla machine learning study to date, we overcome this limitation by capitalizing on the spatial resolution of fMRI to predict local neurotransmitters in the PFC. Critically, we found that prefrontal glutamate could be robustly and exclusively predicted from the functional connectivity fingerprint of one of two anatomically and functionally defined areas that form the pregenual anterior cingulate cortex. Our approach provides greater spatial specificity on neurotransmitter levels, potentially improving the understanding of altered functional connectivity in mental disorders.

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

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

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.

A novel method to measure T1 -relaxation times of macromolecules and quantification of the macromolecular resonances

PURPOSE

Macromolecular peaks underlying metabolite spectra influence the quantification of metabolites. Therefore, it is important to understand the extent of contribution from macromolecules (MMs) in metabolite quantification. However, to model MMs more accurately in spectral fitting, differences in T₁ relaxation times among individual MM peaks must be considered. Characterization of T₁ -relaxation times for all individual MM peaks using a single inversion recovery technique is difficult due to eventual contributions from metabolites. On the contrary, a double inversion recovery (DIR) technique provided flexibility to acquire MM spectra spanning a range of longitudinal magnetizations with minimal metabolite influence. Thus, a novel method to determine T₁ -relaxation times of individual MM peaks is reported in this work.

METHODS

Extensive Bloch simulations were performed to determine inversion time combinations for a DIR technique that yielded adequate MM signal with varying longitudinal magnetizations while minimizing metabolite contributions. MM spectra were acquired using DIR-metabolite-cycled semi-LASER sequence. LCModel concentrations were fitted to the DIR signal equation to calculate T₁ -relaxation times.

RESULTS

T₁ -relaxation times of MMs range from 204 to 510 ms and 253 to 564 ms in gray- and white-matter rich voxels respectively at 9.4T. Additionally, concentrations of 13 MM peaks are reported.

CONCLUSION

A novel DIR method is reported in this work to calculate T₁ -relaxation times of MMs in the human brain. T₁ -relaxation times and relaxation time corrected concentrations of individual MMs are reported in gray- and white-matter rich voxels for the first time at 9.4T.

In vivo characterization of downfield peaks at 9.4 T: T2 relaxation times, quantification, pH estimation, and assignments

PURPOSE

Relaxation times are a valuable asset when determining spectral assignments. In this study, apparent T₂ relaxation times ( T 2 app ) of downfield peaks are reported in the human brain at 9.4 T and are used to guide spectral assignments of some downfield metabolite peaks.

METHODS

Echo time series of downfield metabolite spectra were acquired at 9.4 T using a metabolite-cycled semi-LASER sequence. Metabolite spectral fitting was performed using LCModel V6.3-1L while fitting a pH sweep to estimate the pH of the homocarnosine (hCs) imidazole ring. T 2 app were calculated by fitting the resulting relative amplitudes of the peaks to a mono-exponential decay across the TE series. Furthermore, estimated tissue concentrations of molecules were calculated using the relaxation times and internal water as a reference.

RESULTS

T 2 app of downfield metabolites are reported within a range from 16 to 32 ms except for homocarnosine with T 2 app of 50 ms. Correcting T 2 app for exchange rates ( T 2 c o r r ) resulted in relaxation times between 20 and 33 ms. The estimated pH values based on hCs imidazole range from 7.07 to 7.12 between subjects. Furthermore, analyzing the linewidths of the downfield peaks and their T 2 app contribution led to possible peak assignments.

CONCLUSION

T 2 app relaxation times were longer for the assigned metabolite peaks compared to the unassigned peaks. Tissue pH estimation in vivo with proton MRS and simultaneous quantification of amide protons at 8.30 ± 0.15 ppm is likely possible. Based on concentration, linewidth, and exchange rates measurements, tentative peak assignments are discussed for adenosine triphosphate (ATP), N-acetylaspartylglutamate (NAAG), and urea.

Metabolite cycled liver 1 H MRS on a 7 T parallel transmit system

INTRODUCTION

Single-voxel ¹ H MRS in body applications often suffers from respiratory and other motion induced phase and frequency shifts, which lead to incoherent averaging and hence to suboptimal results.

METHODS

Here we show the application of metabolite cycling (MC) for liver STEAM-localized ¹ H MRS on a 7 T parallel transmit system, using eight transmit-receive fractionated dipole antennas with 16 additional, integrated receive loops. MC-STEAM measurements were made in six healthy, lean subjects and compared with STEAM measurements using VAPOR water suppression. Measurements were performed during free breathing and during synchronized breathing, for which the subjects did breathe in between the MRS acquisitions. Both intra-session repeatability and inter-session reproducibility of liver lipid quantification with MC-STEAM and VAPOR-STEAM were determined.

RESULTS

The preserved water signal in MC-STEAM allowed for robust phase and frequency correction of individual acquisitions before averaging, which resulted in in vivo liver spectra that were of equal quality when measurements were made with free breathing or synchronized breathing. Intra-session repeatability and inter-session reproducibility of liver lipid quantification were better for MC-STEAM than for VAPOR-STEAM. This may also be explained by the more robust phase and frequency correction of the individual MC-STEAM acquisitions as compared with the VAPOR-STEAM acquisitions, for which the low-signal-to-noise ratio lipid signals had to be used for the corrections.

CONCLUSION

Non-water-suppressed MC-STEAM on a 7 T system with parallel transmit is a promising approach for ¹ H MRS applications in the body that are affected by motion, such as in the liver, and yields better repeatability and reproducibility compared with water-suppressed measurements.

Non-water-excitation MR spectroscopy techniques to explore exchanging protons in human brain at 3 T

PURPOSE

To develop localization sequences for in vivo MR spectroscopy (MRS) on clinical scanners of 3 T to record spectra that are not influenced by magnetization transfer from water.

METHODS

Image-selected in vivo spectroscopy (ISIS) localization and chemical-shift-selective excitation (termed I-CSE) was combined in two ways: first, full ISIS localization plus a frequency-selective spin-echo and second, two-dimensional (2D) ISIS plus a frequency-selective excitation and slice-selective refocusing. The techniques were evaluated at 3 T in phantoms and human subjects in comparison to standard techniques with water presaturation or metabolite-cycling. ISIS included gradient-modulated offset-independent adiabatic (GOIA)-type adiabatic inversion pulses; echo times were 8-10 ms.

RESULTS

The novel 2D and 3D I-CSE methods yield upfield spectra that are comparable to those from standard MRS, except for shorter echo times and a limited frequency range. On the downfield/high-frequency side, they yield much more signal for exchangeable protons when compared to MRS with water presaturation or metabolite-cycling and longer echo times.

CONCLUSION

Novel non-water-excitation MRS sequences offer substantial benefits for the detection of metabolite signals that are otherwise suppressed by saturation transfer from water. Avoiding water saturation and using very short echo times allows direct observation of faster exchanging moieties than was previously possible at 3 T and additionally makes the methods less susceptible to fast T₂ relaxation.

Improving MRI's slice selectivity in the presence of strong, metal-derived inhomogeneities

PURPOSE

To develop schemes that deliver faithful 2D slices near field heterogeneities of the kind arising from non-ferromagnetic metal implants, with reduced artifacts and shorter scan times.

METHODS

An excitation scheme relying on cross-term spatio-temporal encoding (xSPEN) was used as basis for developing the new inhomogeneity-insensitive, slice-selective pulse scheme. The resulting Fully refOCUSED cross-term SPatiotemporal ENcoding (FOCUSED-xSPEN) approach involved four adiabatic sweeps. The method was evaluated in silico, in vitro and in vivo using mice models, and compared against a number of existing and of novel alternatives based on both conventional and swept RF pulses, including an analogous method based on LASER's selectivity spatial selectivity.

RESULTS

Calculations and experiments confirmed that multi-sweep derivatives of xSPEN and LASER can deliver localized excitation profiles, centered at the intended positions and endowed with enhanced immunity to B₀ and B₁ distortions. This, however, is achieved at the expense of higher SAR than non-swept counterparts. Furthermore, single-shot FOCUSED-xSPEN and LASER profiles covered limited off-resonance ranges. This could be extended to bands covering arbitrary off-resonance values with uniform slice widths, by looping the experiments over a number of scans possessing suitable transmission and reception offsets.

CONCLUSIONS

A series of novel approaches were introduced to select slices near metals, delivering robustness against B_o and B₁⁺ field inhomogeneities.

General technical remarks on 1HMRS translational research in 7T

Purpose

The aim of the work was to share the practical experience of preclinical and clinical proton magnetic resonance spectroscopy (¹HMRS) studies conducted using a 7-Tesla magnetic field strength scanner, taking into account the specificity of both settings in the context of translational research.

Material and methods

¹HMRS volunteer studies conducted using a Discovery 950 GE 7T scanner, were carried out with PRESS sequence, and a VOI measuring 2.0 × 2.0 × 2.0 cm³ placed in the white matter at the parietal occipital lobe. Rodent spectra obtained using a 7T Bruker were measured with PRESS, with a VOI 2.0 × 2.0 × 5.5 mm³ placed over the hippocampus.

Results

¹HMRS data from humans and rats show that the brain spectra obtained in the same field are characterised by a similar neurochemical structure and spectral resolution. Spectra obtained from rats demonstrate the following metabolites: NAA, Glu, Gln, Ins, Cho, Cr, PCr, Tau, GABA, Lac, NAAG, and Asp. In turn, spectra from humans allowed estimation of the following metabolites: Ala, NAA, Glu, Gln, Ins, Cho, Cr, PCr, Tau, GABA, Lac, NAAG, and Asp. Signals from Gln, Glu with chemical shift around 2.4 ppm, from Cr, PCr, and GABA at 3 ppm, and signals from Cho and Tau at approximately 3.2 ppm, can be properly separated and estimated both in humans and in rats.

Conclusions

These results are promising in terms of broadening the knowledge of many neurological diseases by inducing them on animal models and then transferring this knowledge to clinical practice. In spite of this, important distinctions in the technical aspects and methodological differences of high-field ¹HMRS in both preclinical and clinical conditions should be taken into account.

Metabolite-cycled density-weighted concentric rings k-space trajectory (DW-CRT) enables high-resolution 1 H magnetic resonance spectroscopic imaging at 3-Tesla

Magnetic resonance spectroscopic imaging (MRSI) is a promising technique in both experimental and clinical settings. However, to date, MRSI has been hampered by prohibitively long acquisition times and artifacts caused by subject motion and hardware-related frequency drift. In the present study, we demonstrate that density weighted concentric ring trajectory (DW-CRT) k-space sampling in combination with semi-LASER excitation and metabolite-cycling enables high-resolution MRSI data to be rapidly acquired at 3 Tesla. Single-slice full-intensity MRSI data (short echo time (TE) semi-LASER TE = 32 ms) were acquired from 6 healthy volunteers with an in-plane resolution of 5 × 5 mm in 13 min 30 sec using this approach. Using LCModel analysis, we found that the acquired spectra allowed for the mapping of total N-acetylaspartate (median Cramer-Rao Lower Bound [CRLB] = 3%), glutamate+glutamine (8%), and glutathione (13%). In addition, we demonstrate potential clinical utility of this technique by optimizing the TE to detect 2-hydroxyglutarate (long TE semi-LASER, TE = 110 ms), to produce relevant high-resolution metabolite maps of grade III IDH-mutant oligodendroglioma in a single patient. This study demonstrates the potential utility of MRSI in the clinical setting at 3 Tesla.

Combination of surface and 'vertical' loop elements improves receive performance of a human head transceiver array at 9.4 T

Ultra-high-field (UHF, ≥7 T) human magnetic resonance imaging (MRI) provides undisputed advantages over low-field MRI (≤3 T), but its development remains challenging because of numerous technical issues, including the low efficiency of transmit (Tx) radiofrequency (RF) coils caused by the increase in tissue power deposition with frequency. Tight-fit human head transceiver (TxRx) arrays improve Tx efficiency in comparison with Tx-only arrays, which are larger in order to fit multi-channel receive (Rx)-only arrays inside. A drawback of the TxRx design is that the number of elements in an array is limited by the number of available high-power RF Tx channels (commonly 8 or 16), which is not sufficient for optimal Rx performance. In this work, as a proof of concept, we developed a method for increasing the number of Rx elements in a human head TxRx surface loop array without the need to move the loops away from a sample, which compromises the array Tx performance. We designed and constructed a prototype 16-channel tight-fit array, which consists of eight TxRx surface loops placed on a cylindrical holder circumscribing a head, and eight Rx-only vertical loops positioned along the central axis (parallel to the magnetic field B₀ ) of each TxRx loop, perpendicular to its surface. We demonstrated both experimentally and numerically that the addition of the vertical loops has no measurable effect on the Tx efficiency of the array. An increase in the maximum local specific absorption rate (SAR), evaluated using two human head voxel models (Duke and Ella), measured 3.4% or less. At the same time, the 16-element array provided 30% improvement of central signal-to-noise ratio (SNR) in vivo relative to a surface loop eight-element array. The novel array design also demonstrated an improvement in the parallel Rx performance in the transversal plane. Thus, using this method, both the Rx and Tx performance of the human head array can be optimized simultaneously.