Multi-compartment relaxometry and diffusion informed myelin water imaging – Promises and challenges of new gradient echo myelin water imaging methods

Quantitative MRI distinguishes different leukodystrophies and correlates with clinical measures

OBJECTIVES

The leukodystrophy "vanishing white matter" (VWM) and "metachromatic leukodystrophy" (MLD) affect the brain's white matter, but have very different underlying pathology. We aim to determine whether quantitative MRI reflects known neuropathological differences and correlates with clinical scores in these leukodystrophies.

METHODS

VWM and MLD patients and controls were prospectively included between 2020 and 2023. Clinical scores were recorded. MRI at 3 T included multi-compartment relaxometry diffusion-informed myelin water imaging (MCR-DIMWI) and multi-echo T2-relaxation imaging with compressed sensing (METRICS) to determine myelin water fractions (MWF). Multi-shell diffusion-weighted data were used for diffusion tensor imaging measures and neurite orientation dispersion and density imaging (NODDI) analysis, which estimates neurite density index, orientation dispersion index, and free water fraction. As quantitative MRI measures are age-dependent, ratios between actual and age-expected MRI measures were calculated. We performed the multilevel analysis with subsequent post-hoc and correlation tests to assess differences between groups and clinico-radiological correlations.

RESULTS

Sixteen control (age range: 2.3-61.3 years, 8 male), 37 VWM (2.4-56.5 years, 20 male), and 14 MLD (2.2-41.7 years, 6 male) subjects were included. Neurite density index and MWF were lower in patients than in controls (p < 0.001). Free water fraction was highest in VWM (p = 0.01), but similar to controls in MLD (p = 0.99). Changes in diffusion tensor imaging measures relative to controls were generally more pronounced in VWM than in MLD. In both patient groups, MCR-DIMWI MWF correlated strongest with clinical measures.

CONCLUSION

Quantitative MRI correlates to clinical measures and yields differential profiles in VWM and MLD, in line with differences in neuropathology.

KEY POINTS

Question Can quantitative MRI reflect known neuropathological differences and correlate with clinical scores for these leukodystrophies? Finding Quantitative MRI measures, e.g., MWF, neurite density index, and free water fraction differ between leukodystrophies and controls, in correspondence to known histological differences. Clinical relevance MRI techniques producing quantitative, biologically-specific, measures regarding the health of myelin and axons deliver more comprehensive information regarding pathological changes in leukodystrophies than current approaches, and are thus viable tools for monitoring patients and providing clinical trial outcome measures.

Whole brain multiparametric mapping in two minutes using a dual-flip-angle stack-of-stars blipped multi-gradient-echo acquisition

A new MRI technique is presented for three-dimensional fast simultaneous whole brain mapping of myelin water fraction (MWF), T1, proton density (PD), R2*, magnetic susceptibility (QSM), and B1 transmit field (B1+). Phantom and human (N = 9) datasets were acquired using a dual-flip-angle blipped multi-gradient-echo (DFA-mGRE) sequence with a stack-of-stars (SOS) trajectory. Images were reconstructed using a subspace-based algorithm with a locally low-rank constraint. A novel joint-sparsity-constrained multicomponent T2*-T1 spectrum estimation (JMSE) algorithm is proposed to correct for the T1 saturation effect and B1+/B1- inhomogeneities in the quantification of MWF. A tissue-prior-based B1+ estimation algorithm was adapted for B1 correction in the mapping of T1 and PD. In the phantom study, measurements obtained at an acceleration factor (R) of 12 using prospectively under-sampled SOS showed good consistency (R2 > 0.997) with Cartesian reference for R2*/T1app/M0app. In the in vivo study, results of retrospectively under-sampled SOS with R = 6, 12, 18, showed good quality (structure similarity index measure > 0.95) compared with those of fully-sampled SOS. Besides, results of prospectively under-sampled SOS with R = 12 showed good consistency (intraclass correlation coefficient > 0.91) with Cartesian reference for T1/PD/B1+/MWF/QSM/R2*, and good reproducibility (coefficient of variation < 7.0 %) in the test-retest analysis for T1/PD/B1+/MWF/R2*. This study has demonstrated the feasibility of simultaneous whole brain multiparametric mapping with a two-minute scan using the DFA-mGRE SOS sequence, which may overcome a major obstacle for neurological applications of multiparametric MRI.

Harnessing myelin water fraction as an imaging biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination: A review

Myelin water fraction (MWF) imaging has emerged as a promising magnetic resonance imaging (MRI) biomarker for investigating brain function and composition. This comprehensive review synthesizes the current state of knowledge on MWF as a biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination. The databases used include Web of Science, Scopus, Science Direct, and PubMed. We begin with a brief discussion of the theoretical foundations of MWF imaging, including its basis in MR physics and the mathematical modeling underlying its calculation, with an overview of the most adopted MRI methods of MWF imaging. Next, we delve into the clinical and research applications that have been explored to date, highlighting its advantages and limitations. Finally, we explore the potential of MWF to serve as a predictive biomarker for neurological disorders and identify future research directions for optimizing MWF imaging protocols and interpreting MWF in various contexts. By harnessing the power of MWF imaging, we may gain new insights into brain health and disease across the human lifespan, ultimately informing novel diagnostic and therapeutic strategies.

Quantification of the in vivo brain ultrashort-T2* component in healthy volunteers

PURPOSE

Recent work has shown MRI is able to measure and quantify signals of phospholipid membrane-bound protons associated with myelin in the human brain. This work seeks to develop an improved technique for characterizing this brain ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component in vivo accounting forT 1 $$ {\mathrm{T}}_1 $$ weighting.

METHODS

Data from ultrashort echo time scans from 16 healthy volunteers with variable flip angles (VFA) were collected and fitted into an advanced regression model to quantify signal fraction, relaxation time, and frequency shift of the ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component.

RESULTS

The fitted components show intra-subject differences of different white matter structures and significantly elevated ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ signal fraction in the corticospinal tracts measured at 0.09 versus 0.06 in other white matter structures and significantly elevated ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ frequency shift in the body of the corpus callosum at- $$ - $$ 1.5 versus- $$ - $$ 2.0 ppm in other white matter structures.

CONCLUSION

The significantly different measured components and measuredT 1 $$ {\mathrm{T}}_1 $$ relaxation time of the ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component suggest that this method is picking up novel signals from phospholipid membrane-bound protons.

Applicability of multiple quantitative magnetic resonance methods in genetic brain white matter disorders

BACKGROUND AND PURPOSE

Magnetic resonance imaging (MRI) measures of tissue microstructure are important for monitoring brain white matter (WM) disorders like leukodystrophies and multiple sclerosis. They should be sensitive to underlying pathological changes. Three whole-brain isotropic quantitative methods were applied and compared within a cohort of controls and leukodystrophy patients: two novel myelin water imaging (MWI) techniques (multi-compartment relaxometry diffusion-informed MWI: MCR-DIMWI, and multi-echo T2 relaxation imaging with compressed sensing: METRICS) and neurite orientation dispersion and density imaging (NODDI).

METHODS

For 9 patients with different leukodystrophies (age range 0.4-62.4 years) and 15 control subjects (2.3-61.3 years), T1-weighted MRI, fluid-attenuated inversion recovery, multi-echo gradient echo with variable flip angles, METRICS, and multi-shell diffusion-weighted imaging were acquired on 3 Tesla. MCR-DIMWI, METRICS, NODDI, and quality control measures were extracted to evaluate differences between patients and controls in WM and deep gray matter (GM) regions of interest (ROIs). Pearson correlations, effect size calculations, and multi-level analyses were performed.

RESULTS

MCR-DIMWI and METRICS-derived myelin water fractions (MWFs) were lower and relaxation times were higher in patients than in controls. Effect sizes of MWF values and relaxation times were large for both techniques. Differences between patients and controls were more pronounced in WM ROIs than in deep GM. MCR-DIMWI-MWFs were more homogeneous within ROIs and more bilaterally symmetrical than METRICS-MWFs. The neurite density index was more sensitive in detecting differences between patients and controls than fractional anisotropy. Most measures obtained from MCR-DIMWI, METRICS, NODDI, and diffusion tensor imaging correlated strongly with each other.

CONCLUSION

This proof-of-concept study shows that MCR-DIMWI, METRICS, and NODDI are sensitive techniques to detect changes in tissue microstructure in WM disorders.

Microstructure-Informed Myelin Mapping (MIMM) from Gradient Echo MRI using Stochastic Matching Pursuit

Quantification of the myelin content of the white matter is important for studying demyelination in neurodegenerative diseases such as Multiple Sclerosis (MS), particularly for longitudinal monitoring. A novel noninvasive MRI method, called Microstructure-Informed Myelin Mapping (MIMM), is developed to quantify the myelin volume fraction (MVF) by utilizing a multi gradient echo sequence (mGRE) and a detailed biophysical model of tissue microstructure. Myelin is modeled as anisotropic negative susceptibility source based on the Hollow Cylindrical Fiber Model (HCFM), and iron as isotropic positive susceptibility source in the extracellular region. Voxels with a range of biophysical parameters are simulated to create a dictionary of MR echo time magnitude signals and total susceptibility values. MRI signals measured using a mGRE sequence are then matched voxel-by-voxel to the created dictionary to obtain the spatial distributions of myelin and iron. Three different MIMM versions are presented to deal with the fiber orientation dependent susceptibility effects of the myelin sheaths: a basic variation, which assumes fiber orientation is an unknown to fit, two orientation informed variations, which assume the fiber orientation distribution is available either from a separate diffusion tensor imaging (DTI) acquisition or from a DTI atlas based fiber orientation map. While all showed a significant linear correlation with the reference method based on T2-relaxometry (p < 0.0001), DTI orientation informed and atlas orientation informed variations reduced overestimation at white matter tracts compared to the basic variation. Finally, the implications and usefulness of attaining an additional iron susceptibility distribution map are discussed.

Highlights

novel stochastic matching pursuit algorithm called microstructure-informed myelin mapping (MIMM) is developed to quantify Myelin Volume Fraction (MVF) using Magnetic Resonance Imaging (MRI) and microstructural modeling.utilizes a detailed biophysical model to capture the susceptibility effects on both magnitude and phase to quantify myelin and iron.matter fiber orientation effects are considered for the improved MVF quantification in the major fiber tracts.acquired myelin and iron maps may be utilized to monitor longitudinal disease progress.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Quantitative MRI in leukodystrophies

Leukodystrophies constitute a large and heterogeneous group of genetic diseases primarily affecting the white matter of the central nervous system. Different disorders target different white matter structural components. Leukodystrophies are most often progressive and fatal. In recent years, novel therapies are emerging and for an increasing number of leukodystrophies trials are being developed. Objective and quantitative metrics are needed to serve as outcome measures in trials. Quantitative MRI yields information on microstructural properties, such as myelin or axonal content and condition, and on the chemical composition of white matter, in a noninvasive fashion. By providing information on white matter microstructural involvement, quantitative MRI may contribute to the evaluation and monitoring of leukodystrophies. Many distinct MR techniques are available at different stages of development. While some are already clinically applicable, others are less far developed and have only or mainly been applied in healthy subjects. In this review, we explore the background, current status, potential and challenges of available quantitative MR techniques in the context of leukodystrophies.

On the performance of multi-compartment relaxometry for myelin water imaging (MCR-MWI) - test-retest repeatability and inter-protocol reproducibility

In this study, we optimized the variable flip angle (VFA) acquisition scheme using numerical simulations to shorten the acquisition time of multicompartment relaxometry for myelin water imaging (MCR-MWI) to a clinically practical range in the absence of advanced image reconstruction methods. As the primary objective of this study, the test-retest repeatability of myelin water fraction (MWF) measurements of MCR-MWI is evaluated on three gradient echo (GRE) sequence settings using the optimized VFA schemes with different echo times and repetition times, emulating various scanner setups. The cross-protocol reproducibility of MCR-MWI and MCR with diffusion-informed myelin water imaging (MCR-DIMWI) is also examined. As a secondary objective, we explore the bundle-specific profiles of various microstructural parameters from MCR-(DI)MWI and their cross-correlations to determine if these parameters possess supplementary microstructure information beyond myelin concentration. Numerical simulations indicate that MCR-MWI can be performed with a minimum of three flip angles covering a wide range of T₁ weightings without adding significant bias. This is supported by the results of an in vivo experiment, allowing whole-brain 1.5 mm isotropic MWF maps to be acquired in 9 min, reducing the total scan time to 40% of the original implementation without significant quality degradation. Good test-retest repeatability is observed for MCR-MWI for all three GRE protocols. While good correlations can also be found in MWF across protocols, systematic differences are observed. Bundle-specific MWF analysis reveals that certain white matter bundles are similar in all participants. We also found that microstructure relaxation parameters have low linear correlations with MWF. MCR-MWI is a reproducible measure of myelin. However, attention should be paid to the protocol related MWF differences when comparing different studies, as the MWF bias up to 0.5% can be observed across the protocols examined in this work.

Mitigating the impact of flip angle and orientation dependence in single compartment R2* estimates via 2-pool modeling

PURPOSE

The effective transverse relaxation rate (R 2 * $$ {\mathrm{R}}_2^{\ast } $$ ) is influenced by biological features that make it a useful means of probing brain microstructure. However, confounding factors such as dependence on flip angle (α) and fiber orientation with respect to the main field (θ $$ \uptheta $$ ) complicate interpretation. The α- andθ $$ \uptheta $$ -dependence stem from the existence of multiple sub-voxel micro-environments (e.g., myelin and non-myelin water compartments). Ordinarily, it is challenging to quantify these sub-compartments; therefore, neuroscientific studies commonly make the simplifying assumption of a mono-exponential decay obtaining a singleR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimate per voxel. In this work, we investigated how the multi-compartment nature of tissue microstructure affects single compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates.

METHODS

We used 2-pool (myelin and non-myelin water) simulations to characterize the bias in single compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates. Based on our numeric observations, we introduced a linear model that partitionsR 2 * $$ {\mathrm{R}}_2^{\ast } $$ into α-dependent and α-independent components and validated this in vivo at 7T. We investigated the dependence of both components on the sub-compartment properties and assessed their robustness, orientation dependence, and reproducibility empirically.

RESULTS

R 2 * $$ {\mathrm{R}}_2^{\ast } $$ increased with myelin water fraction and residency time leading to a linear dependence on α. We observed excellent agreement between our numeric and empirical results. Furthermore, the α-independent component of the proposed linear model was robust to the choice of α and reduced dependence on fiber orientation, although it suffered from marginally higher noise sensitivity.

CONCLUSION

We have demonstrated and validated a simple approach that mitigates flip angle and orientation biases in single-compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

Magnetic Resonance Imaging and Spectroscopy Analysis in a Pelizaeus-Merzbacher Disease Rat Model

Pelizaeus−Merzbacher disease (PMD) is an X-linked recessive disorder of the central nervous system. We performed 7 Tesla magnetic resonance imaging of the brain in Tama rats, a rodent PMD model, and control rats, as well as evaluated the diagnostic values. In the white matter of the Tama rats, the T2 values were prolonged, which is similar to that observed in patients with PMD (60.7 ± 1.8 ms vs. 51.6 ± 1.3 ms, p < 0.0001). The apparent diffusion coefficient values in the white matter of the Tama rats were higher than those of the control rats (0.68 ± 0.03 × 10−3 mm2/s vs. 0.64 ± 0.03 × 10−3 mm2/s, p < 0.05). In proton magnetic resonance spectroscopy, the N-acetylaspartate (6.97 ± 0.12 mM vs. 5.98 ± 0.25 mM, p < 0.01) and N-acetylaspartate + N-acetylaspartylglutamate values of the Tama rats were higher (8.22 ± 0.17 mM vs. 7.14 ± 0.35 mM, p < 0.01) than those of the control rats. The glycerophosphocholine + phosphocholine values of the Tama rats were lower than those of the control rats (1.04 ± 0.09 mM vs. 1.45 ± 0.04 mM, p < 0.001). By using Luxol fast blue staining, we confirmed dysmyelination in the Tama rats. These results are similar to those of patients with PMD and other PMD animal models.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

Decoding the microstructural properties of white matter using realistic models

Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed. In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

Multi-compartment analysis of the complex gradient-echo signal quantifies myelin breakdown in premanifest Huntington's disease

White matter (WM) alterations have been identified as a relevant pathological feature of Huntington's disease (HD). Increasing evidence suggests that WM changes in this disorder are due to alterations in myelin-associated biological processes. Multi-compartmental analysis of the complex gradient-echo MRI signal evolution in WM has been shown to quantify myelin in vivo, therefore pointing to the potential of this technique for the study of WM myelin changes in health and disease. This study first characterized the reproducibility of metrics derived from the complex multi-echo gradient-recalled echo (mGRE) signal across the corpus callosum in healthy participants, finding highest reproducibility in the posterior callosal segment. Subsequently, the same analysis pipeline was applied in this callosal region in a sample of premanifest HD patients (n = 19) and age, sex and education matched healthy controls (n = 21). In particular, we focused on two myelin-associated derivatives: i. the myelin water signal fraction (fm), a parameter dependent on myelin content; and ii. The difference in frequency between myelin and intra-axonal water pools (Δω), a parameter dependent on the ratio between the inner and the outer axonal radii. fm was found to be lower in HD patients (β = -0.13, p = 0.03), while Δω did not show a group effect. Performance in tests of working memory, executive function, social cognition and movement was also assessed, and a greater age-related decline in executive function was detected in HD patients (β = -0.06, p = 0.006), replicating previous evidence of executive dysfunction in HD. Finally, the correlation between fm, executive function, and proximity to disease onset was explored in patients, and a positive correlation between executive function and fm was detected (r = 0.542; p = 0.02). This study emphasises the potential of complex mGRE signal analysis for aiding understanding of HD pathogenesis and progression. Moreover, expanding on evidence from pathology and animal studies, it provides novel in vivo evidence supporting myelin breakdown as an early feature of HD.

Imperfect spoiling in variable flip angle T1 mapping at 7T: Quantifying and minimizing impact

Purpose

The variable flip angle (VFA) approach to T₁ mapping assumes perfectly spoiled transverse magnetisation at the end of each repetition time (TR). Despite radiofrequency (RF) and gradient spoiling, this condition is rarely met, leading to erroneous T₁ estimates (T1app). Theoretical corrections can be applied but make assumptions about tissue properties, for example, a global T₂ time. Here, we investigate the effect of imperfect spoiling at 7T and the interaction between the RF and gradient spoiling conditions, additionally accounting for diffusion. We provide guidance on the optimal approach to maximise the accuracy of the T₁ estimate in the context of 3D multi‐echo acquisitions.

Methods

The impact of the spoiling regime was investigated through numerical simulations, phantom and invivo experiments.

Results

The predicted dependence of T1app on tissue properties, system settings, and spoiling conditions was observed in both phantom and in vivo experiments. Diffusion effects modulated the dependence of T1app on both B1+ efficiency and T₂ times.

Conclusion

Error in T1app can be minimized by using an RF spoiling increment and gradient spoiler moment combination that minimizes T₂‐dependence and safeguards image quality. Although the diffusion effect was comparatively small at 7T, correction factors accounting for this effect are recommended.

Click here for author‐reader discussions

Variable flip angle echo planar time-resolved imaging (vFA-EPTI) for fast high-resolution gradient echo myelin water imaging

Myelin water imaging techniques based on multi-compartment relaxometry have been developed as an important tool to measure myelin concentration in vivo, but are limited by the long scan time of multi-contrast multi-echo acquisition. In this work, a fast imaging technique, termed variable flip angle Echo Planar Time-Resolved Imaging (vFA-EPTI), is developed to acquire multi-echo and multi-flip-angle gradient-echo data with significantly reduced acquisition time, providing rich information for multi-compartment analysis of gradient-echo myelin water imaging (GRE-MWI). The proposed vFA-EPTI method achieved 26 folds acceleration with good accuracy by utilizing an efficient continuous readout, optimized spatiotemporal encoding across echoes and flip angles, as well as a joint subspace reconstruction. An approach to estimate off-resonance field changes between different flip-angle acquisitions was also developed to ensure high-quality joint reconstruction across flip angles. The accuracy of myelin water fraction (MWF) estimate under high acceleration was first validated by a retrospective undersampling experiment using a lengthy fully-sampled data as reference. Prospective experiments were then performed where whole-brain MWF and multi-compartment quantitative maps were obtained in 5 min at 1.5 mm isotropic resolution and 24 min at 1 mm isotropic resolution at 3T. Additionally, ultra-high resolution data at 600 μm isotropic resolution were acquired at 7T, which show detailed structures within the cortex such as the line of Gennari, demonstrating the ability of the proposed method for submillimeter GRE-MWI that can be used to study cortical myeloarchitecture in vivo.

Probing the myelin water compartment with a saturation-recovery, multi-echo gradient-recalled echo sequence

PURPOSE

To investigate the effect of varying levels of T 1 -weighting on the evolution of the complex signal from white matter in a multi-echo gradient-recalled echo (mGRE) saturation-recovery sequence.

THEORY AND METHODS

Analysis of the complex signal evolution in an mGRE sequence allows the contributions from short- and long- T 2 ∗ components to be separated, thus providing a measure of the relative strength of signals from the myelin water, and the external and intra-axonal compartments. Here we evaluated the effect of different levels of T 1 -weighting on these signals, expecting that the previously reported, short T 1 of the myelin water would lead to a relative enhancement of the myelin water signal in the presence of signal saturation. Complex, saturation-recovery mGRE data from the splenium of the corpus callosum from 5 healthy volunteers were preprocessed using a frequency difference mapping (FDM) approach and analyzed using the 3-pool model of complex signal evolution in white matter.

RESULTS

An increase in the apparent T 1 as a function of echo time was demonstrated, but this increase was an order of magnitude smaller than that expected from previously reported myelin water T 1 -values. This suggests the presence of magnetization transfer and exchange effects which counteract the T 1 -weighting.

CONCLUSION

Variation of the B 1 + amplitude in a saturation-recovery mGRE sequence can be used to modulate the relative strength of signals from the different compartments in white matter, but the modulation is less than predicted from previously reported T 1 -values.