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

Deciphering adiabatic rotating frame relaxometry in biological tissues

PURPOSE

This work aims to unravel the intricacies of adiabatic rotating frame relaxometry in biological tissues.

THEORY AND METHODS

The classical formalisms of dipolar relaxationR 1 ρ $$ {R}_{1\rho } $$ andR 2 ρ $$ {R}_{2\rho } $$ were systematically analyzed for water molecules reorienting on "fast" and "slow" timescales. These two timescales are, respectively, responsible for the absence and presence ofR 1 ρ $$ {R}_{1\rho } $$ dispersion. A time-averagedR 1 ρ $$ {R}_{1\rho } $$ orR 2 ρ $$ {R}_{2\rho } $$ over an adiabatic pulse duration was recast into a sum ofR 1 $$ {R}_1 $$ andR 2 $$ {R}_2 $$ , but with different weightings. These weightings depend on the specific modulations of adiabatic pulse waveforms. In this context, stretched hyperbolic secant (HSn $$ HSn $$ ) pulses were characterized. Previously publishedH S 1 $$ HS1 $$ R 1 ρ $$ {R}_{1\rho } $$ , continuous-wave (CW)R 1 ρ $$ {R}_{1\rho } $$ , andR 1 $$ {R}_1 $$ measures from 12 agarose phantoms were used to validate the theoretical predictions. A similar validation was also performed on previously publishedHSn $$ HSn $$ R 1 ρ $$ {R}_{1\rho } $$ (n $$ n $$ =1, 4, 8) andHS 1 $$ HS1 $$ R 2 ρ $$ {R}_{2\rho } $$ from bovine cartilage specimens.

RESULTS

Longitudinal relaxation weighting decreases forHSn $$ HSn $$ pulses asn $$ n $$ increases. Predicted CWR 1 ρ cal $$ {R}_{1\rho}^{cal} $$ values from agarose phantoms align well with the measured CWR 1 ρ exp $$ {R}_{1\rho}^{exp} $$ values, as indicated by a linear regression function:R 1 ρ cal = 1.04 * R 1 ρ exp - 1.96 $$ {R}_{1\rho}^{cal}={1.04}^{\ast }{R}_{1\rho}^{exp}-1.96 $$ . The predicted adiabaticR 1 ρ $$ {R}_{1\rho } $$ andR 2 ρ $$ {R}_{2\rho } $$ from cartilage specimens are consistent with those previously measured, as quantified by:R 1 ρ , 2 ρ cal = 1.10 * R 1 ρ , 2 ρ exp - 0.41 $$ {R}_{1\rho, 2\rho}^{cal}={1.10}^{\ast }{R}_{1\rho, 2\rho}^{exp}-0.41 $$ .

CONCLUSION

This work has theoretically and experimentally demonstrated that adiabaticR 1 ρ $$ {R}_{1\rho } $$ andR 2 ρ $$ {R}_{2\rho } $$ can be recast into a sum ofR 1 $$ {R}_1 $$ andR 2 $$ {R}_2 $$ , with varying weightings. Therefore, any suggestions that adiabatic rotating frame relaxometry in biological tissues could provide more information than the standardR 1 $$ {R}_1 $$ andR 2 $$ {R}_2 $$ warrant closer scrutiny.

Characterization of the orientation dependence of magnetization transfer measures in single and crossing-fiber white matter

PURPOSE

To fully characterize the orientation dependence of magnetization transfer (MT) and inhomogeneous MT (ihMT) measures in the whole white matter (WM), for both single-fiber and crossing-fiber voxels.

METHODS

A characterization method was developed using the fiber orientation obtained from diffusion MRI (dMRI) with diffusion tensor imaging (DTI) and constrained spherical deconvolution. This allowed for characterization of the orientation dependence of measures in all of WM, regardless of the number of fiber orientation in a voxel. Furthermore, the orientation dependence inside 31 different WM bundles was characterized to evaluate the homogeneity of the effect. Variation of the results within and between-subject was assessed from a 12-subject dataset.

RESULTS

Previous results for single-fiber voxels were reproduced and a novel characterization was produced in voxels of crossing fibers, which seems to follow trends consistent with single-fiber results. Heterogeneity of the orientation dependence across bundles was observed, but homogeneity within similar bundles was also highlighted. Differences in behavior between MT and ihMT measures, as well as the ratio and saturation versions of these, were noted.

CONCLUSION

Orientation dependence characterization was proven possible over the entirety of WM. The vast range of effects and subtleties of the orientation dependence on MT measures showed the need for, but also the challenges of, a correction method.

T1 relaxation and axon fibre configuration in human white matter

Understanding the effects of white matter (WM) axon fibre microstructure on T1 relaxation is important for neuroimaging. Here, we have studied the interrelationship between T1 and axon fibre configurations at 3T and 7T. T1 and S0 (=signal intensity at zero TI) were computed from MP2RAGE images acquired with six inversion recovery times. Multishell diffusion MRI images were analysed for fractional anisotropy (FA); MD; V1; the volume fractions for the first (f1), second (f2) and third (f3) fibre configuration; and fibre density cross-section images for the first (fdc1), second (fdc2) and third (fdc3) fibres. T1 values were plotted as a function of FA, f1, f2, f3, fdc1, fdc2 and fdc3 to examine interrelationships between the longitudinal relaxation and the diffusion MRI microstructural measures. T1 values decreased with increasing FA, f1 and f2 in a nonlinear fashion. At low FA values (from 0.2 to 0.4), a steep shortening of T1 was followed by a shallow shortening by 6%-10% at both fields. The steep shortening was associated with decreasing S0 and MD. T1 also decreased with increasing fdc1 values in a nonlinear fashion. Instead, only a small T1 change as a function of either f3 or fdc3 was observed. In WM areas selected by fdc1 only masks, T1 was shorter than in those with fdc2/fdc3. In WM areas with high single fibre populations, as delineated by f1/fdc1 masks, T1 was shorter than in tissue with high complex fibre configurations, as segmented by f2/fdc2 or f3/fdc3 masks. T1 differences between these WM areas are attributable to combined effects by T1 anisotropy and lowered FA. The current data show strong interrelationships between T1, axon fibre configuration and orientation in healthy WM. It is concluded that diffusion MRI microstructural measures are essential in the effort to interpret quantitative T1 images in terms of tissue state in health and disease.

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.

A straightforward procedure to build a non-toxic relaxometry phantom with desired T1 and T2 times at 3T

OBJECTIVE

To prepare and analyze soy-lecithin-agar gels for non-toxic relaxometry phantoms with tissue-like relaxation times at 3T.

METHODS

Phantoms mimicking the relaxation times of various tissues (gray and white matter, kidney cortex and medulla, spleen, muscle, liver) were built and tested with a clinical 3T whole-body MR scanner. Simple equations were derived to calculate the appropriate concentrations of soy lecithin and agar in aqueous solutions to achieve the desired relaxation times. Phantoms were tested for correspondence between measurements and calculated T1 and T2 values, reproducibility, spatial homogeneity, and temporal stability. T1 and T2 mapping techniques and a 3D T1-weighted sequence with high spatial resolution were applied.

RESULTS

Except for the liver relaxation phantom, all phantoms were successfully and reproducibly produced. Good agreement was found between the targeted and measured relaxation times. The percentage deviations from the targeted relaxation times were less than 3% for T1 and less than 6.5% for T2. In addition, the phantoms were homogeneous and had little to no air bubbles. However, the phantoms were unstable over time: after a storage period of 4 weeks, mold growth and also changes in relaxation times were detected in almost all phantoms.

CONCLUSION

Soy-lecithin-agar gels are a non-toxic material for the construction of relaxometry phantoms with tissue-like relaxation times. They are easy to prepare, inexpensive and allow independent adjustment of T1 and T2. However, there is still work to be done to improve the long-term stability of the phantoms.

B1 inhomogeneity-corrected T1 mapping and quantitative magnetization transfer imaging via simultaneously estimating Bloch-Siegert shift and magnetization transfer effects

PURPOSE

To introduce a method of inducing Bloch-Siegert shift and magnetization Transfer Simultaneously (BTS) and demonstrate its utilization for measuring binary spin-bath model parameters free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular fraction (f $$ f $$ ), magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) and local transmit field (B 1 + $$ {B}_1^{+} $$ ).

THEORY AND METHODS

Bloch-Siegert shift and magnetization transfer is simultaneously induced through the application of off-resonance irradiation in between excitation and acquisition of an RF-spoiled gradient-echo scheme. Applying the binary spin-bath model, an analytical signal equation is derived and verified through Bloch simulations. Monte Carlo simulations were performed to analyze the method's performance. The estimation of the binary spin-bath parameters withB 1 + $$ {B}_1^{+} $$ compensation was further investigated through experiments, both ex vivo and in vivo.

RESULTS

Comparing BTS with existing methods, simulations showed that existing methods can significantly biasT 1 $$ {T}_1 $$ estimation when not accounting for transmitB 1 $$ {B}_1 $$ heterogeneity and MT effects that are present. Phantom experiments further showed that the degree of this bias increases with increasing macromolecular proton fraction. Multi-parameter fit results from an in vivo brain study generated values in agreement with previous literature. Based on these studies, we confirmed that BTS is a robust method for estimating the binary spin-bath parameters in macromolecule-rich environments, even in the presence ofB 1 + $$ {B}_1^{+} $$ inhomogeneity.

CONCLUSION

A method of estimating Bloch-Siegert shift and magnetization transfer effect has been developed and validated. Both simulations and experiments confirmed that BTS can estimate spin-bath parameters (T 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ) that are free fromB 1 + $$ {B}_1^{+} $$ bias.

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

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.

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.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

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.

Reliability of quantitative multiparameter maps is high for magnetization transfer and proton density but attenuated for R1 and R2 * in healthy young adults

We investigate the reliability of individual differences of four quantities measured by magnetic resonance imaging‐based multiparameter mapping (MPM): magnetization transfer saturation (MT), proton density (PD), longitudinal relaxation rate (R₁), and effective transverse relaxation rate (R₂*). Four MPM datasets, two on each of two consecutive days, were acquired in healthy young adults. On Day 1, no repositioning occurred and on Day 2, participants were repositioned between MPM datasets. Using intraclass correlation effect decomposition (ICED), we assessed the contributions of session‐specific, day‐specific, and residual sources of measurement error. For whole‐brain gray and white matter, all four MPM parameters showed high reproducibility and high reliability, as indexed by the coefficient of variation (CoV) and the intraclass correlation (ICC). However, MT, PD, R₁, and R₂* differed markedly in the extent to which reliability varied across brain regions. MT and PD showed high reliability in almost all regions. In contrast, R₁ and R₂* showed low reliability in some regions outside the basal ganglia, such that the sum of the measurement error estimates in our structural equation model was higher than estimates of between‐person differences. In addition, in this sample of healthy young adults, the four MPM parameters showed very little variability over four measurements but differed in how well they could assess between‐person differences. We conclude that R₁ and R₂* might carry only limited person‐specific information in some regions of the brain in healthy young adults, and, by implication, might be of restricted utility for studying associations to between‐person differences in behavior in those regions.

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.

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.