Rapid whole brain myelin water content mapping without an external water standard at 1.5 T

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.

Association of brain tissue cerebrospinal fluid fraction with age in healthy cognitively normal adults

Background and purpose

Our objective was to apply multi-compartment T2 relaxometry in cognitively normal individuals aged 20-80 years to study the effect of aging on the parenchymal CSF fraction (CSFF), a potential measure of the subvoxel CSF space.

Materials and methods

A total of 60 volunteers (age range, 22-80 years) were enrolled. Voxel-wise maps of short-T2 myelin water fraction (MWF), intermediate-T2 intra/extra-cellular water fraction (IEWF), and long-T2 CSFF were obtained using fast acquisition with spiral trajectory and adiabatic T2prep (FAST-T2) sequence and three-pool non-linear least squares fitting. Multiple linear regression analyses were performed to study the association between age and regional MWF, IEWF, and CSFF measurements, adjusting for sex and region of interest (ROI) volume. ROIs include the cerebral white matter (WM), cerebral cortex, and subcortical deep gray matter (GM). In each model, a quadratic term for age was tested using an ANOVA test. A Spearman's correlation between the normalized lateral ventricle volume, a measure of organ-level CSF space, and the regional CSFF, a measure of tissue-level CSF space, was computed.

Results

Regression analyses showed that there was a statistically significant quadratic relationship with age for CSFF in the cortex (p = 0.018), MWF in the cerebral WM (p = 0.033), deep GM (p = 0.017) and cortex (p = 0.029); and IEWF in the deep GM (p = 0.033). There was a statistically highly significant positive linear relationship between age and regional CSFF in the cerebral WM (p < 0.001) and deep GM (p < 0.001). In addition, there was a statistically significant negative linear association between IEWF and age in the cerebral WM (p = 0.017) and cortex (p < 0.001). In the univariate correlation analysis, the normalized lateral ventricle volume correlated with the regional CSFF measurement in the cerebral WM (ρ = 0.64, p < 0.001), cortex (ρ = 0.62, p < 0.001), and deep GM (ρ = 0.66, p < 0.001).

Conclusion

Our cross-sectional data demonstrate that brain tissue water in different compartments shows complex age-dependent patterns. Parenchymal CSFF, a measure of subvoxel CSF-like water in the brain tissue, is quadratically associated with age in the cerebral cortex and linearly associated with age in the cerebral deep GM and WM.

Magnetic Susceptibility Source Separation Solely from Gradient Echo Data: Histological Validation

Quantitative susceptibility mapping (QSM) facilitates mapping of the bulk magnetic susceptibility of tissue from the phase of complex gradient echo (GRE) MRI data. QSM phase processing combined with an R2* model of magnitude of multiecho gradient echo data (R2*QSM) allows separation of dia- and para-magnetic components (e.g., myelin and iron) that contribute constructively to R2* value but destructively to the QSM value of a voxel. This R2*QSM technique is validated against quantitative histology—optical density of myelin basic protein and Perls’ iron histological stains of rim and core of 10 ex vivo multiple sclerosis lesions, as well as neighboring normal appearing white matter. We found that R2*QSM source maps are in good qualitative agreement with histology, e.g., showing increased iron concentration at the edge of the rim+ lesions and myelin loss in the lesions’ core. Furthermore, our results indicate statistically significant correlation between paramagnetic and diamagnetic tissue components estimated with R2*QSM and optical densities of Perls’ and MPB stains. These findings provide direct support for the use of R2*QSM magnetic source separation based solely on GRE complex data to characterize MS lesion composition.

Susceptibility source separation from gradient echo data using magnitude decay modeling

BACKGROUND AND PURPOSE

The objective is to demonstrate feasibility of separating magnetic sources in quantitative susceptibility mapping (QSM) by incorporating magnitude decay rates R 2 ∗ $R_2^{\rm{*}}$ in gradient echo (GRE) MRI.

METHODS

Magnetic susceptibility source separation was developed using R 2 ∗ $R_2^{\rm{*}}$ and compared with a prior method using R 2 ' = R 2 ∗ - R 2 ${R^{\prime}_2} = R_2^* - {R_2}$ that required an additional sequence to measure the transverse relaxation rate R₂ . Both susceptibility separation methods were compared in multiple sclerosis (MS) patients (n = 17). Susceptibility values of negative sources estimated with R 2 ∗ $R_2^{\rm{*}}$ -based source separation in a set of enhancing MS lesions (n = 44) were correlated against longitudinal myelin water fraction (MWF) changes.

RESULTS

In in vivo data, linear regression of the estimated χ + ${\chi}^{+}$ and χ - ${\chi}^{-}$ susceptibility values between the R 2 ∗ $R_2^*$ - and the R 2 ' ${R^{\prime}_2}$ -based separation methods performed across 182 segmented lesions revealed correlation coefficient r = .96 and slope close .99. Correlation analysis in enhancing lesions revealed a significant positive association between the χ - ${\chi}^{-}$ increase at 1-year post-onset relative to 0 year and the MWF increase at 1 year relative to 0 year (β = -0.144, 95% confidence interval: [-0.199, -0.1], p = .0008) and good agreement between R 2 ' ${R^{\prime}_2}$ and R 2 ∗ $R_2^*$ methods (r = .79, slope = .95).

CONCLUSIONS

Separation of magnetic sources based solely on GRE complex data is feasible by combining magnitude decay rate modeling and phase-based QSM and χ - ${\chi}^{-}$ change may serve as a biomarker for myelin recovery or damage in acute MS lesions.

Quantitative Water Permeability Mapping of Blood-Brain-Barrier Dysfunction in Aging

Blood-brain-barrier (BBB) dysfunction is a hallmark of aging and aging-related disorders, including cerebral small vessel disease and Alzheimer's disease. An emerging biomarker of BBB dysfunction is BBB water exchange rate (k_W) as measured by diffusion-weighted arterial spin labeling (DW-ASL) MRI. We developed an improved DW-ASL sequence for Quantitative Permeability Mapping and evaluated whole brain and region-specific k_W in a cohort of 30 adults without dementia across the age spectrum. In this cross-sectional study, we found higher k_W values in the cerebral cortex (mean = 81.51 min^-1, SD = 15.54) compared to cerebral white matter (mean = 75.19 min^-1, SD = 13.85) (p < 0.0001). We found a similar relationship for cerebral blood flow (CBF), concordant with previously published studies. Multiple linear regression analysis with k_W as an outcome showed that age was statistically significant in the cerebral cortex (p = 0.013), cerebral white matter (p = 0.033), hippocampi (p = 0.043), orbitofrontal cortices (p = 0.042), and precunei cortices (p = 0.009), after adjusting for sex and number of vascular risk factors. With CBF as an outcome, age was statistically significant only in the cerebral cortex (p = 0.026) and precunei cortices (p = 0.020). We further found moderate negative correlations between white matter hyperintensity (WMH) k_W and WMH volume (r = -0.51, p = 0.02), and normal-appearing white matter (NAWM) and WMH volume (r = -0.44, p = 0.05). This work illuminates the relationship between BBB water exchange and aging and may serve as the basis for BBB-targeted therapies for aging-related brain disorders.

Feasibility of Magnetic Resonance Fingerprinting on Aging MRI Hardware

The purpose of this work is to evaluate the feasibility of performing magnetic resonance fingerprinting (MRF) on older and lower-performance MRI hardware as a means to bring advanced imaging to the aging MRI install base. Phantom and in vivo experiments were performed on a 1.5T Siemens Aera (installed 2015) and 1.5T Siemens Symphony (installed 2002). A 2D spiral MRF sequence for simultaneous T₁/T₂/M₀ mapping was implemented on both scanners with different gradient trajectories to accommodate system specifications. In phantom, for T₁/T₂ values in a physiologically relevant range (T₁: 195-1539 ms; T₂: 20-267 ms), scanners had strong correlation (R² > 0.999) with average absolute percent difference of 8.1% and 10.1%, respectively. Comparison of the two trajectories on the newer scanner showed differences of 2.6% (T₁) and 10.9% (T₂), suggesting a partial explanation of the observed inter-scanner bias. Inter-scanner agreement was better when the same trajectory was used, with differences of 6.0% (T₁) and 4.0% (T₂). Intra-scanner coefficient of variation (CV) of T₁ and T₂ estimates in phantom were <2.0% and in vivo were ≤3.5%. In vivo inter-scanner white matter CV was 4.8% (T₁) and 5.1% (T₂). White matter measurements on the aging scanner after two months were consistent, with differences of 1.9% (T₁) and 3.9% (T₂). In conclusion, MRF is feasible on an aging MRI scanner and required only changes to the gradient trajectory.

Water content changes in new multiple sclerosis lesions have a minimal effect on the determination of myelin water fraction values

BACKGROUND AND PURPOSE

Myelin water fraction (MWF) is a histopathologically validated in vivo myelin marker. As MWF is the proportion of water with a short T₂ relative to the total water, increases in water from edema and inflammation may confound MWF determination in multiple sclerosis (MS) lesions. Total water content (TWC) measurement enables calculation of absolute myelin water content (MWC) and can be used to distinguish edema/inflammation from demyelination. We assessed what influence changes in total water might have on MWF by calculating MWC values in new MS lesions.

METHODS

3T 32-echo T₂ relaxation data were collected monthly for 6 months from six relapsing-remitting MS participants. TWC was determined and multiplied with MWF images to calculate corrected MWC images. The effect of this water content correction was examined in 20 new lesions by comparing mean MWF and MWC over time.

RESULTS

On average, at lesion first appearance, lesion TWC increased by 6.4% (p = .003; range: -1% to +21%), MWF decreased by 24% (p = .006; range: -70% to +12%), and MWC decreased by 20% (p = .026; range: -68% to +21%), relative to prelesion values. Average TWC in lesions then gradually decreased, whereas MWF and MWC remained low. The shape of the MWF and MWC lesion evolution curves was nearly identical, differing only by an offset.

CONCLUSION

MWF mirrors MWC and is able to monitor myelin in new lesions. Even after taking into account water content increases, MWC still decreased at lesion first appearance attributed to demyelination.

Impact of Lesion Location on Longitudinal Myelin Water Fraction Change in Chronic Multiple Sclerosis Lesions

BACKGROUND AND PURPOSE

To examine the impact of lesion location on longitudinal myelin water fraction (MWF) changes in chronic multiple sclerosis (MS) lesions. Relative hypoxia, due to vascular watershed regions of the cerebrum, has been implicated in lesion development but impact on ongoing demyelination is unknown.

METHODS

Forty-eight patients with relapsing-remitting and secondary progressive MS had two MWF scans with fast acquisition, spiral trajectory, and T2prep (FAST-T2) sequence, at an interval of 2.0 (±.3) years. Lesion location was identified based upon cerebral lobe and relation to the ventricles. Change in MWF was assessed using a mixed effects model, controlling for lesion location and patient covariates.

RESULTS

Average age was 42.3 (±12) years, mean disease duration was 9.7 (±9.1) years, and median Expanded Disability Status Score (EDSS) was 2.5 (±2.3). The majority of 512 chronic lesions was located in the frontal and parietal lobes (75.6%) and more often periventricular (44.7%). All occipital lesions were periventricular. The average lesion MWF decreased from baseline (.07 ± .03) to 2 years (.06 ±.03) P < .01. Lesions within the occipital lobe showed a significant reduction in MWF as compared to other lobes.

CONCLUSIONS

Chronic lesions in the occipital lobe showed the greatest reduction in MWF. Neuroanatomical localization of lesions to the occipital horns of the lateral ventricles, a watershed region, may contribute to ongoing demyelination in this lesion type.

Imaging outcome measures of neuroprotection and repair in MS: A consensus statement from NAIMS

OBJECTIVE

To summarize current and emerging imaging techniques that can be used to assess neuroprotection and repair in multiple sclerosis (MS), and to provide a consensus opinion on the potential utility of each technique in clinical trial settings.

METHODS

Clinicians and scientists with expertise in the use of MRI in MS convened in Toronto, Canada, in November 2016 at a North American Imaging in Multiple Sclerosis (NAIMS) Cooperative workshop meeting. The discussion was compiled into a manuscript and circulated to all NAIMS members in attendance. Edits and feedback were incorporated until all authors were in agreement.

RESULTS

A wide spectrum of imaging techniques and analysis methods in the context of specific study designs were discussed, with a focus on the utility and limitations of applying each technique to assess neuroprotection and repair. Techniques were discussed under specific themes, and included conventional imaging, magnetization transfer ratio, diffusion tensor imaging, susceptibility-weighted imaging, imaging cortical lesions, magnetic resonance spectroscopy, PET, advanced diffusion imaging, sodium imaging, multimodal techniques, imaging of special regions, statistical considerations, and study design.

CONCLUSIONS

Imaging biomarkers of neuroprotection and repair are an unmet need in MS. There are a number of promising techniques with different strengths and limitations, and selection of a specific technique will depend on a number of factors, notably the question the trial seeks to answer. Ongoing collaborative efforts will enable further refinement and improved methods to image the effect of novel therapeutic agents that exert benefit in MS predominately through neuroprotective and reparative mechanisms.

Yet more evidence that myelin protons can be directly imaged with UTE sequences on a clinical 3T scanner: Bicomponent T2* analysis of native and deuterated ovine brain specimens

PURPOSE

UTE sequences with a minimal nominal TE of 8 µs have shown promise for direct imaging of myelin protons (T2 , < 1 ms). However, there is still debate about the efficiency of 2D slice-selective UTE sequences in exciting myelin protons because the half excitation pulses used in these sequences have a relatively long duration (e.g., 0.3-0.6 ms). Here, we compared UTE and inversion-recovery (IR) UTE sequences used with either hard or half excitation pulses (durations 32 µs or 472 µs, respectively) for imaging myelin in native and deuterated ovine brain at 3T.

METHODS

Freshly frozen ovine brains were dissected into ∼2 mm-thick pure white matter and ∼3 to 8 mm-thick cerebral hemisphere specimens, which were imaged before and/or after different immersion time in deuterium oxide.

RESULTS

Bicomponent T2* analysis of UTE signals obtained with hard excitation pulses detected an ultrashort T2 component (STC) fraction (fS ) of 0% to 10% in native specimens, and up to ∼86% in heavily deuterated specimens. fS values were significantly affected by the TIs used in IR-UTE sequences with either hard or half excitation pulses in native specimens but not in heavily deuterated specimens. The STC T2* was in the range of 150 to 400 µs in all UTE and IR-UTE measurements obtained with either hard or half excitation pulses.

CONCLUSION

Our results further support myelin protons as the major source of the ultrashort T2* signals seen on IR-UTE images and demonstrate the potential of IR-UTE sequences with half excitation pulses for directly imaging myelin using clinical scanners. Magn Reson Med 80:538-547, 2017. © 2017 International Society for Magnetic Resonance in Medicine.