Assessment of microstructural signal compartments across the corpus callosum using multi-echo gradient recalled echo at 7 T

A predictor-corrector phase unwrapping algorithm for temporally undersampled gradient-echo MRI

PURPOSE

To develop a method for unwrapping temporally undersampled and nonlinear gradient recalled echo (GRE) phase.

THEORY AND METHODS

Temporal unwrapping is performed as a sequential one step prediction of the echo phase, followed by a correction to the nearest integer wrap-count. A spatio-temporal extension of the 1D predictor corrector unwrapping (PCU) algorithm improves the prediction accuracy, and thereby maintains spatial continuity. The proposed method is evaluated using numerical phantom, physical phantom, and in vivo brain data at both 3 T and 9.4 T. The unwrapping performance is compared with the state-of-the-art temporal and spatial unwrapping algorithms, and the spatio-temporal iterative virtual-echo based Nyquist sampled (iVENyS) algorithm.

RESULTS

Simulation results showed significant reduction in unwrapping errors at higher echoes compared with the state-of-the-art algorithms. Similar to the iVENyS algorithm, the PCU algorithm was able to generate spatially smooth phase images for in vivo data acquired at 3 T and 9.4 T, bypassing the use of additional spatial unwrapping step. A key advantage over iVENyS algorithm is the superior performance of PCU algorithm at higher echoes.

CONCLUSION

PCU algorithm serves as a robust phase unwrapping method for temporally undersampled and nonlinear GRE phase, particularly in the presence of high field gradients.

Novel insights into axon diameter and myelin content in late childhood and adolescence

White matter microstructural development in late childhood and adolescence is driven predominantly by increasing axon density and myelin thickness. Ex vivo studies suggest that the increase in axon diameter drives developmental increases in axon density observed with pubertal onset. In this cross-sectional study, 50 typically developing participants aged 8-18 years were scanned using an ultra-strong gradient magnetic resonance imaging scanner. Microstructural properties, including apparent axon diameter $({d}_a)$, myelin content, and g-ratio, were estimated in regions of the corpus callosum. We observed age-related differences in ${d}_a$, myelin content, and g-ratio. In early puberty, males had larger ${d}_a$ in the splenium and lower myelin content in the genu and body of the corpus callosum, compared with females. Overall, this work provides novel insights into developmental, pubertal, and cognitive correlates of individual differences in apparent axon diameter and myelin content in the developing human brain.

Multi-compartment analysis of the complex gradient-echo signal quantifies myelin breakdown in premanifest Huntington's disease

White matter (WM) alterations have been identified as a relevant pathological feature of Huntington's disease (HD). Increasing evidence suggests that WM changes in this disorder are due to alterations in myelin-associated biological processes. Multi-compartmental analysis of the complex gradient-echo MRI signal evolution in WM has been shown to quantify myelin in vivo, therefore pointing to the potential of this technique for the study of WM myelin changes in health and disease. This study first characterized the reproducibility of metrics derived from the complex multi-echo gradient-recalled echo (mGRE) signal across the corpus callosum in healthy participants, finding highest reproducibility in the posterior callosal segment. Subsequently, the same analysis pipeline was applied in this callosal region in a sample of premanifest HD patients (n = 19) and age, sex and education matched healthy controls (n = 21). In particular, we focused on two myelin-associated derivatives: i. the myelin water signal fraction (fm), a parameter dependent on myelin content; and ii. The difference in frequency between myelin and intra-axonal water pools (Δω), a parameter dependent on the ratio between the inner and the outer axonal radii. fm was found to be lower in HD patients (β = -0.13, p = 0.03), while Δω did not show a group effect. Performance in tests of working memory, executive function, social cognition and movement was also assessed, and a greater age-related decline in executive function was detected in HD patients (β = -0.06, p = 0.006), replicating previous evidence of executive dysfunction in HD. Finally, the correlation between fm, executive function, and proximity to disease onset was explored in patients, and a positive correlation between executive function and fm was detected (r = 0.542; p = 0.02). This study emphasises the potential of complex mGRE signal analysis for aiding understanding of HD pathogenesis and progression. Moreover, expanding on evidence from pathology and animal studies, it provides novel in vivo evidence supporting myelin breakdown as an early feature of HD.

Probing the myelin water compartment with a saturation-recovery, multi-echo gradient-recalled echo sequence

PURPOSE

To investigate the effect of varying levels of T 1 -weighting on the evolution of the complex signal from white matter in a multi-echo gradient-recalled echo (mGRE) saturation-recovery sequence.

THEORY AND METHODS

Analysis of the complex signal evolution in an mGRE sequence allows the contributions from short- and long- T 2 ∗ components to be separated, thus providing a measure of the relative strength of signals from the myelin water, and the external and intra-axonal compartments. Here we evaluated the effect of different levels of T 1 -weighting on these signals, expecting that the previously reported, short T 1 of the myelin water would lead to a relative enhancement of the myelin water signal in the presence of signal saturation. Complex, saturation-recovery mGRE data from the splenium of the corpus callosum from 5 healthy volunteers were preprocessed using a frequency difference mapping (FDM) approach and analyzed using the 3-pool model of complex signal evolution in white matter.

RESULTS

An increase in the apparent T 1 as a function of echo time was demonstrated, but this increase was an order of magnitude smaller than that expected from previously reported myelin water T 1 -values. This suggests the presence of magnetization transfer and exchange effects which counteract the T 1 -weighting.

CONCLUSION

Variation of the B 1 + amplitude in a saturation-recovery mGRE sequence can be used to modulate the relative strength of signals from the different compartments in white matter, but the modulation is less than predicted from previously reported T 1 -values.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

Field strength influences on gradient recalled echo MRI signal compartment frequency shifts

Different trends of echo time dependent gradient recalled echo MRI signals in different brain regions have been attributed to signal compartments in image voxels. It remains unclear how variations in gradient recalled echo MRI signals change as a function of MRI field strength, and how data processing may impact signal compartment parameters. We used two popular quantitative susceptibility mapping methods of processing raw phase images (Laplacian and path-based unwrapping with V-SHARP) and expressed values in the form of induced frequency shifts (in Hz) in six specific brain regions at 3T and 7T. We found the frequency shift curves to vary with echo time, and a good overlap between 3T and 7T mean frequency shift curves was present. However, the amount of variation across participants was greater at 3T, and we were able to obtain better compartment model fits of the signal at 7T. We also found the temporal trends in the signal and compartment frequency shifts to change with the method used to process images. The inter-participant averaged trends were consistent between 3T and 7T for each quantitative susceptibility pipeline. However, signal compartment frequency shifts generated using different pipelines may not be comparable.

Influence of 7T GRE-MRI Signal Compartment Model Choice on Tissue Parameters

Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders. Ultra-high field magnetic resonance imaging (MRI) data obtained using a multi-echo gradient echo sequence have been shown to contain information on myelin, axonal, and extracellular compartments in tissue. Quantitative assessment of water fraction, relaxation time (T₂*), and frequency shift using multi-compartment models has been shown to be useful in studying white matter properties via specific tissue parameters. It remains unclear how tissue parameters vary with model selection based on 7T multiple echo time gradient-recalled echo (GRE) MRI data. We applied existing signal compartment models to the corpus callosum and investigated whether a three-compartment model can be reduced to two compartments and still resolve white matter parameters [i.e., myelin water fraction (MWF) and g-ratio]. We show that MWF should be computed using a three-compartment model in the corpus callosum, and the g-ratios obtained using three compartment models are consistent with previous reports. We provide results for other parameters, such as signal compartment frequency shifts.

Strong diffusion gradients allow the separation of intra- and extra-axonal gradient-echo signals in the human brain

The quantification of brain white matter properties is a key area of application of Magnetic Resonance Imaging (MRI), with much effort focused on using MR techniques to quantify tissue microstructure. While diffusion MRI probes white matter (WM) microstructure by characterising the sensitivity of Brownian motion of water molecules to anisotropic structures, susceptibility-based techniques probe the tissue microstructure by observing the effect of interaction between the tissue and the magnetic field. Here, we unify these two complementary approaches by combining ultra-strong (300 mT/m) gradients with a novel Diffusion-Filtered Asymmetric Spin Echo (D-FASE) technique. Using D-FASE we can separately assess the evolution of the intra- and extra-axonal signals under the action of susceptibility effects, revealing differences in the behaviour in different fibre tracts. We observed that the effective relaxation rate of the ASE signal in the corpus callosum decreases with increasing b-value in all subjects (from 17.1 ± 0.7 s⁻¹ at b = 0 s/mm² to 14.6 ± 0.7 s⁻¹ at b = 4800 s/mm²), while this dependence on b in the corticospinal tract is less pronounced (from 12.0± 1.1 s⁻¹ at b = 0s/mm² to 10.7 ± 0.5 s⁻¹ at b = 4800 s/mm²). Voxelwise analysis of the signal evolution with respect to b-factor and acquisition delay using a microscopic model demonstrated differences in gradient echo signal evolution between the intra- and extra-axonal pools.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 ​ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 ​T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Correlations of quantitative MRI metrics with myelin basic protein (MBP) staining in a murine model of demyelination

Myelin imaging in the central nervous system is essential for monitoring pathologies involving white matter alterations. Various quantitative MRI protocols relying on the modeling of the interactions of water protons with myelinated tissues have shown sensitivities in case of myelin disruption. Some extracted model parameters are more sensitive to demyelination, such as the bound pool fraction (f) in quantitative magnetization transfer imaging (qMTI), the radial diffusivity in diffusion tensor imaging (DTI), and the myelin water fraction (MWF) in myelin water imaging (MWI). A 3D ultrashort echo time (UTE) sequence within an appropriate water suppression condition (Diff-UTE) is also considered for the direct visualization of the myelin semi-solid matrix (Diff-UTE normalized signal; rSPF). In this paper, we aimed at assessing the sensitivities and correlations of the parameters mentioned above to an immuno-histological study of the myelin basic protein (MBP) in a murine model of demyelination at 7 T. We demonstrated a high sensitivity of the MRI metrics to demyelination, and strong Spearman correlations in the corpus callosum between histology, macromolecular proton fraction (ρ>0.87) and Diff-UTE signal (ρ>0.76), but moderate ones with radial diffusivity and MWF (|ρ|<0.70).

Frequency difference mapping applied to the corpus callosum at 7T

Purpose

Frequency difference mapping (FDM) is a phase processing technique which characterizes the nonlinear temporal evolution of the phase of gradient echo (GE) signals. Here, a novel FDM‐processing algorithm is introduced, which is shown to reveal information about white matter microstructure. Unlike some other phase‐processing techniques, the FDM algorithm presented here does not require the use of phase unwrapping or sophisticated image processing. It uses a series of scaled complex divisions to unwrap phase and remove background fields.

Methods

Ten healthy subjects underwent a series of single‐slice, sagittal multi‐echo GE scans at 7T with the slice positioned at the midline. Phase data were processed with the novel FDM algorithm, and the temporal evolution of the magnitude signal and frequency difference was examined in 5 regions of the corpus callosum (CC; genu, anterior body, middle body, posterior body, and splenium).

Results

Consistent frequency difference contrast relative to surrounding tissue was observed in all subjects in the CC and in other white matter regions where the nerve fibers run perpendicular to B0, such as the superior cerebellar peduncle. Examination of the frequency difference curves shows distinct variations over the CC, with the genu and splenium displaying larger frequency differences than the other regions (in addition to a faster decay of signal magnitude).

Conclusion

The novel FDM algorithm presented here yields images sensitive to tissue microstructure and microstructural differences over the CC in a simple manner, without the requirement for phase unwrapping or sophisticated image processing.

High-resolution imaging of distinct human corpus callosum microstructure and topography of structural connectivity to cortices at high field

Characterization of the microstructural properties and topography of the human corpus callosum (CC) is key to understanding interhemispheric neural communication and brain function. In this work, we tested the hypothesis that high-resolution T₁ relaxometry at high field has adequate sensitivity and specificity for characterizing microstructural properties of the human CC, and elucidating the structural connectivity of the callosal fibers to the cortices of origin. The high-resolution parametric T₁ images acquired from healthy subjects (N = 16) at 7 T clearly showed a consistent T₁ distribution among individuals with substantial variation along the human CC axis, which is highly similar to the spatial patterns of myelin density and myelinated axon size based on the histology study. Compared to the anterior part of the CC, the posterior midbody and splenium had significantly higher T₁ values. In conjunction with T₁-based classification method, the splenial T₁ values were decoded more reliably compared to a conventional partitioning method, showing a much higher T₁ value in the inferior splenium than in the middle/superior splenium. Moreover, the T₁ profile of the callosal subdivision represented the topology of the fiber connectivity to the projected cortical regions: the fibers in the posterior midbody and inferior splenium with a higher T₁ (inferring a larger axon size) were mainly connected to motor-sensory and visual cortical areas, respectively; in contrast, the fibers in the anterior/posterior CC with a lower T₁ (inferring a smaller axon size) were primarily connected to the frontal/parietal-temporal areas. These findings indicate that high-resolution T₁ relaxometry imaging could provide a complementary and robust neuroimaging tool, useful for exploring the complex tissue properties and topographic organization of the human corpus callosum.