Group statistics of DTI fiber bundles using spatial functions of tensor measures

pre-printWe present a framework for hypothesis testing of differences between groups of DTI ber tracts. An anatomical, tract-oriented coordinate system provides a basis for estimating the distribution of diffusion properties. The parametrization of sampled, smooth functions is normalized across a population using DTI atlas building. Functional data analysis, an extension of multivariate statistics to continuous functions is applied to the problem of hypothesis testing and discrimination. B-spline models of fractional anisotropy (FA) and Frobenius norm measures are analyzed jointly. Plots of the discrimination direction provide a clinical interpretation of the group differences. The methodology is tested on a pediatric study of subjects aged one and two years

DTI denoising for data with low signal to noise ratios

Low signal to noise ratio (SNR) experiments in diffusion tensor imaging (DTI) give key information about tracking and anisotropy, e. g., by measurements with small voxel sizes or with high b values. However, due to the complicated and dominating impact of thermal noise such data are still seldom analysed. In this paper Monte Carlo simulations are presented which investigate the distributions of noise for different DTI variables in low SNR situations. Based on this study a strategy for the application of spatial smoothing is derived. Optimal prerequisites for spatial filters are unbiased, bell shaped distributions with uniform variance, but, only few variables have a statistics close to that. To construct a convenient filter a chain of nonlinear Gaussian filters is adapted to peculiarities of DTI and a bias correction is introduced. This edge preserving three dimensional filter is then validated via a quasi realistic model. Further, it is shown that for small sample sizes the filter is as effective as a maximum likelihood estimator and produces reliable results down to a local SNR of approximately 1. The filter is finally applied to very recent data with isotropic voxels of the size 1×1×1mm^3 which corresponds to a spatially mean SNR of 2.5. This application demonstrates the statistical robustness of the filter method. Though the Rician noise model is only approximately realized in the data, the gain of information by spatial smoothing is considerable

Spatial Smoothing for Diffusion Tensor Imaging with low Signal to Noise Ratios

Though low signal to noise ratio (SNR) experiments in DTI give key information about tracking and anisotropy, e.g. by measurements with very small voxel sizes, due to the complicated impact of thermal noise such experiments are up to now seldom analysed. In this paper Monte Carlo simulations are presented which investigate the random fields of noise for different DTI variables in low SNR situations. Based on this study a strategy for spatial smoothing, which demands essentially uniform noise, is derived. To construct a convenient filter the weights of the nonlinear Aurich chain are adapted to DTI. This edge preserving three dimensional filter is then validated in different variants via a quasi realistic model and is applied to very new data with isotropic voxels of the size 1x1x1 mm3 which correspond to a spatial mean SNR of approximately 3

Microstructural Correlates of Resilience against Major Depressive Disorder: Epigenetic Mechanisms?

Mental disorders are a major cause of long-term disability and are a direct cause of mortality, with approximately 800.000 individuals dying from suicide every year worldwide - a high proportion of them related to major depressive disorder (MDD)^1^. Healthy relatives of patients with major depressive disorder (MDD) are at risk to develop the disease. This higher vulnerability is associated with structural^2-4^ and functional brain changes^5^. However, we found using high angular resolution diffusion imaging (HARDI) with 61 diffusion directions that neuron tracts between frontal cortices and limbic as well as temporal and parietal brain regions are characterized by better diffusion coefficients in unaffected relatives (UHR), who managed to stay healthy, compared to healthy volunteers without any family history for a psychiatric disease (HC). Moreover, those UHR with stronger fibre connections better managed incidences of adversity in early life without later developing depression, while in HC axonal connections were found to be decreased when they had early-life adversity. Altogether these findings indicate the presence of stronger neural fibre connections in UHR, which seem to be associated with resilience against environmental stressors, which we suggest occur through epigenetic mechanisms

Neuroimaging of structural pathology and connectomics in traumatic brain injury: Toward personalized outcome prediction.

Recent contributions to the body of knowledge on traumatic brain injury (TBI) favor the view that multimodal neuroimaging using structural and functional magnetic resonance imaging (MRI and fMRI, respectively) as well as diffusion tensor imaging (DTI) has excellent potential to identify novel biomarkers and predictors of TBI outcome. This is particularly the case when such methods are appropriately combined with volumetric/morphometric analysis of brain structures and with the exploration of TBI-related changes in brain network properties at the level of the connectome. In this context, our present review summarizes recent developments on the roles of these two techniques in the search for novel structural neuroimaging biomarkers that have TBI outcome prognostication value. The themes being explored cover notable trends in this area of research, including (1) the role of advanced MRI processing methods in the analysis of structural pathology, (2) the use of brain connectomics and network analysis to identify outcome biomarkers, and (3) the application of multivariate statistics to predict outcome using neuroimaging metrics. The goal of the review is to draw the community's attention to these recent advances on TBI outcome prediction methods and to encourage the development of new methodologies whereby structural neuroimaging can be used to identify biomarkers of TBI outcome

Unsupervised White Matter Fiber Clustering and Tract Probability Map Generation: Applications of a Gaussian Process framework for White Matter Fibers

With the increasing importance of fiber tracking in diffusion tensor images for clinical needs, there has been a growing demand for an objective mathematical framework to perform quantitative analysis of white matter fiber bundles incorporating their underlying physical significance. This paper presents such a novel mathematical framework that facilitates mathematical operations between tracts using an inner product based on Gaussian processes, between fibers which span a metric space. This metric facilitates combination of fiber tracts, rendering operations like tract membership to a bundle or bundle similarity simple. Based on this framework, we have designed an automated unsupervised atlas-based clustering method that does not require manual initialization nor an a priori knowledge of the number of clusters. Quantitative analysis can now be performed on the clustered tract volumes across subjects thereby avoiding the need for point parametrization of these fibers, or the use of medial or envelope representations as in previous work. Experiments on synthetic data demonstrate the mathematical operations. Subsequently, the applicability of the unsupervised clustering framework has been demonstrated on a 21 subject dataset

Development of High Angular Resolution Diffusion Imaging Analysis Paradigms for the Investigation of Neuropathology

Diffusion weighted magnetic resonance imaging (DW-MRI), provides unique insight into the microstructure of neural white matter tissue, allowing researchers to more fully investigate white matter disorders. The abundance of clinical research projects incorporating DW-MRI into their acquisition protocols speaks to the value this information lends to the study of neurological disease. However, the most widespread DW-MRI technique, diffusion tensor imaging (DTI), possesses serious limitations which restrict its utility in regions of complex white matter. Fueled by advances in DW-MRI acquisition protocols and technologies, a group of exciting new DW-MRI models, developed to address these concerns, are now becoming available to clinical researchers. The emergence of these new imaging techniques, categorized as high angular resolution diffusion imaging (HARDI), has generated the need for sophisticated computational neuroanatomic techniques able to account for the high dimensionality and structure of HARDI data. The goal of this thesis is the development of such techniques utilizing prominent HARDI data models. Specifically, methodologies for spatial normalization, population atlas building and structural connectivity have been developed and validated. These methods form the core of a comprehensive analysis paradigm allowing the investigation of local white matter microarcitecture, as well as, systemic properties of neuronal connectivity. The application of this framework to the study of schizophrenia and the autism spectrum disorders demonstrate its sensitivity sublte differences in white matter organization, as well as, its applicability to large population DW-MRI studies