Dynamic multi-coil technique (DYNAMITE) shimming for echo-planar imaging of the human brain at 7 Tesla

High-quality lipid suppression and B0 shimming for human brain 1H MRSI

Magnetic Resonance Spectroscopic Imaging (MRSI) is a powerful technique that can map the metabolic profile in the brain non-invasively. Extracranial lipid contamination and insufficient B0 homogeneity however hampers robustness, and as a result has hindered widespread use of MRSI in clinical and research settings. Over the last six years we have developed highly effective extracranial lipid suppression methods with a second order gradient insert (ECLIPSE) utilizing inner volume selection (IVS) and outer volume suppression (OVS) methods. While ECLIPSE provides > 100-fold in lipid suppression with modest radio frequency (RF) power requirements and immunity to B1+ field variations, axial coverage is reduced for non-elliptical head shapes. In this work we detail the design, construction, and utility of MC-ECLIPSE, a pulsed second order gradient coil with Z2 and X2Y2 fields, combined with a 54-channel multi-coil (MC) array. The MC-ECLIPSE platform allows arbitrary region of interest (ROI) shaped OVS for full-axial slice coverage, in addition to MC-based B0 field shimming, for robust human brain proton MRSI. In vivo experiments demonstrate that MC-ECLIPSE allows axial brain coverage of 92-95 % is achieved following arbitrary ROI shaped OVS for various head shapes. The standard deviation (SD) of the residual B0 field following SH2 and MC shimming were 25 ± 9 Hz and 18 ± 8 Hz over a 5 cm slab, and 18 ± 5 Hz and 14 ± 6 Hz over a 1.5 cm slab, respectively. These results demonstrate that B0 magnetic field shimming with the MC array supersedes second order harmonic capabilities available on standard MRI systems for both restricted and large ROIs. Furthermore, MC based B0 shimming provides comparable shimming performance to an unrestricted SH5 shim set for both restricted, and 5-cm slab shim challenges. Phantom experiments demonstrate the high level of localization performance achievable with MC-ECLIPSE, with ROI edge chemical shift displacements ranging from 1-3 mm with a median value of 2 mm, and transition width metrics ranging from 1-2.5 mm throughout the ROI edge. Furthermore, MC based B0 shimming is comparable to performance following a full set of unrestricted spherical harmonic fields up to order 5. Short echo time MRSI and GABA-edited MRSI acquisitions in the human brain following MC-shimming and arbitrary ROI shaping demonstrate full-axial slice coverage and extracranial lipid artifact free spectra. MC-ECLIPSE allows full-axial coverage and robust MRSI acquisitions, while allowing interrogation of cortical tissue proximal to the skull, which has significant value in a wide range of neurological and psychiatric conditions.

Advancements in MR hardware systems and magnetic field control: B0 shimming, RF coils, and gradient techniques for enhancing magnetic resonance imaging and spectroscopy

High magnetic field homogeneity is critical for magnetic resonance imaging (MRI), functional MRI, and magnetic resonance spectroscopy (MRS) applications. B0 inhomogeneity during MR scans is a long-standing problem resulting from magnet imperfections and site conditions, with the main issue being the inhomogeneity across the human body caused by differences in magnetic susceptibilities between tissues, resulting in signal loss, image distortion, and poor spectral resolution. Through a combination of passive and active shim techniques, as well as technological advances employing multi-coil techniques, optimal coil design, motion tracking, and real-time modifications, improved field homogeneity and image quality have been achieved in MRI/MRS. The integration of RF and shim coils brings a high shim efficiency due to the proximity of participants. This technique will potentially be applied to high-density RF coils with a high-density shim array for improved B0 homogeneity. Simultaneous shimming and image encoding can be achieved using multi-coil array, which also enables the development of novel encoding methods using advanced magnetic field control. Field monitoring enables the capture and real-time compensation for dynamic field perturbance beyond the static background inhomogeneity. These advancements have the potential to better use the scanner performance to enhance diagnostic capabilities and broaden applications of MRI/MRS in a variety of clinical and research settings. The purpose of this paper is to provide an overview of the latest advances in B0 magnetic field shimming and magnetic field control techniques as well as MR hardware, and to emphasize their significance and potential impact on improving the data quality of MRI/MRS.

Artificial intelligence for neuro MRI acquisition: a review

OBJECT

To review recent advances of artificial intelligence (AI) in enhancing the efficiency and throughput of the MRI acquisition workflow in neuroimaging, including planning, sequence design, and correction of acquisition artifacts.

MATERIALS AND METHODS

A comprehensive analysis was conducted on recent AI-based methods in neuro MRI acquisition. The study focused on key technological advances, their impact on clinical practice, and potential risks associated with these methods.

RESULTS

The findings indicate that AI-based algorithms have a substantial positive impact on the MRI acquisition process, improving both efficiency and throughput. Specific algorithms were identified as particularly effective in optimizing acquisition steps, with reported improvements in workflow efficiency.

DISCUSSION

The review highlights the transformative potential of AI in neuro MRI acquisition, emphasizing the technological advances and clinical benefits. However, it also discusses potential risks and challenges, suggesting areas for future research to mitigate these concerns and further enhance AI integration in MRI acquisition.

High resolution simulation and measurement demonstrate oscillatory spatiotemporal B0 fluctuations across the human cardiac cycle

PURPOSE

Functional cardiac MRI scans employing balanced steady-state free precession sequences suffer from dark band artifacts in the myocardium due to B0 inhomogeneity. We recently introduced a novel method for the theoretical derivation of B0 distributions in the human heart. This study aims to simulate the B0 distributions in the heart across the cardiac cycle using structural MR images and validate the simulations via in vivo measured cardiac phase-specific B0 maps on the same subjects at 3T.

METHODS

Cardiac phase-specific B0 field maps were acquired from eight healthy subjects at 3T. B0 conditions were simulated based on tissue masks created from the cardiac-phase specific structural images from the in vivo B0 map scan and anatomical images from a thoracic MRI scan, adopting our recently published approach. The simulations and in vivo measurements were compared by calculating the spatial correlation of their B0 distributions and temporal correlation of the derived spherical harmonic coefficients throughout the cardiac cycle.

RESULTS

The spatial comparison of B0 maps between the simulation and in vivo measurement indicates an overall average correlation coefficient of 0.91 across the cardiac cycle in all subjects. Both groups show consistent high-level B0 patterns. Temporal variations of B0 conditions exhibit sinusoidal characteristics and are strongly correlated between simulation and in vivo.

CONCLUSION

Theoretical simulations employing regional anatomical features were validated by direct in vivo B0 mapping in the same subjects. The spatial B0 condition throughout the cardiac cycle exhibits oscillatory characteristics due to structural distortions of cardiac motion.

3D FRONSAC with PSF reconstruction

PURPOSE

This study extends the Fast ROtary Nonlinear Spatial ACquisition (FRONSAC) method to include 3D acquisitions and reconstructions. It uses a transform domain reconstruction which is needed to make 3D reconstructions practical and provides new insights into how parallel imaging performance is enhanced by FRONSAC encoding.

METHODS

This work developed the first examples of FRONSAC incorporated into a 3D acquisition. 3D FRONSAC was tested on human subjects with both simple gradient echo and MPRAGE Cartesian acquisitions. The quality of the 3D FRONSAC images was evaluated using structural similarity index measure (SSIM), and normalized root mean square error (NRMSE).

RESULTS

FRONSAC encoding did not significantly modify the contrast obtained in either sequence, but it substantially improves the image quality of undersampled reconstruction. FRONSAC images have reduced undersampling ghosts and consistently improved SSIM and NRMSE.

CONCLUSIONS

Acquisition and reconstruction of 3D FRONSAC images are feasible, and the additional FRONSAC encoding improves image quality in highly undersampled images.

In vivo cortical glutathione response to oral fumarate therapy in relapsing-remitting multiple sclerosis: A single-arm open-label phase IV trial using 7-Tesla 1H MRS

BACKGROUND

This is an open-label, single-arm, single-center pilot study using 7-Tesla in vivo proton magnetic resonance spectroscopy (1H MRS) to measure brain cortical glutathione concentration at baseline before and during the use of oral fumarates as a disease-modifying therapy for multiple sclerosis. The primary endpoint of this research was the change in prefrontal cortex glutathione concentration relative to a therapy-naïve baseline after one year of oral fumarate therapy.

METHODS

Brain glutathione concentrations were examined by 1H MRS in single prefrontal and occipital cortex cubic voxels (2.5 × 2.5 × 2.5 cm3) before and during initiation of oral fumarate therapy (120 mg b.i.d. for 7 days and 240 mg b.i.d. thereafter). Additional measurements of related metabolites glutamate, glutamine, myoinositol, total N-acetyl aspartate, and total choline were also acquired in voxels centered on the same regions. Seven relapsing-remitting multiple sclerosis patients (4 f / 3 m, age range 28-50 years, mean age 40 years) naïve to fumarate therapy were scanned at pre-therapy baseline and after 1, 3, 6 and 12 months of therapy. A group of 8 healthy volunteers (4 f / 4 m, age range 33-48 years, mean age 41 years) was also scanned at baseline and Month 6 to characterize 1H-MRS measurement reproducibility over a comparable time frame.

RESULTS

In the multiple sclerosis cohort, general linear models demonstrated a significant positive linear relationship between prefrontal glutathione and time either linearly across all time points (+0.05 ± 0.02 mM/month, t(27) = 2.6, p = 0.02) or specifically for factor variable Month 12 (+0.6 ± 0.3 mM/12 months, t(24) = 2.2, p = 0.04) relative to baseline. No such effects of time on glutathione concentration were demonstrated in the occipital cortex or in the healthy volunteer group. Changes in occipital total choline were further observed in the multiple sclerosis cohort as well as prefrontal total choline and occipital glutamine and myoinositol in the control cohort throughout the study duration.

CONCLUSIONS

While the open-label single-arm pilot study design and abbreviated control series cannot support firm conclusions about the influence of oral fumarate therapy independent of test-retest factors or normal biological variation in a state of either health or disease, these results do justify further investigation at a larger scale into the potential relationship between prefrontal cortex glutathione increases and oral fumarate therapy in relapsing-remitting multiple sclerosis.

Naturalistic Hyperscanning with Wearable Magnetoencephalography

The evolution of human cognitive function is reliant on complex social interactions which form the behavioural foundation of who we are. These social capacities are subject to dramatic change in disease and injury; yet their supporting neural substrates remain poorly understood. Hyperscanning employs functional neuroimaging to simultaneously assess brain activity in two individuals and offers the best means to understand the neural basis of social interaction. However, present technologies are limited, either by poor performance (low spatial/temporal precision) or an unnatural scanning environment (claustrophobic scanners, with interactions via video). Here, we describe hyperscanning using wearable magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs). We demonstrate our approach by simultaneously measuring brain activity in two subjects undertaking two separate tasks-an interactive touching task and a ball game. Despite large and unpredictable subject motion, sensorimotor brain activity was delineated clearly, and the correlation of the envelope of neuronal oscillations between the two subjects was demonstrated. Our results show that unlike existing modalities, OPM-MEG combines high-fidelity data acquisition and a naturalistic setting and thus presents significant potential to investigate neural correlates of social interaction.

Hybrid active and passive local shimming (HAPLS) for two-region MRI

PURPOSE

An MRI scanner is equipped with global shim systems for shimming one region of interest (ROI) only. However, it often fails to reach state-of-the-art when shimming two isolated regions of interest simultaneously, even though the two-area shimming can be essential in scan scenarios, such as bilateral breasts or dyadic brains. To address these challenges, a hybrid active and passive local shimming technique is proposed to simultaneously shim two isolated region-of-interest areas within the whole FOV.

METHODS

A local passive shimming system is constructed by optimized bilateral ferromagnetic chip arrays to compensate for the magnet's significant high-order B0 inhomogeneities at the boundary of the manufacturer's specified homogeneous volume, thus locally improving the available FOV. The local active shimming consists of 40-channel DC loops powered by 64-channel current amplifiers. With the optimized current distribution, active shimming can correct the residual low-order B0 inhomogeneities and subject-specific field inhomogeneities. In addition, active shimming is used to homogenize the center frequencies of the two regions.

RESULTS

With the implementation of the hybrid active and passive local shimming, the 95% peak-to-peak was reduced from 1.92 to 1.12 ppm by 41.7%, and RMS decreased from 0.473 to 0.255 ppm by 46.1% in a two-phantom experiment. The volume ratio containing MR voxels within a 0.5-ppm frequency span increased from 64.3% to 81.3% by 26.3%.

CONCLUSION

The proposed hybrid active and passive local shimming technique uses both passive and active local shimming, and it can efficiently shim two areas simultaneously, which is an unmet need for a commercial MRI scanner.

Multiple sclerosis diagnosis and phenotype identification by multivariate classification of in vivo frontal cortex metabolite profiles

Multiple sclerosis (MS) is a heterogeneous autoimmune disease for which diagnosis continues to rely on subjective clinical judgment over a battery of tests. Proton magnetic resonance spectroscopy (¹H MRS) enables the noninvasive in vivo detection of multiple small-molecule metabolites and is therefore in principle a promising means of gathering information sufficient for multiple sclerosis diagnosis and subtype classification. Here we show that supervised classification using ¹H-MRS-visible normal-appearing frontal cortex small-molecule metabolites alone can indeed differentiate individuals with progressive MS from control (held-out validation sensitivity 79% and specificity 68%), as well as between relapsing and progressive MS phenotypes (held-out validation sensitivity 84% and specificity 74%). Post hoc assessment demonstrated the disproportionate contributions of glutamate and glutamine to identifying MS status and phenotype, respectively. Our finding establishes ¹H MRS as a viable means of characterizing progressive multiple sclerosis disease status and paves the way for continued refinement of this method as an auxiliary or mainstay of multiple sclerosis diagnostics.

A lightweight magnetically shielded room with active shielding

Magnetically shielded rooms (MSRs) use multiple layers of materials such as MuMetal to screen external magnetic fields that would otherwise interfere with high precision magnetic field measurements such as magnetoencephalography (MEG). Optically pumped magnetometers (OPMs) have enabled the development of wearable MEG systems which have the potential to provide a motion tolerant functional brain imaging system with high spatiotemporal resolution. Despite significant promise, OPMs impose stringent magnetic shielding requirements, operating around a zero magnetic field resonance within a dynamic range of ± 5 nT. MSRs developed for OPM-MEG must therefore effectively shield external sources and provide a low remnant magnetic field inside the enclosure. Existing MSRs optimised for OPM-MEG are expensive, heavy, and difficult to site. Electromagnetic coils are used to further cancel the remnant field inside the MSR enabling participant movements during OPM-MEG, but present coil systems are challenging to engineer and occupy space in the MSR limiting participant movements and negatively impacting patient experience. Here we present a lightweight MSR design (30% reduction in weight and 40-60% reduction in external dimensions compared to a standard OPM-optimised MSR) which takes significant steps towards addressing these barriers. We also designed a 'window coil' active shielding system, featuring a series of simple rectangular coils placed directly onto the walls of the MSR. By mapping the remnant magnetic field inside the MSR, and the magnetic field produced by the coils, we can identify optimal coil currents and cancel the remnant magnetic field over the central cubic metre to just |B|= 670 ± 160 pT. These advances reduce the cost, installation time and siting restrictions of MSRs which will be essential for the widespread deployment of OPM-MEG.

High-resolution simulation of B0 field conditions in the human heart from segmented computed tomography images

B₀ inhomogeneity leads to imaging artifacts in cardiac magnetic resonance imaging (MRI), in particular dark band artifacts with steady-state free precession pulse sequences. The limited spatial resolution of MR-derived in vivo B₀ maps and the lack of population data prevent systematic analysis of the problem at hand and the development of optimized B₀ shim strategies. We used readily available clinical computed tomography (CT) images to simulate the B₀ conditions in the human heart at high spatial resolution. Calculated B₀ fields showed consistency with MRI-based B₀ measurements. The B₀ maps for both the simulations and in vivo measurements showed local field inhomogeneities in the vicinity of lung tips with dominant Z3 spherical harmonic terms in the field distribution. The presented simulation approach allows for the derivation of B₀ field conditions at high spatial resolution from CT images and enables the development of subject- and population-specific B₀ shim strategies for the human heart.

Fourier-based decomposition for simultaneous 2-voxel MRS acquisition with 2SPECIAL

PURPOSE

To simultaneously acquire spectroscopic signals from two MRS voxels using a multi-banded 2 spin-echo, full-intensity acquired localized (2SPECIAL) sequence, and to decompose the signal to their respective regions by a novel voxel-GRAPPA (vGRAPPA) decomposition approach for in vivo brain applications at 7 T.

METHODS

A wideband, uniform rate, smooth truncation (WURST) multi-banded pulse was incorporated into SPECIAL to implement 2SPECIAL for simultaneous multi-voxel spectroscopy (sMVS). To decompose the acquired data, the voxel-GRAPPA decomposition algorithm is introduced, and its performance is compared to the SENSE-based decomposition. Furthermore, the limitations of two-voxel excitation concerning the multi-banded adiabatic inversion pulse, as well as of the combined B₀ shim and B₁ ⁺ adjustments, are evaluated.

RESULTS

It was successfully shown that the 2SPECIAL sequence enables sMVS without a significant loss in SNR while reducing the total scan time by 21.6% compared to two consecutive acquisitions. The proposed voxel-GRAPPA algorithm properly reassigns the signal components to their respective origin region and shows no significant differences to the well-established SENSE-based algorithm in terms of leakage (both <10%) or Cramér-Rao lower bounds (CRLB) for in vivo applications, while not requiring the acquisition of additional sensitivity maps and thus decreasing motion sensitivity.

CONCLUSION

The use of 2SPECIAL in combination with the novel voxel-GRAPPA decomposition technique allows a substantial reduction of measurement time compared to the consecutive acquisition of two single voxels without a significant decrease in spectral quality or metabolite quantification accuracy and thus provides a new option for multiple-voxel applications.

Shim Coils Tailored for Correcting B0 Inhomogeneity in the Human Brain (SCOTCH): Design Methodology and 48-Channel Prototype Assessment in 7-Tesla MRI

Increased static field inhomogeneities are a burden for human brain MRI at Ultra-High-Field. In particular they cause enhanced Echo-Planar image distortions and signal losses due to magnetic susceptibility gradients at air-tissue interfaces in the subject's head. In the past decade, Multi-Coil Arrays (MCA) have been proposed to shim the field in the brain better than the 2nd or 3rd order Spherical Harmonic (SH) coils usually offered by MRI manufacturers. Here we present a novel MCA, named SCOTCH, optimized for whole brain shimming. Based on a cylindrical structure, it features several layers of small coils whose shape, size and location are found from a principal component analysis of ideal stream functions computed from an internal 100-brain fieldmap database. From an Open-Access external database of 126 brains, our SCOTCH implementation is shown to be equivalent to a partial 7th-order SH system with unlimited power, outperforming all known existing MCA prototypes. This result is further confirmed by a low-cost  30-cm diameter SCOTCH prototype built with 48 coils on 3 layers, and tested on 7 volunteers at 7T with a parallel-transmit RF coil made to be inserted in SCOTCH. Echo-Planar images of the subject brains before and after SCOTCH shimming show large signal recoveries, especially in the prefrontal cortex.

Physical limits to human brain B0 shimming with spherical harmonics, engineering implications thereof

OBJECTIVE

As the MRI main magnetic field rises for improved signal-to-noise ratio, susceptibility-induced B0-inhomogeneity increases proportionally, aggravating related artifacts. Considering only susceptibility disparities between air and biological tissue, we explore the topological conditions for which perfect shimming could be performed in a Region of Interest (ROI) such as the human brain or part thereof.

MATERIALS AND METHODS

After theoretical considerations for perfect shimming, spherical harmonic (SH) shimming simulations of very high degree are performed, based on a 100-subject database of 1.7-mm-resolved brain fieldmaps acquired at 3T . In addition to the whole brain, shimmed ROIs include slabs targeting the prefrontal cortex, both or single temporal lobes, or spheres in the frontal brain above the nasal sinus.

RESULTS AND DISCUSSION

We show "perfect" SH shimming is possible only if the ROI can be contained in a sphere that does not enclose sources of magnetic field inhomogeneity, which are gathered at the air-tissue interface. We establish a [Formula: see text]Hz inhomogeneity hard shim limit at 7T for whole brain SH shimming, that can only be attained at shimming degree higher than 90. On the other hand, under limited power and SH degree resources, 3D region-specific shimming is shown to greatly improve homogeneity in critical zones such as the prefrontal cortex and around ear canals.

A 31-channel integrated "AC/DC" B0 shim and radiofrequency receive array coil for improved 7T MRI

PURPOSE

To test an integrated "AC/DC" array approach at 7T, where B0 inhomogeneity poses an obstacle for functional imaging, diffusion-weighted MRI, MR spectroscopy, and other applications.

METHODS

A close-fitting 7T 31-channel (31-ch) brain array was constructed and tested using combined Rx and ΔB0 shim channels driven by a set of rapidly switchable current amplifiers. The coil was compared to a shape-matched 31-ch reference receive-only array for RF safety, signal-to-noise ratio (SNR), and inter-element noise correlation. We characterize the coil array's ability to provide global and dynamic (slice-optimized) shimming using ΔB0 field maps and echo planar imaging (EPI) acquisitions.

RESULTS

The SNR and average noise correlation were similar to the 31-ch reference array. Global and slice-optimized shimming provide 11% and 40% improvements respectively compared to baseline second-order spherical harmonic shimming. Birdcage transmit coil efficiency was similar for the reference and AC/DC array setups.

CONCLUSION

Adding ΔB0 shim capability to a 31-ch 7T receive array can significantly boost 7T brain B0 homogeneity without sacrificing the array's rdiofrequency performance, potentially improving ultra-high field neuroimaging applications that are vulnerable to off-resonance effects.

Emerging methods and applications of ultra-high field MR spectroscopic imaging in the human brain

Magnetic Resonance Spectroscopic Imaging (MRSI) of the brain enables insights into the metabolic changes and fluxes in diseases such as tumors, multiple sclerosis, epilepsy, or hepatic encephalopathy, as well as insights into general brain functionality. However, the routine application of MRSI is mostly hampered by very low signal-to-noise ratios (SNR) due to the low concentrations of metabolites, about 10000 times lower than water. Furthermore, MRSI spectra have a dense information content with many overlapping metabolite resonances, especially for proton MRSI. MRI scanners at ultra-high field strengths, like 7 T or above, offer the opportunity to increase SNR, as well as the separation between resonances, thus promising to solve both challenges. Yet, MRSI at ultra-high field strengths is challenged by decreased B0- and B1-homogeneity, shorter T2 relaxation times, stronger chemical shift displacement errors, and aggravated lipid contamination. Therefore, to capitalize on the advantages of ultra-high field strengths, these challenges must be overcome. This review focuses on the challenges MRSI of the human brain faces at ultra-high field strength, as well as the possible applications to this date.

In vivo evidence of differential frontal cortex metabolic abnormalities in progressive and relapsing-remitting multiple sclerosis

The pathophysiology of progressive multiple sclerosis remains elusive, significantly limiting available disease-modifying therapies. Proton MRS (¹ H-MRS) enables in vivo measurement of small molecules implicated in multiple sclerosis, but its application to key metabolites glutamate, γ-aminobutyric acid (GABA), and glutathione has been sparse. We employed, at 7 T, a previously validated ¹ H-MRS protocol to measure glutamate, GABA, and glutathione, as well as glutamine, N-acetyl aspartate, choline, and myoinositol, in the frontal cortex of individuals with relapsing-remitting (N = 26) or progressive (N = 21) multiple sclerosis or healthy control adults (N = 25) in a cross-sectional analysis. Only individuals with progressive multiple sclerosis demonstrated reduced glutamate (F_2,65  = 3.424, p = 0.04; 12.40 ± 0.62 mM versus control 13.17 ± 0.95 mM, p = 0.03) but not glutamine (F_2,65  = 0.352, p = 0.7; 4.71 ± 0.35 mM versus control 4.84 ± 0.42 mM), reduced GABA (F_2,65  = 3.89, p = 0.03; 1.29 ± 0.23 mM versus control 1.47 ± 0.25 mM, p = 0.05), and possibly reduced glutathione (F_2,65  = 0.352, p = 0.056; 2.23 ± 0.46 mM versus control 2.51 ± 0.48 mM, p < 0.1). As a group, multiple sclerosis patients demonstrated significant negative correlations between disease duration and glutamate or GABA (ρ = -0.4, p = 0.02) but not glutamine or glutathione. Alone, only relapsing-remitting multiple sclerosis patients exhibited a significant negative correlation between disease duration and GABA (ρ = -0.5, p = 0.03). Taken together, these results indicate that frontal cortex metabolism is differentially disturbed in progressive and relapsing-remitting multiple sclerosis.

Application of an integrated radio-frequency/shim coil technology for signal recovery in fMRI

PURPOSE

Gradient-echo echo-planar imaging (EPI), which is typically used for blood oxygenation level-dependent (BOLD) functional MRI (fMRI), suffers from distortions and signal loss caused by localized B₀ inhomogeneities. Such artifacts cannot be effectively corrected for with the low-order spherical harmonic (SH) shim coils available on most scanners. The integrated parallel reception, excitation, and shimming (iPRES) coil technology allows radiofrequency (RF) and direct currents to flow on each coil element, enabling imaging and localized B₀ shimming with one coil array. iPRES was previously used to correct for distortions in spin-echo EPI and is further developed here to also recover signal loss in gradient-echo EPI.

METHODS

The cost function in the shim optimization, which typically uses a single term representing the B₀ inhomogeneity, was modified to include a second term representing the signal loss, with an adjustable weight to optimize the trade-off between distortion correction and signal recovery. Simulations and experiments were performed to investigate the shimming performance.

RESULTS

Slice-optimized shimming with iPRES and the proposed cost function substantially reduced the signal loss in the inferior frontal and temporal brain regions compared to shimming with iPRES and the original cost function or 2^nd -order SH shimming with either cost function. In breath-holding fMRI experiments, the ΔB₀ and signal loss root-mean-square errors decreased by -34.3% and -56.2%, whereas the EPI signal intensity and number of activated voxels increased by 60.3% and 174.0% in the inferior frontal brain region.

CONCLUSION

iPRES can recover signal loss in gradient-echo EPI, which is expected to improve BOLD fMRI studies in brain regions suffering from signal loss.

High-resolution fMRI at 7 Tesla: challenges, promises and recent developments for individual-focused fMRI studies

Limited detection power has been a bottleneck for subject-specific functional MRI (fMRI) studies, however the higher signal-to-noise ratio afforded by ultra-high magnetic fields (≥ 7 Tesla) provides levels of sensitivity and resolution needed to study individual subjects. What may be surprising is that higher imaging resolution may provide both higher specificity and sensitivity due to reductions in partial volume effects and reduced physiological noise. However, challenges remain to ensure high data quality and to reduce variability in ultra-high field fMRI. We discuss session-specific biases including those caused by factors related to instrumentation, anatomy, and physiology-which can translate into variability across sessions-and how to minimize these to help ultra-high field fMRI reach its full potential for individual-focused studies.

MRI exploiting frequency-modulated pulses and their nonlinear phase

Frequency-modulated (FM) pulses can provide several advantages over conventional amplitude-modulated pulses in the field of MRI; however, the manner in which spins are manipulated imprints a quadratic phase on the resulting magnetization. Historically this was considered a hindrance and slowed the widespread adoption of FM pulses. This article seeks to provide a historical perspective of the different techniques that researchers have used to exploit the benefits of FM pulses and to compensate for the nonlinear phase created by this class of pulses in MRI. Expanding on existing techniques, a new method of phase compensation is presented that utilizes nonlinear gradients to mitigate the undesirable phase imparted by this class of pulses.

B0 shimming of the human heart at 7T

PURPOSE

Inhomogeneities of the static magnetic B₀ field are a major limiting factor in cardiac MRI at ultrahigh field (≥ 7T), as they result in signal loss and image distortions. Different magnetic susceptibilities of the myocardium and surrounding tissue in combination with cardiac motion lead to strong spatio-temporal B₀ -field inhomogeneities, and their homogenization (B₀ shimming) is a prerequisite. Limitations of state-of-the-art shimming are described, regional B₀ variations are measured, and a methodology for spherical harmonics shimming of the B₀ field within the human myocardium is proposed.

METHODS

The spatial B₀ -field distribution in the heart was analyzed as well as temporal B₀ -field variations in the myocardium over the cardiac cycle. Different shim region-of-interest selections were compared, and hardware limitations of spherical harmonics B₀ shimming were evaluated by calibration-based B₀ -field modeling. The role of third-order spherical harmonics terms was analyzed as well as potential benefits from cardiac phase-specific shimming.

RESULTS

The strongest B₀ -field inhomogeneities were observed in localized spots within the left-ventricular and right-ventricular myocardium and varied between systolic and diastolic cardiac phases. An anatomy-driven shim region-of-interest selection allowed for improved B₀ -field homogeneity compared with a standard shim region-of-interest cuboid. Third-order spherical harmonics terms were demonstrated to be beneficial for shimming of these myocardial B₀ -field inhomogeneities. Initial results from the in vivo implementation of a potential shim strategy were obtained. Simulated cardiac phase-specific shimming was performed, and a shim term-by-term analysis revealed periodic variations of required currents.

CONCLUSION

Challenges in state-of-the-art B₀ shimming of the human heart at 7 T were described. Cardiac phase-specific shimming strategies were found to be superior to vendor-supplied shimming.

Optimal echo times for multi-gradient echo-based B0 field-mapping

B₀ field maps are used ubiquitously in neuroimaging, in disciplines ranging from magnetic resonance spectroscopy to temperature mapping and susceptibility-weighted imaging. Most B₀ maps are acquired using standard gradient-echo-based vendor-provided sequences, often comprised of two echoes spaced a few milliseconds apart. Herein, we analyze the optimal spacing of echo times, defined as those maximizing precision-minimizing the standard deviation-for a fixed total acquisition time. Field estimation is carried out using a weighted least squares estimator. The standard deviation is shown to be approximately inversely proportional to the total acquisition time, suggesting a law of diminishing returns, whereby substantial gains are obtained up to a certain point, with little improvement beyond that point. Validations are provided in a phantom and a group of volunteers. Multi-gradient echo sequences are readily available on all manufacturer platforms, making our recommendations straightforward to implement on any modern scanner.

Dynamic multicoil technique (DYNAMITE) MRI on human brain

PURPOSE

Spatial encoding for MRI is generally based on linear x, y, and z magnetic field gradients generated by a set of dedicated gradient coils. We recently introduced the dynamic multicoil technique (DYNAMITE) for B₀ field control and demonstrated DYNAMITE MRI in a preclinical MR environment. In this study, we report the first realization of DYNAMITE MRI of the in vivo human head.

METHODS

Gradient fields for DYNAMITE MRI were generated with a 28-channel multicoil hardware arranged in 4 rows of 7 coils on a cylindrical surface (length 359 mm, diameter 344 mm, maximum 5 A per coil). DYNAMITE MRIs of a resolution phantom and in vivo human heads were acquired with multislice gradient-echo, multislice spin-echo, and 3D gradient-echo sequences. The resultant image fidelity was compared to that obtained with conventional gradient coil technology.

RESULTS

DYNAMITE field control enabled the realization of all imaging sequences with average gradient errors ≤ 1%. DYNAMITE MRI provided image quality and sensitivity comparable to conventional gradient technology without any obvious artifacts. Some minor geometric deformations were noticed primarily in the image periphery as the result of regional field imperfections. The imperfections can be readily approximated theoretically through numerical integration of the Biot-Savart law and removed through image distortion correction.

CONCLUSION

The first realization of DYNAMITE MRI of the in vivo human head has been presented. The obtained image fidelity is comparable to MRI with conventional gradient coils, paving the way for full-fledged DYNAMITE MRI and B₀ shim systems for human applications.

Elevated homocarnosine and GABA in subject on isoniazid as assessed through 1H MRS at 7T

Typical magnetic resonance spectroscopy J-editing methods designed to quantify GABA suffer from contamination of both overlapping macromolecules and homocarnosine signal, introducing potential confounds. The aim of this study was to develop a novel method to assess accurately both the relative concentrations of homocarnosine as well as GABA free from overlapping creatine, homocarnosine and macromolecule signal. A novel method which utilized the combination of echo time STEAM and MEGA-sLASER magnetic resonance spectroscopy experiments at 7T were used to quantify the concentration of GABA and homocarnsoine independently, which are typically quantified in tandem. The metabolites GABA and homocarnosine were measured in brain of 6 healthy control subjects, and in a single subject medicated with isoniazid. It was found that (16.6±10.2)% of the supposed GABA signal in the brain originated from homocarnosine, and that isoniazid caused significantly elevated concentration of GABA and homocarnosine in a single subject compared to controls.

An orthogonal shim coil for 3T brain imaging

PURPOSE

We designed and implemented an orthogonal shim array consisting of shim coils with their planes perpendicular to the planes of neighboring RF coils. This shim coil improves the magnetic field homogeneity by minimizing the interference to RF coils.

METHODS

Using realistic off-resonance maps of the human brain, we first evaluated the performance of shim coils in different orientations. Based on simulations, we developed a 7-channel orthogonal shim array, whose coil plan was perpendicular to neighboring RF coils, at the forehead. A programmable open-source current driver supplied shim currents.

RESULTS

The 7-channel orthogonal shim array caused only marginal SNR loss to the integrated 32-channel RF-shim array. The 7-channel orthogonal shim array itself improved the magnetic field homogeneity by 30% in slice-optimized shimming, comparable to the baseline shimming offered by the scanner's 2nd order spherical harmonic shimming.

CONCLUSION

Orthogonal shim coils can improve the field homogeneity while maintaining high image SNR.

Switching Circuit Optimization for Matrix Gradient Coils

Matrix gradient coils with up to 84 coil elements were recently introduced for magnetic resonance imaging. Ideally, each element is driven by a dedicated amplifier, which may be technically and financially infeasible. Instead, several elements can be connected in series (called a “cluster”) and driven by a single amplifier. In previous works, a set of clusters, called a “configuration,” was sought to approximate a target field shape. Because a magnetic resonance pulse sequence requires several distinct field shapes, a mechanism to switch between configurations is needed. This can be achieved by a hypothetical switching circuit connecting all terminals of all elements with each other and with the amplifiers. For a predefined set of configurations, a switching circuit can be designed to require only a limited amount of switches. Here we introduce an algorithm to minimize the number of switches without affecting the ability of the configurations to accurately create the desired fields. The problem is modeled using graph theory and split into 2 sequential combinatorial optimization problems that are solved using simulated annealing. For the investigated cases, the results show that compared to unoptimized switching circuits, the reduction of switches in optimized circuits ranges from 8% to up to 44% (average of 31%). This substantial reduction is achieved without impeding circuit functionality. This study shows how technical effort associated with implementation and operation of a matrix gradient coil is related to different hardware setups and how to reduce this effort.

High-fidelity, high-isotropic-resolution diffusion imaging through gSlider acquisition with B 1 + and T1 corrections and integrated ΔB0 /Rx shim array

PURPOSE

B 1 + and T₁ corrections and dynamic multicoil shimming approaches were proposed to improve the fidelity of high-isotropic-resolution generalized slice-dithered enhanced resolution (gSlider) diffusion imaging.

METHODS

An extended reconstruction incorporating B 1 + inhomogeneity and T₁ recovery information was developed to mitigate slab-boundary artifacts in short-repetition time (TR) gSlider acquisitions. Slab-by-slab dynamic B₀ shimming using a multicoil integrated ΔB₀ /Rx shim array and high in-plane acceleration (R_inplane = 4) achieved with virtual-coil GRAPPA were also incorporated into a 1-mm isotropic resolution gSlider acquisition/reconstruction framework to achieve a significant reduction in geometric distortion compared to single-shot echo planar imaging (EPI).

RESULTS

The slab-boundary artifacts were alleviated by the proposed B 1 + and T₁ corrections compared to the standard gSlider reconstruction pipeline for short-TR acquisitions. Dynamic shimming provided >50% reduction in geometric distortion compared to conventional global second-order shimming. One-millimeter isotropic resolution diffusion data show that the typically problematic temporal and frontal lobes of the brain can be imaged with high geometric fidelity using dynamic shimming.

CONCLUSIONS

The proposed B 1 + and T₁ corrections and local-field control substantially improved the fidelity of high-isotropic-resolution diffusion imaging, with reduced slab-boundary artifacts and geometric distortion compared to conventional gSlider acquisition and reconstruction. This enabled high-fidelity whole-brain 1-mm isotropic diffusion imaging with 64 diffusion directions in 20 min using a 3T clinical scanner.

Ultralong TE In Vivo 1 H MR Spectroscopy of Omega-3 Fatty Acids in Subcutaneous Adipose Tissue at 7 T

BACKGROUND

Omega-3 (n-3) fatty acids (FA) play and important role in neural development and other metabolic diseases such as obesity and diabetes. The knowledge about the in vivo content and distribution of n-3 FA in human body tissues is not well established and the standard quantification of FA is invasive and costly.

PURPOSE

To detect omega-3 (n-3 CH3 ) and non-omega-3 (CH3 ) methyl group resonance lines with echo times up to 1200 msec, in oils, for the assessment of n-3 FA content, and the n-3 FA fraction in adipose tissue in vivo.

STUDY TYPE

Prospective technical development.

POPULATION

Three oils with different n-3 FA content and 24 healthy subjects.

FIELD STRENGTH/SEQUENCE

Single-voxel MR spectroscopy (SVS) with a point-resolved spectroscopy (PRESS) sequence with an echo time (TE) of 1000 msec at 7 T.

ASSESSMENT

Knowledge about the J-coupling evolution of both CH3 resonances was used for the optimal detection of the n-3 CH3 resonance line at a TE of 1000 msec. The accuracy of the method in oils and in vivo was validated from a biopsy sample with gas chromatography analysis.

STATISTICAL TESTS

SVS data were compared to gas chromatography with the Pearson correlation coefficient.

RESULTS

T2 relaxation times in oils were assessed as follows: CH2 , 65 ± 22 msec; CH3 , 325 ± 7 msec; and n-3 CH3 , 628 ± 34 msec. The n-3 FA fractions from oil phantom experiments (n = 3) were in agreement with chromatography analysis and the comparison of in vivo obtained data with the results of chromatography analysis (n = 5) yielded a significant correlation (P = 0.029).

DATA CONCLUSION

PRESS with ultralong-TE can detect and quantify the n-3 CH3 signal in vivo at 7 T.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:71-82.

Driving mutually coupled gradient array coils in magnetic resonance imaging

PURPOSE

In contrast to conventional linear gradients, gradient coil arrays with arbitrary spatial dependency might experience strong mutual coupling. Although conventional gradient power amplifiers with feedback loop might compensate the effect of coupling, required voltages for the compensation are generally unknown and has to be considered beforehand to ensure that amplifier voltage limits are not exceeded. A first-order circuit model is proposed to be used as a feedforward model which enables analytical formulas of required voltages to drive the mutually coupled gradient coil arrays.

THEORY AND METHODS

A first-order circuit model including the mutual couplings is provided to analytically calculate the input voltages and minimum achievable rise times for a given set of gradient array currents and amplifier limitations. Previously designed 9-channel Z-gradient coil array and home-built gradient amplifiers (50 V and 20 A) are used in the experiments. Three sets of currents optimized for linear Z-gradient, second-order Z2, and third-order Z3 fields are used in the bench-top experiments. The current weightings for the linear Z-gradient are also used as the readout gradient in the 3T MRI experiments.

RESULTS

Current measurements for the example magnetic field profiles with minimum rise times are demonstrated for the simultaneous use of 9-channel gradient coils and amplifiers. MRI experiments verify that a linear Z-gradient field with a desired time waveform can be generated using a mutually coupled array coils.

CONCLUSION

Bench-top and MRI experiments demonstrate the feasibility of the proposed circuit model and analytical formulas to drive the mutually coupled gradient coils.

Spin echo based cardiac diffusion imaging at 7T: An ex vivo study of the porcine heart at 7T and 3T

Purpose of this work was to assess feasibility of cardiac diffusion tensor imaging (cDTI) at 7 T in a set of healthy, unfixed, porcine hearts using various parallel imaging acceleration factors and to compare SNR and derived cDTI metrics to a reference measured at 3 T. Magnetic resonance imaging was performed on 7T and 3T whole body systems using a spin echo diffusion encoding sequence with echo planar imaging readout. Five reference (b = 0 s/mm2) images and 30 diffusion directions (b = 700 s/mm2) were acquired at both 7 T and 3 T using a GRAPPA acceleration factor R = 1. Scans at 7 T were repeated using R = 2, R = 3, and R = 4. SNR evaluation was based on 30 reference (b = 0 s/mm2) images of 30 slices of the left ventricle and cardiac DTI metrics were compared within AHA segmentation. The number of hearts scanned at 7 T and 3 T was n = 11. No statistically significant differences were found for evaluated helix angle, secondary eigenvector angle, fractional anisotropy and apparent diffusion coefficient at the different field strengths, given sufficiently high SNR and geometrically undistorted images. R≥3 was needed to reduce susceptibility induced geometric distortions to an acceptable amount. On average SNR in myocardium of the left ventricle was increased from 29±3 to 44±6 in the reference image (b = 0 s/mm2) when switching from 3 T to 7 T. Our study demonstrates that high resolution, ex vivo cDTI is feasible at 7 T using commercial hardware.

Distinguishing Lipid Subtypes by Amplifying Contrast from J-Coupling

Previous work has highlighted the complicated and distinctive dynamics that set signal evolution during a train of spin echoes, especially with nonuniform echo spacing applied to complex molecules like fats. The work presented here regards those signal patterns as codes that can be used as a contrast mechanism, capable of distinguishing mixtures of molecules with an imaging sequence, sidestepping many challenges of spectroscopy. For particular arrays of echo spacings, non-monotonic and distinctive signal evolution can be enhanced to improve contrast between target species. This work presents simulations that show how contrast between two molecules: (a) depends on the specific sequence of echo spacing, (b) is directly linked to the presence of J-coupling, and (c) can be relatively insensitive to variations in B0, T2 and B1. Imaging studies with oils demonstrate this phenomenon experimentally and also show that spin echo codes can be used for quantification. Finally, preliminary experiments apply the method to human liver in vivo, verifying that the presence of fat can lead to nonmonotonic codes like those seen in vitro. In summary, nonuniformly spaced echo trains introduce a new approach to molecular imaging of J-coupled species, such as lipids, which may have implications diagnosing metabolic diseases.

Combined imaging and shimming with the dynamic multi-coil technique

PURPOSE

Spatial encoding and shimming in MRI have traditionally been performed using dedicated coils that generate orthogonal spherical harmonic fields. The recently introduced multi-coil hardware has proven that MRI-relevant magnetic fields can also be created by a generic set of localized coils producing non-orthogonal fields. As a step towards establishing a purely multi-coil-based MRI field generation system, the feasibility of performing conventional Cartesian k-space encoding and echo-planar imaging (EPI), as well as concurrent encoding and shimming is demonstrated in this study.

METHODS

We report the use of Dynamic Multi-Coil Technique (DYNAMITE) for combined Cartesian encoding and shimming, and EPI using a 48-channel multi-coil system. Experiments were performed on phantom objects and biological specimens in a 9.4 T pre-clinical scanner. Cartesian Fourier-encoded MRI and EPI were implemented whereby the magnetic fields required for encoding of the three orthogonal spatial dimensions were entirely based on linear combinations of multi-coil fields. Furthermore, DYNAMITE imaging was augmented by concurrent DYNAMITE shimming with the same hardware.

RESULTS

DYNAMITE-based MR and echo-planar images were indistinguishable from those acquired with the conventional linear imaging gradients provided by the scanner. In experiments with concurrent DYNAMITE shimming and imaging, shim challenges that would result in extreme spatial distortion and signal loss were corrected very effectively with more than 92% signal recovery in case of extreme Z² shim challenge that resulted in complete signal dephasing in most slices.

CONCLUSIONS

We demonstrate the first successful implementation of combined DYNAMITE imaging and shimming and show the feasibility of performing EPI with DYNAMITE hardware. Our results substantiate the potential of multi-coil hardware as a full-fledged imaging and shimming system, with additional potential benefits of reduced echo-time and risk of peripheral nerve stimulation while performing EPI.

Improved methods for MRI-compatible implants in nonhuman primates

•

We created custom PEEK implants coated with hydroxylapatite to promote osseointegration.

•

Headposts and chambers were implanted with ceramic screws and the surgical incision was closed subcutaneously.

•

We prevented the animal from picking and scratching after surgery by using a head cap that protects the wound margin.

•

Implants integrated with the skull and remained robust after a year with no growth of granulation tissue.

•

MRI signal and contrast-to-noise ratio improved in NHPs implanted with our new methods.

Background

Neuroscientists commonly use permanently implanted headposts to stabilize the head of nonhuman primates (NHPs) during electrophysiology and functional magnetic resonance imaging (fMRI). Here, we present improved methodology for MRI-compatible implants without the use of acrylic for head stabilization in NHPs.

New method

MRI is used to obtain a 3D-reconstruction of NHP skulls, which are used to create customized implants by modeling intersections with the bone. Implants are manufactured from PEEK using computer numerical control machining and coated with hydroxyapatite to promote osseointegration. Surgically, implants are attached to the skull with ceramic screws, while the skin flap is pulled over the implant and closed subcutaneously.

Results

Quality of blood oxygen level dependent (BOLD) fMRI signal is improved in animals implanted with our method as compared to traditional acrylic implants. Additionally, implants are well-integrated with the skull, remain robust for more than a year and without granulation tissue around the skin margin.

Comparison with existing method(s)

Previous improvements on NHP implants (Chen et al., 2017; McAndrew et al., 2012; Mulliken et al., 2015; Overton et al., 2017) lacked fMRI-compatibility, as they relied on titanium headposts and/or titanium screws. Thus, most fMRI studies in NHPs today still rely on the use of acrylic-based headposts for stabilization and the use of contrast-enhanced agents to improve MRI signal.

Conclusions

Our method preserves fMRI-compatibility and results in measurable improvement in BOLD signal without the use of contrast-enhanced agents. Furthermore, the long-term stability of our implants contributes positively to the wellbeing of NHPs in neuroscience research.

In vivo B0 field shimming methods for MRI at 7T

Functional MRI (fMRI) at 7T and above provides improved Signal-to-Noise Ratio and Contrast-to-Noise Ratio compared to 3T acquisitions. In addition to the beneficial effects on spin polarization and magnetization of deoxyhemoglobin, the increased applied field also further magnetizes air and tissue. While the magnets themselves typically provide a static B0 field with sufficient spatial homogeneity, the diamagnetism of tissue and the paramagnetism of air causes local field deviations inside the human head. These spatially-varying field offsets (ΔB0) cause image artifacts, especially in single shot EPI, including geometric distortion, signal dropout, and blurring. These effects are particularly strong near air-tissue interfaces such as the frontal sinus, and ear canals. Furthermore, if the field offsets are dynamically modulated through physiological processes such as respiration or motion, then the effect on the image time-series can be even more problematic. While post-processing methods have been developed to mitigate these effects, the ideal solution would be to reduce the ΔB0 variations at their source. Typically 7T scanners contain 2nd and some 3rd order spherical harmonic shim coil terms to cancel static ΔB0 variations of low spatial order. In this article, we will motivate the need for improved, higher-order compensation for B0 inhomogeneity and potentially add dynamic control of these fields. We discuss and compare several promising hardware approaches for static and dynamic B0 shimming using either higher-order spherical harmonic shim coils or multi-coil shim arrays as well as passive shimming approaches, and active variants such and adaptive current networks.

Comparision of new element designs for combined RF-Shim arrays at 7 T

PURPOSE

To identify novel concepts for RF-shim loop architectures suitable for 7T made of two concentric conducting loops fulfilling RF and DC functions, respectively, and to determine their relative SNR performance. The goal is to minimize interference between the two systems while making efficient use of the space closest to the body.

THEORY

We show by means of theoretical derivation of the frequency spectrum that the proposed two-loop structure exhibits an anti-resonant null and a resonant peak in the frequency domain.

METHODS

The proposed structure is comprised of two concentric wire loops either arranged as nested loops or in the form of a coaxial cable, in which the two conductors carry the RF and shim signals, respectively. We use theory, simulation, and phantom measurements to obtain frequency spectra and SNR maps for the proposed structures.

RESULTS

Retained SNR is found to be 75% in the coaxial loop and ranges from 57% to 67% in three different coaxial configurations. We have found both implementations to be a viable concept for the use in RF-shim devices if remaining SNR limitations can be overcome.

CONCLUSIONS

We have investigated two new design modalities in 7T RF-shim coil design that separate the RF and shim conductors such that the previously proposed toroidal chokes are eliminated - thereby removing undesired additional loss, bulk, and design complexity.

Spontaneous activity forms a foundation for odor-evoked activation maps in the rat olfactory bulb

Fluctuations in spontaneous activity have been observed by many neuroimaging techniques, but because these resting-state changes are not evoked by stimuli, it is difficult to determine how they relate to task-evoked activations. We conducted multi-modal neuroimaging scans of the rat olfactory bulb, both with and without odor, to examine interaction between spontaneous and evoked activities. Independent component analysis of spontaneous fluctuations revealed resting-state networks, and odor-evoked changes revealed activation maps. We constructed simulated activation maps using resting-state networks that were highly correlated to evoked activation maps. Simulated activation maps derived by intrinsic optical signal (IOS), which covers the dorsal portion of the glomerular sheet, significantly differentiated one odor's evoked activation map from the other two. To test the hypothesis that spontaneous activity of the entire glomerular sheet is relevant for representing odor-evoked activations, we used functional magnetic resonance imaging (fMRI) to map the entire glomerular sheet. In contrast to the IOS results, the fMRI-derived simulated activation maps significantly differentiated all three odors' evoked activation maps. Importantly, no evoked activation maps could be significantly differentiated using simulated activation maps produced using phase-randomized resting-state networks. Given that some highly organized resting-state networks did not correlate with any odors' evoked activation maps, we posit that these resting-state networks may characterize evoked activation maps associated with odors not studied. These results emphasize that fluctuations in spontaneous activity form a foundation for active processing, signifying the relevance of resting-state mapping to functional neuroimaging.

Optimization of Coil Element Configurations for a Matrix Gradient Coil

Recently, matrix gradient coils (also termed multi-coils or multi-coil arrays) were introduced for imaging and B₀ shimming with 24, 48, and even 84 coil elements. However, in imaging applications, providing one amplifier per coil element is not always feasible due to high cost and technical complexity. In this simulation study, we show that an 84-channel matrix gradient coil (head insert for brain imaging) is able to create a wide variety of field shapes even if the number of amplifiers is reduced. An optimization algorithm was implemented that obtains groups of coil elements, such that a desired target field can be created by driving each group with an amplifier. This limits the number of amplifiers to the number of coil element groups. Simulated annealing is used due to the NP-hard combinatorial nature of the given problem. A spherical harmonic basis set up to the full third order within a sphere of 20-cm diameter in the center of the coil was investigated as target fields. We show that the median normalized least squares error for all target fields is below approximately 5% for 12 or more amplifiers. At the same time, the dissipated power stays within reasonable limits. With a relatively small set of amplifiers, switches can be used to sequentially generate spherical harmonics up to third order. The costs associated with a matrix gradient coil can be lowered, which increases the practical utility of matrix gradient coils.

A z-gradient array for simultaneous multi-slice excitation with a single-band RF pulse

PURPOSE

Multi-slice radiofrequency (RF) pulses have higher specific absorption rates, more peak RF power, and longer pulse durations than single-slice RF pulses. Gradient field design techniques using a z-gradient array are investigated for exciting multiple slices with a single-band RF pulse.

THEORY AND METHODS

Two different field design methods are formulated to solve for the required current values of the gradient array elements for the given slice locations. The method requirements are specified, optimization problems are formulated for the minimum current norm and an analytical solution is provided. A 9-channel z-gradient coil array driven by independent, custom-designed gradient amplifiers is used to validate the theory.

RESULTS

Performance measures such as normalized slice thickness error, gradient strength per unit norm current, power dissipation, and maximum amplitude of the magnetic field are provided for various slice locations and numbers of slices. Two and 3 slices are excited by a single-band RF pulse in simulations and phantom experiments.

CONCLUSION

The possibility of multi-slice excitation with a single-band RF pulse using a z-gradient array is validated in simulations and phantom experiments. Magn Reson Med 80:400-412, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

The public multi-coil information (PUMCIN) policy

PURPOSE

Multi-coil (MC) magnetic field modeling has emerged as a viable alternative to conventional field generation based on spherical harmonic shapes, and an active MC community is forming. Although all MC applications share the same modeling concept, the specific MC designs can largely differ as a result of disparities in region of interest (eg, human versus rodent), intended MR application (eg, B0 shimming versus spatial encoding), or other experimental constraints (eg, available bore space or integration with radiofrequency technology). To date, a lack of detailed information on existing MC designs complicates the assessment and precludes a meaningful comparison.

METHODS

Here, we suggest that future publications involving the MC technique not only report the benefits for the application at hand, but also include an explicit description of the MC wire pattern used.

RESULTS

This public multi-coil information (PUMCIN) policy represents a voluntary commitment to promoting free public access to the details necessary for reproducing and benefiting from MC research.

CONCLUSIONS

The PUMCIN policy is expected to initiate a paradigm shift with respect to the way MC innovation is reported. By setting an example, we hope to encourage the evolving MC community to maximize the benefits for science and society by embracing it. Magn Reson Med 78:2042-2047, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Simultaneous use of linear and nonlinear gradients for B1+ inhomogeneity correction

The simultaneous use of linear spatial encoding magnetic fields (L-SEMs) and nonlinear spatial encoding magnetic fields (N-SEMs) in B₁⁺ inhomogeneity problems is formulated and demonstrated with both simulations and experiments. Independent excitation k-space variables for N-SEMs are formulated for the simultaneous use of L-SEMs and N-SEMs by assuming a small tip angle. The formulation shows that, when N-SEMs are considered as an independent excitation k-space variable, numerous different k-space trajectories and frequency weightings differing in dimension, length, and energy can be designed for a given target transverse magnetization distribution. The advantage of simultaneous use of L-SEMs and N-SEMs is demonstrated by B₁⁺ inhomogeneity correction with spoke excitation. To fully utilize the independent k-space formulations, global optimizations are performed for 1D, 2D RF power limited, and 2D RF power unlimited simulations and experiments. Three different cases are compared: L-SEMs alone, N-SEMs alone, and both used simultaneously. In all cases, the simultaneous use of L-SEMs and N-SEMs leads to a decreased standard deviation in the ROI compared with using only L-SEMs or N-SEMs. The simultaneous use of L-SEMs and N-SEMs results in better B₁⁺ inhomogeneity correction than using only L-SEMs or N-SEMs due to the increased number of degrees of freedom.

Modeling real shim fields for very high degree (and order) B0 shimming of the human brain at 9.4 T

PURPOSE

To describe the process of calibrating a B₀ shim system using high-degree (or high order) spherical harmonic models of the measured shim fields, to provide a method that considers amplitude dependency of these models, and to show the advantage of very high-degree B₀ shimming for whole-brain and single-slice applications at 9.4 Tesla (T).

METHODS

An insert shim with up to fourth and partial fifth/sixth degree (order) spherical harmonics was used with a Siemens 9.4T scanner. Each shim field was measured and modeled as input for the shimming algorithm. Optimal shim currents can therefore be calculated in a single iteration. A range of shim currents was used in the modeling to account for possible amplitude nonlinearities. The modeled shim fields were used to compare different degrees of whole-brain B₀ shimming on healthy subjects.

RESULTS

The ideal shim fields did not correctly shim the subject brains. However, using the modeled shim fields improved the B₀ homogeneity from 55.1 (second degree) to 44.68 Hz (partial fifth/sixth degree) on the whole brains of 9 healthy volunteers, with a total applied current of 0.77 and 6.8 A, respectively.

CONCLUSIONS

The necessity of calibrating the shim system was shown. Better B₀ homogeneity drastically reduces signal dropout and distortions for echo-planar imaging, and significantly improves the linewidths of MR spectroscopy imaging. Magn Reson Med 79:529-540, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Gradient-echo EPI using a high-degree shim insert coil at 7 T: Implications for BOLD fMRI

Purpose

To quantitatively assess the effects of high degree and order (1^st–4^th+) relative to 1^st–2^nd degree B₀ shimming at 7 Tesla (T) on gradient‐echo echo planar imaging (GE‐EPI) and blood‐oxygen‐level dependent (BOLD) activation.

Methods

Simulations and GE‐EPI were performed at (2mm)³ and (3mm)³ resolution, evaluating the temporal signal‐to‐noise ratio (tSNR), transverse relaxivity ( R2*), BOLD % signal change and activated pixel counts in a breath‐hold task.

Results

Comparing the 1^st–4^th+ degree with 1^st–2^nd degree shimmed B₀ maps generated spatially varying regions of Δ|B0|=|B01−2|−|B01−4+|. As binned in 10‐Hz intervals, the two center Δ|B₀| (±10 Hz) bins maintained the B₀ offset of 48.6% of gray‐matter pixels. In the positive Δ|B₀| bins greater than 10 Hz, the 1^st–4^th+degree shimming improved the B₀ offset in 41.1%; in negative Δ|B₀| bins less than −10 Hz, the offset worsened in 10.2% of the pixels. In the positive Δ|B₀| bins, we found variable but significant increases in BOLD sensitivity; the negative Δ|B₀| bins showed significant decreases. In the breath‐hold studies, positive bins showed significantly increased activated pixel numbers (+5–29%), whereas negative bins showed −18 to 0% decline.

Conclusion

1^st–4^th+ degree shimming maintained B₀ homogeneity over central brain regions while improving most of the other regions, including the inferior frontal lobe. Magn Reson Med 78:1734–1745, 2017. © 2016 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 NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

The technological future of 7 T MRI hardware

In this article we present our projections of future hardware developments on 7 T human MRI systems. These include compact cryogen-light magnets, improved gradient performance, integrated RF-receive and direct current shimming coil arrays, new RF technology with adaptive impedance matching, patient-specific specific absorption rate estimation and monitoring, and increased integration of physiological monitoring systems. Copyright © 2015 John Wiley & Sons, Ltd.

A 24-channel shim array for the human spinal cord: Design, evaluation, and application

PURPOSE

A novel multichannel shim array is introduced to improve MRI and spectroscopic studies of the human spinal cord.

METHODS

Twenty-four-channel shim and 8-channel transceiver arrays were designed to insert into the patient bed table to lie in close proximity to the subject's spine. The reference field patterns of each of the shim channels (Hz/A) were determined empirically via gradient echo field mapping and subsequently used to demonstrate shim performance at 3 Tesla using an ex vivo phantom, which incorporated a fixed human spine. The shim was further demonstrated on five healthy volunteers.

RESULTS

Application of the shim to the ex vivo phantom reduced the standard deviation of the field over the spinal volume of interest (123.4 cm³ ) from an original 51.3 Hz down to 32.5 Hz, amounting to an improvement in field homogeneity of 36.6%. In vivo, the spine shim resulted in an average improvement in field homogeneity of 63.8 ± 15.4%.

CONCLUSION

The localized spine shim offers a promising new means of correcting magnetic field distortion in the spinal cord. Magn Reson Med 76:1604-1611, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Integrated RF-shim coil allowing two degrees of freedom shim current

High-quality magnetic resonance imaging and spectroscopic measurements require a highly homogeneous magnetic field. Different from global shimming, localized off-resonance can be corrected by using multi-coil shimming. Previously, integrated RF and shimming coils have been used to implement multi-coil shimming. Such coils share the same conductor for RF signal reception and shim field generation. Here we propose a new design of the integrated RF-shim coil at 3-tesla, where two independent shim current paths are allowed in each coil. This coil permits a higher degree of freedom in shim current distribution design. We use both phantom experiments and simulations to demonstrate the feasibility of this new design.

Reproducibility measurement of glutathione, GABA, and glutamate: Towards in vivo neurochemical profiling of multiple sclerosis with MR spectroscopy at 7T

PURPOSE

To determine the reproducibility of a comprehensive single-session measurement of glutathione (GSH), γ-aminobutyric acid (GABA), glutamate, and other biochemicals implicated in the pathophysiology of multiple sclerosis (MS) in the human brain with ¹ H magnetic resonance spectroscopy (MRS).

MATERIALS AND METHODS

Five healthy subjects were studied twice in separate 1-hour sessions at 7T. One MS patient was also scanned once. GSH and GABA were measured with J-difference editing using a semilocalized by adiabatic selective refocusing sequence (semi-LASER, TE = 72 msec). A stimulated echo acquisition mode sequence (STEAM, TE = 10 msec) was used to detect glutamate along with the overall biochemical profile. Spectra were quantified with LCModel. Quantification accuracy was assessed through Cramer-Rao lower bounds (CRLB). Reproducibility of the metabolite quantification was tested using coefficients of variation (CoV).

RESULTS

CRLB were ≤7% for GSH, GABA, and glutamate and average CoV of 7.8 ± 3.2%, 9.5 ± 7.0%, and 3.2 ± 1.7% were achieved, respectively. The average test/retest concentration differences at this measurement reproducibility and quantification accuracy were smaller for GABA and glutamate than intersubject variations in metabolite content with CoV ratios of 0.6 and 0.8, respectively. As proof of principle, GSH, GABA, and glutamate were also detected in an MS patient.

CONCLUSION

GSH, GABA, glutamate, and other metabolites relevant in MS can be quantified at 7T with high accuracy and reproducibility in a single 1-hour session. This methodology might serve as a clinical research tool to investigate biochemical markers associated with MS.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:187-198.