SCI Publications

Many mental illnesses are thought to have their origins in early stages of development, encouraging increased research efforts related to early neurodevelopment. Magnetic resonance imaging (MRI) has provided us with an unprecedented view of the brain in vivo. More recently, diffusion tensor imaging (DTI/DT-MRI), a magnetic resonance imaging technique, has enabled the characterization of the microstrucutral organization of tissue in vivo. As the brain develops, the water content in the brain decreases while protein and fat content increases due to processes such as myelination and axonal organization. Changes of signal intensity in structural MRI and diffusion parameters of DTI reflect these underlying biological changes.

Longitudinal neuroimaging studies provide a unique opportunity for understanding brain maturation by taking repeated scans over a time course within individuals. Despite the availability of detailed images of the brain, there has been little progress in accurate modeling of brain development or creating predictive models of structure that could help identify early signs of illness. We have developed methodologies for the nonlinear parametric modeling of longitudinal structural MRI and DTI changes over the neurodevelopmental period to address this gap. This research provides a normative model of early brain growth trajectory as is represented in structural MRI and DTI data, which will be crucial to understanding the timing and potential mechanisms of atypical development. Growth trajectories are described via intuitive parameters related to delay, rate of growth and expected asymptotic values, all descriptive measures that can answer clinical questions related to quantitative analysis of growth patterns. We demonstrate the potential of the framework on two clinical studies: healthy controls (singletons and twins) and children at risk of autism. Our framework is designed not only to provide qualitative comparisons, but also to give researchers and clinicians quantitative parameters and a statistical testing scheme. Moreover, the method includes modeling of growth trajectories of individuals, resulting in personalized profiles. The statistical framework also allows for prediction and prediction intervals for subject-specific growth trajectories, which will be crucial for efforts to improve diagnosis for individuals and personalized treatment.

Keywords: autism, brain development, image analysis

Working memory emerges in infancy and plays a privileged role in subsequent adaptive cognitive development. The neural networks important for the development of working memory during infancy remain unknown. We used diffusion tensor imaging (DTI) and deterministic fiber tracking to characterize the microstructure of white matter fiber bundles hypothesized to support working memory in 12-month-old infants (n=73). Here we show robust associations between infants' visuospatial working memory performance and microstructural characteristics of widespread white matter. Significant associations were found for white matter tracts that connect brain regions known to support working memory in older children and adults (genu, anterior and superior thalamic radiations, anterior cingulum, arcuate fasciculus, and the temporal-parietal segment). Better working memory scores were associated with higher FA and lower RD values in these selected white matter tracts. These tract-specific brain-behavior relationships accounted for a significant amount of individual variation above and beyond infants' gestational age and developmental level, as measured with the Mullen Scales of Early Learning. Working memory was not associated with global measures of brain volume, as expected, and few associations were found between working memory and control white matter tracts. To our knowledge, this study is among the first demonstrations of brain-behavior associations in infants using quantitative tractography. The ability to characterize subtle individual differences in infant brain development associated with complex cognitive functions holds promise for improving our understanding of normative development, biomarkers of risk, experience-dependent learning and neuro-cognitive periods of developmental plasticity.

Background: Increasing evidence suggests that autism is a disorder of distributed neural networks that may exhibit abnormal developmental trajectories. Characterization of white matter early in the developmental course of the disorder is critical to understanding these aberrant trajectories.

Methods: A cross-sectional study of 2-6 year old children with autism was conducted using diffusion tensor imaging combined with a novel statistical approach employing fractional anisotropy distributions. 58 children aged 18-79 months were imaged: 33 were diagnosed with autism, 8 with general developmental delay (DD), and 17 were typically developing (TD). Fractional anisotropy values within global white matter, cortical lobes, and the cerebellum were measured and transformed to random F distributions for each subject. Each distribution of values for a region was summarized by estimating delta, the estimated mean and standard deviation of the approximating F for each distribution.

Results: The estimated delta parameter, delta-hat, was significantly decreased in individuals with autism compared to the combined control group. This was true in all cortical lobes, as well as in the cerebellum, but differences were strongest in the temporal lobe. Predicted developmental trajectories of delta-hat across the age range in the sample showed patterns that partially distinguished the groups. Exploratory analyses suggested that the variability, rather than the central tendency, component of delta-hat was the driving force behind these results. Conclusions: White matter in young children with autism appears to be abnormally homogeneous, which may reflect poorly organized or differentiated pathways, particularly in the temporal lobe, which is important for social and emotional cognition.

The evaluation of analysis methods for diffusion tensor imaging (DTI) remains challenging due to the lack of gold standards and validation frameworks. Significant work remains in developing metrics for comparing fiber bundles generated from streamline tractography. We propose a set of volumetric and tract oriented measures for evaluating tract differences. The different methods developed for this assessment work are: an overlap measurement, a point cloud distance and a quantification of the diffusion properties at similar locations between fiber bundles. The application of the measures in this paper is a comparison of atlas generated tractography to tractography generated in individual images. For the validation we used a database of 37 subject DTIs, and applied the measurements on five specific fiber bundles: uncinate, cingulum (left and right for both bundles) and genu. Each measurments is interesting for specific use: the overlap measure presents a simple and comprehensive metric but is sensitive to partial voluming and does not give consistent values depending on the bundle geometry. The point cloud distance associated with a quantile interpretation of the distribution gives a good intuition of how close and similar the bundles are. Finally, the functional difference is useful for a comparison of the diffusion properties since it is the focus of many DTI analysis to compare scalar invariants. The comparison demonstrated reasonable similarity of results. The tract difference measures are also applicable to comparison of tractography algorithms, quality control, reproducibility studies, and other validation problems.

This paper presents a method for normalization of diffusion tensor images (DTI) to a fixed DTI template, a pre-processing step to improve the performance of full tensor based registration methods. The proposed method maps the individual tensors of the subject image in to the template space based on matching the cumulative distribution function and the fractional anisotrophy values. The method aims to determine a more accurate deformation field from any full tensor registration method by applying the registration algorithm on the normalized DTI rather than the original DTI. The deformation field applied to the original tensor images are compared to the deformed image without normalization for 11 different cases of mapping seven subjects (neonate through 2 years) to two different atlases. The method shows an improvement in DTI registration based on comparing the normalized fractional anisotropy values of major fiber tracts in the brain.

Recent contributions to the body of knowledge on traumatic brain injury (TBI) favor the view that multimodal neuroimaging using structural and functional magnetic resonance imaging (MRI and fMRI, respectively) as well as diffusion tensor imaging (DTI) has excellent potential to identify novel biomarkers and predictors of TBI outcome. This is particularly the case when such methods are appropriately combined with volumetric/morphometric analysis of brain structures and with the exploration of TBI]related changes in brain network properties at the level of the connectome. In this context, our present review summarizes recent developments on the roles of these two techniques in the search for novel structural neuroimaging biomarkers that have TBI outcome prognostication value. The themes being explored cover notable trends in this area of research, including (1) the role of advanced MRI processing methods in the analysis of structural pathology, (2) the use of brain connectomics and network analysis to identify outcome biomarkers, and (3) the application of multivariate statistics to predict outcome using neuroimaging metrics. The goal of the review is to draw the communityfs attention to these recent advances on TBI outcome prediction methods and to encourage the development of new methodologies whereby structural neuroimaging can be used to identify biomarkers of TBI outcome.

Computer simulation and visualization are having a substantial impact on biomedicine and other areas of science and engineering. Advanced simulation and data acquisition techniques allow biomedical researchers to investigate increasingly sophisticated biological function and structure. A continuing trend in all computational science and engineering applications is the increasing size of resulting datasets. This trend is also evident in data acquisition, especially in image acquisition in biology and medical image databases.

For example, in a collaboration between neuroscientist Robert Marc and our research team at the University of Utah's Scientific Computing and Imaging (SCI) Institute (www.sci.utah.edu), we're creating datasets of brain electron microscopy (EM) mosaics that are 16 terabytes in size. However, while there's no foreseeable end to the increase in our ability to produce simulation data or record observational data, our ability to use this data in meaningful ways is inhibited by current data analysis capabilities, which already lag far behind. Indeed, as the NIH-NSF Visualization Research Challenges report notes, to effectively understand and make use of the vast amounts of data researchers are producing is one of the greatest scientific challenges of the 21st century.

Visual data analysis involves creating images that convey salient information about underlying data and processes, enabling the detection and validation of expected results while leading to unexpected discoveries in science. This allows for the validation of new theoretical models, provides comparison between models and datasets, enables quantitative and qualitative querying, improves interpretation of data, and facilitates decision making. Scientists can use visual data analysis systems to explore \"what if\" scenarios, define hypotheses, and examine data under multiple perspectives and assumptions. In addition, they can identify connections between numerous attributes and quantitatively assess the reliability of hypotheses. In essence, visual data analysis is an integral part of scientific problem solving and discovery.

As applied to biomedical systems, visualization plays a crucial role in our ability to comprehend large and complex data-data that, in two, three, or more dimensions, convey insight into many diverse biomedical applications, including understanding neural connectivity within the brain, interpreting bioelectric currents within the heart, characterizing white-matter tracts by diffusion tensor imaging, and understanding morphology differences among different genetic mice phenotypes.

Keywords: kaust

A population growth model that represents the growth trajectories of individual subjects is critical to study and understand neurodevelopment. This paper presents a framework for jointly estimating and modeling individual and population growth trajectories, and determining significant regional differences in growth pattern characteristics applied to longitudinal neuroimaging data. We use non-linear mixed effect modeling where temporal change is modeled by the Gompertz function. The Gompertz function uses intuitive parameters related to delay, rate of change, and expected asymptotic value; all descriptive measures which can answer clinical questions related to growth. Our proposed framework combines nonlinear modeling of individual trajectories, population analysis, and testing for regional differences. We apply this framework to the study of early maturation in white matter regions as measured with diffusion tensor imaging (DTI). Regional differences between anatomical regions of interest that are known to mature differently are analyzed and quantified. Experiments with image data from a large ongoing clinical study show that our framework provides descriptive, quantitative information on growth trajectories that can be directly interpreted by clinicians. To our knowledge, this is the first longitudinal analysis of growth functions to explain the trajectory of early brain maturation as it is represented in DTI.

Objective: Evidence from prospective studies of high-risk infants suggests that early symptoms of autism usually emerge late in the first or early in the second year of life after a period of relatively typical development. The authors prospectively examined white matter fiber tract organization from 6 to 24 months in high-risk infants who developed autism spectrum disorders (ASDs) by 24 months.

Method: The participants were 92 highrisk infant siblings from an ongoing imaging study of autism. All participants had diffusion tensor imaging at 6 months and behavioral assessments at 24 months; a majority contributed additional imaging data at 12 and/or 24 months. At 24 months, 28 infants met criteria for ASDs and 64 infants did not. Microstructural properties of white matter fiber tracts reported to be associated with ASDs or related behaviors were characterized by fractional anisotropy and radial and axial diffusivity.

Results: The fractional anisotropy trajectories for 12 of 15 fiber tracts differed significantly between the infants who developed ASDs and those who did not. Development for most fiber tracts in the infants with ASDs was characterized by higher fractional anisotropy values at 6 months followed by slower change over time relative to infants without ASDs. Thus, by 24 months of age, those with ASDs had lower values.

Conclusions: These results suggest that aberrant development of white matter pathways may precede the manifestation of autistic symptoms in the first year of life. Longitudinal data are critical to characterizing the dynamic age-related brain and behavior changes underlying this neurodevelopmental disorder.

A population growth model that represents the growth trajectories of individual subjects is critical to study and understand neurodevelopment. This paper presents a framework for jointly estimating and modeling individual and population growth trajectories, and determining significant regional differences in growth pattern characteristics applied to longitudinal neuroimaging data. We use non-linear mixed effect modeling where temporal change is modeled by the Gompertz function. The Gompertz function uses intuitive parameters related to delay, rate of change, and expected asymptotic value; all descriptive measures which can answer clinical questions related to growth. Our proposed framework combines nonlinear modeling of individual trajectories, population analysis, and testing for regional differences. We apply this framework to the study of early maturation in white matter regions as measured with diffusion tensor imaging (DTI). Regional differences between anatomical regions of interest that are known to mature differently are analyzed and quantified. Experiments with image data from a large ongoing clinical study show that our framework provides descriptive, quantitative information on growth trajectories that can be directly interpreted by clinicians. To our knowledge, this is the first longitudinal analysis of growth functions to explain the trajectory of early brain maturation as it is represented in DTI.

Keywords: namic

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) has become the surgical therapy of choice for medically intractable Parkinson's disease. However, quantitative understanding of the interaction between the electric field generated by DBS and the underlying neural tissue is limited. Recently, computational models of varying levels of complexity have been used to study the neural response to DBS. The goal of this study was to evaluate the quantitative impact of incrementally incorporating increasing levels of complexity into computer models of STN DBS. Our analysis focused on the direct activation of experimentally measureable fiber pathways within the internal capsule (IC). Our model system was customized to an STN DBS patient and stimulation thresholds for activation of IC axons were calculated with electric field models that ranged from an electrostatic, homogenous, isotropic model to one that explicitly incorporated the voltage-drop and capacitance of the electrode-electrolyte interface, tissue encapsulation of the electrode, and diffusion-tensor based 3D tissue anisotropy and inhomogeneity. The model predictions were compared to experimental IC activation defined from electromyographic (EMG) recordings from eight different muscle groups in the contralateral arm and leg of the STN DBS patient. Coupled evaluation of the model and experimental data showed that the most realistic predictions of axonal thresholds were achieved with the most detailed model. Furthermore, the more simplistic neurostimulation models substantially overestimated the spatial extent of neural activation.

Keywords: Action Potentials, Action Potentials: physiology, Computer Simulation, Deep Brain Stimulation, Deep Brain Stimulation: instrumentation, Deep Brain Stimulation: methods, Humans, Male, Middle Aged, Models, Neurological, Parkinson Disease, Parkinson Disease: therapy, Subthalamic Nucleus, Subthalamic Nucleus: physiology

Diffusion tensor imaging (DTI) provides a unique source of information about the underlying tissue structure of brain white matter in vivo including both the geometry of major fiber bundles as well as quantitative information about tissue properties represented by derived tensor measures. This paper presents a method for statistical comparison of fiber bundle diffusion properties between populations of diffusion tensor images. Unbiased diffeomorphic atlas building is used to compute a normalized coordinate system for populations of diffusion images. The diffeomorphic transformations between each subject and the atlas provide spatial normalization for the comparison of tract statistics. Diffusion properties, such as fractional anisotropy (FA) and tensor norm, along fiber tracts are modeled as multivariate functions of arc length. Hypothesis testing is performed non-parametrically using permutation testing based on the Hotelling T(2) statistic. The linear discriminant embedded in the T(2) metric provides an intuitive, localized interpretation of detected differences. The proposed methodology was tested on two clinical studies of neurodevelopment. In a study of 1 and 2 year old subjects, a significant increase in FA and a correlated decrease in Frobenius norm was found in several tracts. Significant differences in neonates were found in the splenium tract between controls and subjects with isolated mild ventriculomegaly (MVM) demonstrating the potential of this method for clinical studies.

OBJECTIVE: Despite the clinical success of deep brain stimulation (DBS) for the treatment of Parkinson's disease (PD), little is known about the electrical spread of the stimulation. The primary goal of this study was to integrate neuroimaging, neurophysiology and neurostimulation data sets from 10 patients with PD, unilaterally implanted with subthalamic nucleus (STN) DBS electrodes, to identify the theoretical volume of tissue activated (VTA) by clinically defined therapeutic stimulation parameters.

METHODS: Each patient specific model was created with a series of five steps: (1) definition of the neurosurgical stereotactic coordinate system within the context of preoperative imaging data; (2) entry of intraoperative microelectrode recording locations from neurophysiologically defined thalamic, subthalamic and substantia nigra neurons into the context of the imaging data; (3) fitting a three dimensional brain atlas to the neuroanatomy and neurophysiology of the patient; (4) positioning the DBS electrode in the documented stereotactic location, verified by postoperative imaging data; and (5) calculation of the VTA using a diffusion tensor based finite element neurostimulation model.

RESULTS: The patient specific models show that therapeutic benefit was achieved with direct stimulation of a wide range of anatomical structures in the subthalamic region. Interestingly, of the five patients exhibiting a greater than 40\% improvement in their Unified PD Rating Scale (UPDRS), all but one had the majority of their VTA outside the atlas defined borders of the STN. Furthermore, of the five patients with less than 40\% UPDRS improvement, all but one had the majority of their VTA inside the STN.

CONCLUSIONS: Our results are consistent with previous studies suggesting that therapeutic benefit is associated with electrode contacts near the dorsal border of the STN, and provide quantitative estimates of the electrical spread of the stimulation in a clinically relevant context.

Keywords: Brain Mapping, Brain Mapping: methods, Cerebral,Cerebral: physiology, Computer-Assisted, Computer-Assisted: methods, Deep Brain Stimulation, Deep Brain Stimulation: methods, Diffusion Magnetic Resonance Imaging, Diffusion Magnetic Resonance Imaging: methods, Dominance, Electrodes, Humans, Image Processing, Imaging, Implanted, Magnetic Resonance Imaging, Magnetic Resonance Imaging: methods, Nerve Net, Nerve Net: physiopathology, Neurologic Examination, Neurons, Neurons: physiology, Parkinson Disease, Parkinson Disease: physiopathology, Parkinson Disease: therapy, Substantia Nigra, Substantia Nigra: physiopathology, Subthalamic Nucleus, Subthalamic Nucleus: physiopathology, Synaptic Transmission, Synaptic Transmission: physiology, Thalamus, Thalamus: physiopathology, Three-Dimensional, Tomography, Treatment Outcome, X-Ray Computed, X-Ray Computed: methods

Despite the clinical success of deep brain stimulation (DBS) for the treatment of movement disorders, many questions remain about its effects on the nervous system. This study presents a methodology to predict the volume of tissue activated (VTA) by DBS on a patient-specific basis. Our goals were to identify the intersection between the VTA and surrounding anatomical structures and to compare activation of these structures with clinical outcomes. The model system consisted of three fundamental components: (1) a 3D anatomical model of the subcortical nuclei and DBS electrode position in the brain, each derived from magnetic resonance imaging (MRI); (2) a finite element model of the DBS electrode and electric field transmitted to the brain, with tissue conductivity properties derived from diffusion tensor MRI; (3) VTA prediction derived from the response of myelinated axons to the applied electric field, which is a function of the stimulation parameters (contact, impedance, voltage, pulse width, frequency). We used this model system to analyze the effects of subthalamic nucleus (STN) DBS in a patient with Parkinson's disease. Quantitative measurements of bradykinesia, rigidity, and corticospinal tract (CST) motor thresholds were evaluated over a range of stimulation parameter settings. Our model predictions showed good agreement with CST thresholds. Additionally, stimulation through electrode contacts that improved bradykinesia and rigidity generated VTAs that overlapped the zona incerta/fields of Forel (ZI/H2). Application of DBS technology to various neurological disorders has preceded scientific characterization of the volume of tissue directly affected by the stimulation. Synergistic integration of clinical analysis, neuroimaging, neuroanatomy, and neurostimulation modeling provides an opportunity to address wide ranging questions on the factors linked with the therapeutic benefits and side effects of DBS.

Keywords: Deep Brain Stimulation, Finite Element Analysis, Humans, Imaging, Three-Dimensional, Magnetic Resonance Imaging, Models, Neurological, Parkinson Disease, Parkinson Disease: therapy, Pyramidal Tracts, Pyramidal Tracts: physiology, Subthalamic Nucleus, Subthalamic Nucleus: physiology