Measures for validation of DTI tractography

pre-printThe evaluation of analysis methods for diffusion tensor imaging (DTI) remains challenging due to the lack of gold standards and validation frameworks. Significant work remains in developing metrics for comparing fiber bundles generated from streamline tractography. We propose a set of volumetric and tract oriented measures for evaluating tract differences. The different methods developed for this assessment work are: an overlap measurement, a point cloud distance and a quantification of the diffusion properties at similar locations between fiber bundles. The application of the measures in this paper is a comparison of atlas generated tractography to tractography generated in individual images. For the validation we used a database of 37 subject DTIs, and applied the measurements on five specific fiber bundles: uncinate, cingulum (left and right for both bundles) and genu. Each measurments is interesting for specific use: the overlap measure presents a simple and comprehensive metric but is sensitive to partial voluming and does not give consistent values depending on the bundle geometry. The point cloud distance associated with a quantile interpretation of the distribution gives a good intuition of how close and similar the bundles are. Finally, the functional difference is useful for a comparison of the diffusion properties since it is the focus of many DTI analysis to compare scalar invariants. The comparison demonstrated reasonable similarity of results. The tract difference measures are also applicable to comparison of tractography algorithms, quality control, reproducibility studies, and other validation problems

Feasibility of diffusion and probabilistic white matter analysis in patients implanted with a deep brain stimulator.

Deep brain stimulation (DBS) for Parkinson\u27s disease (PD) is an established advanced therapy that produces therapeutic effects through high frequency stimulation. Although this therapeutic option leads to improved clinical outcomes, the mechanisms of the underlying efficacy of this treatment are not well understood. Therefore, investigation of DBS and its postoperative effects on brain architecture is of great interest. Diffusion weighted imaging (DWI) is an advanced imaging technique, which has the ability to estimate the structure of white matter fibers; however, clinical application of DWI after DBS implantation is challenging due to the strong susceptibility artifacts caused by implanted devices. This study aims to evaluate the feasibility of generating meaningful white matter reconstructions after DBS implantation; and to subsequently quantify the degree to which these tracts are affected by post-operative device-related artifacts. DWI was safely performed before and after implanting electrodes for DBS in 9 PD patients. Differences within each subject between pre- and post-implantation FA, MD, and RD values for 123 regions of interest (ROIs) were calculated. While differences were noted globally, they were larger in regions directly affected by the artifact. White matter tracts were generated from each ROI with probabilistic tractography, revealing significant differences in the reconstruction of several white matter structures after DBS. Tracts pertinent to PD, such as regions of the substantia nigra and nigrostriatal tracts, were largely unaffected. The aim of this study was to demonstrate the feasibility and clinical applicability of acquiring and processing DWI post-operatively in PD patients after DBS implantation. The presence of global differences provides an impetus for acquiring DWI shortly after implantation to establish a new baseline against which longitudinal changes in brain connectivity in DBS patients can be compared. Understanding that post-operative fiber tracking in patients is feasible on a clinically-relevant scale has significant implications for increasing our current understanding of the pathophysiology of movement disorders, and may provide insights into better defining the pathophysiology and therapeutic effects of DBS

Nonlinear tube-fitting for the analysis of anatomical and functional structures

We are concerned with the estimation of the exterior surface and interior summaries of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT) and to a white-matter tract image produced using diffusion tensor imaging (DTI).Comment: Published in at http://dx.doi.org/10.1214/10-AOAS384 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Fuzzy Fibers: Uncertainty in dMRI Tractography

Fiber tracking based on diffusion weighted Magnetic Resonance Imaging (dMRI) allows for noninvasive reconstruction of fiber bundles in the human brain. In this chapter, we discuss sources of error and uncertainty in this technique, and review strategies that afford a more reliable interpretation of the results. This includes methods for computing and rendering probabilistic tractograms, which estimate precision in the face of measurement noise and artifacts. However, we also address aspects that have received less attention so far, such as model selection, partial voluming, and the impact of parameters, both in preprocessing and in fiber tracking itself. We conclude by giving impulses for future research

Disconnection of network hubs and cognitive impairment after traumatic brain injury.

Traumatic brain injury affects brain connectivity by producing traumatic axonal injury. This disrupts the function of large-scale networks that support cognition. The best way to describe this relationship is unclear, but one elegant approach is to view networks as graphs. Brain regions become nodes in the graph, and white matter tracts the connections. The overall effect of an injury can then be estimated by calculating graph metrics of network structure and function. Here we test which graph metrics best predict the presence of traumatic axonal injury, as well as which are most highly associated with cognitive impairment. A comprehensive range of graph metrics was calculated from structural connectivity measures for 52 patients with traumatic brain injury, 21 of whom had microbleed evidence of traumatic axonal injury, and 25 age-matched controls. White matter connections between 165 grey matter brain regions were defined using tractography, and structural connectivity matrices calculated from skeletonized diffusion tensor imaging data. This technique estimates injury at the centre of tract, but is insensitive to damage at tract edges. Graph metrics were calculated from the resulting connectivity matrices and machine-learning techniques used to select the metrics that best predicted the presence of traumatic brain injury. In addition, we used regularization and variable selection via the elastic net to predict patient behaviour on tests of information processing speed, executive function and associative memory. Support vector machines trained with graph metrics of white matter connectivity matrices from the microbleed group were able to identify patients with a history of traumatic brain injury with 93.4% accuracy, a result robust to different ways of sampling the data. Graph metrics were significantly associated with cognitive performance: information processing speed (R(2) = 0.64), executive function (R(2) = 0.56) and associative memory (R(2) = 0.25). These results were then replicated in a separate group of patients without microbleeds. The most influential graph metrics were betweenness centrality and eigenvector centrality, which provide measures of the extent to which a given brain region connects other regions in the network. Reductions in betweenness centrality and eigenvector centrality were particularly evident within hub regions including the cingulate cortex and caudate. Our results demonstrate that betweenness centrality and eigenvector centrality are reduced within network hubs, due to the impact of traumatic axonal injury on network connections. The dominance of betweenness centrality and eigenvector centrality suggests that cognitive impairment after traumatic brain injury results from the disconnection of network hubs by traumatic axonal injury

Vessel tractography using an intensity based tensor model with branch detection

In this paper, we present a tubular structure seg- mentation method that utilizes a second order tensor constructed from directional intensity measurements, which is inspired from diffusion tensor image (DTI) modeling. The constructed anisotropic tensor which is fit inside a vessel drives the segmen- tation analogously to a tractography approach in DTI. Our model is initialized at a single seed point and is capable of capturing whole vessel trees by an automatic branch detection algorithm developed in the same framework. The centerline of the vessel as well as its thickness is extracted. Performance results within the Rotterdam Coronary Artery Algorithm Evaluation framework are provided for comparison with existing techniques. 96.4% average overlap with ground truth delineated by experts is obtained in addition to other measures reported in the paper. Moreover, we demonstrate further quantitative results over synthetic vascular datasets, and we provide quantitative experiments for branch detection on patient Computed Tomography Angiography (CTA) volumes, as well as qualitative evaluations on the same CTA datasets, from visual scores by a cardiologist expert

Age related diffusion and tractography changes in typically developing pediatric cervical and thoracic spinal cord

Background and objective: Diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) are two techniques that can measure white matter integrity of the spinal cord. Recently, DTI indices have been shown to change with age. The purpose of this study is (a) to evaluate the maturational states of the entire pediatric spinal cord using DTI and DTT indices including fractional anisotropy (FA), mean diffusivity (MD), mean length of white matter fiber tracts and tract density and (b) to analyze the DTI and DTT parameters along the entire spinal cord as a function of spinal cord levels and age. Method: A total of 23 typically developing (TD) pediatric subjects ranging in age from 6 to 16 years old (11.94 ± 3.26 (mean ± standard deviation), 13 females and 10 males) were recruited, and scanned using 3.0 T MR scanner. Reduced FOV diffusion tensor images were acquired axially in the same anatomical location prescribed for the T2-weighted images to cover the entire spinal cord (C1-mid L1 levels). To mitigate motion induced artifacts, diffusion directional images were aligned with the reference image (b0) using a rigid body registration algorithm performed by in-house software developed in Matlab (MathWorks, Natick, Massachusetts). Diffusion tensor maps (FA and MD) and streamline deterministic tractography were then generated from the motion corrected DTI dataset. DTI and DTT parameters were calculated by using ROIs drawn to encapsulate the whole cord along the entire spinal cord by an independent board certified neuroradiologist. These indices then were compared between two age groups (age group A = 6–11 years (n = 11) and age group B = 12–16 years (n = 12)) based on similar standards and age definitions used for reporting spinal cord injury in the pediatric population. Standard least squared linear regression based on a restricted maximum likelihood (REML) method was used to evaluate the relationship between age and DTI and DTT parameters. Results: An increase in FA (group A = 0.42 ± 0.097, group B = 0.49 ± 0.116), white matter tract density (group A = 368.01 ± 236.88, group B = 440.13 ± 245.24) and mean length of fiber tracts (group A = 48.16 ± 20.48 mm, group B = 60.28 ± 23.87 mm) and a decrease in MD (group A = 1.06 ± 0.23 × 10−3 mm2/s, group B = 0.82 ± 0.24 × 10−3 mm2/s) were observed with age along the entire spinal cord. Statistically significant increases have been shown in FA (p = 0.004, R2 = 0.57), tract density (p = 0.0004, R2 = 0.58), mean length of fiber tracts (p \u3c 0.001, R2 = 0.5) and a significant decrease has been shown in MD (p = 0.002, R2 = 0.59) between group A and group B. Also, it has been shown DTI and DTT parameters vary along the spinal cord as a function of intervertebral disk and mid-vertebral body level. Conclusion: This study provides an initial understanding of age related changes of DTI values as well as DTT metrics of the spinal cord. The results show significant differences in DTI and DTT parameters which may result from decreasing water content, myelination of fiber tracts, and the thickening diameter of fiber tracts during the maturation process. Consequently, when quantitative DTI and DTT of the spinal cord is undertaken in the pediatric population an age and level matched normative dataset should be used to accurately interpret the quantitative results. © 201

White matter differences between healthy young ApoE4 carriers and non-carriers identified with tractography and support vector machines.

The apolipoprotein E4 (ApoE4) is an established risk factor for Alzheimer's disease (AD). Previous work has shown that this allele is associated with functional (fMRI) changes as well structural grey matter (GM) changes in healthy young, middle-aged and older subjects. Here, we assess the diffusion characteristics and the white matter (WM) tracts of healthy young (20-38 years) ApoE4 carriers and non-carriers. No significant differences in diffusion indices were found between young carriers (ApoE4+) and non-carriers (ApoE4-). There were also no significant differences between the groups in terms of normalised GM or WM volume. A feature selection algorithm (ReliefF) was used to select the most salient voxels from the diffusion data for subsequent classification with support vector machines (SVMs). SVMs were capable of classifying ApoE4 carrier and non-carrier groups with an extremely high level of accuracy. The top 500 voxels selected by ReliefF were then used as seeds for tractography which identified a WM network that included regions of the parietal lobe, the cingulum bundle and the dorsolateral frontal lobe. There was a non-significant decrease in volume of this WM network in the ApoE4 carrier group. Our results indicate that there are subtle WM differences between healthy young ApoE4 carriers and non-carriers and that the WM network identified may be particularly vulnerable to further degeneration in ApoE4 carriers as they enter middle and old age

Test-retest reliability of structural brain networks from diffusion MRI

Structural brain networks constructed from diffusion MRI (dMRI) and tractography have been demonstrated in healthy volunteers and more recently in various disorders affecting brain connectivity. However, few studies have addressed the reproducibility of the resulting networks. We measured the test–retest properties of such networks by varying several factors affecting network construction using ten healthy volunteers who underwent a dMRI protocol at 1.5 T on two separate occasions. Each T1-weighted brain was parcellated into 84 regions-of-interest and network connections were identified using dMRI and two alternative tractography algorithms, two alternative seeding strategies, a white matter waypoint constraint and three alternative network weightings. In each case, four common graph-theoretic measures were obtained. Network properties were assessed both node-wise and per network in terms of the intraclass correlation coefficient (ICC) and by comparing within- and between-subject differences. Our findings suggest that test–retest performance was improved when: 1) seeding from white matter, rather than grey; and 2) using probabilistic tractography with a two-fibre model and sufficient streamlines, rather than deterministic tensor tractography. In terms of network weighting, a measure of streamline density produced better test–retest performance than tract-averaged diffusion anisotropy, although it remains unclear which is a more accurate representation of the underlying connectivity. For the best performing configuration, the global within-subject differences were between 3.2% and 11.9% with ICCs between 0.62 and 0.76. The mean nodal within-subject differences were between 5.2% and 24.2% with mean ICCs between 0.46 and 0.62. For 83.3% (70/84) of nodes, the within-subject differences were smaller than between-subject differences. Overall, these findings suggest that whilst current techniques produce networks capable of characterising the genuine between-subject differences in connectivity, future work must be undertaken to improve network reliability