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

An unconstrained four pool model analysis of proton relaxation and magnetization transfer in ex vivo white matter

Understanding proton relaxation in the brain's white matter remains an active field of magnetic resonance imaging research. Models of varying complexity have been proposed to link measurements to tissue composition/microstructure, in particular myelination. Although the presence of multiple aqueous and nonaqueous proton pools is well established experimentally, so-called "quantitative MRI" is usually based on simpler models due to the large number of model parameters. In this work, a comprehensive set of parameters characterizing a four-pool model is obtained. A piece of fixed porcine spinal-cord WM was investigated at 3 T and temperatures between 21 and 35 °C. Measurements included a wide range of preparations of the spin system in combination with long echo trains to achieve sensitivity to all model parameters. The results allow the extraction of all intrinsic relaxation and exchange rates as well as assigning them to specific dynamic processes involving tissue water. A critical assessment indicates that simpler models often lack specificity to myelin.

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.

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.

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.

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.

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.

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.