Faculty Bibliography

Tissue characterization in brain tumors and, in particular, in high-grade gliomas is challenging as a result of the co-existence of several intra-tumoral tissue types within the same region and the high spatial heterogeneity. This study presents a method for the detection of the relevant tumor substructures (i.e. viable tumor, necrosis and edema), which could be of added value for the diagnosis, treatment planning and follow-up of individual patients. Twenty-four patients with glioma [10 low-grade gliomas (LGGs), 14 high-grade gliomas (HGGs)] underwent a multi-parametric MRI (MP-MRI) scheme, including conventional MRI (cMRI), perfusion-weighted imaging (PWI), diffusion kurtosis imaging (DKI) and short-TE (1)H MRSI. MP-MRI parameters were derived: T2, T1 + contrast, fluid-attenuated inversion recovery (FLAIR), relative cerebral blood volume (rCBV), mean diffusivity (MD), fractional anisotropy (FA), mean kurtosis (MK) and the principal metabolites lipids (Lip), lactate (Lac), N-acetyl-aspartate (NAA), total choline (Cho), etc. Hierarchical non-negative matrix factorization (hNMF) was applied to the MP-MRI parameters, providing tissue characterization on a patient-by-patient and voxel-by-voxel basis. Tissue-specific patterns were obtained and the spatial distribution of each tissue type was visualized by means of abundance maps. Dice scores were calculated by comparing tissue segmentation derived from hNMF with the manual segmentation by a radiologist. Correlation coefficients were calculated between each pathologic tissue source and the average feature vector within the corresponding tissue region. For the patients with HGG, mean Dice scores of 78%, 85% and 83% were obtained for viable tumor, the tumor core and the complete tumor region. The mean correlation coefficients were 0.91 for tumor, 0.97 for necrosis and 0.96 for edema. For the patients with LGG, a mean Dice score of 85% and mean correlation coefficient of 0.95 were found for the tumor region. hNMF was also applied to reduced MRI datasets, showing the added value of individual MRI modalities.

The purpose of this explorative study was to investigate whether diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) parameter changes are reliable measures of white matter integrity changes in Alzheimer's disease (AD) patients using a whole brain voxel-based analysis (VBA). Therefore, age- and gender-matched patients with mild cognitive impairment (MCI) due to AD (n = 18), dementia due to AD (n = 19), and age-matched cognitively healthy controls (n = 14) were prospectively included. The magnetic resonance imaging protocol included routine structural brain imaging and DKI. Datasets were transformed to a population-specific atlas space. Groups were compared using VBA. Differences in diffusion and mean kurtosis measures between MCI and AD patients and controls were shown, and were mainly found in the splenium of the corpus callosum and the corona radiata. Hence, DTI and DKI parameter changes are suggestive of white matter changes in AD.

PURPOSE/OBJECTIVE:Diffusion-weighted magnetic resonance imaging suffers from physiological noise, such as artifacts caused by motion or system instabilities. Therefore, there is a need for robust diffusion parameter estimation techniques. In the past, several techniques have been proposed, including RESTORE and iRESTORE (Chang et al. Magn Reson Med 2005; 53:1088-1095; Chang et al. Magn Reson Med 2012; 68:1654-1663). However, these techniques are based on nonlinear estimators and are consequently computationally intensive. METHOD/METHODS:In this work, we present a new, robust, iteratively reweighted linear least squares (IRLLS) estimator. IRLLS performs a voxel-wise identification of outliers in diffusion-weighted magnetic resonance images, where it exploits the natural skewness of the data distribution to become more sensitive to both signal hyperintensities and signal dropouts. RESULTS:Both simulations and real data experiments were conducted to compare IRLLS with other state-of-the-art techniques. While IRLLS showed no significant loss in accuracy or precision, it proved to be substantially faster than both RESTORE and iRESTORE. In addition, IRLLS proved to be even more robust when considering the overestimation of the noise level or when the signal-to-noise ratio is low. CONCLUSION/CONCLUSIONS:The substantially shortened calculation time in combination with the increased robustness and accuracy, make IRLLS a practical and reliable alternative to current state-of-the-art techniques for the robust estimation of diffusion-weighted magnetic resonance parameters.

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.

BACKGROUND:We assessed the diagnostic accuracy of diffusion kurtosis imaging (DKI), dynamic susceptibility-weighted contrast-enhanced (DSC) MRI, and short echo time chemical shift imaging (CSI) for grading gliomas. METHODS:In this prospective study, 35 patients with cerebral gliomas underwent DKI, DSC, and CSI on a 3 T MR scanner. Diffusion parameters were mean diffusivity (MD), fractional anisotropy, and mean kurtosis (MK). Perfusion parameters were mean relative regional cerebral blood volume (rrCBV), mean relative regional cerebral blood flow (rrCBF), mean transit time, and relative decrease ratio (rDR). The diffusion and perfusion parameters along with 12 CSI metabolite ratios were compared among 22 high-grade gliomas and 14 low-grade gliomas (Mann-Whitney U-test, P < .05). Classification accuracy was determined with a linear discriminant analysis for each MR modality independently. Furthermore, the performance of a multimodal analysis is reported, using a decision-tree rule combining the statistically significant DKI, DSC-MRI, and CSI parameters with the lowest P-value. The proposed classifiers were validated on a set of subsequently acquired data from 19 clinical patients. RESULTS:Statistically significant differences among tumor grades were shown for MK, MD, mean rrCBV, mean rrCBF, rDR, lipids over total choline, lipids over creatine, sum of myo-inositol, and sum of creatine. DSC-MRI proved to be the modality with the best performance when comparing modalities individually, while the multimodal decision tree proved to be most accurate in predicting tumor grade, with a performance of 86%. CONCLUSIONS:Combining information from DKI, DSC-MRI, and CSI increases diagnostic accuracy to differentiate low- from high-grade gliomas, possibly providing diagnosis for the individual patient.

OBJECTIVE:To evaluate the effect of monophasic combined oral contraceptive pill (COCP) and menstrual cycle phase in healthy young women on white matter (WM) organization using diffusion tensor imaging (DTI). METHODS:Thirty young women were included in the study; 15 women used COCP and 15 women had a natural cycle. All subjects underwent DTI magnetic resonance imaging during the follicular and luteal phase of their cycle, or in different COCP cycle phases. DTI parameters were obtained in different WM structures by performing diffusion tensor fibre tractography. Fractional anisotropy and mean diffusivity were calculated for different WM structures. Hormonal plasma concentrations were measured in peripheral venous blood samples and correlated with the DTI findings. RESULTS:We found a significant difference in mean diffusivity in the fornix between the COCP and the natural cycle group. Mean diffusivity values in the fornix were negatively correlated with luteinizing hormone and estradiol blood concentrations. CONCLUSION/CONCLUSIONS:An important part in the limbic system, the fornix, regulates emotional processes. Differences in diffusion parameters in the fornix may contribute to behavioural alternations related to COCP use. This finding also suggests that the use of oral contraceptives needs to be taken into account when designing DTI group studies.

Rodent models of Huntington disease (HD) are valuable tools for investigating HD pathophysiology and evaluating new therapeutic approaches. Non-invasive characterization of HD-related phenotype changes is important for monitoring progression of pathological processes and possible effects of interventions. The first transgenic rat model for HD exhibits progressive late-onset affective, cognitive, and motor impairments, as well as neuropathological features reflecting observations from HD patients. In this report, we contribute to the anatomical phenotyping of this model by comparing high-resolution ex vivo DTI measurements obtained in aged transgenic HD rats and wild-type controls. By region of interest analysis supplemented by voxel-based statistics, we find little evidence of atrophy in basal ganglia regions, but demonstrate altered DTI measurements in the dorsal and ventral striatum, globus pallidus, entopeduncular nucleus, substantia nigra, and hippocampus. These changes are largely compatible with DTI findings in preclinical and clinical HD patients. We confirm earlier reports that HD rats express a moderate neuropathological phenotype, and provide evidence of altered DTI measures in specific HD-related brain regions, in the absence of pronounced morphometric changes.

PURPOSE/OBJECTIVE:Linear least squares estimators are widely used in diffusion MRI for the estimation of diffusion parameters. Although adding proper weights is necessary to increase the precision of these linear estimators, there is no consensus on how to practically define them. In this study, the impact of the commonly used weighting strategies on the accuracy and precision of linear diffusion parameter estimators is evaluated and compared with the nonlinear least squares estimation approach. METHODS:Simulation and real data experiments were done to study the performance of the weighted linear least squares estimators with weights defined by (a) the squares of the respective noisy diffusion-weighted signals; and (b) the squares of the predicted signals, which are reconstructed from a previous estimate of the diffusion model parameters. RESULTS:The negative effect of weighting strategy (a) on the accuracy of the estimator was surprisingly high. Multi-step weighting strategies yield better performance and, in some cases, even outperformed the nonlinear least squares estimator. CONCLUSION/CONCLUSIONS:If proper weighting strategies are applied, the weighted linear least squares approach shows high performance characteristics in terms of accuracy/precision and may even be preferred over nonlinear estimation methods.

PURPOSE/OBJECTIVE:Amyloid deposition in the brain is considered an initial event in the progression of Alzheimer's disease. We hypothesized that the presence of amyloid plaques in the brain of APP/presenilin 1 mice leads to higher diffusion kurtosis measures due to increased microstructural complexity. As such, our purpose was to provide an in vivo proof of principle for detection of amyloidosis by diffusion kurtosis imaging (DKI). METHODS:APPKM670/671NL /presenilin 1 L166P mice (n = 5) and wild-type littermates (n = 5) underwent DKI at the age of 16 months. Averaged diffusion and diffusion kurtosis parameters were obtained for multiple regions (hippocampus-cortex-thalamus-cerebellum). After DKI, mice were sacrificed for amyloid staining. RESULTS:Histograms of the frequency distribution of the DKI parameters tended to shift to higher values. After normalization of absolute values to the cerebellum, a nearly plaque-free region, mean, radial, and axial diffusion kurtosis were significantly higher in APP/presenilin 1 mice as compared to wild-type in the cortex and thalamus, regions demonstrating substantial amyloid staining. CONCLUSION/CONCLUSIONS:The current study, although small-scale, suggests increased DKI metrics, in the absence of alterations in diffusion tensor imaging metrics in the cortex and thalamus of APP/presenilin 1 mice with established amyloidosis. These results warrant further investigations on the potential of DKI as a sensitive marker for Alzheimer's disease.

During the last decade, many approaches have been proposed for improving the estimation of diffusion measures. These techniques have already shown an increase in accuracy based on theoretical considerations, such as incorporating prior knowledge of the data distribution. The increased accuracy of diffusion metric estimators is typically observed in well-defined simulations, where the assumptions regarding properties of the data distribution are known to be valid. In practice, however, correcting for subject motion and geometric eddy current deformations alters the data distribution tremendously such that it can no longer be expressed in a closed form. The image processing steps that precede the model fitting will render several assumptions on the data distribution invalid, potentially nullifying the benefit of applying more advanced diffusion estimators. In this work, we present a generic diffusion model fitting framework that considers some statistics of diffusion MRI data. A central role in the framework is played by the conditional least squares estimator. We demonstrate that the accuracy of that particular estimator can generally be preserved, regardless the applied preprocessing steps, if the noise parameter is known a priori. To fulfill that condition, we also propose an approach for the estimation of spatially varying noise levels.