Multiexponential T2, magnetization transfer, and quantitative histology in white matter tracts of rat spinal cord

Dual Contrast - Magnetic Resonance Fingerprinting (DC-MRF): A Platform for Simultaneous Quantification of Multiple MRI Contrast Agents

Injectable Magnetic Resonance Imaging (MRI) contrast agents have been widely used to provide critical assessments of disease for both clinical and basic science imaging research studies. The scope of available MRI contrast agents has expanded over the years with the emergence of molecular imaging contrast agents specifically targeted to biological markers. Unfortunately, synergistic application of more than a single molecular contrast agent has been limited by MRI's ability to only dynamically measure a single agent at a time. In this study, a new Dual Contrast - Magnetic Resonance Fingerprinting (DC - MRF) methodology is described that can detect and independently quantify the local concentration of multiple MRI contrast agents following simultaneous administration. This "multi-color" MRI methodology provides the opportunity to monitor multiple molecular species simultaneously and provides a practical, quantitative imaging framework for the eventual clinical translation of molecular imaging contrast agents.

Lysophosphatidyl Choline Induced Demyelination in Rat Probed by Relaxation along a Fictitious Field in High Rank Rotating Frame

In this work a new MRI modality entitled Relaxation Along a Fictitious Field in the rotating frame of rank 4 (RAFF4) was evaluated in its ability to detect lower myelin content in lysophosphatidyl choline (LPC)-induced demyelinating lesions. The lesions were induced in two areas of the rat brain with either uniform or complex fiber orientations, i.e., in the corpus callosum (cc) and dorsal tegmental tract (dtg), respectively. RAFF4 showed excellent ability to detect demyelinated lesions and good correlation with myelin content in both brain areas. In comparison, diffusion tensor imaging metrices, fractional anisotropy, mean diffusivity and axonal and radial diffusivity, and magnetization transfer (MT) metrices, longitudinal relaxation during off-resonance irradiation and MT ratio, either failed to detect demyelination in dtg or showed lower correlation with myelin density quantified from gold chloride stained histological sections. Good specifity of RAFF4 to myelin was confirmed by its low correlation with cell density assesed from Nissl stained sections as well as its lack of sensitivity to pH changes in the physiological range as tested in heat denaturated bovine serum albumin phantoms. The excellent ability of RAFF4 to detect myelin content and its insensitivity to fiber orientation distribution, gliosis and pH, together with low specific absorption rate, demonstrates the promise of rotating frame of rank n (RAFFn) as a valuable MRI technique for non-invasive imaging of demyelinating lesions.

Evaluating g-ratio weighted changes in the corpus callosum as a function of age and sex

Recent years have seen a growing interest in relating MRI measurements to the structural-biophysical properties of white matter fibers. The fiber g-ratio, defined as the ratio between the inner and outer radii of the axon myelin sheath, is an important structural property of white matter, affecting signal conduction. Recently proposed modeling methods that use a combination of quantitative-MRI signals, enable a measurement of the fiber g-ratio in vivo. Here we use an MRI-based g-ratio estimation to observe the variance of the g-ratio within the corpus callosum, and evaluate sex and age related differences. To estimate the g-ratio we used a model (Stikov et al., 2011; Duval et al., 2017) based on two different WM microstructure parameters: the relative amounts of myelin (myelin volume fraction, MVF) and fibers (fiber volume fraction, FVF) in a voxel. We derived the FVF from the fractional anisotropy (FA), and estimated the MVF by using the lipid and macromolecular tissue volume (MTV), calculated from the proton density (Mezer et al., 2013). In comparison to other methods of estimating the MVF, MTV represents a stable parameter with a straightforward route of acquisition. To establish our model, we first compared histological MVF measurements (West et al., 2016) with the MRI derived MTV. We then implemented our model on a large database of 92 subjects (44 males), aged 7 to 81, in order to evaluate age and sex related changes within the corpus callosum. Our results show that the MTV provides a good estimation of MVF for calculating g-ratio, and produced values from the corpus callosum that correspond to those found in animals ex vivo and are close to the theoretical optimum, as well as to published in vivo data. Our results demonstrate that the MTV derived g-ratio provides a simple and reliable in vivo g-ratio-weighted (GR*) measurement in humans. In agreement with theoretical predictions, and unlike other tissue parameters measured with MRI, the g-ratio estimations were found to be relatively stable with age, and we found no support for a significant sexual dimorphism with age.

Effect of myelin water exchange on DTI-derived parameters in diffusion MRI: Elucidation of TE dependence

PURPOSE

Water exchange exists between different neuronal compartments of brain tissue but is often ignored in most diffusion models. The goal of the current study was to demonstrate the dependence of diffusion measurements on echo time (TE) in the human brain and to investigate the underlying effects of myelin water exchange.

METHODS

Five healthy subjects were examined with single-shot pulsed-gradient spin-echo echo-planar imaging with fixed duration (δ) and separation (Δ) of diffusion gradient pulses and a set of varying TEs. The effects of water exchange and intrinsic T₂ difference in cellular environments were investigated with Monte Carlo simulations.

RESULTS

Both in vivo measurements and simulations showed that fractional anisotropy (FA) and axial diffusivity (AD) had positive correlations with TE, while radial diffusivity (RD) showed a negative correlation, which is consistent with a previous study. The simulation results further indicated the sensitivity of TE dependence to the change of g-ratio.

CONCLUSION

The exchange between myelin and intra/extra-axonal water pools often plays a non-negligible role in the observed TE dependence of diffusion parameters, which may accompany or alter the effect of intrinsic T₂ in causing such dependence. The TE dependence may potentially serve as a biomarker for demyelination processes (e.g., in multiple sclerosis and Alzheimer's disease). Magn Reson Med 79:1650-1660, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Quantitative in vivo T2 mapping using fast spin echo techniques - A linear correction procedure

A method is presented for correcting the effects of stimulated and indirect echoes on quantitative T2 mapping data acquired with multiple spin echo techniques, such as turbo spin echo. In contrast to similar correction techniques proposed in the literature, the method does not require a priori knowledge of the radio frequency (RF) pulse profiles. In a first step, for the T2 mapping protocol under investigation, signal decay curves S(TE) are simulated for a range of different RF pulse profiles. The actual signal decay S(TE) is then measured on a phantom with known T2, so the approximate RF pulse profiles can be derived via comparison with the simulated decay curves. In a second step, with the RF pulses obtained from step one, signal decay curves S(TE) are simulated for different T2 values and fitted mono-exponentially, thus allowing to deduce the relationship between true T2 and the apparent T2 (T2app) values. Results show that this relationship is approximately linear, allowing for a direct correction of T2app maps. If the amplitude of the transmitted RF field (B1) does not exceed the nominal value by more than 10%, it is shown that a B1-independent correction of T2app maps yields sufficiently accurate results for T2. A B1-dependent version is also presented. The method is tested in vitro on a phantom with different T2 values and in vivo on healthy subjects.

Optimization of inhomogeneous magnetization transfer (ihMT) MRI contrast for preclinical studies using dipolar relaxation time (T1D ) filtering

A pulsed inhomogeneous magnetization transfer (ihMT)-prepared fast imaging sequence was implemented at 11.75 T for preclinical studies on mouse central nervous system. A strategy based on filtering the ihMT signal originating from short dipolar relaxation time (T_1D ) components is proposed. It involves increasing the repetition time of consecutive radiofrequency (RF) pulses of the dual saturation and allows improved signal specificity for long T_1D myelinated structures. Furthermore, frequency offset, power and timing saturation parameters were adjusted to optimize the ihMT sensitivity. The optimization of the ihMT sensitivity, whilst preserving the strong specificity for the long T_1D component of myelinated tissues, allowed measurements of ihMT ratios on the order of 4-5% in white matter (WM), 2.5% in gray matter (GM) and 1-1.3% in muscle. This led to high relative ihMT contrasts between myelinated tissues and others (~3-4 between WM and muscle, and ≥2 between GM and muscle). Conversely, higher ihMT ratios (~6-7% in WM) could be obtained using minimal T_1D filtering achieved with short saturation pulse repetition time or cosine-modulated pulses for the dual-frequency saturation. This study represents a first stage in the process of validating ihMT as a myelin biomarker by providing optimized ihMT preclinical sequences, directly transposable and applicable to other preclinical magnetic fields and scanners. Finally, ihMT ratios measured in various central nervous system areas are provided for future reference.

Impact of transcytolemmal water exchange on estimates of tissue microstructural properties derived from diffusion MRI

PURPOSE

To investigate the influence of transcytolemmal water exchange on estimates of tissue microstructural parameters derived from diffusion MRI using conventional PGSE and IMPULSED methods.

METHODS

Computer simulations were performed to incorporate a broad range of intracellular water life times τin (50-∞ ms), cell diameters d (5-15 μm), and intrinsic diffusion coefficient Din (0.6-2 μm2 /ms) for different values of signal-to-noise ratio (SNR) (10 to 50). For experiments, murine erythroleukemia (MEL) cancer cells were cultured and treated with saponin to selectively change cell membrane permeability. All fitted microstructural parameters from simulations and experiments in vitro were compared with ground-truth values.

RESULTS

Simulations showed that, for both PGSE and IMPULSED methods, cell diameter d can be reliably fit with sufficient SNR (≥ 50), whereas intracellular volume fraction fin is intrinsically underestimated due to transcytolemmal water exchange. Din can be reliably fit only with sufficient SNR and using the IMPULSED method with short diffusion times. These results were confirmed with those obtained in the cell culture experiments in vitro.

CONCLUSION

For the sequences and models considered in this study, transcytolemmal water exchange has minor effects on the fittings of d and Din with physiologically relevant membrane permeabilities if the SNR is sufficient (> 50), but fin is intrinsically underestimated. Magn Reson Med 77:2239-2249, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Histological validation of fast macromolecular proton fraction mapping as a quantitative myelin imaging method in the cuprizone demyelination model

Cuprizone-induced demyelination in mice is a frequently used model in preclinical multiple sclerosis research. A recent quantitative clinically-targeted MRI method, fast macromolecular proton fraction (MPF) mapping demonstrated a promise as a myelin biomarker in human and animal studies with a particular advantage of sensitivity to both white matter (WM) and gray matter (GM) demyelination. This study aimed to histologically validate the capability of MPF mapping to quantify myelin loss in brain tissues using the cuprizone demyelination model. Whole-brain MPF maps were obtained in vivo on an 11.7T animal MRI scanner from 7 cuprizone-treated and 7 control С57BL/6 mice using the fast single-point synthetic-reference method. Brain sections were histologically stained with Luxol Fast Blue (LFB) for myelin quantification. Significant (p < 0.05) demyelination in cuprizone-treated animals was found according to both LFB staining and MPF in all anatomical structures (corpus callosum, anterior commissure, internal capsule, thalamus, caudoputamen, and cortex). MPF strongly correlated with quantitative histology in all animals (r = 0.95, p < 0.001) as well as in treatment and control groups taken separately (r = 0.96, p = 0.002 and r = 0.93, p = 0.007, respectively). Close agreement between histological myelin staining and MPF suggests that fast MPF mapping enables robust and accurate quantitative assessment of demyelination in both WM and GM.

Modeling white matter microstructure

Quantitative magnetic resonance imaging can be combined with advanced biophysical models to measure microstructural features of white matter. Non-invasive microstructural imaging has the potential to revolutionize neuroscience, and acquiring these measures in clinically feasible times would greatly improve patient monitoring and clinical studies of drug efficacy. However, a good understanding of microstructural imaging techniques is essential to set realistic expectations and to prevent over-interpretation of results. This review explains the methodology behind microstructural modeling and imaging, and gives an overview of the breakthroughs and challenges associated with it.

Magnetic resonance imaging assessment of degenerative cervical myelopathy: a review of structural changes and measurement techniques

Degenerative cervical myelopathy encompasses a spectrum of age-related structural changes of the cervical spine that result in static and dynamic injury to the spinal cord and collectively represent the most common cause of myelopathy in adults. Although cervical myelopathy is determined clinically, the diagnosis requires confirmation via imaging, and MRI is the preferred modality. Because of the heterogeneity of the condition and evolution of MRI technology, multiple techniques have been developed over the years in an attempt to quantify the degree of baseline severity and potential for neurological recovery. In this review, these techniques are categorized anatomically into those that focus on bone, ligaments, discs, and the spinal cord. In addition, measurements for the cervical spine canal size and sagittal alignment are also described briefly. These tools have resulted collectively in the identification of numerous useful parameters. However, the development of multiple techniques for assessing the same feature, such as cord compression, has also resulted in a number of challenges, including introducing ambiguity in terms of which methods to use and hindering effective comparisons of analysis in the literature. In addition, newer techniques that use advanced MRI are emerging and providing exciting new tools for assessing the spinal cord in patients with degenerative cervical myelopathy.

New rapid, accurate T2 quantification detects pathology in normal-appearing brain regions of relapsing-remitting MS patients

Introduction

Quantitative T₂ mapping may provide an objective biomarker for occult nervous tissue pathology in relapsing-remitting multiple sclerosis (RRMS). We applied a novel echo modulation curve (EMC) algorithm to identify T₂ changes in normal-appearing brain regions of subjects with RRMS (N = 27) compared to age-matched controls (N = 38).

Methods

The EMC algorithm uses Bloch simulations to model T₂ decay curves in multi-spin-echo MRI sequences, independent of scanner, and scan-settings. T₂ values were extracted from normal-appearing white and gray matter brain regions using both expert manual regions-of-interest and user-independent FreeSurfer segmentation.

Results

Compared to conventional exponential T₂ modeling, EMC fitting provided more accurate estimations of T₂ with less variance across scans, MRI systems, and healthy individuals. Thalamic T₂ was increased 8.5% in RRMS subjects (p < 0.001) and could be used to discriminate RRMS from healthy controls well (AUC = 0.913). Manual segmentation detected both statistically significant increases (corpus callosum & temporal stem) and decreases (posterior limb internal capsule) in T₂ associated with RRMS diagnosis (all p < 0.05). In healthy controls, we also observed statistically significant T₂ differences for different white and gray matter structures.

Conclusions

The EMC algorithm precisely characterizes T₂ values, and is able to detect subtle T₂ changes in normal-appearing brain regions of RRMS patients. These presumably capture both axon and myelin changes from inflammation and neurodegeneration. Further, T₂ variations between different brain regions of healthy controls may correlate with distinct nervous tissue environments that differ from one another at a mesoscopic length-scale.

•

EMC technique provides accurate and scanner-invariant T₂ mapping in MS subjects.

•

Thalamus T₂ differences distinguish relapsing-remitting MS subjects from controls.

•

Normal-appearing brain regions demonstrate T₂ changes in MS patients compared to controls.

•

T₂ values reflect anatomic and function-specific differences in healthy controls.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

g-Ratio weighted imaging of the human spinal cord in vivo

The fiber g-ratio is defined as the ratio of the inner to the outer diameter of the myelin sheath. This ratio provides a measure of the myelin thickness that complements axon morphology (diameter and density) for assessment of demyelination in diseases such as multiple sclerosis. Previous work has shown that an aggregate g-ratio map can be computed using a formula that combines axon and myelin density measured with quantitative MRI. In this work, we computed g-ratio weighted maps in the cervical spinal cord of nine healthy subjects. We utilized the 300mT/m gradients from the CONNECTOM scanner to estimate the fraction of restricted water (fr) with high accuracy, using the CHARMED model. Myelin density was estimated using the lipid and macromolecular tissue volume (MTV) method, derived from normalized proton density (PD) mapping. The variability across spinal level, laterality and subject were assessed using a three-way ANOVA. The average g-ratio value obtained in the white matter was 0.76+/-0.03, consistent with previous histology work. Coefficients of variation of fr and MTV were respectively 4.3% and 13.7%. fr and myelin density were significantly different across spinal tracts (p=3×10-7 and 0.004 respectively) and were positively correlated in the white matter (r=0.42), suggesting shared microstructural information. The aggregate g-ratio did not show significant differences across tracts (p=0.6). This study suggests that fr and myelin density can be measured in vivo with high precision and that they can be combined to produce a g-ratio-weighted map robust to free water pool contamination from cerebrospinal fluid or veins. Potential applications include the study of early demyelination in multiple sclerosis, and the quantitative assessment of remyelination drugs.

Test-retest reliability and concurrent validity of in vivo myelin content indices: Myelin water fraction and calibrated T1 w/T2 w image ratio

In an age-heterogeneous sample of healthy adults, we examined test-retest reliability (with and without participant repositioning) of two popular MRI methods of estimating myelin content: modeling the short spin-spin (T₂ ) relaxation component of multi-echo imaging data and computing the ratio of T₁ -weighted and T₂ -weighted images (T₁ w/T₂ w). Taking the myelin water fraction (MWF) index of myelin content derived from the multi-component T₂ relaxation data as a standard, we evaluate the concurrent and differential validity of T₁ w/T₂ w ratio images. The results revealed high reliability of MWF and T₁ w/T₂ w ratio. However, we found significant correlations of low to moderate magnitude between MWF and the T₁ w/T₂ w ratio in only two of six examined regions of the cerebral white matter. Notably, significant correlations of the same or greater magnitude were observed for T₁ w/T₂ w ratio and the intermediate T₂ relaxation time constant, which is believed to reflect differences in the mobility of water between the intracellular and extracellular compartments. We conclude that although both methods are highly reliable and thus well-suited for longitudinal studies, T₁ w/T₂ w ratio has low criterion validity and may be not an optimal index of subcortical myelin content. Hum Brain Mapp 38:1780-1790, 2017. © 2017 Wiley Periodicals, Inc.

Magnetic Resonance of Myelin Water: An in vivo Marker for Myelin

Myelin is critical for healthy brain function. An accurate in vivo measure of myelin content has important implications for understanding brain plasticity and neurodegenerative diseases. Myelin water imaging is a magnetic resonance imaging method which can be used to visualize myelination in the brain and spinal cord in vivo. This review presents an overview of myelin water imaging data acquisition and analysis, post-mortem validation work, findings in both animal and human studies and a brief discussion about other MR techniques purported to provide in vivo myelin content. Multi-echo T₂ relaxation approaches continue to undergo development and whole-brain imaging time now takes less than 10 minutes; the standard analysis method for this type of data acquisition is a non-negative least squares approach. Alternate methods including the multi-flip angle gradient echo mcDESPOT are also being used for myelin water imaging. Histological validation studies in animal and human brain and spinal cord tissue demonstrate high specificity of myelin water imaging for myelin. Potential confounding factors for in vivo myelin water fraction measurement include the presence of myelin debris and magnetization exchange processes. Myelin water imaging has successfully been used to study animal models of injury, applied in healthy human controls and can be used to assess damage and injury in conditions such as multiple sclerosis, neuromyelitis optica, schizophrenia, phenylketonuria, neurofibromatosis, niemann pick’s disease, stroke and concussion. Other quantitative magnetic resonance approaches that are sensitive to, but not specific for, myelin exist including magnetization transfer, diffusion tensor imaging and T₁ weighted imaging.

Characterizing axonal myelination within the healthy population: a tract-by-tract mapping of effects of age and gender on the fiber g-ratio

The g-ratio, equal to the ratio of the inner-to-outer diameter of a myelinated axon, is associated with the speed of conduction, and thus reflects axonal function and integrity. It is now possible to estimate an “aggregate” g-ratio in vivo using MRI. The aim of this study was to assess the variation of the MRI-derived fiber g-ratio in the brain of healthy individuals, and to characterize its variation across the lifespan. Thirty-eight healthy participants, aged between 20 and 76, were recruited. Whole-brain g-ratio maps were computed and analyzed voxel-wise. Median tract g-ratio values were also extracted. No significant effect of gender was found, whereas age was found to be significantly associated with the g-ratio within the white matter. The tract-specific analysis showed this relationship to follow a nearly-linear increase, although the slope appears to slow down slightly after the 6th decade of life. The most likely interpretation is a subtle but consistent reduction in myelin throughout adulthood, with the density of axons beginning to decrease between the 4th and 5th decade.

AxonSeg: Open Source Software for Axon and Myelin Segmentation and Morphometric Analysis

Segmenting axon and myelin from microscopic images is relevant for studying the peripheral and central nervous system and for validating new MRI techniques that aim at quantifying tissue microstructure. While several software packages have been proposed, their interface is sometimes limited and/or they are designed to work with a specific modality (e.g., scanning electron microscopy (SEM) only). Here we introduce AxonSeg, which allows to perform automatic axon and myelin segmentation on histology images, and to extract relevant morphometric information, such as axon diameter distribution, axon density and the myelin g-ratio. AxonSeg includes a simple and intuitive MATLAB-based graphical user interface (GUI) and can easily be adapted to a variety of imaging modalities. The main steps of AxonSeg consist of: (i) image pre-processing; (ii) pre-segmentation of axons over a cropped image and discriminant analysis (DA) to select the best parameters based on axon shape and intensity information; (iii) automatic axon and myelin segmentation over the full image; and (iv) atlas-based statistics to extract morphometric information. Segmentation results from standard optical microscopy (OM), SEM and coherent anti-Stokes Raman scattering (CARS) microscopy are presented, along with validation against manual segmentations. Being fully-automatic after a quick manual intervention on a cropped image, we believe AxonSeg will be useful to researchers interested in large throughput histology. AxonSeg is open source and freely available at: https://github.com/neuropoly/axonseg.

Assessment of the myelin water fraction in rodent spinal cord using T2-prepared ultrashort echo time MRI

OBJECTIVE

Multi-component T2 relaxation allows for assessing the myelin water fraction in nervous tissue, providing a surrogate marker for demyelination. The assessment of the number and distribution of different T2 components for devising exact models of tissue relaxation has been limited by T2 sampling with conventional MR methods.

MATERIALS AND METHODS

A T2-prepared UTE sequence was used to assess multicomponent T2 relaxation at 9.4 T of fixed mouse and rat spinal cord samples and of mouse spinal cord in vivo. For in vivo scans, a cryogenically cooled probe allowed for 78-µm resolution in 1-mm slices. Voxel-wise non-negative least square analysis was used to assess the number of myelin water-associated T2 components.

RESULTS

More than one myelin water-associated T2 component was detected in only 12 % of analyzed voxels in rat spinal cords and 6 % in mouse spinal cords, both in vivo and in vitro. However, myelin water-associated T2 values of individual voxels varied between 0.1 and 20 ms. While in fixed samples almost no components below 1 ms were identified, in vivo, these contributed 14 % of the T2 spectrum. No significant differences in MWF were observed in mouse spinal cord in vivo versus ex vivo measurements.

CONCLUSION

Voxel-wise analysis methods using relaxation models with one myelin water-associated T2 component are appropriate for assessing myelin content of nervous tissue.

Simulations on the influence of myelin water in diffusion-weighted imaging

While myelinated axons present an important barrier to water diffusion, many models used to interpret DWI signal neglect other potential influences of myelin. In this work, Monte Carlo simulations were used to test the sensitivity of DWI results to the diffusive properties of water within myelin. Within these simulations, the apparent diffusion coefficient (D app) varied slowly over several orders of magnitude of the coefficient of myelin water diffusion (D m), but exhibited important differences compared to D app values simulated that neglect D m (=0). Compared to D app, the apparent diffusion kurtosis (K app) was generally more sensitive to D m. Simulations also tested the sensitivity of D app and K app to the amount of myelin present. Unique variations in D app and K app caused by differences in the myelin volume fraction were diminished when myelin water diffusion was included. Also, expected trends in D app and K app with experimental echo time were reduced or inverted when accounting for myelin water diffusion, and these reduced/inverted trends were seen experimentally in ex vivo rat brain DWI experiments. In general, myelin water has the potential to subtly influence DWI results and bias models of DWI that neglect these components of white matter.

In vivo imaging of neuroinflammation in schizophrenia

In recent years evidence has accumulated to suggest that neuroinflammation might be an early pathology of schizophrenia that later leads to neurodegeneration, yet the exact role in the etiology, as well as the source of neuroinflammation, are still not known. The hypothesis of neuroinflammation involvement in schizophrenia is quickly gaining popularity, and thus it is imperative that we have reliable and reproducible tools and measures that are both sensitive, and, most importantly, specific to neuroinflammation. The development and use of appropriate human in vivo imaging methods can help in our understanding of the location and extent of neuroinflammation in different stages of the disorder, its natural time-course, and its relation to neurodegeneration. Thus far, there is little in vivo evidence derived from neuroimaging methods. This is likely the case because the methods that are specific and sensitive to neuroinflammation are relatively new or only just being developed. This paper provides a methodological review of both existing and emerging positron emission tomography and magnetic resonance imaging techniques that identify and characterize neuroinflammation. We describe \how these methods have been used in schizophrenia research. We also outline the shortcomings of existing methods, and we highlight promising future techniques that will likely improve state-of-the-art neuroimaging as a more refined approach for investigating neuroinflammation in schizophrenia.

Reprint of "Hypomyelinating disorders: An MRI approach

In recent years, the concept of hypomyelinating disorders has been proposed as a group of disorders with varying systemic manifestations that are identified by MR findings of absence or near absence of the T2 hypointensity that develops in white matter as a result of myelination. Initially proposed as a separate group because they were the largest single category of undiagnosed leukodystrophies, their separation as a distinct group that can be recognized by looking for a specific MRI feature has resulted in a marked increase in their diagnosis and a better understanding of the different causes of hypomyelination. This review will discuss the clinical presentations, imaging findings on standard MRI, and new MRI-related techniques that allow a better understanding of these disorders and proposed methods for quantifying the myelination as a potential means of assessing disease course and the effects of proposed treatments. Disorders with hypomyelination of white matter, or hypomyelinating disorders (HMDs), represent the single largest category among undiagnosed genetic leukoencephalopathies (Schiffmann and van der Knaap, 2009; Steenweg et al., 2010). This group of inborn errors of metabolism is characterized by a magnetic resonance imaging (MRI) appearance of reduced or absent myelin development: delay in the development of T2 hypointensity and, often, T1 hyperintensity in the white matter of the brain. The concept of hypomyelination was first conceptualized by (Schiffmann and van der Knaap, 2009; Steenweg et al., 2010; Schiffmann et al., 1994) in a series of papers that showed that these MRI characteristics were easily recognized, were different from the MRI characteristics of dysmyelinating and demyelinating disorders, and that the combination of these imaging findings with specific other clinical and imaging features could be used to make diagnoses with some confidence. In this manuscript, we will discuss the physiologic and genetic bases of hypomyelinating disorders, as well as their classification, clinical manifestations and imaging characteristics.

Mapping an index of the myelin g-ratio in infants using magnetic resonance imaging

Optimal myelination of neuronal axons is essential for effective brain and cognitive function. The ratio of the axon diameter to the outer fiber diameter, known as the g-ratio, is a reliable measure to assess axonal myelination and is an important index reflecting the efficiency and maximal conduction velocity of white matter pathways. Although advanced neuroimaging techniques including multicomponent relaxometry (MCR) and diffusion tensor imaging afford insight into the microstructural characteristics of brain tissue, by themselves they do not allow direct analysis of the myelin g-ratio. Here, we show that by combining myelin content information (obtained with mcDESPOT MCR) with neurite density information (obtained through NODDI diffusion imaging) an index of the myelin g-ratio may be estimated. Using this framework, we present the first quantitative study of myelin g-ratio index changes across childhood, examining 18 typically developing children 3months to 7.5years of age. We report a spatio-temporal pattern of maturation that is consistent with histological and developmental MRI studies, as well as theoretical studies of the myelin g-ratio. This work represents the first ever in vivo visualization of the evolution of white matter g-ratio indices throughout early childhood.

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.

PGSE, OGSE, and sensitivity to axon diameter in diffusion MRI: Insight from a simulation study

PURPOSE

To identify optimal pulsed gradient spin-echo (PGSE) and oscillating gradient spin-echo (OGSE) sequence settings for maximizing sensitivity to axon diameter in idealized and practical conditions.

METHODS

Simulations on a simple two-compartment white matter model (with nonpermeable cylinders) are used to investigate a wide space of clinically plausible PGSE and OGSE sequence parameters with trapezoidal diffusion gradient waveforms. Signal sensitivity is measured as a derivative of the signal with respect to axon diameter. Models of parallel and dispersed fibers are investigated separately to represent idealized and practical conditions.

RESULTS

Simulations show that, for the simple case of gradients perfectly perpendicular to straight parallel fibers, PGSE always gives maximum sensitivity. However, in real-world scenarios where fibers have unknown and dispersed orientation, low-frequency OGSE provides higher sensitivity. Maximum sensitivity results show that on current clinical scanners (Gmax  = 60 mT/m, signal to noise ratio (SNR) = 20) axon diameters below 6 µm are indistinguishable from zero. Scanners with stronger gradient systems such as the Massachusetts General Hospital (MGH) Connectom scanner (Gmax  = 300 mT/m) can extend this sensitivity limit down to 2-3 µm, probing a much greater proportion of the underlying axon diameter distribution.

CONCLUSION

Low-frequency OGSE provides additional sensitivity to PGSE in practical situations. OGSE is particularly advantageous for systems with high performance gradients.

Hypomyelinating disorders: An MRI approach

In recent years, the concept of hypomyelinating disorders has been proposed as a group of disorders with varying systemic manifestations that are identified by MR findings of absence or near absence of the T2 hypointensity that develops in white matter as a result of myelination. Initially proposed as a separate group because they were the largest single category of undiagnosed leukodystrophies, their separation as a distinct group that can be recognized by looking for a specific MRI feature has resulted in a marked increase in their diagnosis and a better understanding of the different causes of hypomyelination. This review will discuss the clinical presentations, imaging findings on standard MRI, and new MRI-related techniques that allow a better understanding of these disorders and proposed methods for quantifying the myelination as a potential means of assessing disease course and the effects of proposed treatments. Disorders with hypomyelination of white matter, or hypomyelinating disorders (HMDs), represent the single largest category among undiagnosed genetic leukoencephalopathies (Schiffmann and van der Knaap, 2009; Steenweg et al., 2010). This group of inborn errors of metabolism is characterized by a magnetic resonance imaging (MRI) appearance of reduced or absent myelin development: delay in the development of T2 hypointensity and, often, T1 hyperintensity in the white matter of the brain. The concept of hypomyelination was first conceptualized by (Schiffmann and van der Knaap, 2009; Steenweg et al., 2010; Schiffmann et al., 1994) in a series of papers that showed that these MRI characteristics were easily recognized, were different from the MRI characteristics of dysmyelinating and demyelinating disorders, and that the combination of these imaging findings with specific other clinical and imaging features could be used to make diagnoses with some confidence. In this manuscript, we will discuss the physiologic and genetic bases of hypomyelinating disorders, as well as their classification, clinical manifestations and imaging characteristics.

Compressed sensing CPMG with group-sparse reconstruction for myelin water imaging

PURPOSE

Myelin content is a marker for nervous system pathology and is quantifiable by myelin water imaging using multi-echo CPMG sequence, which is inherently slow. One way to accelerate the scan is to utilize compressed sensing. However, reconstructing the images piecemeal by standard compressed sensing methods is not the optimal solution, because it only exploits intraimage spatial redundancy. It does not recognize that the different T2 weighted images are scans of the same anatomical volume and hence correlated. The purpose of this work is to test the feasibility of compressed sensed CPMG with group-sparsity promoting optimization for myelin water imaging.

METHODS

Group-sparse reconstruction was performed at various simulated and actual undersampling factors for an electronic phantom, ex vivo rat spinal cord, and in vivo rat spinal cord. Normalized mean square error was used as the metric for comparison.

RESULTS

For both simulated undersampling and the actual undersampling, the method was found to minimally impact myelin water fraction map quality (normalized mean square error < 0.25) when acceleration factor was below two.

CONCLUSION

Compressed sensed CPMG with group-sparse reconstruction is useful for achieving a shorter scan time than traditionally possible.

Analysis and correction of biases in cross-relaxation MRI due to biexponential longitudinal relaxation

PURPOSE

Cross-relaxation imaging (CRI) is a family of quantitative magnetization transfer techniques that utilize images obtained with off-resonance saturation and longitudinal relaxation rate (R1) maps reconstructed by the variable flip angle (VFA) method. It was demonstrated recently that a significant bias in an apparent VFA R1 estimation occurs in macromolecule-rich tissues due to magnetization transfer (MT)-induced biexponential behavior of longitudinal relaxation of water protons. The purpose of this article is to characterize theoretically and experimentally the resulting bias in the CRI maps and propose methods to correct it.

THEORY

The modified CRI algorithm is proposed, which corrects for such biases and yields accurate parametric bound pool fraction f, cross-relaxation rate k, and R1 maps. Additionally, an analytical correction procedure is introduced to recalculate previously obtained parameter values.

RESULTS

The systematic errors due to unaccounted MT-induced biexponential relaxation can be characterized as an overestimation of R1, f, and k, with a relative bias comparable with the magnitude of f. The phantom and human in vivo experiments demonstrate that both proposed modified CRI and analytical correction approaches significantly improve the accuracy of the CRI method.

CONCLUSION

The accuracy of the CRI method can be considerably improved by taking into account the contribution of MT-induced biexponential longitudinal relaxation into variable flip angle R1 measurements.

Quantitative MRI in a non-surgical model of cervical spinal cord injury

Quantitative T2 (qT2), diffusion tensor imaging (DTI), and histology were used to investigate a cervical model of spinal cord injury (SCI) in the rat. While quantitative MRI can significantly increase the specificity in the presence of pathology, it must be validated for each type of injury or disease. In the case of traumatic SCI most models are difficult to image, either due to the location of the injury, or as a result of damage to surrounding tissues resulting from invasive surgical procedures. In this study a non-surgical cervical model of SCI, produced using a combination of focused ultrasound and microbubbles, was used to produce pathology similar to that seen in models of contusive and compressive injuries. qT2 and DTI were performed at 24 h and 1 and 2 weeks following injury, and compared with H&E and luxol fast blue histology. In the injured spinal cord, in addition to intra/extracellular (I/E) water and myelin water in white matter, qT2 revealed a large component with very short T2 of about 3 ms, which was highly correlated with the presence of hemorrhage in both gray and white matter at 24 h, and with the presence of hemosiderin in gray matter at 2 weeks following injury. The T2 of the I/E water peak was also elevated at 24 h in both gray and white matter, which was correlated with the presence of vacuolation/edema on histology. Cystic cavities were only seen at the 1 or 2 week timepoints, and were correlated with the presence of a water peak with T2  > 250 ms. No significant changes in diffusivity parameters were observed. Pathologies were often co-occurring, with opposite effects on the average T2 in a given voxel, reducing the visibility of injured tissue on standard T2 -weighted MR images.

In vivo mapping of human spinal cord microstructure at 300mT/m

The ability to characterize white matter microstructure non-invasively has important applications for the diagnosis and follow-up of several neurological diseases. There exists a family of diffusion MRI techniques, such as AxCaliber, that provide indices of axon microstructure, such as axon diameter and density. However, to obtain accurate measurements of axons with small diameters (<5μm), these techniques require strong gradients, i.e. an order of magnitude higher than the 40-80mT/m currently available in clinical systems. In this study we acquired AxCaliber diffusion data at a variety of different q-values and diffusion times in the spinal cord of five healthy subjects using a 300mT/m whole body gradient system. Acquisition and processing were optimized using state-of-the-art methods (e.g., 64-channel coil, template-based analysis). Results consistently show an average axon diameter of 4.5+/-1.1μm in the spinal cord white matter. Diameters ranged from 3.0μm (gracilis) to 5.9μm (spinocerebellar tracts). Values were similar across laterality (left-right), but statistically different across spinal cord pathways (p<10(-5)). The observed trends are similar to those observed in animal histology. This study shows, for the first time, in vivo mapping of axon diameter in the spinal cord at 300mT/m, thus creating opportunities for applications in spinal cord diseases.

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.

PGSE, OGSE, and sensitivity to axon diameter in diffusion MRI: Insight from a simulation study

Purpose

To identify optimal pulsed gradient spin‐echo (PGSE) and oscillating gradient spin‐echo (OGSE) sequence settings for maximizing sensitivity to axon diameter in idealized and practical conditions.

Methods

Simulations on a simple two‐compartment white matter model (with nonpermeable cylinders) are used to investigate a wide space of clinically plausible PGSE and OGSE sequence parameters with trapezoidal diffusion gradient waveforms. Signal sensitivity is measured as a derivative of the signal with respect to axon diameter. Models of parallel and dispersed fibers are investigated separately to represent idealized and practical conditions.

Results

Simulations show that, for the simple case of gradients perfectly perpendicular to straight parallel fibers, PGSE always gives maximum sensitivity. However, in real‐world scenarios where fibers have unknown and dispersed orientation, low‐frequency OGSE provides higher sensitivity. Maximum sensitivity results show that on current clinical scanners (G_max = 60 mT/m, signal to noise ratio (SNR) = 20) axon diameters below 6 µm are indistinguishable from zero. Scanners with stronger gradient systems such as the Massachusetts General Hospital (MGH) Connectom scanner (G_max = 300 mT/m) can extend this sensitivity limit down to 2–3 µm, probing a much greater proportion of the underlying axon diameter distribution.

Conclusion

Low‐frequency OGSE provides additional sensitivity to PGSE in practical situations. OGSE is particularly advantageous for systems with high performance gradients. Magn Reson Med 75:688–700, 2016. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

In vivo quantitative magnetization transfer imaging correlates with histology during de- and remyelination in cuprizone-treated mice

The pool size ratio measured by quantitative magnetization transfer MRI is hypothesized to closely reflect myelin density, but their relationship has so far been confirmed mostly in ex vivo conditions. We investigate the correspondence between this parameter measured in vivo at 7.0 T, with Black Gold II staining for myelin fibres, and with myelin basic protein and beta-tubulin immunofluorescence in a hybrid longitudinal study of C57BL/6 and SJL/J mice treated with cuprizone, a neurotoxicant causing relatively selective myelin loss followed by spontaneous remyelination upon treatment suspension. Our results confirm that pool size ratio measurements correlate with myelin content, with the correlation coefficient depending on strain and staining method, and demonstrate the in vivo applicability of this MRI technique to experimental mouse models of multiple sclerosis.

Microstructural origins of gadolinium-enhanced susceptibility contrast and anisotropy

PURPOSE

MR histology based on magnetic susceptibility can be used to visualize diamagnetic myelin (and its deterioration) in the central nervous system and is facilitated by the application of high magnetic field strengths and paramagnetic contrast agents. Characterizing the effect of these tools will aid in assessing white matter myelin content and microstructure.

METHODS

Image data from six gadolinium-perfused mouse brain specimens were acquired at 2.0, 7.0, and 9.4 Tesla. Magnetic susceptibility contrast was analyzed for its dependence on field strength, gadolinium concentration, and white matter fiber orientation. A model for this contrast is presented based on the three-pool model for white matter.

RESULTS

The specimen data illustrate that white-gray matter susceptibility contrast is field strength independent. White-gray matter contrast improves significantly as a function of gadolinium contrast agent in the tissue, i.e., white matter appears increasingly more diamagnetic relative to gray matter. The simulated data from the model suggest that susceptibility anisotropy of white matter fiber bundles increases nonlinearly as a function of gadolinium concentration due to contrast agent compartmentalization into the extracellular white matter water pool.

CONCLUSION

Using contrast agents in MR histology facilitates white-gray matter susceptibility contrast modulation and the probing of white matter microstructure and orientation.

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.

Removal of cerebrospinal fluid partial volume effects in quantitative magnetization transfer imaging using a three-pool model with nonexchanging water component

PURPOSE

Parameters of the two-pool model describing magnetization transfer (MT) in macromolecule-rich tissues may be significantly biased in partial volume (PV) voxels containing cerebrospinal fluid (CSF). The purpose of this study was to develop a quantitative MT (qMT) method that provides indices insensitive to CSF PV averaging.

THEORY AND METHODS

We propose a three-pool MT model, in which PV macro-compartment is modeled as an additional nonexchanging water pool. We demonstrate the feasibility of model parameter estimation from several MT-weighted spoiled gradient echo datasets. We validated the three-pool model in numerical, phantom, and in vivo studies.

RESULTS

PV averaging with the free water compartment reduces all qMT parameters, most significantly affecting macromolecular proton fraction (MPF) and cross-relaxation rate. Monte-Carlo simulations confirmed stability of the three-pool model fit. Unlike the standard two-pool model, the three-pool model qMT parameters were not affected by PV averaging in simulations and phantom studies. The three-pool model fit allowed CSF PV correction in brain PV voxels and resulted in good correlation with standard two-pool model parameters in non-PV voxels.

CONCLUSION

Quantitative MT imaging based on a three-pool model with a non-exchanging water component yields a set of CSF-insensitive qMT parameters, which may improve MPF-based assessment of myelination in structures strongly affected by CSF PV averaging such as brain gray matter.

A framework for accurate determination of the T₂ distribution from multiple echo magnitude MRI images

Measurement of the T2 distribution in tissues provides biologically relevant information about normal and abnormal microstructure and organization. Typically, the T2 distribution is obtained by fitting the magnitude MR images acquired by a multi-echo MRI pulse sequence using an inverse Laplace transform (ILT) algorithm. It is well known that the ideal magnitude MR signal follows a Rician distribution. Unfortunately, studies attempting to establish the validity and efficacy of the ILT algorithm assume that these input signals are Gaussian distributed. Violation of the normality (or Gaussian) assumption introduces unexpected artifacts, including spurious cerebrospinal fluid (CSF)-like long T2 components; bias of the true geometric mean T2 values and in the relative fractions of various components; and blurring of nearby T2 peaks in the T2 distribution. Here we apply and extend our previously proposed magnitude signal transformation framework to map noisy Rician-distributed magnitude multi-echo MRI signals into Gaussian-distributed signals with high accuracy and precision. We then perform an ILT on the transformed data to obtain an accurate T2 distribution. Additionally, we demonstrate, by simulations and experiments, that this approach corrects the aforementioned artifacts in magnitude multi-echo MR images over a large range of signal-to-noise ratios.

Multi-parametric MRI characterization of inflammation in murine skeletal muscle

Myopathies often display a common set of complex pathologies that include muscle weakness, inflammation, compromised membrane integrity, fat deposition, and fibrosis. Multi-parametric, quantitative, non-invasive imaging approaches may be able to resolve these individual pathological components. The goal of this study was to use multi-parametric MRI to investigate inflammation as an isolated pathological feature. Proton relaxation, diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT-MRI), and dynamic contrast enhanced (DCE-MRI) parameters were calculated from data acquired in a single imaging session conducted 6-8 hours following the injection of λ-carrageenan, a local inflammatory agent. T2 increased in the inflamed muscle and transitioned to bi-exponential behavior. In diffusion measurements, all three eigenvalues and the apparent diffusion coefficient increased, but λ3 had the largest relative change. Analysis of the qMT data revealed that the T1 of the free pool and the observed T1 both increased in the inflamed tissue, while the ratio of exchanging spins in the solid pool to those in the free water pool (the pool size ratio) significantly decreased. DCE-MRI data also supported observations of an increase in extracellular volume. These findings enriched the understanding of the relation between multiple quantitative MRI parameters and an isolated inflammatory pathology, and may potentially be employed for other single or complex myopathy models.

Effects of diffusion on high-resolution quantitative T2 MRI

Carr-Purcell-Meiboom-Gill-based sequences are often assumed to be insensitive to diffusion. However, imaging gradients always contribute some degree of diffusion weighting which increases with resolution. This may cause an apparent decrease in T2 when using a multi-echo sequence, such as quantitative T2 (qT2) at high resolution. This study investigated the impact of diffusion on the qT2 sequence. An equation was developed relating the diffusion factor associated with each echo (bqT2 ) to the underestimation of T2 , which was strongly dependent on both the actual T2 and the apparent diffusion coefficient of the tissue. The diffusion dependence of the measured T2 was demonstrated in rat spinal cord. The measured T2 was independent of the imaging plane in gray matter, where diffusion was isotropic, and orientation dependent in white matter, where diffusion was strongly anisotropic. The dependence of the measured T2 on the actual T2 value was also demonstrated in MnCl2 phantoms. The relationship between the resolution and underestimation of T2 was investigated both theoretically and experimentally for the original readout and a fully refocused readout. The fully refocused readout increased the resolution at which diffusion effects could be neglected whilst measuring T2 . To avoid the misidentification of cerebrospinal fluid when applying qT2 in the brain or spinal cord, a minimum in-plane voxel dimension of 0.2 mm was suggested for the standard qT2 sequence and 0.1 mm for the refocused readout. Simulations of myelin water fraction measurement indicated that signal-to-noise ratio requirements were increased in the presence of diffusion. Finally, the use of decreasing spoiler gradients to attenuate stimulated echoes should be avoided, as it was found to distort the T2 distribution when the slice thickness was less than 1 mm.

Effect of treatment with interferon beta-1a on changes in voxel-wise magnetization transfer ratio in normal appearing brain tissue and lesions of patients with relapsing-remitting multiple sclerosis: a 24-week, controlled pilot study

Background

This pilot study investigated changes in remyelinating and demyelinating activity in normal appearing brain tissue (NABT) and lesions, by using voxel-wise magnetization transfer ratio (VW-MTR), in patients with relapsing–remitting multiple sclerosis (RRMS) receiving interferon beta-1a 44 mcg subcutaneously (IFN β-1a SC) three times weekly versus healthy controls (HCs) (NCT01085318).

Methods

Increasing (suggestive of remyelination) and decreasing (suggestive of demyelination) VW-MTR changes in NABT and in T2, T1 and gadolinium (Gd)-enhancing lesion volume were measured over 24 weeks in 23 patients treated with IFN β-1a SC and in 15 HCs (where applicable). VW-MTR changes were tested using the Wilcoxon signed–rank or Wilcoxon rank–sum test.

Results

A trend for greater volume of NABT with increasing VW-MTR at 24 weeks was observed for patients versus HCs (median [range] 1206 [0–15278]; 342 [0–951] mm³; p = 0.061). NABT volume with increasing VW-MTR at 12 weeks was significantly greater in patients than in HCs (852 [6–11577]; 360 [0–1755] mm³; p = 0.028). Similar findings were detected for lesion volumes. Two patients with notably high numbers of Gd-enhancing lesions at baseline had a markedly greater volume of tissue with increasing VW-MTR compared with other patients. Volume of NABT tissue with decreasing VW-MTR was significantly greater in patients versus HCs at 24 weeks (942 [0–6141]; 297 [0–852] mm³; p<0.001).

Conclusions

The significant change in NABT volume with increasing VW-MTR at 12 weeks suggests that active remyelination in patients with RRMS may occur during treatment with IFN β-1a SC. Findings from two patients with the highest number of Gd-enhancing lesions at baseline suggest that extensive remyelination in NABT may occur in patients with high disease activity. Tissue volume with decreasing VW-MTR was greater in patients than in HCs, despite treatment, validating the sensitivity of this technique for detecting MS disease activity.

Trial Registration

ClinicalTrials.gov NCT01085318.

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.

Viable and fixed white matter: diffusion magnetic resonance comparisons and contrasts at physiological temperature

PURPOSE

Fixed samples have been used extensively in diffusion MRI (dMRI) studies. However, fixation causes significant structural changes in tissue. The purpose of this study was to evaluate fixed white matter as a surrogate for viable white matter during development and validation of dMRI methods.

METHODS

dMRI data was acquired from fixed and viable rat optic nerves maintained in identical conditions in a viable isolated tissue (VIT) chamber. The chamber preserves tissue integrity for 10 h at 37°C. Diffusion tensors (DT) and multi-compartment white matter signal models were fitted to the data.

RESULTS

When comparing VIT and fixed tissue, DT parameters demonstrated that fixation causes significant reductions in axial diffusivity and increases in radial diffusivity. However, both tissues exhibited similar responses to changes in diffusion times and gradient strengths. Multicompartment models demonstrated differences in parameter estimates (e.g., directional diffusivities) that were analogous to differences in DT parameters. Similarities in multi-compartment model rankings suggested that tissue water populations were broadly maintained postfixation.

CONCLUSIONS

The data demonstrate that fixed tissue, while maintaining the broad water environment of viable tissue, differs significantly in diffusion parameters. Results from dMRI experiments on fixed tissue may correlate with-but will not directly translate into-results from viable tissue.

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.

The radial diffusivity and magnetization transfer pool size ratio are sensitive markers for demyelination in a rat model of type III multiple sclerosis (MS) lesions

Determining biophysical sensitivity and specificity of quantitative magnetic resonance imaging is essential to develop effective imaging metrics of neurodegeneration. Among these metrics, apparent pool size ratio (PSR) from quantitative magnetization transfer (qMT) imaging and radial diffusivity (RD) from diffusion tensor imaging (DTI) are both known to relate to histological measure of myelin density and integrity. However their relative sensitivities towards quantitative myelin detection are unknown. In this study, we correlated high-resolution quantitative magnetic resonance imaging measures of subvoxel tissue structures with corresponding quantitative myelin histology in a lipopolysaccharide (LPS) mediated animal model of MS. Specifically, we acquired quantitative magnetization transfer (qMT) and diffusion tensor imaging (DTI) metrics (on the same tissue sample) in an animal model system of type III oligodendrogliopathy which lacked prominent lymphocytic infiltration, a system that had not been previously examined with quantitative MRI. We find that the qMT measured apparent pool size ratio (PSR) showed the strongest correlation with a histological measure of myelin content. DTI measured RD showed the next strongest correlation, and other DTI and relaxation parameters (such as the longitudinal relaxation rate (R1f) or fractional anisotropy (FA)) showed considerably weaker correlations with myelin content.

Multiexponential T(2) analyses in a murine model of hepatic fibrosis at 11.7 T MRI

This study evaluated the effects of hepatic fibrosis on the multiexponential T(2) (MET(2) ) relaxation of ex vivo murine liver specimens using an 11.7 T MRI. This animal study was approved by the Institutional Animal Care and Use Committee. Eighteen male C57BL/6 mice were divided into control (n = 3) and experimental (n = 15) groups; the latter group was fed a 3,5-dicarbethoxy-1,4-dihydrocollidine-supplemented diet to induce hepatic fibrosis. Ex vivo liver specimens were imaged using an 11.7 T MRI scanner. A multi-echo spin-echo sequence was utilized for subsequent MET(2) analysis. Degrees of fibrosis were determined by a pathologist, as well as by digital image analysis. Scatterplot graphs comparing various features of the MET(2) signal decay with the degrees of fibrosis were generated, and correlation coefficients were calculated. Two distinct peaks of the MET(2) signal decay were identified in all liver specimens: a short T(2) component with a geometric mean T(2) (GMT(2) ) approximating 30 ms; and a long T(2) component with GMT(2) approximating 400 ms. Strong correlation was found between the degree of hepatic fibrosis and the amplitude of the short T(2) component, with a higher degrees of fibrosis associated with a lower amplitude. Moderate correlation was also found between hepatic fibrosis and the GMT(2) values of the long T(2) component, with higher degrees of fibrosis associated with lower GMT(2) values. The study of hepatic microenvironments using MET(2) analyses offers potential utility in the ongoing development of the noninvasive assessment of hepatic fibrosis using MRI.