Design of calibration-free RF pulses for T 2 $$ {}_2 $$ -weighted single-slab 3D turbo-spin-echo sequences at 7T utilizing parallel transmission.

PURPOSE: T 2 $$ {}_2 $$ -weighted turbo-spin-echo (TSE) sequences are a fundamental technique in brain imaging but suffer from field inhomogeneities at ultra-high fields. Several methods have been proposed to mitigate the problem, but were limited so far to nonselective three-dimensional (3D) measurements, making short acquisitions difficult to achieve when targeting very high resolution images, or needed additional calibration procedures, thus complicating their application. METHODS: Slab-selective excitation pulses were designed for flexible placement utilizing the concept of k T $$ {}_T $$ -spokes. Phase-coherent refocusing universal pulses were subsequently optimized with the Gradient Ascent Pulse Engineering algorithm and tested in vivo for improved signal homogeneity. RESULTS: Implemented within a 3D variable flip angle TSE sequence, these pulses led to a signal-to-noise ratio (SNR) improvement ranging from 10% to 30% compared to a two-dimensional (2D) T2w TSE sequence employing B 1 + $$ {\mathrm{B}}_1^{+} $$ -shimmed pulses. B 1 + $$ {\mathrm{B}}_1^{+} $$ field inhomogeneities could be mitigated and artifacts from B 0 $$ {\mathrm{B}}_0 $$ deviations reduced. The concept of universal pulses was successfully applied. CONCLUSION: We present a pulse design method which provides a set of calibration-free universal pulses (UPs) for slab-selective excitation and phase-coherent refocusing in slab-selective TSE sequences.

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.

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.

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.

Optimized interferometric encoding of presaturated TurboFLASH B1 mapping for parallel transmission MRI at 7 T: Preliminary application for quantitative T1 mapping in the spinal cord.

PURPOSE: The acquisition of accurate B1 maps is critical for parallel transmit techniques (pTx). The presaturated turboFLASH (satTFL) method has been widely used in combination with interferometric encoding to provide robust and fast B1 maps. However, typical encodings, mostly evaluated on brain, do not necessarily fit all coils and organs. In this work, we evaluated and improved the accuracy of the satTFL for cervical spine at 7 T, proposing a novel interferometric encoding optimization. The benefits of such improvements were investigated in an exploratory study of quantitative T1 mapping with pTx-MP2RAGE. METHODS: Global optimization of interferometric encoding was implemented by simulating the ability of the satTFL to reconstruct B1 maps, with varying encoding and inclusion of complex noise, inside a region of interest covering the cervical spine. The performance of satTFL before and after optimization was compared to actual flip angle imaging. Optimized and non-optimized B1 maps were then used to calculate pTx pulses for MP2RAGE T1 mapping. RESULTS: Interferometric encoding optimization resulted in satTFL closer to actual flip angle imaging, with substantial gain of signal in regions where non-optimized satTFL could fail. T1 maps measured with non-adiabatic pTx pulses were closer to standard non-pTx results (which used adiabatic pulses) when using optimized-satTFL, with substantially lower specific absorption rate. CONCLUSION: Optimization of the satTFL interferometric encoding improves B1 maps in the spinal cord, in particular in low SNR regions. A linear correction of the satTFL was additionally shown to be required. The method was successfully used for quantitative phantom and in vivo T1 mapping, showing improved results compared to non-optimized satTFL thanks to improved pTx-pulse generation.

Improved detection of juxtacortical lesions using highly accelerated double inversion-recovery MRI in patients with multiple sclerosis.

PURPOSE: The purpose of this study was to compare a highly-accelerated double inversion recovery (fast-DIR) sequence using a recent parallel imaging technique (CAIPIRINHA) with a conventional DIR (conv-DIR) sequence for image quality and the detection of juxtacortical and infratentorial multiple sclerosis (MS) lesions. MATERIALS AND METHODS: A total of 38 patients with MS who underwent brain MRI at 3 T between 2020 and 2021 were included. There were 27 women and 12 men with a mean age of 40 ± 12.8 (standard deviation) years (range: 20-59 years). All patients underwent conv-DIR sequence and fast-DIR sequence. Fast-DIR was obtained with a T2-preparation module to improve contrast and an iterative denoising algorithm to compensate noise enhancement. Two blinded readers reported the number of juxtacortical and infratentorial MS lesions for fast-DIR and conv-DIR, confirmed by further consensus reading that was used as the standard of reference. Image quality and contrast were evaluated for fast-DIR and conv-DIR sequences. Comparisons between fast-DIR and conv-DIR sequences were performed using Wilcoxon test and Lin concordance correlation coefficient. RESULTS: Thirty-eight patients were analyzed. Fast-DIR imaging allowed detection of 289 juxtacortical lesions vs. 238 with conv-DIR, corresponding to a significant improved detection rate with fast-DIR (P < 0.001). Conversely, 117 infratentorial lesions were detected with conv-DIR sequence vs. 80 with fast-DIR sequence (P < 0.001). Inter-observer agreement for lesion detection with fast-DIR and conv-DIR was very high (Lin concordance correlation coefficient ranging between 0.86 and 0.96). CONCLUSION: Fast-DIR improves the detection of juxtacortical MS lesions, but is limited for the detection of infratentorial MS lesions.

Electromagnetic and RF pulse design simulation based optimization of an eight-channel loop array for 11.7T brain imaging.

PURPOSE: Optimization of transmit array performance is crucial in ultra-high-field MRI scanners such as 11.7T because of the increased RF losses and RF nonuniformity. This work presents a new workflow to investigate and minimize RF coil losses, and to choose the optimum coil configuration for imaging. METHODS: An 8-channel transceiver loop-array was simulated to analyze its loss mechanism at 499.415 MHz. A folded-end RF shield was developed to limit radiation loss and improve the B 1 + $$ {B}_1^{+} $$ efficiency. The coil element length, and the shield diameter and length were further optimized using electromagnetic (EM) simulations. The generated EM fields were used to perform RF pulse design (RFPD) simulations under realistic constraints. The chosen coil design was constructed to demonstrate performance equivalence in bench and scanner measurements. RESULTS: The use of conventional RF shields at 11.7T resulted in significantly high radiation losses of 18.4%. Folding the ends of the RF shield combined with optimizing the shield diameter and length increased the absorbed power in biological tissue and reduced the radiation loss to 2.4%. The peak B 1 + $$ {B}_1^{+} $$ of the optimal array was 42% more than the reference array. Phantom measurements validated the numerical simulations with a close match of within 4% of the predicted B 1 + $$ {B}_1^{+} $$ . CONCLUSION: A workflow that combines EM and RFPD simulations to numerically optimize transmit arrays was developed. Results have been validated using phantom measurements. Our findings demonstrate the need for optimizing the RF shield in conjunction with array element design to achieve efficient excitation at 11.7T.

Accurate Diagnosis of Cortical and Infratentorial Lesions in Multiple Sclerosis Using Accelerated Fluid and White Matter Suppression Imaging.

OBJECTIVES: The precise location of multiple sclerosis (MS) cortical lesions can be very challenging at 3 T, yet distinguishing them from subcortical lesions is essential for the diagnosis and prognosis of the disease. Compressed sensing-accelerated fluid and white matter suppression imaging (CS-FLAWS) is a new magnetic resonance imaging sequence derived from magnetization-prepared 2 rapid acquisition gradient echo with promising features for the detection and classification of MS lesions. The objective of this study was to compare the diagnostic performances of CS-FLAWS (evaluated imaging) and phase sensitive inversion recovery (PSIR; reference imaging) for classification of cortical lesions (primary objective) and infratentorial lesions (secondary objective) in MS, in combination with 3-dimensional (3D) double inversion recovery (DIR). MATERIALS AND METHODS: Prospective 3 T scans (MS first diagnosis or follow-up) acquired between March and August 2021 were retrospectively analyzed. All underwent 3D CS-FLAWS, axial 2D PSIR, and 3D DIR. Double-blinded reading sessions exclusively in axial plane and final consensual reading were performed to assess the number of cortical and infratentorial lesions. Wilcoxon test was used to compare the 2 imaging datasets (FLAWS + DIR and PSIR + DIR), and intraobserver and interobserver agreement was assessed using the intraclass correlation coefficient. RESULTS: Forty-two patients were analyzed (38 with relapsing-remitting MS, 29 women, 42.7 ± 12.6 years old). Compressed sensing-accelerated FLAWS allowed the identification of 263 cortical lesions versus 251 with PSIR ( P = 0.74) and 123 infratentorial lesions versus 109 with PSIR ( P = 0.63), corresponding to a nonsignificant difference between the 2 sequences. Compressed sensing-accelerated FLAWS exhibited fewer false-negative findings than PSIR either for cortical lesions (1 vs 13; P < 0.01) or infratentorial lesions (1 vs 15; P < 0.01). No false-positive findings were found with any of the 2 sequences. Diagnostic confidence was high for each contrast. CONCLUSION: Three-dimensional CS-FLAWS is as accurate as 2D PSIR imaging for classification of cortical and infratentorial MS lesions, with fewer false-negative findings, opening the way to a reliable full brain MS exploration in a clinically acceptable duration (5 minutes 15 seconds).

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.

Fast T2-weighted liver MRI: Image quality and solid focal lesions conspicuity using a deep learning accelerated single breath-hold HASTE fat-suppressed sequence.

PURPOSE: Acceleration of MRI acquisitions and especially of T2-weighted sequences is essential to reduce the duration of MRI examinations but also kinetic artifacts in liver imaging. The purpose of this study was to compare the acquisition time and the image quality of a single-shot fat-suppressed turbo spin-echo (TSE) T2-weighted sequence with deep learning reconstruction (HASTEDL) with that of a fat-suppressed T2-weighted BLADE TSE sequence in patients with focal liver lesions. MATERIALS AND METHODS: Ninety-five patients (52 men, 43 women; mean age: 61 ± 14 [SD]; age range: 28-87 years) with 42 focal liver lesions (17 hepatocellular carcinomas, 10 sarcoidosis lesions, 9 myeloma lesions, 3 liver metastases and 3 focal nodular hyperplasias) who underwent liver MRI at 1.5 T including HASTEDL and BLADE sequences were retrospectively included. Overall image quality, noise level in the liver, lesion conspicuity and sharpness of liver lesion contours were assessed by two independent readers. Liver signal-to-noise ratio (SNR) and lesion contrast-to-noise ratio (CNR) were measured and compared between the two sequences, as well as the mean duration of the sequences (Student t-test or Wilcoxon test for paired data). RESULTS: Median overall quality on HASTEDL images (3; IQR: 3, 3) was significantly greater than that on BLADE images (2; IQR: 1, 3) (P < 0.001). Median noise level in the liver on HASTEDL images (0; IQR: 0, 0.5) was significantly lower than that on BLADE images (1; IQR: 1, 2) (P < 0.001). On HASTEDL images, mean liver SNR (107.3 ± 39.7 [SD]) and mean focal liver lesion CNR (87.0 ± 76.6 [SD]) were significantly greater than those on BLADE images (67.1 ± 23.8 [SD], P < 0.001 and 48.6 ± 43.9 [SD], P = 0.027, respectively). Acquisition time was significantly shorter with the HASTEDL sequence (18 ± [0] s; range: 18-18 s) compared to BLADE sequence (152 ± 47 [SD] s; range: 87-263 s) (P < 0.001). CONCLUSION: By comparison with the BLADE sequence, HASTEDL sequence significantly reduces acquisition time while improving image quality, liver SNR and focal liver lesions CNR.

3-Dimensional Fluid and White Matter Suppression Magnetic Resonance Imaging Sequence Accelerated With Compressed Sensing Improves Multiple Sclerosis Cervical Spinal Cord Lesion Detection Compared With Standard 2-Dimensional Imaging.

OBJECTIVES: Fluid and white matter suppression (FLAWS) is a recently proposed magnetic resonance sequence derived from magnetization-prepared 2 rapid acquisition gradient-echo providing 2 coregistered datasets with white matter- and cerebrospinal fluid-suppressed signal, enabling synthetic imaging with amplified contrast. Although these features are high potential for brain multiple sclerosis (MS) imaging, spinal cord has never been evaluated with this sequence to date. The objective of this work was therefore to assess diagnostic performance and self-confidence provided by compressed-sensing (CS) 3-dimensional (3D) FLAWS for cervical MS lesion detection on a head scan that includes the cervical cord without changing standard procedures. MATERIALS AND METHODS: Prospective 3 T scans (MS first diagnosis or follow-up) acquired between 2019 and 2020 were retrospectively analyzed. All patients underwent 3D CS-FLAWS (duration: 5 minutes 40 seconds), axial T 2 turbo spin echo covering cervical spine from cervicomedullary junction to the same inferior level as FLAWS, and sagittal cervical T 2 /short tau inversion recovery imaging. Two readers performed a 2-stage double-blind reading, followed by consensus reading. Wilcoxon tests were used to compare the number of detected spinal cord lesions and the reader's diagnostic self-confidence when using FLAWS versus the reference 2D T 2 -weighted imaging. RESULTS: Fifty-eight patients were included (mean age, 40 ± 13 years, 46 women, 7 ± 6 years mean disease duration). The CS-FLAWS detected significantly more lesions than the reference T 2 -weighted imaging (197 vs 152 detected lesions, P < 0.001), with a sensitivity of 98% (T 2 -weighted imaging sensitivity: 90%) after consensual reading. Considering the subgroup of patients who underwent sagittal T2 + short tau inversion recovery imaging (Magnetic Resonance Imaging for Multiple Sclerosis subgroup), +250% lesions were detected with FLAWS (63 vs 25 lesions detected, P < 0.001). Mean reading self-confidence was significantly better with CS-FLAWS (median, 5 [interquartile range, 1] [no doubt for diagnosis] vs 4 [interquartile range, 1] [high confidence]; P < 0.001). CONCLUSIONS: Imaging with CS-FLAWS provides an improved cervical spinal cord exploration for MS with increased self-confidence compared with conventional T 2 -weighted imaging, in a clinically acceptable time.

Optimizing Full 3D SPARKLING Trajectories for High-Resolution Magnetic Resonance Imaging.

The Spreading Projection Algorithm for Rapid K-space sampLING, or SPARKLING, is an optimization-driven method that has been recently introduced for accelerated 2D MRI using compressed sensing. It has then been extended to address 3D imaging using either stacks of 2D sampling patterns or a local 3D strategy that optimizes a single sampling trajectory at a time. 2D SPARKLING actually performs variable density sampling (VDS) along a prescribed target density while maximizing sampling efficiency and meeting the gradient-based hardware constraints. However, 3D SPARKLING has remained limited in terms of acceleration factors along the third dimension if one wants to preserve a peaky point spread function (PSF) and thus good image quality. In this paper, in order to achieve higher acceleration factors in 3D imaging while preserving image quality, we propose a new efficient algorithm that performs optimization on full 3D SPARKLING. The proposed implementation based on fast multipole methods (FMM) allows us to design sampling patterns with up to 107 k-space samples, thus opening the door to 3D VDS. We compare multi-CPU and GPU implementations and demonstrate that the latter is optimal for 3D imaging in the high-resolution acquisition regime ( 600µ m isotropic). Finally, we show that this novel optimization for full 3D SPARKLING outperforms stacking strategies or 3D twisted projection imaging through retrospective and prospective studies on NIST phantom and in vivo brain scans at 3 Tesla taking the particular case of T2 *-w imaging. Overall the proposed method allows for 2.5-3.75x shorter scan times compared to GRAPPA-4 parallel imaging acquisition at 3 Tesla without compromising image quality.

Standardized universal pulse: A fast RF calibration approach to improve flip angle accuracy in parallel transmission.

PURPOSE: In parallel transmission (pTX), subject-tailored RF pulses allow achieving excellent flip angle (FA) accuracy but often require computationally extensive online optimizations, precise characterization of the static field ( ΔB0 ), and the transmit RF field ( B1+ ) distributions. This costs time and requires expertise from the MR user. Universal pulses (UPs) have been proposed to reduce this burden, yet, with a penalty in FA accuracy. This study introduces the concept of standardized universal pulses (SUPs), where pulses are designed offline and adjusted to the subject through a fast online calibration scan. METHODS: A SUP is designed offline using a so-called standardized database, wherein each B1+ map has been normalized to a reference transmit RF field distribution. When scanning a new subject, a 3-slice B1+ acquisition (scan time <10  s) is performed and used to adjust the SUP to the subject through a linear transform. SUP performance was assessed at 7T with simulations by computing the FA-normalized root mean square error (FA-NRMSE) and the FA pattern stability as measured by the average and coefficient of variation of the FA across 15 control subjects, along with in vivo experiments using an MP2RAGE sequence implementing the SUP variant for the FLASH readout. RESULTS: Adjusted SUP improved the FA-NRMSE (8.8 % for UP vs. 7.1 % for adjusted SUP). Experimentally in vivo, this translated in an improved signal homogeneity and more accurate T1 quantification using MP2RAGE. CONCLUSION: The proposed SUP approach improves excitation accuracy (FA-NRMSE) while preserving the same offline pulse design principle as offered by UPs.

Automatic segmentation of deep grey nuclei using a high-resolution 7T magnetic resonance imaging atlas-Quantification of T1 values in healthy volunteers.

We present a new consensus atlas of deep grey nuclei obtained by shape-based averaging of manual segmentation of two experienced neuroradiologists and optimized from 7T MP2RAGE images acquired at (.6 mm)3 in 60 healthy subjects. A group-wise normalization method was used to build a high-contrast and high-resolution T1 -weighted brain template (.5 mm)3 using data from 30 out of the 60 controls. Delineation of 24 deep grey nuclei per hemisphere, including the claustrum and 12 thalamic nuclei, was then performed by two expert neuroradiologists and reviewed by a third neuroradiologist according to tissue contrast and external references based on the Morel atlas. Corresponding deep grey matter structures were also extracted from the Morel and CIT168 atlases. The data-derived, Morel and CIT168 atlases were all applied at the individual level using non-linear registration to fit the subject reference and to extract absolute mean quantitative T1 values derived from the 3D-MP2RAGE volumes, after correction for residual B1+ biases. Three metrics (the Dice and the volumetric similarity coefficients and a novel Hausdorff distance) were used to estimate the inter-rater agreement of manual MRI segmentation and inter-atlas variability, and these metrics were measured to quantify biases due to image registration, and their impact on the measurements of the quantitative T1 values was highlighted. This represents a fully automated segmentation process permitting the extraction of unbiased normative T1 values in a population of young healthy controls as a reference for characterizing subtle structural alterations of deep grey nuclei relevant to a range of neurological diseases.

T1-Based Synthetic Magnetic Resonance Contrasts Improve Multiple Sclerosis and Focal Epilepsy Imaging at 7 T.

OBJECTIVES: Ultra-high field magnetic resonance imaging (MRI) (≥7 T) is a unique opportunity to improve the clinical diagnosis of brain pathologies, such as multiple sclerosis or focal epilepsy. However, several shortcomings of 7 T MRI, such as radiofrequency field inhomogeneities, could degrade image quality and hinder radiological interpretation. To address these challenges, an original synthetic MRI method based on T1 mapping achieved with the magnetization-prepared 2 rapid acquisition gradient echo (MP2RAGE) sequence was developed. The radiological quality of on-demand T1-based contrasts generated by this technique was evaluated in multiple sclerosis and focal epilepsy imaging at 7 T. MATERIALS AND METHODS: This retrospective study was carried out from October 2017 to September 2019 and included 21 patients with different phenotypes of multiple sclerosis and 35 patients with focal epilepsy who underwent MRI brain examinations using a whole-body investigative 7 T magnetic resonance system. The quality of 2 proposed synthetic contrast images were assessed and compared with conventional images acquired at 7 T using the MP2RAGE sequence by 4 radiologists, evaluating 3 qualitative criteria: signal homogeneity, contrast intensity, and lesion visualization. Statistical analyses were performed on reported quality scores using Wilcoxon rank tests and further multiple comparisons tests. Intraobserver and interobserver reliabilities were calculated as well. RESULTS: Radiological quality scores were reported higher for synthetic images when compared with original images, regardless of contrast, pathologies, or raters considered, with significant differences found for all 3 criteria (P < 0.0001, Wilcoxon rank test). None of the 4 radiologists ever rated a synthetic image "markedly worse" than an original image. Synthetic images were rated slightly less satisfying for only 3 epileptic patients, without precluding lesion identification. CONCLUSION: T1-based synthetic MRI with the MP2RAGE sequence provided on-demand contrasts and high-quality images to the radiologist, facilitating lesion visualization in multiple sclerosis and focal epilepsy, while reducing the magnetic resonance examination total duration by removing an additional sequence.

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.

Intravoxel Incoherent Motion at 7 Tesla to quantify human spinal cord perfusion: limitations and promises.

PURPOSE: To develop a noninvasive technique to map human spinal cord (SC) perfusion in vivo. More specifically, to implement an intravoxel incoherent motion (IVIM) protocol at ultrahigh field for the human SC and assess parameters estimation errors. METHODS: Monte-Carlo simulations were conducted to assess estimation errors of 2 standard IVIM fitting approaches (two-step versus one-step fit) over the range of IVIM values reported for the human brain and for typical SC diffusivities. Required signal-to-noise ratio (SNR) was inferred for estimation of the parameters product, fIVIM D* (microvascular fraction times pseudo-diffusion coefficient), within 10% error margins. In-vivo IVIM imaging of the SC was performed at 7T in 6 volunteers. An image processing pipeline is proposed to generate IVIM maps and register them for an atlas-based region-wise analysis. RESULTS: Required b = 0 SNRs for 10% error estimation on fIVIM D* with the one-step fit were 159 and 185 for diffusion-encoding perpendicular and parallel to the SC axis, respectively. Average in vivo b = 0 SNR within cord was 141 ± 79, corresponding to estimation errors of 12.7% and 14.7% according to numerical simulations. Slice- and group-averaging reduced noise in IVIM maps, highlighting the difference in perfusion between gray and white matter. Mean ± standard deviation fIVIM and D* values across subjects within gray (respectively white) matter were 16.0 ± 1.7 (15.0 ± 1.6)% and 11.4 ± 2.9 (11.5 ± 2.4) × 10-3 mm2 /s. CONCLUSION: Single-subject data SNR at 7T was insufficient for reliable perfusion estimation. However, atlas-averaged IVIM maps highlighted the higher microvascular fraction of gray matter compared to white matter, providing first results of healthy human SC perfusion mapping with MRI.

Anterior fissure, central canal, posterior septum and more: New insights into the cervical spinal cord gray and white matter regional organization using T1 mapping at 7T.

T1 mapping lacks specificity toward a single particular biological feature, however it has the potential to discriminate spinal cord regional tissue organization and characterize tissue microstructural impairments occurring in neurodegenerative pathologies. In this exploratory work, T1 mapping of the cervical spinal cord with a 300-µm in-plane resolution was performed on fourteen healthy subjects at 7T, using the MP2RAGE sequence. Individual images from C1 to C7 vertebral levels provided a clear delineation of spinal cord anatomical details and substructures including motor columns within gray matter (GM) horns, anterior median fissure, central canal, ventral, lateral and dorsal white matter (WM) fasciculi, and posterior median septum. Group studies highlighted regional T1 differences between regions of interest so far hardly visible at lower spatial resolution. Two-dimensional averaged T1 maps and manual parcellation of GM and WM substructures were built based on these data. Benefiting from the very high spatial resolution achievable at ultra-high field for T1 mapping, this work contributes to improve the in vivo characterization of the cervical spinal cord. By allowing investigation within a wider range of functional regions, it also opens new perspectives for pathology diagnosis such as motor neuron disease, neuropathic pain or refined investigation of neurodegeneration.

Regional T1 mapping of the whole cervical spinal cord using an optimized MP2RAGE sequence.

The recently-proposed MP2RAGE sequence was purposely optimized for cervical spinal cord imaging at 3T. Sequence parameters were chosen to optimize gray/white matter T1 contrast with sub-millimetric resolution and scan-time < 10 min while preserving reliable T1 determination with minimal B1+ variation effects within a range of values compatible with pathologies and surrounding structures. Results showed good agreements with IR-based measurements, high MP2RAGE-based T1 reproducibility and preliminary evidences of age- and tract-related T1 variations in the healthy spinal cord.

While T1 measurements present multiple challenges (robustness, acquisition time), the recently proposed MP2RAGE sequence (magnetization-prepared two rapid acquisition gradient echoes) has opened new perspectives to characterize tissue microstructure changes occurring in a pathological or developmental context. Extensively used for brain studies, it was herein adapted to investigate the cervical spinal cord (SC) at 3 T. By integrating Bloch equations, the MP2RAGE sequence parameters were chosen to optimize SC gray matter/white matter (GM/WM) T1 contrast with sub-millimetric resolution, a scan time less than 10 min and a reliable T1 determination with minimal B1+ variation effect, within a range of values compatible with different pathologies and surrounding structures. The residual B1+ effect on T1 values was corrected using a look-up-table approach and B1+ mapping. The accuracy of B1+ -corrected T1 measurements was assessed on a phantom with respect to conventional inversion recovery. In vivo MP2RAGE acquisitions were performed on five young (28.8 ± 4.3 years old) and five elderly (60.2 ± 2.9 years old) volunteers and analyzed using a template-based approach. Phantom experiments led to high agreements between inversion-recovery spin-echo and MP2RAGE-based T1 values (R2  = 0.997). In vivo T1 values for cervical WM, anterior GM (aGM), posterior sensory tracts (PSTs) and lateral motor tracts (LMTs) were 917 ± 29 s, 934 ± 33 ms, 920 ± 37 ms and 877 ± 35 ms, respectively, with all subjects and cervical levels considered. Significant differences were observed between aGM and LMTs, and between LMTs and PSTs, in agreement with the literature. Repeated T1 measurements demonstrated high reproducibility of the MP2RAGE in the SC (variation coefficient < 5% in all regions of interest). Finally, preliminary assessment of age-related SC tissue microstructure variation additionally showed evidence of SC atrophy and slight trends of T1 decrease with age in all regions. Overall, this study shows that fast, robust and accurate sub-millimetric resolution T1 mapping in the cervical SC using the MP2RAGE sequence is possible, paving the way for future multi-centric and longitudinal clinical studies investigating the pathological cord.