Characterization of the NMR behavior of white matter in bovine brain

Multicompartment imaging of the brain using a comprehensive MR imaging protocol

In this study, we describe a comprehensive 3D magnetic resonance imaging (MRI) protocol designed to assess major tissue and fluid components in the brain. The protocol comprises four different sequences: 1) magnetization transfer prepared Cones (MT-Cones) for two-pool MT modeling to quantify macromolecular content; 2) short-TR adiabatic inversion-recovery prepared Cones (STAIR-Cones) for myelin water imaging; 3) proton-density weighted Cones (PDw-Cones) for total water imaging; and 4) highly T2 weighted Cones (T2w-Cones) for free water imaging. By integrating these techniques, we successfully mapped key brain components-namely macromolecules, myelin water, intra/extracellular water, and free water-in ten healthy volunteers and five patients with multiple sclerosis (MS) using a 3T clinical scanner. Brain macromolecular proton fraction (MMPF), myelin water proton fraction (MWPF), intra/extracellular water proton fraction (IEWPF), and free water proton fraction (FWPF) values were generated in white matter (WM), grey matter (GM), and MS lesions. Excellent repeatability of the protocol was demonstrated with high intra-class correlation coefficient (ICC) values. In MS patients, the MMPF and MWPF values of the lesions and normal-appearing WM (NAWM) were significantly lower than those in normal WM (NWM) in healthy volunteers. Moreover, we observed significantly higher FWPF values in MS lesions compared to those in NWM and NAWM regions. This study demonstrates the capability of our technique to volumetrically map major brain components. The technique may have particular value in providing a comprehensive assessment of neuroinflammatory and neurodegenerative diseases of the brain.

Anisotropic longitudinal water proton relaxation in white matter investigated ex vivo in porcine spinal cord with sample rotation

A variation of the longitudinal relaxation time T 1 in brain regions that differ in their main fiber direction has been occasionally reported, however, with inconsistent results. Goal of the present study was to clarify such inconsistencies, and the origin of potential T 1 orientation dependence, by applying direct sample rotation and comparing the results from different approaches to measure T 1 . A section of fixed porcine spinal cord white matter was investigated at 3 T with variation of the fiber-to-field angle θ FB . The experiments included one-dimensional inversion-recovery, MP2RAGE, and variable flip-angle T 1 measurements at 22 °C and 36 °C as well as magnetization-transfer (MT) and diffusion-weighted acquisitions. Depending on the technique, different degrees of T 1 anisotropy (between 2 and 10%) were observed as well as different dependencies on θ FB (monotonic variation or T 1 maximum at 30-40°). More pronounced anisotropy was obtained with techniques that are more sensitive to MT effects. Furthermore, strong correlations of θ FB -dependent MT saturation and T 1 were found. A comprehensive analysis based on the binary spin-bath model for MT revealed an interplay of several orientation-dependent parameters, including the transverse relaxation times of the macromolecular and the water pool as well as the longitudinal relaxation time of the macromolecular pool.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics

Inversion pulses are commonly employed in MRI for T 1 $$ {T}_1 $$ -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 $$ {B}_0 $$ and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 $$ {T}_1 $$ relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous ¹ H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 $$ {T}_1 $$ relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

High-resolution magnetization-transfer imaging of post-mortem marmoset brain: Comparisons with relaxometry and histology

Cell membranes and macromolecules or paramagnetic compounds interact with water proton spins, which modulates magnetic resonance imaging (MRI) contrast providing information on tissue composition. For a further investigation, quantitative magnetization transfer (qMT) parameters (at 3T), including the ratio of the macromolecular and water proton pools, F, and the exchange-rate constant as well as the (observed) longitudinal and the effective transverse relaxation rates (at 3T and 7T), R₁^obs and R₂^*, respectively, were measured at high spatial resolution (200 µm) in a slice of fixed marmoset brain and compared to histology results obtained with Gallyas' myelin stain and Perls' iron stain. R₁^obs and R₂^* were linearly correlated with the iron content for the entire slice, whereas distinct differences were obtained between gray and white matter for correlations of relaxometry and qMT parameters with myelin content. The combined results suggest that the macromolecular pool interacting with water consists of myelin and (less efficient) non-myelin contributions. Despite strong correlation of F and R₁^obs, none of these parameters was uniquely specific to myelination. Due to additional sensitivity to iron stores, R₁^obs and R₂^* were more sensitive for depicting microstructural differences between cortical layers than F.

Synthetic MR: Physical principles, clinical implementation, and new developments

Current clinical MR imaging practices rely on the qualitative assessment of images for diagnosis and treatment planning. While contrast in MR images is dependent on the spin parameters of the imaged tissue, pixel values on MR images are relative and are not scaled to represent any tissue properties. Synthetic MR is a fully featured imaging workflow consisting of efficient multiparameter mapping acquisition, synthetic image generation, and volume quantitation of brain tissues. As the application becomes more widely available on multiple vendors and scanner platforms, it has also gained widespread adoption as clinicians begin to recognize the benefits of rapid quantitation. This review will provide details about the sequence with a focus on the physical principles behind its relaxometry mechanisms. It will present an overview of the products in their current form and some potential issues when implementing it in the clinic. It will conclude by highlighting some recent advances of the technique, including a 3D mapping method and its associated applications. This article is protected by copyright. All rights reserved.

Impact of the acquisition protocol on the sensitivity to demyelination and axonal loss of clinically feasible DWI techniques: a simulation study

OBJECTIVE

To evaluate: (a) the specific effect that the demyelination and axonal loss have on the DW signal, and (b) the impact of the sequence parameters on the sensitivity to damage of two clinically feasible DWI techniques, i.e. DKI and NODDI.

METHODS

We performed a Monte Carlo simulation of water diffusion inside a novel synthetic model of white matter in the presence of axonal loss and demyelination, with three compartments with permeable boundaries between them. We compared DKI and NODDI in their ability to detect and assess the damage, using several acquisition protocols. We used the F test statistic as an index of the sensitivity for each DWI parameter to axonal loss and demyelination, respectively.

RESULTS

DKI parameters significantly changed with increasing axonal loss, but, in most cases, not with demyelination; all the NODDI parameters showed sensitivity to both the damage processes (at p < 0.01). However, the acquisition protocol strongly affected the sensitivity to damage of both the DKI and NODDI parameters and, especially for NODDI, the parameter absolute values also.

DISCUSSION

This work is expected to impact future choices for investigating white matter microstructure in focusing on specific stages of the disease, and for selecting the appropriate experimental framework to obtain optimal data quality given the purpose of the experiment.

Deuterium double quantum-filtered NMR studies of peripheral and optic nerves

OBJECTIVE

Characterization of the nerve components by deuterium double quantum-filtered magnetization transfer (DQF-MT) NMR.

METHODS

Nerves were equilibrated in deuterated saline and ²H single-pulse and ²H DQF-MT NMR spectra were measured, enabling the separation of the different water compartments, according to their quadrupolar splittings.

RESULTS

Rat sciatic and brachial nerves and porcine optic nerve immersed in deuterated saline yielded ²H DQF spectra composed of three pairs of quadrupolar-split signals assigned to the water in the collagenous compartments and the myelin bilayer and one narrow signal assigned to the axonal water. Stretching of the nerves, application of osmotic stress and incubation in collagenase did not affect the quadrupolar splitting of the myelin water. The signals of myelin and axonal water were shown to decay during Wallerian degeneration and to rise during maturation. The chemical exchange between the myelin and the intra-axonal water was measured for optic nerve during maturation. The quadrupolar splitting of the signal of myelin water was not sensitive to its orientation relative to the magnetic field. This resembles liquid crystalline behavior, but leaves its mechanism open for interpretation.

CONCLUSIONS

²H DQF-MT NMR characterizes the different components of nerves, the water exchange between them and their changes during processes such as nerve maturation and Wallerian degeneration.

Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain

PURPOSE

To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types.

METHODS

A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles.

RESULTS

Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup.

CONCLUSIONS

Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.

Temperature dependence and histological correlation of inhomogeneous magnetization transfer and myelin water imaging in ex vivo brain

PURPOSE

The promise of inhomogeneous magnetization transfer (ihMT) as a new myelin imaging method was studied in ex vivo human brain tissue and in relation to myelin water fraction (MWF). The temperature dependence of both methods was characterized, as well as their correspondence with a histological measure of myelin content. Unfiltered and filtered ihMT protocols were studied by adjusting the saturation scheme to preserve or attenuate signal from tissue with short dipolar relaxation time T_1D.

METHODS

ihMT ratio (ihMTR) and MWF maps were acquired at 7 T from formalin-fixed human brain samples at 22.5 °C, 30 °C and 37 °C. The impact of temperature on unfiltered ihMTR, filtered ihMTR and MWF was investigated and compared to myelin basic protein staining.

RESULTS

Unfiltered ihMTR exhibited no temperature dependence, whereas filtered ihMTR increased with increasing temperature. MWF decreased at higher temperature, with an increasing prevalence of areas where the myelin water signal was unreliably determined, likely related to a reduction in T₂ peak separability at higher temperatures ex vivo. MWF and ihMTR showed similar per-sample correlation with myelin staining at room temperature. At 37 ^°C, filtered ihMTR was more strongly correlated with myelin staining and had increased dynamic range compared to unfiltered ihMTR.

CONCLUSIONS

Given the temperature dependence of filtered ihMT, increased dynamic range, and strong myelin specificity that persists at higher temperatures, we recommend carefully controlled temperatures close to 37 °C for filtered ihMT acquisitions. Unfiltered ihMT may also be useful, due to its independence from temperature, higher amplitude values, and sensitivity to short T_1D components. Ex vivo myelin water imaging should be performed at room temperature, to avoid fitting issues found at higher temperatures.

Identification of water compartments in spinal cords by 2 H double quantum filtered NMR

In ² H double quantum filtered (DQF) NMR, the various water compartments are characterized by their different residual quadrupolar interactions. The spectral separation between the different signals enables the measurement of the relaxation of each compartment and the magnetization transfer (MT) between them. In the current study, five water compartments were identified in the ² H DQF spectra of porcine spinal cord. The most prominent signal was the pair of satellites with a quadrupolar splitting of about 550 Hz. ² H DQF MRI optimized for the 550 Hz quadrupolar splitting indicated that this signal originated mainly from the white matter and it was assigned to the myelin water. This splitting does not change upon changing the orientation of the spinal cord relative to the magnetic field, indicating a liquid crystalline nature. Another site exhibiting splitting of about 1500 Hz was assigned to collagenous connective tissue. The narrow central peak was assigned to a combination of intra- and inter-axonal water. The assignment of the other two sites is not certain and requires further study. The rates of MT between the various sites were recorded.

Understanding aqueous and non-aqueous proton T1 relaxation in brain

A full picture of longitudinal relaxation in complex heterogeneous environments like white matter brain tissue remains elusive. In tissue, successive approximations, from the solvation layer model to the two pool model, have highlighted how longitudinal magnetization evolution depends on both inter-compartmental exchange and spin-lattice relaxation. In white matter, however, these models fail to capture the behaviour of the two distinct aqueous pools, myelin water and intra/extra-cellular water. A challenge with testing more comprehensive multi-pool models lies in directly observing all pools, both aqueous and non-aqueous. In this work, we advance these efforts by integrating three main experimental and analytical elements: direct observation of the longitudinal relaxation of both the aqueous and the non-aqueous protons in white matter, a wide range of different initial conditions, and application of an analysis pipeline which includes lineshape, CPMG, and fitting of a four pool model. An eigenvector interpretation of the four pool model highlights how longitudinal relaxation in white matter depends on initial conditions. We find that a single set of model parameters is able to describe the entire range of relaxation behaviour observed in all the separable aqueous and non-aqueous pools in experiments involving six different initial conditions. Understanding of the nature and connectedness of the tissue components is crucial in the design and interpretation of many MRI measurements, especially those based on magnetization transfer and longitudinal relaxation. In particular, the dependency of relaxation behaviour on initial conditions is likely the basis for understanding method-dependent discrepancies in in vivo T₁.

Comparison of multi echo T2 relaxation and steady state approaches for myelin imaging in the central nervous system

The traditional approach for measuring myelin-associated water with quantitative magnetic resonance imaging (MRI) uses multi-echo T₂ relaxation data to calculate the myelin water fraction (MWF). A fundamentally different approach, abbreviated “mcDESPOT”, uses a more efficient steady-state acquisition to generate an equivalent metric (f_M). Although previous studies have demonstrated inherent instability and bias in the complex mcDESPOT analysis procedure, f_M has often been used as a surrogate for MWF. We produced and compared multivariate atlases of MWF and f_M in healthy human brain and cervical spinal cord (available online) and compared their ability to detect multiple sclerosis pathology. A significant bias was found in all regions (p < 10^–5), albeit reversed for spinal cord (f_M-MWF =  − 3.4%) compared to brain (+ 6.2%). MWF and f_M followed an approximately linear relationship for regions with MWF <  ~ 10%. For MWF >  ~ 10%, the relationship broke down and f_M no longer increased in tandem with MWF. For multiple sclerosis patients, MWF and f_M Z score maps showed overlapping areas of low Z score and similar trends between patients and brain regions, although those of f_M generally had greater spatial extent and magnitude of severity. These results will guide future choice of myelin-sensitive quantitative MRI and improve interpretation of studies using either myelin imaging approach.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Stability of the Meat Protein Type I Collagen: Influence of pH, Ionic Strength, and Phenolic Antioxidant

The water-holding capacity (WHC) is among the key factors in determining the quality of meat and its value, which is strongly influenced by the content and quality of the connective tissue proteins like collagen. Therefore, the factors that influence the proteins’ stability, e.g., pH, ionic strength, and the antioxidants which are used to increase the meat shelf-life, also affect the WHC. The interaction of collagen, whose structure is strongly influenced by the interaction with water molecules, can be studied following the behavior of water diffusion by low-resolution ¹H NMR experiments. The present study is addressed to study the collagen stability as a function of pH, ionic strength, and phenolic antioxidants like catechin. The experimental study demonstrated how the ¹H NMR time domain (TD) experiments are able to evaluate the hydration properties of collagen, not only as a function of ionic strength and pH, but also in determining the ability of catechin to interact both on the surface of the collagen fibrils and inside the fibrillar domain.

Origin of orientation-dependent R1 (=1/T1 ) relaxation in white matter

PURPOSE

In a recent MRI study, it was shown that the longitudinal relaxation rate, R₁ , in white matter (WM) is influenced by the relative orientation of nerve fibers with respect to the main magnetic field (B₀ ). Even though the exact nature of this R₁ orientation dependency is still unclear, it can be assumed that the origin of the phenomenon can be attributed to the anisotropic and unique molecular environment within the myelin sheath surrounding the axons. The current work investigates the contribution of dipolar induced R₁ relaxation of the myelin associated hydrogen nuclei theoretically and compares the results with the experimentally observed R₁ orientation dependency.

METHODS

Atomistic molecular dynamics simulations were employed and the R₁ relaxation rate of hydrogen nuclei of a myelin-alike molecular environment was calculated for various orientations of the trajectory sets relative to the B₀ -field. Based on the calculated relaxation rates, the observable R₁ relaxation was simulated for various fiber orientations and fitted to the experimental data using a suitable signal weighting-scheme.

RESULTS

The results obtained show that the R₁ relaxation rate of both solid myelin (SM) and myelin water (MW) depends on the fiber orientation relative to the main B₀ -field. Moreover, employing a realistic signal weighing scheme and tissue characteristics, the theoretically investigated R₁ orientation dependency matches the experimental data well.

CONCLUSION

The good agreement between theoretical and experimental findings indicates that the R₁ orientation dependency in WM mainly originates from anisotropic dipole-dipole interactions between hydrogen nuclei located within the myelin sheath.

Dipolar induced spin-lattice relaxation in the myelin sheath: A molecular dynamics study

Interactions between hydrogen protons of water molecules and macromolecules within the myelin sheath surrounding the axons are a major factor influencing the magnetic resonance (MR) contrast in white matter (WM) regions. In past decades, several studies have investigated the underlying effects and reported a wide range of R1 rates for the myelin associated compartments at different field strengths. However, it was also shown that the experimental quantification of the compartment-specific R1 rates is associated with large uncertainties. The current study therefore investigates the longitudinal relaxation rates within the myelin sheath using a molecular dynamic (MD) simulation. For this purpose, a realistic molecular model of the myelin sheath was employed to determine the dipole-dipole induced R1 relaxation rate of the hydrogen protons at clinically relevant field strengths. The results obtained clearly reflect the spatial heterogeneity of R1 with a increased relaxivity of myelin water due to a reduced molecular mobility near the membrane surface. Moreover, the calculated R1 rates for both myelin water and macromolecules are in excellent agreement with experimental findings from the literature at different field strengths.

Iron, Myelin, and the Brain: Neuroimaging Meets Neurobiology

Although iron is crucial for neuronal functioning, many aspects of cerebral iron biology await clarification. The ability to quantify specific iron forms in the living brain would open new avenues for diagnosis, therapeutic monitoring, and understanding pathogenesis of diseases. A modality that allows assessment of brain tissue composition in vivo, in particular of iron deposits or myelin content on a submillimeter spatial scale, is magnetic resonance imaging (MRI). Multimodal strategies combining MRI with complementary analytical techniques ex vivo have emerged, which may lead to improved specificity. Interdisciplinary collaborations will be key to advance beyond simple correlative analyses in the biological interpretation of MRI data and to gain deeper insights into key factors leading to iron accumulation and/or redistribution associated with neurodegeneration.

White matter intercompartmental water exchange rates determined from detailed modeling of the myelin sheath

PURPOSE

Magnetization exchange (ME) between hydrogen protons of water and large molecules (semisolids [SS]) in lipid bilayers is an important factor in MRI signal generation and can be exploited to study white matter pathology. Current models used to quantify ME in white matter generally consider water to reside in 1 or 2 distinct compartments, ignoring the complexities of the myelin sheath's multicompartment structure of alternating myelin SS and myelin water (MW) layers. Here, we investigated the effect of this by fitting ME data obtained from human brain at 7 T with a multilayer model of myelin.

METHODS

A multi-echo acquisition for a T₂ ^* -based separation of MW from other water signals was combined with various preparation pulses to change the (relative) state of the SS and water pools and analyzed by fitting with a multilayer exchange model.

RESULTS

The estimated lifetime within a single MW layer was 260 µs, corresponding to a lipid bilayer permeability of 6.7 µm/s. The magnetization lifetime of the aggregate of all MW was estimated at 13 ms, shorter than previously reported values in the range of 40 to 140 ms.

CONCLUSION

Contrary to expectations and previous reports, ME between protons in myelin SS and water is not limited by the myelin sheath but rather by the exchange between SS and water protons. The analysis of ME contrast should account for the relatively short MW lifetime and affects the interpretation of tissue compartmentalization from MRI contrasts such as T₁ - and diffusion-weighting.

SyMRI of the Brain: Rapid Quantification of Relaxation Rates and Proton Density, With Synthetic MRI, Automatic Brain Segmentation, and Myelin Measurement

Conventional magnetic resonance images are usually evaluated using the image signal contrast between tissues and not based on their absolute signal intensities. Quantification of tissue parameters, such as relaxation rates and proton density, would provide an absolute scale; however, these methods have mainly been performed in a research setting. The development of rapid quantification, with scan times in the order of 6 minutes for full head coverage, has provided the prerequisites for clinical use. The aim of this review article was to introduce a specific quantification method and synthesis of contrast-weighted images based on the acquired absolute values, and to present automatic segmentation of brain tissues and measurement of myelin based on the quantitative values, along with application of these techniques to various brain diseases. The entire technique is referred to as “SyMRI” in this review. SyMRI has shown promising results in previous studies when used for multiple sclerosis, brain metastases, Sturge-Weber syndrome, idiopathic normal pressure hydrocephalus, meningitis, and postmortem imaging.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

The role of magnetization transfer in the observed contrast in T1 weighted imaging under clinical setups

In T₁ weighted magnetic resonance imaging of brain and spinal cord in the clinical setting, the white matter (WM) appears with greater intensity than the gray matter (GM). This contrast has been assigned to differences in T₁ values. In these experiments the RF pulses are too long to excite both the water and the species with restricted motion of the protons (SRMP). In in vitro studies using short RF pulses, the contrast is reversed, with greater intensity for the GM. These results raise the question of whether magnetization transfer (MT) plays a role in the contrast observed in the T₁ weighting experiments. In the present work we implemented selective saturation recovery alone and together with the conventional magnetization transfer contrast (MTC) method. The results confirm that a major factor that determines the characteristic WM/GM averaged intensity ratio observed in T₁ weighted imaging under clinical conditions is MT between the SRMP and water. When selective saturation recovery is combined with MTC, the SRMP yields spectral widths ranging from a few to tens of kilohertz, indicating that more than one type of SRMP is involved in the MT. The z-spectrum obtained with this combination is free of the effect of direct saturation of the water peak. Selective saturation recovery enables an independent measurement of the exchange time and T₁ , while the combination with MTC with complete saturation of the SRMP enables measurement of T₁ without the effect of MT. The latter measurement can be carried out on a timescale much shorter than T_1.

Brain microstructure by multi-modal MRI: Is the whole greater than the sum of its parts?

The MRI signal is dependent upon a number of sub-voxel properties of tissue, which makes it potentially able to detect changes occurring at a scale much smaller than the image resolution. This "microstructural imaging" has become one of the main branches of quantitative MRI. Despite the exciting promise of unique insight beyond the resolution of the acquired images, its widespread application is limited by the relatively modest ability of each microstructural imaging technique to distinguish between differing microscopic substrates. This is mainly due to the fact that MRI provides a very indirect measure of the tissue properties in which we are interested. A strategy to overcome this limitation lies in the combination of more than one technique, to exploit the relative contributions of differing physiological and pathological substrates to selected MRI contrasts. This forms the basis of multi-modal MRI, a broad concept that refers to many different ways of effectively combining information from more than one MRI contrast. This paper will review a range of methods that have been proposed to maximise the output of this combination, primarily falling into one of two approaches. The first one relies on data-driven methods, exploiting multivariate analysis tools able to capture overlapping and complementary information. The second approach, which we call "model-driven", aims at combining parameters extracted by existing biophysical or signal models to obtain new parameters, which are believed to be more accurate or more specific than the original ones. This paper will attempt to provide an overview of the advantages and limitations of these two philosophies.

Inversion recovery ultrashort echo time imaging of ultrashort T2 tissue components in ovine brain at 3 T: a sequential D2 O exchange study

Inversion recovery ultrashort echo time (IR-UTE) imaging holds the potential to directly characterize MR signals from ultrashort T2 tissue components (STCs), such as collagen in cartilage and myelin in brain. The application of IR-UTE for myelin imaging has been challenging because of the high water content in brain and the possibility that the ultrashort T2 * signals are contaminated by water protons, including those associated with myelin sheaths. This study investigated such a possibility in an ovine brain D2 O exchange model and explored the potential of IR-UTE imaging for the quantification of ultrashort T2 * signals in both white and gray matter at 3 T. Six specimens were examined before and after sequential immersion in 99.9% D2 O. Long T2 MR signals were measured using a clinical proton density-weighted fast spin echo (PD-FSE) sequence. IR-UTE images were first acquired with different inversion times to determine the optimal inversion time to null the long T2 signals (TInull ). Then, at this TInull , images with echo times (TEs) of 0.01-4 ms were acquired to measure the T2 * values of STCs. The PD-FSE signal dropped to near zero after 24 h of immersion in D2 O. A wide range of TInull values were used at different time points (240-330 ms for white matter and 320-350 ms for gray matter at TR = 1000 ms) because the T1 values of the long T2 tissue components changed significantly. The T2 * values of STCs were 200-300 μs in both white and gray matter (comparable with the values obtained from myelin powder and its mixture with D2 O or H2 O), and showed minimal changes after sequential immersion. The ultrashort T2 * signals seen on IR-UTE images are unlikely to be from water protons as they are exchangeable with deuterons in D2 O. The source is more likely to be myelin itself in white matter, and might also be associated with other membranous structures in gray matter.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.

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.

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.

Utility of a Multiparametric Quantitative MRI Model That Assesses Myelin and Edema for Evaluating Plaques, Periplaque White Matter, and Normal-Appearing White Matter in Patients with Multiple Sclerosis: A Feasibility Study

BACKGROUND AND PURPOSE

T1 and T2 values and proton density can now be quantified on the basis of a single MR acquisition. The myelin and edema in a voxel can also be estimated from these values. The purpose of this study was to evaluate a multiparametric quantitative MR imaging model that assesses myelin and edema for characterizing plaques, periplaque white matter, and normal-appearing white matter in patients with MS.

MATERIALS AND METHODS

We examined 3T quantitative MR imaging data from 21 patients with MS. The myelin partial volume, excess parenchymal water partial volume, the inverse of T1 and transverse T2 relaxation times (R1, R2), and proton density were compared among plaques, periplaque white matter, and normal-appearing white matter.

RESULTS

All metrics differed significantly across the 3 groups (P < .001). Those in plaques differed most from those in normal-appearing white matter. The percentage changes of the metrics in plaques and periplaque white matter relative to normal-appearing white matter were significantly more different from zero for myelin partial volume (mean, -61.59 ± 20.28% [plaque relative to normal-appearing white matter], and mean, -10.51 ± 11.41% [periplaque white matter relative to normal-appearing white matter]), and excess parenchymal water partial volume (13.82 × 10³ ± 49.47 × 10³% and 51.33 × 10² ± 155.31 × 10²%) than for R1 (-35.23 ± 13.93% and -6.08 ± 8.66%), R2 (-21.06 ± 11.39% and -4.79 ± 6.79%), and proton density (23.37 ± 10.30% and 3.37 ± 4.24%).

CONCLUSIONS

Multiparametric quantitative MR imaging captures white matter damage in MS. Myelin partial volume and excess parenchymal water partial volume are more sensitive to the MS disease process than R1, R2, and proton density.

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.

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.

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.

Modeling the Presence of Myelin and Edema in the Brain Based on Multi-Parametric Quantitative MRI

The aim of this study was to present a model that uses multi-parametric quantitative MRI to estimate the presence of myelin and edema in the brain. The model relates simultaneous measurement of R₁ and R₂ relaxation rates and proton density to four partial volume compartments, consisting of myelin partial volume, cellular partial volume, free water partial volume, and excess parenchymal water partial volume. The model parameters were obtained using spatially normalized brain images of a group of 20 healthy controls. The pathological brain was modeled in terms of the reduction of myelin content and presence of excess parenchymal water, which indicates the degree of edema. The method was tested on spatially normalized brain images of a group of 20 age-matched multiple sclerosis (MS) patients. Clear differences were observed with respect to the healthy controls: the MS group had a 79 mL smaller brain volume (1069 vs. 1148 mL), a 38 mL smaller myelin volume (119 vs. 157 mL), and a 21 mL larger excess parenchymal water volume (78 vs. 57 mL). Template regions of interest of various brain structures indicated that the myelin partial volume in the MS group was 1.6 ± 1.5% lower for gray matter (GM) structures and 2.8 ± 1.0% lower for white matter (WM) structures. The excess parenchymal water partial volume was 9 ± 10% larger for GM and 5 ± 2% larger for WM. Manually placed ROIs indicated that the results using the template ROIs may have suffered from loss of anatomical detail due to the spatial normalization process. Examples of the application of the method on high-resolution images are provided for three individual subjects: a 45-year-old healthy subject, a 72-year-old healthy subject, and a 45-year-old MS patient. The observed results agreed with the expected behavior considering both age and disease. In conclusion, the proposed model may provide clinically important parameters, such as the total brain volume, degree of myelination, and degree of edema, based on a single qMRI acquisition with a clinically acceptable scan time.

In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter

The presence of micrometer-level restrictions leads to a decrease of diffusion coefficient with diffusion time. Here we investigate this effect in human white matter in vivo. We focus on a broad range of diffusion times, up to 600 ms, covering diffusion length scales up to about 30 μm. We perform stimulated echo diffusion tensor imaging on 5 healthy volunteers and observe a relatively weak time-dependence in diffusion transverse to major fiber tracts. Remarkably, we also find notable time-dependence in the longitudinal direction. Comparing models of diffusion in ordered, confined and disordered media, we argue that the time-dependence in both directions can arise due to structural disorder, such as axonal beads in the longitudinal direction, and the random packing geometry of fibers within a bundle in the transverse direction. These time-dependent effects extend beyond a simple picture of Gaussian compartments, and may lead to novel markers that are specific to neuronal fiber geometry at the micrometer scale.

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.

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.

Influence of water compartmentation and heterogeneous relaxation on quantitative magnetization transfer imaging in rodent brain tumors

PURPOSE

The goal of this study was to investigate the influence of water compartmentation and heterogeneous relaxation properties on quantitative magnetization transfer (qMT) imaging in tissues, and in particular whether a two-pool model is sufficient to describe qMT data in brain tumors.

METHODS

Computer simulations and in vivo experiments with a series of qMT measurements before and after injection of Gd-DTPA were performed. Both off-resonance pulsed saturation (pulsed) and on-resonance selective inversion recovery (SIR) qMT methods were used, and all data were fit with a two-pool model only.

RESULTS

Simulations indicated that a two-pool fitting of four-pool data yielded accurate measures of pool size ratio (PSR) of macromolecular versus free water protons when there were fast transcytolemmal exchange and slow R1 recovery. The fitted in vivo PSR of both pulsed and SIR qMT methods showed no dependence on R1 variations caused by different concentrations of Gd-DTPA during wash-out, whereas the fitted kex (magnetization transfer exchange rate) changed significantly with R1 .

CONCLUSION

A two-pool model provides reproducible estimates of PSR in brain tumors independent of relaxation properties in the presence of relatively fast transcytolemmal exchange, whereas estimates of kex are biased by relaxation variations. In addition, estimates of PSR in brain tumors using the pulsed and SIR qMT methods agree well with one another. Magn Reson Med 76:635-644, 2016. © 2015 Wiley Periodicals, Inc.

Structure predicts function: combining non-invasive electrophysiology with in-vivo histology

We present an approach for combining high resolution MRI-based myelin mapping with functional information from electroencephalography (EEG) or magnetoencephalography (MEG). The main contribution to the primary currents detectable with EEG and MEG comes from ionic currents in the apical dendrites of cortical pyramidal cells, aligned perpendicularly to the local cortical surface. We provide evidence from an in-vivo experiment that the variation in MRI-based myeloarchitecture measures across the cortex predicts the variation of the current density over individuals and thus is of functional relevance. Equivalent current dipole locations and moments due to pitch onset evoked response fields (ERFs) were estimated by means of a variational Bayesian algorithm. The myeloarchitecture was estimated indirectly from individual high resolution quantitative multi-parameter maps (MPMs) acquired at 800μm isotropic resolution. Myelin estimates across cortical areas correlated positively with dipole magnitude. This correlation was spatially specific: regions of interest in the auditory cortex provided significantly better models than those covering whole hemispheres. Based on the MPM data we identified the auditory cortical area TE1.2 as the most likely origin of the pitch ERFs measured by MEG. We can now proceed to exploit the higher spatial resolution of quantitative MPMs to identify the cortical origin of M/EEG signals, inform M/EEG source reconstruction and explore structure-function relationships at a fine structural level in the living human brain.

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.

Hypomyelinating leukodystrophies: translational research progress and prospects

Hypomyelinating leukodystrophies represent a genetically heterogeneous but clinically overlapping group of heritable disorders. Current management approaches in the care of the patient with a hypomyelinating leukodystrophy include use of serial magnetic resonance imaging (MRI) to establish and monitor hypomyelination, molecular diagnostics to determine a specific etiology, and equally importantly, careful attention to neurologic complications over time. Emerging research in oligodendrocyte biology and neuroradiology with bedside applications may result in the possibility of clinical trials in the near term, yet there are significant gaps in knowledge in disease classification, characterization, and outcome measures in this group of disorders. Here we review the biological background of myelination, the clinical and genetic variability in hypomyelinating leukodystrophies, and the insights that can be obtained from current MRI techniques. In addition, we discuss ongoing research approaches to define potential outcome markers for future clinical trials.

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.

Rapid multi-echo measurement of brain metabolite T₂ values at 7 T using a single-shot spectroscopic Carr-Purcell-Meiboom-Gill sequence and prior information

We present a method for the robust and accurate estimation of brain metabolite transverse relaxation times (T2 ) from multiple spin-echo data acquired with a single-shot Carr-Purcell-Meiboom-Gill (CPMG) spectroscopic sequence. Each acquired echo consists of a small number of complex time-domain data points. The amplitudes of the spectral components in each echo are calculated by solving a set of linear equations in which previously estimated frequencies and linewidths serve as prior information. These priors are obtained from a short MRS experiment in which a large number of time-domain data points are acquired, and are subsequently estimated using linear prediction with singular value decomposition (LPSVD) processing. We show that this process can be used to accurately and rapidly measure the T2 values for the main singlet resonances in single-volume MRS measurements in the brain. The proposed method can be generalized to any set of MRS experiments comprising repeated measurements of amplitude changes, e.g. as a function of an experimental parameter, such as TE, inversion time or diffusion weighting.

Temporal phase correction of multiple echo T2 magnetic resonance images

Typically, magnetic resonance imaging (MRI) analysis is performed on magnitude data, and multiple echo T2 data consist of numerous images of the same slice taken with different echo spacing, giving voxel-wise temporal sampling of the noise as the signals decay according to T2 relaxation. Magnitude T2 decay data has Rician distributed noise which is characterized by a change in the noise distribution from Gaussian, through a transitional region, to Rayleigh as the signal to noise ratio decreases with increasing echo time. Non-Gaussian noise distributions may produce errors in the commonly applied non-negative least squares (NNLS) algorithm that is used to assess multiple echo decays for compartmentalized water environments through the creation of T2 distributions. Typically, Gaussian noise is sought by performing spatial-based phase correction on the MRI data however, these methods cannot capitalize on the temporal information available from multiple echo T2 acquisitions. Here we describe a temporal phase correction (TPC) algorithm that utilizes the temporal noise information available in multiple echo T2 acquisitions to put the relevant decay information in the Real portion of the decay data and leave only noise in the Imaginary portion. We apply this TPC algorithm to create real-valued multiple echo T2 data from human subjects measured at 1.5 T. We show that applying TPC causes changes in the T2 distribution estimates; notably the possible resolution of separate extracellular and intracellular water environments, and the disappearance of the commonly labeled cerebrospinal fluid peak, which might be an artefact observed in many previously published multiple echo T2 analyses.

Matrix-algebra-based calculations of the time evolution of the binary spin-bath model for magnetization transfer

Quantification of magnetization-transfer (MT) experiments are typically based on the assumption of the binary spin-bath model. This model allows for the extraction of up to six parameters (relative pool sizes, relaxation times, and exchange rate constants) for the characterization of macromolecules, which are coupled via exchange processes to the water in tissues. Here, an approach is presented for estimating MT parameters acquired with arbitrary saturation schemes and imaging pulse sequences. It uses matrix algebra to solve the Bloch-McConnell equations without unwarranted simplifications, such as assuming steady-state conditions for pulsed saturation schemes or neglecting imaging pulses. The algorithm achieves sufficient efficiency for voxel-by-voxel MT parameter estimations by using a polynomial interpolation technique. Simulations, as well as experiments in agar gels with continuous-wave and pulsed MT preparation, were performed for validation and for assessing approximations in previous modeling approaches. In vivo experiments in the normal human brain yielded results that were consistent with published data.