Extended phase graphs: Dephasing, RF pulses, and echoes - pure and simple: Extended Phase Graphs

Efficient performance analysis and optimization of transient-state sequences for multiparametric magnetic resonance imaging

In transient-state multiparametric MRI sequences such as Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT), MR fingerprinting, or hybrid-state imaging, the flip angle pattern of the RF excitation varies over the sequence. This gives considerable freedom to choose an optimal pattern of flip angles. For pragmatic reasons, most optimization methodologies choose for a single-voxel approach (i.e., without taking the spatial encoding scheme into account). Particularly in MR-STAT, the context of spatial encoding is important. In the current study, we present a methodology, called BLock Analysis of a K-space-domain Jacobian (BLAKJac), which is sufficiently fast to optimize a sequence in the context of a predetermined phase-encoding pattern. Based on MR-STAT acquisitions and reconstructions, we show that sequences optimized using BLAKJac are more reliable in terms of actually achieved precision than conventional single-voxel-optimized sequences. In addition, BLAKJac provides analytical tools that give insights into the performance of the sequence in a very limited computation time. Our experiments are based on MR-STAT, but the theory is equally valid for other transient-state multiparametric methods.

A Framework for Simulating Cardiac MR Images With Varying Anatomy and Contrast

One of the limiting factors for the development and adoption of novel deep-learning (DL) based medical image analysis methods is the scarcity of labeled medical images. Medical image simulation and synthesis can provide solutions by generating ample training data with corresponding ground truth labels. Despite recent advances, generated images demonstrate limited realism and diversity. In this work, we develop a flexible framework for simulating cardiac magnetic resonance (MR) images with variable anatomical and imaging characteristics for the purpose of creating a diversified virtual population. We advance previous works on both cardiac MR image simulation and anatomical modeling to increase the realism in terms of both image appearance and underlying anatomy. To diversify the generated images, we define parameters: 1)to alter the anatomy, 2) to assign MR tissue properties to various tissue types, and 3) to manipulate the image contrast via acquisition parameters. The proposed framework is optimized to generate a substantial number of cardiac MR images with ground truth labels suitable for downstream supervised tasks. A database of virtual subjects is simulated and its usefulness for aiding a DL segmentation method is evaluated. Our experiments show that training completely with simulated images can perform comparable with a model trained with real images for heart cavity segmentation in mid-ventricular slices. Moreover, such data can be used in addition to classical augmentation for boosting the performance when training data is limited, particularly by increasing the contrast and anatomical variation, leading to better regularization and generalization. The database is publicly available at https://osf.io/bkzhm/ and the simulation code will be available at https://github.com/sinaamirrajab/CMRI.

Fast-spin-echo versus rapid gradient-echo for 3D magnetization-prepared acquisitions: Application to inhomogeneous magnetization transfer

PURPOSE

To evaluate the benefits of fast spin echo (FSE) imaging over rapid gradient-echo (RAGE) for magnetization-prepared inhomogeneous magnetization transfer (ihMT) imaging.

METHODS

A 3D FSE sequence was modified to include an ihMT preparation (ihMT-FSE) with an optional CSF suppression based on an inversion-recovery (ihMT-FLAIR). After numeric simulations assessing SNR benefits of FSE and the potential impact of an additional inversion-recovery, ihMT-RAGE, ihMT-FSE, and ihMT-FLAIR sequences were compared in a group of six healthy volunteers, evaluating image quality, thermal, and physiological noise as well as quantification using an ihMT saturation (ihMTsat) approach. A preliminary exploration in the cervical spinal cord was also conducted in a group of three healthy volunteers.

RESULTS

Several fold improvements in thermal SNR were observed with ihMT-FSE in agreement with numerical simulations. However, we observed significantly higher physiological noise in ihMT-FSE compared to ihMT-RAGE that was mitigated in ihMT-FLAIR, which provided the best total SNR (+74% and +49% compared to ihMT-RAGE in the white and gray matter, P ≤ 0.004). IhMTsat quantification was successful in all cases with strong correlation between all sequences (r2  > 0.75). Early experiments showed potential for spinal cord imaging.

CONCLUSIONS

FSE generally offers higher SNR compared to gradient-echo based acquisitions for magnetization-prepared contrasts as illustrated here in the case of ihMT. However, physiological noise has a significant effect, but an inversion-recovery-based CSF suppression was shown to be efficient in mitigating effects of CSF motion.

A method for measuring B0 field inhomogeneity using quantitative double-echo in steady-state

PURPOSE

To develop and validate a method forB 0 $$ {B}_0 $$ mapping for knee imaging using the quantitative Double-Echo in Steady-State (qDESS) exploiting the phase difference (Δ θ $$ \Delta \theta $$ ) between the two echoes acquired. Contrary to a two-gradient-echo (2-GRE) method,Δ θ $$ \Delta \theta $$ depends only on the first echo time.

METHODS

Bloch simulations were applied to investigate robustness to noise of the proposed methodology and all imaging studies were validated with phantoms and in vivo simultaneous bilateral knee acquisitions. Two phantoms and five healthy subjects were scanned using qDESS, water saturation shift referencing (WASSR), and multi-GRE sequences.Δ B 0 $$ \Delta {B}_0 $$ maps were calculated with the qDESS and the 2-GRE methods and compared against those obtained with WASSR. The comparison was quantitatively assessed exploiting pixel-wise difference maps, Bland-Altman (BA) analysis, and Lin's concordance coefficient (ρ c $$ {\rho}_c $$ ). For in vivo subjects, the comparison was assessed in cartilage using average values in six subregions.

RESULTS

The proposed method for measuringΔ B 0 $$ \Delta {B}_0 $$ inhomogeneities from a qDESS acquisition providedΔ B 0 $$ \Delta {B}_0 $$ maps that were in good agreement with those obtained using WASSR.Δ B 0 $$ \Delta {B}_0 $$ ρ c $$ {\rho}_c $$ values were≥ $$ \ge $$ 0.98 and 0.90 in phantoms and in vivo, respectively. The agreement between qDESS and WASSR was comparable to that of a 2-GRE method.

CONCLUSION

The proposed method may allow B0 correction for qDESST 2 $$ {T}_2 $$ mapping using an inherently co-registeredΔ B 0 $$ \Delta {B}_0 $$ map without requiring an additional B0 measurement sequence. More generally, the method may help shorten knee imaging protocols that require an auxiliaryΔ B 0 $$ \Delta {B}_0 $$ map by exploiting a qDESS acquisition that also providesT 2 $$ {T}_2 $$ measurements and high-quality morphological imaging.

Quantitative MRI in cardiometabolic disease: From conventional cardiac and liver tissue mapping techniques to multi-parametric approaches

Cardiometabolic disease refers to the spectrum of chronic conditions that include diabetes, hypertension, atheromatosis, non-alcoholic fatty liver disease, and their long-term impact on cardiovascular health. Histological studies have confirmed several modifications at the tissue level in cardiometabolic disease. Recently, quantitative MR methods have enabled non-invasive myocardial and liver tissue characterization. MR relaxation mapping techniques such as T₁, T_1ρ, T₂ and T₂* provide a pixel-by-pixel representation of the corresponding tissue specific relaxation times, which have been shown to correlate with fibrosis, altered tissue perfusion, oedema and iron levels. Proton density fat fraction mapping approaches allow measurement of lipid tissue in the organ of interest. Several studies have demonstrated their utility as early diagnostic biomarkers and their potential to bear prognostic implications. Conventionally, the quantification of these parameters by MRI relies on the acquisition of sequential scans, encoding and mapping only one parameter per scan. However, this methodology is time inefficient and suffers from the confounding effects of the relaxation parameters in each single map, limiting wider clinical and research applications. To address these limitations, several novel approaches have been proposed that encode multiple tissue parameters simultaneously, providing co-registered multiparametric information of the tissues of interest. This review aims to describe the multi-faceted myocardial and hepatic tissue alterations in cardiometabolic disease and to motivate the application of relaxometry and proton-density cardiac and liver tissue mapping techniques. Current approaches in myocardial and liver tissue characterization as well as latest technical developments in multiparametric quantitative MRI are included. Limitations and challenges of these novel approaches, and recommendations to facilitate clinical validation are also discussed.

Confounder-corrected T1 mapping in the liver through simultaneous estimation of T1 , PDFF, R 2 * , and B 1 + in a single breath-hold acquisition

PURPOSE

Quantitative volumetric T₁ mapping in the liver has the potential to aid in the detection, diagnosis, and quantification of liver fibrosis, inflammation, and spatially resolved liver function. However, accurate measurement of hepatic T₁ is confounded by the presence of fat and inhomogeneous B 1 + $$ {B}_1^{+} $$ excitation. Furthermore, scan time constraints related to respiratory motion require tradeoffs of reduced volumetric coverage and/or increased acquisition time. This work presents a novel 3D acquisition and estimation method for confounder-corrected T₁ measurement over the entire liver within a single breath-hold through simultaneous estimation of T₁ , fat and B 1 + $$ {B}_1^{+} $$ .

THEORY AND METHODS

The proposed method combines chemical shift encoded MRI and variable flip angle MRI with a B 1 + $$ {B}_1^{+} $$ mapping technique to enable confounder-corrected T₁ mapping. The method was evaluated theoretically and demonstrated in both phantom and in vivo acquisitions at 1.5 and 3.0T. At 1.5T, the method was evaluated both pre- and post- contrast enhancement in healthy volunteers.

RESULTS

The proposed method demonstrated excellent linear agreement with reference inversion-recovery spin-echo based T₁ in phantom acquisitions at both 1.5 and 3.0T, with minimal bias (5.2 and 45 ms, respectively) over T₁ ranging from 200-1200 ms. In vivo results were in general agreement with reference saturation-recovery based 2D T₁ maps (SMART₁ Map, GE Healthcare).

CONCLUSION

The proposed 3D T₁ mapping method accounts for fat and B 1 + $$ {B}_1^{+} $$ confounders through simultaneous estimation of T₁ , B 1 + $$ {B}_1^{+} $$ , PDFF and R 2 * $$ {R}_2^{\ast } $$ . It demonstrates strong linear agreement with reference T₁ measurements, with low bias and high precision, and can achieve full liver coverage in a single breath-hold.

Muscle quantitative MRI as a novel biomarker in hereditary transthyretin amyloidosis with polyneuropathy: a cross-sectional study

BACKGROUND

The development of reproducible and sensitive outcome measures has been challenging in hereditary transthyretin (ATTRv) amyloidosis. Recently, quantification of intramuscular fat by magnetic resonance imaging (MRI) has proven as a sensitive marker in patients with other genetic neuropathies. The aim of this study was to investigate the role of muscle quantitative MRI (qMRI) as an outcome measure in ATTRv.

METHODS

Calf- and thigh-centered multi-echo T2-weighted spin-echo and gradient-echo sequences were obtained in patients with ATTRv amyloidosis with polyneuropathy (n = 24) and healthy controls (n = 12). Water T2 (wT2) and fat fraction (FF) were calculated. Neurological assessment was performed in all ATTRv subjects. Quantitative MRI parameters were correlated with clinical and neurophysiological measures of disease severity.

RESULTS

Quantitative imaging revealed significantly higher FF in lower limb muscles in patients with ATTRv amyloidosis compared to controls. In addition, wT2 was significantly higher in ATTRv patients. There was prominent involvement of the posterior compartment of the thighs. Noticeably, FF and wT2 did not exhibit a length-dependent pattern in ATTRv patients. MRI biomarkers correlated with previously validated clinical outcome measures, Polyneuropathy Disability scoring system, Neuropathy Impairment Score (NIS) and NIS-lower limb, and neurophysiological parameters of axonal damage regardless of age, sex, treatment and TTR mutation.

CONCLUSIONS

Muscle qMRI revealed significant difference between ATTRv and healthy controls. MRI biomarkers showed high correlation with clinical and neurophysiological measures of disease severity making qMRI as a promising tool to be further investigated in longitudinal studies to assess its role at monitoring onset, progression, and therapy efficacy for future clinical trials on this treatable condition.

Radiomics and machine learning applied to STIR sequence for prediction of quantitative parameters in facioscapulohumeral disease

Purpose

Quantitative Muscle MRI (qMRI) is a valuable and non-invasive tool to assess disease involvement and progression in neuromuscular disorders being able to detect even subtle changes in muscle pathology. The aim of this study is to evaluate the feasibility of using a conventional short-tau inversion recovery (STIR) sequence to predict fat fraction (FF) and water T2 (wT2) in skeletal muscle introducing a radiomic workflow with standardized feature extraction combined with machine learning algorithms.

Methods

Twenty-five patients with facioscapulohumeral muscular dystrophy (FSHD) were scanned at calf level using conventional STIR sequence and qMRI techniques. We applied and compared three different radiomics workflows (WF1, WF2, WF3), combined with seven Machine Learning regression algorithms (linear, ridge and lasso regression, tree, random forest, k-nearest neighbor and support vector machine), on conventional STIR images to predict FF and wT2 for six calf muscles.

Results

The combination of WF3 and K-nearest neighbor resulted to be the best predictor model of qMRI parameters with a mean absolute error about ± 5 pp for FF and ± 1.8 ms for wT2.

Conclusion

This pilot study demonstrated the possibility to predict qMRI parameters in a cohort of FSHD subjects starting from conventional STIR sequence.

Optimisation of data acquisition towards continuous cardiac Magnetic Resonance Fingerprinting applications

PURPOSE

Assess and optimise acquisition parameters for continuous cardiac Magnetic Resonance Fingerprinting (MRF).

METHODS

Different acquisition schemes (flip angle amplitude, lobe size, T2-preparation pulses) for cardiac MRF were assessed in simulations and phantom and demonstrated in one healthy volunteer. Three different experimental designs were evaluated using central composite and fractional factorial designs. Relative errors for T1 and T2 were calculated for a wide range of realistic T1 and T2 value combinations. The effect of different designs on the accuracy of T1 and T2 was assessed using response surface modelling and Cohen's f calculations.

RESULTS

Larger flip angle amplitudes lead to an improvement of T2 accuracy and precision for simulations and phantom experiments. Similar effects could also be shown qualitatively in in-vivo scans. Accuracy and precision of T1 were robust to different design parameters with improved values for faster flip angle variation. Cohen's f showed that T2-preparation pulses influence the accuracy of T2. The number of pulses used is the most important parameter. Without T2-preparation pulses, RMSE were 3.0 ± 8.09 % for T1 and 16.24 ± 14.47 % for T2. Using those pulses reduced the RMSE to 2.3 ± 8.4 % for T1 and 14.11 ± 13.46 % for T2. Nonetheless, even if the improvement is significant, RMSE are still too high for reliable quantification.

CONCLUSION

In contrast to previous study using triggered MRF sequences using < 30° flip angles, large flip angle amplitudes led to better results for continuous cardiac MRF sequences. T2-preparation pulse can improve the accuracy of T2 estimation but lead to longer scan times.

3D magnetic resonance fingerprinting for rapid simultaneous T1, T2, and T1ρ volumetric mapping of human articular cartilage at 3 T

Quantitative MRI can detect early biochemical changes in cartilage; however, the conventional techniques only measure one parameter (e.g., T1 , T2 , and T1ρ ) at a time while also being comparatively slow. We implemented a 3D magnetic resonance fingerprinting (3D-MRF) technique for simultaneous, volumetric mapping of T1 , T2 , and T1ρ in knee articular cartilage in under 9 min. It is evaluated on 11 healthy volunteers (mean age: 53 ± 9 years), five mild knee osteoarthritis (OA) patients (Kellgren-Lawrence (KL) score: 2, mean age: 60 ± 4 years), and the National Institute of Standards and Technology (NIST)/International Society for Magnetic Resonance in Medicine (ISMRM) system phantom. Proton density image, and T1 , T2, T1ρ relaxation times, and B1 + were estimated in the NIST/ISMRM system phantom as well as in the human knee medial and lateral femur, medial and lateral tibia, and patellar cartilage. The repeatability and reproducibility of the proposed technique were assessed in the phantom using analysis of the Bland-Altman plots. The intrasubject repeatability was assessed with the coefficient of variation (CV) and root mean square CV (rmsCV). The Mann-Whitney U test was used to assess the difference between healthy subjects and mild knee OA patients. The Bland-Altman plots in the NIST/ISMRM phantom demonstrated an average difference of 0.001% ± 015%, 1.2% ± 7.1%, and 0.47% ± 3% between two scans from the same 3-T scanner (repeatability), and 0.002% ± 015%, 0.62% ± 10.5%, and 0.97% ± 14% between the scans acquired on two different 3-T scanners (reproducibility) for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with rmsCV less than 1%, 2%, and 1% for T1 , T2 , and T1ρ , respectively. T1ρ relaxation time in the mild knee OA patients was significantly higher (p < 0.05) than in healthy subjects. The proposed 3D-MRF sequence is fast, reproducible, robust to B1 + inhomogeneity, and can simultaneously measure the T1 , T2 , T1ρ , and B1 + volumetric maps of the knee joint in a single scan within a clinically feasible scan time.

Choroid plexus tissue perfusion and blood to CSF barrier function in rats measured with continuous arterial spin labeling

The choroid plexus (ChP) of the cerebral ventricles is a source of cerebrospinal fluid (CSF) production and also plays a key role in immune surveillance at the level of blood-to-CSF-barrier (BCSFB). In this study, we quantify ChP blood perfusion and BCSFB mediated water exchange from arterial blood into ventricular CSF using non-invasive continuous arterial spin labelling magnetic resonance imaging (CASL-MRI). Systemic administration of anti-diuretic hormone (vasopressin) was used to validate BCSFB water flow as a metric of choroidal CSF secretory function. To further investigate the coupling between ChP blood perfusion and BCSFB water flow, we characterized the effects of two anesthetic regimens known to have large-scale differential effects on cerebral blood flow. For quantification of ChP blood perfusion a multi-compartment perfusion model was employed, and we discovered that partial volume correction improved measurement accuracy. Vasopressin significantly reduced both ChP blood perfusion and BCSFB water flow. ChP blood perfusion was significantly higher with pure isoflurane anesthesia (2-2.5%) when compared to a balanced anesthesia with dexmedetomidine and low-dose isoflurane (1.0 %), and significant correlation between ChP blood perfusion and BCSFB water flow was observed, however there was no significant difference in BCSFB water flow. In summary, here we introduce a non-invasive, robust, and spatially resolved in vivo imaging platform to quantify ChP blood perfusion as well as BCSFB water flow which can be applied to study coupling of these two key parameters in future clinical translational studies.

Acceleration Strategies for MR-STAT: Achieving High-Resolution Reconstructions on a Desktop PC Within 3 Minutes

MR-STAT is an emerging quantitative magnetic resonance imaging technique which aims at obtaining multi-parametric tissue parameter maps from single short scans. It describes the relationship between the spatial-domain tissue parameters and the time-domain measured signal by using a comprehensive, volumetric forward model. The MR-STAT reconstruction solves a large-scale nonlinear problem, thus is very computationally challenging. In previous work, MR-STAT reconstruction using Cartesian readout data was accelerated by approximating the Hessian matrix with sparse, banded blocks, and can be done on high performance CPU clusters with tens of minutes. In the current work, we propose an accelerated Cartesian MR-STAT algorithm incorporating two different strategies: firstly, a neural network is trained as a fast surrogate to learn the magnetization signal not only in the full time-domain but also in the compressed low-rank domain; secondly, based on the surrogate model, the Cartesian MR-STAT problem is re-formulated and split into smaller sub-problems by the alternating direction method of multipliers. The proposed method substantially reduces the computational requirements for runtime and memory. Simulated and in-vivo balanced MR-STAT experiments show similar reconstruction results using the proposed algorithm compared to the previous sparse Hessian method, and the reconstruction times are at least 40 times shorter. Incorporating sensitivity encoding and regularization terms is straightforward, and allows for better image quality with a negligible increase in reconstruction time. The proposed algorithm could reconstruct both balanced and gradient-spoiled in-vivo data within 3 minutes on a desktop PC, and could thereby facilitate the translation of MR-STAT in clinical settings.

Single-shot multi-parametric mapping based on multiple overlapping-echo detachment (MOLED) imaging

Multi-parametric quantitative magnetic resonance imaging (mqMRI) allows the characterization of multiple tissue properties non-invasively and has shown great potential to enhance the sensitivity of MRI measurements. However, real-time mqMRI during dynamic physiological processes or general motions remains challenging. To overcome this bottleneck, we propose a novel mqMRI technique based on multiple overlapping-echo detachment (MOLED) imaging, termed MQMOLED, to enable mqMRI in a single shot. In the data acquisition of MQMOLED, multiple MR echo signals with different multi-parametric weightings and phase modulations are generated and acquired in the same k-space. The k-space data is Fourier transformed and fed into a well-trained neural network for the reconstruction of multi-parametric maps. We demonstrated the accuracy and repeatability of MQMOLED in simultaneous mapping apparent proton density (APD) and any two parameters among T₂, T₂*, and apparent diffusion coefficient (ADC) in 130-170 ms. The abundant information delivered by the multiple overlapping-echo signals in MQMOLED makes the technique potentially robust to system imperfections, such as inhomogeneity of static magnetic field or radiofrequency field. Benefitting from the single-shot feature, MQMOLED exhibits a strong motion tolerance to the continuous movements of subjects. For the first time, it captured the synchronous changes of ADC, T₂, and T₁-weighted APD in contrast-enhanced perfusion imaging on patients with brain tumors, providing additional information about vascular density to the hemodynamic parametric maps. We expect that MQMOLED would promote the development of mqMRI technology and greatly benefit the applications of mqMRI, including therapeutics and analysis of metabolic/functional processes.

Artificial intelligence in cardiac magnetic resonance fingerprinting

Magnetic resonance fingerprinting (MRF) is a fast MRI-based technique that allows for multiparametric quantitative characterization of the tissues of interest in a single acquisition. In particular, it has gained attention in the field of cardiac imaging due to its ability to provide simultaneous and co-registered myocardial T₁ and T₂ mapping in a single breath-held cardiac MRF scan, in addition to other parameters. Initial results in small healthy subject groups and clinical studies have demonstrated the feasibility and potential of MRF imaging. Ongoing research is being conducted to improve the accuracy, efficiency, and robustness of cardiac MRF. However, these improvements usually increase the complexity of image reconstruction and dictionary generation and introduce the need for sequence optimization. Each of these steps increase the computational demand and processing time of MRF. The latest advances in artificial intelligence (AI), including progress in deep learning and the development of neural networks for MRI, now present an opportunity to efficiently address these issues. Artificial intelligence can be used to optimize candidate sequences and reduce the memory demand and computational time required for reconstruction and post-processing. Recently, proposed machine learning-based approaches have been shown to reduce dictionary generation and reconstruction times by several orders of magnitude. Such applications of AI should help to remove these bottlenecks and speed up cardiac MRF, improving its practical utility and allowing for its potential inclusion in clinical routine. This review aims to summarize the latest developments in artificial intelligence applied to cardiac MRF. Particularly, we focus on the application of machine learning at different steps of the MRF process, such as sequence optimization, dictionary generation and image reconstruction.

Accuracy and repeatability of joint sparsity multi-component estimation in MR Fingerprinting

MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estimation for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmentations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL- and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray matter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter.

Cerebrovascular activity is a major factor in the cerebrospinal fluid flow dynamics

Cerebrospinal fluid (CSF) provides physical protection to the central nervous system as well as an essential homeostatic environment for the normal functioning of neurons. Additionally, it has been proposed that the pulsatile movement of CSF may assist in glymphatic clearance of brain metabolic waste products implicated in neurodegeneration. In awake humans, CSF flow dynamics are thought to be driven primarily by cerebral blood volume fluctuations resulting from a number of mechanisms, including a passive vascular response to blood pressure variations associated with cardiac and respiratory cycles. Recent research has shown that mechanisms that rely on the action of vascular smooth muscle cells ("cerebrovascular activity") such as neuronal activity, changes in intravascular CO2, and autonomic activation from the brainstem, may lead to CSF pulsations as well. Nevertheless, the relative contribution of these mechanisms to CSF flow remains unclear. To investigate this further, we developed an MRI approach capable of disentangling and quantifying CSF flow components of different time scales associated with these mechanisms. This approach was evaluated on human control subjects (n = 12) performing intermittent voluntary deep inspirations, by determining peak flow velocities and displaced volumes between these mechanisms in the fourth ventricle. We found that peak flow velocities were similar between the different mechanisms, while displaced volumes per cycle were about a magnitude larger for deep inspirations. CSF flow velocity peaked at around 10.4 s (range 7.1-14.8 s, n = 12) following deep inspiration, consistent with known cerebrovascular activation delays for this autonomic challenge. These findings point to an important role of cerebrovascular activity in the genesis of CSF pulsations. Other regulatory triggers for cerebral blood flow such as autonomic arousal and orthostatic challenges may create major CSF pulsatile movement as well. Future quantitative comparison of these and possibly additional types of CSF pulsations with the proposed approach may help clarify the conditions that affect CSF flow dynamics.

Repeatability of simultaneous 3D 1H MRF/23Na MRI in brain at 7 T

Proton MRI can provide detailed morphological images, but it reveals little information about cell homeostasis. On the other hand, sodium MRI can provide metabolic information but cannot resolve fine structures. The complementary nature of proton and sodium MRI raises the prospect of their combined use in a single experiment. In this work, we assessed the repeatability of normalized proton density (PD), T₁, T₂, and normalized sodium density-weighted quantification measured with simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T, from ten healthy volunteers who were scanned three times each. The coefficients of variation (CV) and the intra-class correlation (ICC) were calculated for the mean and standard deviation (SD) of these 4 parameters in grey matter, white matter, and cerebrospinal fluid. As result, the CVs were lower than 3.3% for the mean values and lower than 6.9% for the SD values. The ICCs were higher than 0.61 in all 24 measurements. We conclude that the measurements of normalized PD, T₁, T₂, and normalized sodium density-weighted from simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T showed high repeatability. We estimate that changes > 6.6% (> 2 CVs) in mean values of both ¹H and ²³Na metrics could be detectable with this method.

Improved myocardial scar visualization with fast free-breathing motion-compensated black-blood T1-rho-prepared late gadolinium enhancement MRI

PURPOSE

Clinical guidelines recommend the use of bright-blood late gadolinium enhancement (BR-LGE) for the detection and quantification of regional myocardial fibrosis and scar. This technique, however, may suffer from poor contrast at the blood-scar interface, particularly in patients with subendocardial myocardial infarction. The purpose of this study was to assess the clinical performance of a two-dimensional black-blood LGE (BL-LGE) sequence, which combines free-breathing T₁-rho-prepared single-shot acquisitions with an advanced non-rigid motion-compensated patch-based reconstruction.

MATERIALS AND METHODS

Extended phase graph simulations and phantom experiments were performed to investigate the performance of the motion-correction algorithm and to assess the black-blood properties of the proposed sequence. Fifty-one patients (37 men, 14 women; mean age, 55 ± 15 [SD] years; age range: 19-81 years) with known or suspected cardiac disease prospectively underwent free-breathing T₁-rho-prepared BL-LGE imaging with inline non-rigid motion-compensated patch-based reconstruction at 1.5T. Conventional breath-held BR-LGE images were acquired for comparison purposes. Acquisition times were recorded. Two readers graded the image quality and relative contrasts were calculated. Presence, location, and extent of LGE were evaluated.

RESULTS

BL-LGE images were acquired with full ventricular coverage in 115 ± 25 (SD) sec (range: 64-160 sec). Image quality was significantly higher on free-breathing BL-LGE imaging than on its breath-held BR-LGE counterpart (3.6 ± 0.7 [SD] [range: 2-4] vs. 3.9 ± 0.2 [SD] [range: 3-4]) (P <0.01) and was graded as diagnostic for 44/51 (86%) patients. The mean scar-to-myocardium and scar-to-blood relative contrasts were significantly higher on BL-LGE images (P < 0.01 for both). The extent of LGE was larger on BL-LGE (median, 5 segments [IQR: 2, 7 segments] vs. median, 4 segments [IQR: 1, 6 segments]) (P < 0.01), the method being particularly sensitive in segments with LGE involving the subendocardium or papillary muscles. In eight patients (16%), BL-LGE could ascertain or rule out a diagnosis otherwise inconclusive on BR-LGE.

CONCLUSION

Free-breathing T₁-rho-prepared BL-LGE imaging with inline motion compensated reconstruction offers a promising diagnostic technology for the non-invasive assessment of myocardial injuries.

Data-driven myelin water imaging based on T1 and T2 relaxometry

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T₁ and T₂ relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T₁ /T₂ data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Feasibility of simultaneous high-resolution anatomical and quantitative magnetic resonance imaging of sciatic nerves in patients with Charcot-Marie-Tooth type 1A (CMT1A) at 7T

INTRODUCTION/AIMS

Magnetic resonance imaging (MRI) of peripheral nerves can provide image-based anatomical information and quantitative measurement. The aim of this pilot study was to investigate the feasibility of high-resolution anatomical and quantitative MRI assessment of sciatic nerve fascicles in patients with Charcot-Marie-Tooth (CMT) 1A using 7T field strength.

METHODS

Six patients with CMT1A underwent imaging on a high-gradient 7T MRI scanner using a 28-channel knee coil. Two high-resolution axial images were simultaneously acquired using a quantitative double-echo in steady-state (DESS) sequence. By comparing the two DESS echoes, T2 and apparent diffusion coefficient (ADC) maps were calculated. The cross-sectional areas and mean T2 and ADC were measured in individual fascicles of the tibial and fibular (peroneal) portions of the sciatic nerve at its bifurcation and 10 mm distally. Disease severity was measured using Charcot-Marie-Tooth Examination Score (CMTES) version 2 and compared to imaging findings.

RESULTS

We demonstrated the feasibility of 7T MRI of the proximal sciatic nerve in patients with CMT1A. Using the higher field, it was possible to measure individual bundles in the tibial and fibular divisions of the sciatic nerve. There was no apparent correlation between diffusion measures and disease severity in this small cohort.

DISCUSSION

This pilot study indicated that high-resolution MRI that allows for combined anatomical and quantitative imaging in one scan is feasible at 7T field strengths and can be used to investigate the microstructure of individual nerve fascicles.

A Fetal Brain magnetic resonance Acquisition Numerical phantom (FaBiAN)

Accurate characterization of in utero human brain maturation is critical as it involves complex and interconnected structural and functional processes that may influence health later in life. Magnetic resonance imaging is a powerful tool to investigate equivocal neurological patterns during fetal development. However, the number of acquisitions of satisfactory quality available in this cohort of sensitive subjects remains scarce, thus hindering the validation of advanced image processing techniques. Numerical phantoms can mitigate these limitations by providing a controlled environment with a known ground truth. In this work, we present FaBiAN, an open-source Fetal Brain magnetic resonance Acquisition Numerical phantom that simulates clinical T2-weighted fast spin echo sequences of the fetal brain. This unique tool is based on a general, flexible and realistic setup that includes stochastic fetal movements, thus providing images of the fetal brain throughout maturation comparable to clinical acquisitions. We demonstrate its value to evaluate the robustness and optimize the accuracy of an algorithm for super-resolution fetal brain magnetic resonance imaging from simulated motion-corrupted 2D low-resolution series compared to a synthetic high-resolution reference volume. We also show that the images generated can complement clinical datasets to support data-intensive deep learning methods for fetal brain tissue segmentation.

Myelin water fraction mapping from multiple echo spin echoes and an independent B1 + map

Purpose

Myelin water fraction (MWF) is often obtained from a multiple echo spin echo (MESE) sequence using multi‐component T₂ fitting with non‐negative least squares. This process fits many unknowns including B₁⁺ to produce a T₂ spectrum for each voxel. Presented is an alternative using a rapid B₁⁺ mapping sequence to supply B₁⁺ for the MWF fitting procedure.

Methods

Effects of B₁⁺ errors on MWF calculations were modeled for 2D and 3D MESE using Bloch and extended phase graph simulations, respectively. Variations in SNR and relative refocusing widths were tested. Human brain experiments at 3 T used 2D MESE and an independent B₁⁺ map. MWF maps were produced with the standard approach and with the use of the independent B₁⁺ map. Differences in B₁⁺ and mean MWF in specific brain regions were compared.

Results

For 2D MESE, MWF with the standard method was strongly affected by B₁⁺ misestimations arising from limited SNR and response asymmetry around 180°, which decreased with increasing relative refocusing width. Using an independent B₁⁺ map increased mean MWF and decreased coefficient of variation. Notable differences in vivo in 2D MESE were in areas of high B₁⁺ such as thalamus and splenium where mean MWF increased by 88% and 31%, respectively (P < 0.001). Simulations also demonstrated the advantages of this approach for 3D MESE when SNR is <500.

Conclusion

For 2D MESE, because of increased complexity of decay curves and limited SNR, supplying B₁⁺ improves MWF results in peripheral and central brain regions where flip angles differ substantially from 180°.

Data‐driven motion‐corrected brain MRI incorporating pose‐dependent B 0 fields

Purpose

To develop a fully data‐driven retrospective intrascan motion‐correction framework for volumetric brain MRI at ultrahigh field (7 Tesla) that includes modeling of pose‐dependent changes in polarizing magnetic (B₀) fields.

Theory and Methods

Tissue susceptibility induces spatially varying B₀ distributions in the head, which change with pose. A physics‐inspired B₀ model has been deployed to model the B₀ variations in the head and was validated in vivo. This model is integrated into a forward parallel imaging model for imaging in the presence of motion. Our proposal minimizes the number of added parameters, enabling the developed framework to estimate dynamic B₀ variations from appropriately acquired data without requiring navigators. The effect on data‐driven motion correction is validated in simulations and in vivo.

Results

The applicability of the physics‐inspired B₀ model was confirmed in vivo. Simulations show the need to include the pose‐dependent B₀ fields in the reconstruction to improve motion‐correction performance and the feasibility of estimating B₀ evolution from the acquired data. The proposed motion and B₀ correction showed improved image quality for strongly corrupted data at 7 Tesla in simulations and in vivo.

Conclusion

We have developed a motion‐correction framework that accounts for and estimates pose‐dependent B₀ fields. The method improves current state‐of‐the‐art data‐driven motion‐correction techniques when B₀ dependencies cannot be neglected. The use of a compact physics‐inspired B₀ model together with leveraging the parallel imaging encoding redundancy and previously proposed optimized sampling patterns enables a purely data‐driven approach.

Simultaneous bilateral T1 , T2 , and T1ρ relaxation mapping of the hip joint with magnetic resonance fingerprinting

Quantitative MRI can detect early biochemical changes in cartilage, but its bilateral use in clinical routines is challenging. The aim of this prospective study was to demonstrate the feasibility of magnetic resonance fingerprinting for bilateral simultaneous T₁ , T₂ , and T_1ρ mapping of the hip joint. The study population consisted of six healthy volunteers with no known trauma or pain in the hip. Monoexponential T₁ , T₂ , and T_1ρ relaxation components were assessed in femoral lateral, superolateral, and superomedial, and inferior, as well as acetabular, superolateral, and superomedial subregions in left and right hip cartilage. Aligned ranked nonparametric factorial analysis was used to assess the side's impact on the subregions. Kruskal-Wallis and Wilcoxon tests were used to compare subregions, and coefficient of variation to assess repeatability. Global averages of T₁ (676.0 ± 45.4 and 687.6 ± 44.5 ms), T₂ (22.5 ± 2.6 and 22.1 ± 2.5 ms), and T_1ρ (38.2 ± 5.5 and 38.2 ± 5.5 ms) were measured in the left and right hip, and articular cartilage, respectively. The Kruskal-Wallis test showed a significant difference between different subregions' relaxation times regardless of the hip side (p < 0.001 for T₁ , p = 0.012 for T₂ , and p < 0.001 for T_1ρ ). The Wilcoxon test showed that T₁ of femoral layers was significantly (p < 0.003) higher than that for acetabular cartilage. The experiments showed excellent repeatability with CV_rms of 1%, 2%, and 4% for T₁ , T₂ , and T_1ρ, respectively. It was concluded that bilateral T₁ , T₂ , and T_1ρ relaxation times, as well as B₁⁺ maps, can be acquired simultaneously from hip joints using the proposed MRF sequence.

Simultaneous 3D acquisition of 1 H MRF and 23 Na MRI

PURPOSE

To develop a 3D MR technique to simultaneously acquire proton multiparametric maps (T₁ , T₂ , and proton density) and sodium density weighted images over the whole brain.

METHODS

We implemented a 3D stack-of-stars MR pulse sequence which consists of interleaved proton (¹ H) and sodium (²³ Na) excitations, tailored slice encoding gradients that can encode the same slice for both nuclei, and simultaneous readout with different radial trajectories (¹ H, full-radial; ²³ Na, center-out radial). The receive chain of our 7T scanner was modified to enable simultaneous acquisition of ¹ H and ²³ Na signal. A heuristically optimized flip angle train was implemented for proton MR fingerprinting (MRF). The SNR and the accuracy of proton T₁ and T₂ were evaluated in phantoms. Finally, in vivo application of the method was demonstrated in five healthy subjects.

RESULTS

The SNR for the simultaneous measurement was almost identical to that for the single-nucleus measurements (<2% change). The proton T₁ and T₂ maps remained similar to the results from a reference 2D MRF technique (normalized RMS error in T₁ ≈ 4.2% and T₂ ≈ 11.3%). Measurements in healthy subjects corroborated these results and demonstrated the feasibility of our method for in vivo application. The in vivo T₁ values measured using our method were lower than the results measured by other conventional techniques.

CONCLUSIONS

With the 3D simultaneous implementation, we were able to acquire sodium and proton density weighted images in addition to proton T₁ , T₂ , and B1+ from ¹ H MRF that covers the whole brain volume within 21 min.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T₁ ), the spin-spin relaxation time (T₂ ), and the spin-lattice relaxation in the rotating frame (T_1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T₁ -, T₂ -, or T_1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

The Digital Brain Bank, an open access platform for post-mortem imaging datasets

Post-mortem magnetic resonance imaging (MRI) provides the opportunity to acquire high-resolution datasets to investigate neuroanatomy and validate the origins of image contrast through microscopy comparisons. We introduce the Digital Brain Bank (open.win.ox.ac.uk/DigitalBrainBank), a data release platform providing open access to curated, multimodal post-mortem neuroimaging datasets. Datasets span three themes-Digital Neuroanatomist: datasets for detailed neuroanatomical investigations; Digital Brain Zoo: datasets for comparative neuroanatomy; and Digital Pathologist: datasets for neuropathology investigations. The first Digital Brain Bank data release includes 21 distinctive whole-brain diffusion MRI datasets for structural connectivity investigations, alongside microscopy and complementary MRI modalities. This includes one of the highest-resolution whole-brain human diffusion MRI datasets ever acquired, whole-brain diffusion MRI in fourteen nonhuman primate species, and one of the largest post-mortem whole-brain cohort imaging studies in neurodegeneration. The Digital Brain Bank is the culmination of our lab's investment into post-mortem MRI methodology and MRI-microscopy analysis techniques. This manuscript provides a detailed overview of our work with post-mortem imaging to date, including the development of diffusion MRI methods to image large post-mortem samples, including whole, human brains. Taken together, the Digital Brain Bank provides cross-scale, cross-species datasets facilitating the incorporation of post-mortem data into neuroimaging studies.

Simultaneous relaxometry and morphometry of human brain structures with 3D magnetic resonance fingerprinting: a multicenter, multiplatform, multifield-strength study

Relaxation times and morphological information are fundamental magnetic resonance imaging-derived metrics of the human brain that reflect the status of the underlying tissue. Magnetic resonance fingerprinting (MRF) enables simultaneous acquisition of T1 and T2 maps inherently aligned to the anatomy, allowing whole-brain relaxometry and morphometry in a single scan. In this study, we revealed the feasibility of 3D MRF for simultaneous brain structure-wise morphometry and relaxometry. Comprehensive test-retest scan analyses using five 1.5-T and three 3.0-T systems from a single vendor including different scanner types across 3 institutions demonstrated that 3D MRF-derived morphological information and relaxation times are highly repeatable at both 1.5 T and 3.0 T. Regional cortical thickness and subcortical volume values showed high agreement and low bias across different field strengths. The ability to acquire a set of regional T1, T2, thickness, and volume measurements of neuroanatomical structures with high repeatability and reproducibility facilitates the ability of longitudinal multicenter imaging studies to quantitatively monitor changes associated with underlying pathologies, disease progression, and treatments.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging

PURPOSE

To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF).

METHODS

Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B₀ inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect.

RESULTS

The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T₁ and T₂ mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures.

CONCLUSIONS

The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T₁, T₂, T₂*, PD and B₁⁺ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

T1-independent exchange rate quantification using saturation- or phase sensitive-water exchange spectroscopy

PURPOSE

Water Exchange Spectroscopy (WEX) is a direct measurement of the exchange rate k_sw of labile protons from a solute to water in which the exchange time is varied. However the useful information can be masked by the T₁-decay of the solvent pool. We propose Saturation-WEX and Phase Sensitive WEX (PS-WEX) as an extension upon the WEX approach to reduce T₁-masking. Additionally PS-WEX takes advantage of the phase information contained in the WEX signal to improve the dynamic range.

METHODS

By introducing an additional RF-pulse and fixing the exchange time delay the T₁-dependence of the signal is reduced. By exploiting the phase sensitivity of the WEX pathway the dynamic range can be increased. This approach is validated using simulations as well as phantom measurements.

RESULTS

The improved dynamic range is demonstrated in measurements. The fixed exchange time reduces the influence of the T₁-decay on the signal curve leading to improved fit quality.

CONCLUSION

Sat-WEX and PS-WEX are an extension to the well established WEX approach with a less complex fit equation and in the case of PS-WEX improved dynamic range, allowing more accurate quantification.

Myocardial T1, T2, T2*, and fat fraction quantification via low-rank motion-corrected cardiac MR fingerprinting

Purpose

Develop a novel 2D cardiac MR fingerprinting (MRF) approach to enable simultaneous T1, T2, T2*, and fat fraction (FF) myocardial tissue characterization in a single breath‐hold scan.

Methods

Simultaneous, co‐registered, multi‐parametric mapping of T1, T2, and FF has been recently achieved with cardiac MRF. Here, we further incorporate T2* quantification within this approach, enabling simultaneous T1, T2, T2*, and FF myocardial tissue characterization in a single breath‐hold scan. T2* quantification is achieved with an eight‐echo readout that requires a long cardiac acquisition window. A novel low‐rank motion‐corrected (LRMC) reconstruction is exploited to correct for cardiac motion within the long acquisition window. The proposed T1/T2/T2*/FF cardiac MRF was evaluated in phantom and in 10 healthy subjects in comparison to conventional mapping techniques.

Results

The proposed approach achieved high quality parametric mapping of T1, T2, T2*, and FF with corresponding normalized RMS error (RMSE) T1 = 5.9%, T2 = 9.6% (T2 values <100 ms), T2* = 3.3% (T2* values <100 ms), and FF = 0.8% observed in phantom scans. In vivo, the proposed approach produced higher left‐ventricular myocardial T1 values than MOLLI (1148 vs 1056 ms), lower T2 values than T2‐GraSE (42.8 vs 50.6 ms), lower T2* values than eight‐echo gradient echo (GRE) (35.0 vs 39.4 ms), and higher FF values than six‐echo GRE (0.8 vs 0.3 %) reference techniques. The proposed approach achieved considerable reduction in motion artifacts compared to cardiac MRF without motion correction, improved spatial uniformity, and statistically higher apparent precision relative to conventional mapping for all parameters.

Conclusion

The proposed cardiac MRF approach enables simultaneous, co‐registered mapping of T1, T2, T2*, and FF in a single breath‐hold for comprehensive myocardial tissue characterization, achieving higher apparent precision than conventional methods.

Comparison of T2 Quantification Strategies in the Abdominal-Pelvic Region for Clinical Use

OBJECTIVES

The aim of the study was to compare different magnetic resonance imaging (MRI) acquisition strategies appropriate for T2 quantification in the abdominal-pelvic area. The different techniques targeted in the study were chosen according to 2 main considerations: performing T2 measurement in an acceptable time for clinical use and preventing/correcting respiratory motion.

MATERIALS AND METHODS

Acquisitions were performed at 3 T. To select sequences for in vivo measurements, a phantom experiment was conducted, for which the T2 values obtained with the different techniques of interest were compared with the criterion standard (single-echo SE sequence, multiple acquisitions with varying echo time). Repeatability and temporal reproducibility studies for the different techniques were also conducted on the phantom. Finally, an in vivo study was conducted on 12 volunteers to compare the techniques that offer acceptable acquisition time for clinical use and either address or correct respiratory motion.

RESULTS

For the phantom study, the DESS and T2-preparation techniques presented the lowest precision (ρ2 = 0.9504 and ρ2 = 0.9849 respectively), and showed a poor repeatability/reproducibility compared with the other techniques. The strategy relying on SE-EPI showed the best precision and accuracy (ρ2 = 0.9994 and Cb = 0.9995). GRAPPATINI exhibited a very good precision (ρ2 = 0.9984). For the technique relying on radial TSE, the precision was not as good as GRAPPATINI (ρ2 = 0.9872). The in vivo study demonstrated good respiratory motion management for all of the selected techniques. It also showed that T2 estimate ranges were different from one method to another. For GRAPPATINI and radial TSE techniques, there were significant differences between all the different types of organs of interest.

CONCLUSIONS

To perform T2 measurement in the abdominal-pelvic region, one should favor a technique with acceptable acquisition time for clinical use, with proper respiratory motion management, with good repeatability, reproducibility, and precision. In this study, the techniques relying respectively on SE-EPI, radial TSE, and GRAPPATINI appeared as good candidates.

A method to remove the influence of fixative concentration on postmortem T2 maps using a kinetic tensor model

Formalin fixation has been shown to substantially reduce T₂ estimates, primarily driven by the presence of fixative in tissue. Prior to scanning, post‐mortem samples are often placed into a fluid that has more favourable imaging properties. This study investigates whether there is evidence for a change in T₂ in regions close to the tissue surface due to fixative outflux into this surrounding fluid. Furthermore, we investigate whether a simulated spatial map of fixative concentration can be used as a confound regressor to reduce T₂ inhomogeneity. To achieve this, T₂ maps and diffusion tensor estimates were obtained in 14 whole, formalin‐fixed post‐mortem brains placed in Fluorinert approximately 48 hr prior to scanning. Seven brains were fixed with 10% formalin and seven brains were fixed with 10% neutral buffered formalin (NBF). Fixative outflux was modelled using a proposed kinetic tensor (KT) model, which incorporates voxelwise diffusion tensor estimates to account for diffusion anisotropy and tissue‐specific diffusion coefficients. Brains fixed with 10% NBF revealed a spatial T₂ pattern consistent with modelled fixative outflux. Confound regression of fixative concentration reduced T₂ inhomogeneity across both white and grey matter, with the greatest reduction attributed to the KT model versus simpler models of fixative outflux. No such effect was observed in brains fixed with 10% formalin. Correlations between the transverse relaxation rate R₂ and ferritin/myelin proteolipid protein (PLP) histology lead to an increased similarity for the relationship between R₂ and PLP for the two fixative types after KT correction.

Laminar perfusion imaging with zoomed arterial spin labeling at 7 Tesla

Laminar fMRI based on BOLD and CBV contrast at ultrahigh magnetic fields has been applied for studying the dynamics of mesoscopic brain networks. However, the quantitative interpretations of BOLD/CBV fMRI results are confounded by different baseline physiology across cortical layers. Here we introduce a novel 3D zoomed pseudo-continuous arterial spin labeling (pCASL) technique at 7T that offers the capability for quantitative measurements of laminar cerebral blood flow (CBF) both at rest and during task activation with high spatial specificity and sensitivity. We found arterial transit time in superficial layers is ∼100 ms shorter than in middle/deep layers revealing the time course of labeled blood flowing from pial arteries to downstream microvasculature. Resting state CBF peaked in the middle layers which is highly consistent with microvascular density measured from human cortex specimens. Finger tapping induced a robust two-peak laminar profile of CBF increases in the superficial (somatosensory and premotor input) and deep (spinal output) layers of M1, while finger brushing task induced a weaker CBF increase in superficial layers (somatosensory input). This observation is highly consistent with reported laminar profiles of CBV activation on M1. We further demonstrated that visuospatial attention induced a predominant CBF increase in deep layers and a smaller CBF increase on top of the lower baseline CBF in superficial layers of V1 (feedback cortical input), while stimulus driven activity peaked in the middle layers (feedforward thalamic input). With the capability for quantitative CBF measurements both at baseline and during task activation, high-resolution ASL perfusion fMRI at 7T provides an important tool for in vivo assessment of neurovascular function and metabolic activities of neural circuits across cortical layers.

Ultra-long-TE arterial spin labeling reveals rapid and brain-wide blood-to-CSF water transport in humans

The study of brain clearance mechanisms is an active area of research. While we know that the cerebrospinal fluid (CSF) plays a central role in one of the main existing clearance pathways, the exact processes for the secretion of CSF and the removal of waste products from tissue are under debate. CSF is thought to be created by the exchange of water and ions from the blood, which is believed to mainly occur in the choroid plexus. This exchange has not been thoroughly studied in vivo.

We propose a modified arterial spin labeling (ASL) MRI sequence and image analysis to track blood water as it is transported to the CSF, and to characterize its exchange from blood to CSF. We acquired six pseudo-continuous ASL sequences with varying labeling duration (LD) and post-labeling delay (PLD) and a segmented 3D-GRASE readout with a long echo train (8 echo times (TE)) which allowed separation of the very long-T₂ CSF signal. ASL signal was observed at long TEs (793 ms and higher), indicating presence of labeled water transported from blood to CSF. This signal appeared both in the CSF proximal to the choroid plexus and in the subarachnoid space surrounding the cortex. ASL signal was separated into its blood, gray matter and CSF components by fitting a triexponential function with T₂s taken from literature. A two-compartment dynamic model was introduced to describe the exchange of water through time and TE. From this, a water exchange time from the blood to the CSF (T_bl->CSF) was mapped, with an order of magnitude of approximately 60 s.

T2 quantification in brain using 3D fast spin-echo imaging with long echo trains

Purpose

Three‐dimensional fast spin‐echo (FSE) sequences commonly use very long echo trains (>64 echoes) and severely reduced refocusing angles. They are increasingly used in brain exams due to high, isotropic resolution and reasonable scan time when using long trains and short interecho spacing. In this study, T₂ quantification in 3D FSE is investigated to achieve increased resolution when comparing with established 2D (proton‐density dual‐echo and multi‐echo spin‐echo) methods.

Methods

The FSE sequence design was explored to use long echo trains while minimizing T₂ fitting error and maintaining typical proton density and T₂‐weighted contrasts. Constant and variable flip angle trains were investigated using extended phase graph and Bloch equation simulations. Optimized parameters were analyzed in phantom experiments and validated in vivo in comparison to 2D methods for eight regions of interest in brain, including deep gray‐matter structures and white‐matter tracts.

Results

Phantom and healthy in vivo brain T₂ measurements showed that optimized variable echo‐train 3D FSE performs similarly to previous 2D methods, while achieving three‐fold‐higher slice resolution, evident visually in the 3D T₂ maps. Optimization resulted in better T₂ fitting and compared well with standard multi‐echo spin echo (within the 8‐ms confidence limits defined based on Bland‐Altman analysis).

Conclusion

T₂ mapping using 3D FSE with long echo trains and variable refocusing angles provides T₂ accuracy in agreement with 2D methods with additional high‐resolution benefits, allowing isotropic views while avoiding incidental magnetization transfer effects. Consequently, optimized 3D sequences should be considered when choosing T₂ mapping methods for high anatomic detail.

Multi-exponential MRI T2 maps: A tool to classify and characterize fruit tissues

The estimation of multi-exponential relaxation time T₂ and their associated amplitudes A₀ at the voxel level has been made possible by recent developments in the field of image processing. These data are of great interest for the characterization of biological tissues, such as fruit tissues. However, they represent a high number of information, not easily interpretable. Moreover, the non-uniformity of the MRI images, which mainly directly impacts A₀, could induce interpretation errors. In this paper, we propose a post-processing scheme that clusters similar voxels according to the multi-exponential relaxation parameters in order to reduce the complexity of the information while avoiding the problems associated with intensity non-uniformity. We also suggest a data representation suitable for the visualization of the multi-T₂ distribution within each tissue. We illustrate this work with results for different fruits, demonstrating the great potential of multi-T₂ information to shed new light on fruit characterization.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares

Multi-echo T₂ magnetic resonance images contain information about the distribution of T₂ relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T₂ spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X²) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X², and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T₂ spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T₂ distributions and MWF maps.

A Spatial Off-Resonance Correction in Spirals for Magnetic Resonance Fingerprinting

In MR Fingerprinting (MRF), balanced Steady-State Free Precession (bSSFP) has advantages over unbalanced SSFP because it retains the spin history achieving a higher signal-to-noise ratio (SNR) and scan efficiency. However, bSSFP-MRF is not frequently used because it is sensitive to off-resonance, producing artifacts and blurring, and affecting the parametric map quality. Here we propose a novel Spatial Off-resonance Correction (SOC) approach for reducing these artifacts in bSSFP-MRF with spiral trajectories. SOC-MRF uses each pixel's Point Spread Function to create system matrices that encode both off-resonance and gridding effects. We iteratively compute the inverse of these matrices to reduce the artifacts. We evaluated the proposed method using brain simulations and actual MRF acquisitions of a standardized T1/T2 phantom and five healthy subjects. The results show that the off-resonance distortions in T1/T2 maps were considerably reduced using SOC-MRF. For T2, the Normalized Root Mean Square Error (NRMSE) was reduced from 17.3 to 8.3% (simulations) and from 35.1 to 14.9% (phantom). For T1, the NRMS was reduced from 14.7 to 7.7% (simulations) and from 17.7 to 6.7% (phantom). For in-vivo, the mean and standard deviation in different ROI in white and gray matter were significantly improved. For example, SOC-MRF estimated an average T2 for white matter of 77ms (the ground truth was 74ms) versus 50 ms of MRF. For the same example the standard deviation was reduced from 18 ms to 6ms. The corrections achieved with the proposed SOC-MRF may expand the potential applications of bSSFP-MRF, taking advantage of its better SNR property.

A New Method for Cartilage Evaluation in Femoroacetabular Impingement Using Quantitative T2 Magnetic Resonance Imaging: Preliminary Validation against Arthroscopic Findings

OBJECTIVE

The outcome of arthroscopic treatment for femoroacetabular impingement (FAI) depends on the preoperative status of the hip cartilage. Quantitative T2 can detect early biochemical cartilage changes, but its routine implementation is challenging. Furthermore, intrinsic T2 variability between patients makes it difficult to define a threshold to identify cartilage lesions. To address this, we propose a normalized T2-index as a new method to evaluate cartilage in FAI.

DESIGN

We retrospectively analyzed magnetic resonance imaging (MRI) data of 18 FAI patients with arthroscopically confirmed cartilage defects. Cartilage T2 maps were reconstructed from multi-spin-echo 3-T data using the echo-modulation-curve (EMC) model-based technique. The central femoral cartilage, assumed healthy in early-stage FAI, was used as the normalization reference to define a T2-index. We investigated the ability of the T2-index to detect surgically confirmed cartilage lesions.

RESULTS

The average T2-index was 1.14 ± 0.1 and 1.13 ± 0.1 for 2 separated segmentations. Using T2-index >1 as the threshold for damaged cartilage, accuracy was 88% and 100% for the 2 segmentations. We found moderate intraobserver repeatability, although separate segmentations yielded comparable accuracy. Damaged cartilage could not be identified using nonnormalized average T2 values.

CONCLUSIONS

This preliminary study confirms the importance of normalizing T2 values to account for interpatient variability and suggests that the T2-index is a promising biomarker for the detection of cartilage lesions in FAI. Future work is needed to confirm that combining T2-index with morphologic MRI and other quantitative biomarkers could improve cartilage assessment in FAI.

B1 inhomogeneity correction of RARE MRI at low SNR: Quantitative in vivo 19 F MRI of mouse neuroinflammation with a cryogenically-cooled transceive surface radiofrequency probe

PURPOSE

Low SNR in fluorine-19 (¹⁹ F) MRI benefits from cryogenically-cooled transceive surface RF probes (CRPs), but strong B₁ inhomogeneities hinder quantification. Rapid acquisition with refocused echoes (RARE) is an SNR-efficient method for MRI of neuroinflammation with perfluorinated compounds but lacks an analytical signal intensity equation to retrospectively correct B₁ inhomogeneity. Here, a workflow was proposed and validated to correct and quantify ¹⁹ F-MR signals from the inflamed mouse brain using a ¹⁹ F-CRP.

METHODS

In vivo ¹⁹ F-MR images were acquired in a neuroinflammation mouse model with a quadrature ¹⁹ F-CRP using an imaging setup including 3D-printed components to acquire co-localized anatomical and ¹⁹ F images. Model-based corrections were validated on a uniform ¹⁹ F phantom and in the neuroinflammatory model. Corrected ¹⁹ F-MR images were benchmarked against reference images and overlaid on in vivo ¹ H-MR images. Computed concentration uncertainty maps using Monte Carlo simulations served as a measure of performance of the B₁ corrections.

RESULTS

Our study reports on the first quantitative in vivo ¹⁹ F-MR images of an inflamed mouse brain using a ¹⁹ F-CRP, including in vivo T₁ calculations for ¹⁹ F-nanoparticles during pathology and B₁ corrections for ¹⁹ F-signal quantification. Model-based corrections markedly improved ¹⁹ F-signal quantification from errors > 50% to < 10% in a uniform phantom (p < 0.001). Concentration uncertainty maps ex vivo and in vivo yielded uncertainties that were generally < 25%. Monte Carlo simulations prescribed SNR ≥ 10.1 to reduce uncertainties < 10%, and SNR ≥ 4.25 to achieve uncertainties < 25%.

CONCLUSION

Our model-based correction method facilitated ¹⁹ F signal quantification in the inflamed mouse brain when using the SNR-boosting ¹⁹ F-CRP technology, paving the way for future low-SNR ¹⁹ F-MRI applications in vivo.

Simultaneous T1 , T2 , and T1ρ cardiac magnetic resonance fingerprinting for contrast agent-free myocardial tissue characterization

PURPOSE

To develop a simultaneous T₁ , T₂ , and T_1ρ cardiac magnetic resonance fingerprinting (MRF) approach to enable comprehensive contrast agent-free myocardial tissue characterization in a single breath-hold scan.

METHODS

A 2D gradient-echo electrocardiogram-triggered cardiac MRF sequence with low flip angles, varying magnetization preparation, and spiral trajectory was acquired at 1.5 T to encode T₁ , T₂ , and T_1⍴ simultaneously. The MRF images were reconstructed using low-rank inversion, regularized with a multicontrast patch-based higher-order reconstruction. Parametric maps were generated and matched in the singular value domain to extended phase graph-based dictionaries. The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional methods in terms of coefficients of determination and best fits for the phantom study, and in terms of Bland-Altman agreement, average values and coefficient of variation of T₁ , T₂ , and T_1⍴ for the healthy subjects study.

RESULTS

The T₁ , T₂ , and T_1⍴ MRF values showed excellent correlation with conventional spin-echo and clinical mapping methods in phantom studies (r² > 0.97). Measured MRF values in myocardial tissue (mean ± SD) were 1133 ± 33 ms, 38.8 ± 3.5 ms, and 52.0 ± 4.0 ms for T₁ , T₂ and T1_⍴ , respectively, against 1053 ± 47 ms, 50.4 ± 3.9 ms, and 55.9 ± 3.3 ms for T₁ modified Look-Locker inversion imaging, T₂ gradient and spin echo, and T_1⍴ turbo field echo, respectively.

CONCLUSION

A cardiac MRF approach for simultaneous quantification of myocardial T₁ , T₂ , and T_1ρ in a single breath-hold MR scan of about 16 seconds has been proposed. The approach has been investigated in phantoms and healthy subjects showing good agreement with reference spin echo measurements and conventional clinical maps.

Simultaneous comprehensive liver T1 , T2 , T 2 ∗ , T1ρ , and fat fraction characterization with MR fingerprinting

PURPOSE

To develop a novel simultaneous co-registered T₁ , T₂ , T 2 ∗ , T_1ρ , and fat fraction abdominal MR fingerprinting (MRF) approach for fully comprehensive liver-tissue characterization in a single breath-hold scan.

METHODS

A gradient-echo liver MRF sequence with low fixed flip angle, multi-echo radial readout, and varying magnetization preparation pulses for multiparametric encoding is performed at 1.5 T. The T 2 ∗ and fat fraction are estimated from a graph/cut water/fat separation method using a six-peak fat model. Water/fat singular images obtained are then matched to an MRF dictionary, estimating water-specific T₁ , T₂ , and T_1ρ . The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional sequences.

RESULTS

For the phantom studies, linear fits show excellent coefficients of determination (r² > 0.9) for every parametric map. For in vivo studies, the average values measured within regions of interest drawn on liver, spleen, muscle, and fat are statistically different from the reference scans (p < 0.05) for T₁ , T₂ , and T_1⍴ but not for T 2 ∗ and fat fraction, whereas correlation between MRF and reference scans is excellent for each parameter (r² > 0.92 for every parameter).

CONCLUSION

The proposed multi-echo inversion-recovery, T₂ , and T_1⍴ prepared liver MRF sequence presented in this work allows for quantitative T₁ , T₂ , T 2 ∗ , T_1⍴ , and fat fraction liver-tissue characterization in a single breath-hold scan of 18 seconds. The approach showed good agreement and correlation with respect to reference clinical maps.

A strategy to reduce the sensitivity of inhomogeneous magnetization transfer (ihMT) imaging to radiofrequency transmit field variations at 3 T

PURPOSE

To minimize the sensitivity of inhomogeneous magnetization transfer gradient-echo (ihMT-GRE) imaging to radiofrequency (RF) transmit field ( B 1 + ) inhomogeneities at 3 T.

METHODS

The ihMT-GRE sequence was optimized by varying the concentration of the RF saturation energy over time, obtained by increasing the saturation pulse power while extending the sequence repetition time (TR). Different protocols were tested using numerical simulations and human in vivo experiments in the brain white matter (WM) of healthy subjects at 3 T. The sensitivity of the ihMT ratio (ihMTR) to B 1 + variations was investigated by comparing measurements obtained at nominal transmitter adjustments and following a 20% global B 1 + drop. The resulting relative variations (δ_ihMTR ) were evaluated voxelwise as a function of the local B 1 + distribution. The reproducibility of the protocol providing minimal B 1 + bias was assessed in a test-retest experiment.

RESULTS

In line with simulations, ihMT-GRE experiments conducted at high concentration of the RF energy over time demonstrated strong reduction of the B 1 + inhomogeneity effects in the human WM. Under the optimal conditions of 350-ms TR and 3-µT root mean square (RMS) saturation power, 73% of all WM voxels presented δ_ihMTR below 10%. Reproducibility analysis yielded a close-to-zero systematic bias (ΔihMTR = -0.081%) and a high correlation (ρ² = 0.977) between test and retest experiments.

CONCLUSION

Concentrating RF saturation energy in ihMT-GRE sequences mitigates the sensitivity of the ihMTR to B 1 + variations and allows for clinical-ready ihMT imaging at 3 T. This feature is of particular interest for high and ultra-high field applications.

Complex B1 + mapping with Carr-Purcell spin echoes and its application to electrical properties tomography

Purpose

To present a new complex‐valued B₁⁺ mapping method for electrical properties tomography using Carr‐Purcell spin echoes.

Methods

A Carr‐Purcell (CP) echo train generates pronounced flip‐angle dependent oscillations that can be used to estimate the magnitude of B₁⁺. To this end, a dictionary is used that takes into account the slice profile as well as T₂ relaxation along the echo train. For validation, the retrieved B₁⁺ map is compared with the actual flip angle imaging (AFI) method in a phantom (79 ε₀, 0.34 S/m). Moreover, the phase of the first echo reflects the transceive phase. Overall, the CP echo train yields an estimate of the complex‐valued B₁⁺, allowing electrical properties tomography with both permittivity and conductivity. The presented method is evaluated in phantom scans as well as for in vivo brain at 3 T.

Results

In the phantom, the obtained magnitude B₁⁺ maps retrieved from the CP echo train and the AFI method show excellent agreement, and both the reconstructed estimated permittivity (79 ± 3) ε₀ and conductivity (0.35 ± 0.04) S/m values are in accordance with expectations. In the brain, the obtained electrical properties are also close to expectations. In addition to the retrieved complex B₁⁺ information, the decay of the CP echo trains also yields an estimate for T₂.

Conclusion

The CP sequence can be used to simultaneously provide both B₁⁺ magnitude and phase estimations, and therefore allows for full reconstruction of the electrical properties.

A unifying view on extended phase graphs and Bloch simulations for quantitative MRI

Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.