Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

Susceptibility-induced internal gradients reveal axon morphology and cause anisotropic effects in the diffusion-weighted MRI signal

Diffusion-weighted MRI is our most promising method for estimating microscopic tissue morphology in vivo. The signal acquisition is based on scanner-generated external magnetic gradients. However, it will also be affected by susceptibility-induced internal magnetic gradients caused by interactions between the tissue and the static magnetic field of the scanner. With 3D in silico experiments, we show how internal gradients cause morphology-, compartment-, and orientation-dependence of spin-echo and pulsed-gradient spin-echo experiments in myelinated axons. These effects surpass those observed with previous 2D modelling corresponding to straight cylinders. For an ex vivo monkey brain, we observe the orientation-dependence generated only when including non-circular cross-sections in the in silico morphological configurations, and find orientation-dependent deviation of up to 17% for diffusion tensor metrics. Interestingly, we find that the orientation-dependence not only biases the signal across different brain regions, but also carries a sensitivity to the morphology of axonal cross-sections which is not attainable by the idealised theoretical diffusion-weighted MRI signal.

The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI

PURPOSE

This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique.

THEORY AND METHODS

A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images.

RESULTS

Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes.

CONCLUSION

The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.

Optimisation of quantitative brain diffusion-relaxation MRI acquisition protocols with physics-informed machine learning

Diffusion-relaxation MRI aims to extract quantitative measures that characterise microstructural tissue properties such as orientation, size, and shape, but long acquisition times are typically required. This work proposes a physics-informed learning framework to extract an optimal subset of diffusion-relaxation MRI measurements for enabling shorter acquisition times, predict non-measured signals, and estimate quantitative parameters. In vivo and synthetic brain 5D-Diffusion-T1-T2∗-weighted MRI data obtained from five healthy subjects were used for training and validation, and from a sixth participant for testing. One fully data-driven and two physics-informed machine learning methods were implemented and compared to two manual selection procedures and Cramér-Rao lower bound optimisation. The physics-informed approaches could identify measurement-subsets that yielded more consistently accurate parameter estimates in simulations than other approaches, with similar signal prediction error. Five-fold shorter protocols yielded error distributions of estimated quantitative parameters with very small effect sizes compared to estimates from the full protocol. Selected subsets commonly included a denser sampling of the shortest and longest inversion time, lowest echo time, and high b-value. The proposed framework combining machine learning and MRI physics offers a promising approach to develop shorter imaging protocols without compromising the quality of parameter estimates and signal predictions.

Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction (f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction (f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation off AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics off A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmentalT 2 and can substantially enhance the comparability between EM- and DWI-based metrics off A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-basedf A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter, we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a 2-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of free water fractions while also capturing compartmental T 2 relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.

Evaluation of tractography-based myelin-weighted connectivity across the lifespan

Introduction

Recent studies showed that the myelin of the brain changes in the life span, and demyelination contributes to the loss of brain plasticity during normal aging. Diffusion-weighted magnetic resonance imaging (dMRI) allows studying brain connectivity in vivo by mapping axons in white matter with tractography algorithms. However, dMRI does not provide insight into myelin; thus, combining tractography with myelin-sensitive maps is necessary to investigate myelin-weighted brain connectivity. Tractometry is designated for this purpose, but it suffers from some serious limitations. Our study assessed the effectiveness of the recently proposed Myelin Streamlines Decomposition (MySD) method in estimating myelin-weighted connectomes and its capacity to detect changes in myelin network architecture during the process of normal aging. This approach opens up new possibilities compared to traditional Tractometry.

Methods

In a group of 85 healthy controls aged between 18 and 68 years, we estimated myelin-weighted connectomes using Tractometry and MySD, and compared their modulation with age by means of three well-known global network metrics.

Results

Following the literature, our results show that myelin development continues until brain maturation (40 years old), after which degeneration begins. In particular, mean connectivity strength and efficiency show an increasing trend up to 40 years, after which the process reverses. Both Tractometry and MySD are sensitive to these changes, but MySD turned out to be more accurate.

Conclusion

After regressing the known predictors, MySD results in lower residual error, indicating that MySD provides more accurate estimates of myelin-weighted connectivity than Tractometry.

MR myelin imaging in multiple sclerosis: A scoping review

The inability of disease-modifying therapies to stop the progression of multiple sclerosis (MS), has led to the development of a new therapeutic strategy focussing on myelin repair. While conventional MRI lacks sensitivity for quantifying myelin damage, advanced MRI techniques are proving effective. The development of targeted therapeutics requires histological validation of myelin imaging results, alongside the crucial task of establishing correlations between myelin imaging results and clinical assessments, so that the effectiveness of therapeutic interventions can be evaluated. The aims of this scoping review were to identify myelin imaging methods - some of which have been histologically validated, and to determine how these approaches correlate with clinical assessments of people with MS (pwMS), thus allowing for effective therapeutic evaluation. A search of two databases was undertaken for publications relating to studies on adults MS using either MRI/MR-histology of the MS brain in the range 1990-to-2022. The myelin imaging methods specified were relaxometry, magnetization transfer, and quantitative susceptibility. Relaxometry was used most frequently, with myelin water fraction (MWF) being the primary metric. Studies conducted on tissue from various regions of the brain showed that MWF was significantly lower in pwMS than in healthy controls. Magnetization transfer ratio indicated that the macromolecular content of lesions was lower than that of normal-appearing tissue. Higher magnetic susceptibility of lesions were indicative of myelin breakdown and iron accumulation. Several myelin imaging metrics were correlated with disability, disease severity and duration. Many studies showed a good correlation between myelin measured histologically and by MR myelin imaging techniques.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

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

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.

Phase-shifted transverse relaxation orientation dependences in human brain white matter

This work aimed to demonstrate an essential phase shift ε 0 $$ {\varepsilon}_0 $$ for better quantifying R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ in human brain white matter (WM), and to further elucidate its origin related to the directional diffusivities from standard diffusion tensor imaging (DTI). ε 0 $$ {\varepsilon}_0 $$ was integrated into a proposed generalized transverse relaxation model for characterizing previously published R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ orientation dependence profiles in brain WM, and then comparisons were made with those without ε 0 $$ {\varepsilon}_0 $$ . It was theorized that anisotropic diffusivity direction ε $$ \varepsilon $$ was collinear with an axon fiber subject to all eigenvalues and eigenvectors from an apparent diffusion tensor. To corroborate the origin of ε 0 $$ {\varepsilon}_0 $$ , R 2 $$ {R}_2 $$ orientation dependences referenced by ε $$ \varepsilon $$ were compared with those referenced by the standard principal diffusivity direction Φ $$ \Phi $$ at b-values of 1000 and 2500 (s/mm² ). These R 2 $$ {R}_2 $$ orientation dependences were obtained from T 2 $$ {T}_2 $$ -weighted images (b = 0) of ultrahigh-resolution Connectome DTI datasets in the public domain. A normalized root-mean-square error ( NRMSE % $$ NRMSE\% $$ ) and an F $$ F $$ -test were used for evaluating curve-fittings, and statistical significance was considered to be a p of 0.05 or less. A phase-shifted model resulted in significantly reduced NRMSE % $$ NRMSE\% $$ compared with that without ε 0 $$ {\varepsilon}_0 $$ in quantifying various R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ profiles, both in vivo and ex vivo at multiple B 0 $$ {B}_0 $$ fields. The R 2 $$ {R}_2 $$ profiles based on Φ $$ \Phi $$ manifested a right-shifted phase ( ε 0 > 0 $$ {\varepsilon}_0>0 $$ ) at two b-values, while those based on ε $$ \varepsilon $$ became free from ε 0 $$ {\varepsilon}_0 $$ . For all phase-shifted R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ profiles, ε 0 $$ {\varepsilon}_0 $$ generally depended on the directional diffusivities by tan - 1 D ⊥ / D ∥ $$ {\tan}^{-1}\left({D}_{\perp }/{D}_{\parallel}\right) $$ , as predicted. In summary, a ubiquitous phase shift ε 0 $$ {\varepsilon}_0 $$ has been demonstrated as a prerequisite for better quantifying transverse relaxation orientation dependences in human brain WM. Furthermore, the origin of ε 0 $$ {\varepsilon}_0 $$ associated with the directional diffusivities from DTI has been elucidated. These findings could have a significant impact on interpretations of prior R 2 $$ {R}_2 $$ and R 2 * $$ {R}_2^{\ast } $$ datasets and on future research.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

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

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

Orientation dependence of R2 relaxation in the newborn brain

In MRI the transverse relaxation rate, R₂ = 1/T₂^, shows dependence on the orientation of ordered tissue relative to the main magnetic field. In previous studies, orientation effects of R₂ relaxation in the mature brain's white matter have been found to be described by a susceptibility-based model of diffusion through local magnetic field inhomogeneities created by the diamagnetic myelin sheaths. Orientation effects in human newborn white matter have not yet been investigated. The newborn brain is known to contain very little myelin and is therefore expected to exhibit a decrease in orientation dependence driven by susceptibility-based effects. We measured R₂ orientation dependence in the white matter of human newborns. R₂ data were acquired with a 3D Gradient and Spin Echo (GRASE) sequence and fiber orientation was mapped with diffusion tensor imaging (DTI). We found orientation dependence in newborn white matter that is not consistent with the susceptibility-based model and is best described by a model of residual dipolar coupling. In the near absence of myelin in the newborn brain, these findings suggest the presence of residual dipolar coupling between rotationally restricted water molecules. This has important implications for quantitative imaging methods such as myelin water imaging, and suggests orientation dependence of R₂ as a potential marker in early brain development.

Magic angle effect on diffusion tensor imaging in ligament and brain

PURPOSE

To evaluate the magic angle effect on diffusion tensor imaging (DTI) measurements in rat ligaments and mouse brains.

METHODS

Three rat knee joints and three mouse brains were scanned at 9.4 T using a modified 3D diffusion-weighted spin echo pulse sequence with the isotropic spatial resolution of 45 μm. The b value was 1000 s/mm2 for rat knee and 4000 s/mm2 for mouse brain. DTI model was used to investigate the quantitative metrics at different orientations with respect to the main magnetic field. The collagen fiber structure of the ligament was validated with polarized light microscopy (PLM) imaging.

RESULTS

The signal intensity, signal-to-noise ratio (SNR), and DTI metrics in the ligament were strongly dependent on the collagen fiber orientation with respect to the main magnetic field from both simulation and actual MRI scans. The variation of fractional anisotropy (FA) was about ~32%, and the variation of mean diffusivity (MD) was ~11%. These findings were further validated with the numerical simulation at different SNRs (~10.0 to 86.0). Compared to the ligament, the DTI metrics showed little orientation dependence in mouse brains.

CONCLUSION

Magic angle effect plays an important role in DTI measurements in the highly ordered collagen-rich tissues, while MD showed less orientation dependence than FA.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T₂-DKI-FWE model that exploits the T₂ relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T₂ of tissue). In our approach, the T₂ of tissue is estimated as an unknown parameter, whereas the T₂ of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T₂ of free water is studied. Next, the improved conditioning of T₂-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T₂-DKI-FWE model is compared to that of the DKI-FWE and T₂-DKI models on both simulated and real datasets. The error due to a biased approximation of the T₂ of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T₂-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T₂-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T₂ relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Bundle myelin fraction (BMF) mapping of different white matter connections using microstructure informed tractography

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, "tractometry" can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named "myelin streamline decomposition" (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T₂ that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T₂, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T₂ based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T₂ estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

MRI of healthy brain aging: A review

We present a review of the characterization of healthy brain aging using MRI with an emphasis on morphology, lesions, and quantitative MR parameters. A scope review found 6612 articles encompassing the keywords "Brain Aging" and "Magnetic Resonance"; papers involving functional MRI or not involving imaging of healthy human brain aging were discarded, leaving 2246 articles. We first consider some of the biogerontological mechanisms of aging, and the consequences of aging in terms of cognition and onset of disease. Morphological changes with aging are reviewed for the whole brain, cerebral cortex, white matter, subcortical gray matter, and other individual structures. In general, volume and cortical thickness decline with age, beginning in mid-life. Prevalent silent lesions such as white matter hyperintensities, microbleeds, and lacunar infarcts are also observed with increasing frequency. The literature regarding quantitative MR parameter changes includes T₁ , T₂ , T₂ *, magnetic susceptibility, spectroscopy, magnetization transfer, diffusion, and blood flow. We summarize the findings on how each of these parameters varies with aging. Finally, we examine how the aforementioned techniques have been used for age prediction. While relatively large in scope, we present a comprehensive review that should provide the reader with sound understanding of what MRI has been able to tell us about how the healthy brain ages.