The influence of white matter fibre orientation on MR signal phase and decay

The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

So You Want to Image Myelin Using MRI: Magnetic Susceptibility Source Separation for Myelin Imaging

In MRI, researchers have long endeavored to effectively visualize myelin distribution in the brain, a pursuit with significant implications for both scientific research and clinical applications. Over time, various methods such as myelin water imaging, magnetization transfer imaging, and relaxometric imaging have been developed, each carrying distinct advantages and limitations. Recently, an innovative technique named as magnetic susceptibility source separation has emerged, introducing a novel surrogate biomarker for myelin in the form of a diamagnetic susceptibility map. This paper comprehensively reviews this cutting-edge method, providing the fundamental concepts of magnetic susceptibility, susceptibility imaging, and the validation of the diamagnetic susceptibility map as a myelin biomarker that indirectly measures myelin content. Additionally, the paper explores essential aspects of data acquisition and processing, offering practical insights for readers. A comparison with established myelin imaging methods is also presented, and both current and prospective clinical and scientific applications are discussed to provide a holistic understanding of the technique. This work aims to serve as a foundational resource for newcomers entering this dynamic and rapidly expanding field.

Anisotropy of the R1/T2* value dependent on white matter fiber orientation with respect to the B0 field

The R1 (1/T1) map divided by the T2* map (R1/T2* map) draws attention as a high-resolution myelin-related map. However, both R1 and R2* (1/T2*) values demonstrate anisotropy dependent on the white matter (WM) fiber orientation with respect to the static magnetic (B0) field. Therefore, this study primarily aimed to investigate the comprehensive impact of these angular-dependent anisotropies on the R1/T2* value. This study enrolled 10 healthy human volunteers (age = 25 ± 1.3) on the 3.0 T MRI system. For R1/T2* map calculation, whole brain R1 and T2* maps were repeatedly obtained in three head tilt positions by magnetization-prepared two rapid gradient echoes and multiple spoiled gradient echo sequences, respectively. Afterward, all maps were spatially normalized and registered to the Johns Hopkins University WM atlas. R1/T2*, R1, and R2* values were binned for fiber orientation related to the B0 field, which was estimated from diffusion-weighted echo-planar imaging data with 3° intervals, to investigate angular-dependent anisotropies in vivo. A larger change in the R1/T2* value in the global WM region as a function of fiber orientation with respect to the B0 field was observed compared to the R1 and R2* values alone. The minimum R1/T2* value at the near magic-angle range was 18.86% lower than the maximum value at the perpendicular angle range. Furthermore, R1/T2* values in the corpus callosum tract and the right and left cingulum cingulate gyrus tracts changed among the three head tilt positions due to fiber orientation changes. In conclusion, the R1/T2* value demonstrates distinctive and complicated angular-dependent anisotropy indicating the trends of both R1 and R2* values and may provide supplemental information for detecting slight changes in the microstructure of myelin and axons.

Quantum dipole interactions and transient hydrogen bond orientation order in cells, cellular membranes and myelin sheath: Implications for MRI signal relaxation, anisotropy, and T1 magnetic field dependence

PURPOSE

Despite significant impact on the study of human brain, MRI lacks a theory of signal formation that integrates quantum interactions involving proton dipoles (a primary MRI signal source) with brain intricate cellular environment. The purpose of the present study is developing such a theory.

METHODS

We introduce the Transient Hydrogen Bond (THB) model, where THB-mediated quantum dipole interactions between water and protons of hydrophilic heads of amphipathic biomolecules forming cells, cellular membranes and myelin sheath serve as a major source of MR signal relaxation.

RESULTS

The THB theory predicts the existence of a hydrogen-bond-driven structural order of dipole-dipole connections within THBs as a primary factor for the anisotropy observed in MRI signal relaxation. We have also demonstrated that the conventional Lorentzian spectral density function decreases too fast at high frequencies to adequately capture the field dependence of brain MRI signal relaxation. To bridge this gap, we introduced a stretched spectral density function that surpasses the limitations of Lorentzian dispersion. In human brain, our findings reveal that at any time point only about 4% to 7% of water protons are engaged in quantum encounters within THBs. These ultra-short (2 to 3 ns), but frequent quantum spin exchanges lead to gradual recovery of magnetization toward thermodynamic equilibrium, that is, relaxation of MRI signal.

CONCLUSION

By incorporating quantum proton interactions involved in brain imaging, the THB approach introduces new insights on the complex relationship between brain tissue cellular structure and MRI measurements, thus offering a promising new tool for better understanding of brain microstructure in health and disease.

Myelin water imaging at 0.55 T using a multigradient-echo sequence

PURPOSE

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T.

METHODS

Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T.

RESULTS

The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate.

CONCLUSION

This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix

PURPOSE

To overcome some limitations of previous proton orientation-dependent transverse relaxation formalisms in human brain white matter (WM) by a generalized magic angle effect function.

METHODS

A cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R₂ and R₂^∗ were divided into isotropic R₂ⁱ and anisotropic parts, R₂^a ∗ f(α,Φ - ε₀), with α denoting an open angle and ε₀ an orientation (Φ) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared to prior models without ε₀ on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3 T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P ≤ 0.05.

RESULTS

Fit A significantly (P < 0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted ε₀ was in good agreement with calculated ε₀ from directional diffusivities. Compared with those from healthy adult, the fitted R₂ⁱ, R₂^a, and α from neonates were substantially reduced but ε₀ increased, consistent with developing myelination. Significant positive (R₂ⁱ) and negative (α and R₂^a) correlations were found with aging (demyelination) in elderly.

CONCLUSION

The developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in the human brain WM, thus shedding new light on myelin microstructural alterations at the molecular level.

Phase-shifted transverse relaxation orientation dependences in human brain white matter

This work aimed to demonstrate an essential phase shift ε 0 $$ {\varepsilon}_0 $$ for better quantifying R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ in human brain white matter (WM), and to further elucidate its origin related to the directional diffusivities from standard diffusion tensor imaging (DTI). ε 0 $$ {\varepsilon}_0 $$ was integrated into a proposed generalized transverse relaxation model for characterizing previously published R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ orientation dependence profiles in brain WM, and then comparisons were made with those without ε 0 $$ {\varepsilon}_0 $$ . It was theorized that anisotropic diffusivity direction ε $$ \varepsilon $$ was collinear with an axon fiber subject to all eigenvalues and eigenvectors from an apparent diffusion tensor. To corroborate the origin of ε 0 $$ {\varepsilon}_0 $$ , R 2 $$ {R}_2 $$ orientation dependences referenced by ε $$ \varepsilon $$ were compared with those referenced by the standard principal diffusivity direction Φ $$ \Phi $$ at b-values of 1000 and 2500 (s/mm² ). These R 2 $$ {R}_2 $$ orientation dependences were obtained from T 2 $$ {T}_2 $$ -weighted images (b = 0) of ultrahigh-resolution Connectome DTI datasets in the public domain. A normalized root-mean-square error ( NRMSE % $$ NRMSE\% $$ ) and an F $$ F $$ -test were used for evaluating curve-fittings, and statistical significance was considered to be a p of 0.05 or less. A phase-shifted model resulted in significantly reduced NRMSE % $$ NRMSE\% $$ compared with that without ε 0 $$ {\varepsilon}_0 $$ in quantifying various R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ profiles, both in vivo and ex vivo at multiple B 0 $$ {B}_0 $$ fields. The R 2 $$ {R}_2 $$ profiles based on Φ $$ \Phi $$ manifested a right-shifted phase ( ε 0 > 0 $$ {\varepsilon}_0>0 $$ ) at two b-values, while those based on ε $$ \varepsilon $$ became free from ε 0 $$ {\varepsilon}_0 $$ . For all phase-shifted R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ profiles, ε 0 $$ {\varepsilon}_0 $$ generally depended on the directional diffusivities by tan - 1 D ⊥ / D ∥ $$ {\tan}^{-1}\left({D}_{\perp }/{D}_{\parallel}\right) $$ , as predicted. In summary, a ubiquitous phase shift ε 0 $$ {\varepsilon}_0 $$ has been demonstrated as a prerequisite for better quantifying transverse relaxation orientation dependences in human brain WM. Furthermore, the origin of ε 0 $$ {\varepsilon}_0 $$ associated with the directional diffusivities from DTI has been elucidated. These findings could have a significant impact on interpretations of prior R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ datasets and on future research.

Comparison of quantitative susceptibility mapping methods for iron-sensitive susceptibility imaging at 7T: An evaluation in healthy subjects and patients with Huntington's disease

Quantitative susceptibility mapping (QSM) is a promising tool for investigating iron dysregulation in neurodegenerative diseases, including Huntington's disease (HD). Many diverse methods have been proposed to generate accurate and robust QSM images. In this study, we evaluated the performance of different dipole inversion algorithms for iron-sensitive susceptibility imaging at 7T on healthy subjects of a large age range and patients with HD. We compared an iterative least-squares-based method (iLSQR), iterative methods that use regularization, single-step approaches, and deep learning-based techniques. Their performance was evaluated by comparing: (1) deviations from a multiple-orientation QSM reference; (2) visual appearance of QSM maps and the presence of artifacts; (3) susceptibility in subcortical brain regions with age; (4) regional brain susceptibility with published postmortem brain iron quantification; and (5) susceptibility in HD-affected basal ganglia regions between HD subjects and healthy controls. We found that single-step QSM methods with either total variation or total generalized variation constraints (SSTV/SSTGV) and the single-step deep learning method iQSM generally provided the best performance in terms of correlation with iron deposition and were better at differentiating between healthy controls and premanifest HD individuals, while deep learning QSM methods trained with multiple-orientation susceptibility data created QSM maps that were most similar to the multiple orientation reference and with the best visual scores.

Orientation dependence of R2 relaxation in the newborn brain

In MRI the transverse relaxation rate, R₂ = 1/T₂^, shows dependence on the orientation of ordered tissue relative to the main magnetic field. In previous studies, orientation effects of R₂ relaxation in the mature brain's white matter have been found to be described by a susceptibility-based model of diffusion through local magnetic field inhomogeneities created by the diamagnetic myelin sheaths. Orientation effects in human newborn white matter have not yet been investigated. The newborn brain is known to contain very little myelin and is therefore expected to exhibit a decrease in orientation dependence driven by susceptibility-based effects. We measured R₂ orientation dependence in the white matter of human newborns. R₂ data were acquired with a 3D Gradient and Spin Echo (GRASE) sequence and fiber orientation was mapped with diffusion tensor imaging (DTI). We found orientation dependence in newborn white matter that is not consistent with the susceptibility-based model and is best described by a model of residual dipolar coupling. In the near absence of myelin in the newborn brain, these findings suggest the presence of residual dipolar coupling between rotationally restricted water molecules. This has important implications for quantitative imaging methods such as myelin water imaging, and suggests orientation dependence of R₂ as a potential marker in early brain development.

Static and dynamic BOLD fMRI components along white matter fibre tracts and their dependence on the orientation of the local diffusion tensor axis relative to the B0-field

Recent studies have reported functional MRI (fMRI) activation within cerebral white matter (WM) using blood-oxygenation-level-dependent (BOLD) contrast. Many blood vessels in WM run parallel to the fibre bundles, and other studies observed dependence of susceptibility contrast-based measures of blood volume on the local orientation of the fibre bundles relative to the magnetic field or B₀ axis. Motivated by this, we characterized the dependence of gradient-echo BOLD fMRI on fibre orientation (estimated by the local diffusion tensor) relative to the B₀ axis to test whether the alignment between bundles and vessels imparts an orientation dependence on resting-state BOLD fluctuations in the WM. We found that the baseline signal level of the T₂*-weighted data is 11% higher in voxels containing fibres parallel to B₀ than those containing perpendicular fibres, consistent with a static influence of either fibre or vessel orientation on local T₂* values. We also found that BOLD fluctuations in most bundles exhibit orientation effects expected from oxygenation changes, with larger amplitudes from voxels containing perpendicular fibres. Different magnitudes of this orientation effect were observed across the major WM bundles, with inferior fasciculus, corpus callosum and optic radiation exhibiting 14-19% higher fluctuations in voxels containing perpendicular compared to parallel fibres.

Relaxation anisotropy of quantitative MRI parameters in biological tissues

Quantitative MR relaxation parameters vary in the sensitivity to the orientation of the tissue in the magnetic field. In this study, the orientation dependence of multiple relaxation parameters was assessed in various tissues. Ex vivo samples of each tissue type were prepared either from bovine knee (tendon, cartilage) or mouse (brain, spinal cord, heart, kidney), and imaged at 9.4 T MRI with T1, T2, continuous wave (CW-) T1ρ, adiabatic T1ρ and T2ρ, and Relaxation along fictitious field (RAFF2-4) sequences at five different orientations with respect to the main magnetic field. Relaxation anisotropy of the measured parameters was quantified and compared. The highly ordered collagenous tissues, i.e. cartilage and tendon, presented the highest relaxation anisotropy for T2, CW-T1ρ with spin-lock power < 1 kHz, Ad-T2ρ and RAFF2-4. Maximally anisotropy was 75% in cartilage and 30% in tendon. T1 and adiabatic T1ρ did not exhibit observable anisotropy. In the other measured tissue types, anisotropy was overall less than 10% for all the parameters. The results confirm that highly ordered collagenous tissues have properties that induce very clearly observable relaxation anisotropy, whereas in other tissues the effect is not as prominent. Quantitative comparison of anisotropy of different relaxation parameters highlights the importance of sequence choice and design in MR imaging.

Regularized Asymmetric Susceptibility Tensor Imaging in the Human Brain in vivo

Susceptibility tensor imaging (STI) is a promising tool for studying orientation-dependent tissue magnetic susceptibility and for mapping white matter fiber orientations complementary to diffusion tensor imaging (DTI). However, the limited head rotation range within modern head coils for data acquisition makes in vivo STI reconstruction ill-conditioned. Conventional STI reconstruction method is usually vulnerable to noise and requires sufficiently large head rotations to solve this ill-conditioned inverse problem. In this study, based on the recently proposed asymmetric STI (aSTI) model, a new method termed aSTI+ was proposed to improve in vivo STI reconstruction by enforcing isotropic susceptibility tensor inside cerebrospinal fluid (CSF) and applying morphology constraint in white matter. Experimental results showed superior performance of the proposed method with reduced noise, improved tissue contrast and better fiber orientation estimation over previous methods. Thus aSTI+ may promote in vivo human brain STI studies on white matter and myelin-related brain diseases.

Magnetic susceptibility anisotropy in normal appearing white matter in multiple sclerosis from single-orientation acquisition

•

Orientation dependence of QSM is studied in a large cohort of MS patients.

•

Apparent magnetic susceptibility anisotropy (MSA) obtained from single-orientation QSM.

•

Apparent MSA found decreased in optic radiation (OR) of MS patients.

•

Apparent MSA decreases with lesion load in OR and with disease duration in splenium.

•

Negative apparent MSA observed in SLF indicates limitations of the proposed method.

Quantitative susceptibility mapping (QSM) has been successfully applied to study changes in deep grey matter nuclei as well as in lesional tissue, but its application to white matter has been complicated by the observed orientation dependence of gradient echo signal. The anisotropic susceptibility tensor is thought to be at the origin of this orientation dependence, and magnetic susceptibility anisotropy (MSA) derived from this tensor has been proposed as a marker of the state and integrity of the myelin sheath and may therefore be of particular interest for the study of demyelinating pathologies such as multiple sclerosis (MS). Reconstruction of the susceptibility tensor, however, requires repeated measurements with multiple head orientations, rendering the approach impractical for clinical applications.

In this study, we combined single-orientation QSM with fibre orientation information to assess apparent MSA in three white matter tracts, i.e., optic radiation (OR), splenium of the corpus callosum (SCC), and superior longitudinal fascicle (SLF), in two cohorts of 64 healthy controls and 89 MS patients. The apparent MSA showed a significant decrease in optic radiation in the MS cohort compared with healthy controls. It decreased in the MS cohort with increasing lesion load in OR and with disease duration in the splenium. All of this suggests demyelination in normal appearing white matter. However, the apparent MSA observed in the SLF pointed to potential systematic issues that require further exploration to realize the full potential of the presented approach. Despite the limitations of such single-orientation ROI-specific estimation, we believe that our clinically feasible approach to study degenerative changes in WM is worthy of further investigation.

Brain pathological changes during neurodegenerative diseases and their identification methods: How does QSM perform in detecting this process?

The presence of iron is essential for many biological processes in the body. But sometimes, for various reasons, the amount of iron deposition in different areas of the brain increases, which leads to problems related to the nervous system. Quantitative susceptibility mapping (QSM) is one of the newest magnetic resonance imaging (MRI)-based methods for assessing iron accumulation in target areas. This Narrative Review article aims to evaluate the performance of QSM compared to other methods of assessing iron deposition in the clinical field. Based on the results, we introduced related basic definitions, some neurodegenerative diseases, methods of examining iron deposition in these diseases, and their advantages and disadvantages. This article states that the QSM method can be introduced as a new, reliable, and non-invasive technique for clinical evaluations.

Learn Less, Infer More: Learning in the Fourier Domain for Quantitative Susceptibility Mapping

Quantitative susceptibility mapping (QSM) aims to evaluate the distribution of magnetic susceptibility from magnetic resonance phase measurements by solving the ill-conditioned dipole inversion problem. Removing the artifacts and preserving the anisotropy of tissue susceptibility simultaneously is still a challenge in QSM. To deal with this issue, a novel k-QSM network is proposed to resolve dipole inversion issues in QSM reconstruction. The k-QSM network converts the results obtained by truncated k-space division (TKD) into the Fourier domain as inputs. After passing through several convolutional and residual blocks, the ill-posed signals of TKD are corrected by making the network output close to the calculation of susceptibility through multiple orientation sampling (COSMOS)-labeled QSM. To evaluate the superiority of k-QSM, comparisons with several state-of-the-art methods are performed in terms of QSM artifacts removing, anisotropy preserving, generalization ability, and clinical applications. Compared to existing methods, the k-QSM achieves a 22.31% lower normalized root mean square error, 10.30% higher peak signal-to-noise ratio (PSNR), 33.10% lower high-frequency error norm, and 1.06% higher structural similarity. In addition, the orientation-dependent susceptibility variation obtained by k-QSM is significant, verifying that k-QSM has the ability to preserve susceptibility anisotropy. When the trained models are tested on the dataset from different centers, our k-QSM shows a strong generalization ability with the highest PSNR. Moreover, by comparing the susceptibility maps between healthy controls and drug addicts with different methods, we found the proposed k-QSM is more sensitive to the susceptibility abnormality in the patients. The proposed k-QSM method learns less—only to fix the ill-posed signals of TKD, but infers more—both COSMOS-like and anisotropy-preserving QSM results. Its generalization ability and great sensitivity to susceptibility changes can make it a potential method for distinguishing some diseases.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

In vivo assessment of anisotropy of apparent magnetic susceptibility in white matter from a single orientation acquisition

Multiple studies have reported a significant dependence of the effective transverse relaxation rate constant (R₂^*) and the phase of gradient-echo based (GRE) signal on the orientation of white matter fibres in the human brain. It has also been hypothesized that magnetic susceptibility, as obtained by single-orientation quantitative susceptibility mapping (QSM), exhibits such a dependence. In this study, we investigated this hypothesized relationship in a cohort of healthy volunteers. We show that R₂^* follows the predicted orientation dependence consistently across white matter regions, whereas the apparent magnetic susceptibility is related differently to fibre orientation across the brain and often in a complex non-monotonic manner. In addition, we explored the effect of fractional anisotropy measured by diffusion-weighted MRI on the strength of the orientation dependence and observed only a limited influence in many regions. However, with careful consideration of such an impact and the limitations imposed by the ill-posed nature of the dipole inversion process, it is possible to study magnetic susceptibility anisotropy in specific brain regions with a single orientation acquisition.

Quantitative magnetic resonance imaging towards clinical application in multiple sclerosis

Quantitative MRI provides biophysical measures of the microstructural integrity of the CNS, which can be compared across CNS regions, patients, and centres. In patients with multiple sclerosis, quantitative MRI techniques such as relaxometry, myelin imaging, magnetization transfer, diffusion MRI, quantitative susceptibility mapping, and perfusion MRI, complement conventional MRI techniques by providing insight into disease mechanisms. These include: (i) presence and extent of diffuse damage in CNS tissue outside lesions (normal-appearing tissue); (ii) heterogeneity of damage and repair in focal lesions; and (iii) specific damage to CNS tissue components. This review summarizes recent technical advances in quantitative MRI, existing pathological validation of quantitative MRI techniques, and emerging applications of quantitative MRI to patients with multiple sclerosis in both research and clinical settings. The current level of clinical maturity of each quantitative MRI technique, especially regarding its integration into clinical routine, is discussed. We aim to provide a better understanding of how quantitative MRI may help clinical practice by improving stratification of patients with multiple sclerosis, and assessment of disease progression, and evaluation of treatment response.

Sensitivity of fiber orientation dependent R 2 ∗ to temperature and post mortem interval

PURPOSE

R 2 ∗ imaging of brain white matter is well known for being sensitive to the orientation of nerve fibers with respect to the B₀ field of the MRI scanner. The goal of this study was to evaluate whether and to which extent fiber orientation dependent R 2 ∗ differs between in vivo and post mortem in situ examinations, and to investigate the influence of varying temperatures and post mortem intervals (PMI).

METHODS

Post mortem in situ and in vivo MRI scans were conducted at 3T. R 2 ∗ was acquired with a multi-echo gradient-echo sequence, and the orientation of white matter fibers was computed using diffusion tensor imaging (DTI). Fitting of the measured fiber orientation dependent R 2 ∗ was performed using three different formulations of a previously proposed model.

RESULTS

R 2 ∗ increased with increasing fiber angle for in vivo and post mortem in situ examinations, whereby the orientation dependency was lower post mortem. The different formulations of the fiber orientation model resulted in an identical fit, but showed large variations of the estimated parameters. The higher order orientation dependent R 2 ∗ components significantly decreased with decreasing temperature, while the orientation independent R 2 ∗ components showed no significant correlation with either temperature or PMI.

CONCLUSION

Although the mean diffusivity is strongly reduced post mortem, we could successfully estimate the fiber angle using DTI. Due to the strong correlation of the higher order orientation dependent R 2 ∗ components with temperature, the decreased R 2 ∗ fiber orientation dependency post mortem in situ might primarily be attributed to the lower brain temperature.

Magnetic Resonance Iron Imaging in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) results in progressive impairment of upper and lower motor neurons. Increasing evidence from both in vivo and ex vivo studies suggest that iron accumulation in the motor cortex is a neuropathological hallmark in ALS. An in vivo neuroimaging marker of iron dysregulation in ALS would be useful in disease diagnosis and prognosis. Magnetic resonance imaging (MRI), with its unique capability to generate a variety of soft tissue contrasts, provides opportunities to image iron distribution in the human brain with millimeter to sub-millimeter anatomical resolution. Conventionally, MRI T1-weighted, T2-weighted, and T2*-weighted images have been used to investigate iron dysregulation in the brain in vivo. Susceptibility weighted imaging has enhanced contrast for para-magnetic materials that provides superior sensitivity to iron in vivo. Recently, the development of quantitative susceptibility mapping (QSM) has realized the possibility of using quantitative assessments of magnetic susceptibility measures in brain tissues as a surrogate measurement of in vivo brain iron. In this review, we provide an overview of MRI techniques that have been used to investigate iron dysregulation in ALS in vivo. The potential uses, strengths, and limitations of these techniques in clinical trials, disease diagnosis, and prognosis are presented and discussed. We recommend further longitudinal studies with appropriate cohort characterization to validate the efficacy of these techniques. We conclude that quantitative iron assessment using recent advances in MRI including QSM holds great potential to be a sensitive diagnostic and prognostic marker in ALS. The use of multimodal neuroimaging markers in combination with iron imaging may also offer improved sensitivity in ALS diagnosis and prognosis that could make a major contribution to clinical care and treatment trials. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 3.

Myelin water imaging depends on white matter fiber orientation in the human brain

Purpose

The multi‐exponential T₂ decay of the MRI signal from cerebral white matter can be separated into short T₂ components related to myelin water and long T₂ components related to intracellular and extracellular water. In this study, we investigated to what degree the apparent myelin water fraction (MWF) depends on the angle between white matter fibers and the main magnetic field.

Methods

Maps of the apparent MWF were acquired using multi‐echo Carr‐Purcell‐Meiboom‐Gill and gradient‐echo spin‐echo sequences. The Carr‐Purcell‐Meiboom‐Gill sequence was acquired with a TR of 1073 ms, 1500 ms, and 2000 ms. The fiber orientation was mapped with DTI. By angle‐wise pooling the voxels across the brain’s white matter, orientation‐dependent apparent MWF curves were generated.

Results

We found that the apparent MWF varied between 25% and 35% across different fiber orientations. Furthermore, the selection of the TR influences the apparent MWF.

Conclusion

White matter fiber orientation induces a strong systematic bias on the estimation of the apparent MWF. This finding has implications for future research and the interpretation of MWI results in previously published studies.

Sensitivity to myelin using model-free analysis of the water resonance line-shape in postmortem mouse brain

PURPOSE

Dysmyelinating diseases are characterized by abnormal myelin formation and function. Such microstructural abnormalities in myelin have been demonstrated to produce measurable effects on the MR signal. This work examines these effects on measurements of voxel-wise, high-resolution water spectra acquired using a 3D echo-planar spectroscopic imaging (EPSI) pulse sequence from both postmortem fixed control mouse brains and a dysmyelination mouse brain model.

METHODS

Perfusion fixed, resected control (n = 5) and shiverer (n = 4) mouse brains were imaged using 3D-EPSI with 100 µm isotropic resolution. The free induction decay (FID) was sampled every 2.74 ms over 192 echoes, for a total sampling duration of 526.08 ms. Voxel-wise FIDs were Fourier transformed to produce water spectra with 1.9 Hz resolution. Spectral asymmetry was computed and compared between the two tissue types.

RESULTS

The water resonance is more asymmetrically broadened in the white matter of control mouse brain compared with dysmyelinated white matter. In control brain, this is modulated by and consistent with previously reported orientationally dependent effects of white matter relative to B₀ . Similar sensitivity to orientation is observed in dysmyelinated white matter as well; however, the magnitude of the resonance asymmetry is much lower across all directions.

CONCLUSION

Results demonstrate that components of the spectra are specifically differentially affected by myelin concentration. This suggests that water proton spectra may be sensitive to the presence of myelin, and as such, could serve as a MRI-based biomarker of dysmyelinating disease, free of mathematical models.

Influence of 7T GRE-MRI Signal Compartment Model Choice on Tissue Parameters

Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders. Ultra-high field magnetic resonance imaging (MRI) data obtained using a multi-echo gradient echo sequence have been shown to contain information on myelin, axonal, and extracellular compartments in tissue. Quantitative assessment of water fraction, relaxation time (T₂*), and frequency shift using multi-compartment models has been shown to be useful in studying white matter properties via specific tissue parameters. It remains unclear how tissue parameters vary with model selection based on 7T multiple echo time gradient-recalled echo (GRE) MRI data. We applied existing signal compartment models to the corpus callosum and investigated whether a three-compartment model can be reduced to two compartments and still resolve white matter parameters [i.e., myelin water fraction (MWF) and g-ratio]. We show that MWF should be computed using a three-compartment model in the corpus callosum, and the g-ratios obtained using three compartment models are consistent with previous reports. We provide results for other parameters, such as signal compartment frequency shifts.

Strong diffusion gradients allow the separation of intra- and extra-axonal gradient-echo signals in the human brain

The quantification of brain white matter properties is a key area of application of Magnetic Resonance Imaging (MRI), with much effort focused on using MR techniques to quantify tissue microstructure. While diffusion MRI probes white matter (WM) microstructure by characterising the sensitivity of Brownian motion of water molecules to anisotropic structures, susceptibility-based techniques probe the tissue microstructure by observing the effect of interaction between the tissue and the magnetic field. Here, we unify these two complementary approaches by combining ultra-strong (300 mT/m) gradients with a novel Diffusion-Filtered Asymmetric Spin Echo (D-FASE) technique. Using D-FASE we can separately assess the evolution of the intra- and extra-axonal signals under the action of susceptibility effects, revealing differences in the behaviour in different fibre tracts. We observed that the effective relaxation rate of the ASE signal in the corpus callosum decreases with increasing b-value in all subjects (from 17.1 ± 0.7 s⁻¹ at b = 0 s/mm² to 14.6 ± 0.7 s⁻¹ at b = 4800 s/mm²), while this dependence on b in the corticospinal tract is less pronounced (from 12.0± 1.1 s⁻¹ at b = 0s/mm² to 10.7 ± 0.5 s⁻¹ at b = 4800 s/mm²). Voxelwise analysis of the signal evolution with respect to b-factor and acquisition delay using a microscopic model demonstrated differences in gradient echo signal evolution between the intra- and extra-axonal pools.

Characterizing white matter fiber orientation effects on multi-parametric quantitative BOLD assessment of oxygen extraction fraction

Relative oxygen extraction fraction (rOEF) is a fundamental indicator of cerebral metabolic function. An easily applicable method for magnetic resonance imaging (MRI) based rOEF mapping is the multi-parametric quantitative blood oxygenation level dependent (mq-BOLD) approach with separate acquisitions of transverse relaxation times   T 2 * and T2 and dynamic susceptibility contrast (DSC) based relative cerebral blood volume (rCBV). Given that transverse relaxation and rCBV in white matter (WM) strongly depend on nerve fiber orientation, mq-BOLD derived rOEF is expected to be affected as well. To investigate fiber orientation related rOEF artefacts, we present a methodological study characterizing anisotropy effects of WM as measured by diffusion tensor imaging (DTI) on mq-BOLD in 30 healthy volunteers. Using a 3T clinical MRI-scanner, we performed a comprehensive correlation of all parameters (  T 2 * , T2, R 2 ' , rCBV, rOEF, where R 2 ' =1/  T 2 * -1/T2) with DTI-derived fiber orientation towards the main magnetic field (B0). Our results confirm strong dependencies of transverse relaxation and rCBV on the nerve fiber orientation towards B0, with anisotropy-driven variations up to 37%. Comparably weak orientation-dependent variations of mq-BOLD derived rOEF (3.8%) demonstrate partially counteracting influences of R 2 ' and rCBV effects, possibly suggesting applicability of rOEF as an oxygenation sensitive biomarker. However, unresolved issues warrant caution when applying mq-BOLD to WM.

Myelin water imaging and R2* mapping in neonates: Investigating R2* dependence on myelin and fibre orientation in whole brain white matter

R₂^* relaxation provides a semiquantitative method of detecting myelin, iron and white matter fibre orientation angles. Compared with standard histogram-based analyses, angle-resolved analysis of R₂^* has previously been shown to substantially improve the detection of subtle differences in the brain between healthy siblings of subjects with multiple sclerosis and unrelated healthy controls. Neonates, who are born with very little myelin and iron, and an underdeveloped connectome, provide researchers with an opportunity to investigate whether R₂^* is intimately linked with fibre-angle or myelin content as it is in adults, which may in future studies be explored as a potential white matter developmental biomarker. Five healthy adult volunteers (mean age [±SD] = 31.2 [±8.3] years; three males) were recruited from Vancouver, Canada. Eight term neonates (mean age = 38.6 ± 1.2 weeks; five males) were recruited from the Children's Hospital of Chongqing Medical University neonatal ward. All subjects were scanned on identical 3 T Philips Achieva scanners equipped with an eight-channel SENSE head coil and underwent a multiecho gradient echo scan, a 32-direction DTI scan and a myelin water imaging scan. For both neonates and adults, bin-averaged R₂^* variation across the brain's white matter was found to be best explained by fibre orientation. For adults, this represented a difference in R₂^* values of 3.5 Hz from parallel to perpendicular fibres with respect to the main magnetic field. In neonates, the fibre orientation dependency displayed a cosine wave shape, with a small R₂^* range of 0.4 Hz. This minor relationship in neonates provides further evidence for the key role myelin probably plays in creating this fibre orientation dependence later in life, but suggests limited clinical application in newborn populations. Future studies should investigate fibre-orientation dependency in infants in the first 5 years, when substantial myelin development occurs.

The role of iron and myelin in orientation dependent R2* of white matter

Brain myelin and iron content are important parameters in neurodegenerative diseases such as multiple sclerosis (MS). Both myelin and iron content influence the brain's R₂^* relaxation rate. However, their quantification based on R₂^* maps requires a realistic tissue model that can be fitted to the measured data. In structures with low myelin content, such as deep gray matter, R₂^* shows a linear increase with increasing iron content. In white matter, R₂^* is not only affected by iron and myelin but also by the orientation of the myelinated axons with respect to the external magnetic field. Here, we propose a numerical model which incorporates iron and myelin, as well as fibre orientation, to simulate R₂^* decay in white matter. Applying our model to fibre orientation-dependent in vivo R₂^* data, we are able to determine a unique solution of myelin and iron content in global white matter. We determine an averaged myelin volume fraction of 16.02 ± 2.07% in non-lesional white matter of patients with MS, 17.32 ± 2.20% in matched healthy controls, and 18.19 ± 2.98% in healthy siblings of patients with MS. Averaged iron content was 35.6 ± 8.9 mg/kg tissue in patients, 43.1 ± 8.3 mg/kg in controls, and 47.8 ± 8.2 mg/kg in siblings. All differences in iron content between groups were significant, while the difference in myelin content between MS patients and the siblings of MS patients was significant. In conclusion, we demonstrate that a model that combines myelin-induced orientation-dependent and iron-induced orientation-independent components is able to fit in vivo R₂^* data.

Three-dimensional MRI sequences in MS diagnosis and research

The most recent guidelines for magnetic resonance imaging (MRI) in multiple sclerosis (MS) recommend three-dimensional (3D) MRI sequences over their two-dimensional (2D) counterparts. This development has been made possible by advances in MRI scanner hardware and software. In this article, we review the 3D versions of conventional sequences, including T1-weighted, T2-weighted and fluid-attenuated inversion recovery (FLAIR), as well as more advanced scans, including double inversion recovery (DIR), FLAIR2, FLAIR*, phase-sensitive inversion recovery, and susceptibility weighted imaging (SWI).

Fibre-orientation dependent R1(=1/T1) relaxation in the brain: The role of susceptibility induced spin-lattice relaxation in the myelin water compartment

We have recently observed a dependence of the longitudinal relaxation rate, R₁, on the orientation of nerve fibres with respect to the main magnetic field. A similar dependence of R₂^∗ is long established and can be well explained by spin-dephasing in an inhomogeneous magnetic field induced by the susceptibility shift between myelin and water protons. The current study investigates if the same effect can also explain the R₁ dependence, neglecting a possible directional dependence of magnetisation transfer between solid myelin and myelin water. A molecular model of the myelin lipid bilayer was employed to simulate the susceptibility induced fields on a microscopic scale for the different nerve fibre orientations. The resulting simulated magnetic fields were used to calculate an orientation dependent relaxation offset, ΔR₁, based on both first-order perturbation theory and a simulation of the spin transition probabilities. Even though both methods yielded consistent orientation dependent relaxation offsets with a distribution that resembles the experimental data, the determined ΔR₁ values are too low to explain the reported R₁ angular dependency. Therefore, unlike R₂^∗, susceptibility induced spin flips can be excluded as a dominant source for the observed R₁ angular dependence.

Flexible proton density (PD) mapping using multi-contrast variable flip angle (VFA) data

Quantitative proton density (PD) maps measure the amount of free water, which is important for non-invasive tissue characterization in pathology and across lifespan. PD mapping requires the estimation and subsequent removal of factors influencing the signal intensity other than PD. These factors include the T1, T2* relaxation effects, transmit field inhomogeneities, receiver coil sensitivity profile (RP) and the spatially invariant factor that is required to scale the data. While the transmit field can be reliably measured, the RP estimation is usually based on image post-processing techniques due to limitations of its measurement at magnetic fields higher than 1.5 T. The post-processing methods are based on unified bias-field/tissue segmentation, fitting the sensitivity profile from images obtained with different coils, or on the linear relationship between T1 and PD. The scaling factor is derived from the signal within a specific tissue compartment or reference object. However, these approaches for calculating the RP and scaling factor have limitations particularly in severe pathology or over a wide age range, restricting their application. We propose a new approach for PD mapping based on a multi-contrast variable flip angle acquisition protocol and a data-driven estimation method for the RP correction and map scaling. By combining all the multi-contrast data acquired at different echo times, we are able to fully correct the MRI signal for T2* relaxation effects and to decrease the variance and the entropy of PD values within tissue class of the final map. The RP is determined from the corrected data applying a non-parametric bias estimation, and the scaling factor is based on the median intensity of an external calibration object. Finally, we compare the signal intensity and homogeneity of the multi-contrast PD map with the well-established effective PD (PD*) mapping, for which the RP is based on concurrent bias field estimation and tissue classification, and the scaling factor is estimated from the mean white matter signal. The multi-contrast PD values homogeneity and accuracy within the cerebrospinal fluid (CSF) and deep brain structures are increased beyond that obtained using PD* maps. We demonstrate that the multi-contrast RP approach is insensitive to anatomical or a priori tissue information by applying it in a patient with extensive brain abnormalities and for whole body PD mapping in post-mortem foetal imaging.

The effects of variations in tissue microstructure from postmortem rat brain on the asymmetry of the water proton resonance

PURPOSE

This work was performed to investigate the effects of tissue microstructure from postmortem rat brain on the shape of the water proton spectrum.

METHODS

Perfusion-fixed, resected rat brains (N = 4) were imaged at 9.4T. 3D DTI and 3D echo-planar spectroscopic imaging (EPSI) data were acquired with 150 μm isotropic resolution. DTI data were acquired over 60 directions with b = 3000 s/mm² . Water spectra were produced from EPSI data acquired over 128 echoes, with 2.9 Hz spectral resolution. A voxel-wise metric reflecting spectral asymmetry about the peak of the resonance was computed and compared with orientation estimates from DTI data by fitting data with the susceptibility anisotropy model.

RESULTS

Asymmetric broadening of the water resonance was computed for mixed populations of grey and/or white matter as determined by thresholding the fractional anisotropy. Asymmetry was shown to be differentially affected by tract orientation relative to B₀ in high FA voxels, whereas low FA voxels exhibited little sensitivity. Anatomic structures in the hippocampus were also found to produce distinct changes in the water resonance.

CONCLUSION

Present results demonstrate that structural variations in tissue architecture cause characteristic, reproducible changes in the water resonance shape. This suggests that water spectra are sensitive to cytoarchitectural variations in brain tissue.

Lorentzian effects in magnetic susceptibility mapping of anisotropic biological tissues

The ultimate goal of MRI is to provide information on biological tissue microstructure and function. Quantitative Susceptibility Mapping (QSM) is one of the newer approaches for studying tissue microstructure by means of measuring phase of Gradient Recalled Echo (GRE) MRI signal. The fundamental question in the heart of this approach is: what is the relationship between the net phase/frequency of the GRE signal from an imaging voxel and the underlying tissue microstructure at the cellular and sub-cellular levels? In the presence of external magnetic field, biological media (e.g. cells, cellular components, blood) become magnetized leading to the MR signal frequency shift that is affected not only by bulk magnetic susceptibility but by the local cellular environment as well. The latter effect is often termed the Lorentzian contribution to the frequency shift. Evaluating the Lorentzian contribution - one of the most intriguing and challenging problems in this field - is the main focus of this review. While the traditional approach to this problem is based on introduction of an imaginary Lorentzian cavity, a more rigorous treatment was proposed recently based on a statistical approach and a direct solution of the Maxwell equations. This approach, termed the Generalized Lorentzian Tensor Approach (GLTA), is especially fruitful for describing anisotropic biological media. The GLTA adequately accounts for two types of anisotropy: anisotropy of magnetic susceptibility and tissue structural anisotropy (e.g., cylindrical axonal bundles in white matter). In the framework of the GLTA the frequency shift due to the local environment is described in terms of the Lorentzian tensor L̂ which can have a substantially different structure than the susceptibility tensor χ̂. While the components of χ̂ are compartmental susceptibilities "weighted" by their volume fractions, the components of L̂ are additionally weighted by specific numerical factors depending on cellular geometrical symmetry. In addition to describing the GLTA that is a phenomenological approach largely based on considering the system symmetry, we also briefly discuss a microscopic approaches to the problem that are based on modeling of the MR signal in different regimes (i.e. static dephasing vs. motion narrowing) and in different cellular environments (e.g., accounting for WM microstructure).

Quantifying MRI frequency shifts due to structures with anisotropic magnetic susceptibility using pyrolytic graphite sheet

Magnetic susceptibility is an important source of contrast in magnetic resonance imaging (MRI), with spatial variations in the susceptibility of tissue affecting both the magnitude and phase of the measured signals. This contrast has generally been interpreted by assuming that tissues have isotropic magnetic susceptibility, but recent work has shown that the anisotropic magnetic susceptibility of ordered biological tissues, such as myelinated nerves and cardiac muscle fibers, gives rise to unexpected image contrast. This behavior occurs because the pattern of field variation generated by microstructural elements formed from material of anisotropic susceptibility can be very different from that predicted by modelling the effects in terms of isotropic susceptibility. In MR images of tissue, such elements are manifested at a sub-voxel length-scale, so the patterns of field variation that they generate cannot be directly visualized. Here, we used pyrolytic graphite sheet which has a large magnetic susceptibility anisotropy to form structures of known geometry with sizes large enough that the pattern of field variation could be mapped directly using MRI. This allowed direct validation of theoretical expressions describing the pattern of field variation from anisotropic structures with biologically relevant shapes (slabs, spherical shells and cylindrical shells).

Studying brain microstructure with magnetic susceptibility contrast at high-field

A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.

Quantitative MRI provides markers of intra-, inter-regional, and age-related differences in young adult cortical microstructure

Measuring the structural composition of the cortex is critical to understanding typical development, yet few investigations in humans have charted markers in vivo that are sensitive to tissue microstructural attributes. Here, we used a well-validated quantitative MR protocol to measure four parameters (R₁, MT, R₂*, PD*) that differ in their sensitivity to facets of the tissue microstructural environment (R₁, MT: myelin, macromolecular content; R₂*: myelin, paramagnetic ions, i.e., iron; PD*: free water content). Mapping these parameters across cortical regions in a young adult cohort (18–39 years, N = 93) revealed expected patterns of increased macromolecular content as well as reduced tissue water content in primary and primary adjacent cortical regions. Mapping across cortical depth within regions showed decreased expression of myelin and related processes – but increased tissue water content – when progressing from the grey/white to the grey/pial boundary, in all regions. Charting developmental change in cortical microstructure cross-sectionally, we found that parameters with sensitivity to tissue myelin (R₁ & MT) showed linear increases with age across frontal and parietal cortex (change 0.5–1.0% per year). Overlap of robust age effects for both parameters emerged in left inferior frontal, right parietal and bilateral pre-central regions. Our findings afford an improved understanding of ontogeny in early adulthood and offer normative quantitative MR data for inter- and intra-cortical composition, which may be used as benchmarks in further studies.

•

We mapped multi-parameter maps (MPMs) across and within cortical regions.

•

We charted age effects (ages 18–39) on myelin and related processes.

•

MPMs sensitive to myelin (R₁, MT) showed elevated values in primary areas over most cortical depths.

•

R₂* map foci tended to overlap MPMs sensitive to myelin (R₁, MT).

•

R₁ and MT increased with age (0.5–1.0% per year) at mid-depth in frontal and parietal cortex.

The impact of white matter fiber orientation in single-acquisition quantitative susceptibility mapping

The aim of this work was to assess the impact of tissue structural orientation on quantitative susceptibility mapping (QSM) reliability, and to provide a criterion to identify voxels in which measures of magnetic susceptibility (χ) are most affected by spatial orientation effects. Four healthy volunteers underwent 7-T magnetic resonance imaging (MRI). Multi-echo, gradient-echo sequences were used to obtain quantitative maps of frequency shift (FS) and χ. Information from diffusion tensor imaging (DTI) was used to investigate the relationship between tissue orientation and FS measures and QSM. After sorting voxels on the basis of their fractional anisotropy (FA), the variations in FS and χ values over tissue orientation were measured. Using a K-means clustering algorithm, voxels were separated into two groups depending on the variability of measures within each FA interval. The consistency of FS and QSM values, observed at low FA, was disrupted for FA > 0.6. The standard deviation of χ measured at high FA (0.0103 ppm) was nearly five times that at low FA (0.0022 ppm). This result was consistent through data across different head positions and for different brain regions considered separately, which confirmed that such behavior does not depend on structures with different bulk susceptibility oriented along particular angles. The reliability of single-orientation QSM anticorrelates with local FA. QSM provides replicable values with little variability in brain regions with FA < 0.6, but QSM should be interpreted cautiously in major and coherent fiber bundles, which are strongly affected by structural anisotropy and magnetic susceptibility anisotropy.

Phase imaging of axonal integrity of cranial corticospinal tract in experimental spinal cord injury at 9.4T

It was known that Spinal cord injury (SCI) leads to neuronal cell death, axonal damage and demyelination, and often results in severe functional loss. Noninvasive imaging of axonal integrity may hold great importance in clinic. This study evaluated the role of phase imaging in detection of axonal integrity of cranial corticospinal tract after spinal cord injury. Quantitative MR phase imaging and diffusion tensor imaging were performed on a group of rats prior to and after spinal cord injury at 9.4 T in this study. The phase contrast decreased mainly in contralateral white matter containing the corticospinal tract (CST), i.e., pyramid by 46%, two weeks postinjury compared to preinjury levels. Diffusion measurements further showed fractional anisotropy (FA) increased and radial diffusivity (λ⊥) decreased significantly, and axial diffusivity (λ∥) remained unchanged in these regions confirmed intact fiber tracts two weeks after injury, precluding the possibility of axonal damage and potential contribution to the shifted susceptibility. It was concluded that phase imaging is a potential endogenous biomarker for brain axonal integrity after spinal cord injury. This novel tool could potentially be applied to diagnose brain pathologies in patients who have spinal cord injury.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.

Effects of biological tissue structural anisotropy and anisotropy of magnetic susceptibility on the gradient echo MRI signal phase: theoretical background

Quantitative susceptibility mapping is a potentially powerful technique for mapping tissue magnetic susceptibility from gradient recalled echo (GRE) MRI signal phase. In this review, we present up-to-date theoretical developments in analyzing the relationships between GRE signal phase and the underlying tissue microstructure and magnetic susceptibility at the cellular level. Two important phenomena contributing to the GRE signal phase are at the focus of this review - tissue structural anisotropy (e.g. cylindrical axonal bundles in white matter) and magnetic susceptibility anisotropy. One of the most intriguing and challenging problems in this field is calculating the so-called Lorentzian contribution to the phase shift induced by the local environment - magnetized tissue structures that have dimensions smaller than the imaging voxel (e.g. cells, cellular components, blood capillaries). In this review, we briefly discuss a "standard" approach to this problem, based on introduction of an imaginary Lorentzian cavity, as well as a more recent method - the generalized Lorentzian tensor approach (GLTA) - that is based on a statistical approach and a direct solution of the magnetostatic Maxwell equations. The latter adequately 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 the components of χ^ are compartmental susceptibilities "weighted" by their volume fractions, the components of L^ are weighted by specific numerical factors depending on tissue geometrical microsymmetry. In multi-compartment structures, the components of the Lorentzian tensor also depend on the compartmental relaxation properties, hence the MR pulse sequence settings. Copyright © 2016 John Wiley & Sons, Ltd.

Heterogeneous anisotropic magnetic susceptibility of the myelin-water layers causes local magnetic field perturbations in axons

One goal of MRI is to determine the myelin water fraction in neural tissue. One approach is to measure the reduction in T₂ * arising from microscopic perturbations in the magnetic field caused by heterogeneities in the magnetic susceptibility of myelin. In this paper, analytic expressions for the induced magnetic field distribution are derived within and around an axon, assuming that the myelin susceptibility is anisotropic. Previous models considered the susceptibility to be piecewise continuous, whereas this model considers a sinusoidally varying susceptibility. Many conclusions are common in both models. When the magnetic field is applied perpendicular to the axon, the magnetic field in the intraaxonal space is uniformly perturbed, the magnetic field in the myelin sheath oscillates between the lipid and water layers, and the magnetic field in the extracellular space just outside the myelin sheath is heterogeneous. These field heterogeneities cause the spins to dephase, shortening T₂ *. When the magnetic field is applied along the axon, the field is homogeneous within water-filled regions, including between lipid layers. Therefore the spins do not dephase and the magnetic susceptibility has no effect on T₂ *. Generally, the response of an axon is given as the superposition of these two contributions. The sinusoidal model uses a different set of approximations compared with the piecewise model, so their common predictions indicate that the models are not too sensitive to the details of the myelin-water distribution. Other predictions, such as the sensitivity to water diffusion between myelin and water layers, may highlight differences between the two approaches. Copyright © 2016 John Wiley & Sons, Ltd.

Susceptibility tensor imaging (STI) of the brain

Susceptibility tensor imaging (STI) is a recently developed MRI technique that allows quantitative determination of orientation-independent magnetic susceptibility parameters from the dependence of gradient echo signal phase on the orientation of biological tissues with respect to the main magnetic field. By modeling the magnetic susceptibility of each voxel as a symmetric rank-2 tensor, individual magnetic susceptibility tensor elements as well as the mean magnetic susceptibility and magnetic susceptibility anisotropy can be determined for brain tissues that would still show orientation dependence after conventional scalar-based quantitative susceptibility mapping to remove such dependence. Similar to diffusion tensor imaging, STI allows mapping of brain white matter fiber orientations and reconstruction of 3D white matter pathways using the principal eigenvectors of the susceptibility tensor. In contrast to diffusion anisotropy, the main determinant factor of the susceptibility anisotropy in brain white matter is myelin. Another unique feature of the susceptibility anisotropy of white matter is its sensitivity to gadolinium-based contrast agents. Mechanistically, MRI-observed susceptibility anisotropy is mainly attributed to the highly ordered lipid molecules in the myelin sheath. STI provides a consistent interpretation of the dependence of phase and susceptibility on orientation at multiple scales. This article reviews the key experimental findings and physical theories that led to the development of STI, its practical implementations, and its applications for brain research. Copyright © 2016 John Wiley & Sons, Ltd.

Mechanisms of T2 * anisotropy and gradient echo myelin water imaging

In MRI, structurally aligned molecular or micro-organization (e.g. axonal fibers) can be a source of substantial signal variations that depend on the structural orientation and the applied magnetic field. This signal anisotropy gives us a unique opportunity to explore information that exists at a resolution several orders of magnitude smaller than that of typical MRI. In this review, one of the signal anisotropies, T₂ * anisotropy in white matter, and a related imaging method, gradient echo myelin water imaging (GRE-MWI), are explored. The T₂ * anisotropy has been attributed to isotropic and anisotropic magnetic susceptibility of myelin and compartmentalized microstructure of white matter fibers (i.e. axonal, myelin, and extracellular space). The susceptibility and microstructure create magnetic frequency shifts that change with the relative orientation of the fiber and the main magnetic field, generating the T₂ * anisotropy. The resulting multi-component magnitude decay and nonlinear phase evolution have been utilized for GRE-MWI, assisting in resolving the signal fraction of the multiple compartments in white matter. The GRE-MWI method has been further improved by signal compensation techniques including physiological noise compensation schemes. The T₂ * anisotropy and GRE-MWI provide microstructural information on a voxel (e.g. fiber orientation and tissue composition), and may serve as sensitive biomarkers for microstructural changes in the brain. Copyright © 2016 John Wiley & Sons, Ltd.

Iron quantification with susceptibility

Iron is an essential trace element involved in a variety of biological mechanisms in the human body. Disturbances of iron homeostasis have been observed in several inflammatory and degenerative diseases, which have raised strong interest in non-invasive iron mapping techniques. Numerous MRI techniques have been proposed so far, mostly based on the field changes induced by the magnetic properties of iron. Each of these approaches has a specific sensitivity for iron and its microstructural environment. Quantitative susceptibility mapping is the latest development and provides a direct measure of bulk susceptibility. However, field changes induced by iron are not always directly related to the concentration of iron, but rather reflect the structure of iron compounds and its cellular distribution. This review provides an overview of the most relevant iron compounds in the human body, their magnetic properties and their cellular distribution. In addition, MRI methods based on direct or indirect susceptibility changes are presented and discussed with respect to technical aspects and clinical applicability. Copyright © 2016 John Wiley & Sons, Ltd.