Whole brain voxel-wise analysis of single-subject serial DTI by permutation testing

Personalized microstructural evaluation using a Mahalanobis-distance based outlier detection strategy on epilepsy patients' DTI data - Theory, simulations and example cases

Quantitative MRI methods have recently gained extensive interest and are seeing substantial developments; however, their application in single patient vs control group comparisons is often limited by inherent statistical difficulties. One such application is detecting malformations of cortical development (MCDs) behind drug resistant epilepsies, a task that, especially when based solely on conventional MR images, may represent a serious challenge. We aimed to develop a novel straightforward voxel-wise evaluation method based on the Mahalanobis-distance, combining quantitative MRI data into a multidimensional parameter space and detecting lesion voxels as outliers. Simulations with standard multivariate Gaussian distribution and resampled DTI-eigenvalue data of 45 healthy control subjects determined the optimal critical value, cluster size threshold, and the expectable lesion detection performance through ROC-analyses. To reduce the effect of false positives emanating from registration artefacts and gyrification differences, an automatic classification method was applied, fine-tuned using a leave-one-out strategy based on diffusion and T1-weighted data of the controls. DWI processing, including thorough corrections and robust tensor fitting was performed with ExploreDTI, spatial coregistration was achieved with the DARTEL tools of SPM12. Additional to simulations, clusters of outlying diffusion profile, concordant with neuroradiological evaluation and independent calculations with the MAP07 toolbox were identified in 12 cases of a 13 patient example population with various types of MCDs. The multidimensional approach proved sufficiently sensitive in pinpointing regions of abnormal tissue microstructure using DTI data both in simulations and in the heterogeneous example population. Inherent limitations posed by registration artefacts, age-related differences, and the different or mixed pathologies limit the generalization of specificity estimation. Nevertheless, the proposed statistical method may aid the everyday examination of individual subjects, ever so more upon extending the framework with quantitative information from other modalities, e.g. susceptibility mapping, relaxometry, or perfusion.

A prospective microstructure imaging study in mixed-martial artists using geometric measures and diffusion tensor imaging: methods and findings

Although diffusion magnetic resonance imaging (dMRI) has been widely used to characterize the effects of repetitive mild traumatic brain injury (rmTBI), to date no studies have investigated how novel geometric models of microstructure relate to more typical diffusion tensor imaging (DTI) sequences. Moreover, few studies have evaluated the sensitivity of different registration pipelines (non-linear, linear and tract-based spatial statistics) for detecting dMRI abnormalities in clinical populations. Results from single-subject analyses in healthy controls (HC) indicated a strong negative relationship between fractional anisotropy (FA) and orientation dispersion index (ODI) in both white and gray matter. Equally important, only moderate relationships existed between all other estimates of free/intracellular water volume fractions and more traditional DTI metrics (FA, mean, axial and radial diffusivity). These findings suggest that geometric measures provide differential information about the cellular microstructure relative to traditional DTI measures. Results also suggest greater sensitivity for non-linear registration pipelines that maximize the anatomical information available in T1-weighted images. Clinically, rmTBI resulted in a pattern of decreased FA and increased ODI, largely overlapping in space, in conjunction with increased intracellular and free water fractions, highlighting the potential role of edema following repeated head trauma. In summary, current results suggest that geometric models of diffusion can provide relatively unique information regarding potential mechanisms of pathology that contribute to long-term neurological damage.

Improved spatial regression analysis of diffusion tensor imaging for lesion detection during longitudinal progression of multiple sclerosis in individual subjects

Subject-specific longitudinal DTI study is vital for investigation of pathological changes of lesions and disease evolution. Spatial Regression Analysis of Diffusion tensor imaging (SPREAD) is a non-parametric permutation-based statistical framework that combines spatial regression and resampling techniques to achieve effective detection of localized longitudinal diffusion changes within the whole brain at individual level without a priori hypotheses. However, boundary blurring and dislocation limit its sensitivity, especially towards detecting lesions of irregular shapes. In the present study, we propose an improved SPREAD (dubbed improved SPREAD, or iSPREAD) method by incorporating a three-dimensional (3D) nonlinear anisotropic diffusion filtering method, which provides edge-preserving image smoothing through a nonlinear scale space approach. The statistical inference based on iSPREAD was evaluated and compared with the original SPREAD method using both simulated and in vivo human brain data. Results demonstrated that the sensitivity and accuracy of the SPREAD method has been improved substantially by adapting nonlinear anisotropic filtering. iSPREAD identifies subject-specific longitudinal changes in the brain with improved sensitivity, accuracy, and enhanced statistical power, especially when the spatial correlation is heterogeneous among neighboring image pixels in DTI.

Spatial regression analysis of serial DTI for subject-specific longitudinal changes of neurodegenerative disease

Quantitative measurement of localized longitudinal changes in brain abnormalities at an individual level may offer critical information for disease diagnosis and treatment. The voxel-wise permutation-based method SPREAD/iSPREAD, which combines resampling and spatial regression of neighboring voxels, provides an effective and robust method for detecting subject-specific longitudinal changes within the whole brain, especially for longitudinal studies with a limited number of scans. As an extension of SPREAD/iSPREAD, we present a general method that facilitates analysis of serial Diffusion Tensor Imaging (DTI) measurements (with more than two time points) for testing localized changes in longitudinal studies. Two types of voxel-level test statistics (model-free test statistics, which measure intra-subject variability across time, and test statistics based on general linear model that incorporate specific lesion evolution models) were estimated and tested against the null hypothesis among groups of DTI data across time. The implementation and utility of the proposed statistical method were demonstrated by both Monte Carlo simulations and applications on clinical DTI data from human brain in vivo. By a design of test statistics based on the disease progression model, it was possible to apportion the true significant voxels attributed to the disease progression and those caused by underlying anatomical differences that cannot be explained by the model, which led to improvement in false positive (FP) control in the results. Extension of the proposed method to include other diseases or drug effect models, as well as the feasibility of global statistics, was discussed. The proposed statistical method can be extended to a broad spectrum of longitudinal studies with carefully designed test statistics, which helps to detect localized changes at the individual level.

•

A nonparametric method for detecting subject-specific localized changes in longitudinal DTI

•

Obtain sufficient statistical power even with limited scans available

•

Various voxel-level test statistics for hypothesis tests among groups across time

•

Excellent at controlling false positive ratio with the model-based test statistics

Abnormal white matter integrity in antipsychotic-naïve first-episode psychosis patients assessed by a DTI principal component analysis

INTRODUCTION

Diffusion tensor imaging (DTI) studies in patients with schizophrenia have shown abnormalities in the microstructure of white matter tracts. Specifically, reduced fractional anisotropy (FA) has been described across multiple white matter tracts, in studies that have mainly included patients treated with antipsychotic medications.

OBJECTIVE

To compare FA in antipsychotic-naïve patients experiencing a first episode of psychosis (FEP) to FA in healthy controls to demonstrate that the variance of FA can be grouped, in a coincidental manner, in four predetermined factors in accordance with a theoretical partition of the white matter tracts, using a principal components analysis (PCA).

METHODS

Thirty-five antipsychotic-naïve FEP patients and 35 age- and gender-matched healthy controls underwent DTI at 3T. Analysis was performed using a tract-based spatial statistics (TBSS) method and exploratory PCA.

RESULTS

DTI analysis showed extensive FA reduction in white matter tracts in FEP patients compared with the control group. The PCA grouped the white matter tracts into four factors explaining 66% of the total variance. Comparison of the FA values within each factor highlighted the differences between FEP patients and controls.

DISCUSSION

Our study confirms extensive white matter tracts anomalies in patients with schizophrenia, more specifically, in drug-naïve FEP patients. The results also indicate that a small number of white matter tracts share common FA anomalies that relate to deficit symptoms in FEP patients. Our study adds to a growing body of literature emphasizing the need for treatments targeting white matter function and structure in FEP patients.

Harnessing graphics processing units for improved neuroimaging statistics

Simple models and algorithms based on restrictive assumptions are often used in the field of neuroimaging for studies involving functional magnetic resonance imaging, voxel based morphometry, and diffusion tensor imaging. Nonparametric statistical methods or flexible Bayesian models can be applied rather easily to yield more trustworthy results. The spatial normalization step required for multisubject studies can also be improved by taking advantage of more robust algorithms for image registration. A common drawback of algorithms based on weaker assumptions, however, is the increase in computational complexity. In this short overview, we will therefore present some examples of how inexpensive PC graphics hardware, normally used for demanding computer games, can be used to enable practical use of more realistic models and accurate algorithms, such that the outcome of neuroimaging studies really can be trusted.

SPatial REgression Analysis of Diffusion tensor imaging (SPREAD) for longitudinal progression of neurodegenerative disease in individual subjects

OBJECTIVES

To develop a novel statistical method for analysis of longitudinal DTI data in individual subjects.

MATERIALS AND METHODS

The proposed SPatial REgression Analysis of Diffusion tensor imaging (SPREAD) method incorporates a spatial regression fitting of DTI data among neighboring voxels and a resampling method among data at different times. Both numerical simulations and real DTI data from healthy volunteers and multiple sclerosis (MS) patients were used in the study to evaluate this method.

RESULTS

Statistical inference based on SPREAD was shown to perform well through both group comparisons among simulated DTI data of individuals (especially when the group size is smaller than 5) and longitudinal comparisons of human DTI data within the same individual.

CONCLUSIONS

When pathological changes of neurodegenerative diseases are heterogeneous in a population, SPREAD provides a unique way to assess abnormality during disease progression at the individual level. Consequently, it has the potential to shed light on how the brain has changed as a result of disease or injury.

Longitudinal change detection: inference on the diffusion tensor along white matter pathways

Diffusion weighted magnetic resonance imaging (DW-MRI) makes it possible to probe brain connections in vivo. This paper presents a change detection framework that relies on white matter pathways with application to neuromyelitis optica (NMO). The objective is to detect local or global fiber diffusion property modifications between two longitudinal DW-MRI acquisitions of a patient. To this end, we develop two frameworks based on statistical tests on tensor eigenvalues to detect local or global changes along the white matter pathways: a pointwise test that compares tensor populations extracted in bundles cross sections and a fiberwise test that compares paired tensors along all the fiber bundles. Experiments on both synthetic and real data highlight the benefit of considering fiber based statistical tests compared to standard voxelwise strategies.

Logical circularity in voxel-based analysis: normalization strategy may induce statistical bias

Recent discussions within the neuroimaging community have highlighted the problematic presence of selection bias in experimental design. Although initially centering on the selection of voxels during the course of fMRI studies, we demonstrate how this bias can potentially corrupt voxel-based analyses. For such studies, template-based registration plays a critical role in which a representative template serves as the normalized space for group alignment. A standard approach maps each subject's image to a representative template before performing statistical comparisons between different groups. We analytically demonstrate that in these scenarios the popular sum of squared difference (SSD) intensity metric, implicitly surrogating as a quantification of anatomical alignment, instead explicitly maximizes effect size--an experimental design flaw referred to as "circularity bias." We illustrate how this selection bias varies in strength with the similarity metric used during registration under the hypothesis that while SSD-related metrics, such as Demons, will manifest similar effects, other metrics which are not formulated based on absolute intensity differences will produce less of an effect. Consequently, given the variability in voxel-based analysis outcomes with similarity metric choice, we caution researchers specifically in the use of SSD and SSD-related measures where normalization and statistical analysis involve the same image set. Instead, we advocate a more cautious approach where normalization of the individual subject images to the reference space occurs through corresponding image sets which are independent of statistical testing. Alternatively, one can use similarity terms that are less sensitive to this bias.

Longitudinal change detection in diffusion MRI using multivariate statistical testing on tensors

This paper presents a longitudinal change detection framework for detecting relevant modifications in diffusion MRI, with application to neuromyelitis optica (NMO) and multiple sclerosis (MS). The core problem is to identify image regions that are significantly different between two scans. The proposed method is based on multivariate statistical testing which was initially introduced for tensor population comparison. We use this method in the context of longitudinal change detection by considering several strategies to build sets of tensors characterizing the variability of each voxel. These strategies make use of the variability existing in the diffusion weighted images (thanks to a bootstrap procedure), or in the spatial neighborhood of the considered voxel, or a combination of both. Results on synthetic evolutions and on real data are presented. Interestingly, experiments on NMO patients highlight the ability of the proposed approach to detect changes in the normal-appearing white matter (according to conventional MRI) that are related with physical status outcome. Experiments on MS patients highlight the ability of the proposed approach to detect changes in evolving and non-evolving lesions (according to conventional MRI). These findings might open promising prospects for the follow-up of NMO and MS pathologies.

Volumetric navigators for real-time motion correction in diffusion tensor imaging

Prospective motion correction methods using an optical system, diffusion-weighted prospective acquisition correction, or a free induction decay navigator have recently been applied to correct for motion in diffusion tensor imaging. These methods have some limitations and drawbacks. This article describes a novel technique using a three-dimensional-echo planar imaging navigator, of which the contrast is independent of the b-value, to perform prospective motion correction in diffusion weighted images, without having to reacquire volumes during which motion occurred, unless motion exceeded some preset thresholds. Water phantom and human brain data were acquired using the standard and navigated diffusion sequences, and the mean and whole brain histogram of the fractional anisotropy and mean diffusivity were analyzed. Our results show that adding the navigator does not influence the diffusion sequence. With head motion, the whole brain histogram-fractional anisotropy shows a shift toward lower anisotropy with a significant decrease in both the mean fractional anisotropy and the fractional anisotropy histogram peak location (P<0.01), whereas the whole brain histogram-mean diffusivity shows a shift toward higher diffusivity with a significant increase in the mean diffusivity (P<0.01), even after retrospective motion correction. These changes in the mean and the shape of the histograms are recovered substantially in the prospective motion corrected data acquired using the navigated sequence.

Longitudinal change detection: inference on the diffusion tensor along white-matter pathways

Diffusion tensor magnetic resonance imaging (DT-MRI) tractography allows to probe brain connections in vivo. This paper presents a change detection framework that relies on white-matter pathways with application to neuromyelitis optica (NMO). The objective is to detect global or local fiber diffusion property modifications between two longitudinal DT-MRI acquisitions of a patient. To this end, estimation and testing tools on tensors along the white-matter pathways are considered. Two tests are implemented: a pointwise test that compares at each sampling point of the fiber bundle the tensor populations of the two exams in the cross section of the bundle and a fiberwise test that compares paired tensors along all the fiber bundle. Experiments on both synthetic and real data highlight the benefit of considering fiber based statistical tests compared to the standard voxelwise strategy.

Diffusion MR imaging: basic principles

From their origin as simple techniques primarily used for detecting acute cerebral ischemia, diffusion MR imaging techniques have rapidly evolved into a versatile set of tools that provide the only noninvasive means of characterizing brain microstructure and connectivity, becoming a mainstay of both clinical and investigational brain MR imaging. In this article, the basic principles required for understanding diffusion MR imaging techniques are reviewed with clinical neuroradiologists in mind.

Change detection in diffusion MRI using multivariate statistical testing on tensors

This paper presents a longitudinal change detection framework for detecting relevant modifications in diffusion MRI, with application to Multiple Sclerosis (MS). The proposed method is based on multivariate statistical testings which were initially introduced for tensor population comparison. We use these methods in the context of longitudinal change detection by considering several strategies to build sets of tensors characterizing the variability of each voxel. These testing tools have been considered either for the comparison of tensor eigenvalues or eigenvectors, thus enabling to differentiate orientation and diffusivity changes. Results on simulated MS lesion evolutions and on real data are presented. Interestingly, experiments on an MS patient highlight the ability of the proposed approach to detect changes in non evolving lesions (according to conventional MRI) and around lesions (in the normal appearing white matter), which might open promising perspectives for the follow-up of the MS pathology.

A statistical method (cross-validation) for bone loss region detection after spaceflight

Astronauts experience bone loss after the long spaceflight missions. Identifying specific regions that undergo the greatest losses (e.g. the proximal femur) could reveal information about the processes of bone loss in disuse and disease. Methods for detecting such regions, however, remains an open problem. This paper focuses on statistical methods to detect such regions. We perform statistical parametric mapping to get t-maps of changes in images, and propose a new cross-validation method to select an optimum suprathreshold for forming clusters of pixels. Once these candidate clusters are formed, we use permutation testing of longitudinal labels to derive significant changes.

Statistical detection of longitudinal changes between apparent diffusion coefficient images: application to multiple sclerosis

The automatic analysis of longitudinal changes between Diffusion Tensor Imaging (DTI) acquisitions is a promising tool for monitoring disease evolution. However, few works address this issue and existing methods are generally limited to the detection of changes between scalar images characterizing diffusion properties, such as Fractional Anisotropy or Mean Diffusivity, while richer information can be exploited from the whole set of Apparent Diffusion Coefficient (ADC) images that can be derived from a DTI acquisition. In this paper, we present a general framework for detecting changes between two sets of ADC images and we investigate the performance of four statistical tests. Results are presented on both simulated and real data in the context of the follow-up of multiple sclerosis lesion evolution.

VOXEL-WISE GROUP ANALYSIS OF DTI

Diffusion tensor MRI (DTI) is now a widely used modality to investigate the fiber tissues in vivo, especially the white matter in brain. An automatic pipeline is described in this paper to conduct a localized voxel-wise multiple-subject group comparison study of DTI. The pipeline consists of 3 steps: 1) Preprocessing, including image format converting, image quality check, eddy-current and motion artifact correction, skull stripping and tensor image estimation, 2) study-specific unbiased DTI atlas computation via affine followed by fluid nonlinear registration and warping of all individual DTI images into the common atlas space to achieve voxel-wise correspondence, 3) voxelwise statistical analysis via heterogeneous linear regression and wild bootstrap technique for correcting for multiple comparisons. This pipeline was applied to process data from a fitness and aging study and preliminary results are presented. The results show that this fully automatic pipeline is suitable for voxel-wise group DTI analysis.

Diffusion tensor MR imaging and fiber tractography: technical considerations

This second article of the 2-part review builds on the theoretic background provided by the first article to cover the major technical factors that affect image quality in diffusion imaging, including the acquisition sequence, magnet field strength, gradient amplitude, and slew rate as well as multichannel radio-frequency coils and parallel imaging. The sources of many common diffusion image artifacts are also explored in detail. The emphasis is on optimizing these technical factors for state-of-the-art diffusion-weighted imaging and diffusion tensor imaging (DTI) based on the best available evidence in the literature. An overview of current methods for quantitative analysis of DTI data and fiber tractography in clinical research is also provided.