Imaging brain microstructure with diffusion MRI: practicality and applications

Leveraging multi-shell diffusion for studies of brain development in youth and young adulthood

•

Multi-shell imaging sequences may improve sensitivity to developmental effects.

•

Models that leverage multi-shell information are often less sensitive to the confounding effects of motion.

•

Multi-shell sequences and models that leverage this data may be of particular utility for studying the developing brain.

Diffusion weighted imaging (DWI) has advanced our understanding of brain microstructure evolution over development. Recently, the use of multi-shell diffusion imaging sequences has coincided with advances in modeling the diffusion signal, such as Neurite Orientation Dispersion and Density Imaging (NODDI) and Laplacian-regularized Mean Apparent Propagator MRI (MAPL). However, the relative utility of recently-developed diffusion models for understanding brain maturation remains sparsely investigated. Additionally, despite evidence that motion artifact is a major confound for studies of development, the vulnerability of metrics derived from contemporary models to in-scanner motion has not been described. Accordingly, in a sample of 120 youth and young adults (ages 12–30) we evaluated metrics derived from diffusion tensor imaging (DTI), NODDI, and MAPL for associations with age and in-scanner head motion at multiple scales. Specifically, we examined mean white matter values, white matter tracts, white matter voxels, and connections in structural brain networks. Our results revealed that multi-shell diffusion imaging data can be leveraged to robustly characterize neurodevelopment, and demonstrate stronger age effects than equivalent single-shell data. Additionally, MAPL-derived metrics were less sensitive to the confounding effects of head motion. Our findings suggest that multi-shell imaging data and contemporary modeling techniques confer important advantages for studies of neurodevelopment.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b ​≤ ​3,000 ​s/mm² (or 3 ​ms/μm²), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm² or 3 ​ms/μm²). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (b_max ​= ​40,000 ​s/mm², or 40 ​ms/μm²) and human (b_max ​= ​10,000 ​s/mm², or 10 ​ms/μm²) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

The role of diffusion tractography in refining glial tumor resection

Primary brain tumors are notoriously hard to resect surgically. Due to their infiltrative nature, finding the optimal resection boundary without damaging healthy tissue can be challenging. One potential tool to help make this decision is diffusion-weighted magnetic resonance imaging (dMRI) tractography. dMRI exploits the diffusion of water molecule along axons to generate a 3D modelization of the white matter bundles in the brain. This feature is particularly useful to visualize how a tumor affects its surrounding white matter and plan a surgical path. This paper reviews the different ways in which dMRI can be used to improve brain tumor resection, its benefits and also its limitations. We expose surgical tools that can be paired with dMRI to improve its impact on surgical outcome, such as loading the 3D tractography in the neuronavigation system and direct electrical stimulation to validate the position of the white matter bundles of interest. We also review articles validating dMRI findings using other anatomical investigation techniques, such as postmortem dissections, manganese-enhanced MRI, electrophysiological stimulations, and phantom studies with known ground truth. We will be discussing the areas of the brain where dMRI performs well and where the future challenges are. We will conclude this review with suggestions and take home messages for neurosurgeons, tractographers, and vendors for advancing the field and on how to benefit from tractography's use in clinical practice.

Double Diffusion Encoding for Probing Radiation-Induced Microstructural Changes in a Tumor Model: A Proof-of-Concept Study With Comparison to the Apparent Diffusion Coefficient and Histology

BACKGROUND

Microstructure analyses are gaining interest in cancer MRI as an alternative to the conventional apparent diffusion coefficient (ADC), of which the determinants remain unclear.

PURPOSE

To assess the sensitivity of parameters calculated from a double diffusion encoding (DDE) sequence to changes in a tumor's microstructure early after radiotherapy and to compare them with ADC and histology.

STUDY TYPE

Cohort study on experimental tumors.

ANIMAL MODEL

Sixteen WAG/Rij rats grafted with one rhabdomyosarcoma fragment in each thigh. Thirty-one were imaged at days 1 and 4, of which 17 tumors received a 20 Gy radiation dose after the first imagery.

FIELD STRENGTH/SEQUENCE

3T. Diffusion-weighted imaging, DDE with flow compensated, and noncompensated measurements.

ASSESSMENTS

1) To compare, after irradiation, DDE-derived parameters (intracellular fraction, cell size, and cell density) to their histological counterparts (fraction of stained area, minimal Feret diameter, and nuclei count, respectively). 2) To compare percentage changes in DDE-derived parameters and ADC. 3) To evaluate the evolution of DDE-derived parameters describing perfusion.

STATISTICAL TESTS

Wilcoxon rank sum test.

RESULTS

1) Intracellular fraction, cell size, and cell density were respectively lower (-24%, P < 0.001), higher (+7.5%, P < 0.001) and lower (-38%, P < 0.001) in treated tumors as compared to controls. Fraction of stained area, minimal Feret diameter, and nuclei count were respectively lower (-20%, P < 0.001), higher (+28%, P < 0.001), and lower (-34%, P < 0.001) in treated tumors. 2) The magnitude of ADC's percentage change due to irradiation (16.4%) was superior to the one of cell size (8.4%, P < 0.01) but inferior to those of intracellular fraction (35.5%, P < 0.001) and cell density (42%, P < 0.001). 3) After treatment, the magnitude of the vascular fraction's decrease was higher than the increase of flow velocity (33.3%, vs. 13.3%, P < 0.001).

DATA CONCLUSION

The DDE sequence allows quantitatively monitoring the effects of radiotherapy on a tumor's microstructure, whereas ADC only reveals global changes.

EVIDENCE LEVEL

2.

TECHNICAL EFFICACY

Stage 4. J. Magn. Reson. Imaging 2020;52:941-951.

Towards unconstrained compartment modeling in white matter using diffusion-relaxation MRI with tensor-valued diffusion encoding

PURPOSE

To optimize diffusion-relaxation MRI with tensor-valued diffusion encoding for precise estimation of compartment-specific fractions, diffusivities, and T₂ values within a two-compartment model of white matter, and to explore the approach in vivo.

METHODS

Sampling protocols featuring different b-values (b), b-tensor shapes (b_Δ ), and echo times (TE) were optimized using Cramér-Rao lower bounds (CRLB). Whole-brain data were acquired in children, adults, and elderly with white matter lesions. Compartment fractions, diffusivities, and T₂ values were estimated in a model featuring two microstructural compartments represented by a "stick" and a "zeppelin."

RESULTS

Precise parameter estimates were enabled by sampling protocols featuring seven or more "shells" with unique b/b_Δ /TE-combinations. Acquisition times were approximately 15 minutes. In white matter of adults, the "stick" compartment had a fraction of approximately 0.5 and, compared with the "zeppelin" compartment, featured lower isotropic diffusivities (0.6 vs. 1.3 μm² /ms) but higher T₂ values (85 vs. 65 ms). Children featured lower "stick" fractions (0.4). White matter lesions exhibited high "zeppelin" isotropic diffusivities (1.7 μm² /ms) and T₂ values (150 ms).

CONCLUSIONS

Diffusion-relaxation MRI with tensor-valued diffusion encoding expands the set of microstructure parameters that can be precisely estimated and therefore increases their specificity to biological quantities.

Self-navigation for 3D multishot EPI with data-reference

PURPOSE

In this study, we sought to develop a self-navigation strategy for improving the reconstruction of diffusion weighted 3D multishot echo planar imaging (EPI). We propose a method for extracting the phase correction information from the acquisition itself, eliminating the need for a 2D navigator, further accelerating the acquisition.

METHODS

In-vivo acquisitions at 3T with 0.9 mm and 1.5 mm isotropic resolutions were used to evaluate the performance of the self-navigation strategy. Sensitivity to motion was tested using a large difference in pitch position of the head. Using a multishell diffusion weighted acquisition, tractography results were obtained at (0.9 mm)³ to validate the quality with conventional acquisition.

RESULTS

The use of 3D multislab EPI with self-navigation enables 3D diffusion-weighted spin echo EPI acquisitions that have the same efficiency as 2D single-shot acquisition. For matched acquisition time the image signal-to-noise ratio (SNR) between 3D and 2D acquisition is shown to be comparable for whole-brain coverage with (1.5 mm)³ resolution and for (0.9 mm)³ resolution the 3D acquisition has higher SNR than what can be obtained with 2D acquisitions using current state-of-art multiband techniques. The self-navigation technique was shown to be stable under inter-volume motion. In tractography analysis, the higher resolution afforded by our technique enabled clear delineation of the tapetum and posterior corona radiata.

CONCLUSION

The proposed self-navigation approach utilized a self-consistent phase in 3D diffusion weighted acquisitions. Its efficiency and stability were demonstrated for a plurality of common acquisitions. The proposed self-navigation approach allows for faster acquisition of 3D multishot EPI desirable for large field of view and/or higher resolution.

Time-dependent diffusion in undulating thin fibers: Impact on axon diameter estimation

Diffusion MRI may enable non-invasive mapping of axonal microstructure. Most approaches infer axon diameters from effects of time-dependent diffusion on the diffusion-weighted MR signal by modeling axons as straight cylinders. Axons do not, however, propagate in straight trajectories, and so far the impact of the axonal trajectory on diameter estimation has been insufficiently investigated. Here, we employ a toy model of axons, which we refer to as the undulating thin fiber model, to analyze the impact of undulating trajectories on the time dependence of diffusion. We study time-dependent diffusion in the frequency domain and characterize the diffusion spectrum by its height, width, and low-frequency behavior (power law exponent). Results show that microscopic orientation dispersion of the thin fibers is the main parameter that determines the characteristics of the diffusion spectra. At lower frequencies (longer diffusion times), straight cylinders and undulating thin fibers can have virtually identical spectra. If the straight-cylinder assumption is used to interpret data from undulating thin axons, the diameter is overestimated by an amount proportional to the undulation amplitude and microscopic orientation dispersion of the fibers. At higher frequencies (shorter diffusion times), spectra from cylinders and undulating thin fibers differ. The low-frequency behavior of the spectra from the undulating thin fibers may also differ from that of cylinders, because the power law exponent of undulating fibers can reach values below 2 for experimentally relevant frequency ranges. In conclusion, we argue that the non-straight nature of axonal trajectories should not be overlooked when analyzing and interpreting diffusion MRI data.

Revisiting double diffusion encoding MRS in the mouse brain at 11.7T: Which microstructural features are we sensitive to?

Brain metabolites, such as N-acetylaspartate or myo-inositol, are constantly probing their local cellular environment under the effect of diffusion. Diffusion-weighted NMR spectroscopy therefore presents unparalleled potential to yield cell-type specific microstructural information. Double diffusion encoding (DDE) consists in applying two diffusion blocks, where gradient's direction in the second block is varied during the course of the experiment. Unlike single diffusion encoding, DDE measurements at long mixing time display some angular modulation of the signal amplitude which reflects microscopic anisotropy (μA), while requiring relatively low gradient strength. This angular dependence has been formerly used to quantify cell fiber diameter using a model of isotropically oriented infinite cylinders. However, how additional features of the cell microstructure (such as cell body diameter, fiber length and branching) may also influence the DDE signal has been little explored. Here, we used a cryoprobe as well as state-of-the-art post-processing to perform DDE acquisitions with high accuracy and precision in the mouse brain at 11.7 ​T. We then compared our results to simulated DDE datasets obtained in various 3D cell models in order to pinpoint which features of cell morphology may influence the most the angular dependence of the DDE signal. While the infinite cylinder model poorly fits our experimental data, we show that incorporating branched fiber structure in our model allows more realistic interpretation of the DDE signal. Lastly, data acquired in the short mixing time regime suggest that some sensitivity to cell body diameter might be retrieved, although additional experiments would be required to further support this statement.

Constant gradient FEXSY: A time-efficient method for measuring exchange

Filter-Exchange NMR Spectroscopy (FEXSY) is a method for measurement of apparent transmembranal water exchange rates. The experiment is comprised of two co-linear sequential pulsed-field gradient (PFG) blocks, separated by a mixing period in which exchange takes place. The first block remains constant and serves as a diffusion-based filter that removes signal coming from fast-diffusing water. The mixing time and the gradient area (q-value) of the second block are varied on repeated iterations to produce a 2D data set that is analyzed using a bi-compartmental model which assumes that intra- and extra-cellular water are slow and fast diffusing, respectively. Here we suggest a variant of the FEXSY method in which measurements for different mixing times are taken at a constant gradient. This Constant Gradient FEXSY (CG-FEXSY) allows for the determination of the exchange rate by using a smaller 1D data set, so that the same information can be gathered during a considerably shorter scan time. Furthermore, in the limit of high diffusion weighting, such that the extra-cellular water signal is removed while the intra-cellular signal is retained, CG-FEXSY also allows for determination of the intra-cellular mean residence time (MRT). The theoretical results are validated on a living yeast cells sample and on a fixed porcine optic nerve, where the values obtained from the two methods are shown to be in agreement.

An improved deep network for tissue microstructure estimation with uncertainty quantification

Deep learning based methods have improved the estimation of tissue microstructure from diffusion magnetic resonance imaging (dMRI) scans acquired with a reduced number of diffusion gradients. These methods learn the mapping from diffusion signals in a voxel or patch to tissue microstructure measures. In particular, it is beneficial to exploit the sparsity of diffusion signals jointly in the spatial and angular domains, and the deep network can be designed by unfolding iterative processes that adaptively incorporate historical information for sparse reconstruction. However, the number of network parameters is huge in such a network design, which could increase the difficulty of network training and limit the estimation performance. In addition, existing deep learning based approaches to tissue microstructure estimation do not provide the important information about the uncertainty of estimates. In this work, we continue the exploration of tissue microstructure estimation using a deep network and seek to address these limitations. First, we explore the sparse representation of diffusion signals with a separable spatial-angular dictionary and design an improved deep network for tissue microstructure estimation. The procedure for updating the sparse code associated with the separable dictionary is derived and unfolded to construct the deep network. Second, with the formulation of sparse representation of diffusion signals, we propose to quantify the uncertainty of network outputs with a residual bootstrap strategy. Specifically, because of the sparsity constraint in the signal representation, we perform a Lasso bootstrap strategy for uncertainty quantification. Experiments were performed on brain dMRI scans with a reduced number of diffusion gradients, where the proposed method was applied to two representative biophysical models for describing tissue microstructure and compared with state-of-the-art methods of tissue microstructure estimation. The results show that our approach compares favorably with the competing methods in terms of estimation accuracy. In addition, the uncertainty measures provided by our method correlate with estimation errors and produce reasonable confidence intervals; these results suggest potential application of the proposed uncertainty quantification method in brain studies.

Mild traumatic brain injury impacts associations between limbic system microstructure and post-traumatic stress disorder symptomatology

BACKGROUND

Post-traumatic stress disorder (PTSD) is a psychiatric disorder that afflicts many individuals, yet the neuropathological mechanisms that contribute to this disorder remain to be fully determined. Moreover, it is unclear how exposure to mild traumatic brain injury (mTBI), a condition that is often comorbid with PTSD, particularly among military personnel, affects the clinical and neurological presentation of PTSD. To address these issues, the present study explores relationships between PTSD symptom severity and the microstructure of limbic and paralimbic gray matter brain regions, as well as the impact of mTBI comorbidity on these relationships.

METHODS

Structural and diffusion MRI data were acquired from 102 male veterans who were diagnosed with current PTSD. Diffusion data were analyzed with free-water imaging to quantify average CSF-corrected fractional anisotropy (FA) and mean diffusivity (MD) in 18 limbic and paralimbic gray matter regions. Associations between PTSD symptom severity and regional average dMRI measures were examined with repeated measures linear mixed models. Associations were studied separately in veterans with PTSD only, and in veterans with PTSD and a history of military mTBI.

RESULTS

Analyses revealed that in the PTSD only cohort, more severe symptoms were associated with higher FA in the right amygdala-hippocampus complex, lower FA in the right cingulate cortex, and lower MD in the left medial orbitofrontal cortex. In the PTSD and mTBI cohort, more severe PTSD symptoms were associated with higher FA bilaterally in the amygdala-hippocampus complex, with higher FA bilaterally in the nucleus accumbens, with lower FA bilaterally in the cingulate cortex, and with higher MD in the right amygdala-hippocampus complex.

CONCLUSIONS

These findings suggest that the microstructure of limbic and paralimbic brain regions may influence PTSD symptomatology. Further, given the additional associations observed between microstructure and symptom severity in veterans with head trauma, we speculate that mTBI may exacerbate the impact of brain microstructure on PTSD symptoms, especially within regions of the brain known to be vulnerable to chronic stress. A heightened sensitivity to the microstructural environment of the brain could partially explain why individuals with PTSD and mTBI comorbidity experience more severe symptoms and poorer illness prognoses than those without a history of brain injury. The relevance of these microstructural findings to the conceptualization of PTSD as being a disorder of stress-induced neuronal connectivity loss is discussed.

Self-Immolative Activation of β-Galactosidase-Responsive Probes for In Vivo MR Imaging in Mouse Models

Our lab has developed a new series of self-immolative MR agents for the rapid detection of enzyme activity in mouse models expressing β-galactosidase (β-gal). We investigated two molecular architectures to create agents that detect β-gal activity by modulating the coordination of water to GdIII . The first is an intermolecular approach, wherein we designed several structural isomers to maximize coordination of endogenous carbonate ions. The second involves an intramolecular mechanism for q modulation. We incorporated a pendant coordinating carboxylate ligand with a 2, 4, 6, or 8 carbon linker to saturate ligand coordination to the GdIII ion. This renders the agent ineffective. We show that one agent in particular (6-C pendant carboxylate) is an extremely effective MR reporter for the detection of enzyme activity in a mouse model expressing β-gal.

Age-Related Changes of Peak Width Skeletonized Mean Diffusivity (PSMD) Across the Adult Lifespan: A Multi-Cohort Study

Parameters of water diffusion in white matter derived from diffusion-weighted imaging (DWI), such as fractional anisotropy (FA), mean, axial, and radial diffusivity (MD, AD, and RD), and more recently, peak width of skeletonized mean diffusivity (PSMD), have been proposed as potential markers of normal and pathological brain ageing. However, their relative evolution over the entire adult lifespan in healthy individuals remains partly unknown during early and late adulthood, and particularly for the PSMD index. Here, we gathered and analyzed cross-sectional diffusion tensor imaging (DTI) data from 10 population-based cohort studies in order to establish the time course of white matter water diffusion phenotypes from post-adolescence to late adulthood. DTI data were obtained from a total of 20,005 individuals aged 18.1 to 92.6 years and analyzed with the same pipeline for computing skeletonized DTI metrics from DTI maps. For each individual, MD, AD, RD, and FA mean values were computed over their FA volume skeleton, PSMD being calculated as the 90% peak width of the MD values distribution across the FA skeleton. Mean values of each DTI metric were found to strongly vary across cohorts, most likely due to major differences in DWI acquisition protocols as well as pre-processing and DTI model fitting. However, age effects on each DTI metric were found to be highly consistent across cohorts. RD, MD, and AD variations with age exhibited the same U-shape pattern, first slowly decreasing during post-adolescence until the age of 30, 40, and 50 years, respectively, then progressively increasing until late life. FA showed a reverse profile, initially increasing then continuously decreasing, slowly until the 70s, then sharply declining thereafter. By contrast, PSMD constantly increased, first slowly until the 60s, then more sharply. These results demonstrate that, in the general population, age affects PSMD in a manner different from that of other DTI metrics. The constant increase in PSMD throughout the entire adult life, including during post-adolescence, indicates that PSMD could be an early marker of the ageing process.

Modified connectivity of vulnerable brain nodes in multiple sclerosis, their impact on cognition and their discriminative value

Brain structural network modifications in multiple sclerosis (MS) seem to be clinically relevant. The discriminative ability of those changes to identify MS patients or their cognitive status remains unknown. Therefore, this study aimed to investigate connectivity changes in MS patients related to their cognitive status, and to define an automatic classification method to classify subjects as patients and healthy volunteers (HV) or as cognitively preserved (CP) and impaired (CI) patients. We analysed structural brain connectivity in 45 HV and 188 MS patients (104 CP and 84 CI). A support vector machine with k-fold cross-validation was built using the graph metrics features that best differentiate the groups (p < 0.05). Local efficiency (LE) and node strength (NS) network properties showed the largest differences: 100% and 69.7% of nodes had reduced LE and NS in CP patients compared to HV. Moreover, 55.3% and 57.9% of nodes had decreased LE and NS in CI compared to CP patients, in associative multimodal areas. The classification method achieved an accuracy of 74.8–77.2% to differentiate patients from HV, and 59.9–60.8% to discriminate CI from CP patients. Structural network integrity is widely reduced and worsens as cognitive function declines. Central network properties of vulnerable nodes can be useful to classify MS patients.

A comprehensive approach for correcting voxel-wise b-value errors in diffusion MRI

PURPOSE

In diffusion MRI, the actual b-value played out on the scanner may deviate from the nominal value due to magnetic field imperfections. A simple image-based correction method for this problem is presented.

METHODS

The apparent diffusion constant (ADC) of a water phantom was measured voxel-wise along 64 diffusion directions at b = 1000 s/mm² . The true diffusion constant of water was estimated, considering the phantom temperature. A voxel-wise correction factor, providing an effective b-value including any magnetic field deviations, was determined for each diffusion direction by relating the measured ADC to the true diffusion constant. To test the method, the measured b-value map was used to calculate the corrected voxel-wise ADC for additionally acquired diffusion data sets on the same water phantom and data sets acquired on a small water phantom at three different positions. Diffusion tensor was estimated by applying the measured b-value map to phantom and in vivo data sets.

RESULTS

The b-value-corrected ADC maps of the phantom showed the expected spatial uniformity as well as a marked improvement in consistency across diffusion directions. The b-value correction for the brain data resulted in a 5.8% and 5.5% decrease in mean diffusivity and angular differences of the primary diffusion direction of 2.71° and 0.73° inside gray and white matter, respectively.

CONCLUSION

The actual b-value deviates significantly from its nominal setting, leading to a spatially variable error in the common diffusion outcome measures. The suggested method measures and corrects these artifacts.

On the need for bundle-specific microstructure kernels in diffusion MRI

Probing microstructure with diffusion magnetic resonance imaging (dMRI) on a scale orders of magnitude below the imaging resolution relies on biophysical modelling of the signal response in the tissue. The vast majority of these biophysical models of diffusion in white matter assume that the measured dMRI signal is the sum of the signals emanating from each of the constituent compartments, each of which exhibits a distinct behaviour in the b-value and/or orientation domain. Many of these models further assume that the dMRI behaviour of the oriented compartments (e.g. the intra-axonal space) is identical between distinct fibre populations, at least at the level of a single voxel. This implicitly assumes that any potential biological differences between fibre populations are negligible, at least as far as is measurable using dMRI. Here, we validate this assumption by means of a voxel-wise, model-free signal decomposition that, under the assumption above and in the absence of noise, is shown to be rank-1. We evaluate the effect size of signal components beyond this rank-1 representation and use permutation testing to assess their significance. We conclude that in the healthy adult brain, the dMRI signal is adequately represented by a rank-1 model, implying that biologically more realistic, but mathematically more complex fascicle-specific microstructure models do not capture statistically significant or anatomically meaningful structure, even in extended high-b diffusion MRI scans.

Single- and double-Diffusion encoding MRI for studying ex vivo apparent axon diameter distribution in spinal cord white matter

Mapping average axon diameter (AAD) and axon diameter distribution (ADD) in neuronal tissues non-invasively is a challenging task that may have a tremendous effect on our understanding of the normal and diseased central nervous system (CNS). Water diffusion is used to probe microstructure in neuronal tissues, however, the different water populations and barriers that are present in these tissues turn this into a complex task. Therefore, it is not surprising that recently we have witnessed a burst in the development of new approaches and models that attempt to obtain, non-invasively, detailed microstructural information in the CNS. In this work, we aim at challenging and comparing the microstructural information obtained from single diffusion encoding (SDE) with double diffusion encoding (DDE) MRI. We first applied SDE and DDE MR spectroscopy (MRS) on microcapillary phantoms and then applied SDE and DDE MRI on an ex vivo porcine spinal cord (SC), using similar experimental conditions. The obtained diffusion MRI data were fitted by the same theoretical model, assuming that the signal in every voxel can be approximated as the superposition of a Gaussian-diffusing component and a series of restricted components having infinite cylindrical geometries. The diffusion MRI results were then compared with histological findings. We found a good agreement between the fittings and the experimental data in white matter (WM) voxels of the SC in both diffusion MRI methods. The microstructural information and apparent AADs extracted from SDE MRI were found to be similar or somewhat larger than those extracted from DDE MRI especially when the diffusion time was set to 40 ms. The apparent ADDs extracted from SDE and DDE MRI show reasonable agreement but somewhat weaker correspondence was observed between the diffusion MRI results and histology. The apparent subtle differences between the microstructural information obtained from SDE and DDE MRI are briefly discussed.

White Matter Abnormalities in Multiple Sclerosis Evaluated by Quantitative Synthetic MRI, Diffusion Tensor Imaging, and Neurite Orientation Dispersion and Density Imaging

BACKGROUND AND PURPOSE

A number of MR-derived quantitative metrics have been suggested to assess the pathophysiology of MS, but the reports about combined analyses of these metrics are scarce. Our aim was to assess the spatial distribution of parameters for white matter myelin and axon integrity in patients with relapsing-remitting MS by multiparametric MR imaging.

MATERIALS AND METHODS

Twenty-four patients with relapsing-remitting MS and 24 age- and sex-matched controls were prospectively scanned by quantitative synthetic and 2-shell diffusion MR imaging. Synthetic MR imaging data were used to retrieve relaxometry parameters (R1 and R2 relaxation rates and proton density) and myelin volume fraction. Diffusion tensor metrics (fractional anisotropy and mean, axial, and radial diffusivity) and neurite orientation and dispersion index metrics (intracellular volume fraction, isotropic volume fraction, and orientation dispersion index) were retrieved from diffusion MR imaging data. These data were analyzed using Tract-Based Spatial Statistics.

RESULTS

Patients with MS showed significantly lower fractional anisotropy and myelin volume fraction and higher isotropic volume fraction in widespread white matter areas. Areas with different isotropic volume fractions were included within areas with lower fractional anisotropy. Myelin volume fraction showed no significant difference in some areas with significantly decreased fractional anisotropy in MS, including in the genu of the corpus callosum and bilateral anterior corona radiata, whereas myelin volume fraction was significantly decreased in some areas where fractional anisotropy showed no significant difference, including the bilateral posterior limb of the internal capsule, external capsule, sagittal striatum, fornix, and uncinate fasciculus.

CONCLUSIONS

We found differences in spatial distribution of abnormality in fractional anisotropy, isotropic volume fraction, and myelin volume fraction distribution in MS, which might be useful for characterizing white matter in patients with MS.

Effect of intravoxel incoherent motion on diffusion parameters in normal brain

At very low diffusion weighting the diffusion MRI signal is affected by intravoxel incoherent motion (IVIM) caused by dephasing of magnetization due to incoherent blood flow in capillaries or other sources of microcirculation. While IVIM measurements at low diffusion weightings have been frequently used to investigate perfusion in the body as well as in malignant tissue, the effect and origin of IVIM in normal brain tissue is not completely established. We investigated the IVIM effect on the brain diffusion MRI signal in a cohort of 137 radiologically-normal patients (62 male; mean age = 50.2 ± 17.8, range = 18 to 94). We compared the diffusion tensor parameters estimated from a mono-exponential fit at b = 0 and 1000 s/mm² versus at b = 250 and 1000 s/mm². The asymptotic fitting method allowed for quantitative assessment of the IVIM signal fraction f* in specific brain tissue and regions. Our results show a mean (median) percent difference in the mean diffusivity of about 4.5 (4.9)% in white matter (WM), about 7.8 (8.7)% in cortical gray matter (GM), and 4.3 (4.2)% in thalamus. Corresponding perfusion fraction f* was estimated to be 0.033 (0.032) in WM, 0.066 (0.065) in cortical GM, and 0.033 (0.030) in the thalamus. The effect of f* with respect to age was found to be significant in cortical GM (Pearson correlation ρ ​= ​0.35, p ​= ​3*10^-5) and the thalamus (Pearson correlation ρ = 0.20, p = 0.022) with an average increase in f* of 5.17*10^-4/year and 3.61*10^-4/year, respectively. Significant correlations between f* and age were not observed for WM, and corollary analysis revealed no effect of gender on f*. Possible origins of the IVIM effect in normal brain tissue are discussed.

In vivo magnetic resonance imaging and spectroscopy. Technological advances and opportunities for applications continue to abound

Over the past decades, the field of in vivo magnetic resonance (MR) has built up an impressive repertoire of data acquisition and analysis technologies for anatomical, functional, physiological, and molecular imaging, the description of which requires many book volumes. As such it is impossible for a few authors to have an authoritative overview of the field and for a brief article to be inclusive. We will therefore focus mainly on data acquisition and attempt to give some insight into the principles underlying current advanced methods in the field and the potential for further innovation. In our view, the foreseeable future is expected to show continued rapid progress, for instance in imaging of microscopic tissue properties in vivo, assessment of functional and anatomical connectivity, higher resolution physiologic and metabolic imaging, and even imaging of receptor binding. In addition, acquisition speed and information content will continue to increase due to the continuous development of approaches for parallel imaging (including simultaneous multi-slice imaging), compressed sensing, and MRI fingerprinting. Finally, artificial intelligence approaches are becoming more realistic and will have a tremendous effect on both acquisition and analysis strategies. Together, these developments will continue to provide opportunity for scientific discovery and, in combination with large data sets from other fields such as genomics, allow the ultimate realization of precision medicine in the clinic.

Associations between modifiable risk factors and white matter of the aging brain: insights from diffusion tensor imaging studies

There is increasing interest in factors that may modulate white matter (WM) breakdown and, consequentially, age-related cognitive and behavioral deficits. Recent diffusion tensor imaging studies have examined the relationship of such factors with WM microstructure. This review summarizes the evidence regarding the relationship between WM microstructure and recognized modifiable factors, including hearing loss, hypertension, diabetes, obesity, smoking, depressive symptoms, physical (in) activity, and social isolation, as well as sleep disturbances, diet, cognitive training, and meditation. Current cross-sectional evidence suggests a clear link between loss of WM integrity (lower fractional anisotropy and higher mean diffusivity) and hypertension, obesity, diabetes, and smoking; a relationship that seems to hold for hearing loss, social isolation, depressive symptoms, and sleep disturbances. Physical activity, cognitive training, diet, and meditation, on the other hand, may protect WM with aging. Preliminary evidence from cross-sectional studies of treated risk factors suggests that modification of factors could slow down negative effects on WM microstructure. Careful intervention studies are needed for this literature to contribute to public health initiatives going forward.

Histological validation of per-bundle water diffusion metrics within a region of fiber crossing following axonal degeneration

Micro-architectural characteristics of white matter can be inferred through analysis of diffusion-weighted magnetic resonance imaging (dMRI). The diffusion-dependent signal can be analyzed through several methods, with the tensor model being the most frequently used due to its straightforward interpretation and low requirements for acquisition parameters. While valuable information can be gained from the tensor-derived metrics in regions of homogeneous tissue organization, this model does not provide reliable microstructural information at crossing fiber regions, which are pervasive throughout human white matter. Several multiple fiber models have been proposed that seem to overcome the limitations of the tensor, with few providing per-bundle dMRI-derived metrics. However, biological interpretations of such metrics are limited by the lack of histological confirmation. To this end, we developed a straightforward biological validation framework. Unilateral retinal ischemia was induced in ten rats, which resulted in axonal (Wallerian) degeneration of the corresponding optic nerve, while the contralateral was left intact; the intact and injured axonal populations meet at the optic chiasm as they cross the midline, generating a fiber crossing region in which each population has different diffusion properties. Five rats served as controls. High-resolution ex vivo dMRI was acquired five weeks after experimental procedures. We correlated and compared histology to per-bundle descriptors derived from three methodologies for dMRI analysis (constrained spherical deconvolution and two multi-tensor representations). We found a tight correlation between axonal density (as evaluated through automatic segmentation of histological sections) with per-bundle apparent fiber density and fractional anisotropy (derived from dMRI). The multi-fiber methods explored were able to correctly identify the damaged fiber populations in a region of fiber crossings (chiasm). Our results provide validation of metrics that bring substantial and clinically useful information about white-matter tissue at crossing fiber regions. Our proposed framework is useful to validate other current and future dMRI methods.

Dimensionality reduction of diffusion MRI measures for improved tractometry of the human brain

Various diffusion MRI (dMRI) measures have been proposed for characterising tissue microstructure over the last 15 years. Despite the growing number of experiments using different dMRI measures in assessments of white matter, there has been limited work on: 1) examining their covariance along specific pathways; and on 2) combining these different measures to study tissue microstructure. Indeed, it quickly becomes intractable for existing analysis pipelines to process multiple measurements at each voxel and at each vertex forming a streamline, highlighting the need for new ways to visualise or analyse such high-dimensional data. In a sample of 36 typically developing children aged 8–18 years, we profiled various commonly used dMRI measures across 22 brain pathways. Using a data-reduction approach, we identified two biologically-interpretable components that capture 80% of the variance in these dMRI measures. The first derived component captures properties related to hindrance and restriction in tissue microstructure, while the second component reflects characteristics related to tissue complexity and orientational dispersion. We then demonstrate that the components generated by this approach preserve the biological relevance of the original measurements by showing age-related effects across developmentally sensitive pathways. In summary, our findings demonstrate that dMRI analyses can benefit from dimensionality reduction techniques, to help disentangling the neurobiological underpinnings of white matter organisation.

White matter alterations in adult with autism spectrum disorder evaluated using diffusion kurtosis imaging

PURPOSE

Autism spectrum disorder (ASD) is related to impairment in various white matter (WM) pathways. Utility of the recently developed two-compartment model of diffusion kurtosis imaging (DKI) to analyse axial diffusivity of WM is restricted by several limitations. The present study aims to validate the utility of model-free DKI in the evaluation of WM alterations in ASD and analyse the potential relationship between DKI-evident WM alterations and personality scales.

METHODS

Overall, 15 participants with ASD and 15 neurotypical (NT) controls were scanned on a 3 T magnetic resonance (MR) scanner, and scores for autism quotient (AQ), systemising quotient (SQ) and empathising quotient (EQ) were obtained for both groups. Multishell diffusion-weighted MR data were acquired using two b-values (1000 and 2000 s/mm²). Differences in mean kurtosis (MK), radial kurtosis (RK) and axial kurtosis (AK) between the groups were evaluated using tract-based spatial statistics (TBSS). Finally, the relationships between the kurtosis indices and personality quotients were examined.

RESULTS

The ASD group demonstrated significantly lower AK in the body and splenium of corpus callosum than the NT group; however, no other significant differences were identified. Negative correlations were found between AK and AQ or SQ, predominantly in WM areas related to social-emotional processing such as uncinate fasciculus, inferior fronto-occipital fasciculus, and inferior and superior longitudinal fasciculi.

CONCLUSIONS

Model-free DKI and its indices may represent a novel, objective method for detecting the disease severity and WM alterations in patients with ASD.

Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling

In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.

Nonparenchymal fluid is the source of increased mean diffusivity in preclinical Alzheimer's disease

Introduction

Although increased mean diffusivity of the white matter has been repeatedly linked to Alzheimer’s disease pathology, the underlying mechanism is not known.

Methods

Here, we used ADNI-3 multishell diffusion magnetic resonance imaging data to separate the diffusion signal of the parenchyma from less hindered fluid pools within the white matter such as perivascular space fluid and fluid-filled cavities.

Results

We found that the source of the pathological increase of the mean diffusivity is the increased nonparenchymal fluid, often found in lacunes and perivascular spaces. In this cohort, the cognitive decline was significantly associated with the fluid increase and not with the microstructural changes of the white matter parenchyma itself. The white matter fluid increase was dominantly observed in the sagittal stratum and anterior thalamic radiation.

Discussion

These findings are positive steps toward understanding the pathophysiology of white matter alteration and its role in the cognitive decline.

Brain networks and their relevance for stroke rehabilitation

Stroke has long been regarded as focal disease with circumscribed damage leading to neurological deficits. However, advances in methods for assessing the human brain and in statistics have enabled new tools for the examination of the consequences of stroke on brain structure and function. Thereby, it has become evident that stroke has impact on the entire brain and its network properties and can therefore be considered as a network disease. The present review first gives an overview of current methodological opportunities and pitfalls for assessing stroke-induced changes and reorganization in the human brain. We then summarize principles of plasticity after stroke that have emerged from the assessment of networks. Thereby, it is shown that neurological deficits do not only arise from focal tissue damage but also from local and remote changes in white-matter tracts and in neural interactions among wide-spread networks. Similarly, plasticity and clinical improvements are associated with specific compensatory structural and functional patterns of neural network interactions. Innovative treatment approaches have started to target such network patterns to enhance recovery. Network assessments to predict treatment response and to individualize rehabilitation is a promising way to enhance specific treatment effects and overall outcome after stroke.

Towards a 'resolution limit' for DW-MRI tumor microstructural models: A simulation study investigating the feasibility of distinguishing between microstructural changes

PURPOSE

To determine the feasibility of extracting sufficiently precise estimates of cell radius, R, and intracellular volume fraction, fi , from DW-MRI data in order to distinguish between specific microstructural changes tissue may undergo, specifically focusing on cell death in tumors.

METHODS

Simulations with optimized and non-optimized clinical acquisitions were performed for a range of microstructures, using a two-compartment model. The ability to distinguish between (i) cell shrinkage with cell density constant, mimicking apoptosis, and (ii) cell size constant with cell density decreasing, mimicking loss of cells, was evaluated based on the precision of simulated parameter estimates. Relationships between parameter precision, SNR, and the magnitude of specific parameter changes, were used to infer SNR requirements for detecting changes.

RESULTS

Accuracy and precision depended on microstructural properties, SNR, and the acquisition protocol. The main benefit of optimized acquisitions tended to be improved accuracy and precision of R, particularly for small cells. In most cases considered, higher SNR was required for detecting changes in R than for changes in fi .

CONCLUSIONS

Given the relative changes in R and fi due to apoptosis, simulations indicate that, for a range of microstructures, detecting changes in R require higher SNR than detecting changes in fi , and that such SNR is typically not achieved in clinical data. This suggests that if apoptotic cell size decreases are to be detected in clinical settings, improved SNR is required. Comparing measurement precision with the magnitude of expected biological changes should form part of the validation process for potential biomarkers.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

Modelling white matter with spherical deconvolution: How and why?

Since the realization that diffusion MRI can probe the microstructural organization and orientation of biological tissue in vivo and non-invasively, a multitude of diffusion imaging methods have been developed and applied to study the living human brain. Diffusion tensor imaging was the first model to be widely adopted in clinical and neuroscience research, but it was also clear from the beginning that it suffered from limitations when mapping complex configurations, such as crossing fibres. In this review, we highlight the main steps that have led the field of diffusion imaging to move from the tensor model to the adoption of diffusion and fibre orientation density functions as a more effective way to describe the complexity of white matter organization within each brain voxel. Among several techniques, spherical deconvolution has emerged today as one of the main approaches to model multiple fibre orientations and for tractography applications. Here we illustrate the main concepts and the reasoning behind this technique, as well as the latest developments in the field. The final part of this review provides practical guidelines and recommendations on how to set up processing and acquisition protocols suitable for spherical deconvolution.

Modeling conduction delays in the corpus callosum using MRI-measured g-ratio

Conduction of action potentials along myelinated axons is affected by their structural features, such as the axonal g-ratio, the ratio between the inner and outer diameters of the myelin sheath surrounding the axon. The effect of g-ratio variance on conduction properties has been quantitatively evaluated using single-axon models. It has recently become possible to estimate a g-ratio weighted measurement in vivo using quantitative MRI. Nevertheless, it is still unclear whether the variance in the g-ratio in the healthy human brain leads to significant differences in conduction velocity. In this work we tested whether the g-ratio MRI measurement can be used to predict conduction delays in the corpus callosum. We present a comprehensive framework in which the structural properties of fibers (i.e. length and g-ratio, measured using MRI), are incorporated in a biophysical model of axon conduction, to model conduction delays of long-range white matter fibers. We applied this framework to the corpus callosum, and found conduction delay estimates that are compatible with previously estimated values of conduction delays. We account for the variance in the velocity given the axon diameter distribution in the splenium, mid-body and genu, to further compare the fibers within the corpus callosum. Conduction delays have been suggested to increase with age. Therefore, we investigated whether there are differences in the g-ratio and the fiber length between young and old adults, and whether this leads to a difference in conduction speed and delays. We found very small differences between the predicted delays of the two groups in the motor fibers of the corpus callosum. We also found that the motor fibers of the corpus callosum have the fastest conduction estimates. Using the axon diameter distributions, we found that the occipital fibers have the slowest estimations, while the frontal and motor fiber tracts have similar estimates. Our study provides a framework for predicting conduction latencies in vivo. The framework could have major implications for future studies of white matter diseases and large range network computations. Our results highlight the need for improving additional in vivo measurements of white matter microstructure.

The Superoanterior Fasciculus (SAF): A Novel White Matter Pathway in the Human Brain?

Fiber tractography (FT) using diffusion magnetic resonance imaging (dMRI) is widely used for investigating microstructural properties of white matter (WM) fiber-bundles and for mapping structural connections of the human brain. While studying the architectural configuration of the brain's circuitry with FT is not without controversy, recent progress in acquisition, processing, modeling, analysis, and visualization of dMRI data pushes forward the reliability in reconstructing WM pathways. Despite being aware of the well-known pitfalls in analyzing dMRI data and several other limitations of FT discussed in recent literature, we present the superoanterior fasciculus (SAF), a novel bilateral fiber tract in the frontal region of the human brain that-to the best of our knowledge-has not been documented. The SAF has a similar shape to the anterior part of the cingulum bundle, but it is located more frontally. To minimize the possibility that these FT findings are based on acquisition or processing artifacts, different dMRI data sets and processing pipelines have been used to describe the SAF. Furthermore, we evaluated the configuration of the SAF with complementary methods, such as polarized light imaging (PLI) and human brain dissections. The FT results of the SAF demonstrate a long pathway, consistent across individuals, while the human dissections indicate fiber pathways connecting the postero-dorsal with the antero-dorsal cortices of the frontal lobe.

Survey of Image Processing Techniques for Brain Pathology Diagnosis: Challenges and Opportunities

In recent years, a number of new products introduced to the global market combine intelligent robotics, artificial intelligence and smart interfaces to provide powerful tools to support professional decision making. However, while brain disease diagnosis from the brain scan images is supported by imaging robotics, the data analysis to form a medical diagnosis is performed solely by highly trained medical professionals. Recent advances in medical imaging techniques, artificial intelligence, machine learning and computer vision present new opportunities to build intelligent decision support tools to aid the diagnostic process, increase the disease detection accuracy, reduce error, automate the monitoring of patient's recovery, and discover new knowledge about the disease cause and its treatment. This article introduces the topic of medical diagnosis of brain diseases from the MRI based images. We describe existing, multi-modal imaging techniques of the brain's soft tissue and describe in detail how are the resulting images are analyzed by a radiologist to form a diagnosis. Several comparisons between the best results of classifying natural scenes and medical image analysis illustrate the challenges of applying existing image processing techniques to the medical image analysis domain. The survey of medical image processing methods also identified several knowledge gaps, the need for automation of image processing analysis, and the identification of the brain structures in the medical images that differentiate healthy tissue from a pathology. This survey is grounded in the cases of brain tumor analysis and the traumatic brain injury diagnoses, as these two case studies illustrate the vastly different approaches needed to define, extract, and synthesize meaningful information from multiple MRI image sets for a diagnosis. Finally, the article summarizes artificial intelligence frameworks that are built as multi-stage, hybrid, hierarchical information processing work-flows and the benefits of applying these models for medical diagnosis to build intelligent physician's aids with knowledge transparency, expert knowledge embedding, and increased analytical quality.

Modeling white matter microstructure with fiber ball imaging

Fiber ball imaging (FBI) provides a means of calculating the fiber orientation density function (fODF) in white matter from diffusion MRI (dMRI) data obtained over a spherical shell with a b-value of about 4000 s/mm2 or higher. By supplementing this FBI-derived fODF with dMRI data acquired for two lower b-value shells, it is shown that several microstructural parameters may be estimated, including the axonal water fraction (AWF) and the intrinsic intra-axonal diffusivity. This fiber ball white matter (FBWM) modeling method is demonstrated for dMRI data acquired from healthy volunteers, and the results are compared with those of the white matter tract integrity (WMTI) method. Both the AWF and the intra-axonal diffusivity obtained with FBWM are found to be significantly larger than for WMTI, with the FBWM values for the intra-axonal diffusivity being more consistent with recent results obtained using isotropic diffusion weighting. An important practical advantage of FBWM is that the only nonlinear fitting required is the minimization of a cost function with just a single free parameter, which facilitates the implementation of efficient and robust numerical routines.

Diffusional kurtosis imaging and white matter microstructure modeling in a clinical study of major depressive disorder

Major depressive disorder (MDD) is a globally prevalent psychiatric disorder that results from disruption of multiple neural circuits involved in emotional regulation. Although previous studies using diffusion tensor imaging (DTI) found smaller values of fractional anisotropy (FA) in the white matter, predominantly in the frontal lobe, of patients with MDD, studies using diffusion kurtosis imaging (DKI) are scarce. Here, we used DKI whole-brain analysis with tract-based spatial statistics (TBSS) to investigate the brain microstructural abnormalities in MDD. Twenty-six patients with MDD and 42 age- and sex-matched control subjects were enrolled. To investigate the microstructural pathology underlying the observations in DKI, a compartment model analysis was conducted focusing on the corpus callosum. In TBSS, the patients with MDD showed significantly smaller values of FA in the genu and frontal portion of the body of the corpus callosum. The patients also had smaller values of mean kurtosis (MK) and radial kurtosis (RK), but MK and RK abnormalities were distributed more widely compared with FA, predominantly in the frontal lobe but also in the parietal, occipital, and temporal lobes. Within the callosum, the regions with smaller MK and RK were located more posteriorly than the region with smaller FA. Model analysis suggested significantly smaller values of intra-neurite signal fraction in the body of the callosum and greater fiber dispersion in the genu, which were compatible with the existing literature of white matter pathology in MDD. Our results show that DKI is capable of demonstrating microstructural alterations in the brains of patients with MDD that cannot be fully depicted by conventional DTI. Though the issues of model validation and parameter estimation still remain, it is suggested that diffusion MRI combined with a biophysical model is a promising approach for investigation of the pathophysiology of MDD.

Characterization of salivary gland tumours with diffusion tensor imaging

OBJECTIVES

To characterize salivary glands tumours with diffusion tensor imaging.

METHODS

This study was conducted upon 53 patients (aged 18-81 years: mean 37 years) with salivary gland tumours that underwent diffusion tensor imaging was obtained using a single-shot echoplanar imaging sequence with parallel imaging at 1.5 T scanner. 48 slices were obtained, with a thickness of 2.5 mm, with no gap and the total scan duration was 7-8 min. The fractional anisotropy (FA) and the mean diffusivity (MD) value of the salivary gland tumours was calculated and correlated with pathological findings. Image analysis was performed by one radiologist. The receiver operating characteristic curve was drawn to detect the cut-off point of FA and MD used to characterize salivary gland tumours.

RESULTS

The mean FA and MD of malignant salivary gland tumours (n = 17) (0.41 ± 0.07 and 0.89 ± 0.15 × 10-3 mm2 s-1) was significantly different (p = 0.001) than that of benign tumours (n = 36) (0.19 ± 0.07 and 1.28 ± 0.42 × 10-3 mm2 s-1), respectively. Combined FA and MD used to differentiate malignant from benign tumours has an area under the curve (AUC) of 0.974, and an accuracy of 86%. There was a significant difference in FA between Warthin tumours and malignant tumours (p = 0.001). Selection FA of 0.35 to differentiate malignant tumours from Warthin tumours revealed AUC of 0.878 and an accuracy of 80%. There was a significant difference in FA and MD of malignant tumours and pleomorphic adenomas (p = 0.001). Combined FA and MD used to differentiate malignant tumours from pleomorphic adenomas revealed AUC of 0.993, and an accuracy of 93%. There was a significant difference in FA and MD of Warthin tumours and pleomorphic adenomas (p = 0.001). Combined FA and MD used to differentiate Warthin tumours from pleomorphic adenomas revealed AUC of 0.978, and an accuracy of 86%.

CONCLUSIONS

Diffusion-weighed imaging is a promising non-invasive method and it may be useful for the characterization and differentiation of benign and malignant salivary gland tumours.

The impact of in-scanner head motion on structural connectivity derived from diffusion MRI

Multiple studies have shown that data quality is a critical confound in the construction of brain networks derived from functional MRI. This problem is particularly relevant for studies of human brain development where important variables (such as participant age) are correlated with data quality. Nevertheless, the impact of head motion on estimates of structural connectivity derived from diffusion tractography methods remains poorly characterized. Here, we evaluated the impact of in-scanner head motion on structural connectivity using a sample of 949 participants (ages 8-23 years old) who passed a rigorous quality assessment protocol for diffusion magnetic resonance imaging (dMRI) acquired as part of the Philadelphia Neurodevelopmental Cohort. Structural brain networks were constructed for each participant using both deterministic and probabilistic tractography. We hypothesized that subtle variation in head motion would systematically bias estimates of structural connectivity and confound developmental inference, as observed in previous studies of functional connectivity. Even following quality assurance and retrospective correction for head motion, eddy currents, and field distortions, in-scanner head motion significantly impacted the strength of structural connectivity in a consistency- and length-dependent manner. Specifically, increased head motion was associated with reduced estimates of structural connectivity for network edges with high inter-subject consistency, which included both short- and long-range connections. In contrast, motion inflated estimates of structural connectivity for low-consistency network edges that were primarily shorter-range. Finally, we demonstrate that age-related differences in head motion can both inflate and obscure developmental inferences on structural connectivity. Taken together, these data delineate the systematic impact of head motion on structural connectivity, and provide a critical context for identifying motion-related confounds in studies of structural brain network development.

On modeling

Mapping tissue microstructure with MRI holds great promise as a noninvasive window into tissue organization at the cellular level. Having originated within the realm of diffusion NMR in the late 1970s, this field is experiencing an exponential growth in the number of publications. At the same time, model-based approaches are also increasingly incorporated into advanced MRI acquisition and reconstruction techniques. However, after about two decades of intellectual and financial investment, microstructural mapping has yet to find a single commonly accepted clinical application. Here, we suggest that slow progress in clinical translation may signify unresolved fundamental problems. We outline such problems and related practical pitfalls, as well as review strategies for developing and validating tissue microstructure models, to provoke a discussion on how to bridge the gap between our scientific aspirations and the clinical reality. We argue for recalibrating the efforts of our community toward a more systematic focus on fundamental research aimed at identifying relevant degrees of freedom affecting the measured MR signal. Such a focus is essential for realizing the truly revolutionary potential of noninvasive three-dimensional in vivo microstructural mapping.