To compare the diagnostic value of DKI and DTI in differentiation of high-grade glioma recurrence and pseudoprogression (PsP).

Diffusion tensor imaging (DTI) is a powerful and noninvasive imaging method for characterizing tissue microstructure and white matter organization in the brain. While it has been applied extensively in research studies of the human brain, DTI studies of non-human primates have been performed only recently. The growing application of DTI in rhesus monkey studies would significantly benefit from a standardized framework to compare findings across different studies. A very common strategy for image analysis is to spatially normalize (co-register) the individual scans to a representative template space. This paper presents the development of a DTI brain template, UWRMAC-DTI271, for adolescent Rhesus Macaque (Macaca mulatta) monkeys. The template was generated from 271 rhesus monkeys, collected as part of a unique brain imaging genetics study. It is the largest number of animals ever used to generate a computational brain template, which enables the generation of a template that has high image quality and accounts for variability in the species. The quality of the template is further ensured with the use of DTI-TK, a well-tested and high-performance DTI spatial normalization method in human studies. We demonstrated its efficacy in monkey studies for the first time by comparing it to other commonly used scalar-methods for DTI normalization. It is anticipated that this template will play an important role in facilitating cross-site voxelwise DTI analyses in Rhesus Macaques. Such analyses are crucial in investigating the role of white matter structure in brain function, development, and other psychopathological disorders for which there are well-validated non-human primate models.

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.

This paper introduces image quality transfer. The aim is to learn the fine structural detail of medical images from high quality data sets acquired with long acquisition times or from bespoke devices and transfer that information to enhance lower quality data sets from standard acquisitions. We propose a framework for solving this problem using random forest regression to relate patches in the low-quality data set to voxel values in the high quality data set. Two examples in diffusion MRI demonstrate the idea. In both cases, we learn from the Human Connectome Project (HCP) data set, which uses an hour of acquisition time per subject, just for diffusion imaging, using custom built scanner hardware and rapid imaging techniques. The first example, super-resolution of diffusion tensor images (DTIs), enhances spatial resolution of standard data sets with information from the high-resolution HCP data. The second, parameter mapping, constructs neurite orientation density and dispersion imaging (NODDI) parameter maps, which usually require specialist data sets with two b-values, from standard single-shell high angular resolution diffusion imaging (HARDI) data sets with b = 1000 smm-2. Experiments quantify the improvement against alternative image reconstructions in comparison to ground truth from the HCP data set in both examples and demonstrate efficacy on a standard data set.

This paper compares a range of compartment models for diffusion MRI data on in vivo human acquisitions from a standard 60mT/m system (Philips 3T Achieva) and a unique 300mT/m system (Siemens Connectom). The key aim is to determine whether both systems support broadly the same models or whether the Connectom higher gradient system supports significantly more complex models. A single volunteer underwent 8h of acquisition on each system to provide uniquely wide and dense sampling of the available space of pulsed-gradient spin-echo (PGSE) measurements. We select a set of promising models from the wide set of possible three-compartment models for in vivo white matter (WM) that previous work and preliminary experiments suggest as strong candidates, but extend them to fit for compartmental T2 and diffusivity. We focus on the corpus callosum where the WM fibre architecture is simplest and compare their ability to explain the measured data, using Akaike's information criterion (AIC), and to predict unseen data, using cross-validation. We also compare the stability of parameter estimates in the presence of i) noise, using bootstrapping, and ii) spatial variation, using visual assessment and comparison with anatomical knowledge. Broadly similar models emerge from the AIC and cross-validation experiments in both data sets. Specifically, a three-compartment model consisting of either a Bingham distribution of sticks or a Cylinder for the intracellular compartment, an anisotropic diffusion tensor (DT) model for the extracellular compartment, as well as an isotropic CSF compartment, performs consistently well. However, various other models also perform well and no single model emerges as clear winner. The WM data (with virtually no CSF contamination) do not support compartmental T2 but partially support compartmental diffusivity. Evaluation of parameter stability favours simpler models than those identified by AIC or cross-validation. They suggest that the level of complexity in models underpinning currently popular microstructure imaging techniques such as NODDI, CHARMED, or ActiveAx, where the number of free parameters is about 4 or 5 rather than 10 or 11, may reflect the level of complexity achievable for a useful technique on current systems, although the 300mT/m data may support more complex models.

Cognitive impairment after traumatic brain injury remains hard to predict. This is partly because axonal injury, which is of fundamental importance, is difficult to measure clinically. Advances in MRI allow axonal injury to be detected after traumatic brain injury, but the most sensitive approach is unclear. Here, we compare the performance of diffusion tensor imaging, neurite orientation dispersion and density-imaging and volumetric measures of brain atrophy in the identification of white-matter abnormalities after traumatic brain injury. Thirty patients with moderate-severe traumatic brain injury in the chronic phase and 20 age-matched controls had T1-weighted and diffusion MRI. Neuropsychological tests of processing speed, executive functioning and memory were used to detect cognitive impairment. Extensive abnormalities in neurite density index and orientation dispersion index were observed, with distinct spatial patterns. Fractional anisotropy and mean diffusivity also indicated widespread abnormalities of white-matter structure. Neurite density index was significantly correlated with processing speed. Slower processing speed was also related to higher mean diffusivity in the corticospinal tracts. Lower white-matter volumes were seen after brain injury with greater effect sizes compared to diffusion metrics; however, volume was not sensitive to changes in cognitive performance. Volume was the most sensitive at detecting change between groups but was not specific for determining relationships with cognition. Abnormalities in fractional anisotropy and mean diffusivity were the most sensitive diffusion measures; however, neurite density index and orientation dispersion index may be more spatially specific. Lower neurite density index may be a useful metric for examining slower processing speed.

This paper presents Contextual Fibre Growth (ConFiG), an approach to generate white matter numerical phantoms by mimicking natural fibre genesis. ConFiG grows fibres one-by-one, following simple rules motivated by real axonal guidance mechanisms. These simple rules enable ConFiG to generate phantoms with tuneable microstructural features by growing fibres while attempting to meet morphological targets such as user-specified density and orientation distribution. We compare ConFiG to the state-of-the-art approach based on packing fibres together by generating phantoms in a range of fibre configurations including crossing fibre bundles and orientation dispersion. Results demonstrate that ConFiG produces phantoms with up to 20% higher densities than the state-of-the-art, particularly in complex configurations with crossing fibres. We additionally show that the microstructural morphology of ConFiG phantoms is comparable to real tissue, producing diameter and orientation distributions close to electron microscopy estimates from real tissue as well as capturing complex fibre cross sections. Signals simulated from ConFiG phantoms match real diffusion MRI data well, showing that ConFiG phantoms can be used to generate realistic diffusion MRI data. This demonstrates the feasibility of ConFiG to generate realistic synthetic diffusion MRI data for developing and validating microstructure modelling approaches.

In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. In the literature, improved alignment is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We apply this framework in three ways. Firstly we assess the most commonly used method for artefact correction, FSL's eddy_correct, and compare it to a recently proposed alternative, eddy. We demonstrate quantitatively that using eddy_correct leads to significant errors in the corrected data, whilst eddy is able to provide much improved correction. Secondly we investigate the datasets required to achieve good correction with eddy, by looking at the minimum number of directions required and comparing the recommended full-sphere acquisitions to equivalent half-sphere protocols. Finally, we investigate the impact of correction quality by examining the fits from microstructure models to real and simulated data.

In this paper we evaluate the three main methods for correcting the susceptibility-induced artefact in diffusion-weighted magnetic-resonance (DW-MR) data, and assess how correction is affected by the susceptibility field's interaction with motion. The susceptibility artefact adversely impacts analysis performed on the data and is typically corrected in post-processing. Correction strategies involve either registration to a structural image, the application of an acquired field-map or the use of additional images acquired with different phase-encoding. Unfortunately, the choice of which method to use is made difficult by the absence of any systematic comparisons of them. In this work we quantitatively evaluate these methods, by extending and employing a recently proposed framework that allows for the simulation of realistic DW-MR datasets with artefacts. Our analysis separately evaluates the ability for methods to correct for geometric distortions and to recover lost information in regions of signal compression. In terms of geometric distortions, we find that registration-based methods offer the poorest correction. Field-mapping techniques are better, but are influenced by noise and partial volume effects, whilst multiple phase-encode methods performed best. We use our simulations to validate a popular surrogate metric of correction quality, the comparison of corrected data acquired with AP and LR phase-encoding, and apply this surrogate to real datasets. Furthermore, we demonstrate that failing to account for the interaction of the susceptibility field with head movement leads to increased errors when analysing DW-MR data. None of the commonly used post-processing methods account for this interaction, and we suggest this may be a valuable area for future methods development.

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/mm2 (or 3 ​ms/μm2), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm2 or 3 ​ms/μm2). 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 (bmax ​= ​40,000 ​s/mm2, or 40 ​ms/μm2) and human (bmax ​= ​10,000 ​s/mm2, or 10 ​ms/μm2) 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.

Most motion correction methods work by aligning a set of volumes together, or to a volume that represents a reference location. These are based on an implicit assumption that the subject remains motionless during the several seconds it takes to acquire all slices in a volume, and that any movement occurs in the brief moment between acquiring the last slice of one volume and the first slice of the next. This is clearly an approximation that can be more or less good depending on how long it takes to acquire one volume and in how rapidly the subject moves. In this paper we present a method that increases the temporal resolution of the motion correction by modelling movement as a piecewise continous function over time. This intra-volume movement correction is implemented within a previously presented framework that simultaneously estimates distortions, movement and movement-induced signal dropout. We validate the method on highly realistic simulated data containing all of these effects. It is demonstrated that we can estimate the true movement with high accuracy, and that scalar parameters derived from the data, such as fractional anisotropy, are estimated with greater fidelity when data has been corrected for intra-volume movement. Importantly, we also show that the difference in fidelity between data affected by different amounts of movement is much reduced when taking intra-volume movement into account. Additional validation was performed on data from a healthy volunteer scanned when lying still and when performing deliberate movements. We show an increased correspondence between the "still" and the "movement" data when the latter is corrected for intra-volume movement. Finally we demonstrate a big reduction in the telltale signs of intra-volume movement in data acquired on elderly subjects.

To compare the added value of diffusion kurtosis imaging (DKI) with the combination of dynamic susceptibility contrast-enhanced (DSC) MRI in differentiating glioma recurrence from pseudoprogression.

To date, numerical simulations of the brain tissue have been limited by their lack of realism and flexibility. The purpose of this work is to propose a controlled and flexible generative model for brain cell morphology and an efficient computational pipeline for the reliable and robust simulation of realistic cellular structures with application to numerical simulation of intra-cellular diffusion-weighted MR (DW-MR) signal features. Inspired by the advances in computational neuroscience for modelling brain cells, we propose a generative model that enables users to simulate molecular diffusion within realistic digital brain cells, such as neurons, in a completely controlled and flexible fashion. We validate our new approach by showing an excellent match between the morphology (no statistically different 3D Sholl metrics, P > 0.05) and simulated intra-cellular DW-MR signal (mean relative difference < 2%) of the generated digital model of brain cells and those of digital reconstruction of real brain cells from available open-access databases. We demonstrate the versatility and potential of the framework by showing a select set of examples of relevance for the DW-MR community. The computational models introduced here are useful for synthesizing intra-cellular DW-MR signals, similar to those one might measure from brain metabolites DW-MRS experiments. They also provide the foundation for a more complete simulation system that will potentially include signals from extra-cellular compartments and exchange processes, necessary for synthesizing DW-MR signals of relevance for DW-MRI experiments.

We introduce a novel image-processing framework for tracking longitudinal changes in white matter microstructure using diffusion tensor imaging (DTI). Charting the trajectory of such temporal changes offers new insight into disease progression but to do so accurately faces a number of challenges. Recent developments have highlighted the importance of processing each subject's data at multiple time points in an unbiased way. In this paper, we aim to highlight a different challenge critical to the processing of longitudinal DTI data, namely the approach to image alignment. Standard approaches in the literature align DTI data by registering the corresponding scalar-valued fractional anisotropy (FA) maps. We propose instead a DTI registration algorithm that leverages full tensor information to drive improved alignment. This proposed pipeline is evaluated against the standard FA-based approach using a DTI dataset from an ongoing study of Alzheimer's disease (AD). The dataset consists of subjects scanned at two time points and at each time point the DTI acquisition consists of two back-to-back repeats in the same scanning session. The repeated scans allow us to evaluate the specificity of each pipeline, using a test-retest design, and assess precision, using bootstrap-based method. The results show that the tensor-based pipeline achieves both higher specificity and precision than the standard FA-based approach. Tensor-based registration for longitudinal processing of DTI data in clinical studies may be of particular value in studies assessing disease progression.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) are widely used models to infer microstructural features in the brain from diffusion-weighted MRI. Several studies have recently applied both models to increase sensitivity to biological changes, however, it remains uncertain how these measures are associated. Here we show that cortical distributions of DTI and NODDI are associated depending on the choice of b-value, a factor reflecting strength of diffusion weighting gradient. We analyzed a combination of high, intermediate and low b-value data of multi-shell diffusion-weighted MRI (dMRI) in healthy 456 subjects of the Human Connectome Project using NODDI, DTI and a mathematical conversion from DTI to NODDI. Cortical distributions of DTI and DTI-derived NODDI metrics were remarkably associated with those in NODDI, particularly when applied highly diffusion-weighted data (b-value = 3000 sec/mm2). This was supported by simulation analysis, which revealed that DTI-derived parameters with lower b-value datasets suffered from errors due to heterogeneity of cerebrospinal fluid fraction and partial volume. These findings suggest that high b-value DTI redundantly parallels with NODDI-based cortical neurite measures, but the conventional low b-value DTI is hard to reasonably characterize cortical microarchitecture.

The histopathological basis of "unidentified bright objects" (UBOs) (hyperintense regions seen on T2-weighted magnetic resonance (MR) brain scans in neurofibromatosis-1 (NF1)) remains unclear. New in vivo MRI-based techniques (multi-exponential T2 relaxation (MET2) and diffusion MR imaging (dMRI)) provide measures relating to microstructural change. We combined these methods and present previously unreported data on in vivo UBO microstructure in NF1.

Previous studies have demonstrated that the accumulation of amyloid-β (Aβ) pathologies has distinctive stage-specific effects on the structural and functional brain networks along the Alzheimer's disease (AD) continuum. A more comprehensive account of both types of brain network may provide a better characterization of the stage-specific effects of Aβ pathologies. A potential candidate for this joint characterization is the coupling between the structural and functional brain networks (SC-FC coupling). We therefore investigated the effect of Aβ accumulation on global SC-FC coupling in patients with mild cognitive impairment (MCI), AD, and healthy controls. Patients with MCI were dichotomized according to their level of Aβ pathology seen in 18F-flutemetamol PET-CT scans-namely, Aβ-negative and Aβ-positive. Our results show that there was no difference in global SC-FC coupling between different cohorts. During the prodromal AD stage, there was a significant negative correlation between the level of Aβ pathology and the global SC-FC coupling of MCI patients with positive Aβ, but no significant correlation for MCI patients with negative Aβ. During the AD dementia stage, the correlation between Aβ pathology and global SC-FC coupling in patients with AD was positive. Our results suggest that Aβ pathology has distinctive stage-specific effects on global coupling between the structural and functional brain networks along the AD continuum.

There is increasing interest in improving understanding of the timing and nature of early neurodegeneration in Alzheimer's disease (AD) and developing methods to measure this in vivo. Autosomal dominant familial Alzheimer's disease (FAD) provides the opportunity for investigation of presymptomatic change. We assessed early microstructural breakdown of cortical grey matter in FAD with diffusion-weighted MRI.

Diffusion MRI (dMRI) studies using the tensor model have identified abnormal white matter development associated with perinatal risk factors in preterm infants studied at term equivalent age (TEA). However, this model is an oversimplification of the underlying neuroanatomy. Fixel-based analysis (FBA) is a novel quantitative framework, which identifies microstructural and macrostructural changes in individual fibre populations within voxels containing crossing fibres. The aim of this study was to apply FBA to investigate the relationship between fixel-based measures of apparent fibre density (FD), fibre bundle cross-section (FC), and fibre density and cross-section (FDC) and perinatal risk factors in preterm infants at TEA. We studied 50 infants (28 male) born at 24.0-32.9 (median 30.4) weeks gestational age (GA) and imaged at 38.6-47.1 (median 42.1) weeks postmenstrual age (PMA). dMRI data were acquired in non-collinear directions with b-value 2500 s/mm2 on a 3 Tesla system sited on the neonatal intensive care unit. FBA was performed to assess the relationship between FD, FC, FDC and PMA at scan, GA at birth, days on mechanical ventilation, days on total parenteral nutrition (TPN), birthweight z-score, and sex. FBA reveals fibre population-specific alterations in FD, FC and FDC associated with clinical risk factors. FD was positively correlated with GA at birth and was negatively correlated with number of days requiring ventilation. FC was positively correlated with GA at birth, birthweight z-scores and was higher in males. FC was negatively correlated with number of days on ventilation and days on TPN. FDC was positively correlated with GA at birth and birthweight z-scores, negatively correlated with days on ventilation and days on TPN and higher in males. We demonstrate that these relationships are fibre-specific even within regions of crossing fibres. These results show that aberrant white matter development involves both microstructural changes and macrostructural alterations.