A comparison of phase imaging and quantitative susceptibility mapping in the imaging of multiple sclerosis lesions at ultrahigh field.

OBJECTIVE: The aim of this study was to compare the use of high-resolution phase and QSM images acquired at ultra-high field in the investigation of multiple sclerosis (MS) lesions with peripheral rings, and to discuss their usefulness for drawing inferences about underlying tissue composition. MATERIALS AND METHODS: Thirty-nine Subjects were scanned at 7 T, using 3D T 2*-weighted and T 1-weighted sequences. Phase images were then unwrapped and filtered, and quantitative susceptibility maps were generated using a thresholded k-space division method. Lesions were compared visually and using a 1D profiling algorithm. RESULTS: Lesions displaying peripheral rings in the phase images were identified in 10 of the 39 subjects. Dipolar projections were apparent in the phase images outside of the extent of several of these lesions; however, QSM images showed peripheral rings without such projections. These projections appeared ring-like in a small number of phase images where no ring was observed in QSM. 1D profiles of six well-isolated example lesions showed that QSM contrast corresponds more closely to the magnitude images than phase contrast. CONCLUSIONS: Phase images contain dipolar projections, which confounds their use in the investigation of tissue composition in MS lesions. Quantitative susceptibility maps correct these projections, providing insight into the composition of MS lesions showing peripheral rings.

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.

Fiber orientation-dependent white matter contrast in gradient echo MRI.

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

Global intravascular and local hyperoxia contrast phase-based blood oxygenation measurements.

The measurement of venous cerebral blood oxygenation (Yv) has potential applications in the study of patient groups where oxygen extraction and/or metabolism are compromised. It is also useful for fMRI studies to assess the stimulus-induced changes in Yv, particularly since basal Yv partially accounts for inter-subject variation in the haemodynamic response to a stimulus. A range of MRI-based methods of measuring Yv have been developed recently. Here, we use a method based on the change in phase in the MR image arising from the field perturbation caused by deoxygenated haemoglobin in veins. We build on the existing phase based approach (Method I), where Yv is measured in a large vein (such as the superior sagittal sinus) based on the field shift inside the vein with assumptions as to the vein's shape and orientation. We demonstrate two novel modifications which address limitations of this method. The first modification (Method II), maps the actual form of the vein, rather than assume a given shape and orientation. The second modification (Method III) uses the intra and perivascular phase change in response to a known change in Yv on hyperoxia to measure normoxic Yv in smaller veins. Method III can be applied to veins whose shape, size and orientation are not accurately known, thus allowing more localised measures of venous oxygenation. Results demonstrate that the use of an overly fine spatial filter caused an overestimation in Yv for Method I, whilst the measurement of Yv using Method II was less sensitive to this bias, giving Yv = 0.62 ± 0.03. Method III was applied to mapping of Yv in local veins across the brain, yielding a distribution of values with a mode of Yv = 0.661 ± 0.008.

Functional quantitative susceptibility mapping (fQSM).

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) is a powerful technique, typically based on the statistical analysis of the magnitude component of the complex time-series. Here, we additionally interrogated the phase data of the fMRI time-series and used quantitative susceptibility mapping (QSM) in order to investigate the potential of functional QSM (fQSM) relative to standard magnitude BOLD fMRI. High spatial resolution data (1mm isotropic) were acquired every 3 seconds using zoomed multi-slice gradient-echo EPI collected at 7 T in single orientation (SO) and multiple orientation (MO) experiments, the latter involving 4 repetitions with the subject's head rotated relative to B0. Statistical parametric maps (SPM) were reconstructed for magnitude, phase and QSM time-series and each was subjected to detailed analysis. Several fQSM pipelines were evaluated and compared based on the relative number of voxels that were coincidentally found to be significant in QSM and magnitude SPMs (common voxels). We found that sensitivity and spatial reliability of fQSM relative to the magnitude data depended strongly on the arbitrary significance threshold defining "activated" voxels in SPMs, and on the efficiency of spatio-temporal filtering of the phase time-series. Sensitivity and spatial reliability depended slightly on whether MO or SO fQSM was performed and on the QSM calculation approach used for SO data. Our results present the potential of fQSM as a quantitative method of mapping BOLD changes. We also critically discuss the technical challenges and issues linked to this intriguing new technique.

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.

Increased iron accumulation occurs in the earliest stages of demyelinating disease: an ultra-high field susceptibility mapping study in Clinically Isolated Syndrome.

OBJECTIVE: To determine, using ultra-high field magnetic resonance imaging (MRI), whether changes in iron content occur in the earliest phases of demyelinating disease, by quantifying the magnetic susceptibility of deep grey matter structures in patients with Clinically Isolated Syndrome (CIS) that is suggestive of multiple sclerosis (MS), as compared with age-matched healthy subjects. METHODS: We compared 19 CIS patients to 20 age-matched, healthy controls. Scanning of the study subjects was performed on a 7T Philips Achieva system, using a 3-dimensional, T2*-weighted gradient echo acquisition. Phase data were first high-pass filtered, using a dipole fitting method, and then inverted to produce magnetic susceptibility maps. Region of interest (ROI) analysis was used to estimate magnetic susceptibility values for deep grey matter structures (caudate nucleus, putamen, globus pallidus, the thalamus and its pulvinar). RESULTS: Significantly increased relative susceptibilities were found in the CIS group, compared with controls, for the caudate nucleus (p = < 0.01), putamen (p < 0.01), globus pallidus (p < 0.01) and pulvinar (p < 0.05). We found no significant nor consistent trends in the relationship between susceptibility and age for either the study controls or CIS patients, in any ROI (r(2) < 0.5; p > 0.05). In CIS patients, the time elapsed since the clinical event and the Expanded Disability Status Scale (EDSS) scores were not correlated with iron levels in any ROI (r(2) < 0.5; p > 0.05); however, a moderate correlation (r(2) = 0.3; p < 0.01) was found between the T1 lesion load and the mean susceptibility of the caudate nucleus. CONCLUSION: CIS patients showed an increased iron accumulation, as measured using susceptibility mapping of the deep grey matter, suggesting that iron changes did occur at the earlier stages of CIS disease.

Calibrated BOLD using direct measurement of changes in venous oxygenation.

Calibration of the BOLD signal is potentially of great value in providing a closer measure of the underlying changes in brain function related to neuronal activity than the BOLD signal alone, but current approaches rely on an assumed relationship between cerebral blood volume (CBV) and cerebral blood flow (CBF). This is poorly characterised in humans and does not reflect the predominantly venous nature of BOLD contrast, whilst this relationship may vary across brain regions and depend on the structure of the local vascular bed. This work demonstrates a new approach to BOLD calibration which does not require an assumption about the relationship between cerebral blood volume and cerebral blood flow. This method involves repeating the same stimulus both at normoxia and hyperoxia, using hyperoxic BOLD contrast to estimate the relative changes in venous blood oxygenation and venous CBV. To do this the effect of hyperoxia on venous blood oxygenation has to be calculated, which requires an estimate of basal oxygen extraction fraction, and this can be estimated from the phase as an alternative to using a literature estimate. Additional measurement of the relative change in CBF, combined with the blood oxygenation change can be used to calculate the relative change in CMRO(2) due to the stimulus. CMRO(2) changes of 18 ± 8% in response to a motor task were measured without requiring the assumption of a CBV/CBF coupling relationship, and are in agreement with previous approaches.

High resolution magnetic susceptibility mapping of the substantia nigra in Parkinson's disease.

PURPOSE: To determine if tissue magnetic susceptibility is a more direct marker of tissue iron content than other MR markers of iron. This study presents the first quantitative, in vivo measurements of the susceptibility of the substantia nigra in patients with Parkinson's disease. MATERIALS AND METHODS: Nine patients and 11 controls were studied at 7 Tesla. Susceptibility maps were created by inverting the filtered phase maps associated with T2* weighted images. RESULTS: On average, patients showed an increase in susceptibility of the pars compacta compared with controls, which correlates with the predicted increase in brain iron in Parkinson's disease. A rostral-caudal gradient in susceptibility was also observed in controls and patients. CONCLUSION: Susceptibility mapping may provide a new tool for studying the development of Parkinson's disease.