No Arabic abstract
Computational models of biophysical tissue properties have been widely used in diffusion MRI (dMRI) research to elucidate the link between microstructural properties and MR signal formation. For brain tissue, the research community has developed the so-called Standard Model (SM) that has been widely used. However, in clinically applicable acquisition protocols, the inverse problem that recovers the SM parameters from a set of MR diffusion measurements using pairs of short pulsed field gradients was shown to be ill-posed. Multidimensional dMRI was shown to solve this problem by combining linear and planar tensor encoding data. Given sufficient measurements, multiple choices of b-tensor sets provide enough information to estimate all SM parameters. However, in the presence of noise, some sets will provide better results. In this work, we develop a framework for optimal experimental design of multidimensional dMRI sequences applicable to the SM. This framework is based on maximising the determinant of the Fisher information matrix, which is averaged over the full SM parameter space. This averaging provides a fairly objective information metric tailored for the expected signal but that only depends on the acquisition configuration. The optimisation of this metric can be further restricted to any subclass of desirable design constraints like, for instance, hardware-specific constraints. In this work, we compute the optimal acquisitions over the set of all b-tensors with fixed eigenvectors.
Diffusion tractography is routinely used to study white matter architecture and brain connectivity in vivo. A key step for successful tractography of neuronal tracts is the correct identification of tract directions in each voxel. Here we propose a fingerprinting-based methodology to identify these fiber directions in Orientation Distribution Functions, dubbed ODF-Fingerprinting (ODF-FP). In ODF-FP, fiber configurations are selected based on the similarity between measured ODFs and elements in a pre-computed library. In noisy ODFs, the library matching algorithm penalizes the more complex fiber configurations. ODF simulations and analysis of bootstrapped partial and whole-brain in vivo datasets show that the ODF-FP approach improves the detection of fiber pairs with small crossing angles while maintaining fiber direction precision, which leads to better tractography results. Rather than focusing on the ODF maxima, the ODF-FP approach uses the whole ODF shape to infer fiber directions to improve the detection of fiber bundles with small crossing angle. The resulting fiber directions aid tractography algorithms in accurately displaying neuronal tracts and calculating brain connectivity.
Purpose: Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, the general Standard Model has recently shown that model parameter estimation from dMRI data is ill-posed unless very strong magnetic gradients are used. We analyse this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from Single Diffusion Encoding (SDE) to Double Diffusion Encoding (DDE) solves the ill-posedness and increases the accuracy of the parameter estimation. Methods: We analyse theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions. Results: We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space. Conclusions: DDE adds additional information for estimating the model parameters, unexplored by SDE, which is enough to solve the degeneracy in the NODDIDA model parameter estimation.
Purpose: Diffusion-weighted steady-state free precession (DW-SSFP) is shown to provide a means to probe non-Gaussian diffusion through manipulation of the flip angle. A framework is presented to define an effective b-value in DW-SSFP. Theory: The DW-SSFP signal is a summation of coherence pathways with different b-values. The relative contribution of each pathway is dictated by the flip angle. This leads to an apparent diffusion coefficient (ADC) estimate that depends on the flip angle in non-Gaussian diffusion regimes. By acquiring DW-SSFP data at multiple flip angles and modelling the variation in ADC for a given form of non-Gaussianity, the ADC can be estimated at a well-defined effective b-value. Methods: A gamma distribution is used to model non-Gaussian diffusion, embedded in the Buxton signal model for DW-SSFP. Monte-Carlo simulations of non-Gaussian diffusion in DW-SSFP and diffusion-weighted spin-echo (DW-SE) sequences are used to verify the proposed framework. Dependence of ADC on flip angle in DW-SSFP is verified with experimental measurements in a whole, human post-mortem brain. Results: Monte-Carlo simulations reveal excellent agreement between ADCs estimated with DW-SE and the proposed framework. Experimental ADC estimates vary as a function of flip angle over the corpus callosum of the postmortem brain, estimating the mean and standard deviation of the gamma distribution as $1.50cdot 10^{-4} mm^2/s$ and $2.10cdot 10^{-4} mm^2/s$. Conclusion: DW-SSFP can be used to investigate non-Gaussian diffusion by varying the flip angle. By fitting a model of non-Gaussian diffusion, the ADC in DW-SSFP can be estimated at an effective b-value, comparable to more conventional diffusion sequences.
In utero diffusion MRI provides unique opportunities to non-invasively study the microstructure of tissue during fetal development. A wide range of developmental processes, such as the growth of white matter tracts in the brain, the maturation of placental villous trees, or the fibres in the fetal heart remain to be studied and understood in detail. Advances in fetal interventions and surgery furthermore increase the need for ever more precise antenatal diagnosis from fetal MRI. However, the specific properties of the in utero environment, such as fetal and maternal motion, increased field-of-view, tissue interfaces and safety considerations, are significant challenges for most MRI techniques, and particularly for diffusion. Recent years have seen major improvements, driven by the development of bespoke techniques adapted to these specific challenges in both acquisition and processing. Fetal diffusion MRI, an emerging research tool, is now adding valuable novel information for both research and clinical questions. This paper will highlight specific challenges, outline strategies to target them, and discuss two main applications: fetal brain connectomics and placental maturation.
Diffusion MRI offers a unique probe into neural microstructure and connectivity in the developing brain. However, analysis of neonatal brain imaging data is complicated by inevitable subject motion, leading to a series of scattered slices that need to be aligned within and across diffusion-weighted contrasts. Here, we develop a reconstruction method for scattered slice multi-shell high angular resolution diffusion imaging (HARDI) data, jointly estimating an uncorrupted data representation and motion parameters at the slice or multiband excitation level. The reconstruction relies on data-driven representation of multi-shell HARDI data using a bespoke spherical harmonics and radial decomposition (SHARD), which avoids imposing model assumptions, thus facilitating to compare various microstructure imaging methods in the reconstructed output. Furthermore, the proposed framework integrates slice-level outlier rejection, distortion correction, and slice profile correction. We evaluate the method in the neonatal cohort of the developing Human Connectome Project (650 scans). Validation experiments demonstrate accurate slice-level motion correction across the age range and across the range of motion in the population. Results in the neonatal data show successful reconstruction even in severely motion-corrupted subjects. In addition, we illustrate how local tissue modelling can extract advanced microstructure features such as orientation distribution functions from the motion-corrected reconstructions.