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

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

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

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.

Making acute ischemic stroke thrombi visible in MRI imaging

Knowledge of thrombus behavior and visualization on MRI in acute ischemic stroke is less than optimal. However, MRI sequences could be enhanced based on the typical T1 and T2 relaxation times of the target tissues, which mainly determine their signal intensities on imaging. We studied the relaxation times of a broad spectrum of clot analogs along with their image characteristics of three sequences analyzed: a T1-weighted turbo inversion-recovery sequence (T1w Turbo IR), a T1-weighted turbo spin echo with fat suppression (T1w TSE SPIR), and a T2-weighted 3D TSE with magnetization refocusing to remove T1 dependence (T2w TSE DRIVE). We compared their imaging behavior with the intensity values of normal brain tissue using the same imaging protocols as for clots. Each histological and biochemical clot component contributed to each of the relaxation times. Overall, histological composition correlated strongly with T1 times, and iron content, specifically, with T2 relaxation time. Using decision trees, fibrin content was selected as the primary biomarker for T1 relaxation times, inducing an increase. Up to four clot subgroups could be defined based on its distinctive T1 relaxation time. Clot signal intensity in the T1 and T2-weighted images varied significantly according to T1 and T2 relaxation times. Moreover, in comparison with normal brain tissue intensity values, T2w DRIVE images depict thrombi according to the principle of the more fibrin, the higher the intensity, and in T1w TSE, the more erythrocytes, the higher the intensity. These findings could facilitate improvements in MRI sequences for clot visualization and indicate that T2w DRIVE and T1w TSE sequences should depict the vast majority of acute ischemic stroke thrombi as more hyperintense than surrounding tissues.

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.

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 T₁ in white matter (WM) is influenced by fiber orientation in B₀ . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T₁ 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 T₁ plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T₁ and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T₁ within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B₀ and T₁ , and diffusion MRI images acquired at 9.4 T. T₁ angular plots were determined at several rotation angles in B₀ .

RESULTS

Angular T₁ plots from global WM provided reference for estimated fiber orientation-linked T₁ 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 T₁ , matching that estimated from WM T₁ data. In CC, where large and giant axons are numerous, the measured T₁ change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T₁ plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B₀ to the T₁ relaxation anisotropy in WM.

White matter microstructure and longitudinal relaxation time anisotropy in human brain at 3 and 7 T

A high degree of structural order by white matter (WM) fibre tracts creates a physicochemical environment where water relaxations are rendered anisotropic. Recently, angularly dependent longitudinal relaxation has been reported in human WM. We have characterised interrelationships between T1 relaxation and diffusion MRI microstructural indices at 3 and 7 T. Eleven volunteers consented to participate in the study. Multishell diffusion MR images were acquired with b-values of 0/1500/3000 and 0/1000/2000 s/mm² at 1.5 and 1.05 mm³ isotropic resolutions at 3 and 7 T, respectively. DTIFIT was used to compute DTI indices; the fibre-to-field angle (θ_FB ) maps were obtained using the principal eigenvector images. The orientations and volume fractions of multiple fibre populations were estimated using BedpostX in FSL, and the orientation dispersion index (ODI) was estimated using the NODDI protocol. MP2RAGE was used to acquire images for T1 maps at 1.0 and 0.9 mm³ isotropic resolutions at 3 and 7 T, respectively. At 3 T, T1 as a function of θ_FB in WM with high fractional anisotropy and one-fibre orientation volume fraction or low ODI shows a broad peak centred at 50^o , but a flat baseline at 0^o and 90^o . The broad peak amounted up to 7% of the mean T1. At 7 T, the broad peak appeared at 40^o and T1 in fibres running parallel to B0 was longer by up to 75 ms (8.3% of the mean T1) than in those perpendicular to the field. The peak at 40^o was approximately 5% of mean T1 (i.e., proportionally smaller than that at 54^o at 3 T). The data demonstrate T1 anisotropy in WM with high microstructural order at both fields. The angular patterns are indicative of the B0-dependency of T1 anisotropy. Thus myelinated WM fibres influence T1 contrast both by acting as a T1 contrast agent and rendering T1 dependent on fibre orientation with B0.

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.

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.

Probing the myelin water compartment with a saturation-recovery, multi-echo gradient-recalled echo sequence

PURPOSE

To investigate the effect of varying levels of T 1 -weighting on the evolution of the complex signal from white matter in a multi-echo gradient-recalled echo (mGRE) saturation-recovery sequence.

THEORY AND METHODS

Analysis of the complex signal evolution in an mGRE sequence allows the contributions from short- and long- T 2 ∗ components to be separated, thus providing a measure of the relative strength of signals from the myelin water, and the external and intra-axonal compartments. Here we evaluated the effect of different levels of T 1 -weighting on these signals, expecting that the previously reported, short T 1 of the myelin water would lead to a relative enhancement of the myelin water signal in the presence of signal saturation. Complex, saturation-recovery mGRE data from the splenium of the corpus callosum from 5 healthy volunteers were preprocessed using a frequency difference mapping (FDM) approach and analyzed using the 3-pool model of complex signal evolution in white matter.

RESULTS

An increase in the apparent T 1 as a function of echo time was demonstrated, but this increase was an order of magnitude smaller than that expected from previously reported myelin water T 1 -values. This suggests the presence of magnetization transfer and exchange effects which counteract the T 1 -weighting.

CONCLUSION

Variation of the B 1 + amplitude in a saturation-recovery mGRE sequence can be used to modulate the relative strength of signals from the different compartments in white matter, but the modulation is less than predicted from previously reported T 1 -values.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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