Methodology for computing white matter nerve fiber orientation in human histological slices.

BACKGROUND: The gold standard for mapping nerve fiber orientation in white matter of the human brain is histological analysis through biopsy. Such mappings are a crucial step in validating non-invasive techniques for assessing nerve fiber orientation in the human brain by using diffusion MRI. However, the manual extraction of nerve fiber directions of histological slices is tedious, time consuming, and prone to human error. NEW METHOD: The presented semi-automated algorithm first creates a binary-segmented mask of the nerve fibers in the histological image, and then extracts an estimate of average directionality of nerve fibers through a Fourier-domain analysis of the masked image. It also generates an uncertainty level for its estimate of average directionality. RESULTS AND COMPARISON WITH EXISTING METHODS: The average orientations of the semi-automatic method were first compared to a qualitative expert opinion based on visual inspection of nerve fibers. A weighted RMS difference between the expert estimate and the algorithmically determined angle (weighted by expert's confidence in his estimate) was 15.4°, dropping to 9.9° when only cases with an expert confidence level of greater than 50% were included. The algorithmically determined angles were then compared with angles extracted using a manual segmentation technique, yielding an RMS difference of 11.2°. CONCLUSION: The presented semi-automated method is in good agreement with both qualitative and quantitative manual expert-based approaches for estimating directionality of nerve fibers in white matter from images of stained histological slices of the human brain.

Mean diffusivity and fractional anisotropy as indicators of disease and genetic liability to schizophrenia.

The goals of this study were to first determine whether the fractional anisotropy (FA) and mean diffusivity (MD) of major white matter pathways associate with schizophrenia, and secondly to characterize the extent to which differences in these metrics might reflect a genetic predisposition to schizophrenia. Differences in FA and MD were identified using a comprehensive atlas-based tract mapping approach using diffusion tensor imaging and high-resolution structural data from 35 patients, 28 unaffected first-degree relatives of patients, 29 community controls, and 14 first-degree relatives of controls. Schizophrenia patients had significantly higher MD in the following tracts compared to controls: the right anterior thalamic radiations, the forceps minor, the bilateral inferior fronto-occipital fasciculus (IFO), the temporal component of the left superior longitudinal fasciculus (tSLF), and the bilateral uncinate. FA showed schizophrenia effects and a linear relationship to genetic liability (represented by schizophrenia patients, first-degree relatives, and controls) for the bilateral IFO, the left inferior longitudinal fasciculus (ILF), and the left tSLF. Diffusion tensor imaging studies have previously identified white matter abnormalities in all three of these tracts in schizophrenia; however, this study is the first to identify a significant genetic liability. Thus, FA of these three tracts may serve as biomarkers for studies seeking to identify how genes influence brain structure predisposing to schizophrenia. However, differences in FA and MD in frontal and temporal white matter pathways may be additionally driven by state variables that involve processes associated with the disease.

Patterns of fractional anisotropy changes in white matter of cerebellar peduncles distinguish spinocerebellar ataxia-1 from multiple system atrophy and other ataxia syndromes.

AIM: To determine prospectively if qualitative and quantitative diffusion tensor imaging (DTI) metrics of white matter integrity are better than conventional magnetic resonance imaging (MRI) metrics for discriminating cerebellar diseases. METHODS: Conventional MRI images from 31 consecutive patients with ataxia and 12 controls were interpreted by a neuroradiologist given only a clinical indication of ataxia. An expert ataxologist, blinded to radiological findings, determined the clinical diagnosis, as well as ataxia severity and asymmetry for each patient. For qualitative analysis, a comparison of the cerebellar white matter in ataxic vs. control patients was made by visual inspection of directionally encoded color (DEC) images. For quantitative analysis, segmentation of the cerebellar white matter in the inferior, middle, and superior cerebellar peduncles (ICP, MCP, and SCP) was attempted using three methods: a region of interest method, a deterministic DTI tractography (DDT) method, and a probabilistic DTI tractography (PDT) method. A statistical comparison of the average fractional anisotropy (FA) in these tracts was made between subject groups, and correlated to clinical diagnosis, severity, and asymmetry. RESULTS: Of the 31 consecutive patients with ataxia, the two largest subgroups had a clinical diagnosis of multiple system atrophy (cerebellar subtype; MSA-C), and spinocerebellar ataxia-1 (SCA1). Conventional MRI features, such as degree of pontocerebellar atrophy, correlated with ataxia severity, but were neither sensitive nor specific for the ataxia subtypes. PDT was the most accurate and least variable method of the three methods used for determining FA, especially in the ICP. Average FA in all ataxic patients was significantly decreased in the MCP, SCP and ICP and this decrease correlated to disease severity. Asymmetric ataxia correlated to proportionately larger contralateral MCP, ICP and SCP FA values. MCP, ICP, and SCP FA difference values formed distinct clusters that distinguished MSA-C from SCA-1, and other ataxia syndromes. CONCLUSIONS: Qualitative and quantitative reductions in DTI metrics of white matter integrity in the cerebellar peduncles correlated better to clinical features of patients with sporadic and hereditary ataxias than conventional structural MRI measures of pontocerebellar atrophy.

A diffusion tensor imaging tractography algorithm based on Navier-Stokes fluid mechanics.

We introduce a fluid mechanics based tractography method for estimating the most likely connection paths between points in diffusion tensor imaging (DTI) volumes. We customize the Navier-Stokes equations to include information from the diffusion tensor and simulate an artificial fluid flow through the DTI image volume. We then estimate the most likely connection paths between points in the DTI volume using a metric derived from the fluid velocity vector field. We validate our algorithm using digital DTI phantoms based on a helical shape. Our method segmented the structure of the phantom with less distortion than was produced using implementations of heat-based partial differential equation (PDE) and streamline based methods. In addition, our method was able to successfully segment divergent and crossing fiber geometries, closely following the ideal path through a digital helical phantom in the presence of multiple crossing tracts. To assess the performance of our algorithm on anatomical data, we applied our method to DTI volumes from normal human subjects. Our method produced paths that were consistent with both known anatomy and directionally encoded color images of the DTI dataset.

Fiber tractography reveals disruption of temporal lobe white matter tracts in schizophrenia.

Diffusion tensor imaging (DTI) studies have demonstrated abnormal anisotropic diffusion in schizophrenia. However, examining data with low spatial resolution and/or a low number of gradient directions and limitations associated with analysis approaches sensitive to registration confounds may have contributed to mixed findings concerning the regional specificity and direction of results. This study examined three major white matter tracts connecting lateral and medial temporal lobe regions with neocortical association regions widely implicated in systems-level functional and structural disturbances in schizophrenia. Using DTIstudio, a previously validated regions of interest tractography method was applied to 30 direction diffusion weighted imaging data collected from demographically similar schizophrenia (n=23) and healthy control subjects (n=22). The diffusion tensor was computed at each voxel after intra-subject registration of diffusion-weighted images. Three-dimensional tract reconstruction was performed using the Fiber Assignment by Continuous Tracking (FACT) algorithm. Tractography results showed reduced fractional anisotropy (FA) of the arcuate fasciculi (AF) and inferior longitudinal fasciculi (ILF) in patients compared to controls. FA changes within the right ILF were negatively correlated with measures of thinking disorder. Reduced volume of the left AF was also observed in patients. These results, which avoid registration issues associated with voxel-based analyses of DTI data, support that fiber pathways connecting lateral and medial temporal lobe regions with neocortical regions are compromised in schizophrenia. Disruptions of connectivity within these pathways may potentially contribute to the disturbances of memory, language, and social cognitive processing that characterize the disorder.

Mean diffusivity: a biomarker for CSF-related disease and genetic liability effects in schizophrenia.

Mean diffusivity (MD), the rotationally invariant magnitude of water diffusion that is greater in cerebrospinal fluid (CSF) and smaller in organized brain tissue, has been suggested to reflect schizophrenia-associated cortical atrophy. Regional changes, associations with CSF, and the effects of genetic predisposition towards schizophrenia, however, remain uncertain. Six-direction diffusion tensor imaging DTI and high-resolution structural images were obtained from 26 schizophrenia patients, 36 unaffected first-degree patient relatives, 20 control subjects and 32 control relatives (N=114). Registration procedures aligned diffusion tensor imaging (DTI) data across imaging modalities. MD was averaged within lobar regions and the cingulate and superior temporal gyri. CSF volume and MD were highly correlated. Significant bilateral temporal, and superior temporal MD increases were observed in schizophrenia compared with unrelated control probands. First-degree relatives of schizophrenia probands showed larger MD measures compared with controls within bilateral superior temporal regions with CSF volume correction. Superior temporal lobe brain tissue deficits and proximal CSF enlargements are widely documented in schizophrenia. Larger MD indices in patients and their relatives may thus reflect similar pathophysiological mechanisms. However, persistence of regional MD effects after controlling for CSF volume, suggests that MD is a sensitive biological marker of disease and genetic liability, characterizing at least partially distinct aspects of brain structural integrity.

Columnar specificity of microvascular oxygenation and volume responses: implications for functional brain mapping.

Cortical neurons with similar properties are grouped in columnar structures and supplied by matching vascular networks. The hemodynamic response to neuronal activation, however, is not well described on a fine spatial scale. We investigated the spatiotemporal characteristics of microvascular responses to neuronal activation in rat barrel cortex using optical intrinsic signal imaging and spectroscopy. Imaging was performed at 570 nm to provide functional maps of cerebral blood volume (CBV) changes and at 610 nm to estimate oxygenation changes. To emphasize parenchymal rather than large vessel contributions to the functional hemodynamic responses, we developed an ANOVA-based statistical analysis technique. Perfusion-based maps were compared with underlying neuroanatomy with cytochrome oxidase staining. Statistically determined CBV responses localized accurately to individually stimulated barrel columns and could resolve neighboring columns with a resolution better than 400 microm. Both CBV and early oxygenation responses extended beyond anatomical boundaries of single columns, but this vascular point spread did not preclude spatial specificity. These results indicate that microvascular flow control structures providing targeted flow increases to metabolically active neuronal columns also produce finely localized changes in CBV. This spatial specificity, along with the high contrast/noise ratio, makes the CBV response an attractive mapping signal. We also found that functional oxygenation changes can achieve submillimeter specificity not only during the transient deoxygenation ("initial dip") but also during the early part of the hyperoxygenation. We, therefore, suggest that to optimize hemodynamic spatial specificity, appropriate response timing (using < or =2-3 sec changes) is more important than etiology (oxygenation or volume).