The restricted SAR protocol: A method to assess MRI coil prototypes in an unconditionally safe manner.

PURPOSE: Testing an RF coil prototype on subjects involves laborious verifications to ensure its safety. In particular, it requires preliminary electromagnetic simulations and their validations on phantoms to accurately predict the specific absorption rate (SAR). For coil design validation with a simpler safety procedure, the restricted SAR (rS) mode is proposed, enabling representative first experiments in vivo. The goal of the developed approach is to accelerate the transition of a custom coil system from prototype to clinical use. METHODS: The restricted specific absorption rate (SAR) (rS) mode imposes a radical limitation on the transmitted RF power based on a worst-case scenario of local RF power absorption. The limitations used are independent of the SAR spatial distribution, making this approach unconditionally safe. The developed rS protocol contains the sequences required for coil evaluation and satisfies the imposed rS conditions. It provides a quantitative characterization of the coil transmission and reception profiles and a qualitative evaluation of the anatomical images. Protocol validation was performed on commercial and pre-industrial prototype coils on a small cohort of healthy volunteers. RESULTS: The proposed rS protocol enables coil evaluation within an acquisition time compatible with common clinical protocol duration. The total time of all evaluation steps does not exceed 17 min. At the same time, the global SAR remains 100 times less than the International Electrotechnical Commission safety limit for played sequences. CONCLUSION: The rS protocol allows characterizing and comparing coil prototypes on volunteers without extensive electromagnetic calculations and phantom validations in an unconditionally safe way.

Pushing MP2RAGE boundaries: Ultimate time-efficient parameterization combined with exhaustive T1 synthetic contrasts.

PURPOSE: MP2RAGE parameter optimization is redefined to allow more time-efficient MR acquisitions, whereas the T1 -based synthetic imaging framework is used to obtain on-demand T1 -weighted contrasts. Our aim was to validate this concept on healthy volunteers and patients with multiple sclerosis, using plug-and-play parallel-transmission brain imaging at 7 T. METHODS: A "time-efficient" MP2RAGE sequence was designed with optimized parameters including TI and TR set as small as possible. Extended phase graph formalism was used to set flip-angle values to maximize the gray-to-white-matter contrast-to-noise ratio (CNR). Several synthetic contrasts (UNI, EDGE, FGATIR, FLAWSMIN , FLAWSHCO ) were generated online based on the acquired T1 maps. Experimental validation was performed on 4 healthy volunteers at various spatial resolutions. Clinical applicability was evaluated on 6 patients with multiple sclerosis, scanned with both time-efficient and conventional MP2RAGE parameterizations. RESULTS: The proposed time-efficient MP2RAGE protocols reduced acquisition time by 40%, 30%, and 19% for brain imaging at (1 mm)3 , (0.80 mm)3 and (0.65 mm)3 , respectively, when compared with conventional parameterizations. They also provided all synthetic contrasts and comparable contrast-to-noise ratio on UNI images. The flexibility in parameter selection allowed us to obtain a whole-brain (0.45 mm)3 acquisition in 19 min 56 s. On patients with multiple sclerosis, a (0.67 mm)3 time-efficient acquisition enhanced cortical lesion visualization compared with a conventional (0.80 mm)3 protocol, while decreasing the scan time by 15%. CONCLUSION: The proposed optimization, associated with T1 -based synthetic contrasts, enabled substantial decrease of the acquisition time or higher spatial resolution scans for a given time budget, while generating all typical brain contrasts derived from MP2RAGE.

Impact of B 0 $$ {\mathrm{B}}_0 $$ field imperfections correction on BOLD sensitivity in 3D-SPARKLING fMRI data.

PURPOSE: Static and dynamic B 0 $$ {\mathrm{B}}_0 $$ field imperfections are detrimental to functional MRI (fMRI) applications, especially at ultra-high magnetic fields (UHF). In this work, a field camera is used to assess the benefits of retrospectively correcting B 0 $$ {\mathrm{B}}_0 $$ field perturbations on Blood Oxygen Level Dependent (BOLD) sensitivity in non-Cartesian three-dimensional (3D)-SPARKLING fMRI acquisitions. METHODS: fMRI data were acquired at 1 mm 3 $$ {}^3 $$ and for a 2.4s-TR while concurrently monitoring in real-time field perturbations using a Skope Clip-on field camera in a novel experimental setting involving a shorter TR than the required minimal TR of the field probes. Measurements of the dynamic field deviations were used along with a static Δ B 0 $$ \Delta {\mathrm{B}}_0 $$ map to retrospectively correct static and dynamic field imperfections, respectively. In order to evaluate the impact of such a correction on fMRI volumes, a comparative study was conducted on healthy volunteers. RESULTS: Correction of B 0 $$ {\mathrm{B}}_0 $$ deviations improved image quality and yielded between 20% and 30% increase in median temporal signal-to-noise ratio (tSNR).Using fMRI data collected during a retinotopic mapping experiment, we demonstrated a significant increase in sensitivity to the BOLD contrast and improved accuracy of the BOLD phase maps: 44% (resp., 159%) more activated voxels were retrieved when using a significance control level based on a p-value of 0.001 without correcting for multiple comparisons (resp., 0.05 with a false discovery rate correction). CONCLUSION: 3D-SPARKLING fMRI hugely benefits from static and dynamic B 0 $$ {\mathrm{B}}_0 $$ imperfections correction. However, the proposed experimental protocol is flexible enough to be deployed on a large spectrum of encoding schemes, including arbitrary non-Cartesian readouts.

Efficiently building receive arrays with electromagnetic simulations and additive manufacturing: A two-layer, 32-channel prototype for 7T brain MRI.

PURPOSE: We propose a comprehensive workflow to design and build fully customized dense receive arrays for MRI, providing prediction of SNR and g-factor. Combined with additive manufacturing, this method allows an efficient implementation for any arbitrary loop configuration. To demonstrate the methodology, an innovative two-layer, 32-channel receive array is proposed. METHODS: The design workflow is based on numerical simulations using a commercial 3D electromagnetic software associated with circuit model co-simulations to provide the most accurate results in an efficient time. A model to compute the noise covariance matrix from circuit model scattering parameters is proposed. A 32-channel receive array at 7 T is simulated and fabricated with a two-layer design made of non-geometrically decoupled loops. Decoupling between loops is achieved using home-built direct high-impedance preamplifiers. The loops are 3D-printed with a new additive manufacturing technique to speed up integration while preserving the detailed geometry as simulated. The SNR and parallel-imaging performances of the proposed design are compared with a commercial coil, and in vivo images are acquired. RESULTS: The comparison of SNR and g-factors showed a good agreement between simulations and measurements. Experimental values are comparable with the ones measured on the commercial coil. Preliminary in vivo images also ensured the absence of any unexpected artifacts. CONCLUSION: A new design and performance analysis workflow is proposed and tested with a non-conventional 32-channel prototype at 7 T. Additive manufacturing of dense arrays of loops for brain imaging at ultrahigh field is validated for clinical use.

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.

Efficient optimization of chemical exchange saturation transfer MRI at 7 T using universal pulses and virtual observation points.

PURPOSE: To optimize the homogeneity of the presaturation module in a Chemical Exchange Saturation Transfer (CEST) acquisition at 7 T using parallel transmission (pTx). THEORY AND METHODS: An optimized pTx-CEST presaturation scheme based on precomputed universal pulses was designed. The optimization was performed by minimizing the L2-norm between the effective B 1 , RMS + $$ {B}_{1,\mathrm{RMS}}^{+} $$ and a given target while imposing energy constraints under virtual observation points (VOPs) supervision. The proposed method was evaluated through simulations and experimentally, both in vitro, on a realistic human head phantom, and in vivo, on healthy volunteers. The results were compared with circular polarization (CP) presaturation and other pTx approaches previously proposed. All experiments were conducted on a 7 T MRI scanner using a commercial 8Tx/32Rx head coil. RESULTS: The simulations show that the proposed pTx strategy boosted with VOPs is superior to the CP mode and existent pTx approaches. While the best results are obtained with subject specific pulses, the gain provided by the use of VOPs renders the universal pulses superior to tailored pulses optimized under vendor provided Specific Absorption Rate (SAR) management. In the phantom, the glucose MTR asym $$ {\mathrm{MTR}}_{\mathrm{asym}} $$ map was significantly more homogeneous than with CP (root mean square error [RMSE] 17% vs. 30%). The efficiency of the method for in vivo hydroxyl, glutamate and rNOE weighted CEST acquisitions was also demonstrated. CONCLUSION: The use of a pTx presaturation scheme based on universal pulses optimized under VOP SAR management is significantly benefiting CEST imaging at high magnetic field.

Improving spreading projection algorithm for rapid k-space sampling trajectories through minimized off-resonance effects and gridding of low frequencies.

PURPOSE: Non-Cartesian MRI with long arbitrary readout directions are susceptible to off-resonance artifacts due to patient induced B 0 $$ {B}_0 $$ inhomogeneities. This results in degraded image quality with strong signal losses and blurring. Current solutions to address this issue involve correcting the off-resonance artifacts during image reconstruction or reducing inhomogeneities through improved shimming. THEORY: The recently developed SPARKLING algorithm is extended to drastically diminish off-resonance artifacts by generating temporally smooth k-space sampling patterns. For doing so, the cost function which is optimized in SPARKLING is modified using a temporal weighting factor. Additionally, oversampling of the center of k-space beyond the Nyquist criteria is prevented through the use of gridded sampling in the region, enforced with affine constraints. METHODS: Prospective k-space data was acquired at 3 T on new trajectories, and we show robustness to B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities through in silico experiments by adding Δ B 0 $$ \Delta {\mathrm{B}}_0 $$  through artificial degradation of system B 0 $$ {\mathrm{B}}_0 $$ shimming. Later on, in vivo experiments were carried out to optimize parameters of the new improvements and benchmark the gain in performance. RESULTS: The improved trajectories allowed for the recovery of signal dropouts observed on original SPARKLING acquisitions at larger B 0 $$ {\mathrm{B}}_0 $$ field inhomogeneities. Furthermore, imposing gridded sampling at the center of k-space provided improved reconstructed image quality with limited artifacts. CONCLUSION: These advancements allowed us for nearly 4 . 62 × $$ 4.62\times $$ shorter scan time compared to GRAPPA-p4x1, allowing us to reach 600 µm isotropic resolution in 3D T 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ -w imaging in just 3.3 min at 3 T with negligible degradation in image quality.

Deep learning-assisted model-based off-resonance correction for non-Cartesian SWI.

PURPOSE: Patient-induced inhomogeneities in the static magnetic field cause distortions and blurring (off-resonance artifacts) during acquisitions with long readouts such as in SWI. Conventional versatile correction methods based on extended Fourier models are too slow for clinical practice in computationally demanding cases such as 3D high-resolution non-Cartesian multi-coil acquisitions. THEORY: Most reconstruction methods can be accelerated when performing off-resonance correction by reducing the number of iterations, compressed coils, and correction components. Recent state-of-the-art unrolled deep learning architectures could help but are generally not adapted to corrupted measurements as they rely on the standard Fourier operator in the data consistency term. The combination of correction models and neural networks is therefore necessary to reduce reconstruction times. METHODS: Hybrid pipelines using UNets were trained stack-by-stack over 99 SWI 3D SPARKLING 20-fold accelerated acquisitions at 0.6 mm isotropic resolution using different off-resonance correction methods. Target images were obtained using slow model-based corrections based on self-estimated Δ B 0 $$ \Delta {B}_0 $$ field maps. The proposed strategies, tested over 11 volumes, are compared to model-only and network-only pipelines. RESULTS: The proposed hybrid pipelines achieved scores competing with two to three times slower baseline methods, and neural networks were observed to contribute both as pre-conditioner and through inter-iteration memory by allowing more degrees of freedom over the model design. CONCLUSION: A combination of model-based and network-based off-resonance correction was proposed to significantly accelerate conventional methods. Different promising synergies were observed between acceleration factors (iterations, coils, correction) and model/network that could be expanded in the future.

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.

Iterative static field map estimation for off-resonance correction in non-Cartesian susceptibility weighted imaging.

PURPOSE: Patient-induced inhomogeneities in the magnetic field cause distortions and blurring during acquisitions with long readouts such as in susceptibility-weighted imaging (SWI). Most correction methods require collecting an additional ΔB0$$ \Delta {\mathrm{B}}_0 $$ field map to remove these artifacts. THEORY: The static ΔB0$$ \Delta {\mathrm{B}}_0 $$ field map can be approximated with an acceptable error directly from a single echo acquisition in SWI. The main component of the observed phase is linearly related to ΔB0$$ \Delta {\mathrm{B}}_0 $$ and the echo time (TE), and the relative impact of non- ΔB0$$ \Delta {\mathrm{B}}_0 $$ terms becomes insignificant with TE$$ \mathrm{TE} $$ >20 ms at 3 T for a well-tuned system. METHODS: The main step is to combine and unfold the multi-channel phase maps wrapped many times, and several competing algorithms are compared for this purpose. Four in vivo brain data sets collected using the recently proposed 3D spreading projection algorithm for rapid k-space sampling (SPARKLING) readouts are used to assess the proposed method. RESULTS: The estimated 3D field maps generated with a 0.6 mm isotropic spatial resolution provide overall similar off-resonance corrections compared to reference corrections based on an external ΔB0$$ \Delta {\mathrm{B}}_0 $$ acquisitions, and even improved for 2 of 4 individuals. Although a small estimation error is expected, no aftermath was observed in the proposed corrections, whereas degradations were observed in the references. CONCLUSION: A static ΔB0$$ \Delta {\mathrm{B}}_0 $$ field map estimation method was proposed to take advantage of acquisitions with long echo times, and outperformed the reference technique based on an external field map. The difference can be attributed to an inherent robustness to mismatches between volumes and external ΔB0$$ \Delta {\mathrm{B}}_0 $$ maps, and diverse other sources investigated.

Measuring radiofrequency field-induced temperature variations in brain MRI exams with motion compensated MR thermometry and field monitoring.

PURPOSE: An MR thermometry (MRT) method with motion and field fluctuation compensation is proposed to measure non-invasively sub-degree brain temperature variations occurring through radiofrequency (RF) power deposition during MR exams. METHODS: MRT at 7T with a multi-slice echo planar imaging (EPI) sequence and concurrent field monitoring was first tested in vitro to assess accuracy in the presence of external field perturbations, an optical probe being used for ground truth. In vivo, this strategy was complemented by a motion compensation scheme based on a dictionary pre-scan, as reported in some previous work, and was adapted to the human brain. Precision reached with this scheme was assessed on eight volunteers with a 5 minute-long low specific absorption rate (SAR) scan. Finally, temperature rise in the brain was measured twice on the same volunteers and with the same strategy, this time by employing a 20-minutes scan at the maximum SAR delivered with a commercial volume head coil. RESULTS: In vitro, the root mean square (RMS) error between optical probe and MRT measurements was 0.02°C with field sensor correction. In vivo, the low SAR scan returned a precision in temperature change measurement with field monitoring and motion compensation of 0.05°C. The 20-minutes maximum SAR scan returned a temperature rise throughout the inner-brain in the range of 0-0.2°C. Brain periphery remained too sensitive with respect to motion to lead to equally conclusive results. CONCLUSION: Sub-degree temperature rise in the inner human brain was characterized experimentally throughout RF exposure. Potential applications include improvement of human thermal models and revision of safety margins.

Magnetic field strength dependent SNR gain at the center of a spherical phantom and up to 11.7T.

PURPOSE: The SNR at the center of a spherical phantom of known electrical properties was measured in quasi-identical experimental conditions as a function of magnetic field strength between 3 T and 11.7 T. METHODS: The SNR was measured at the center of a spherical water saline phantom with a gradient-recalled echo sequence. Measurements were performed at NeuroSpin at 3, 7, and 11.7 T. The phantom was then shipped to Maastricht University and then to the University of Minnesota for additional data points at 7, 9.4, and 10.5 T. Experiments were carried out with the exact same type of birdcage volume coil (except at 3 T, where a similar coil was used) to attempt at isolating the evolution of SNR with field strength alone. Phantom electrical properties were characterized over the corresponding frequency range. RESULTS: Electrical properties were found to barely vary over the frequency range. Removing the influence of the flip-angle excitation inhomogeneity was crucial, as expected. After such correction, measurements revealed a gain of SNR growing as B0 1.94 ± 0.16 compared with B0 2.13 according to ultimate intrinsic SNR theory. CONCLUSIONS: By using quasi-identical experimental setups (RF volume coil, phantom, electrical properties, and protocol), this work reports experimental data between 3 T and 11.7 T, enabling the comparison with SNR theories in which conductivity and permittivity can be assumed to be constant with respect to field strength. According to ultimate SNR theory, these results can be reasonably extrapolated to the performance of receive arrays with greater than about 32 elements for central SNR in the same spherical phantom.

Open birdcage coil for head imaging at 7T.

PURPOSE: To theoretically describe, design, and test the new geometry of the birdcage coil for 7 Tesla anatomical brain imaging, which includes a large window on top, without deliberately jeopardizing its homogeneity and efficiency. This opencage will not only improve patient comfort but also enable the volunteer to follow functional MRI stimuli. This design could also facilitate the tracking of patient compliance and enable better correction of the movement. METHODS: Via the transfer matrix approach, a birdcage-like coil with a nonperiodic distribution of rungs is constructed with optimized currents in the coil rungs. Subsequently, the coil is adjusted in full-wave simulations. Then, the coil is assembled, fine-tuned, and matched on the bench. Finally, these results are confirmed experimentally on a phantom and in vivo. RESULTS: Indeed, the computed isolation of -14.9 dB between the feeding ports of the coil and the symmetry of the circular polarized mode pattern transmit RF magnetic field ( B1+ ) showed that the coil was properly optimized. An experimental assessment of the developed coil showed competitive transmit efficiency and coverage compared with the conventional birdcage coil of similar size. CONCLUSION: The proposed opencage coil can be designed and work without a dramatic drop of performance in terms of the B1+ field homogeneity, transmit efficiency ( B1+ / Pref ), peak local specific absorption rate ( SAR10g ) and SAR efficiency ( B1+ / SAR10g ).

RF heating measurement using MR thermometry and field monitoring: Methodological considerations and first in vivo results.

PURPOSE: A MR thermometry (MRT) method with field monitoring is proposed to improve the measurement of small temperature variations induced in brain MRI exams. METHODS: MR thermometry experiments were performed at 7 Tesla with concurrent field monitoring and RF heating. Images were reconstructed with nominal k-space trajectories and with first-order spherical harmonics correction. Experiments were performed in vitro with deliberate field disturbances and on an anesthetized macaque in 2 different specific absorption rate regimes, that is, at 50% and 100% of the maximal specific absorption rate level allowed in the International Electrotechnical Commission normal mode of operation. Repeatability was assessed by running a second separate session on the same animal. RESULTS: Inclusion of magnetic field fluctuations in the reconstruction improved temperature measurement accuracy in vitro down to 0.02°C. Measurement precision in vivo was on the order of 0.15°C in areas little affected by motion. In the same region, temperature increase reached 0.5 to 0.8°C after 20 min of heating at 100% specific absorption rates and followed a rough factor of 2 with the 50% specific absorption rate scans. A horizontal temperature plateau, as predicted by Pennes bioheat model with thermal constants from the literature and constant blood temperature assumption, was not observed. CONCLUSION: Inclusion of field fluctuations in image reconstruction was beneficial for the measurement of small temperature rises encountered in standard brain exams. More work is needed to correct for motion-induced field disturbances to extract reliable temperature maps.

What can 7T sodium MRI tell us about cellular energy depletion and neurotransmission in Alzheimer's disease?

The pathophysiological processes underlying the development and progression of Alzheimer's disease (AD) on the neuronal level are still unclear. Previous research has hinted at metabolic energy deficits and altered sodium homeostasis with impaired neuronal function as a potential metabolic marker relevant for neurotransmission in AD. Using sodium (23 Na) magnetic resonance (MR) imaging on an ultra-high-field 7 Tesla MR scanner, we found increased cerebral tissue sodium concentration (TSC) in 17 biomarker-defined AD patients compared to 22 age-matched control subjects in vivo. TSC was highly discriminative between controls and early AD stages and was predictive for cognitive state, and associated with regional tau load assessed with flortaucipir-positron emission tomography as a possible mediator of TSC-associated neurodegeneration. TSC could therefore serve as a non-invasive, stage-dependent, metabolic imaging marker. Setting a focus on cellular metabolism and potentially disturbed interneuronal communication due to energy-dependent altered cell homeostasis could hamper progressive cognitive decline by targeting these processes in future interventions.

Mapping subcomponents of numerical cognition in relation to functional and anatomical landmarks of human parietal cortex.

Human functional imaging has identified the middle part of the intraparietal sulcus (IPS) as an important brain substrate for different types of numerical tasks. This area is often equated with the macaque ventral intraparietal area (VIP) where neuronal selectivity for non-symbolic numerical stimuli (sets of items) is found. However, the low spatial resolution and whole-brain averaging analysis performed in most fMRI studies limit the extent to which an exact correspondence of activations in different numerical tasks with specific sub-regions of the IPS can be established. Here we acquired high-resolution 7T fMRI data in a group of human adults and related the activations in several numerical contrasts (implying different numerical stimuli and tasks) to anatomical and functional landmarks on the cortical surface. Our results reveal a functional heterogeneity within human intraparietal cortex where the retinotopic visual field maps in superior/medial parts of the IPS and superior parietal gyrus respond preferentially to the visual processing of concrete sets of items (over single Arabic numerals), whereas lateral/inferior parts of the IPS are predominantly recruited during numerical operations such as calculation and quantitative comparison. Since calculation and comparison-related activity fell mainly outside the retinotopic visual field maps considered the human functional equivalent of the monkey VIP/LIP complex, the areas most activated during such numerical operations in humans are likely different from VIP.

Simultaneous proton density, T1 , T2 , and flip-angle mapping of the brain at 7 T using multiparametric 3D SSFP imaging and parallel-transmission universal pulses.

PURPOSE: Performing simultaneous quantitative MRI at ultrahigh field is challenging, as B0 and B1+ heterogeneities as well as specific absorption rate increase. Too large deviations of flip angle from the target can induce biases and impair SNR in the quantification process. In this work, we use calibration-free parallel transmission, a dedicated pulse-sequence parameter optimization and signal fitting to recover 3D proton density, flip angle, T1 , and T2 maps over the whole brain, in a clinically suitable time. METHODS: Eleven optimized contrasts were acquired with an unbalanced SSFP sequence by varying flip-angle amplitude and RF phase-cycling increment, at a 1.0 × 1.0 × 3.0 mm3 resolution. Acquisition time was 16 minutes 36 seconds for the whole brain. Parallel transmission and universal pulses were used to mitigate B1+ heterogeneity, to improve the results' reliability over 6 healthy volunteers (3 females/3 males, age 22.6 ± 2.7 years old). Quantification of the physical parameters was performed by fitting the acquired contrasts to the simulated ones using the Bloch-Torrey equations with a realistic diffusion coefficient. RESULTS: Whole-brain 3D maps of effective flip angle, proton density, and relaxation times were estimated. Parallel transmission improved the robustness of the results at 7 T. Results were in accordance with literature and with measurements from standard methods. CONCLUSION: These preliminary results show robust proton density, flip angle, T1 , and T2 map retrieval. Other parameters, such as ADC, could be assessed. With further optimization in the acquisition, scan time could be reduced and spatial resolution increased to bring this multiparametric quantification method to clinical research routine at 7 T.

3D variable-density SPARKLING trajectories for high-resolution T2*-weighted magnetic resonance imaging.

We have recently proposed a new optimization algorithm called SPARKLING (Spreading Projection Algorithm for Rapid K-space sampLING) to design efficient compressive sampling patterns for magnetic resonance imaging (MRI). This method has a few advantages over conventional non-Cartesian trajectories such as radial lines or spirals: i) it allows to sample the k-space along any arbitrary density while the other two are restricted to radial densities and ii) it optimizes the gradient waveforms for a given readout time. Here, we introduce an extension of the SPARKLING method for 3D imaging by considering both stacks-of-SPARKLING and fully 3D SPARKLING trajectories. Our method allowed to achieve an isotropic resolution of 600 µm in just 45 seconds for T2∗-weighted ex vivo brain imaging at 7 Tesla over a field-of-view of 200 × 200 × 140 mm3 . Preliminary in vivo human brain data shows that a stack-of-SPARKLING is less subject to off-resonance artifacts than a stack-of-spirals.

Post-contrast 3D T1-weighted TSE MR sequences (SPACE, CUBE, VISTA/BRAINVIEW, isoFSE, 3D MVOX): Technical aspects and clinical applications.

Post-contrast three-dimensional T1-weighted imaging of the brain is widely used for a broad range of vascular, inflammatory or tumoral diseases. The variable flip angle 3D TSE sequence is now available from several manufacturers (CUBE, General Electric; SPACE, Siemens; VISTA/BRAINVIEW, Philips; isoFSE, Itachi; 3D MVOX, Canon). Compared to gradient-echo (GRE) techniques, 3D TSE offers the advantages of useful image contrasts and reduction of artifacts from static field inhomogeneity. However, the respective role of 3D TSE and GRE MR sequences remains to be elucidated, particularly in the setting of post-contrast imaging. The purpose of this review was (1) to describe the technical aspects of 3D TSE sequences, (2) to illustrate the main clinical applications of the post-contrast 3D T1-w TSE sequence through clinical cases, (3) to discuss the respective role of post-contrast 3D TSE and GRE imaging in the field of neuroimaging.

SPARKLING: variable-density k-space filling curves for accelerated T2* -weighted MRI.

PURPOSE: To present a new optimition-driven design of optimal k-space trajectories in the context of compressed sensing: Spreading Projection Algorithm for Rapid K-space sampLING (SPARKLING). THEORY: The SPARKLING algorithm is a versatile method inspired from stippling techniques that automatically generates optimized sampling patterns compatible with MR hardware constraints on maximum gradient amplitude and slew rate. These non-Cartesian sampling curves are designed to comply with key criteria for optimal sampling: a controlled distribution of samples (e.g., variable density) and a locally uniform k-space coverage. METHODS: Ex vivo and in vivo prospective T2* -weighted acquisitions were performed on a 7-Tesla scanner using the SPARKLING trajectories for various setups and target densities. Our method was compared to radial and variable-density spiral trajectories for high-resolution imaging. RESULTS: Combining sampling efficiency with compressed sensing, the proposed sampling patterns allowed up to 20-fold reductions in MR scan time (compared to fully sampled Cartesian acquisitions) for two-dimensional T2* -weighted imaging without deterioration of image quality, as demonstrated by our experimental results at 7 Tesla on in vivo human brains for a high in-plane resolution of 390 µm. In comparison to existing non-Cartesian sampling strategies, the proposed technique also yielded superior image quality. CONCLUSIONS: The proposed optimization-driven design of k-space trajectories is a versatile framework that is able to enhance MR sampling performance in the context of compressed sensing.