Submillimeter balanced SSFP BOLD-functional MRI accelerated with 3D stack-of-spirals at 9.4 T.

PURPOSE: This work aims to improve the speed of balanced SSFP (bSSFP) acquisition with segmented 3D stack-of-spirals for functional brain studies at ultrahigh field. METHODS: Functional experiments were performed with an accelerated 3D stack-of-spirals sequence with water excitation for fat suppression. The resulting data were reconstructed using an iterative algorithm with corrections for system imperfections such as trajectory deviations and B0 inhomogeneity. In the first set of experiments, we evaluated the signal change and stability with respect to echo and TR for a full-field checkerboard stimulus. To demonstrate the high spatio-temporal resolution of the developed method, the results of three optimized protocols at submillimeter resolution (0.6-mm isotropic and 0.8-mm isotropic) and at 1.2 mm isotropic resolution for whole-brain coverage were shown. RESULTS: Water excitation and the model-based iterative reconstruction improved image quality. The BOLD-related signal changes increased with longer TE and longer TR. We observed an increase in thermal noise performance at lower TE and higher TR. However, signal stability deteriorates at higher TE and TR. Therefore, optimized protocols used shorter TE and moderately long TR to maximize the sensitivity and speed. Reproducible activations were detected along the gray-matter gyri in the submillimeter protocols with a median signal change of approximately 4% across subjects. CONCLUSIONS: Three-dimensional stack-of-spirals enables passband balanced SSFP functional imaging at a much higher spatial and temporal scale, compared with conventional spoiled gradient-echo train sequences.

Intra-scan RF power amplifier drift correction.

PURPOSE: The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated. METHODS: RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan. RESULTS: RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability. CONCLUSION: Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Improving MR axon radius estimation in human white matter using spiral acquisition and field monitoring.

PURPOSE: To compare MR axon radius estimation in human white matter using a multiband spiral sequence combined with field monitoring to the current state-of-the-art echo-planar imaging (EPI)-based approach. METHODS: A custom multiband spiral sequence was used for diffusion-weighted imaging at ultra-high b $$ b $$ -values. Field monitoring and higher order image reconstruction were employed to greatly reduce artifacts in spiral images. Diffusion weighting parameters were chosen to match a state-of-the art EPI-based axon radius mapping protocol. The spiral approach was compared to the EPI approach by comparing the image signal-to-noise ratio (SNR) and performing a test-retest study to assess the respective variability and repeatability of axon radius mapping. Effective axon radius estimates were compared over white matter voxels and along the left corticospinal tract. RESULTS: Increased SNR and reduced artifacts in spiral images led to reduced variability in resulting axon radius maps, especially in low-SNR regions. Test-retest variability was reduced by a factor of approximately 1.5 using the spiral approach. Reduced repeatability due to significant bias was found for some subjects in both spiral and EPI approaches, and attributed to scanner instability, pointing to a previously unknown limitation of the state-of-the-art approach. CONCLUSION: Combining spiral readouts with field monitoring improved mapping of the effective axon radius compared to the conventional EPI approach.

The possible influence of third-order shim coils on gradient-magnet interactions: an inter-field and inter-site study.

OBJECTIVE: To assess the possible influence of third-order shim coils on the behavior of the gradient field and in gradient-magnet interactions at 7 T and above. MATERIALS AND METHODS: Gradient impulse response function measurements were performed at 5 sites spanning field strengths from 7 to 11.7 T, all of them sharing the same exact whole-body gradient coil design. Mechanical fixation and boundary conditions of the gradient coil were altered in several ways at one site to study the impact of mechanical coupling with the magnet on the field perturbations. Vibrations, power deposition in the He bath, and field dynamics were characterized at 11.7 T with the third-order shim coils connected and disconnected inside the Faraday cage. RESULTS: For the same whole-body gradient coil design, all measurements differed greatly based on the third-order shim coil configuration (connected or not). Vibrations and gradient transfer function peaks could be affected by a factor of 2 or more, depending on the resonances. Disconnecting the third-order shim coils at 11.7 T also suppressed almost completely power deposition peaks at some frequencies. DISCUSSION: Third-order shim coil configurations can have major impact in gradient-magnet interactions with consequences on potential hardware damage, magnet heating, and image quality going beyond EPI acquisitions.

Clinically Meaningful Magnetic Resonance Endpoints Sensitive to Preataxic Spinocerebellar Ataxia Types 1 and 3.

OBJECTIVE: This study was undertaken to identify magnetic resonance (MR) metrics that are most sensitive to early changes in the brain in spinocerebellar ataxia type 1 (SCA1) and type 3 (SCA3) using an advanced multimodal MR imaging (MRI) protocol in the multisite trial setting. METHODS: SCA1 or SCA3 mutation carriers and controls (n = 107) underwent MR scanning in the US-European READISCA study to obtain structural, diffusion MRI, and MR spectroscopy data using an advanced protocol at 3T. Morphometric, microstructural, and neurochemical metrics were analyzed blinded to diagnosis and compared between preataxic SCA (n = 11 SCA1, n = 28 SCA3), ataxic SCA (n = 14 SCA1, n = 37 SCA3), and control (n = 17) groups using nonparametric testing accounting for multiple comparisons. MR metrics that were most sensitive to preataxic abnormalities were identified using receiver operating characteristic (ROC) analyses. RESULTS: Atrophy and microstructural damage in the brainstem and cerebellar peduncles and neurochemical abnormalities in the pons were prominent in both preataxic groups, when patients did not differ from controls clinically. MR metrics were strongly associated with ataxia symptoms, activities of daily living, and estimated ataxia duration. A neurochemical measure was the most sensitive metric to preataxic changes in SCA1 (ROC area under the curve [AUC] = 0.95), and a microstructural metric was the most sensitive metric to preataxic changes in SCA3 (AUC = 0.92). INTERPRETATION: Changes in cerebellar afferent and efferent pathways underlie the earliest symptoms of both SCAs. MR metrics collected with a harmonized advanced protocol in the multisite trial setting allow detection of disease effects in individuals before ataxia onset with potential clinical trial utility for subject stratification. ANN NEUROL 2022.

Reproducibility of rapid multi-parameter mapping at 3T and 7T with highly segmented and accelerated 3D-EPI.

PURPOSE: Quantitative multi-parameter mapping (MPM) has been shown to provide good longitudinal and cross-sectional reproducibility for clinical research. Unfortunately, acquisition times (TAs) are typically infeasible for routine scanning at high resolutions. METHODS: A fast whole-brain MPM protocol based on interleaved multi-shot 3D-EPI with controlled aliasing (SC-EPI) at 3T and 7T is proposed and compared with MPM using a standard spoiled gradient echo (FLASH) sequence. Four parameters (R1 , PD, R 2 * $$ {R}_2^{\ast } $$ , and MTsat) were measured in less than 3 min at 1 mm isotropic resolution. Five subjects went through the same scanning sessions twice at each scanner. The intra-subject coefficient of variation (scan-rescan) (CoV) was estimated for each protocol and scanner to assess the longitudinal reproducibility. RESULTS: At 3T, the CoV of SC-EPI ranged between 1.2%-4.8% for PD and R1 , 2.8%-10.6% for R 2 * $$ {R}_2^{\ast } $$ and MTsat, which was comparable with FLASH (0.6%-4.9% for PD and R1 , 2.6%-11.3% for R 2 * $$ {R}_2^{\ast } $$ and MTsat). At 7T, where the SC-EPI TA was reduced to ∼2 min, the CoV of SC-EPI (1.4%-10.6% for PD, R1 , and R 2 * $$ {R}_2^{\ast } $$ ) was 1.2-2.4 times larger than the CoV of FLASH (1.0%-15%) and MTsat showed much higher variability across subjects. The SC-EPI-MPM protocol at 3T showed high reproducibility and yielded stable quantitative maps at a clinically feasible resolution and scan time, whereas at 7T, MT saturation homogeneity needs to be improved. CONCLUSION: SC-EPI-based MPM is feasible as an additional MRI modality in clinical or population studies where the parameters offer great potential as biomarkers.

Open-source MR imaging and reconstruction workflow.

PURPOSE: This work presents an end-to-end open-source MR imaging workflow. It is highly flexible in rapid prototyping across the whole imaging process and integrates vendor-independent openly available tools. The whole workflow can be shared and executed on different MR platforms. It is also integrated in the JEMRIS simulation framework, which makes it possible to generate simulated data from the same sequence that runs on the MRI scanner using the same pipeline for image reconstruction. METHODS: MRI sequences can be designed in Python or JEMRIS using the Pulseq framework, allowing simplified integration of new sequence design tools. During the sequence design process, acquisition metadata required for reconstruction is stored in the MR raw data format. Data acquisition is possible on MRI scanners supported by Pulseq and in simulations through JEMRIS. An image reconstruction and postprocessing pipeline was implemented into a Python server that allows real-time processing of data as it is being acquired. The Berkeley Advanced Reconstruction Toolbox is integrated into this framework for image reconstruction. The reconstruction pipeline supports online integration through a vendor-dependent interface. RESULTS: The flexibility of the workflow is demonstrated with different examples, containing 3D parallel imaging with controlled aliasing in volumetric parallel imaging (CAIPIRINHA) acceleration, spiral imaging, and B0 mapping. All sequences, data, and the corresponding processing pipelines are publicly available. CONCLUSION: The proposed workflow is highly flexible and allows integration of advanced tools at all stages of the imaging process. All parts of this workflow are open-source, simplifying collaboration across different MR platforms or sites and improving reproducibility of results.

Combining navigator and optical prospective motion correction for high-quality 500 µm resolution quantitative multi-parameter mapping at 7T.

PURPOSE: High-resolution quantitative multi-parameter mapping shows promise for non-invasively characterizing human brain microstructure but is limited by physiological artifacts. We implemented corrections for rigid head movement and respiration-related B0-fluctuations and evaluated them in healthy volunteers and dementia patients. METHODS: Camera-based optical prospective motion correction (PMC) and FID navigator correction were implemented in a gradient and RF-spoiled multi-echo 3D gradient echo sequence for mapping proton density (PD), longitudinal relaxation rate (R1) and effective transverse relaxation rate (R2*). We studied their effectiveness separately and in concert in young volunteers and then evaluated the navigator correction (NAVcor) with PMC in a group of elderly volunteers and dementia patients. We used spatial homogeneity within white matter (WM) and gray matter (GM) and scan-rescan measures as quality metrics. RESULTS: NAVcor and PMC reduced artifacts and improved the homogeneity and reproducibility of parameter maps. In elderly participants, NAVcor improved scan-rescan reproducibility of parameter maps (coefficient of variation decreased by 14.7% and 11.9% within WM and GM respectively). Spurious inhomogeneities within WM were reduced more in the elderly than in the young cohort (by 9% vs. 2%). PMC increased regional GM/WM contrast and was especially important in the elderly cohort, which moved twice as much as the young cohort. We did not find a significant interaction between the two corrections. CONCLUSION: Navigator correction and PMC significantly improved the quality of PD, R1, and R2* maps, particularly in less compliant elderly volunteers and dementia patients.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile.

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination Ï´) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination Ï´ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.

Automated olfactory bulb segmentation on high resolutional T2-weighted MRI.

The neuroimage analysis community has neglected the automated segmentation of the olfactory bulb (OB) despite its crucial role in olfactory function. The lack of an automatic processing method for the OB can be explained by its challenging properties (small size, location, and poor visibility on traditional MRI scans). Nonetheless, recent advances in MRI acquisition techniques and resolution have allowed raters to generate more reliable manual annotations. Furthermore, the high accuracy of deep learning methods for solving semantic segmentation problems provides us with an option to reliably assess even small structures. In this work, we introduce a novel, fast, and fully automated deep learning pipeline to accurately segment OB tissue on sub-millimeter T2-weighted (T2w) whole-brain MR images. To this end, we designed a three-stage pipeline: (1) Localization of a region containing both OBs using FastSurferCNN, (2) Segmentation of OB tissue within the localized region through four independent AttFastSurferCNN - a novel deep learning architecture with a self-attention mechanism to improve modeling of contextual information, and (3) Ensemble of the predicted label maps. For this work, both OBs were manually annotated in a total of 620 T2w images for training (n=357) and testing. The OB pipeline exhibits high performance in terms of boundary delineation, OB localization, and volume estimation across a wide range of ages in 203 participants of the Rhineland Study (Dice Score (Dice): 0.852, Volume Similarity (VS): 0.910, and Average Hausdorff Distance (AVD): 0.215 mm). Moreover, it also generalizes to scans of an independent dataset never encountered during training, the Human Connectome Project (HCP), with different acquisition parameters and demographics, evaluated in 30 cases at the native 0.7 mm HCP resolution (Dice: 0.738, VS: 0.790, and AVD: 0.340 mm), and the default 0.8 mm pipeline resolution (Dice: 0.782, VS: 0.858, and AVD: 0.268 mm). We extensively validated our pipeline not only with respect to segmentation accuracy but also to known OB volume effects, where it can sensitively replicate age effects (ß=-0.232, p<.01). Furthermore, our method can analyze a 3D volume in less than a minute (GPU) in an end-to-end fashion, providing a validated, efficient, and scalable solution for automatically assessing OB volumes.

The traveling heads 2.0: Multicenter reproducibility of quantitative imaging methods at 7 Tesla.

OBJECT: This study evaluates inter-site and intra-site reproducibility at ten different 7 T sites for quantitative brain imaging. MATERIAL AND METHODS: Two subjects - termed the "traveling heads" - were imaged at ten different 7 T sites with a harmonized quantitative brain MR imaging protocol. In conjunction with the system calibration, MP2RAGE, QSM, CEST and multi-parametric mapping/relaxometry were examined. RESULTS: Quantitative measurements with MP2RAGE showed very high reproducibility across sites and subjects, and errors were in concordance with previous results and other field strengths. QSM had high inter-site reproducibility for relevant subcortical volumes. CEST imaging revealed systematic differences between the sites, but reproducibility was comparable to results in the literature. Relaxometry had also very high agreement between sites, but due to the high sensitivity, differences caused by different applications of the B1 calibration of the two RF coil types used were observed. CONCLUSION: Our results show that quantitative brain imaging can be performed with high reproducibility at 7 T and with similar reliability as found at 3 T for multicenter studies of the supratentorial brain.

Impact of prospective motion correction, distortion correction methods and large vein bias on the spatial accuracy of cortical laminar fMRI at 9.4 Tesla.

Functional imaging with sub-millimeter spatial resolution is a basic requirement for assessing functional MRI (fMRI) responses across different cortical depths and is used extensively in the emerging field of laminar fMRI. Such studies seek to investigate the detailed functional organization of the brain and may develop to a new powerful tool for human neuroscience. However, several studies have shown that measurement of laminar fMRI responses can be biased by the image acquisition and data processing strategies. In this work, measurements with three different gradient-echo EPI BOLD fMRI protocols with a voxel size down to 650 â€‹µm isotropic were performed at 9.4 â€‹T. We estimated how prospective motion correction can help to improve spatial accuracy by reducing the number of spatial resampling steps in postprocessing. In addition, we demonstrate key requirements for accurate geometric distortion correction to ensure that distortion correction maps are properly aligned to the functional data and that strong variations of distortions near large veins can lead to signal overlays which cannot be corrected for during postprocessing. Furthermore, this study illustrates the spatial extent of bias induced by pial and other larger veins in laminar BOLD experiments. Since these issues under investigation affect studies performed with more conventional spatial resolutions, the methods applied in this work may also help to improve the understanding of the BOLD signal more broadly.

Whole brain snapshot CEST at 3T using 3D-EPI: Aiming for speed, volume, and homogeneity.

PURPOSE: CEST MRI enables imaging of distributions of low-concentrated metabolites as well as proteins and peptides and their alterations in diseases. CEST examinations often suffer from low spatial resolution, long acquisition times, and concomitant motion artifacts. This work aims to maximize both resolution and volume coverage with a 3D-EPI snapshot CEST approach at 3T, allowing for fast and robust whole-brain CEST MRI. METHODS: Resolution and temporal SNR of 3D-EPI examinations with nonselective excitation were optimized at a clinical 3T MR scanner in five healthy subjects using a clinical head/neck coil. A CEST presaturation module for low power relayed nuclear Overhauser enhancement and amide proton transfer contrast was applied as an example. The suggested postprocessing included motion correction, dynamic B0 correction, denoising, and B1 correction and was compared to an established 3D-gradient echo-based sequence. RESULTS: CEST examinations were performed at 1.8 mm nominal isotropic resolution in 4.3 s per presaturation offset. In contrast to slab-selective 3D or multislice approaches, the whole brain was covered. Repeated examinations at three different B1 values took 13 minutes for 58 presaturation offsets with temporal SNR around 75. The resulting CEST effects revealed significant gray and white matter contrast and were of similar quality across the whole brain. Coefficient of variation across three healthy subjects was below 9%. CONCLUSION: The suggested protocol enables whole brain coverage at 1.8 mm isotropic resolution and fast acquisition of 4.3 s per presaturation offset. For the fitted CEST amplitudes, high reproducibility was proven, increasing the opportunities of quantitative CEST investigations at 3T significantly.

Structure or Exchange? On the Feasibility of Chemical Exchange Detection with Balanced Steady-State Free Precession in Tissue - An In Vitro Study.

Balanced steady-state free precession imaging has recently been suggested for chemical exchange detection (bSSFPX). The objective of this work is to investigate the contributions of microstructural, chemical shift and chemical exchange effects to the asymmetry of the bSSFP profile at field strengths of 3 T and 9.4 T. To this end, in vitro bSSFPX experiments are performed for a range of repetition times and flip angles in glucose water solutions with different MnCl2 concentrations and tissue homogenates obtained from the brainstem of pig brains. The experimental results are compared to multi-pool Bloch-McConnell simulations. Additionally, the influence of white matter tract geometry is analyzed ex vivo in pig brain hemispheres measured at two different angles with respect to B0 . The detectable bSSFP profile asymmetry in glucose solutions with tissue-like relaxation times and white matter homogenates was consistent with Bloch-McConnell simulations but relatively small. In intact white matter tracts, the asymmetry was dominated by structural effects with a strong dependency on tract orientation relative to B0 . In tracts perpendicular to B0 , the asymmetry was ≈ 3-4 times higher than in the homogenates, thus barely affected by chemical exchange effects. In conclusion, chemical exchange-related bSSFP profile asymmetries are detectable in tissue homogenates, however, the observed asymmetry level is generally low and prone to confounding structural effects rendering in vivo chemical exchange detection with bSSFP challenging in the brain.

Whole-brain B1 -mapping using three-dimensional DREAM.

PURPOSE: Demonstration of a 3D version of the DREAM sequence (3DREAM) for rapid 3D flip angle and B1+ mapping of the human brain. METHODS: A rectangular non-selective STE preparation is followed by a 3D readout with a Cartesian center-out spiral phase encoding order. This enables parallel imaging acceleration in both phase encoding dimensions as well as early capture of the prepared magnetization. RESULTS: B1+ mapping of the whole human head is demonstrated on a 7T system at a nominal resolution of 5 mm with and without parallel imaging acceleration. Artifacts caused by the different signal decay of the FID and STE signal during the long imaging train is suppressed by appropriate filtering of the FID image. Remaining blurring can be controlled by adjusting the echo train length and readout flip angle. CONCLUSIONS: 3DREAM provides a whole-brain flip angle map in a few seconds or individual maps for an 8-channel array in about a minute.

Whole-brain snapshot CEST imaging at 7 T using 3D-EPI.

PURPOSE: The aim of this work is to develop a fast and robust CEST sequence in order to allow the acquisition of a whole-brain imaging volume after a single preparation block (snapshot acquisition). METHODS: A 3D-CEST sequence with an optimized 3D-EPI readout module was developed, which acquires the complete k-space data following a single CEST preparation for 1 saturation offset. Whole-brain mapping of the Z-spectrum with 2 mm isotropic resolution is achieved at 68 saturation frequencies in 5 minutes (4.33 s per offset). We analyzed the B1 distribution in order to optimize B1 correction and to provide accurate CEST quantification across the whole brain. RESULTS: We obtained maps for 3 different CEST contrasts from 4 healthy subjects. Based on our B1 distribution analysis, we conclude that 3 B1 sampling points allow for sufficient compensation of B1 variations across most of the brain. Two brain regions, the cerebellum and the temporal lobes, are difficult to quantify at 7 T due to very low B1 that was achieved in these regions. CONCLUSIONS: The proposed sequence enables robust acquisition of 2 mm isotropic whole-brain CEST maps at 7 Tesla within a total scan time of 16 minutes.

Assessment of frequency drift on CEST MRI and dynamic correction: application to gagCEST at 7 T.

PURPOSE: To investigate the effect of a frequency drift of the static magnetic field on 3D CEST MRI based on glycosaminoglycans (GAGs) of articular cartilage at 7 T and to introduce a retrospective correction method that uses the phase images of the gradient-echo readout. METHODS: Repeated gagCEST and B0 measurements were performed in a glucose model solution and in vivo in the knee joint of 3 healthy volunteers at 7 T. Phase images of the modified 3D rectangular spiral centric-reordered gradient-echo CEST sequence were used to quantify and compensate the apparent frequency drift in repeated gagCEST measurements. RESULTS: The frequency drift of the MRI scanner strongly influences the gagCEST signal in the articular cartilage of the human knee joint. The gagCEST signal in the articular cartilage is changed by 0.18%/Hz while an average drift of 0.7 ± 0.2 Hz/minute was observed. The proposed correction method can be applied retrospectively without the need of additional measurements and provides improved comparability and reproducibility for gagCEST studies. This correction method may also be of interest for other applications of CEST MRI. CONCLUSION: Prospective or retrospective correction of the frequency drift of the MRI scanner is essential for reproducible gagCEST measurements. The proposed retrospective correction method fulfills this requirement without the need of additional measurements.

DeepCEST: 9.4 T Chemical exchange saturation transfer MRI contrast predicted from 3 T data - a proof of concept study.

PURPOSE: To determine the feasibility of employing the prior knowledge of well-separated chemical exchange saturation transfer (CEST) signals in the 9.4 T Z-spectrum to separate overlapping CEST signals acquired at 3 T, using a deep learning approach trained with 3 T and 9.4 T CEST spectral data from brains of the same subjects. METHODS: Highly spectrally resolved Z-spectra from the same volunteer were acquired by 3D-snapshot CEST MRI at 3 T and 9.4 T at low saturation power of B1 = 0.6 µT. The volume-registered 3 T Z-spectra-stack was then used as input data for a three layer deep neural network with the volume-registered 9.4 T fitted parameter stack as target data. RESULTS: An optimized neural net architecture could be found and verified in healthy volunteers. The gray-/white-matter contrast of the different CEST effects was predicted with only small deviations (Pearson R = 0.89). The 9.4 T prediction was less noisy compared to the directly measured CEST maps, although at the cost of slightly lower tissue contrast. Application to an unseen tumor patient measured at 3 T and 9.4 T revealed that tumorous tissue Z-spectra and corresponding hyper-/hypointensities of different CEST effects can also be predicted (Pearson R = 0.84). CONCLUSION: The 9.4 T CEST signals acquired at low saturation power can be accurately estimated from CEST imaging at 3 T using a neural network trained with coregistered 3 T and 9.4 T data of healthy subjects. The deepCEST approach generalizes to Z-spectra of tumor areas and might indicate whether additional ultrahigh-field (UHF) scans will be beneficial.

Human in-vivo brain magnetic resonance current density imaging (MRCDI).

Magnetic resonance current density imaging (MRCDI) and MR electrical impedance tomography (MREIT) are two emerging modalities, which combine weak time-varying currents injected via surface electrodes with magnetic resonance imaging (MRI) to acquire information about the current flow and ohmic conductivity distribution at high spatial resolution. The injected current flow creates a magnetic field in the head, and the component of the induced magnetic field ΔBz,c parallel to the main scanner field causes small shifts in the precession frequency of the magnetization. The measured MRI signal is modulated by these shifts, allowing to determine ΔBz,c for the reconstruction of the current flow and ohmic conductivity. Here, we demonstrate reliable ΔBz,c measurements in-vivo in the human brain based on multi-echo spin echo (MESE) and steady-state free precession free induction decay (SSFP-FID) sequences. In a series of experiments, we optimize their robustness for in-vivo measurements while maintaining a good sensitivity to the current-induced fields. We validate both methods by assessing the linearity of the measured ΔBz,c with respect to the current strength. For the more efficient SSFP-FID measurements, we demonstrate a strong influence of magnetic stray fields on the ΔBz,c images, caused by non-ideal paths of the electrode cables, and validate a correction method. Finally, we perform measurements with two different current injection profiles in five subjects. We demonstrate reliable recordings of ΔBz,c fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔBz,c measurements with simulated ΔBz,c images based on FEM calculations and individualized head models reveals significant linear correlations in all subjects, but only for the stray field-corrected data. As final step, we reconstruct current density distributions from the measured and simulated ΔBz,c data. Reconstructions from non-corrected ΔBz,c measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔBz,c measurements and from simulated ΔBz,c images reveals an average coefficient of determination R2 of 71%. In addition, it shows that the simulations underestimated the current strength on average by 24%. Our results open up the possibility of using MRI to systematically validate and optimize numerical field simulations that play an important role in several neuroscience applications, such as transcranial brain stimulation, and electro- and magnetoencephalography.

Chemical exchange saturation transfer MRI contrast in the human brain at 9.4 T.

PURPOSE: The high chemical shift separation at 9.4 T allows for selective saturation of proton pools in exchange with water protons. For the first time, highly selective and comprehensive chemical exchange saturation transfer (CEST) experiments were performed in the human brain at 9.4 T. This work provides insight into CEST signals in the human brain in comparison with existing animal studies, as well as with CEST effects in vivo at lower field strengths. METHODS: A novel snapshot-CEST method for human brain scans at 9.4 T was optimized and employed for highly-spectrally-resolved (95 offsets) CEST measurements in healthy subjects and one brain tumor patient. Reproducibility and stability between scans was verified in grey and white matter after B0, B1, and motion correction of the acquired 3D CEST volumes. Two-step Lorentzian fitting was used to further improve separation of spectrally discernible signals to create known and novel CEST contrast maps at 9.4 T. RESULTS: At a saturation power of B1 = 0.5 µT most selective CEST effects could be obtained in the human brain with high inter-scan reproducibility. While contrast behavior of previously measured signals at lower field, namely amide-, guanidyl- and NOE-CEST effects, could be reproduced, novel signals at 2.7 ppm, and -1.6 ppm could be verified in healthy subjects and in a brain tumor patient for the first time. CONCLUSION: High spectral resolution chemical exchange saturation transfer at 9.4 T allows deeper insights into the Z-spectrum structure of the human brain, and provides many different contrasts showing different correlations in healthy tissue and in tumor-affected areas of the brain, generating hypotheses for future investigations of in-vivo-CEST at UHF.