Publications

Diffusion MRI has been successful in identifying the existence of white matter abnormalities in schizophrenia in vivo. However, the role of these abnormalities in the etiology of schizophrenia is not well understood. Accumulating evidence from imaging, histological, genetic, and immunochemical studies support the involvement of axonal degeneration and neuroinflammation—ubiquitous components of neurodegenerative disorders—as the underlying pathologies of these abnormalities. Nevertheless, the current imaging modalities cannot distinguish neuroinflammation from axonal degeneration, and therefore provide little specificity with respect to the pathophysiology progression and whether it is related to a neurodegenerative process. Free-water imaging is a new methodology that is sensitive to water molecules diffusing in the extracellular space. Excessive extracellular volume is a surrogate biomarker for neuroinflammation and can be separated out to reveal abnormalities such as axonal degeneration that affect diffusion characteristics in the tissue. We applied free-water imaging on diffusion MRI data acquired from schizophrenia-diagnosed human subjects with a first psychotic episode. We found a significant increase in the extracellular volume in both white and gray matter. In contrast, significant signs of axonal degeneration were limited to focal areas in the frontal lobe white matter. Our findings demonstrate that neuroinflammation is more prominent than axonal degeneration in the early stage of schizophrenia, revealing a pattern shared by many neurodegenerative disorders, in which prolonged inflammation leads to axonal degeneration. These findings promote anti-inflammatory treatment for early diagnosed schizophrenia patients.

The properties of carbon nanotube (CNT)/polymer composites are strongly dependent on the dispersion and orientation of CNTs in the host matrix. Quantification of the dispersion and orientation of CNTs by means of microstructure observation and image analysis has been demonstrated as a useful way to understand the structure-property relationship of CNT/polymer composites. However, due to the various morphologies and large amount of CNTs in one image, automatic and accurate identification of CNTs has become the bottleneck for dispersion/orientation analysis. To solve this problem, shape identification is performed for each pixel in the filler identification step, so that individual CNTs can be extracted from images automatically. The improved filler identification enables more accurate analysis of CNT dispersion and orientation. The dispersion index and orientation index obtained for both synthetic and real images from model compounds correspond well with the observations. Moreover, these indices help to explain the electrical properties of CNT/silicone composite, which is used as a model compound. This method can also be extended to other polymer composites with high-aspect-ratio fillers.

The registration of a pair of point sets as well as the estimation of their pointwise correspondences is a challenging and important task in computer vision. In this paper, we present a method to estimate the diffeomorphic deformation, together with the pointwise correspondences, between a pair of point sets. Many of the registration problems are iteratively solved by estimating the correspondence, locally optimizing certain cost functionals over the rigid or similarity or affine transformation group, then estimating the correspondence again, and so on. This type of approach, however, is well-known to be susceptible to suboptimal local solutions. In this paper, we first adopt the perspective of treating the registration as a posterior estimation optimization problem and solve it accordingly via a particle-filtering framework. Second, within such a framework, the diffeomorphic registration is performed to correct the nonlinear deformation of the points. In doing so, we provide a solution less susceptible to local minima. We provide the experimental results, which include challenging medical data sets where the two point sets differ by 180 (°) rotation as well as local deformations, to highlight the algorithm’s capability of robustly finding the more globally optimal solution for the registration task.

BACKGROUND: Knowledge of the individual course of the optic radiations (ORs) is important to avoid postoperative visual deficits. Cadaveric studies of the visual pathways are limited because it has not been possible to separate the OR from neighboring tracts accurately and results may not apply to individual patients. Diffusion tensor imaging studies may be able to demonstrate the relationships between the OR and neighboring fibers in vivo in individual subjects. OBJECTIVE: To use diffusion tensor imaging tractography to study the OR and the Meyer loop (ML) anatomy in vivo. METHODS: Ten healthy subjects underwent magnetic resonance imaging with diffusion imaging at 3 T. With the use of a fiducial-based diffusion tensor imaging tractography tool (Slicer 3.3), seeds were placed near the lateral geniculate nucleus to reconstruct individual visual pathways and neighboring tracts. Projections of the ORs onto 3-dimensional brain models were shown individually to quantify relationships to key landmarks. RESULTS: Two patterns of visual pathways were found. The OR ran more commonly deep in the whole superior and middle temporal gyri and superior temporal sulcus. The OR was closely surrounded in all cases by an inferior longitudinal fascicle and a parieto/occipito/temporo-pontine fascicle. The mean left and right distances between the tip of the OR and temporal pole were 39.8 ± 3.8 and 40.6 ± 5.7 mm, respectively. CONCLUSION: Diffusion tensor imaging tractography provides a practical complementary method to study the OR and the Meyer loop anatomy in vivo with reference to individual 3-dimensional brain anatomy.

The desire to understand complex mental processes using functional MRI drives development of imaging techniques that scan the whole human brain at a high spatial and temporal resolution. In this work, an accelerated multishot three-dimensional echo-planar imaging sequence is proposed to increase the temporal resolution of these studies. A combination of two modern acceleration techniques, UNFOLD and GRAPPA is used in the secondary phase encoding direction to reduce the scan time effectively. The sequence (repetition time of 1.02 s) was compared with standard two-dimensional echo-planar imaging (3 s) and multishot three-dimensional echo-planar imaging (3 s) sequences with both block design and event-related functional MRI paradigms. With the same experimental setup and imaging time, the temporal resolution improvement with our sequence yields similar activation regions in the block design functional MRI paradigm with slightly increased t-scores. Moreover, additional information on the timing of rapid dynamic changes was extracted from accelerated images for the case of the event related complex mental paradigm.

Extracting anatomical and functional significant structures renders one of the important tasks for both the theoretical study of the medical image analysis, and the clinical and practical community. In the past, much work has been dedicated only to the algorithmic development. Nevertheless, for clinical end users, a well designed algorithm with an interactive software is necessary for an algorithm to be utilized in their daily work. Furthermore, the software would better be open sourced in order to be used and validated by not only the authors but also the entire community. Therefore, the contribution of the present work is twofolds: first, we propose a new robust statistics based conformal metric and the conformal area driven multiple active contour framework, to simultaneously extract multiple targets from MR and CT medical imagery in 3D. Second, an open source graphically interactive 3D segmentation tool based on the aforementioned contour evolution is implemented and is publicly available for end users on multiple platforms. In using this software for the segmentation task, the process is initiated by the user drawn strokes (seeds) in the target region in the image. Then, the local robust statistics are used to describe the object features, and such features are learned adaptively from the seeds under a non-parametric estimation scheme. Subsequently, several active contours evolve simultaneously with their interactions being motivated by the principles of action and reaction-this not only guarantees mutual exclusiveness among the contours, but also no longer relies upon the assumption that the multiple objects fill the entire image domain, which was tacitly or explicitly assumed in many previous works. In doing so, the contours interact and converge to equilibrium at the desired positions of the desired multiple objects. Furthermore, with the aim of not only validating the algorithm and the software, but also demonstrating how the tool is to be used, we provide the reader reproducible experiments that demonstrate the capability of the proposed segmentation tool on several public available data sets.

Thalamo-cortical feedback loops play a key role in the processing and coordination of processing and integration of perceptual inputs and outputs, and disruption in this connection has long been hypothesized to contribute significantly to neuropsychological disturbances in schizophrenia. To test this hypothesis, we applied diffusion tensor tractography to 18 patients suffering schizophrenia and 20 control subjects. Fractional anisotropy (FA) was evaluated in the bilateral anterior and posterior limbs of the internal capsule, and correlated with clinical and neurocognitive measures. Patients diagnosed with schizophrenia showed significantly reduced FA bilaterally in the anterior but not the posterior limb of the internal capsule, compared with healthy control subjects. Lower FA correlated with lower scores on tests of declarative episodic memory in the patient group only. These findings suggest that disruptions, bilaterally, in thalamo-cortical connections in schizophrenia may contribute to disease-related impairment in the coordination of mnemonic processes of encoding and retrieval that are vital for efficient learning of new information.

PURPOSE: Diffusion magnetic resonance imaging (MRI) in combination with functional MRI promises a whole new vista for scientists to investigate noninvasively the structural and functional connectivity of the human brain-the human connectome, which had heretofore been out of reach. As with other imaging modalities, diffusion MRI data are inherently noisy and its acquisition time-consuming. Further, a faithful representation of the human connectome that can serve as a predictive model requires a robust and accurate data-analytic pipeline. The focus of this paper is on one of the key segments of this pipeline-in particular, the development of a sparse and optimal acquisition (SOA) design for diffusion MRI multiple-shell acquisition and beyond. METHODS: The authors propose a novel optimality criterion for sparse multiple-shell acquisition and quasimultiple-shell designs in diffusion MRI and a novel and effective semistochastic and moderately greedy combinatorial search strategy with simulated annealing to locate the optimum design or configuration. The goal of the optimality criteria is threefold: first, to maximize uniformity of the diffusion measurements in each shell, which is equivalent to maximal incoherence in angular measurements; second, to maximize coverage of the diffusion measurements around each radial line to achieve maximal incoherence in radial measurements for multiple-shell acquisition; and finally, to ensure maximum uniformity of diffusion measurement directions in the limiting case when all the shells are coincidental as in the case of a single-shell acquisition. The approach taken in evaluating the stability of various acquisition designs is based on the condition number and the A-optimal measure of the design matrix. RESULTS: Even though the number of distinct configurations for a given set of diffusion gradient directions is very large in general-e.g., in the order of 10(232) for a set of 144 diffusion gradient directions, the proposed search strategy was found to be effective in finding the optimum configuration. It was found that the square design is the most robust (i.e., with stable condition numbers and A-optimal measures under varying experimental conditions) among many other possible designs of the same sample size. Under the same performance evaluation, the square design was found to be more robust than the widely used sampling schemes similar to that of 3D radial MRI and of diffusion spectrum imaging (DSI). CONCLUSIONS: A novel optimality criterion for sparse multiple-shell acquisition and quasimultiple-shell designs in diffusion MRI and an effective search strategy for finding the best configuration have been developed. The results are very promising, interesting, and practical for diffusion MRI acquisitions.

Features of the diffusion-time dependence of the diffusion-weighted magnetic resonance imaging (MRI) signal provide a new contrast that could be altered by numerous biological processes and pathologies in tissue at microscopic length scales. An anomalous diffusion model, based on the theory of Brownian motion in fractal and disordered media, is used to characterize the temporal scaling (TS) characteristics of diffusion-related quantities, such as moments of the displacement and zero-displacement probabilities, in excised rat hippocampus specimens. To reduce the effect of noise in magnitude-valued MRI data, a novel numerical procedure was employed to yield accurate estimation of these quantities even when the signal falls below the noise floor. The power-law dependencies characterize the TS behavior in all regions of the rat hippocampus, providing unique information about its microscopic architecture. The relationship between the TS characteristics and diffusion anisotropy is investigated by examining the anisotropy of TS, and conversely, the TS of anisotropy. The findings suggest the robustness of the technique as well as the reproducibility of estimates. TS characteristics of the diffusion-weighted signals could be used as a new and useful marker of tissue microstructure.

The overall goal of this research is the design of statistical atlas models that can be created from normal subjects, but may generalize to be applicable to abnormal brains. We present a new style of joint modeling of fMRI, DTI, and structural MRI. Motivated by the fact that a white matter tract and related cortical areas are likely to displace together in the presence of a mass lesion (brain tumor), in this work we propose a rotation and translation invariant model that represents the spatial relationship between fiber tracts and anatomic and functional landmarks. This landmark distance model provides a new basis for representation of fiber tracts and can be used for detection and prediction of fiber tracts based on landmarks. Our results indicate that the measured model is consistent across normal subjects, and thus suitable for atlas building. Our experiments demonstrate that the model is robust to displacement and missing data, and can be successfully applied to a small group of patients with mass lesions.