Origin of B0 orientation dependent R2(*) (=1/T2(*)) in white matter

χ-sepnet: Deep Neural Network for Magnetic Susceptibility Source Separation

Magnetic susceptibility source separation (χ-separation), an advanced quantitative susceptibility mapping (QSM) method, enables the separate estimation of paramagnetic and diamagnetic susceptibility source distributions in the brain. Similar to QSM, it requires solving the ill-conditioned problem of dipole inversion, suffering from so-called streaking artifacts. Additionally, the method utilizes reversible transverse relaxation (R 2 ' = R 2 * - R 2 $$ {R}_2^{\prime }={R}_2^{\ast }-{R}_2 $$ ) to complement frequency shift information for estimating susceptibility source concentrations, requiring time-consuming data acquisition forR 2 $$ {R}_2 $$ (e.g., multi-echo spin-echo) in addition to multi-echo GRE data forR 2 * $$ {R}_2^{\ast } $$ . To address these challenges, we develop a new deep learning network, χ-sepnet, and propose two deep learning-based susceptibility source separation pipelines, χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ for inputs with multi-echo GRE and multi-echo spin-echo (or turbo spin-echo) and χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ for input with multi-echo GRE only. The neural network is trained using multiple head orientation data that provide streaking artifact-free labels, generating high-quality χ-separation maps. The evaluation of the pipelines encompasses both qualitative and quantitative assessments in healthy subjects, and visual inspection of lesion characteristics in multiple sclerosis patients. The susceptibility source-separated maps of the proposed pipelines delineate detailed brain structures with substantially reduced artifacts compared to those from the conventional regularization-based reconstruction methods. In quantitative analysis, χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ achieves the best outcomes followed by χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ , outperforming the conventional methods. When the lesions of multiple sclerosis patients are classified into subtypes, most lesions are identified as the same subtype in the maps from χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ and χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ (paramagnetic susceptibility: 99.6% and diamagnetic susceptibility: 98.4%; both out of 250 lesions). The χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ pipeline, which only requires multi-echo GRE data, has demonstrated its potential to offer broad clinical and scientific applications, although further evaluations for various diseases and pathological conditions are necessary.

In-vivo high-resolution χ-separation at 7T

A recently introduced quantitative susceptibility mapping (QSM) technique, χ-separation, offers the capability to separate paramagnetic (χpara) and diamagnetic (χdia) susceptibility distribution within the brain. In-vivo high-resolution mapping of iron and myelin distribution, estimated by χ-separation, could provide a deeper understanding of brain substructures, assisting the investigation of their functions and alterations. This can be achieved using 7T MRI, which benefits from a high signal-to-noise ratio and susceptibility effects. However, applying χ-separation at 7T presents difficulties due to the requirement of an R2 map, coupled with issues such as high specific absorption rate (SAR), large B1 transmit field inhomogeneities, and prolonged scan time. To address these challenges, we developed a novel deep neural network, R2PRIMEnet7T, designed to convert a 7T R2* map into a 3T R2' map. Building on this development, we present a new pipeline for χ-separation at 7T, enabling us to generate high-resolution χ-separation maps from multi-echo gradient-echo data. The proposed method is compared with alternative pipelines, such as an end-to-end network and linearly-scaled R2', and is validated against χ-separation maps at 3T, demonstrating its accuracy. The 7T χ-separation maps generated by the proposed method exhibit similar contrasts to those from 3T, while 7T high-resolution maps offer enhanced clarity and detail. Quantitative analysis confirms that the proposed method surpasses the alternative pipelines. The proposed method results well delineate the detailed brain structures associated with iron and myelin. This new pipeline holds promise for analyzing iron and myelin concentration changes in various neurodegenerative diseases through precise structural examination.

Histological Validation of Multi-Echo Gradient Echo (MGRE)-Derived Myelin Water Fraction (MWF) at 9.4 T and the Influence of Orientation on Quantification

Myelin is essential in the nervous system of mammals. As the location and degree of myelin loss can reflect varied pathophysiological status, noninvasive measurement of myelin is of high importance. The magnetic resonance imaging (MRI) technique of myelin water fraction (MWF) derived from multi-echo gradient echo (MGRE) sequence is a promising tool for the quantification of myelin content due to the low specific absorption rate (SAR) compared with the spin-echo sequence, time efficiency, and wide availability. Yet to our knowledge, MGRE-derived MWF has never been quantitatively validated with histology. The main objective of this study was to quantitatively validate the MRI findings by referencing the myelin histology using a rat model. As a second objective, we investigated how the orientation of white matter fibers with respect to the static B0 field impacted both the apparent transverse relaxation rate (R2* = 1/T2*) and the derived MWF. Moreover, MWF is known to change with age; thus, we compared rat brains of different ages. The orientation effect of MWF in a clinical setting was studied using 3 T human data. Twenty ex vivo rat brains with different ages and three healthy volunteers were scanned on a 9.4 T Bruker and 3.0 T Siemens systems, respectively. The 3D MGRE and diffusion tensor imaging (DTI) data were acquired. Our results showed a highly significant correlation between MGRE-derived MWF and histological stain of myelin, and susceptibility and diffusivity also demonstrated a significant association with myelin. Both MWF and R2* (R2* = 1/T2*) values changed as a function of orientation, and the function varied with age. Furthermore, MWF and R2* were more sensitive to age than DTI. In vivo 3 T human MWF also changed substantially with the orientation as well. Our results support that MGRE-derived MWF can be used to assess the myelin content quantitatively.

Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal

Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.

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.

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).

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo

PURPOSE

Recent studies indicate that T1 in white matter (WM) is influenced by fiber orientation in B0 . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T1 relaxation time in humans in vivo as well as in rat brain ex vivo.

METHODS

Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T1 plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T1 and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T1 within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B0 and T1 , and diffusion MRI images acquired at 9.4 T. T1 angular plots were determined at several rotation angles in B0 .

RESULTS

Angular T1 plots from global WM provided reference for estimated fiber orientation-linked T1 changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T1 , matching that estimated from WM T1 data. In CC, where large and giant axons are numerous, the measured T1 change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T1 plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B0 to the T1 relaxation anisotropy in WM.

Myelin biomarkers in the healthy adult brain: Correlation, reproducibility, and the effect of fiber orientation

PURPOSE

We investigated the correlation, reproducibility, and effect of white matter fiber orientation for three myelin-sensitive MRI techniques: magnetization transfer ratio (MTR), inhomogeneous magnetization transfer ratio (ihMTR), and gradient and spin echo-derived myelin water fraction (MWF).

METHODS

We measured the three metrics in 17 white and three deep grey matter regions in 17 healthy adults at 3 T.

RESULTS

We found a strong correlation between ihMTR and MTR (r = 0.70, p < 0.001) and ihMTR and MWF (r = 0.79, p < 0.001), and a weaker correlation between MTR and MWF (r = 0.54, p < 0.001). The dynamic range in white matter was greatest for MWF (2.0%-27.5%), followed by MTR (14.4%-23.2%) and then ihMTR (1.2%-5.4%). The average scan-rescan coefficient of variation for white matter regions was 0.6% MTR, 0.3% ihMTR, and 0.7% MWF in metric units; however, when adjusted by the dynamic range, these became 6.3%, 6.1% and 2.8%, respectively. All three metrics varied with fiber direction: MWF and ihMTR were lower in white matter fibers perpendicular to B₀ by 6% and 1%, respectively, compared with those parallel, whereas MTR was lower by 0.5% at about 40°, with the highest values at 90°. However, separating the apparent orientation dependence by white matter region revealed large dissimilarities in the trends, suggesting that real differences in myelination between regions are confounding the apparent orientation dependence measured using this method.

CONCLUSION

The strong correlation between ihMTR and MWF suggests that these techniques are measuring the same myelination; however, the larger dynamic range of MWF may provide more power to detect small differences in myelin.

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.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

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.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

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.

Introduction of a new simple dynamic phantom for physical BOLD effect simulation

OBJECTIVE

To propose a new method of simulating the BOLD contrast using a dynamic, easy to construct and operate, low-cost physical phantom.

MATERIALS AND METHODS

A structure of thin pipelines passing through a gel volume was used to simulate blood vessels in human tissue. Quantitative T2*, R2* measurements were used to study the signal change of the phantom. BOLD fMRI experiments and analysis were performed to evaluate its potential use as an fMRI simulator.

RESULTS

Experimental T2*, R2* measurements showed similar behavior with published references. BOLD contrast was successfully achieved with the proposed method. In addition, there were several proposed parameters, like the angle of the phantom relative to B0, which can easily adjust the signal change and the activation area. Coefficients of variation showed good reproducibility within a month period. Statistical t maps were produced with in-house software for the BOLD measurements.

DISCUSSION

T2*maps and BOLD images confirm the potential use of this phantom as an fMRI simulator and also as a tool for studying sensitivity and specificity of BOLD sequences/algorithms.

Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data

PURPOSE

A method named DECOMPOSE-QSM is developed to decompose bulk susceptibility measured with QSM into sub-voxel paramagnetic and diamagnetic components based on a three-pool complex signal model.

METHODS

Multi-echo gradient echo signal is modeled as a summation of three weighted exponentials corresponding to three types of susceptibility sources: reference susceptibility, diamagnetic and paramagnetic susceptibility relative to the reference. Paramagnetic component susceptibility (PCS) and diamagnetic component susceptibility (DCS) maps are constructed to represent the sub-voxel compartments by solving for linear and nonlinear parameters in the model.

RESULTS

Numerical forward simulation and phantom validation confirmed the ability of DECOMPOSE-QSM to separate the mixture of paramagnetic and diamagnetic components. The PCS obtained from temperature-variant brainstem imaging follows the Curie's Law, which further validated the model and the solver. Initial in vivo investigation of human brain images showed the ability to extract sub-voxel PCS and DCS sources that produce visually enhanced contrast between brain structures comparing to threshold QSM.

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.

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.

Myelin quantification with MRI: A systematic review of accuracy and reproducibility

OBJECTIVES

Currently, multiple sclerosis is treated with anti-inflammatory therapies, but these treatments lack efficacy in progressive disease. New treatment strategies aim to repair myelin damage and efficacy evaluation of such new therapies would benefit from validated myelin imaging techniques. Several MRI methods for quantification of myelin density are available now. This systematic review aims to analyse the performance of these MRI methods.

METHODS

Studies comparing myelin quantification by MRI with histology, the current gold standard, or assessing reproducibility were retrieved from PubMed/MEDLINE and Embase (until December 2019). Included studies assessed both myelin histology and MRI quantitatively. Correlation or variance measurements were extracted from the studies. Non-parametric tests were used to analyse differences in study methodologies.

RESULTS

The search yielded 1348 unique articles. Twenty-two animal studies and 13 human studies correlated myelin MRI with histology. Eighteen clinical studies analysed the reproducibility. Overall bias risk was low or unclear. All MRI methods performed comparably, with a mean correlation between MRI and histology of R²=0.54 (SD=0.30) for animal studies, and R²=0.54 (SD=0.18) for human studies. Reproducibility for the MRI methods was good (ICC=0.75-0.93, R²=0.90-0.98, COV=1.3-27%), except for MTR (ICC=0.05-0.51).

CONCLUSIONS

Overall, MRI-based myelin imaging methods show a fairly good correlation with histology and a good reproducibility. However, the amount of validation data is too limited and the variability in performance between studies is too large to select the optimal MRI method for myelin quantification yet.

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.

Superficial white matter imaging: Contrast mechanisms and whole-brain in vivo mapping

Superficial white matter (SWM) contains the most cortico-cortical white matter connections in the human brain encompassing the short U-shaped association fibers. Despite its importance for brain connectivity, very little is known about SWM in humans, mainly due to the lack of noninvasive imaging methods. Here, we lay the groundwork for systematic in vivo SWM mapping using ultrahigh resolution 7 T magnetic resonance imaging. Using biophysical modeling informed by quantitative ion beam microscopy on postmortem brain tissue, we demonstrate that MR contrast in SWM is driven by iron and can be linked to the microscopic iron distribution. Higher SWM iron concentrations were observed in U-fiber-rich frontal, temporal, and parietal areas, potentially reflecting high fiber density or late myelination in these areas. Our SWM mapping approach provides the foundation for systematic studies of interindividual differences, plasticity, and pathologies of this crucial structure for cortico-cortical connectivity in humans.

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.

The influence of iron oxidation state on quantitative MRI parameters in post mortem human brain

A variety of Magnetic Resonance Imaging (MRI) techniques are known to be sensitive to brain iron content. In principle, iron sensitive MRI techniques are based on local magnetic field variations caused by iron particles in tissue. The purpose of this study was to investigate the sensitivity of MR relaxation and magnetization transfer parameters to changes in iron oxidation state compared to changes in iron concentration. Therefore, quantitative MRI parameters including R₁, R₂, R₂∗, quantitative susceptibility maps (QSM) and magnetization transfer ratio (MTR) of post mortem human brain tissue were acquired prior and after chemical iron reduction to change the iron oxidation state and chemical iron extraction to decrease the total iron concentration. All assessed parameters were shown to be sensitive to changes in iron concentration whereas only R₂, R₂∗ and QSM were also sensitive to changes in iron oxidation state. Mass spectrometry confirmed that iron accumulated in the extraction solution but not in the reduction solution. R₂∗ and QSM are often used as markers for iron content. Changes in these parameters do not necessarily reflect variations in iron content but may also be a result of changes in the iron's oxygenation state from ferric towards more ferrous iron or vice versa.

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.

Dependence of resting-state fMRI fluctuation amplitudes on cerebral cortical orientation relative to the direction of B0 and anatomical axes

Functional magnetic resonance imaging (fMRI) is now capable of sub-millimetre scale measurements over the entire human brain, however with such high resolutions each voxel is influenced by the local fine-scale details of the cerebral cortical vascular anatomy. The cortical vasculature is structured with the pial vessels lying tangentially along the grey matter surface, intracortical diving arterioles and ascending venules running perpendicularly to the surface, and a randomly oriented capillary network within the parenchyma. It is well-known that the amplitude of the blood-oxygenation level dependent (BOLD) signal emanating from a vessel depends on its orientation relative to the B0-field. Thus the vascular geometric hierarchy will impart an orientation dependence to the BOLD signal amplitudes and amplitude differences due to orientation differences constitute a bias for interpreting neuronal activity. Here, we demonstrate a clear effect of cortical orientation to B0 in the resting-state BOLD-fMRI amplitude (quantified as the coefficient of temporal signal variation) for 1.1 mm isotropic data at 7T and 2 mm isotropic at 3T. The maximum bias, i.e. the fluctuation amplitude difference between regions where cortex is perpendicular to vs. parallel to B0, is about +70% at the pial surface at 7T and +11% at 3T. The B0 orientation bias declines with cortical depth, becomes progressively smaller closer to the white matter surface, but then increases again to a local maximum within the white matter just beneath the cortical grey matter, suggesting a distinct tangential network of white matter vessels that also generate a BOLD orientation effect. We further found significant (negative) biases with the cortex orientation to the anterior-posterior anatomical axis of the head: a maximum negative bias of about -30% at the pial surface at 7T and about -13% at 3T. The amount of signal variance explained by the low frequency drift, motion and the respiratory cycle also showed a cortical orientation dependence; only the cardiac cycle induced signal variance was independent of cortical orientation, suggesting that the cardiac induced component of the image time-series fluctuations is not related to a significant change in susceptibility. Although these orientation effects represent a signal bias, and are likely to be a nuisance in high-resolution analyses, they may help characterize the vascular influences on candidate fMRI acquisitions and, thereby, may be exploited to improve the neuronal specificity of fMRI.

Biophysically motivated efficient estimation of the spatially isotropic R 2 * component from a single gradient-recalled echo measurement

Purpose

To propose and validate an efficient method, based on a biophysically motivated signal model, for removing the orientation‐dependent part of R2* using a single gradient‐recalled echo (GRE) measurement.

Methods

The proposed method utilized a temporal second‐order approximation of the hollow‐cylinder‐fiber model, in which the parameter describing the linear signal decay corresponded to the orientation‐independent part of R2*. The estimated parameters were compared to the classical, mono‐exponential decay model for R2* in a sample of an ex vivo human optic chiasm (OC). The OC was measured at 16 distinct orientations relative to the external magnetic field using GRE at 7T. To show that the proposed signal model can remove the orientation dependence of R2*, it was compared to the established phenomenological method for separating R2* into orientation‐dependent and ‐independent parts.

Results

Using the phenomenological method on the classical signal model, the well‐known separation of R2* into orientation‐dependent and ‐independent parts was verified. For the proposed model, no significant orientation dependence in the linear signal decay parameter was observed.

Conclusions

Since the proposed second‐order model features orientation‐dependent and ‐independent components at distinct temporal orders, it can be used to remove the orientation dependence of R2* using only a single GRE measurement.

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.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

The influence of brain iron on myelin water imaging

With myelin playing a vital role in normal brain integrity and function and thus in various neurological disorders, myelin sensitive magnetic resonance imaging (MRI) techniques are of great importance. In particular, multi-exponential T₂ relaxation was shown to be highly sensitive to myelin. The myelin water imaging (MWI) technique allows to separate the T₂ decay into short components, specific to myelin water, and long components reflecting the intra- and extracellular water. The myelin water fraction (MWF) is the ratio of the short components to all components. In the brain's white matter (WM), myelin and iron are closely linked via the presence of iron in the myelin generating oligodendrocytes. Iron is known to decrease T₂ relaxation times and may therefore mimic myelin. In this study, we investigated if variations in WM iron content can lead to apparent MWF changes. We performed MWI in post mortem human brain tissue prior and after chemical iron extraction. Histology for iron and myelin confirmed a decrease in iron content and no change in myelin content after iron extraction. In MRI, iron extraction lead to a decrease in MWF by 26%-28% in WM. Thus, a change in MWF does not necessarily reflect a change in myelin content. This observation has important implications for the interpretation of MWI findings in previously published studies and future research.

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.

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.

Untangling the R2* contrast in multiple sclerosis: A combined MRI-histology study at 7.0 Tesla

T2*-weighted multi-echo gradient-echo magnetic resonance imaging and its reciprocal R2* are used in brain imaging due to their sensitivity to iron content. In patients with multiple sclerosis who display pathological alterations in iron and myelin contents, the use of R2* may offer a unique way to untangle mechanisms of disease. Coronal slices from 8 brains of deceased multiple sclerosis patients were imaged using a whole-body 7.0 Tesla MRI scanner. The scanning protocol included three-dimensional (3D) T2*-w multi-echo gradient-echo and 2D T2-w turbo spin echo (TSE) sequences. Histopathological analyses of myelin and iron content were done using Luxol fast blue and proteolipid myelin staining and 3,3′-diaminobenzidine tetrahydrochloride enhanced Turnbull blue staining. Quantification of R2*, myelin and iron intensity were obtained. Variations in R2* were found to be affected differently by myelin and iron content in different regions of multiple sclerosis brains. The data shall inform clinical investigators in addressing the role of T2*/R2* variations as a biomarker of tissue integrity in brains of MS patients, in vivo.

Assessment of ferritin content in multiple sclerosis brains using temperature-induced R*2 changes

PURPOSE

Current MRI techniques cannot reliably assess iron content in white matter due to the confounding diamagnetic effect of myelin. The purpose of this study was to validate with histology a novel iron mapping technique that uses the temperature dependency of the paramagnetic susceptibility in multiple sclerosis (MS) brains, where white matter has been reported to show significant variations in iron content.

METHODS

We investigated post mortem brain tissue from three MS patients and one control subject. Temperature-dependent R2* relaxometry was performed between 4°C and 37°C. The resulting temperature coefficient ( TcR2*) maps were compared with immunohistochemical stains for ferritin light chain.

RESULTS

Good agreement between TcR2* maps and ferritin staining was found by way of visual comparison and quantitative analysis. The highest iron concentrations were detected at the edge of MS lesions and in the basal ganglia. For all regions, except the subcortical U-fibers, there was a significant negative correlation between the TcR2* values and the ferritin count.

CONCLUSION

This study provides further evidence that TcR2* may be a reliable measure of white matter iron content due to the elimination of myelin-induced susceptibility changes and is well suited for further research into neurological diseases with distortions of the iron homeostasis. Magn Reson Med 79:1609-1615, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Internal gradient distributions: A susceptibility-derived tensor delivering morphologies by magnetic resonance

Nuclear magnetic resonance is a powerful tool for probing the structures of chemical and biological systems. Combined with field gradients it leads to NMR imaging (MRI), a widespread tool in non-invasive examinations. Sensitivity usually limits MRI’s spatial resolution to tens of micrometers, but other sources of information like those delivered by constrained diffusion processes, enable one extract morphological information down to micron and sub-micron scales. We report here on a new method that also exploits diffusion – isotropic or anisotropic– to sense morphological parameters in the nm-mm range, based on distributions of susceptibility-induced magnetic field gradients. A theoretical framework is developed to define this source of information, leading to the proposition of internal gradient-distribution tensors. Gradient-based spin-echo sequences are designed to measure these new observables. These methods can be used to map orientations even when dealing with unconstrained diffusion, as is here demonstrated with studies of structured systems, including tissues.

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.

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.