An analytical model of restricted diffusion in bovine optic nerve

AxonPacking: An Open-Source Software to Simulate Arrangements of Axons in White Matter

HIGHLIGHTS

AxonPacking: Open-source software for simulating white matter microstructure.

Validation on a theoretical disk packing problem.

Reproducible and stable for various densities and diameter distributions.

Can be used to study interplay between myelin/fiber density and restricted fraction.

Quantitative Magnetic Resonance Imaging (MRI) can provide parameters that describe white matter microstructure, such as the fiber volume fraction (FVF), the myelin volume fraction (MVF) or the axon volume fraction (AVF) via the fraction of restricted water (fr). While already being used for clinical application, the complex interplay between these parameters requires thorough validation via simulations. These simulations required a realistic, controlled and adaptable model of the white matter axons with the surrounding myelin sheath. While there already exist useful algorithms to perform this task, none of them combine optimisation of axon packing, presence of myelin sheath and availability as free and open source software. Here, we introduce a novel disk packing algorithm that addresses these issues. The performance of the algorithm is tested in term of reproducibility over 50 runs, resulting density, and stability over iterations. This tool was then used to derive multiple values of FVF and to study the impact of this parameter on fr and MVF in light of the known microstructure based on histology sample. The standard deviation of the axon density over runs was lower than 10⁻³ and the expected hexagonal packing for monodisperse disks was obtained with a density close to the optimal density (obtained: 0.892, theoretical: 0.907). Using an FVF ranging within [0.58, 0.82] and a mean inter-axon gap ranging within [0.1, 1.1] μm, MVF ranged within [0.32, 0.44] and fr ranged within [0.39, 0.71], which is consistent with the histology. The proposed algorithm is implemented in the open-source software AxonPacking (https://github.com/neuropoly/axonpacking) and can be useful for validating diffusion models as well as for enabling researchers to study the interplay between microstructure parameters when evaluating qMRI methods.

Microstructure Imaging of Crossing (MIX) White Matter Fibers from diffusion MRI

Diffusion MRI (dMRI) reveals microstructural features of the brain white matter by quantifying the anisotropic diffusion of water molecules within axonal bundles. Yet, identifying features such as axonal orientation dispersion, density, diameter, etc., in complex white matter fiber configurations (e.g. crossings) has proved challenging. Besides optimized data acquisition and advanced biophysical models, computational procedures to fit such models to the data are critical. However, these procedures have been largely overlooked by the dMRI microstructure community and new, more versatile, approaches are needed to solve complex biophysical model fitting problems. Existing methods are limited to models assuming single fiber orientation, relevant to limited brain areas like the corpus callosum, or multiple orientations but without the ability to extract detailed microstructural features. Here, we introduce a new and versatile optimization technique (MIX), which enables microstructure imaging of crossing white matter fibers. We provide a MATLAB implementation of MIX, and demonstrate its applicability to general microstructure models in fiber crossings using synthetic as well as ex-vivo and in-vivo brain data.

Characterization of spinal cord white matter by suppressing signal from hindered space. A Monte Carlo simulation and an ex vivo ultrahigh-b diffusion-weighted imaging study

Signal measured from white matter in diffusion-weighted imaging is difficult to interpret because of the heterogeneous structure of white matter. Characterization of the white matter will be straightforward if the signal contributed from the hindered space is suppressed and purely restricted signal is analyzed. In this study, a Monte Carlo simulation (MCS) of water diffusion in white matter was performed to understand the behavior of the diffusion-weighted signal in white matter. The signal originating from the hindered space of an excised pig cervical spinal cord white matter was suppressed using the ultrahigh-b radial diffusion-weighted imaging. A light microscopy image of a section of white matter was obtained from the excised pig cervical spinal cord for the MCS. The radial diffusion-weighted signals originating from each of the intra-axonal, extra-axonal, and total spaces were studied using the MCS. The MCS predicted that the radial diffusion-weighted signal remains almost constant in the intra-axonal space and decreases gradually to about 2% of its initial value in the extra-axonal space when the b-value is increased to 30,000s/mm². The MCS also revealed that the diffusion-weighted signal for a b-value greater than 20,000s/mm² is mostly from the intra-axonal space. The decaying behavior of the signal-b curve obtained from ultrahigh-b diffusion-weighted imaging (b_max∼30,000s/mm²) of the excised pig cord was very similar to the decaying behavior of the total signal-b curve synthesized in the MCS. A mono-exponential plus constant fitting of the signal-b curve obtained from a white matter pixel estimated the values of constant fraction and apparent diffusion coefficient of decaying fraction as 0.32±0.05 and (0.16±0.01)×10^-3mm²/s, respectively, which agreed well with the results of the MCS. The signal measured in the ultrahigh-b region (b>20,000s/mm²) is mostly from the restricted (intra-axonal) space. Integrity and intactness of the axons can be evaluated by assessing the remaining signal in the ultrahigh-b region.

Investigating white matter fibre density and morphology using fixel-based analysis

Voxel-based analysis of diffusion MRI data is increasingly popular. However, most white matter voxels contain contributions from multiple fibre populations (often referred to as crossing fibres), and therefore voxel-averaged quantitative measures (e.g. fractional anisotropy) are not fibre-specific and have poor interpretability. Using higher-order diffusion models, parameters related to fibre density can be extracted for individual fibre populations within each voxel ('fixels'), and recent advances in statistics enable the multi-subject analysis of such data. However, investigating within-voxel microscopic fibre density alone does not account for macroscopic differences in the white matter morphology (e.g. the calibre of a fibre bundle). In this work, we introduce a novel method to investigate the latter, which we call fixel-based morphometry (FBM). To obtain a more complete measure related to the total number of white matter axons, information from both within-voxel microscopic fibre density and macroscopic morphology must be combined. We therefore present the FBM method as an integral piece within a comprehensive fixel-based analysis framework to investigate measures of fibre density, fibre-bundle morphology (cross-section), and a combined measure of fibre density and cross-section. We performed simulations to demonstrate the proposed measures using various transformations of a numerical fibre bundle phantom. Finally, we provide an example of such an analysis by comparing a clinical patient group to a healthy control group, which demonstrates that all three measures provide distinct and complementary information. By capturing information from both sources, the combined fibre density and cross-section measure is likely to be more sensitive to certain pathologies and more directly interpretable.

Demonstration of Non-Gaussian Restricted Diffusion in Tumor Cells Using Diffusion Time-Dependent Diffusion-Weighted Magnetic Resonance Imaging Contrast

The diffusion-weighted magnetic resonance imaging (DWI) technique enables quantification of water mobility for probing microstructural properties of biological tissue and has become an effective tool for collecting information about the underlying pathology of cancerous tissue. Measurements using multiple b-values have indicated biexponential signal attenuation, ascribed to “fast” (high ADC) and “slow” (low ADC) diffusion components. In this empirical study, we investigate the properties of the diffusion time (Δ)-dependent components of the diffusion-weighted (DW) signal in a constant b-value experiment. A xenograft gliobastoma mouse was imaged using Δ = 11 ms, 20 ms, 40 ms, 60 ms, and b = 500–4000 s/mm² in intervals of 500 s/mm². Data were corrected for EPI distortions, and the Δ-dependence on the DW-signal was measured within three regions of interest [intermediate- and high-density tumor regions and normal-appearing brain (NAB) tissue regions]. In this study, we verify the assumption that the slow decaying component of the DW-signal is non-Gaussian and dependent on Δ, consistent with restricted diffusion of the intracellular space. As the DW-signal is a function of Δ and is specific to restricted diffusion, manipulating Δ at constant b-value (cb) provides a complementary and direct approach for separating the restricted from the hindered diffusion component. We found that Δ-dependence is specific to the tumor tissue signal. Based on an extended biexponential model, we verified the interpretation of the diffusion time-dependent contrast and successfully estimated the intracellular restricted ADC, signal volume fraction, and cell size within each ROI.

Assessing the accuracy of using oscillating gradient spin echo sequences with AxCaliber to infer micron-sized axon diameters

OBJECTIVE

Current magnetic resonance imaging (MRI) axon diameter measurements rely on the pulsed gradient spin-echo sequence, which is unable to provide diffusion times short enough to measure small axon diameters. This study combines the AxCaliber axon diameter fitting method with data generated from Monte Carlo simulations of oscillating gradient spin-echo sequences (OGSE) to infer micron-sized axon diameters, in order to determine the feasibility of using MRI to infer smaller axon diameters in brain tissue.

MATERIALS AND METHODS

Monte Carlo computer simulation data were synthesized from tissue geometries of cylinders of different diameters using a range of gradient frequencies in the cosine OGSE sequence . Data were fitted to the AxCaliber method modified to allow the new pulse sequence. Intra- and extra-axonal water were studied separately and together.

RESULTS

The simulations revealed the extra-axonal model to be problematic. Rather than change the model, we found that restricting the range of gradient frequencies such that the measured apparent diffusion coefficient was constant over that range resulted in more accurate fitted diameters. Thus a careful selection of frequency ranges is needed for the AxCaliber method to correctly model extra-axonal water, or adaptations to the method are needed. This restriction helped reduce the necessary gradient strengths for measurements that could be performed with parameters feasible for a Bruker BG6 gradient set. For these experiments, the simulations inferred diameters as small as 0.5 μm on square-packed and randomly packed cylinders. The accuracy of the inferred diameters was found to be dependent on the signal-to-noise ratio (SNR), with smaller diameters more affected by noise, although all diameter distributions were distinguishable from one another for all SNRs tested.

CONCLUSION

The results of this study indicate the feasibility of using MRI with OGSE on preclinical scanners to infer small axon diameters.

Fusion in diffusion MRI for improved fibre orientation estimation: An application to the 3T and 7T data of the Human Connectome Project

Determining the acquisition parameters in diffusion magnetic resonance imaging (dMRI) is governed by a series of trade-offs. Images of lower resolution have less spatial specificity but higher signal to noise ratio (SNR). At the same time higher angular contrast, important for resolving complex fibre patterns, also yields lower SNR. Considering these trade-offs, the Human Connectome Project (HCP) acquires high quality dMRI data for the same subjects at different field strengths (3T and 7T), which are publically released. Due to differences in the signal behavior and in the underlying scanner hardware, the HCP 3T and 7T data have complementary features in k- and q-space. The 3T dMRI has higher angular contrast and resolution, while the 7T dMRI has higher spatial resolution. Given the availability of these datasets, we explore the idea of fusing them together with the aim of combining their benefits. We extend a previously proposed data-fusion framework and apply it to integrate both datasets from the same subject into a single joint analysis. We use a generative model for performing parametric spherical deconvolution and estimate fibre orientations by simultaneously using data acquired under different protocols. We illustrate unique features from each dataset and how they are retained after fusion. We further show that this allows us to complement benefits and improve brain connectivity analysis compared to analyzing each of the datasets individually.

Simulations on the influence of myelin water in diffusion-weighted imaging

While myelinated axons present an important barrier to water diffusion, many models used to interpret DWI signal neglect other potential influences of myelin. In this work, Monte Carlo simulations were used to test the sensitivity of DWI results to the diffusive properties of water within myelin. Within these simulations, the apparent diffusion coefficient (D app) varied slowly over several orders of magnitude of the coefficient of myelin water diffusion (D m), but exhibited important differences compared to D app values simulated that neglect D m (=0). Compared to D app, the apparent diffusion kurtosis (K app) was generally more sensitive to D m. Simulations also tested the sensitivity of D app and K app to the amount of myelin present. Unique variations in D app and K app caused by differences in the myelin volume fraction were diminished when myelin water diffusion was included. Also, expected trends in D app and K app with experimental echo time were reduced or inverted when accounting for myelin water diffusion, and these reduced/inverted trends were seen experimentally in ex vivo rat brain DWI experiments. In general, myelin water has the potential to subtly influence DWI results and bias models of DWI that neglect these components of white matter.

Multi-compartment microscopic diffusion imaging

This paper introduces a multi-compartment model for microscopic diffusion anisotropy imaging. The aim is to estimate microscopic features specific to the intra- and extra-neurite compartments in nervous tissue unconfounded by the effects of fibre crossings and orientation dispersion, which are ubiquitous in the brain. The proposed MRI method is based on the Spherical Mean Technique (SMT), which factors out the neurite orientation distribution and thus provides direct estimates of the microscopic tissue structure. This technique can be immediately used in the clinic for the assessment of various neurological conditions, as it requires only a widely available off-the-shelf sequence with two b-shells and high-angular gradient resolution achievable within clinically feasible scan times. To demonstrate the developed method, we use high-quality diffusion data acquired with a bespoke scanner system from the Human Connectome Project. This study establishes the normative values of the new biomarkers for a large cohort of healthy young adults, which may then support clinical diagnostics in patients. Moreover, we show that the microscopic diffusion indices offer direct sensitivity to pathological tissue alterations, exemplified in a preclinical animal model of Tuberous Sclerosis Complex (TSC), a genetic multi-organ disorder which impacts brain microstructure and hence may lead to neurological manifestations such as autism, epilepsy and developmental delay.

Parametric Probability Distribution Functions for Axon Diameters of Corpus Callosum

Axon diameter is an important neuroanatomical characteristic of the nervous system that alters in the course of neurological disorders such as multiple sclerosis. Axon diameters vary, even within a fiber bundle, and are not normally distributed. An accurate distribution function is therefore beneficial, either to describe axon diameters that are obtained from a direct measurement technique (e.g., microscopy), or to infer them indirectly (e.g., using diffusion-weighted MRI). The gamma distribution is a common choice for this purpose (particularly for the inferential approach) because it resembles the distribution profile of measured axon diameters which has been consistently shown to be non-negative and right-skewed. In this study we compared a wide range of parametric probability distribution functions against empirical data obtained from electron microscopy images. We observed that the gamma distribution fails to accurately describe the main characteristics of the axon diameter distribution, such as location and scale of the mode and the profile of distribution tails. We also found that the generalized extreme value distribution consistently fitted the measured distribution better than other distribution functions. This suggests that there may be distinct subpopulations of axons in the corpus callosum, each with their own distribution profiles. In addition, we observed that several other distributions outperformed the gamma distribution, yet had the same number of unknown parameters; these were the inverse Gaussian, log normal, log logistic and Birnbaum-Saunders distributions.

White matter microstructure from nonparametric axon diameter distribution mapping

We report the development of a double diffusion encoding (DDE) MRI method to estimate and map the axon diameter distribution (ADD) within an imaging volume. A variety of biological processes, ranging from development to disease and trauma, may lead to changes in the ADD in the central and peripheral nervous systems. Unlike previously proposed methods, this ADD experimental design and estimation framework employs a more general, nonparametric approach, without a priori assumptions about the underlying form of the ADD, making it suitable to analyze abnormal tissue. In the current study, this framework was used on an ex vivo ferret spinal cord, while emphasizing the way in which the ADD can be weighted by either the number or the volume of the axons. The different weightings, which result in different spatial contrasts, were considered throughout this work. DDE data were analyzed to derive spatially resolved maps of average axon diameter, ADD variance, and extra-axonal volume fraction, along with a novel sub-micron restricted structures map. The morphological information contained in these maps was then used to segment white matter into distinct domains by using a proposed k-means clustering algorithm with spatial contiguity and left-right symmetry constraints, resulting in identifiable white matter tracks. The method was validated by comparing histological measures to the estimated ADDs using a quantitative similarity metric, resulting in good agreement. With further acquisition acceleration and experimental parameters adjustments, this ADD estimation framework could be first used preclinically, and eventually clinically, enabling a wide range of neuroimaging applications for improved understanding of neurodegenerative pathologies and assessing microstructural changes resulting from trauma.

Clarifying Human White Matter

Progress in magnetic resonance imaging (MRI) now makes it possible to identify the major white matter tracts in the living human brain. These tracts are important because they carry many of the signals communicated between different brain regions. MRI methods coupled with biophysical modeling can measure the tissue properties and structural features of the tracts that impact our ability to think, feel, and perceive. This review describes the fundamental ideas of the MRI methods used to identify the major white matter tracts in the living human brain.

Towards higher sensitivity and stability of axon diameter estimation with diffusion-weighted MRI

Diffusion-weighted MRI is an important tool for in vivo and non-invasive axon morphometry. The ActiveAx technique utilises an optimised acquisition protocol to infer orientationally invariant indices of axon diameter and density by fitting a model of white matter to the acquired data. In this study, we investigated the factors that influence the sensitivity to small-diameter axons, namely the gradient strength of the acquisition protocol and the model fitting routine. Diffusion-weighted ex. vivo images of the mouse brain were acquired using 16.4-T MRI with high (Gmax of 300 mT/m) and ultra-high (Gmax of 1350 mT/m) gradient strength acquisitions. The estimated axon diameter indices of the mid-sagittal corpus callosum were validated using electron microscopy. In addition, a dictionary-based fitting routine was employed and evaluated. Axon diameter indices were closer to electron microscopy measures when higher gradient strengths were employed. Despite the improvement, estimated axon diameter indices (a lower bound of ~ 1.8 μm) remained higher than the measurements obtained using electron microscopy (~1.2 μm). We further observed that limitations of pulsed gradient spin echo (PGSE) acquisition sequences and axonal dispersion could also influence the sensitivity with which axon diameter indices could be estimated. Our results highlight the influence of acquisition protocol, tissue model and model fitting, in addition to gradient strength, on advanced microstructural diffusion-weighted imaging techniques. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Two-Dimensional Correlation of Isotropic and Directional Diffusion Using NMR

Diffusion nuclear magnetic resonance (NMR) is a powerful technique for studying porous media, but yields ambiguous results when the sample comprises multiple regions with different pore sizes, shapes, and orientations. Inspired by solid-state NMR techniques for correlating isotropic and anisotropic chemical shifts, we propose a diffusion NMR method to resolve said ambiguity. Numerical data inversion relies on sparse representation of the data in a basis of radial and axial diffusivities. Experiments are performed on a composite sample with a cell suspension and a liquid crystal.

Biexponential diffusion alterations in the normal-appearing white matter of glioma patients might indicate the presence of global vasogenic edema

PURPOSE

To investigate normal-appearing white matter (NAWM) microstructure of glioma patients with biexponential diffusion analysis in order to reveal the nature of diffusion abnormalities and to assess whether they are region-specific or global.

MATERIALS AND METHODS

Twenty-four newly diagnosed glioma patients (grade II-IV) and 24 matched control subjects underwent diffusion-weighted imaging at 3T. Diffusion parameters were calculated using monoexponential and biexponential models. Apparent diffusion coefficient (ADC) values were measured in the entire NAWM of the hemisphere contralateral and ipsilateral to the tumor. In the contralateral NAWM, regional ADC values were assessed in the frontal, parietal, occipital, and temporal NAWM.

RESULTS

ADCmono and ADCfast were significantly higher than control values in all investigated regions except the temporal NAWM (P < 0.04). ADCslow was significantly increased in the total contralateral, frontal, and parietal NAWM (P < 0.03), while pslow was decreased in both total hemispheric NAWM and the parietal NAWM of glioma patients compared to controls (P < 0.04). ADCmono , ADCfast , ADCslow , and pslow were significantly different among the NAWM of the four lobes of the contralateral hemisphere in both groups (P < 0.0001), and these regional differences were similar in patients and controls (P > 0.05). Hemispheric ADCmono and pslow differences were different between groups (P < 0.05).

CONCLUSION

Globally altered diffusion parameters suggest the presence of global vasogenic edema in the NAWM of glioma patients, which is further supported by the finding that regional differences in patients follow those found in controls. Alternatively, some tumor infiltration might contribute to diffusion abnormalities in the NAWM, especially in the tumor-affected hemisphere. J. Magn. Reson. Imaging 2016;44:633-641.

Modern Methods for Interrogating the Human Connectome

OBJECTIVES

Connectionist theories of brain function took hold with the seminal contributions of Norman Geschwind a half century ago. Modern neuroimaging techniques have expanded the scientific interest in the study of brain connectivity to include the intact as well as disordered brain.

METHODS

In this review, we describe the most common techniques used to measure functional and structural connectivity, including resting state functional MRI, diffusion MRI, and electroencephalography and magnetoencephalography coherence. We also review the most common analytical approaches used for examining brain interconnectivity associated with these various imaging methods.

RESULTS

This review presents a critical analysis of the assumptions, as well as methodological limitations, of each imaging and analysis approach.

CONCLUSIONS

The overall goal of this review is to provide the reader with an introduction to evaluating the scientific methods underlying investigations that probe the human connectome.

Heritability of complex white matter diffusion traits assessed in a population isolate

INTRODUCTION

Diffusion weighted imaging (DWI) methods can noninvasively ascertain cerebral microstructure by examining pattern and directions of water diffusion in the brain. We calculated heritability for DWI parameters in cerebral white (WM) and gray matter (GM) to study the genetic contribution to the diffusion signals across tissue boundaries.

METHODS

Using Old Order Amish (OOA) population isolate with large family pedigrees and high environmental homogeneity, we compared the heritability of measures derived from three representative DWI methods targeting the corpus callosum WM and cingulate gyrus GM: diffusion tensor imaging (DTI), the permeability-diffusivity (PD) model, and the neurite orientation dispersion and density imaging (NODDI) model. These successively more complex models represent the diffusion signal modeling using one, two, and three diffusion compartments, respectively.

RESULTS

We replicated the high heritability of the DTI-based fractional anisotropy (h(2)  = 0.67) and radial diffusivity (h(2)  = 0.72) in WM. High heritability in both WM and GM tissues were observed for the permeability-diffusivity index from the PD model (h(2)  = 0.64 and 0.84), and the neurite density from the NODDI model (h(2)  = 0.70 and 0.55). The orientation dispersion index from the NODDI model was only significantly heritable in GM (h(2)  = 0.68).

CONCLUSION

DWI measures from multicompartmental models were significantly heritable in WM and GM. DWI can offer valuable phenotypes for genetic research; and genes thus identified may reveal mechanisms contributing to mental and neurological disorders in which diffusion imaging anomalies are consistently found. Hum Brain Mapp 37:525-535, 2016. © 2015 Wiley Periodicals, Inc.

In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter

The presence of micrometer-level restrictions leads to a decrease of diffusion coefficient with diffusion time. Here we investigate this effect in human white matter in vivo. We focus on a broad range of diffusion times, up to 600 ms, covering diffusion length scales up to about 30 μm. We perform stimulated echo diffusion tensor imaging on 5 healthy volunteers and observe a relatively weak time-dependence in diffusion transverse to major fiber tracts. Remarkably, we also find notable time-dependence in the longitudinal direction. Comparing models of diffusion in ordered, confined and disordered media, we argue that the time-dependence in both directions can arise due to structural disorder, such as axonal beads in the longitudinal direction, and the random packing geometry of fibers within a bundle in the transverse direction. These time-dependent effects extend beyond a simple picture of Gaussian compartments, and may lead to novel markers that are specific to neuronal fiber geometry at the micrometer scale.

Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter.

Fiber Orientation and Compartment Parameter Estimation From Multi-Shell Diffusion Imaging

Diffusion MRI offers the unique opportunity of assessing the structural connections of human brains in vivo. With the advance of diffusion MRI technology, multi-shell imaging methods are becoming increasingly practical for large scale studies and clinical application. In this work, we propose a novel method for the analysis of multi-shell diffusion imaging data by incorporating compartment models into a spherical deconvolution framework for fiber orientation distribution (FOD) reconstruction. For numerical implementation, we develop an adaptively constrained energy minimization approach to efficiently compute the solution. On simulated and real data from Human Connectome Project (HCP), we show that our method not only reconstructs sharp and clean FODs for the modeling of fiber crossings, but also generates reliable estimation of compartment parameters with great potential for clinical research of neurological diseases. In comparisons with publicly available DSI-Studio and BEDPOSTX of FSL, we demonstrate that our method reconstructs sharper FODs with more precise estimation of fiber directions. By applying probabilistic tractography to the FODs computed by our method, we show that more complete reconstruction of the corpus callosum bundle can be achieved. On a clinical, two-shell diffusion imaging data, we also demonstrate the feasibility of our method in analyzing white matter lesions.

Simulation of changes in diffusion related to different pathologies at cellular level after traumatic brain injury

PURPOSE

The goal of the current study was to investigate tissue pathology at the cellular level in traumatic brain injury (TBI) as revealed by Monte Carlo simulation of diffusion tensor imaging (DTI)-derived parameters and elucidate the possible sources of conflicting findings of DTI abnormalities as reported in the TBI literature.

METHODS

A model with three compartments separated by permeable membranes was employed to represent the diffusion environment of water molecules in brain white matter. The dynamic diffusion process was simulated with a Monte Carlo method using adjustable parameters of intra-axonal diffusivity, axon separation, glial cell volume fraction, and myelin sheath permeability. The effects of tissue pathology on DTI parameters were investigated by adjusting the parameters of the model corresponding to different stages of brain injury.

RESULTS

The results suggest that the model is appropriate and the DTI-derived parameters simulate the predominant cellular pathology after TBI. Our results further indicate that when edema is not prevalent, axial and radial diffusivity have better sensitivity to axonal injury and demyelination than other DTI parameters.

CONCLUSION

DTI is a promising biomarker to detect and stage tissue injury after TBI. The observed inconsistencies among previous studies are likely due to scanning at different stages of tissue injury after TBI. Magn Reson Med 76:290-300, 2016. © 2015 Wiley Periodicals, Inc.

Effects of pore-size and shape distributions on diffusion pore imaging by nuclear magnetic resonance

In medical imaging and porous media research, NMR diffusion measurements are extensively used to investigate the structure of diffusion restrictions such as cell membranes. Recently, several methods have been proposed to unambiguously determine the shape of arbitrary closed pores or cells filled with an NMR-visible medium by diffusion experiments. The first approach uses a combination of a long and a short diffusion-weighting gradient pulse, while the other techniques employ short gradient pulses only. While the eventual aim of these methods is to determine pore-size and shape distributions, the focus has been so far on identical pores. Thus, the aim of this work is to investigate the ability of these different methods to resolve pore-size and orientation distributions. Simulations were performed comparing the various pore imaging techniques employing different distributions of pore size and orientation and varying timing parameters. The long-narrow gradient profile is most advantageous to investigate pore distributions, because average pore images can be directly obtained. The short-gradient methods suppress larger pores or induce a considerable blurring. Moreover, pore-shape-specific artifacts occur; for example, the central part of a distribution of cylinders may be largely underestimated. Depending on the actual pore distribution, short-gradient methods may nonetheless yield good approximations of the average pore shape. Furthermore, the application of short-gradient methods can be advantageous to differentiate whether pore-size distributions or intensity distributions, e.g., due to surface relaxation, are predominant.

Time-Dependent Influence of Cell Membrane Permeability on MR Diffusion Measurements

PURPOSE

To investigate the influence of cell membrane permeability on diffusion measurements over a broad range of diffusion times.

METHODS

Human myelogenous leukemia K562 cells were cultured and treated with saponin to selectively alter cell membrane permeability, resulting in a broad physiologically relevant range of 0.011-0.044 μm/ms. Apparent diffusion coefficient (ADC) values were acquired with the effective diffusion time (Δeff ) ranging from 0.42 to 3000 ms. Cosine-modulated oscillating gradient spin echo (OGSE) measurements were performed to achieve short Δeff from 0.42 to 5 ms, while stimulated echo acquisitions were used to achieve long Δeff from 11 to 2999 ms. Computer simulations were also performed to support the experimental results.

RESULTS

Both computer simulations and experiments in vitro showed that the influence of membrane permeability on diffusion MR measurements is highly dependent on the choice of diffusion time, and it is negligible only when the diffusion time is at least one order of magnitude smaller than the intracellular exchange lifetime.

CONCLUSION

The influence of cell membrane permeability on the measured ADCs is negligible in OGSE measurements at moderately high frequencies. By contrast, cell membrane permeability has a significant influence on ADC and quantitative diffusion measurements at low frequencies such as those sampled using conventional pulsed gradient methods.

Unraveling the multiscale structural organization and connectivity of the human brain: the role of diffusion MRI

The structural architecture and the anatomical connectivity of the human brain show different organizational principles at distinct spatial scales. Histological staining and light microscopy techniques have been widely used in classical neuroanatomical studies to unravel brain organization. Using such techniques is a laborious task performed on 2-dimensional histological sections by skilled anatomists possibly aided by semi-automated algorithms. With the recent advent of modern magnetic resonance imaging (MRI) contrast mechanisms, cortical layers and columns can now be reliably identified and their structural properties quantified post-mortem. These developments are allowing the investigation of neuroanatomical features of the brain at a spatial resolution that could be interfaced with that of histology. Diffusion MRI and tractography techniques, in particular, have been used to probe the architecture of both white and gray matter in three dimensions. Combined with mathematical network analysis, these techniques are increasingly influential in the investigation of the macro-, meso-, and microscopic organization of brain connectivity and anatomy, both in vivo and ex vivo. Diffusion MRI-based techniques in combination with histology approaches can therefore support the endeavor of creating multimodal atlases that take into account the different spatial scales or levels on which the brain is organized. The aim of this review is to illustrate and discuss the structural architecture and the anatomical connectivity of the human brain at different spatial scales and how recently developed diffusion MRI techniques can help investigate these.

Microstructural parameter estimation in vivo using diffusion MRI and structured prior information

Purpose

Diffusion MRI has recently been used with detailed models to probe tissue microstructure. Much of this work has been performed ex vivo with powerful scanner hardware, to gain sensitivity to parameters such as axon radius. By contrast, performing microstructure imaging on clinical scanners is extremely challenging.

Methods

We use an optimized dual spin‐echo diffusion protocol, and a Bayesian fitting approach, to obtain reproducible contrast (histogram overlap of up to 92%) in estimated maps of axon radius index in healthy adults at a modest, widely‐available gradient strength (35 mT m −1). A key innovation is the use of influential priors.

Results

We demonstrate that our priors can improve precision in axon radius estimates—a 7‐fold reduction in voxelwise coefficient of variation in vivo—without significant bias. Our results may reflect true axon radius differences between white matter regions, but this interpretation should be treated with caution due to the complexity of the tissue relative to our model.

Conclusions

Some sensitivity to relatively large axons (3–15 μm) may be available at clinical field and gradient strengths. Future applications at higher gradient strength will benefit from the favorable eddy current properties of the dual spin‐echo sequence, and greater precision available with suitable priors. Magn Reson Med, 2015. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. Magn Reson Med 75:1787–1796, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance.

Parameter estimation using macroscopic diffusion MRI signal models

Macroscopic models of the diffusion MRI (dMRI) signal can be helpful to understanding the relationship between the tissue microstructure and the dMRI signal. We study the least squares problem associated with estimating tissue parameters such as the cellular volume fraction, the residence times and the effective diffusion coefficients using a recently developed macroscopic model of the dMRI signal called the Finite Pulse Kärger model that generalizes the original Kärger model to non-narrow gradient pulses. In order to analyze the quality of the estimation in a controlled way, we generated synthetic noisy dMRI signals by including the effect of noise on the exact signal produced by the Finite Pulse Kärger model. The noisy signals were then fitted using the macroscopic model. Minimizing the least squares, we estimated the model parameters. The bias and standard deviations of the estimated model parameters as a function of the signal to noise ratio (SNR) were obtained. We discuss the choice of the b-values, the least square weights, the extension to experimentally obtained dMRI data as well noise correction.

Numerical study of a cylinder model of the diffusion MRI signal for neuronal dendrite trees

We study numerically how the neuronal dendrite tree structure can affect the diffusion magnetic resonance imaging (dMRI) signal in brain tissue. For a large set of randomly generated dendrite trees, synthetic dMRI signals are computed and fitted to a cylinder model to estimate the effective longitudinal diffusivity D(L) in the direction of neurites. When the dendrite branches are short compared to the diffusion length, D(L) depends significantly on the ratio between the average branch length and the diffusion length. In turn, D(L) has very weak dependence on the distribution of branch lengths and orientations of a dendrite tree, and the number of branches per node. We conclude that the cylinder model which ignores the connectivity of the dendrite tree, can still be adapted to describe the apparent diffusion coefficient in brain tissue.

Longitudinal measurement of the developing grey matter in preterm subjects using multi-modal MRI

Preterm birth is a major public health concern, with the severity and occurrence of adverse outcome increasing with earlier delivery. Being born preterm disrupts a time of rapid brain development: in addition to volumetric growth, the cortex folds, myelination is occurring and there are changes on the cellular level. These neurological events have been imaged non-invasively using diffusion-weighted (DW) MRI. In this population, there has been a focus on examining diffusion in the white matter, but the grey matter is also critically important for neurological health. We acquired multi-shell high-resolution diffusion data on 12 infants born at ≤ 28 weeks of gestational age at two time-points: once when stable after birth, and again at term-equivalent age. We used the Neurite Orientation Dispersion and Density Imaging model (NODDI) (Zhang et al., 2012) to analyse the changes in the cerebral cortex and the thalamus, both grey matter regions. We showed region-dependent changes in NODDI parameters over the preterm period, highlighting underlying changes specific to the microstructure. This work is the first time that NODDI parameters have been evaluated in both the cortical and the thalamic grey matter as a function of age in preterm infants, offering a unique insight into neuro-development in this at-risk population.

Ground Truth for Diffusion MRI in Cancer: A Model-Based Investigation of a Novel Tissue-Mimetic Material

This work presents preliminary results on the development, characterisation, and use of a novel physical phantom designed as a simple mimic of tumour cellular structure, for diffusion-weighted magnetic resonance imaging (DW-MRI) applications. The phantom consists of a collection of roughly spherical, micron-sized core-shell polymer 'cells', providing a system whose ground truth microstructural properties can be determined and compared with those obtained from modelling the DW-MRI signal. A two-compartment analytic model combining restricted diffusion inside a sphere with hindered extracellular diffusion was initially investigated through Monte Carlo diffusion simulations, allowing a comparison between analytic and simulated signals. The model was then fitted to DW-MRI data acquired from the phantom over a range of gradient strengths and diffusion times, yielding estimates of 'cell' size, intracellular volume fraction and the free diffusion coefficient. An initial assessment of the accuracy and precision of these estimates is provided, using independent scanning electron microscope measurements and bootstrap-style simulations. Such phantoms may be useful for testing microstructural models relevant to the characterisation of tumour tissue.

The impact of gradient strength on in vivo diffusion MRI estimates of axon diameter

Diffusion magnetic resonance imaging (MRI) methods for axon diameter mapping benefit from higher maximum gradient strengths than are currently available on commercial human scanners. Using a dedicated high-gradient 3T human MRI scanner with a maximum gradient strength of 300 mT/m, we systematically studied the effect of gradient strength on in vivo axon diameter and density estimates in the human corpus callosum. Pulsed gradient spin echo experiments were performed in a single scan session lasting approximately 2h on each of three human subjects. The data were then divided into subsets with maximum gradient strengths of 77, 145, 212, and 293 mT/m and diffusion times encompassing short (16 and 25 ms) and long (60 and 94 ms) diffusion time regimes. A three-compartment model of intra-axonal diffusion, extra-axonal diffusion, and free diffusion in cerebrospinal fluid was fitted to the data using a Markov chain Monte Carlo approach. For the acquisition parameters, model, and fitting routine used in our study, it was found that higher maximum gradient strengths decreased the mean axon diameter estimates by two to three fold and decreased the uncertainty in axon diameter estimates by more than half across the corpus callosum. The exclusive use of longer diffusion times resulted in axon diameter estimates that were up to two times larger than those obtained with shorter diffusion times. Axon diameter and density maps appeared less noisy and showed improved contrast between different regions of the corpus callosum with higher maximum gradient strength. Known differences in axon diameter and density between the genu, body, and splenium of the corpus callosum were preserved and became more reproducible at higher maximum gradient strengths. Our results suggest that an optimal q-space sampling scheme for estimating in vivo axon diameters should incorporate the highest possible gradient strength. The improvement in axon diameter and density estimates that we demonstrate from increasing maximum gradient strength will inform protocol development and encourage the adoption of higher maximum gradient strengths for use in commercial human scanners.

One diffusion acquisition and different white matter models: how does microstructure change in human early development based on WMTI and NODDI?

White matter microstructural changes during the first three years of healthy brain development are characterized using two different models developed for limited clinical diffusion data: White Matter Tract Integrity (WMTI) metrics from Diffusional Kurtosis Imaging (DKI) and Neurite Orientation Dispersion and Density Imaging (NODDI). Both models reveal a non-linear increase in intra-axonal water fraction and in tortuosity of the extra-axonal space as a function of age, in the genu and splenium of the corpus callosum and the posterior limb of the internal capsule. The changes are consistent with expected behavior related to myelination and asynchrony of fiber development. The intra- and extracellular axial diffusivities as estimated with WMTI do not change appreciably in normal brain development. The quantitative differences in parameter estimates between models are examined and explained in the light of each model's assumptions and consequent biases, as highlighted in simulations. Finally, we discuss the feasibility of a model with fewer assumptions.

Mapping mean axon diameter and axonal volume fraction by MRI using temporal diffusion spectroscopy

Mapping mean axon diameter and intra-axonal volume fraction may have significant clinical potential because nerve conduction velocity is directly dependent on axon diameter, and several neurodegenerative diseases affect axons of specific sizes and alter axon counts. Diffusion-weighted MRI methods based on the pulsed gradient spin echo (PGSE) sequence have been reported to be able to assess axon diameter and volume fraction non-invasively. However, due to the relatively long diffusion times used, e.g. >20ms, the sensitivity to small axons (diameter<2μm) is low, and the derived mean axon diameter has been reported to be overestimated. In the current study, oscillating gradient spin echo (OGSE) diffusion sequences with variable frequency gradients were used to assess rat spinal white matter tracts with relatively short effective diffusion times (1-5ms). In contrast to previous PGSE-based methods, the extra-axonal diffusion cannot be modeled as hindered (Gaussian) diffusion when short diffusion times are used. Appropriate frequency-dependent rates are therefore incorporated into our analysis and validated by histology-based computer simulation of water diffusion. OGSE data were analyzed to derive mean axon diameters and intra-axonal volume fractions of rat spinal white matter tracts (mean axon diameter of ~1.27-5.54μm). The estimated values were in good agreement with histology, including the small axon diameters (<2.5μm). This study establishes a framework for the quantification of nerve morphology using the OGSE method with high sensitivity to small axons.

Exploring diffusion across permeable barriers at high gradients. II. Localization regime

We present an analytical solution of the one-dimensional Bloch-Torrey equation for diffusion across multiple semi-permeable barrier. This solution generalizes the seminal work by Stoller, Happer, and Dyson, in which the non-Gaussian stretched-exponential behavior of the pulsed-gradient spin-echo (PGSE) signal was first predicted at high gradients in the so-called localization regime. We investigate how the diffusive exchange across a semi-permeable barrier modifies this asymptotic behavior, and explore the transition between the localization regime at low permeability and the Gaussian regime at high permeability. High gradients are suitable to spatially localize the contribution of the nuclei near the barrier and to enhance the sensitivity of the PGSE signal to the barrier permeability. The emergence of the localization regime for three-dimensional domains is discussed.

Numerical study of a macroscopic finite pulse model of the diffusion MRI signal

Diffusion magnetic resonance imaging (dMRI) is an imaging modality that probes the diffusion characteristics of a sample via the application of magnetic field gradient pulses. The dMRI signal from a heterogeneous sample includes the contribution of the water proton magnetization from all spatial positions in a voxel. If the voxel can be spatially divided into different Gaussian diffusion compartments with inter-compartment exchange governed by linear kinetics, then the dMRI signal can be approximated using the macroscopic Karger model, which is a system of coupled ordinary differential equations (ODEs), under the assumption that the duration of the diffusion-encoding gradient pulses is short compared to the diffusion time (the narrow pulse assumption). Recently, a new macroscopic model of the dMRI signal, without the narrow pulse restriction, was derived from the Bloch-Torrey partial differential equation (PDE) using periodic homogenization techniques. When restricted to narrow pulses, this new homogenized model has the same form as the Karger model. We conduct a numerical study of the new homogenized model for voxels that are made up of periodic copies of a representative volume that contains spherical and cylindrical cells of various sizes and orientations and show that the signal predicted by the new model approaches the reference signal obtained by solving the full Bloch-Torrey PDE in O(ε(2)), where ε is the ratio between the size of the representative volume and a measure of the diffusion length. When the narrow gradient pulse assumption is not satisfied, the new homogenized model offers a much better approximation of the full PDE signal than the Karger model. Finally, preliminary results of applying the new model to a voxel that is not made up of periodic copies of a representative volume are shown and discussed.

Exploring diffusion across permeable barriers at high gradients. I. Narrow pulse approximation

The adaptive variation of the gradient intensity with the diffusion time at a constant optimal b-value is proposed to enhance the contribution of the nuclei diffusing across permeable barriers, to the pulsed-gradient spin-echo (PGSE) signal. An exact simple formula the PGSE signal is derived under the narrow pulse approximation in the case of one-dimensional diffusion across a single permeable barrier. The barrier contribution to the signal is shown to be maximal at a particular b-value. The exact formula is then extended to multiple permeable barriers, while the PGSE signal is shown to be sensitive to the permeability and to the inter-barrier distance. Potential applications of the protocol to survey diffusion in three-dimensional domains with permeable membranes are illustrated through numerical simulations.

Accelerated Microstructure Imaging via Convex Optimization (AMICO) from diffusion MRI data

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.

High angular resolution diffusion imaging with stimulated echoes: compensation and correction in experiment design and analysis

Stimulated echo acquisition mode (STEAM) diffusion MRI can be advantageous over pulsed-gradient spin-echo (PGSE) for diffusion times that are long compared with T2 . It therefore has potential for biomedical diffusion imaging applications at 7T and above where T2 is short. However, gradient pulses other than the diffusion gradients in the STEAM sequence contribute much greater diffusion weighting than in PGSE and lead to a disrupted experimental design. Here, we introduce a simple compensation to the STEAM acquisition that avoids the orientational bias and disrupted experiment design that these gradient pulses can otherwise produce. The compensation is simple to implement by adjusting the gradient vectors in the diffusion pulses of the STEAM sequence, so that the net effective gradient vector including contributions from diffusion and other gradient pulses is as the experiment intends. High angular resolution diffusion imaging (HARDI) data were acquired with and without the proposed compensation. The data were processed to derive standard diffusion tensor imaging (DTI) maps, which highlight the need for the compensation. Ignoring the other gradient pulses, a bias in DTI parameters from STEAM acquisition is found, due both to confounds in the analysis and the experiment design. Retrospectively correcting the analysis with a calculation of the full B matrix can partly correct for these confounds, but an acquisition that is compensated as proposed is needed to remove the effect entirely.

A model for extra-axonal diffusion spectra with frequency-dependent restriction

PURPOSE

In the brain, there is growing interest in using the temporal diffusion spectrum to characterize axonal geometry in white matter because of the potential to be more sensitive to small pores compared to conventional time-dependent diffusion. However, analytical expressions for the diffusion spectrum of particles have only been derived for simple, restricting geometries such as cylinders, which are often used as a model for intra-axonal diffusion. The extra-axonal space is more complex, but the diffusion spectrum has largely not been modeled. We propose a model for the extra-axonal space, which can be used for interpretation of experimental data.

THEORY AND METHODS

An empirical model describing the extra-axonal space diffusion spectrum was compared with simulated spectra. Spectra were simulated using Monte Carlo methods for idealized, regularly and randomly packed axons over a wide range of packing densities and spatial scales. The model parameters are related to the microstructural properties of tortuosity, axonal radius, and separation for regularly packed axons and pore size for randomly packed axons.

RESULTS

Forward model predictions closely matched simulations. The model fitted the simulated spectra well and accurately estimated microstructural properties.

CONCLUSIONS

This simple model provides expressions that relate the diffusion spectrum to biologically relevant microstructural properties.

Quantifying white matter tract diffusion parameters in the presence of increased extra-fiber cellularity and vasogenic edema

The effect of extra-fiber structural and pathological components confounding diffusion tensor imaging (DTI) computation was quantitatively investigated using data generated by both Monte-Carlo simulations and tissue phantoms. Increased extent of vasogenic edema, by addition of various amount of gel to fixed normal mouse trigeminal nerves or by increasing non-restricted isotropic diffusion tensor components in Monte-Carlo simulations, significantly decreased fractional anisotropy (FA) and increased radial diffusivity, while less significantly increased axial diffusivity derived by DTI. Increased cellularity, mimicked by graded increase of the restricted isotropic diffusion tensor component in Monte-Carlo simulations, significantly decreased FA and axial diffusivity with limited impact on radial diffusivity derived by DTI. The MC simulation and tissue phantom data were also analyzed by the recently developed diffusion basis spectrum imaging (DBSI) to simultaneously distinguish and quantify the axon/myelin integrity and extra-fiber diffusion components. Results showed that increased cellularity or vasogenic edema did not affect the DBSI-derived fiber FA, axial or radial diffusivity. Importantly, the extent of extra-fiber cellularity and edema estimated by DBSI correlated with experimentally added gel and Monte-Carlo simulations. We also examined the feasibility of applying 25-direction diffusion encoding scheme for DBSI analysis on coherent white matter tracts. Results from both phantom experiments and simulations suggested that the 25-direction diffusion scheme provided comparable DBSI estimation of both fiber diffusion parameters and extra-fiber cellularity/edema extent as those by 99-direction scheme. An in vivo 25-direction DBSI analysis was performed on experimental autoimmune encephalomyelitis (EAE, an animal model of human multiple sclerosis) optic nerve as an example to examine the validity of derived DBSI parameters with post-imaging immunohistochemistry verification. Results support that in vivo DBSI using 25-direction diffusion scheme correctly reflect the underlying axonal injury, demyelination, and inflammation of optic nerves in EAE mice.

Functional correlates of central white matter maturation in perinatal period in rabbits

Anisotropy indices derived from diffusion tensor imaging (DTI) are being increasingly used as biomarkers of central WM structural maturation, myelination and even functional development. Our hypothesis was that the rate of functional changes in central WM tracts directly reflects rate of changes in structural development as determined by DTI indices. We examined structural and functional development of four major central WM tracts with different maturational trajectories, including internal capsule (IC), corpus callosum (CC), fimbria hippocampi (FH) and anterior commissure (AC). Rabbits were chosen due to perinatal brain development being similar to humans, and four time points were studied: P1, P11, P18 and adults. Imaging parameters of structural maturation included fractional anisotropy (FA), mean and directional diffusivities derived from DTI, and T2 relaxation time. Axonal composition and degree of myelination were confirmed on electron microscopy. To assess functional maturation, conduction velocity was measured in myelinated and non-myelinated fibers by electrophysiological recordings of compound action potential in perfused brain slices. Diffusion indices and T2 relaxation time in rabbits followed a sigmoid curve during development similar to that for humans, with active changes even at premyelination stage. The shape of the developmental curve was different between the fiber tracts, with later onset but steeper rapid phase of development in IC and FH than in CC. The structural development was not directly related to myelination or to functional development. Functional properties in projection (IC) and limbic tracts (FH) matured earlier than in associative and commissural tracts (CC and AC). The rapid phase of changes in diffusion anisotropy and T2 relaxation time coincided with the development of functional responses and myelination in IC and FH between the second and third weeks of postnatal development in rabbits. In these two tracts, MRI indices could serve as surrogate markers of the early stage of myelination. However, the discordance between developmental change of diffusion indices, myelination and functional properties in CC and AC cautions against equating DTI index changes as biomarkers for myelination in all tracts.

Can T2-spectroscopy resolve submicrometer axon diameters?

The microscopic geometry of white matter carries rich information about brain function in health and disease. A key challenge for medical imaging is to estimate microstructural features noninvasively. One important parameter is the axon diameter, which correlates with the conduction time delay of action potentials and is affected by various neurological disorders. Diffusion magnetic resonance (MR) experiments are the method of choice today when we aim to recover the axon diameter distribution, although the technique requires very high gradient strengths in order to assess nerve fibers with one micrometer or less in diameter. In practice in-vivo brain imaging is only sensitive to the largest axons, not least due to limitations in the human physiology which tolerates only moderate gradient strengths. This work studies, from a theoretical perspective, the feasibility of T2-spectroscopy to resolve submicrometer tissue structures. Exploiting the surface relaxation effect, we formulate a plausible biophysical model relating the axon diameter distribution to the T2-weighted signal, which is based on a surface-to-volume ratio approximation of the Bloch-Torrey equation. Under a certain regime of bulk and surface relaxation coefficients, our simulation results suggest that it might be possible to reveal axons smaller than one micrometer in diameter.

An analytical model for estimating water exchange rate in white matter using diffusion MRI

Substantial effort is being expended on using micro-structural modeling of the white matter, with the goal of relating diffusion weighted magnetic resonance imaging (DWMRI) to the underlying structure of the tissue, such as axonal density. However, one of the important parameters affecting diffusion is the water exchange rate between the intra- and extra-axonal space, which has not been fully investigated and is a crucial marker of brain injury such as multiple sclerosis (MS), stroke, and traumatic brain injury (TBI). To our knowledge, there is no diffusion analytical model which includes the Water eXchange Rate (WXR) without the requirement of short gradient pulse (SGP) approximation. We therefore propose a new analytical model by deriving the diffusion signal for a permeable cylinder, assuming a clinically feasible pulse gradient spin echo (PGSE) sequence. Simulations based on Markov Random Walk confirm that the exchange parameter included in our model has a linear correlation (R²>0.88) with the actual WXR. Moreover, increasing WXR causes the estimated values of diameter and volume fraction of the cylinders to increase and decrease, respectively, which is consistent with our findings from histology measurements in tissues near TBI regions. This model was also applied to the diffusion signal acquired from ex vivo brains of 14 male (10 TBI and 4 normal) rats using hybrid diffusion imaging. The estimated values of axon diameter and axonal volume fraction are in agreement with their corresponding histological measurements in normal brains, with 0.96 intra-class correlation coefficient value resulting from consistency analysis. Moreover, a significant increase (p = 0.001) in WXR and diameter and decrease in axonal volume fraction in the TBI boundary were detected in the TBI rats compared with the normal rats.