Effect of intercompartmental water exchange on the apparent myelin water fraction in multiexponential T2 measurements of rat spinal cord

Rapid multicomponent relaxometry in steady state with correction of magnetization transfer effects

PURPOSE

To study the effects of magnetization transfer (MT) on multicomponent T2 parameters obtained using mcDESPOT in macromolecule-rich tissues and to propose a new method called mcRISE to correct MT-induced biases.

METHODS

The two-pool mcDESPOT model was modified by the addition of an exchanging macromolecule proton pool to model the MT effect in cartilage. The mcRISE acquisition scheme was developed to provide sensitivity to all pools. An incremental fitting was applied to estimate MT and relaxometry parameters with minimized coupling. The interaction between MT and relaxometry parameters, efficacy of MT correction, and feasibility of mcRISE in vivo were investigated in simulations and in healthy volunteers.

RESULTS

The MT effect caused significant errors in multicomponent T1/T2 values and in fast-relaxing water fraction fF , which is consistent with previous experimental observations. fF increased significantly with macromolecule content if MT was ignored. mcRISE resulted in a multifold reduction of MT biases and yielded decoupled multicomponent T1/T2 relaxometry and quantitative MT parameters.

CONCLUSION

mcRISE is an efficient approach for correcting MT biases in multicomponent relaxometry based on steady state sequences. Improved specificity of mcRISE may help to elucidate the sources of the previously described high sensitivity of noncorrected mcDESPOT parameters to disease-related changes in cartilage and the brain.

Fast and simplified mapping of mean axon diameter using temporal diffusion spectroscopy

Mapping axon diameter is of interest for the potential diagnosis and monitoring of various neuronal pathologies. Advanced diffusion-weighted MRI methods have been developed to measure mean axon diameters non-invasively, but suffer major drawbacks that prevent their direct translation into clinical practice, such as complex non-linear data fitting and, more importantly, long scanning times that are usually not tolerable for most human subjects. In the current study, temporal diffusion spectroscopy using oscillating diffusion gradients was used to measure mean axon diameters with high sensitivity to small axons in the central nervous system. Axon diameters have been found to be correlated with a novel metric, DDR⊥ (the rate of dispersion of the perpendicular diffusion coefficient with gradient frequency), which is a model-free quantity that does not require complex data analyses and can be obtained from two diffusion coefficient measurements in clinically relevant times with conventional MRI machines. A comprehensive investigation including computer simulations and animal experiments ex vivo showed that measurements of DDR⊥ agree closely with histological data. In humans in vivo, DDR⊥ was also found to correlate well with reported mean axon diameters in human corpus callosum, and the total scan time was only about 8 min. In conclusion, DDR⊥ may have potential to serve as a fast, simple and model-free approach to map the mean axon diameter of white matter in clinics for assessing axon diameter changes.

Modeling T(1) and T(2) relaxation in bovine white matter

The fundamental basis of T1 and T2 contrast in brain MRI is not well understood; recent literature contains conflicting views on the nature of relaxation in white matter (WM). We investigated the effects of inversion pulse bandwidth on measurements of T1 and T2 in WM. Hybrid inversion-recovery/Carr-Purcell-Meiboom-Gill experiments with broad or narrow bandwidth inversion pulses were applied to bovine WM in vitro. Data were analysed with the commonly used 1D-non-negative least squares (NNLS) algorithm, a 2D-NNLS algorithm, and a four-pool model which was based upon microscopically distinguishable WM compartments (myelin non-aqueous protons, myelin water, non-myelin non-aqueous protons and intra/extracellular water) and incorporated magnetization exchange between adjacent compartments. 1D-NNLS showed that different T2 components had different T1 behaviours and yielded dissimilar results for the two inversion conditions. 2D-NNLS revealed significantly more complicated T1/T2 distributions for narrow bandwidth than for broad bandwidth inversion pulses. The four-pool model fits allow physical interpretation of the parameters, fit better than the NNLS techniques, and fits results from both inversion conditions using the same parameters. The results demonstrate that exchange cannot be neglected when analysing experimental inversion recovery data from WM, in part because it can introduce exponential components having negative amplitude coefficients that cannot be correctly modeled with nonnegative fitting techniques. While assignment of an individual T1 to one particular pool is not possible, the results suggest that under carefully controlled experimental conditions the amplitude of an apparent short T1 component might be used to quantify myelin water.

Simulation of changes in diffusion related to different pathologies at cellular level after traumatic brain injury

PURPOSE

The goal of the current study was to investigate tissue pathology at the cellular level in traumatic brain injury (TBI) as revealed by Monte Carlo simulation of diffusion tensor imaging (DTI)-derived parameters and elucidate the possible sources of conflicting findings of DTI abnormalities as reported in the TBI literature.

METHODS

A model with three compartments separated by permeable membranes was employed to represent the diffusion environment of water molecules in brain white matter. The dynamic diffusion process was simulated with a Monte Carlo method using adjustable parameters of intra-axonal diffusivity, axon separation, glial cell volume fraction, and myelin sheath permeability. The effects of tissue pathology on DTI parameters were investigated by adjusting the parameters of the model corresponding to different stages of brain injury.

RESULTS

The results suggest that the model is appropriate and the DTI-derived parameters simulate the predominant cellular pathology after TBI. Our results further indicate that when edema is not prevalent, axial and radial diffusivity have better sensitivity to axonal injury and demyelination than other DTI parameters.

CONCLUSION

DTI is a promising biomarker to detect and stage tissue injury after TBI. The observed inconsistencies among previous studies are likely due to scanning at different stages of tissue injury after TBI. Magn Reson Med 76:290-300, 2016. © 2015 Wiley Periodicals, Inc.

Physiological noise compensation in gradient-echo myelin water imaging

In MRI, physiological noise which originates from cardiac and respiratory functions can induce substantial errors in detecting small signals in the brain. In this work, we explored the effects of the physiological noise and their compensation methods in gradient-echo myelin water imaging (GRE-MWI). To reduce the cardiac function induced inflow noise, flow saturation RF pulses were applied to the inferior portion of the head, saturating inflow blood signals. For the respiratory function induced B0 fluctuation compensation, a navigator echo was acquired, and respiration induced phase errors were corrected during reconstruction. After the compensations, the resulting myelin water images show substantially improved image quality and reproducibility. These improvements confirm the importance and usefulness of the physiological noise compensations in GRE-MWI.

The microstructural correlates of T1 in white matter

PURPOSE

Several studies have shown strong correlations between myelin content and T1 within the brain, and have even suggested that T1 can be used to estimate myelin content. However, other micro-anatomical features such as compartment size are known to affect longitudinal relaxation rates, similar to compartment size effects in porous media.

METHODS

T1 measurements were compared with measured or otherwise published axon size measurements in white matter tracts of the rat spinal cord, rat brain, and human brain.

RESULTS

In both ex vivo and in vivo studies, correlations were present between the relaxation rate 1/T1 and axon size across regions of rat spinal cord with nearly equal myelin content.

CONCLUSION

While myelination is likely the dominant determinant of T1 in white matter, variations in white matter microstructure, independent of myelin volume fraction, may also be reflected in T1 differences between regions or subjects.

Mapping mean axon diameter and axonal volume fraction by MRI using temporal diffusion spectroscopy

Mapping mean axon diameter and intra-axonal volume fraction may have significant clinical potential because nerve conduction velocity is directly dependent on axon diameter, and several neurodegenerative diseases affect axons of specific sizes and alter axon counts. Diffusion-weighted MRI methods based on the pulsed gradient spin echo (PGSE) sequence have been reported to be able to assess axon diameter and volume fraction non-invasively. However, due to the relatively long diffusion times used, e.g. >20ms, the sensitivity to small axons (diameter<2μm) is low, and the derived mean axon diameter has been reported to be overestimated. In the current study, oscillating gradient spin echo (OGSE) diffusion sequences with variable frequency gradients were used to assess rat spinal white matter tracts with relatively short effective diffusion times (1-5ms). In contrast to previous PGSE-based methods, the extra-axonal diffusion cannot be modeled as hindered (Gaussian) diffusion when short diffusion times are used. Appropriate frequency-dependent rates are therefore incorporated into our analysis and validated by histology-based computer simulation of water diffusion. OGSE data were analyzed to derive mean axon diameters and intra-axonal volume fractions of rat spinal white matter tracts (mean axon diameter of ~1.27-5.54μm). The estimated values were in good agreement with histology, including the small axon diameters (<2.5μm). This study establishes a framework for the quantification of nerve morphology using the OGSE method with high sensitivity to small axons.

Probing signal phase in direct visualization of short transverse relaxation time component (ViSTa)

PURPOSE

To demonstrate the phase evolutions of direct visualization of short transverse relaxation time component (ViSTa) matches with those of myelin water.

METHOD

Myelin water imaging (MWI) measures short transverse signals and has been suggested as a biomarker for myelin. Recently, a new approach, ViSTa, has been proposed to acquire short T2* signals by suppressing long T1 signals. This method does not require any ill-conditioned data processing and therefore provides high-quality images. In this study, the phase of the ViSTa signal was compared with the phase of myelin water simulated by the magnetic susceptibility model of hollow cylinder.

RESULTS

The phase evolutions of the ViSTa signal were similar to the simulated myelin water phase evolutions. When fiber orientation was perpendicular relative to the main magnetic field, both the ViSTa and the simulated myelin water phase showed large positive frequency shifts, whereas the gradient echo phase showed a slightly negative frequency shift. Additionally, the myelin water phase map generated using diffusion tensor imaging (DTI) information revealed a good match with the ViSTa phase image.

CONCLUSION

The results of this study support the origin of ViSTa signal as myelin water. ViSTa phase may potentially provide sensitivity to demyelination.

MRI-based myelin water imaging: A technical review

Multiexponential T2 relaxation time measurement in the central nervous system shows a component that originates from water trapped between the lipid bilayers of myelin. This myelin water component is of significant interest as it provides a myelin-specific MRI signal of value in assessing myelin changes in cerebral white matter in vivo. In this article, the various acquisition and analysis strategies proposed to date for myelin water imaging are reviewed and research conducted into their validity and clinical applicability is presented. Comparisons between the imaging methods are made with a discussion regarding potential difficulties and model limitations.

In vivo longitudinal Myelin Water Imaging in rat spinal cord following dorsal column transection injury

Longitudinal Myelin Water Imaging was carried out in vivo to characterize white matter damage following dorsal column transection (DC Tx) injury at the lumbar level L1 of rat spinal cords. A transmit-receive implantable coil system was used to acquire multiple spin-echo (MSE) quantitative T2 data from the lumbar spinal cords of 16 rats at one week pre-injury as well as 3 and 8weeks post-injury (117 microns in-plane resolution and 1.5mm slice thickness). In addition, ex vivo MSE and DTI data were acquired from cords fixed and excised at 3 or 8weeks post injury using a solenoid coil. The MSE data were used to generate Myelin Water Fractions (MWFs) as a surrogate measure of myelin content, while DTI data were acquired to study damage to the axons. Myelin damage was assessed histologically with Eriochrome cyanine (EC) and Myelin Basic Protein in degenerated myelin (dgen-MBP) staining, and axonal damage was assessed by neurofilament-H in combination with neuron specific beta-III-tubulin (NF/Tub) staining. These MRI and histological measures of injury were studied in the dorsal column at 5mm cranial and 5mm caudal to injury epicenter. MWF increased significantly at 3weeks post-injury at both the cranial and caudal sites, relative to baseline. The values on the cranial side of injury returned to baseline at 8weeks post-injury but remained elevated on the caudal side. This trend was found in both in vivo and ex vivo data. This MWF increase was likely due to the presence of myelin debris, which were cleared by 8 weeks on the cranial, but not the caudal, side. Both EC and dgen-MBP stains displayed similar trends. MWF showed significant correlation with EC staining (R=0.63, p=0.005 in vivo and R=0.74, p=0.0001 ex vivo). MWF also correlated strongly with the dgen-MBP stain, but only on the cranial side (R=0.64, p=0.05 in vivo; R=0.63, p=0.038 ex vivo). This study demonstrates that longitudinal MWI in vivo can accurately characterize white matter damage in DC Tx model of injury in the rat spinal cord.

Quantitative MRI and ultrastructural examination of the cuprizone mouse model of demyelination

The cuprizone mouse model of demyelination was used to investigate the influence that white matter changes have on different magnetic resonance imaging results. In vivo T2 -weighted and magnetization transfer images (MTIs) were acquired weekly in control (n = 5) and cuprizone-fed (n = 5) mice, with significant increases in signal intensity in T2 -weighted images (p < 0.001) and lower magnetization transfer ratio (p < 0.001) in the corpus callosum of the cuprizone-fed mice starting at 3 weeks and peaking at 4 and 5 weeks, respectively. Diffusion tensor imaging (DTI), quantitative MTI (qMTI), and T1/T2 measurements were used to analyze freshly excised tissue after 6 weeks of cuprizone administration. In multicomponent T2 analysis with 10 ms echo spacing, there was no visible myelin water component associated with the short T2 value. Quantitative MTI metrics showed significant differences in the corpus callosum and external capsule of the cuprizone-fed mice, similar to previous studies of multiple sclerosis in humans and animal models of demyelination. Fractional anisotropy was significantly lower and mean, axial, and radial diffusivity were significantly higher in the cuprizone-fed mice. Cellular distributions measured in electron micrographs of the corpus callosum correlated strongly to several different quantitative MRI metrics. The largest Spearman correlation coefficient varied depending on cellular type: T1 versus the myelinated axon fraction (ρ = -0.90), the bound pool fraction (ƒ) versus the myelin sheath fraction (ρ = 0.93), and axial diffusivity versus the non-myelinated cell fraction (ρ = 0.92). Using Pearson's correlation coefficient, ƒ was strongly correlated to the myelin sheath fraction (r = 0.98) with a linear equation predicting myelin content (5.37ƒ - 0.25). Of the calculated MRI metrics, ƒ was the strongest indicator of myelin content, while longitudinal relaxation rates and diffusivity measurements were the strongest indicators of changes in tissue structure.

In-vivo multi-exponential T2, magnetization transfer and quantitative histology in a rat model of intramyelinic edema

Two MRI methods, multi-exponential analysis of transverse relaxation (MET₂) and quantitative magnetization transfer (qMT), were used along with quantitative evaluation of histology in a study of intra-myelinic edema in rat spinal white matter. The results showed a strong linear correlation between a distinct long-T₂ signal from MET₂ analysis and the edema water volume fraction as measured by histology, although this analysis overestimated the edema water content by ≈ 100% relative to quantitative histological measurements. This overestimation was reasoned to result from the effects of inter-compartmental water exchange on observed transverse relaxation. Commonly studied MRI markers for myelin, the myelin water fraction (from MET₂ analysis) and the macromolecular pool size ratio (from qMT analysis) produced results that could not be explained purely by changes in myelin content. The results demonstrate the potential for MET₂ analysis as well as the limits of putative myelin markers for characterizing white matter abnormalities involving intra-myelinic edema.

•

We studied a rat model of intra-myelinic edema induced by hexachlorophene ingestion.

•

We used multi-exponential T₂ (MET₂) and quantitative magnetization transfer MRI.

•

Histology was quantitatively evaluated to measure edema volume and myelin content.

•

MET₂ provides a measure that correlates but overestimates with edema volume fraction.

•

MET₂ measure of edema is affected by microscopic water dynamics.

MR susceptibility imaging

This work reviews recent developments in the use of magnetic susceptibility contrast for human MRI, with a focus on the study of brain anatomy. The increase in susceptibility contrast with modern high field scanners has led to novel applications and insights into the sources and mechanism contributing to this contrast in brain tissues. Dedicated experiments have demonstrated that in most of healthy brain, iron and myelin dominate tissue susceptibility variations, although their relative contribution varies substantially. Local variations in these compounds can affect both amplitude and frequency of the MRI signal. In white matter, the myelin sheath introduces an anisotropic susceptibility that has distinct effects on the water compartments inside the axons, between the myelin sheath, and the axonal space, and renders their signals dependent on the angle between the axon and the magnetic field. This offers opportunities to derive tissue properties specific to these cellular compartments.

Micro-compartment specific T2* relaxation in the brain

MRI at high field can be sensitized to the magnetic properties of tissues, which introduces a signal dependence on the orientation of white matter (WM) fiber bundles relative to the magnetic field. In addition, study of the NMR relaxation properties of this signal has indicated contributions from compartmentalized water environments inside and outside the myelin sheath that may be separable. Here we further investigated the effects of water compartmentalization on the MRI signal with the goal of extracting compartment-specific information. By comparing MRI measurements of human and marmoset brain at 7T with magnetic field modeling, we show that: (1) water between the myelin lipid bilayers, in the axonal, and in the interstitial space each experience characteristic magnetic field effects that depend on fiber orientation (2) these field effects result in characteristic relaxation properties and frequency shifts for these compartments; and (3) compartmental contributions may be separated by multi-component fitting of the MRI signal relaxation (i.e. decay) curve. We further show the potential application of these findings to the direct mapping of myelin content and assessment of WM fiber integrity with high field MRI.

Characterizing inter-compartmental water exchange in myelinated tissue using relaxation exchange spectroscopy

PURPOSE

To investigate inter-compartmental water exchange in two model myelinated tissues ex vivo using relaxation exchange spectroscopy.

METHODS

Building upon a previously developed theoretical framework, a three-compartment (myelin, intra-axonal, and extra-axonal water) model of the inversion-recovery prepared relaxation exchange spectroscopy signal was applied in excised rat optic nerve and frog sciatic nerve samples to estimate the water residence time constants in myelin (τmyelin ).

RESULTS

In the rat optic nerve samples, τmyelin = 138 ± 15 ms (mean ± standard deviation) was estimated. In sciatic nerve, which possesses thicker myelin sheaths than optic nerve, a much longer τmyelin = 2046 ± 140 ms was observed.

CONCLUSION

Consistent with previous studies in rat spinal cord, the extrapolation of exchange rates in optic nerve to in vivo conditions indicates that τmyelin < 100 ms. This suggests that there is a significant effect of inter-compartmental water exchange on the transverse relaxation of water protons in white matter. The much longer τmyelin values in sciatic nerve supports the postulate that the inter-compartmental water exchange rate is mediated by myelin thickness. Together, these findings point to the potential for MRI methods to probe variations in myelin thickness in white matter.

Segmentation of magnetization transfer ratio lesions for longitudinal analysis of demyelination and remyelination in multiple sclerosis

We demonstrate a new technique to quantify longitudinal changes in magnetization transfer ratio (MTR) magnetic resonance imaging (MRI). These changes are indicative of demyelination and remyelination. This technique comprises a definition of ΔMTR lesions, which are identified directly from the MTR images, and an automatic procedure for segmenting these lesions. We used this technique to analyze MTR changes in lesions of subjects with rapidly progressing multiple sclerosis before and after treatment with immunoablation and autologous stem cell transplant. Subjects who experienced clinical improvement after treatment showed significantly improved MTR recovery in lesions that were recovering during treatment (p<0.0001) while those who were clinically stable after treatment showed significantly poorer MTR recovery (p=0.002). The statistical power of this technique to detect treatment effects on MTR recovery was shown to be considerably better than previous methods. These results suggest that longitudinal measurements of MTR in ΔMTR lesions may be an important technique for the assessment of treatment effects on remyelination in clinical trials.

Fiber orientation-dependent white matter contrast in gradient echo MRI

Recent studies have shown that there is a direct link between the orientation of the nerve fibers in white matter (WM) and the contrast observed in magnitude and phase images acquired using gradient echo MRI. Understanding the origin of this link is of great interest because it could offer access to a new diagnostic tool for investigating tissue microstructure. Since it has been suggested that myelin is the dominant source of this contrast, creating an accurate model for characterizing the effect of the myelin sheath on the evolution of the NMR signal is an essential step toward fully understanding WM contrast. In this study, we show by comparison of the results of simulations and experiments carried out on human subjects at 7T, that the magnitude and phase of signals acquired from WM in vivo can be accurately characterized by (i) modeling the myelin sheath as a hollow cylinder composed of material having an anisotropic magnetic susceptibility that is described by a tensor with a radially oriented principal axis, and (ii) adopting a two-pool model in which the water in the sheath has a reduced T(2) relaxation time and spin density relative to its surroundings, and also undergoes exchange. The accuracy and intrinsic simplicity of the hollow cylinder model provides a versatile framework for future exploitation of the effect of WM microstructure on gradient echo contrast in clinical MRI.

On the inherent precision of mcDESPOT

A statistical analysis of the mcDESPOT protocol for characterizing two exchanging water proton pools--a seven-dimensional problem that fits to multiple flip angle measurements of both spoiled and refocused gradient echoes--is presented. Theoretical calculations of the Cramér-Rao lower bounds of the variance of fitted model parameters were made using a variety of model system parameters, meant to mimic those expected in human white matter. The results, validated by Monte Carlo simulations, indicated that mcDESPOT signals acquired at feasibly attainable signal-to-noise ratios cannot provide parameter estimates with useful levels of precision. Precision can be greatly improved by constraining solutions with a priori model information, although this will generally lead to biased parameter estimates with less specificity. These results indicate that previous, apparently successful applications of mcDESPOT to human white matter may have used data fitting methods that implicitly constrained parameter solutions, or that the two-pool model of white matter may not be sufficient to describe the observed water proton signal in mcDESPOT acquisitions. In either case, mcDESPOT-derived estimates of two-pool model parameters cannot yet be unambiguously related to specific tissue characteristics.