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.

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.

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.

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 2023;93:686-701.

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 B₀ 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.

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 (R₁ , 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 R₁ , 2.8%-10.6% for R 2 * $$ {R}_2^{\ast } $$ and MTsat, which was comparable with FLASH (0.6%-4.9% for PD and R₁ , 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, R₁ , 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.

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.

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.

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.

Temporal SNR optimization through RF coil combination in fMRI: The more, the better?

For functional MRI with a multi-channel receiver RF coil, images are often reconstructed channel by channel, resulting into multiple images per time frame. The final image to analyze usually is the result of the covariance Sum-of-Squares (covSoS) combination across these channels. Although this reconstruction is quasi-optimal in SNR, it is not necessarily the case in terms of temporal SNR (tSNR) of the time series, which is yet a more relevant metric for fMRI data quality. In this work, we investigated tSNR optimality through voxel-wise RF coil combination and its effects on BOLD sensitivity. An analytical solution for an optimal RF coil combination is described, which is somewhat tied to the extended Krueger-Glover model involving both thermal and physiological noise covariance matrices. Compared experimentally to covSOS on four volunteers at 7T, the method yielded great improvement of tSNR but, surprisingly, did not result into higher BOLD sensitivity. Solutions to improve the method such as for example the t-score for the mean recently proposed are also explored, but result into similar observations once the statistics are corrected properly. Overall, the work shows that data-driven RF coil combinations based on tSNR considerations alone should be avoided unless additional and unbiased assumptions can be made.

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.

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 MnCl₂ 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 B₀ . 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 B₀ . In tracts perpendicular to B₀ , 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.

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 B₀ 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.

Snapshot-CEST: Optimizing spiral-centric-reordered gradient echo acquisition for fast and robust 3D CEST MRI at 9.4 T

Gradient echo (GRE)-based acquisition provides a robust readout method for chemical exchange saturation transfer (CEST) at ultrahigh field (UHF). To develop a snapshot-CEST approach, the transient GRE signal and point spread function were investigated in detail, leading to optimized measurement parameters and reordering schemes for fast and robust volumetric CEST imaging. Simulation of the transient GRE signal was used to determine the optimal sequence parameters and the maximum feasible number of k-space lines. Point spread function analysis provided an insight into the induced k-space filtering and the performance of different rectangular reordering schemes in terms of blurring, signal-to-noise ratio (SNR) and relaxation dependence. Simulation results were confirmed in magnetic resonance imaging (MRI) measurements of healthy subjects. Minimal repetition time (TR) is beneficial for snapshot-GRE readout. At 9.4 T, for TR = 4 ms and optimal flip angle close to the Ernst angle, a maximum of 562 k-space lines can be acquired after a single presaturation, providing decent SNR with high image quality. For spiral-centric reordered k-space acquisition, the image quality can be further improved using a rectangular spiral reordering scheme adjusted to the field of view. Application of the derived snapshot-CEST sequence for fast imaging acquisition in the human brain at 9.4 T shows excellent image quality in amide and nuclear Overhauser enhancement (NOE), and enables guanidyl CEST detection. The proposed snapshot-CEST establishes a fast and robust volumetric CEST approach ready for the imaging of known and novel exchange-weighted contrasts at UHF.

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 ΔB_z,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 ΔB_z,c for the reconstruction of the current flow and ohmic conductivity. Here, we demonstrate reliable ΔB_z,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 ΔB_z,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 ΔB_z,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 ΔB_z,c fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔB_z,c measurements with simulated ΔB_z,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 ΔB_z,c data. Reconstructions from non-corrected ΔB_z,c measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔB_z,c measurements and from simulated ΔB_z,c images reveals an average coefficient of determination R² 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.

SAR and scan-time optimized 3D whole-brain double inversion recovery imaging at 7T

PURPOSE

The aim of this project was to implement an ultra-high field (UHF) optimized double inversion recovery (DIR) sequence for gray matter (GM) imaging, enabling whole brain coverage in short acquisition times ( ≈5 min, image resolution 1 mm³ ).

METHODS

A 3D variable flip angle DIR turbo spin echo (TSE) sequence was optimized for UHF application. We implemented an improved, fast, and specific absorption rate (SAR) efficient TSE imaging module, utilizing improved reordering. The DIR preparation was tailored to UHF application. Additionally, fat artifacts were minimized by employing water excitation instead of fat saturation.

RESULTS

GM images, covering the whole brain, were acquired in 7 min scan time at 1 mm isotropic resolution. SAR issues were overcome by using a dedicated flip angle calculation considering SAR and SNR efficiency. Furthermore, UHF related artifacts were minimized.

CONCLUSION

The suggested sequence is suitable to generate GM images with whole-brain coverage at UHF. Due to the short total acquisition times and overall robustness, this approach can potentially enable DIR application in a routine setting and enhance lesion detection in neurological diseases. Magn Reson Med 79:2620-2628, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Multiline balanced SSFP for rapid functional imaging at ultrahigh field

PURPOSE

The goal of this study is to develop and evaluate a multiline balanced steady-state free-precession (bSSFP) sequence for passband functional MRI at ultrahigh field.

METHODS

Passband bSSFP functional MRI experiments using a visual task were performed on a 9.4 T system with echo trains ranging from one up to seven echoes. We analyze the acquisition efficiency, temporal and thermal signal-to-noise ratio, as well as the observed blood oxygen-level-dependent (BOLD) signal changes.

RESULTS

With increasing repetition time and echo train length, the BOLD-related signal change as well as the thermal and temporal noise were improved. Activation patterns and signal changes were stable and reproducible across subjects.

CONCLUSIONS

We propose a multiline bSSFP for functional BOLD imaging that approaches the speed of echo-planar imaging and that shows an increased BOLD sensitivity compared with single-line bSSFP. Magn Reson Med 79:994-1000, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Depth-dependence of visual signals in the human superior colliculus at 9.4 T

The superior colliculus (SC) is a layered structure located in the midbrain. We exploited the improved spatial resolution and BOLD signal strength available at 9.4 T to investigate the depth profile of visual BOLD responses in the human SC based on distortion-corrected EPI data with a 1 mm isotropic resolution. We used high resolution (350 µm in-plane) anatomical images to determine regions-of-interest of the SC and applied a semi-automated method to segment it into superficial, intermediate, and deep zones. A greater than linear increase in sensitivity of the functional signal at 9.4 T allowed us to detect a statistically significant depth pattern in a group analysis with a 20 min stimulation paradigm. Descriptive data showed consistent depth profiles also in single individuals. The highest signals were localized to the superficial layers of the right and left SC during contralateral stimulation, which was in good agreement with its functional architecture known from non-human primates. This study thus demonstrates the potential of 9.4 T MRI for functional neuroimaging even in deeply located, particularly challenging brain structures such as the SC. Hum Brain Mapp 38:574-587, 2017. © 2016 Wiley Periodicals, Inc.

Desynchronization of Cartesian k-space sampling and periodic motion for improved retrospectively self-gated 3D lung MRI using quasi-random numbers

PURPOSE

To demonstrate that desynchronization between Cartesian k-space sampling and periodic motion in free-breathing lung MRI improves the robustness and efficiency of retrospective respiratory self-gating.

METHODS

Desynchronization was accomplished by reordering the phase (k_y ) and partition (k_z ) encoding of a three-dimensional FLASH sequence according to two-dimensional, quasi-random (QR) numbers. For retrospective respiratory self-gating, the k-space center signal (DC signal) was acquired separately after each encoded k-space line. QR sampling results in a uniform distribution of k-space lines after gating. Missing lines resulting from the gating process were reconstructed using iterative GRAPPA. Volunteer measurements were performed to compare quasi-random with conventional sampling. Patient measurements were performed to demonstrate the feasibility of QR sampling in a clinical setting.

RESULTS

The uniformly sampled k-space after retrospective gating allows for a more stable iterative GRAPPA reconstruction and improved ghost artifact reduction compared with conventional sampling. It is shown that this stability can either be used to reduce the total scan time or to reconstruct artifact-free data sets in different respiratory phases, both resulting in an improved efficiency of retrospective respiratory self-gating.

CONCLUSION

QR sampling leads to desynchronization between repeated data acquisition and periodic respiratory motion. This results in an improved motion artifact reduction in shorter scan time. Magn Reson Med 77:787-793, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Dynamically phase-cycled radial balanced SSFP imaging for efficient banding removal

PURPOSE

Balanced steady-state free precession (bSSFP) imaging suffers from banding artifacts due to its inherent sensitivity to inhomogeneities in the main magnetic field. These artifacts can be removed by the acquisition of multiple images at different frequency offsets. However, conventional phase-cycling is hindered by a long scan time. The purpose of this work is to present a novel approach for efficient banding removal in bSSFP imaging.

THEORY AND METHODS

To this end, the phase-cycle during a single-shot radial acquisition of an image was dynamically changed. Thus, each projection is acquired with a different frequency offset. Using conventional radial gridding, an artifact-free image can be reconstructed out of this dataset.

RESULTS

The approach is validated at clinical field strength [3.0 Tesla (T)] as well as at ultrahigh field (9.4T). Robust elimination of banding artifacts was obtained for different imaging regions, including brain imaging at ultrahigh field with an in-plane resolution of 0.25 × 0.25 mm(2). Besides banding artifact-free imaging, the applicability of the proposed technique for fat-water separation is demonstrated.

CONCLUSION

Dynamically phase-cycled radial bSSFP has the potential for banding-free bSSFP imaging in a short scan time, in the presence of severe field inhomogeneities and at high resolution.

Three-dimensional through-time radial GRAPPA for renal MR angiography

PURPOSE

To achieve high temporal and spatial resolution for contrast-enhanced time-resolved MR angiography exams (trMRAs), fast imaging techniques such as non-Cartesian parallel imaging must be used. In this study, the three-dimensional (3D) through-time radial generalized autocalibrating partially parallel acquisition (GRAPPA) method is used to reconstruct highly accelerated stack-of-stars data for time-resolved renal MRAs.

MATERIALS AND METHODS

Through-time radial GRAPPA has been recently introduced as a method for non-Cartesian GRAPPA weight calibration, and a similar concept can also be used in 3D acquisitions. By combining different sources of calibration information, acquisition time can be reduced. Here, different GRAPPA weight calibration schemes are explored in simulation, and the results are applied to reconstruct undersampled stack-of-stars data.

RESULTS

Simulations demonstrate that an accurate and efficient approach to 3D calibration is to combine a small number of central partitions with as many temporal repetitions as exam time permits. These findings were used to reconstruct renal trMRA data with an in-plane acceleration factor as high as 12.6 with respect to the Nyquist sampling criterion, where the lowest root mean squared error value of 16.4% was achieved when using a calibration scheme with 8 partitions, 16 repetitions, and a 4 projection × 8 read point segment size.

CONCLUSION

3D through-time radial GRAPPA can be used to successfully reconstruct highly accelerated non-Cartesian data. By using in-plane radial undersampling, a trMRA can be acquired with a temporal footprint less than 4s/frame with a spatial resolution of approximately 1.5 mm × 1.5 mm × 3 mm.

Dynamic MRI of the TMJ under physical load

OBJECTIVES

The objective of this study was to examine the kinematics of structures of the temporomandibular joint (TMJ) under physiological load while masticating.

METHODS

Radial MRI was chosen as a fast imaging method to dynamically capture the motions of the joint's anatomy. The technique included a golden ratio-based increment angle and a sliding window reconstruction. The measurements were performed on 22 subjects with and without deformation/displacement of the intra-articular disc while they were biting on a cooled caramel toffee.

RESULTS

The reconstructed dynamic images provided sufficient information about the size and localization of the disc as well as the change of the intra-articular distance with and without loading.

CONCLUSIONS

The feasibility of the golden ratio-based radial MRI technique to dynamically capture the anatomy of the TMJ under physical load was demonstrated in this initial study.

DC-gated high resolution three-dimensional lung imaging during free-breathing

PURPOSE

To use the acquisition of the k-space center signal (DC signal) implemented into a Cartesian three-dimensional (3D) FLASH sequence for retrospective respiratory self-gating and, thus, for the examination of the whole human lung in high spatial resolution during free breathing.

MATERIALS AND METHODS

Volunteer as well as patient measurements were performed under free breathing conditions. The DC signal is acquired after the actual image data acquisition within each excitation of a 3D FLASH sequence. The DC signal is then used to track respiratory motion for retrospective respiratory gating.

RESULTS

It is shown that the acquisition of the DC signal after the imaging module can be used in a 3D FLASH sequence to extract respiratory motion information for retrospective respiratory self-gating and allows for shorter echo times (TE) and therefore increased lung parenchyma SNR.

CONCLUSION

The acquisition of the DC signal after image signal acquisition allows successful retrospective gating, enabling the reconstruction of high resolution images of the whole human lung under free breathing conditions.

Improved temporal resolution in cardiac imaging using through-time spiral GRAPPA

Previous work has shown that the use of radial GRAPPA for the reconstruction of undersampled real-time free-breathing cardiac data allows for frame rates of up to 30 images/s. It is well known that the spiral trajectory offers a higher scan efficiency compared to radial trajectories. For this reason, we have developed a novel through-time spiral GRAPPA method and demonstrate its application to real-time cardiac imaging. By moving from the radial trajectory to the spiral trajectory, the temporal resolution can be further improved at lower acceleration factors compared to radial GRAPPA. In addition, the image quality is improved compared to those generated using the radial trajectory due to the lower acceleration factor. Here, we show that 2D frame rates of up to 56 images/s can be achieved using this parallel imaging method with the spiral trajectory.

Improved radial GRAPPA calibration for real-time free-breathing cardiac imaging

To generate real-time, nongated, free-breathing cardiac images, the undersampled radial trajectory combined with parallel imaging in the form of radial GRAPPA has shown promise. However, this method starts to fail at high undersampling factors due to the assumptions that must be made for the purposes of calibrating the GRAPPA weight sets. In this manuscript, a novel through-time radial GRAPPA calibration scheme is proposed which greatly improves image quality for the high acceleration factors required for real-time cardiac imaging. This through-time calibration method offers better image quality than standard radial GRAPPA, but it requires many additional calibration frames to be acquired. By combining the through-time calibration method proposed here with the standard through-k-space radial GRAPPA calibration method, images with high acceleration factors can be reconstructed using few calibration frames. Both the through-time and the hybrid through-time/through-k-space methods are investigated to determine the most advantageous calibration parameters for an R = 6 in vivo short-axis cardiac image. Once the calibration parameters have been established, they are then used to reconstruct several in vivo real-time, free-breathing cardiac datasets with temporal resolutions better than 45 msec, including one with a temporal resolution of 35 msec and an in-plane resolution of 1.56 mm(2) .

Using the GRAPPA operator and the generalized sampling theorem to reconstruct undersampled non-Cartesian data

As expected from the generalized sampling theorem of Papoulis, the use of a bunched sampling acquisition scheme in conjunction with a conjugate gradient (CG) reconstruction algorithm can decrease scan time by reducing the number of phase-encoding lines needed to generate an unaliased image at a given resolution. However, the acquisition of such bunched data requires both modified pulse sequences and high gradient performance. A novel method of generating the "bunched" data using self-calibrating GRAPPA operator gridding (GROG), a parallel imaging method that shifts data points by small distances in k-space (with Deltak usually less than 1.0, depending on the receiver coil) using the GRAPPA operator, is presented here. With the CG reconstruction method, these additional "bunched" points can then be used to reconstruct an image with reduced artifacts from undersampled data. This method is referred to as GROG-facilitated bunched phase encoding (BPE), or GROG-BPE. To better understand how the patterns of bunched points, maximal blip size, and number of bunched points affect the reconstruction quality, a number of simulations were performed using the GROG-BPE approach. Finally, to demonstrate that this method can be combined with a variety of trajectories, examples of images with reduced artifacts reconstructed from undersampled in vivo radial, spiral, and rosette data are shown.

Solution structure of the Legionella pneumophila Mip-rapamycin complex

Background

Legionella pneumphila is the causative agent of Legionnaires' disease. A major virulence factor of the pathogen is the homodimeric surface protein Mip. It shows peptidyl-prolyl cis/trans isomerase activty and is a receptor of FK506 and rapamycin, which both inhibit its enzymatic function. Insight into the binding process may be used for the design of novel Mip inhibitors as potential drugs against Legionnaires' disease.

Results

We have solved the solution structure of free Mip^77–213 and the Mip^77–213-rapamycin complex by NMR spectroscopy. Mip^77–213 showed the typical FKBP-fold and only minor rearrangements upon binding of rapamycin. Apart from the configuration of a flexible hairpin loop, which is partly stabilized upon binding, the solution structure confirms the crystal structure. Comparisons to the structures of free FKBP12 and the FKBP12-rapamycin complex suggested an identical binding mode for both proteins.

Conclusion

The structural similarity of the Mip-rapamycin and FKBP12-rapamycin complexes suggests that FKBP12 ligands may be promising starting points for the design of novel Mip inhibitors. The search for a novel drug against Legionnaires' disease may therefore benefit from the large variety of known FKBP12 inhibitors.

Domain motions of the Mip protein from Legionella pneumophila

The homodimeric 45.6 kDa (total mass) Mip protein, a virulence factor from Legionella pneumophila, was investigated with solution NMR spectroscopy and molecular dynamics (MD) simulations. Two Mip monomers are dimerized via an N-terminal helix bundle that is connected via a long alpha-helix to a C-terminal FKBP domain in each subunit. More than 85% of the amino acids were identified in triple-resonance NMR spectra. (15)N relaxation analysis showed a bimodal distribution of R(1)/R(2) values, with the lower ratio in the N-terminal domain. Relaxation dispersion measurements confirmed that these reduced ratios did not originate from conformational exchange. Thus, two different correlation times (tau(c)) can be deduced, reflecting partly uncoupled motions of both domains. Relaxation data of a Mip(77)(-)(213) monomer mutant were similar to those observed in the dimer, corroborating that the FKBP domain, including part of the connecting helix, behaves as one dynamic entity. MD simulations (18 ns) of the Mip dimer also yielded two different correlation times for the two domains and thus confirm the independence of the domain motions. Principal component analysis of the dihedral space covariance matrix calculated from the MD trajectory suggests a flexible region in the long connecting helix that acts as a hinge between the two domains. Such motion provides a possible explanation of how Mip can bind to complex molecular components of the extracellular matrix and mediate alveolar damage and bacterial spread in the lung.