Publications by Year: 2011

This paper proposes a method for the registration of white matter tract bundles traced from diffusion images and its extension to atlas generation, Our framework is based on a Gaussian process representation of tract density maps. Such a representation avoids the need for point-to-point correspondences, is robust to tract interruptions and reconnections and seamlessly handles the comparison and combination of white matter tract bundles. Moreover, being a parametric model, this approach has the potential to be defined in the Gaussian processes’ parameter space, without the need for resampling the fiber bundles during the registration process. We use the similarity measure of our Gaussian process framework, which is in fact an inner product, to drive a diffeomorphic registration algorithm between two sets of homologous bundles which is not biased by point-to-point correspondences or the parametrization of the tracts. We estimate a dense deformation of the underlying white matter using the bundles as anatomical landmarks and obtain a population atlas of those fiber bundles. Finally we test our results in several different bundles obtained from in-vivo data.

BACKGROUND: Structural abnormalities in the callosal fibers connecting the heteromodal association areas of the prefrontal and temporoparietal cortices bilaterally have been suggested to play a role in the etiology of schizophrenia. AIMS: To investigate for geometric abnormalities in these callosal fibers in schizophrenia patients by using a novel Diffusion-Tensor Imaging (DTI) metric of fiber geometry named Shape-Normalized Dispersion (SHD). METHODS: DTIs (3T, 51 gradient directions, 1.7mm isotropic voxels) were acquired from 26 schizophrenia patients and 23 matched healthy controls. The prefrontal and temporoparietal fibers of the corpus callosum were extracted by means of whole-brain tractography, and their mean SHD calculated. RESULTS: The schizophrenia patients exhibited subnormal levels of SHD in the prefrontal callosal fibers when controlling for between-group differences in Fractional Anisotropy. Reduced SHD could reflect either irregularly turbulent or inhomogeneously distributed fiber trajectories in the corpus callosum. CONCLUSIONS: The results suggest that the transcallosal misconnectivity thought to be associated with schizophrenia could reflect abnormalities in fiber geometry. These abnormalities in fiber geometry could potentially be underpinned by neurodevelopmental irregularities.

Voxel-based morphometry (VBM) is a hypothesis-free, whole-brain, voxel-by-voxel analytic method that attempts to compare imaging data between populations. Schizophrenia studies have utilized this method to localize differences in diffusion tensor imaging (DTI) derived fractional anisotropy (FA), a measure of white matter integrity, between patients and healthy controls. The number of publications has grown, although it is unclear how reliable and reproducible this method is, given the subtle white matter abnormalities expected in schizophrenia. Here we analyze and combine results from 23 studies published to date that use VBM to study schizophrenia in order to evaluate the reproducibility of this method in DTI analysis. Coordinates of each region reported in DTI VBM studies published thus far in schizophrenia were plotted onto a Montreal Neurological Institute atlas, and their anatomical locations were recorded. Results indicated that the reductions of FA in patients with schizophrenia were scattered across the brain. Moreover, even the most consistently reported regions were reported independently in less than 35% of the articles studied. Other instances of reduced FA were replicated at an even lower rate. Our findings demonstrate striking inconsistency, with none of the regions reported in much more than a third of the published articles. This poor replication rate suggests that the application of VBM to DTI data may not be the optimal way for finding the subtle microstructural abnormalities suggested in schizophrenia.

BACKGROUND: Abnormalities within language-related anatomical structures have been associated with clinical symptoms and with language and memory deficits in schizophrenia. Recent studies suggest disruptions in functional connectivity within the Inferior Frontal Gyrus (IFG) network in schizophrenia. However, due to technical challenges, anatomical connectivity abnormalities within this network and their involvement in clinical and cognitive deficits have not been studied. MATERIAL AND METHODS: Diffusion and anatomical scans were obtained from 23 chronic schizophrenia patients and 23 matched controls. The IFG was automatically segmented, and its white matter connections extracted and measured with newly-developed stochastic tractography tools. Correlations between anatomical structures and measures of semantic processing were also performed. RESULTS: White Matter connections between the IFG and posterior brain regions followed two distinct pathways: dorsal and ventral. Both demonstrated left lateralization, but ventral pathway abnormalities were only found in schizophrenia. IFG volumes also showed left lateralization and abnormalities in schizophrenia. Further, despite similar laterality and abnormality patterns, IFG volumes and white matter connectivity were not correlated with each other in either group. Interestingly, measures of semantic processing correlated with white matter connectivity in schizophrenia and with gray matter volumes in controls. Finally, hallucinations were best predicted by both gray matter and white matter measures together. CONCLUSIONS: Our results suggest abnormalities within the ventral IFG network in schizophrenia, with white matter abnormalities better predicting semantic deficits. The lack of a statistical relationship between coexisting gray and white matter deficits might suggest their different origin and the necessity for a multimodal approach in future schizophrenia studies.

The vervet monkey is an important nonhuman primate model that allows the study of isolated environmental factors in a controlled environment. Analysis of monkey MRI often suffers from lower quality images compared with human MRI because clinical equipment is typically used to image the smaller monkey brain and higher spatial resolution is required. This, together with the anatomical differences of the monkey brains, complicates the use of neuroimage analysis pipelines tuned for human MRI analysis. In this paper we developed an open source image analysis framework based on the tools available within the 3D Slicer software to support a biological study that investigates the effect of chronic ethanol exposure on brain morphometry in a longitudinally followed population of male vervets. We first developed a computerized atlas of vervet monkey brain MRI, which was used to encode the typical appearance of the individual brain structures in MRI and their spatial distribution. The atlas was then used as a spatial prior during automatic segmentation to process two longitudinal scans per subject. Our evaluation confirms the consistency and reliability of the automatic segmentation. The comparison of atlas construction strategies reveals that the use of a population-specific atlas leads to improved accuracy of the segmentation for subcortical brain structures. The contribution of this work is twofold. First, we describe an image processing workflow specifically tuned towards the analysis of vervet MRI that consists solely of the open source software tools. Second, we develop a digital atlas of vervet monkey brain MRIs to enable similar studies that rely on the vervet model.

The characterization of the distribution of noise in the magnitude MR image is a very important problem within image processing algorithms. The Rician noise assumed in single-coil acquisitions has been the keystone for signal-to-noise ratio estimation, image filtering, or diffusion tensor estimation for years. With the advent of parallel protocols such as sensitivity encoding or Generalized Autocalibrated Partially Parallel Acquisitions that allow accelerated acquisitions, this noise model no longer holds. Since Generalized Autocalibrated Partially Parallel Acquisitions reconstructions yield the combination of the squared signals recovered at each receiving coil, noncentral Chi statistics have been previously proposed to model the distribution of noise. However, we prove in this article that this is a weak model due to several artifacts in the acquisition scheme, mainly the correlation existing between the signals obtained at each coil. Alternatively, we propose to model such correlations with a reduction in the number of degrees of freedom of the signal, which translates in an equivalent nonaccelerated system with a minor number of independent receiving coils and, consequently, a lower signal-to-noise ratio. With this model, a noncentral Chi distribution can be assumed for all pixels in the image, whose effective number of coils and effective variance of noise can be explicitly computed in a closed form from the Generalized Autocalibrated Partially Parallel Acquisitions interpolation coefficients. Extensive experiments over both synthetic and in vivo data sets have been performed to show the goodness of fit of out model.

This work proposes a new method to detect abnormalities in fiber bundles of first-episode (FE) schizophrenia patients. Existing methods have either examined a particular region of interest or used voxel-based morphometry or used tracts generated using the single tensor model for locating statistically different fiber bundles. Further, a two-sample t test, which assumes a Gaussian distribution for each population, is the most widely used statistical hypothesis testing algorithm. In this study, we use the unscented Kalman filter based two-tensor tractography algorithm for tracing neural fiber bundles of the brain that connect 105 different cortical and subcortical regions. Next, fiber bundles with significant connectivity across the entire population were determined. Several diffusion measures derived from the two-tensor model were computed and used as features in the subsequent analysis. For each fiber bundle, an affine-invariant descriptor was computed, thus obviating the need for precise registration of patients to an atlas. A kernel-based statistical hypothesis testing algorithm, which makes no assumption regarding the distribution of the underlying population, was then used to determine the abnormal diffusion properties of all fiber bundles for 20 FE patients and 20 age-matched healthy controls. Of the 1254 fiber bundles with significant connectivity, 740 fiber bundles were found to be significantly different in at least one diffusion measure after correcting for multiple comparisons. Thus, the changes affecting first-episode patients seem to be global in nature (spread throughout the brain).

Arterial spin labeling is a noninvasive technique that can quantitatively measure cerebral blood flow. While traditionally arterial spin labeling employs 2D echo planar imaging or spiral acquisition trajectories, single-shot 3D gradient echo and spin echo (GRASE) is gaining popularity in arterial spin labeling due to inherent signal-to-noise ratio advantage and spatial coverage. However, a major limitation of 3D GRASE is through-plane blurring caused by T(2) decay. A novel technique combining 3D GRASE and a periodically rotated overlapping parallel lines with enhanced reconstruction trajectory (PROPELLER) is presented to minimize through-plane blurring without sacrificing perfusion sensitivity or increasing total scan time. Full brain perfusion images were acquired at a 3 × 3 × 5 mm(3) nominal voxel size with pulsed arterial spin labeling preparation sequence. Data from five healthy subjects was acquired on a GE 1.5T scanner in less than 4 minutes per subject. While showing good agreement in cerebral blood flow quantification with 3D gradient echo and spin echo, 3D GRASE PROPELLER demonstrated reduced through-plane blurring, improved anatomical details, high repeatability and robustness against motion, making it suitable for routine clinical use.

Parallel imaging methods are routinely used to accelerate the image acquisition process in cardiac cine imaging. The addition of a temporal acceleration method, whereby k-space is sampled differently for different time frames, has been shown in prior work to improve image quality as compared to parallel imaging by itself. However, such temporal acceleration strategies prove difficult to combine with retrospectively gated cine imaging. The only currently published method to feature such combination, by Hansen et al. [Magn Reson Med 55 (2006) 85-91] tends to be associated with prohibitively long reconstruction times. The goal of the present work was to develop a retrospectively gated cardiac cine method that features both parallel imaging and temporal acceleration, capable of achieving significant acceleration factors on commonly available hardware and associated with reconstruction times short enough for practical use in a clinical context. Seven cardiac patients and a healthy volunteer were recruited and imaged, with acceleration factors of 3.5 or 4.5, using an eight-channel product cardiac array on a 1.5-T system. The prescribed FOV value proved slightly too small in three patients, and one of the patients had a bigemini condition. Despite these additional challenges, good-quality results were obtained for all slices and all patients, with a reconstruction time of 0.98±0.07 s per frame, or about 20 s for a 20-frame slice, using a single processor on a single PC. As compared to using parallel imaging by itself, the addition of a temporal acceleration strategy provided much resistance to artifacts.

BACKGROUND: Activity-induced structural remodeling of dendritic spines and glial cells was recently proposed as an important factor in neuroplasticity and suggested to accompany the induction of long-term potentiation (LTP). Although T1 and diffusion MRI have been used to study structural changes resulting from long-term training, the cellular basis of the findings obtained and their relationship to neuroplasticity are poorly understood. METHODOLOGY/PRINCIPAL FINDING: Here we used diffusion tensor imaging (DTI) to examine the microstructural manifestations of neuroplasticity in rats that performed a spatial navigation task. We found that DTI can be used to define the selective localization of neuroplasticity induced by different tasks and that this process is age-dependent in cingulate cortex and corpus callosum and age-independent in the dentate gyrus. CONCLUSION/SIGNIFICANCE: We relate the observed DTI changes to the structural plasticity that occurs in astrocytes and discuss the potential of MRI for probing structural neuroplasticity and hence indirectly localizing LTP.