Automatic clustering of white matter fibers in brain diffusion mri with an application to genetics

To understand factors that affect brain connectivity and integrity, it is beneficial to automatically cluster white matter (WM) fibers into anatomically recognizable tracts. Whole brain tractography, based on diffusion-weighted MRI, generates vast sets of fibers throughout the brain; clustering them into consistent and recognizable bundles can be difficult as there are wide individual variations in the trajectory and shape of WM pathways. Here we introduce a novel automated tract clustering algorithm based on label fusion - a concept from traditional intensity-based segmentation. Streamline tractography generates many incorrect fibers, so our top-down approach extracts tracts consistent with known anatomy, by mapping multiple hand-labeled atlases into a new dataset. We fuse clustering results from different atlases, using a mean distance fusion scheme. We reliably extracted the major tracts from 105-gradient high angular resolution diffusion images (HARDI) of 198 young normal twins. To compute population statistics, we use a pointwise correspondence method to match, compare, and average WM tracts across subjects. We illustrate our method in a genetic study of white matter tract heritability in twins.

Automatic population HARDI white matter tract clustering by label fusion of multiple tract atlases

Automatic labeling of white matter fibres in diffusion-weighted brain MRI is vital for comparing brain integrity and connectivity across populations, but is challenging. Whole brain tractography generates a vast set of fibres throughout the brain, but it is hard to cluster them into anatomically meaningful tracts, due to wide individual variations in the trajectory and shape of white matter pathways. We propose a novel automatic tract labeling algorithm that fuses information from tractography and multiple hand-labeled fibre tract atlases. As streamline tractography can generate a large number of false positive fibres, we developed a top-down approach to extract tracts consistent with known anatomy, based on a distance metric to multiple hand-labeled atlases. Clustering results from different atlases were fused, using a multi-stage fusion scheme. Our "label fusion" method reliably extracted the major tracts from 105-gradient HARDI scans of 100 young normal adults. © 2012 Springer-Verlag.

Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: Comparing meta and megaanalytical approaches for data pooling

Combining datasets across independent studies can boost statistical power by increasing the numbers of observations and can achieve more accurate estimates of effect sizes. This is especially important for genetic studies where a large number of observations are required to obtain sufficient power to detect and replicate genetic effects. There is a need to develop and evaluate methods for joint-analytical analyses of rich datasets collected in imaging genetics studies. The ENIGMA-DTI consortium is developing and evaluating approaches for obtaining pooled estimates of heritability through meta-and mega-genetic analytical approaches, to estimate the general additive genetic contributions to the intersubject variance in fractional anisotropy (FA) measured from diffusion tensor imaging (DTI). We used the ENIGMA-DTI data harmonization protocol for uniform processing of DTI data from multiple sites. We evaluated this protocol in five family-based cohorts providing data from a total of 2248 children and adults (ages: 9-85) collected with various imaging protocols. We used the imaging genetics analysis tool, SOLAR-Eclipse, to combine twin and family data from Dutch, Australian and Mexican-American cohorts into one large "mega-family". We showed that heritability estimates may vary from one cohort to another. We used two meta-analytical (the sample-size and standard-error weighted) approaches and a mega-genetic analysis to calculate heritability estimates across-population. We performed leave-one-out analysis of the joint estimates of heritability, removing a different cohort each time to understand the estimate variability. Overall, meta- and mega-genetic analyses of heritability produced robust estimates of heritability.

Predicting white matter integrity from multiple common genetic variants

Several common genetic variants have recently been discovered that appear to influence white matter microstructure, as measured by diffusion tensor imaging (DTI). Each genetic variant explains only a small proportion of the variance in brain microstructure, so we set out to explore their combined effect on the white matter integrity of the corpus callosum. We measured six common candidate single-nucleotide polymorphisms (SNPs) in the COMT, NTRK1, BDNF, ErbB4, CLU, and HFE genes, and investigated their individual and aggregate effects on white matter structure in 395 healthy adult twins and siblings (age: 20-30 years). All subjects were scanned with 4-tesla 94-direction high angular resolution diffusion imaging. When combined using mixed-effects linear regression, a joint model based on five of the candidate SNPs (COMT, NTRK1, ErbB4, CLU, and HFE) explained ∼ 6% of the variance in the average fractional anisotropy (FA) of the corpus callosum. This predictive model had detectable effects on FA at 82% of the corpus callosum voxels, including the genu, body, and splenium. Predicting the brain's fiber microstructure from genotypes may ultimately help in early risk assessment, and eventually, in personalized treatment for neuropsychiatric disorders in which brain integrity and connectivity are affected.

White matter integrity measured by fractional anisotropy correlates poorly with actual individual fiber anisotropy

Fractional anisotropy (FA), a very widely used measure of fiber integrity based on diffusion tensor imaging (DTI), is a problematic concept as it is influenced by several quantities including the number of dominant fiber directions within each voxel, each fiber's anisotropy, and partial volume effects from neighboring gray matter. High-angular resolution diffusion imaging (HARDI) can resolve more complex diffusion geometries than standard DTI, including fibers crossing or mixing. The tensor distribution function (TDF) can be used to reconstruct multiple underlying fibers per voxel, representing the diffusion profile as a probabilistic mixture of tensors. Here we found that DTIderived mean diffusivity (MD) correlates well with actual individual fiber MD, but DTI-derived FA correlates poorly with actual individual fiber anisotropy, and may be suboptimal when used to detect disease processes that affect myelination. Analysis of the TDFs revealed that almost 40% of voxels in the white matter had more than one dominant fiber present. To more accurately assess fiber integrity in these cases, we here propose the differential diffusivity (DD), which measures the average anisotropy based on all dominant directions in each voxel.

Scalar connectivity measures from fast-marching tractography reveal heritability of white matter architecture

Recent advances in diffusion-weighted MRI (DWI) have enabled studies of complex white matter tissue architecture in vivo. To date, the underlying influence of genetic and environmental factors in determining central nervous system connectivity has not been widely studied. In this work, we introduce new scalar connectivity measures based on a computationally-efficient fast-marching algorithm for quantitative tractography. We then calculate connectivity maps for a DTI dataset from 92 healthy adult twins and decompose the genetic and environmental contributions to the variance in these metrics using structural equation models. By combining these techniques, we generate the first maps to directly examine genetic and environmental contributions to brain connectivity in humans. Our approach is capable of extracting statistically significant measures of genetic and environmental contributions to neural connectivity.

Automatic clustering and population analysis of white matter tracts using maximum density paths

We introduce a framework for population analysis of white matter tracts based on diffusion-weighted images of the brain. The framework enables extraction of fibers from high angular resolution diffusion images (HARDI); clustering of the fibers based partly on prior knowledge from an atlas; representation of the fiber bundles compactly using a path following points of highest density (maximum density path; MDP); and registration of these paths together using geodesic curve matching to find local correspondences across a population. We demonstrate our method on 4-Tesla HARDI scans from 565 young adults to compute localized statistics across 50 white matter tracts based on fractional anisotropy (FA). Experimental results show increased sensitivity in the determination of genetic influences on principal fiber tracts compared to the tract-based spatial statistics (TBSS) method. Our results show that the MDP representation reveals important parts of the white matter structure and considerably reduces the dimensionality over comparable fiber matching approaches.

Dynamics of gray matter loss in Alzheimer's disease

We detected and mapped a dynamically spreading wave of gray matter loss in the brains of patients with Alzheimer's disease (AD). The loss pattern was visualized in four dimensions as it spread over time from temporal and limbic cortices into frontal and occipital brain regions, sparing sensorimotor cortices. The shifting deficits were asymmetric (left hemisphere > right hemisphere) and correlated with progressively declining cognitive status (p < 0.0006). Novel brain mapping methods allowed us to visualize dynamic patterns of atrophy in 52 high-resolution magnetic resonance image scans of 12 patients with AD (age 68.4 ± 1.9 years) and 14 elderly matched controls (age 71.4 ± 0.9 years) scanned longitudinally (two scans; interscan interval 2.1 ± 0.4 years). A cortical pattern matching technique encoded changes in brain shape and tissue distribution across subjects and time. Cortical atrophy occurred in a well defined sequence as the disease progressed, mirroring the sequence of neurofibrillary tangle accumulation observed in cross sections at autopsy. Advancing deficits were visualized as dynamic maps that change over time. Frontal regions, spared early in the disease, showed pervasive deficits later (< 15% loss). The maps distinguished different phases of AD and differentiated AD from normal aging. Local gray matter loss rates (5.3 ± 2.3% per year in AD v 0.9 ± 0.9% per year in controls) were faster in the left hemisphere (p < 0.029) than the right. Transient barriers to disease progression appeared at limbic/frontal boundaries. This degenerative sequence, observed in vivo as it developed, provides the first quantitative, dynamic visualization of cortical atrophic rates in normal elderly populations and in those with dementia.

Changes in white matter connectivity following therapy for anomia post stroke

Background. The majority of studies investigating the neural mechanisms underlying treatment-induced recovery in aphasia have focused on the cortical regions associated with language processing. However, the integrity of the white matter connecting these regions may also be crucial to understanding treatment mechanisms. Objective. This study investigated the integrity of the arcuate fasciculus (AF) and uncinate fasciculus (UF) before and after treatment for anomia in people with aphasia. Method. Eight people with aphasia received 12 treatment sessions to improve naming; alternating between phonologically-based and semantic-based tasks, with high angular resolution diffusion imaging conducted pre and post treatment. The mean generalized fractional anisotropy (GFA), a measure of fiber integrity, and number of fibers in the AF and UF were compared pre and post treatment, as well as with a group of 14 healthy older controls. Results. Pre treatment, participants with aphasia had significantly fewer fibers and lower mean GFA in the left AF compared with controls. Post treatment, mean GFA increased in the left AF to be statistically equivalent to controls. Additionally, mean GFA in the left AF pre and post treatment positively correlated with maintenance of the phonologically based treatment. No differences were found in the right AF, or the UF in either hemisphere, between participants with aphasia and controls, and no changes were observed in these tracts following treatment. Conclusions. Anomia treatments may improve the integrity of the white matter connecting cortical language regions. These preliminary results add to the understanding of the mechanisms underlying treatment outcomes in people with aphasia post stroke.

Serum cholesterol and variant in cholesterol-related gene CETP predict white matter microstructure

Several common genetic variants influence cholesterol levels, which play a key role in overall health. Myelin synthesis and maintenance are highly sensitive to cholesterol concentrations, and abnormal cholesterol levels increase the risk for various brain diseases, including Alzheimer's disease. We report significant associations between higher serum cholesterol (CHOL) and high-density lipoprotein levels and higher fractional anisotropy in 403 young adults (23.8 ± 2.4years) scanned with diffusion imaging and anatomic magnetic resonance imaging at 4Tesla. By fitting a multi-locus genetic model within white matter areas associated with CHOL, we found that a set of 18 cholesterol-related, single-nucleotide polymorphisms implicated in Alzheimer's disease risk predicted fractional anisotropy. We focused on the single-nucleotide polymorphism with the largest individual effects, CETP (rs5882), and found that increased G-allele dosage was associated with higher fractional anisotropy and lower radial and mean diffusivities in voxel-wise analyses of the whole brain. A follow-up analysis detected white matter associations with rs5882 in the opposite direction in 78 older individuals (74.3 ± 7.3years). Cholesterol levels may influence white matter integrity, and cholesterol-related genes may exert age-dependent effects on the brain.

Identifying candidate gene effects by restricting search space in a multivariate genetic analysis of white matter microstructure

Several genetic variants are thought to influence white matter (WM) integrity, measured with diffusion tensor imaging (DTI). Voxel based methods can test genetic associations, but heavy multiple comparisons corrections are required to adjust for searching the whole brain and for all genetic variants analyzed. Thus, genetic associations are hard to detect even in large studies. Using a recently developed multi-SNP analysis, we examined the joint predictive power of a group of 18 cholesterol-related single nucleotide polymorphisms (SNPs) on WM integrity, measured by fractional anisotropy. To boost power, we limited the analysis to brain voxels that showed significant associations with total serum cholesterol levels. From this space, we identified two genes with effects that replicated in individual voxel-wise analyses of the whole brain. Multivariate analyses of genetic variants on a reduced anatomical search space may help to identify SNPs with strongest effects on the brain from a broad panel of genes.

Stabilization of stormwater biofilters: Impacts of wetting and drying phases and the addition of organic matter to filter media

Ripening period refers to a phase of stabilization in sand filters in water treatment systems that follows a new installation or cleaning of the filter. Intermittent wetting and drying, a unique property of stormwater biofilters, would similarly be subjected to a phase of stabilization. Suspended solids, is an important parameter that is often used to monitor the stabilization of sand filters in water treatment systems. Stormwater biofilters however, contain organic material that is added to the filter layer to enhance nitrate removal, the dynamics of which is seldom analysed in stabilization of stormwater biofilters. Therefore, in this study of stormwater biofiltration in addition to suspended solids (Turbidity), organic matter (TOC, DOC, TN and TKN) was also monitored as a parameter for stabilization of the stormwater biofilter. One Perspex bioretention column (94 mm internal diameter) was fabricated with filter layer that contained 8% organic material and fed with tapwater with different antecedent dry days (0 – 40 day) at 100 mL/min. Samples were collected from the outflow at different time intervals between 2 – 150 minutes and were tested for Total Organic Carbon, Dissolved Organic Carbon, Total Nitrogen, Total Kjeldhal Nitrogen and Turbidity. The column was observed to experience two phases of stabilization, one at the beginning of each event that lasted for 30 minutes while the other phase was observed across subsequent events that related to the age of filter.

Effects of biodiesel chemical composition on the chemical and physical properties of the primary and secondary diesel particulate matter

This thesis improves our insight towards the effects of using biodiesels on the particulate matter emission of diesel engines and contributes to our understanding of their potential adverse health effects. The novelty of this project is the use of biodiesel fuel with controlled chemical composition that enables us to relate changes of physiochemical properties of particles to specific properties of the biodiesel. For the first time, the possibility of a correlation of the volatility and the Reactive Oxygen Species concentration of the particles is investigated versus the saturation, oxygen content and carbon chain length of the fuel.

Can’t see the trees for the forest? Why IS-ServQual items matter

Despite longstanding concern with the dimensionality of the service quality construct as measured by ServQual and IS-ServQual instruments, variations on the IS-ServQual instrument have been enduringly prominent in both academic research and practice in the field of IS. We explain the continuing popularity of the instrument based on the salience of the item set for predicting overall customer satisfaction, suggesting that the preoccupation with the dimensions has been a distraction. The implicit mutual exclusivity of the items suggests a more appropriate conceptualization of IS-ServQual as a formative index. This conceptualization resolves the paradox in IS-ServQual research, that of how an instrument with such well-known and well-documented weaknesses continue to be very influential and widely used by academics and practitioners. A formative conceptualization acknowledges and addresses the criticisms of IS-ServQual, while simultaneously explaining its enduring salience by focusing on the items rather than the “dimensions.” By employing an opportunistic sample and adopting the most recent IS-ServQual instrument published in a leading IS journal (virtually, any valid IS- ServQual sample in combination with a previously tested instrument variant would suffice for study purposes), we demonstrate that when re-specified as both first-order and second-order formatives, IS-ServQual has good model quality metrics and high predictive power on customer satisfaction. We conclude that this formative specification has higher practical use and is more defensible theoretically.

A preconditioned numerical solver for stiff nonlinear reaction-diffusion equations with fractional Laplacians that avoids dense matrices

The numerical solution of fractional partial differential equations poses significant computational challenges in regard to efficiency as a result of the spatial nonlocality of the fractional differential operators. The dense coefficient matrices that arise from spatial discretisation of these operators mean that even one-dimensional problems can be difficult to solve using standard methods on grids comprising thousands of nodes or more. In this work we address this issue of efficiency for one-dimensional, nonlinear space-fractional reaction–diffusion equations with fractional Laplacian operators. We apply variable-order, variable-stepsize backward differentiation formulas in a Jacobian-free Newton–Krylov framework to advance the solution in time. A key advantage of this approach is the elimination of any requirement to form the dense matrix representation of the fractional Laplacian operator. We show how a banded approximation to this matrix, which can be formed and factorised efficiently, can be used as part of an effective preconditioner that accelerates convergence of the Krylov subspace iterative solver. Our approach also captures the full contribution from the nonlinear reaction term in the preconditioner, which is crucial for problems that exhibit stiff reactions. Numerical examples are presented to illustrate the overall effectiveness of the solver.