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

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.

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.

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.

Radiation damping at clinical field strength: Characterization and compensation in quantitative measurements

PURPOSE

In any MR experiment, the bulk magnetization acts on itself, caused by the induced current in the RF receiver circuit that generates an oscillating damping field. This effect, known as "radiation damping" (RD), is usually weak and, therefore, unconsidered in MRI, but can affect quantitative studies performed with dedicated coils that provide a high SNR. The current work examined RD in a setup for investigations of small tissue specimens including a quantitative characterization of the spin-coil system.

THEORY AND METHODS

A custom-made Helmholtz coil (radius and spacing 16 mm) was interfaced to a transmit-receive (Tx/Rx) switch with integrated passive feedback for modulation or suppression of RD similar to preamplifier decoupling. Pulse sequences included pulse-width arrays to demonstrate the absence/ presence of RD and difference techniques employing gradient pulses or composite RF pulses to quantify RD effects during free precession and transmission, respectively. Experiments were performed at 3T in small samples of MnCl2 solution.

RESULTS

Significant RD effects may impact RF pulse application and evolution periods. Effective damping time constants were comparable to typical T2 * times or echo spacings in multi-echo sequences. Measurements of the phase relation showed that deviations from the commonly assumed 90° angle between the damping field and the transverse magnetization may occur.

CONCLUSION

Radiation damping may affect the accuracy of quantitative MR measurements performed with dedicated RF coils. Efficient mitigation can be achieved hardware-based or by appropriate consideration in the pulse sequence.