Micro-compartment specific T2* relaxation in the brain

Generalized Lorentzian Tensor Approach (GLTA) as a biophysical background for quantitative susceptibility mapping

PURPOSE

Quantitative susceptibility mapping (QSM) is a potentially powerful technique for mapping tissue magnetic susceptibility from gradient recalled echo (GRE) MRI. Herein we aim to derive the relationships between GRE signal phase and the underlying tissue microstructure and magnetic susceptibility at the cellular level.

METHODS

We use Maxwell's equations and a statistical approach to derive the expression for the magnetic-susceptibility-induced MR signal frequency shift of the GRE signal in single- and multicompartment systems, in which inhomogeneous magnetic field is induced by the cellular constituents (proteins, lipids, iron, etc.) distributed in intra- and extracellular spaces.

RESULTS

We introduce the Generalized Lorentzian Tensor Approach (GLTA) that accounts for both types of anisotropy: the anisotropy of magnetic susceptibility and the structural tissue anisotropy. In the GLTA the frequency shift due to the local environment is characterized by the Lorentzian tensor L⁁ which has a substantially different structure than the susceptibility tensor χ⁁. While components of χ⁁ are simply compartmental susceptibilities "weighted" by their relative volumes, the components of L⁁ are weighted by specific numerical factors depending on tissue micro-symmetry and parameters related to the MR pulse sequence. We also provide equations bridging phenomenological and microscopic considerations.

CONCLUSION

The GLTA provides a consistent background for deciphering phase data.

The role of iron in brain ageing and neurodegenerative disorders

SUMMARY

In the CNS, iron in several proteins is involved in many important processes such as oxygen transportation, oxidative phosphorylation, myelin production, and the synthesis and metabolism of neurotransmitters. Abnormal iron homoeostasis can induce cellular damage through hydroxyl radical production, which can cause the oxidation and modification of lipids, proteins, carbohydrates, and DNA. During ageing, different iron complexes accumulate in brain regions associated with motor and cognitive impairment. In various neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, changes in iron homoeostasis result in altered cellular iron distribution and accumulation. MRI can often identify these changes, thus providing a potential diagnostic biomarker of neurodegenerative diseases. An important avenue to reduce iron accumulation is the use of iron chelators that are able to cross the blood-brain barrier, penetrate cells, and reduce excessive iron accumulation, thereby affording neuroprotection.

Reproducibility of T2 * mapping in the human cerebral cortex in vivo at 7 tesla MRI

BACKGROUND

To assess the test-retest reproducibility of cortical mapping of T2 * relaxation rates at 7 Tesla (T) MRI. T2 * maps have been used for studying cortical myelo-architecture patterns in vivo and for characterizing conditions associated with changes in iron and/or myelin concentration.

METHODS

T2 * maps were calculated from 7T multi-echo T2 *-weighted images acquired during separate scanning sessions on 8 healthy subjects. The reproducibility of surface-based cortical T2 * mapping was assessed at different depths of the cortex; from pial surface (0% depth) towards gray/white matter boundary (100% depth), across cortical regions and hemispheres, using coefficients of variation (COVs = SD/mean) between each couple (scan-rescan) of average T2 * measurements.

RESULTS

Average cortical T2 * was significantly different among 25%, 50%, and 75% depths (analysis of variance, P < 0.001). Coefficient of variations were very low within cortical regions, and whole cortex (average COV = 0.83-1.79%), indicating a high degree of reproducibility in T2 * measures.

CONCLUSION

Surface-based mapping of T2 * relaxation rates as a function of cortical depth is reproducible and could prove useful for studying the laminar architecture of the cerebral cortex in vivo, and for investigating physiological and pathological states associated with changes in iron and/or myelin concentration.

A torque balance measurement of anisotropy of the magnetic susceptibility in white matter

PURPOSE

Recent MRI studies have suggested that the magnetic susceptibility of white matter (WM) in the human brain is anisotropic, providing a new contrast mechanism for the visualization of fiber bundles and allowing the extraction of cellular compartment-specific information. This study provides an independent confirmation and quantification of this anisotropy.

METHODS

Anisotropic magnetic susceptibility results in a torque exerted on WM when placed in a uniform magnetic field, tending to align the WM fibers with the field. To quantify the effect, excised spinal cord samples were placed in a torque balance inside the magnet of a 7 T MRI system and the magnetic torque was measured as function of orientation.

RESULTS

All tissue samples (n = 5) showed orienting effects, confirming the presence of anisotropic susceptibility. Analysis of the magnetic torque resulted in reproducible values for the WM volume anisotropy that ranged from 13.6 to 19.2 ppb.

CONCLUSION

The independently determined anisotropy values confirm estimates inferred from MRI experiments and validate the use of anisotropy to extract novel information about brain fiber structure and myelination.

Quantitative susceptibility mapping: current status and future directions

Quantitative susceptibility mapping (QSM) is a new technique for quantifying magnetic susceptibility. It has already found various applications in quantifying in vivo iron content, calcifications and changes in venous oxygen saturation. The accuracy of susceptibility mapping is dependent on several factors. In this review, we evaluate the entire process of QSM from data acquisition to individual data processing steps. We also show preliminary results of several new concepts introduced in this review in an attempt to improve the quality and accuracy for certain steps. The uncertainties in estimating susceptibility differences using susceptibility maps, phase images, and T2* maps are analyzed and compared. Finally, example clinical applications are presented. We conclude that QSM holds great promise in quantifying iron and becoming a standard clinical tool.

The role of iron in brain ageing and neurodegenerative disorders

In the CNS, iron in several proteins is involved in many important processes such as oxygen transportation, oxidative phosphorylation, myelin production, and the synthesis and metabolism of neurotransmitters. Abnormal iron homoeostasis can induce cellular damage through hydroxyl radical production, which can cause the oxidation and modification of lipids, proteins, carbohydrates, and DNA. During ageing, different iron complexes accumulate in brain regions associated with motor and cognitive impairment. In various neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, changes in iron homoeostasis result in altered cellular iron distribution and accumulation. MRI can often identify these changes, thus providing a potential diagnostic biomarker of neurodegenerative diseases. An important avenue to reduce iron accumulation is the use of iron chelators that are able to cross the blood-brain barrier, penetrate cells, and reduce excessive iron accumulation, thereby affording neuroprotection.

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.

Mean magnetic susceptibility regularized susceptibility tensor imaging (MMSR-STI) for estimating orientations of white matter fibers in human brain

PURPOSE

An increasing number of studies show that magnetic susceptibility in white matter fibers is anisotropic and may be described by a tensor. However, the limited head rotation possible for in vivo human studies leads to an ill-conditioned inverse problem in susceptibility tensor imaging (STI). Here we suggest the combined use of limiting the susceptibility anisotropy to white matter and imposing morphology constraints on the mean magnetic susceptibility (MMS) for regularizing the STI inverse problem.

METHODS

The proposed MMS regularized STI (MMSR-STI) method was tested using computer simulations and in vivo human data collected at 3T. The fiber orientation estimated from both the STI and MMSR-STI methods was compared to that from diffusion tensor imaging (DTI).

RESULTS

Computer simulations show that the MMSR-STI method provides a more accurate estimation of the susceptibility tensor than the conventional STI approach. Similarly, in vivo data show that use of the MMSR-STI method leads to a smaller difference between the fiber orientation estimated from STI and DTI for most selected white matter fibers.

CONCLUSION

The proposed regularization strategy for STI can improve estimation of the susceptibility tensor in white matter.

Animal models and high field imaging and spectroscopy

A plethora of magnetic resonance (MR) techniques developed in the last two decades provide unique and noninvasive measurement capabilities for studies of basic brain function and brain diseases in humans. Animal model experiments have been an indispensible part of this development. MR imaging and spectroscopy measurements have been employed in animal models, either by themselves or in combination with complementary and often invasive techniques, to enlighten us about the information content of such MR methods and/or verify observations made in the human brain. They have also been employed, with or independently of human efforts, to examine mechanisms underlying pathological developments in the brain, exploiting the wealth of animal models available for such studies. In this endeavor, the desire to push for ever-higher spatial and/or spectral resolution, better signal-to-noise ratio, and unique image contrast has inevitably led to the introduction of increasingly higher magnetic fields. As a result, today, animal model studies are starting to be conducted at magnetic fields ranging from ~ 11 to 17 Tesla, significantly enhancing the armamentarium of tools available for the probing brain function and brain pathologies.

Comparison of R2' measurement methods in the normal brain at 3 Tesla

PURPOSE

R2', the reversible component of transverse relaxation, is an important susceptibility measurement for studies of brain physiology and pathologies. In existing literature, different R2' measurement methods are used with assumption of equivalency. This study explores the choice of measurement method in healthy, young subjects at 3T.

METHODS

In this study, a modified gradient-echo sampling of free induction decay and echo (GESFIDE) sequence was used to compare four standard R2' measurement methods: asymmetric spin echo (ASE), standard GESFIDE, gradient echo sampling of the spin echo (GESSE), and separate R2 and R2* mapping.

RESULTS

GESSE returned lower R2' measurements than other methods (P < 0.05). Intersubject mean R2' in gray matter was found to be 2.7 s(-1) using standard GESFIDE and GESSE, versus 3.4-3.8 s(-1) using other methods. In white matter, mean R2' from GESSE was 2.3 s(-1) while other methods produced 3.7-4.3 s(-1) . R2 correction was applied to partially reduce the discrepancies between the methods, but significant differences remained, likely due to violation of the fundamental assumption of a single-compartmental tissue model, and hence monoexponential decay.

CONCLUSION

R2' measurements are influenced significantly by the choice of method. Awareness of this issue is important when designing and interpreting studies that involve R2' measurements.

Effects of white matter microstructure on phase and susceptibility maps

Purpose

To investigate the effects on quantitative susceptibility mapping (QSM) and susceptibility tensor imaging (STI) of the frequency variation produced by the microstructure of white matter (WM).

Methods

The frequency offsets in a WM tissue sample that are not explained by the effect of bulk isotropic or anisotropic magnetic susceptibility, but rather result from the local microstructure, were characterized for the first time. QSM and STI were then applied to simulated frequency maps that were calculated using a digitized whole-brain, WM model formed from anatomical and diffusion tensor imaging data acquired from a volunteer. In this model, the magnitudes of the frequency contributions due to anisotropy and microstructure were derived from the results of the tissue experiments.

Results

The simulations suggest that the frequency contribution of microstructure is much larger than that due to bulk effects of anisotropic magnetic susceptibility. In QSM, the microstructure contribution introduced artificial WM heterogeneity. For the STI processing, the microstructure contribution caused the susceptibility anisotropy to be significantly overestimated.

Conclusion

Microstructure-related phase offsets in WM yield artifacts in the calculated susceptibility maps. If susceptibility mapping is to become a robust MRI technique, further research should be carried out to reduce the confounding effects of microstructure-related frequency contributions. Magn Reson Med 73:1258–1269, 2015. © 2014 Wiley Periodicals, Inc.

Quantitative magnetic susceptibility mapping without phase unwrapping using WASSR

The magnetic susceptibility of tissue within and around an image voxel affects the magnetic field and thus the local frequency in that voxel. Recently, it has been shown that spatial maps of frequency can be used to quantify local susceptibility if the contributions of surrounding tissue can be deconvolved. Currently, such quantitative susceptibility mapping (QSM) methods employ gradient recalled echo (GRE) imaging to measure spatial differences in the signal phase evolution as a function of echo time, from which frequencies can be deduced. Analysis of these phase images, however, is complicated by phase wraps, despite the availability and usage of various phase unwrapping algorithms. In addition, lengthy high-resolution GRE scanning often heats the magnet bore, causing the magnetic field to drift over several Hertz, which is on the order of the frequency differences between tissues. Here, we explore the feasibility of applying the WAter Saturation Shift Referencing (WASSR) method for 3D whole brain susceptibility imaging. WASSR uses direct saturation of water protons as a function of frequency irradiation offset to generate frequency maps without phase wraps, which can be combined with any image or spectroscopy acquisition. By utilizing a series of fast short-echo-time direct saturation images with multiple radiofrequency offsets, a frequency correction for field drift can be applied based on the individual image phases. Regions of interest were delineated with an automated atlas-based method, and the average magnetic susceptibilities calculated from frequency maps obtained from WASSR correlated well with those from the phase-based multi-echo GRE approach at 3T.

Longitudinal MR imaging of iron in multiple sclerosis: an imaging marker of disease

PURPOSE

To investigate the relationship between magnetic resonance (MR) imaging markers of iron content and disease severity in patients with multiple sclerosis (MS) over a 2-year period.

MATERIALS AND METHODS

This prospective study was approved by the local ethics committee, and written informed consent was obtained from all participants. Seventeen patients with MS and 17 control subjects were examined twice, 2 years apart, by using phase imaging and transverse relaxation (R2*) mapping at 4.7 T. Quantitative differences in iron content in deep gray matter between patients and control subjects were evaluated with repeated-measures multivariate analysis of variance separately for R2* mapping and phase imaging. Multiple regression analysis was used to evaluate correlations of MR imaging measures, both 2-year-difference and single-time measurements, to baseline disease severity.

RESULTS

R2* mapping using 2-year-difference measurements had the highest correlation to disease severity (r = 0.905, P < .001) compared with R2* mapping using single-time measurements (r = 0.560, P = .019) and phase imaging by using either single-time (r = 0.539, P = .026) or 2-year-difference (r = 0.644, P = .005) measurements. Significant increases in R2* occur during 2 years in the substantia nigra (P < .001) and globus pallidus (P = .035), which are both predictors of disease in regression analysis, in patients compared with control subjects. There were group differences in the substantia nigra, globus pallidus, pulvinar thalamus, thalamus, and caudate nucleus, compared with control subjects with R2* mapping (P < .05), and group differences in the caudate nucleus and pulvinar thalamus, compared with control subjects with phase imaging (P < .05).

CONCLUSION

There are significant changes in deep gray matter iron content in MS during 2 years measured with MR imaging, changes that are strongly related to physical disability. Longitudinal measurements may produce a higher correlation to disease severity compared with single-time measurements because baseline iron content of deep gray matter is variable among subjects.

Comparison of Magnetic Susceptibility Tensor and Diffusion Tensor of the Brain

Susceptibility tensor imaging (STI) provides a novel approach for noninvasive assessment of the white matter pathways of the brain. Using mouse brain ex vivo, we compared STI with diffusion tensor imaging (DTI), in terms of tensor values, principal tensor values, anisotropy values, and tensor orientations. Despite the completely different biophysical underpinnings, magnetic susceptibility tensors and diffusion tensors show many similarities in the tensor and principal tensor images, for example, the tensors perpendicular to the fiber direction have the highest gray-white matter contrast, and the largest principal tensor is along the fiber direction. Comparison to DTI fractional anisotropy, the susceptibility anisotropy provides much higher sensitivity to the chemical composition of the white matter, especially myelin. The high sensitivity can be further enhanced with the perfusion of ProHance, a gadolinium-based contrast agent. Regarding the tensor orientations, the direction of the largest principal susceptibility tensor agrees with that of diffusion tensors in major white matter fiber bundles. The STI fiber tractography can reconstruct the fiber pathways for the whole corpus callosum and for white matter fiber bundles that are in close contact but in different orientations. There are some differences between susceptibility and diffusion tensor orientations, which are likely due to the limitations in the current STI reconstruction. With the development of more accurate reconstruction methods, STI holds the promise for probing the white matter micro-architectures with more anatomical details and higher chemical sensitivity.

Gradient echo based fiber orientation mapping using R2* and frequency difference measurements

Fiber orientation mapping through diffusion tensor imaging (DTI) is a powerful MRI-based technique for visualising white matter (WM) microstructure in the brain. Although DTI provides a robust way to measure fiber orientation, it has some limitations linked to the use of EPI read-outs and long diffusion encoding periods, including relatively low spatial resolution. Development of alternative MRI-based methods for fiber orientation mapping is therefore valuable, in part to allow validation of DTI results. In this study, we used the orientation dependence of R2* (1/T2*) and frequency difference measurements to generate three dimensional maps of the fiber orientation in WM from multi-echo gradient-echo (GE) images acquired from post mortem brain tissue samples oriented at multiple angles to B0. Through analytical derivation and numerical simulation, the relationships connecting variations in R2* and frequency difference values to the angle between the underlying WM fiber orientation and the direction of B0 were characterised. High resolution 3D fiber orientation maps (FOM) were then formed by comparing R2* and frequency difference data, acquired with the sample at multiple orientations to the field, to generalised models based on the derived expressions for the angular dependence of each parameter. By comparing the resulting GE-based FOM with DTI-based FOM from the same tissue sample, we demonstrate that fiber orientation mapping based on gradient echo MRI has the potential to become an important tool for investigating microstructure in brain tissue.