Publications by Year: 2010

Kahn T, Jolesz FA, Lewin JS. Proceedings of the 8th Interventional MRI Symposium. 8th Interventional MRI Symposium. 2010;8 :1-340. 2010 iMRI Symposium Proceedings
Kapur T, Tempany CM, Jolesz FA. Proceedings of the 3rd Image Guided Therapy Workshop. Image Guided Therapy Workshop. 2010;3 :1-87. 2010 IGT Workshop Proceedings
Whitford TJ, Kubicki M, Schneiderman JS, O'Donnell LJ, King R, Alvarado JL, Khan U, Markant D, Nestor PG, Niznikiewicz M, et al. Corpus Callosum Abnormalities and their Association with Psychotic Symptoms in Patients with Schizophrenia. Biol Psychiatry. 2010;68 (1) :70-7.Abstract

BACKGROUND: While the neuroanatomical underpinnings of the functional brain disconnectivity observed in patients with schizophrenia (SZ) remain elusive, white matter fiber bundles of the brain are a likely candidate, given that they represent the infrastructure for long-distance neural communication. METHODS: This study investigated for diffusion abnormalities in 19 patients with chronic SZ, relative to 19 matched control subjects, across tractography-defined segments of the corpus callosum. Diffusion-weighted images were acquired with 51 noncollinear gradients on a 3T scanner (1.7 mm isotropic voxels). The corpus callosum was extracted by means of whole-brain tractography and automated fiber clustering and was parcelled into six segments on the basis of fiber trajectories. The diffusion indexes of fractional anisotropy (FA) and mode were calculated for each segment. RESULTS: Relative to the healthy control subjects, the SZ patients exhibited mode increases in the parietal fibers, suggesting a relative absence of crossing fibers. Schizophrenia patients also exhibited FA reductions in the frontal fibers, which were underpinned by increases in radial diffusivity, consistent with myelin abnormalities. Significant correlations were observed between patients' degree of reality distortion and their FA and radial diffusivity, such that the most severely psychotic patients were the least abnormal in terms of their frontal fiber diffusivity. CONCLUSIONS: The SZ patients exhibited a variety of diffusion abnormalities in the corpus callosum, which were related to the severity of their psychotic symptoms. To the extent that diffusion abnormalities influence axonal transmission velocities, these results provide support for those theories that emphasize neural timing abnormalities in the etiology of schizophrenia.

Voineskos AN, Lobaugh NJ, Bouix S, Rajji TK, Miranda D, Kennedy JL, Mulsant BH, Pollock BG, Shenton ME. Diffusion Tensor Tractography Findings in Schizophrenia across the Adult Lifespan. Brain. 2010;133 (Pt 5) :1494-504.Abstract

In healthy adult individuals, late life is a dynamic time of change with respect to the microstructural integrity of white matter tracts. Yet, elderly individuals are generally excluded from diffusion tensor imaging studies in schizophrenia. Therefore, we examined microstructural integrity of frontotemporal and interhemispheric white matter tracts in schizophrenia across the adult lifespan. Diffusion tensor imaging data from 25 younger schizophrenic patients (< or = 55 years), 25 younger controls, 25 older schizophrenic patients (> or = 56 years) and 25 older controls were analysed. Patients with schizophrenia in each group were individually matched to controls. Whole-brain tractography and clustering segmentation were employed to isolate white matter tracts. Groups were compared using repeated measures analysis of variance with 12 within-group measures of fractional anisotropy: (left and right) uncinate fasciculus, arcuate fasciculus, inferior longitudinal fasciculus, inferior occipito-frontal fasciculus, cingulum bundle, and genu and splenium of the corpus callosum. For each white matter tract, fractional anisotropy was then regressed against age in patients and controls, and correlation coefficients compared. The main effect of group (F(3,92) = 12.2, P < 0.001), and group by tract interactions (F(26,832) = 1.68, P = 0.018) were evident for fractional anisotropy values. Younger patients had significantly lower fractional anisotropy than younger controls (Bonferroni-corrected alpha = 0.0042) in the left uncinate fasciculus (t(48) = 3.7, P = 0.001) and right cingulum bundle (t(48) = 3.6, P = 0.001), with considerable effect size, but the older groups did not differ. Schizophrenic patients did not demonstrate accelerated age-related decline compared with healthy controls in any white matter tract. To our knowledge, this is the first study to examine the microstructural integrity of frontotemporal white matter tracts across the adult lifespan in schizophrenia. The left uncinate fasciculus and right cingulum bundle are disrupted in younger chronic patients with schizophrenia compared with matched controls, suggesting that these white matter tracts are related to frontotemporal disconnectivity. The absence of accelerated age-related decline, or differences between older community-dwelling patients and controls, suggests that these patients may possess resilience to white matter disruption.

Tokuda J, Fischer GS, Dimaio SP, Gobbi DG, Csoma C, Mewes PW, Fichtinger G, Tempany CM, Hata N. Integrated navigation and control software system for MRI-guided robotic prostate interventions. Comput Med Imaging Graph. 2010;34 (1) :3-8.Abstract
A software system to provide intuitive navigation for MRI-guided robotic transperineal prostate therapy is presented. In the system, the robot control unit, the MRI scanner, and the open-source navigation software are connected together via Ethernet to exchange commands, coordinates, and images using an open network communication protocol, OpenIGTLink. The system has six states called "workphases" that provide the necessary synchronization of all components during each stage of the clinical workflow, and the user interface guides the operator linearly through these workphases. On top of this framework, the software provides the following features for needle guidance: interactive target planning; 3D image visualization with current needle position; treatment monitoring through real-time MR images of needle trajectories in the prostate. These features are supported by calibration of robot and image coordinates by fiducial-based registration. Performance tests show that the registration error of the system was 2.6mm within the prostate volume. Registered real-time 2D images were displayed 1.97 s after the image location is specified.
Yeo TBT, Sabuncu MR, Vercauteren T, Ayache N, Fischl B, Golland P. Spherical demons: fast diffeomorphic landmark-free surface registration. IEEE Trans Med Imaging. 2010;29 (3) :650-68.Abstract
We present the Spherical Demons algorithm for registering two spherical images. By exploiting spherical vector spline interpolation theory, we show that a large class of regularizors for the modified Demons objective function can be efficiently approximated on the sphere using iterative smoothing. Based on one parameter subgroups of diffeomorphisms, the resulting registration is diffeomorphic and fast. The Spherical Demons algorithm can also be modified to register a given spherical image to a probabilistic atlas. We demonstrate two variants of the algorithm corresponding to warping the atlas or warping the subject. Registration of a cortical surface mesh to an atlas mesh, both with more than 160 k nodes requires less than 5 min when warping the atlas and less than 3 min when warping the subject on a Xeon 3.2 GHz single processor machine. This is comparable to the fastest nondiffeomorphic landmark-free surface registration algorithms. Furthermore, the accuracy of our method compares favorably to the popular FreeSurfer registration algorithm. We validate the technique in two different applications that use registration to transfer segmentation labels onto a new image 1) parcellation of in vivo cortical surfaces and 2) Brodmann area localization in ex vivo cortical surfaces.
Vaswani N, Rathi Y, Yezzi A, Tannenbaum A. Deform PF-MT: particle filter with mode tracker for tracking nonaffine contour deformations. IEEE Trans Image Process. 2010;19 (4) :841-57.Abstract
We propose algorithms for tracking the boundary contour of a deforming object from an image sequence, when the nonaffine (local) deformation over consecutive frames is large and there is overlapping clutter, occlusions, low contrast, or outlier imagery. When the object is arbitrarily deforming, each, or at least most, contour points can move independently. Contour deformation then forms an infinite (in practice, very large), dimensional space. Direct application of particle filters (PF) for large dimensional problems is impractically expensive. However, in most real problems, at any given time, most of the contour deformation occurs in a small number of dimensions ("effective basis space") while the residual deformation in the rest of the state space ("residual space") is small. This property enables us to apply the particle filtering with mode tracking (PF-MT) idea that was proposed for such large dimensional problems in recent work. Since most contour deformation is low spatial frequency, we propose to use the space of deformation at a subsampled set of locations as the effective basis space. The resulting algorithm is called deform PF-MT. It requires significant modifications compared to the original PF-MT because the space of contours is a non-Euclidean infinite dimensional space.
Toews M, Wells W, Collins LD, Arbel T. Feature-based morphometry: discovering group-related anatomical patterns. Neuroimage. 2010;49 (3) :2318-27.Abstract
This paper presents feature-based morphometry (FBM), a new fully data-driven technique for discovering patterns of group-related anatomical structure in volumetric imagery. In contrast to most morphometry methods which assume one-to-one correspondence between subjects, FBM explicitly aims to identify distinctive anatomical patterns that may only be present in subsets of subjects, due to disease or anatomical variability. The image is modeled as a collage of generic, localized image features that need not be present in all subjects. Scale-space theory is applied to analyze image features at the characteristic scale of underlying anatomical structures, instead of at arbitrary scales such as global or voxel-level. A probabilistic model describes features in terms of their appearance, geometry, and relationship to subject groups, and is automatically learned from a set of subject images and group labels. Features resulting from learning correspond to group-related anatomical structures that can potentially be used as image biomarkers of disease or as a basis for computer-aided diagnosis. The relationship between features and groups is quantified by the likelihood of feature occurrence within a specific group vs. the rest of the population, and feature significance is quantified in terms of the false discovery rate. Experiments validate FBM clinically in the analysis of normal (NC) and Alzheimer's (AD) brain images using the freely available OASIS database. FBM automatically identifies known structural differences between NC and AD subjects in a fully data-driven fashion, and an equal error classification rate of 0.80 is achieved for subjects aged 60-80 years exhibiting mild AD (CDR=1).
Malcolm JG, Michailovich O, Bouix S, Westin C-F, Shenton ME, Rathi Y. A filtered approach to neural tractography using the Watson directional function. Med Image Anal. 2010;14 (1) :58-69.Abstract
We propose a technique to simultaneously estimate the local fiber orientations and perform multi-fiber tractography. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the measured signal or estimated fiber orientation. Further, to overcome noise, many algorithms use a filter as a post-processing step to obtain a smooth trajectory. We formulate fiber tracking as causal estimation: at each step of tracing the fiber, the current estimate of the signal is guided by the previous. To do this, we model the signal as a discrete mixture of Watson directional functions and perform tractography within a filtering framework. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the signal and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides an accurate estimate of the local structure at each point along the fiber. We choose the Watson function since it provides a compact representation of the signal parameterized by the principal diffusion direction and a scaling parameter describing anisotropy, and also allows analytic reconstruction of the oriented diffusion function from those parameters. Using a mixture of two and three components (corresponding to two-fiber and three-fiber models) we demonstrate in synthetic experiments that this approach reduces signal reconstruction error and significantly improves the angular resolution at crossings and branchings. In vivo experiments examine the corpus callosum and internal capsule and confirm the ability to trace through regions known to contain such crossing and branching while providing inherent path regularization.
Savadjiev P, Kindlmann GL, Bouix S, Shenton ME, Westin C-F. Local white matter geometry from diffusion tensor gradients. Neuroimage. 2010;49 (4) :3175-86.Abstract
We introduce a mathematical framework for computing geometrical properties of white matter fibers directly from diffusion tensor fields. The key idea is to isolate the portion of the gradient of the tensor field corresponding to local variation in tensor orientation, and to project it onto a coordinate frame of tensor eigenvectors. The resulting eigenframe-centered representation then makes it possible to define scalar indices (or measures) that describe the local white matter geometry directly from the diffusion tensor field and its gradient, without requiring prior tractography. We derive new scalar indices of (1) fiber dispersion and (2) fiber curving, and we demonstrate them on synthetic and in vivo data. Finally, we illustrate their applicability to a group study on schizophrenia.
Berlow YA, Wells WM, Ellison JM, Sung YH, Renshaw PF, Harper DG. Neuropsychiatric correlates of white matter hyperintensities in Alzheimer's disease. Int J Geriatr Psychiatry. 2010;25 (8) :780-8.Abstract
OBJECTIVE: To investigate the association of behavioral and psychological symptoms of dementia (BPSD) in Alzheimer's disease (AD) and magnetic resonance imaging (MRI) measures of brain atrophy and white matter hyperintensities (WMH). METHODS: Thirty-seven patients with probable AD received the Neuropsychiatric Inventory (NPI), the Mini Mental Status Exam (MMSE), and an MRI scan as part of their initial evaluation at the Outpatient Memory Diagnostic Clinic at McLean Hospital. MRI-based volumetric measurements of whole brain atrophy, hippocampal volumes, and WMH were obtained. Analysis of covariance models, using age as a covariate and the presence of specific BPSD as independent variables, were used to test for differences in whole brain volumes, hippocampal volumes and WMH volumes. RESULTS: Increased WMH were associated with symptoms of anxiety, aberrant motor behavior, and night time disturbance, while symptoms of disinhibition were linked to lower WMH volume. No associations were found for whole brain or hippocampal volumes and BPSD. CONCLUSIONS: These findings suggest that white matter changes are associated with the presence of BPSD in AD.
Lashkari D, Sridharan R, Golland P. Categories and Functional Units: An Infinite Hierarchical Model for Brain Activations. Adv Neural Inf Process Syst. 2010;23 :1252-1260.Abstract
We present a model that describes the structure in the responses of different brain areas to a set of stimuli in terms of stimulus categories (clusters of stimuli) and functional units (clusters of voxels). We assume that voxels within a unit respond similarly to all stimuli from the same category, and design a nonparametric hierarchical model to capture inter-subject variability among the units. The model explicitly encodes the relationship between brain activations and fMRI time courses. A variational inference algorithm derived based on the model learns categories, units, and a set of unit-category activation probabilities from data. When applied to data from an fMRI study of object recognition, the method finds meaningful and consistent clusterings of stimuli into categories and voxels into units.
Lashkari D, Vul E, Kanwisher N, Golland P. Discovering structure in the space of fMRI selectivity profiles. Neuroimage. 2010;50 (3) :1085-98.Abstract
We present a method for discovering patterns of selectivity in fMRI data for experiments with multiple stimuli/tasks. We introduce a representation of the data as profiles of selectivity using linear regression estimates, and employ mixture model density estimation to identify functional systems with distinct types of selectivity. The method characterizes these systems by their selectivity patterns and spatial maps, both estimated simultaneously via the EM algorithm. We demonstrate a corresponding method for group analysis that avoids the need for spatial correspondence among subjects. Consistency of the selectivity profiles across subjects provides a way to assess the validity of the discovered systems. We validate this model in the context of category selectivity in visual cortex, demonstrating good agreement with the findings based on prior hypothesis-driven methods.
Langs G, Golland P, Tie Y, Rigolo L, Golby AJ. Functional Geometry Alignment and Localization of Brain Areas. Adv Neural Inf Process Syst. 2010;1 :1225-1233.Abstract
Matching functional brain regions across individuals is a challenging task, largely due to the variability in their location and extent. It is particularly difficult, but highly relevant, for patients with pathologies such as brain tumors, which can cause substantial reorganization of functional systems. In such cases spatial registration based on anatomical data is only of limited value if the goal is to establish correspondences of functional areas among different individuals, or to localize potentially displaced active regions. Rather than rely on spatial alignment, we propose to perform registration in an alternative space whose geometry is governed by the functional interaction patterns in the brain. We first embed each brain into a functional map that reflects connectivity patterns during a fMRI experiment. The resulting functional maps are then registered, and the obtained correspondences are propagated back to the two brains. In application to a language fMRI experiment, our preliminary results suggest that the proposed method yields improved functional correspondences across subjects. This advantage is pronounced for subjects with tumors that affect the language areas and thus cause spatial reorganization of the functional regions.
Agar NYR, Malcolm JG, Mohan V, Yang HW, Johnson MD, Tannenbaum A, Agar JN, Black PM. Imaging of meningioma progression by matrix-assisted laser desorption ionization time-of-flight mass spectrometry. Anal Chem. 2010;82 (7) :2621-5.Abstract
Often considered benign, meningiomas represent 32% of intracranial tumors with three grades of malignancy defined by the World Health Organization (WHO) histology based classification. Malignant meningiomas are associated with less than 2 years median survival. The inability to predict recurrence and progression of meningiomas induces significant anxiety for patients and limits physicians in implementing prophylactic treatment approaches. This report presents an analytical approach to tissue characterization based on matrix-assisted laser desorption ionization time-of-flight (MALDI TOF) mass spectrometry imaging (MSI) which is introduced in an attempt to develop a reference database for predictive classification of brain tumors. This pilot study was designed to evaluate the potential of such an approach and to begin to address limitations of the current methodology. Five recurrent and progressive meningiomas for which surgical specimens were available from the original and progressed grades were selected and tested against nonprogressive high-grade meningiomas, high-grade gliomas, and nontumor brain specimens. The common profiling approach of data acquisition was compared to imaging and revealed significant benefits in spatially resolved acquisition for improved spectral definition. A preliminary classifier based on the support vector machine showed the ability to distinguish meningioma image spectra from the nontumor brain and from gliomas, a different type of brain tumor, and to enable class imaging of surgical tissue. Although the development of classifiers was shown to be sensitive to data preparation parameters such as recalibration and peak picking criteria, it also suggested the potential for maturing into a predictive algorithm if provided with a larger series of well-defined cases.
Lee JW, Wen PY, Hurwitz S, Black P, Kesari S, Drappatz J, Golby AJ, Wells WM, Warfield SK, Kikinis R, et al. Morphological characteristics of brain tumors causing seizures. Arch Neurol. 2010;67 (3) :336-42.Abstract
OBJECTIVE: To quantify size and localization differences between tumors presenting with seizures vs nonseizure neurological symptoms. DESIGN: Retrospective imaging survey. We performed magnetic resonance imaging-based morphometric analysis and nonparametric mapping in patients with brain tumors. SETTING: University-affiliated teaching hospital. PATIENTS OR OTHER PARTICIPANTS: One hundred twenty-four patients with newly diagnosed supratentorial glial tumors. MAIN OUTCOME MEASURES: Volumetric and mapping methods were used to evaluate differences in size and location of the tumors in patients who presented with seizures as compared with patients who presented with other symptoms. RESULTS: In high-grade gliomas, tumors presenting with seizures were smaller than tumors presenting with other neurological symptoms, whereas in low-grade gliomas, tumors presenting with seizures were larger. Tumor location maps revealed that in high-grade gliomas, deep-seated tumors in the pericallosal regions were more likely to present with nonseizure neurological symptoms. In low-grade gliomas, tumors of the temporal lobe as well as the insular region were more likely to present with seizures. CONCLUSIONS: The influence of size and location of the tumors on their propensity to cause seizures varies with the grade of the tumor. In high-grade gliomas, rapidly growing tumors, particularly those situated in deeper structures, present with non-seizure-related symptoms. In low-grade gliomas, lesions in the temporal lobe or the insula grow large without other symptoms and eventually cause seizures. Quantitative image analysis allows for the mapping of regions in each group that are more or less susceptible to seizures.
Ou W, Nummenmaa A, Ahveninen J, Belliveau JW, Hämäläinen MS, Golland P. Multimodal functional imaging using fMRI-informed regional EEG/MEG source estimation. Neuroimage. 2010;52 (1) :97-108.Abstract
We propose a novel method, fMRI-Informed Regional Estimation (FIRE), which utilizes information from fMRI in E/MEG source reconstruction. FIRE takes advantage of the spatial alignment between the neural and the vascular activities, while allowing for substantial differences in their dynamics. Furthermore, with a region-based approach, FIRE estimates the model parameters for each region independently. Hence, it can be efficiently applied on a dense grid of source locations. The optimization procedure at the core of FIRE is related to the re-weighted minimum-norm algorithms. The weights in the proposed approach are computed from both the current source estimates and fMRI data, leading to robust estimates in the presence of silent sources in either fMRI or E/MEG measurements. We employ a Monte Carlo evaluation procedure to compare the proposed method to several other joint E/MEG-fMRI algorithms. Our results show that FIRE provides the best trade-off in estimation accuracy between the spatial and the temporal accuracy. Analysis using human E/MEG-fMRI data reveals that FIRE significantly reduces the ambiguities in source localization present in the minimum-norm estimates, and that it accurately captures activation timing in adjacent functional regions.
Gholami B, Haddad WM, Tannenbaum AR. Relevance vector machine learning for neonate pain intensity assessment using digital imaging. IEEE Trans Biomed Eng. 2010;57 (6) :1457-66.Abstract
Pain assessment in patients who are unable to verbally communicate is a challenging problem. The fundamental limitations in pain assessment in neonates stem from subjective assessment criteria, rather than quantifiable and measurable data. This often results in poor quality and inconsistent treatment of patient pain management. Recent advancements in pattern recognition techniques using relevance vector machine (RVM) learning techniques can assist medical staff in assessing pain by constantly monitoring the patient and providing the clinician with quantifiable data for pain management. The RVM classification technique is a Bayesian extension of the support vector machine (SVM) algorithm, which achieves comparable performance to SVM while providing posterior probabilities for class memberships and a sparser model. If classes represent "pure" facial expressions (i.e., extreme expressions that an observer can identify with a high degree of confidence), then the posterior probability of the membership of some intermediate facial expression to a class can provide an estimate of the intensity of such an expression. In this paper, we use the RVM classification technique to distinguish pain from nonpain in neonates as well as assess their pain intensity levels. We also correlate our results with the pain intensity assessed by expert and nonexpert human examiners.
Depa M, Sabuncu MR, Holmvang G, Nezafat R, Schmidt EJ, Golland P. Robust Atlas-Based Segmentation of Highly Variable Anatomy: Left Atrium Segmentation. Stat Atlases Comput Models Heart. 2010;6364 :85-94.Abstract
Automatic segmentation of the heart's left atrium offers great benefits for planning and outcome evaluation of atrial ablation procedures. However, the high anatomical variability of the left atrium presents significant challenges for atlas-guided segmentation. In this paper, we demonstrate an automatic method for left atrium segmentation using weighted voting label fusion and a variant of the demons registration algorithm adapted to handle images with different intensity distributions. We achieve accurate automatic segmentation that is robust to the high anatomical variations in the shape of the left atrium in a clinical dataset of MRA images.
Venkataraman A, Kubicki M, Westin C-F, Golland P. Robust Feature Selection in Resting-State fMRI Connectivity Based on Population Studies. Conf Comput Vis Pattern Recognit Workshops. 2010 :63-70.Abstract
We propose an alternative to univariate statistics for identifying population differences in functional connectivity. Our feature selection method is based on a procedure that searches across subsets of the data to isolate a set of robust, predictive functional connections. The metric, known as the Gini Importance, also summarizes multivariate patterns of interaction, which cannot be captured by univariate techniques. We compare the Gini Importance with univariate statistical tests to evaluate functional connectivity changes induced by schizophrenia. Our empirical results indicate that univariate features vary dramatically across subsets of the data and have little classification power. In contrast, relevant features based on Gini Importance are considerably more stable and allow us to accurately predict the diagnosis of a test subject.