Understanding aqueous and non-aqueous proton T1 relaxation in brain

Validation of an ultrahigh contrast divided subtracted inversion recovery technique using a standard T1 phantom

The divided subtracted inversion recovery (dSIR) is a high T1 contrast technique that shows changes in white matter in patients with traumatic brain injury and hypoxic injury. The changes can be explained by small differences in T1; however, to date, there has been no independent validation of the technique using a standard reference. The present study develops the theory of the dSIR signal and performs validation using the NIST/ISMRM T1 phantom. Non-idealities are explored, including the influence of noise bias and finite repetition time (TR), which leads to the introduction of an optimally efficient TR for inversion recovery acquisitions. Results show excellent agreement with theoretical calculations.

Multicompartment imaging of the brain using a comprehensive MR imaging protocol

In this study, we describe a comprehensive 3D magnetic resonance imaging (MRI) protocol designed to assess major tissue and fluid components in the brain. The protocol comprises four different sequences: 1) magnetization transfer prepared Cones (MT-Cones) for two-pool MT modeling to quantify macromolecular content; 2) short-TR adiabatic inversion-recovery prepared Cones (STAIR-Cones) for myelin water imaging; 3) proton-density weighted Cones (PDw-Cones) for total water imaging; and 4) highly T2 weighted Cones (T2w-Cones) for free water imaging. By integrating these techniques, we successfully mapped key brain components-namely macromolecules, myelin water, intra/extracellular water, and free water-in ten healthy volunteers and five patients with multiple sclerosis (MS) using a 3T clinical scanner. Brain macromolecular proton fraction (MMPF), myelin water proton fraction (MWPF), intra/extracellular water proton fraction (IEWPF), and free water proton fraction (FWPF) values were generated in white matter (WM), grey matter (GM), and MS lesions. Excellent repeatability of the protocol was demonstrated with high intra-class correlation coefficient (ICC) values. In MS patients, the MMPF and MWPF values of the lesions and normal-appearing WM (NAWM) were significantly lower than those in normal WM (NWM) in healthy volunteers. Moreover, we observed significantly higher FWPF values in MS lesions compared to those in NWM and NAWM regions. This study demonstrates the capability of our technique to volumetrically map major brain components. The technique may have particular value in providing a comprehensive assessment of neuroinflammatory and neurodegenerative diseases of the brain.

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.

Yet more evidence that non-aqueous myelin lipids can be directly imaged with ultrashort echo time (UTE) MRI on a clinical 3T scanner: a lyophilized red blood cell membrane lipid study

Direct imaging of semi-solid lipids, such as myelin, is of great interest as a noninvasive biomarker of neurodegenerative diseases. Yet, the short T2 relaxation times of semi-solid lipid protons hamper direct detection through conventional magnetic resonance imaging (MRI) pulse sequences. In this study, we examined whether a three-dimensional ultrashort echo time (3D UTE) sequence can directly acquire signals from membrane lipids. Membrane lipids from red blood cells (RBC) were collected from commercially available blood as a general model of the myelin lipid bilayer and subjected to D2O exchange and freeze-drying for complete water removal. Sufficiently high MR signals were detected with the 3D UTE sequence, which showed an ultrashort T2* of ∼77-271 µs and a short T1 of ∼189 ms for semi-solid RBC membrane lipids. These measurements can guide designing UTE-based sequences for direct in vivo imaging of membrane lipids.

Sensitivity of unconstrained quantitative magnetization transfer MRI to Amyloid burden in preclinical Alzheimer's disease

Magnetization transfer MRI is sensitive to semi-solid macromolecules, including amyloid beta, and has previously been used to discriminate Alzheimer's disease (AD) patients from controls. Here, we fit an unconstrained 2-pool quantitative MT (qMT) model, i.e., without constraints on the longitudinal relaxation rateR 1 s of semi-solids, and investigate the sensitivity of the estimated parameters to amyloid accumulation in preclinical subjects. We scanned 15 cognitively normal volunteers, of which 9 were amyloid positive by [18F]Florbetaben PET. A 12 min hybrid-state qMT scan with an effective resolution of 1.24 mm isotropic and whole-brain coverage was acquired to estimate the unconstrained 2-pool qMT parameters. Group comparisons and correlations with Florbetaben PET standardized uptake value ratios were analyzed at the lobar level. We find that the exchange rate and semi-solid pool'sR 1 s were sensitive to the amyloid concentration, while morphometric measures of cortical thickness derived from structural MRI were not. Changes in the exchange rate are consistent with previous reports in clinical AD, while changes inR 1 s have not been reported previously as its value is typically constrained in the literature. Our results demonstrate that qMT MRI may be a promising surrogate marker of amyloid beta without the need for contrast agents or radiotracers.

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.

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.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

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.

Optimization of acquisition parameters for cortical inhomogeneous magnetization transfer (ihMT) imaging using a rapid gradient echo readout

PURPOSE

Imaging biomarkers with increased myelin specificity are needed to better understand the complex progression of neurological disorders. Inhomogeneous magnetization transfer (ihMT) imaging is an emergent technique that has a high degree of specificity for myelin content but suffers from low signal to-noise ratio (SNR). This study used simulations to determine optimal sequence parameters for ihMT imaging for use in high-resolution cortical mapping.

METHODS

MT-weighted cortical image intensity and ihMT SNR were simulated using modified Bloch equations for a range of sequence parameters. The acquisition time was limited to 4.5 min/volume. A custom MT-weighted RAGE sequence with center-out k-space encoding was used to enhance SNR at 3 T. Pulsed MT imaging was studied over a range of saturation parameters, and the impact of the turbo factor on the effective ihMT resolution was investigated. 1 mm isotropic ihMTsat maps were generated in 25 healthy adults.

RESULTS

Greater SNR was observed for larger number of bursts consisting of 6-8 saturation pulses each, combined with a high readout turbo factor. However, that protocol suffered from a point spread function that was more than twice the nominal resolution. For high-resolution cortical imaging, we selected a protocol with a higher effective resolution at the cost of a lower SNR. We present the first group-average ihMTsat whole-brain map at 1 mm isotropic resolution.

CONCLUSION

This study presents the impact of saturation and excitation parameters on ihMTsat SNR and resolution. We demonstrate the feasibility of high-resolution cortical myelin imaging using ihMTsat in less than 20 min.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Quantitative magnetic resonance mapping of the myelin bilayer reflects pathology in multiple sclerosis brain tissue

Multiple sclerosis (MS) is a neuroinflammatory disease characterized by loss of myelin (demyelination) and, to a certain extent, subsequent myelin repair (remyelination). To better understand the pathomechanisms underlying de- and remyelination and to monitor the efficacy of treatments aimed at regenerating myelin, techniques offering noninvasive visualizations of myelin are warranted. Magnetic resonance (MR) imaging has long been at the forefront of efforts to visualize myelin, but it has only recently become feasible to access the rapidly decaying resonance signals stemming from the myelin lipid-protein bilayer itself. Here, we show that direct MR mapping of the bilayer yields highly specific myelin maps in brain tissue from patients with MS. Furthermore, examination of the bilayer signal behavior is found to reveal pathological alterations in normal-appearing white and gray matter. These results indicate promise for in vivo implementations of the myelin bilayer mapping technique, with prospective applications in basic research, diagnostics, disease monitoring, and drug development.

Inhomogeneous magnetization transfer imaging: Concepts and directions for further development

Off-resonance radio frequency irradiation can induce the ordering of proton spins in the dipolar fields of their neighbors, in molecules with restricted mobility. This dipolar order decays with a characteristic relaxation time, T_1D , that is very different from the T₁ and T₂ relaxation of the nuclear alignment with the main magnetic field. Inhomogeneous magnetization transfer (ihMT) imaging is a refinement of magnetization transfer (MT) imaging that isolates the MT signal dependence on dipolar order relaxation times within motion-constrained molecules. Because T_1D relaxation is a unique contrast mechanism, ihMT may enable improved characterization of tissue. Initial work has stressed the high correlation between ihMT signal and myelin density. Dipolar order relaxation appears to be much longer in membrane lipids than other molecules. Recent work has shown, however, that ihMT acquisitions may also be adjusted to emphasize different ranges of T_1D . These newer approaches may be sensitive to other microstructural components of tissue. Here, we review the concepts and history of ihMT and outline the requirements for further development to realize its full potential.

Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions

PURPOSE

To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases.

METHODS

The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T_1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T_2,f , T_2,b , R, and free pool T₁ ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T₁ .

RESULTS

Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T₁ (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT₁  = -19.0 ms). The single-point methodology yielded consistent MPF and T₁ values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter.

CONCLUSION

The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics

Inversion pulses are commonly employed in MRI for T 1 $$ {T}_1 $$ -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 $$ {B}_0 $$ and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 $$ {T}_1 $$ relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous ¹ H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 $$ {T}_1 $$ relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

High-resolution magnetization-transfer imaging of post-mortem marmoset brain: Comparisons with relaxometry and histology

Cell membranes and macromolecules or paramagnetic compounds interact with water proton spins, which modulates magnetic resonance imaging (MRI) contrast providing information on tissue composition. For a further investigation, quantitative magnetization transfer (qMT) parameters (at 3T), including the ratio of the macromolecular and water proton pools, F, and the exchange-rate constant as well as the (observed) longitudinal and the effective transverse relaxation rates (at 3T and 7T), R₁^obs and R₂^*, respectively, were measured at high spatial resolution (200 µm) in a slice of fixed marmoset brain and compared to histology results obtained with Gallyas' myelin stain and Perls' iron stain. R₁^obs and R₂^* were linearly correlated with the iron content for the entire slice, whereas distinct differences were obtained between gray and white matter for correlations of relaxometry and qMT parameters with myelin content. The combined results suggest that the macromolecular pool interacting with water consists of myelin and (less efficient) non-myelin contributions. Despite strong correlation of F and R₁^obs, none of these parameters was uniquely specific to myelination. Due to additional sensitivity to iron stores, R₁^obs and R₂^* were more sensitive for depicting microstructural differences between cortical layers than F.

Mapping the myelin bilayer with short-T2 MRI: Methods validation and reference data for healthy human brain

PURPOSE

To explore the properties of short-T₂ signals in human brain, investigate the impact of various experimental procedures on these properties and evaluate the performance of three-component analysis.

METHODS

Eight samples of non-pathological human brain tissue were subjected to different combinations of experimental procedures including D₂ O exchange and frozen storage. Short-T₂ imaging techniques were employed to acquire multi-TE (33-2067 μs) data, to which a three-component complex model was fitted in two steps to recover the properties of the underlying signal components and produce amplitude maps of each component. For validation of the component amplitude maps, the samples underwent immunohistochemical myelin staining.

RESULTS

The signal component representing the myelin bilayer exhibited super-exponential decay with T_2,min of 5.48 μs and a chemical shift of 1.07 ppm, and its amplitude could be successfully mapped in both white and gray matter in all samples. These myelin maps corresponded well to myelin-stained tissue sections. Gray matter signals exhibited somewhat different components than white matter signals, but both tissue types were well represented by the signal model. Frozen tissue storage did not alter the signal components but influenced component amplitudes. D₂ O exchange was necessary to characterize the non-aqueous signal components, but component amplitude mapping could be reliably performed also in the presence of H₂ O signals.

CONCLUSIONS

The myelin mapping approach explored here produced reasonable and stable results for all samples. The extensive tissue and methodological investigations performed in this work form a basis for signal interpretation in future studies both ex vivo and in vivo.

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.

Inhomogeneous Magnetization Transfer (ihMT) imaging in the acute cuprizone mouse model of demyelination/remyelination

BACKGROUND

To investigate the association of ihMT (inhom signals with the demyelination and remyelination phases of the acute cuprizone mouse model in comparison with histology, and to assess the extent of tissue damage and repair from MRI data.

METHODS

Acute demyelination by feeding 0.2% cuprizone for five weeks, followed by a four-week remyelination period was applied on genetically modified plp-GFP mice. Animals were scanned at different time points of the demyelination and remyelination phases of the cuprizone model using a multimodal MRI protocol, including ihMT T_1D-filters, MPF (Macromolecular Proton Fraction) and R₁ (longitudinal relaxation rate). For histology, plp-GFP (proteolipid protein - Green Fluorescent Protein) microscopy and LFB (Luxol Fast Blue) staining were employed as references for the myelin content. Comparison of MRI with histology was performed in the medial corpus callosum (mCC) and cerebral cortex (CTX) at two brain levels whereas ROI-wise and voxel-based analyses of the MRI metrics allowed investigating in vivo the spatial extent of myelin alterations.

RESULTS

IhMT high-pass T_1D-filters, targeted toward long T_1D components, showed significant temporal variations in the mCC consistent with the effects induced by the cuprizone toxin. In addition, the corresponding signals correlated strongly and significantly with the myelin content assessed by GFP fluorescence and LFB staining over the demyelination and the remyelination phases. The signal of the band-pass T_1D-filter, which isolates short T_1D components, showed changes over time that were poorly correlated with histology, hence suggesting a sensitivity to pathological processes possibly not related to myelin. Although MPF was also highly correlated to histology, ihMT high-pass T_1D-filters showed better capability to characterize the spatial-temporal patterns during the demyelination and remyelination phases of the acute cuprizone model (e.g., rostro-caudal gradient of demyelination in the mCC previously described in the literature).

CONCLUSIONS

IhMT sequences selective for long T_1D components are specific and sensitive in vivo markers of demyelination and remyelination and have successfully captured the spatially heterogeneous pattern of the demyelination and remyelination mechanisms in the cuprizone model. Interestingly, differences in signal variations between the ihMT high-pass and band-pass T_1D-filter, suggest a sensitivity of the ihMT sequences targeted to short T_1Ds to alterations other than those of myelin. Future studies will need to further address these differences by examining more closely the origin of the short T_1D components and the variation of each T_1D component in pathology.

Primary Multiparametric Quantitative Brain MRI: State-of-the-Art Relaxometric and Proton Density Mapping Techniques

This review on brain multiparametric quantitative MRI (MP-qMRI) focuses on the primary subset of quantitative MRI (qMRI) parameters that represent the mobile ("free") and bound ("motion-restricted") proton pools. Such primary parameters are the proton densities, relaxation times, and magnetization transfer parameters. Diffusion qMRI is also included because of its wide implementation in complete clinical MP-qMRI application. MP-qMRI advances were reviewed over the past 2 decades, with substantial progress observed toward accelerating image acquisition and increasing mapping accuracy. Areas that need further investigation and refinement are identified as follows: (a) the biologic underpinnings of qMRI parameter values and their changes with age and/or disease and (b) the theoretical limitations implicitly built into most qMRI mapping algorithms that do not distinguish between the different spatial scales of voxels versus spin packets, the central physical object of the Bloch theory. With rapidly improving image processing techniques and continuous advances in computer hardware, MP-qMRI has the potential for implementation in a wide range of clinical applications. Currently, three emerging MP-qMRI applications are synthetic MRI, macrostructural qMRI, and microstructural tissue modeling.

Generalized Bloch model: A theory for pulsed magnetization transfer

PURPOSE

The paper introduces a classical model to describe the dynamics of large spin-1/2 ensembles associated with nuclei bound in large molecule structures, commonly referred to as the semi-solid spin pool, and their magnetization transfer (MT) to spins of nuclei in water.

THEORY AND METHODS

Like quantum-mechanical descriptions of spin dynamics and like the original Bloch equations, but unlike existing MT models, the proposed model is based on the algebra of angular momentum in the sense that it explicitly models the rotations induced by radiofrequency (RF) pulses. It generalizes the original Bloch model to non-exponential decays, which are, for example, observed for semi-solid spin pools. The combination of rotations with non-exponential decays is facilitated by describing the latter as Green's functions, comprised in an integro-differential equation.

RESULTS

Our model describes the data of an inversion-recovery magnetization-transfer experiment with varying durations of the inversion pulse substantially better than established models. We made this observation for all measured data, but in particular for pulse durations smaller than 300 μs. Furthermore, we provide a linear approximation of the generalized Bloch model that reduces the simulation time by approximately a factor 15,000, enabling simulation of the spin dynamics caused by a rectangular RF-pulse in roughly 2 μs.

CONCLUSION

The proposed theory unifies the original Bloch model, Henkelman's steady-state theory for MT, and the commonly assumed rotation induced by hard pulses (i.e., strong and infinitesimally short applications of RF-fields) and describes experimental data better than previous models.

T1D -weighted ihMT imaging - Part II. Investigating the long- and short-T1D components correlation with myelin content. Comparison with R1 and the macromolecular proton fraction

PURPOSE

To investigate the long- and short-T_1D components correlation with myelin content using inhomogeneous magnetization transfer (ihMT) high-pass and band-pass T_1D -filters and to compare ihMT, R₁ , and the macromolecular proton fraction (MPF) for myelin specific imaging.

METHODS

The 3D ihMT rapid gradient echo (ihMTRAGE) sequences with increasing switching times (Δt) were used to derive ihMT high-pass T_1D -filters with increasing T_1D cutoff values and an ihMT band-pass T_1D -filter for components in the 100 µs to 1 ms range. 3D spoiled gradient echo quantitative MT (SPGR-qMT) protocols were used to derive R₁ and MPF maps. The specificity of R₁ , MPF, and ihMT T_1D -filters was evaluated by comparison with two histological reference techniques for myelin imaging.

RESULTS

The higher contribution of long-T_1D s as compared to the short components as Δt got longer led to an increase in the specificity to myelination. In contrast, focusing on the signal originating from a narrow range of short-T_1D s (< 1 ms) as isolated by the band-pass T_1D -filter led to lower specificity. In addition, the significantly lower r² correlation coefficient of the band-pass T_1D -filter suggests that the origin of short-T_1D components is mostly associated with non-myelin protons. Also, the important contribution of short-T_1D s to the estimated MPF, explains its low specificity to myelination as compared to the ihMT high-pass T_1D -filters.

CONCLUSION

Long-T_1D components imaging by means of ihMT high-pass T_1D -filters is proposed as an MRI biomarker for myelin content. Future studies should enable the investigation of the sensitivity of ihMT T_1D -filters for demyelinating processes.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.