A Semi-Automated Workflow Solution for Multimodal Neuroimaging

Total Page:16

File Type:pdf, Size:1020Kb

A Semi-Automated Workflow Solution for Multimodal Neuroimaging Brain Informatics (2016) 3:1–15 DOI 10.1007/s40708-015-0026-y A semi-automated workflow solution for multimodal neuroimaging: application to patients with traumatic brain injury Koon-Pong Wong . Marvin Bergsneider . Thomas C. Glenn . Vladimir Kepe . Jorge R. Barrio . David A. Hovda . Paul M. Vespa . Sung-Cheng Huang Received: 2 November 2015 / Accepted: 16 November 2015 / Published online: 1 December 2015 Ó The Author(s) 2015. This article is published with open access at Springerlink.com Abstract Traumatic brain injury (TBI) is a major cause 15O-water (0–10 min) and 18F-FDDNP (0–6 min) for of mortality and morbidity, placing a significant financial measuring cerebral blood flow in patients with TBI. burden on the healthcare system worldwide. Non-invasive neuroimaging technologies have been playing a pivotal Keywords Traumatic brain injury (TBI) Á Cerebral blood role in the study of TBI, providing important information flow (CBF) Á Magnetic resonance imaging (MRI) Á for surgical planning and patient management. Advances in 15O-water Á 18F-FDDNP Á Positron emission tomography understanding the basic mechanisms and pathophysiology (PET) Á Spatial normalization of the brain following TBI are hindered by a lack of reli- able image analysis methods for accurate quantitative assessment of TBI-induced structural and pathophysiolog- 1 Introduction ical changes seen on anatomical and functional images obtained from multiple imaging modalities. Conventional Traumatic brain injury (TBI) is an important public health region-of-interest (ROI) analysis based on manual labeling and socio-economic problem throughout the world. It is of brain regions is time-consuming and the results could be one of the most common causes of death and long-term inconsistent within and among investigators. In this study, disability in adolescents, young adults, and the elderly. In we propose a workflow solution framework that combined the United States, it was estimated that 1.7 million people the use of non-linear spatial normalization of structural sustain a TBI annually [1]. Of these people, approximately brain images and template-based anatomical labeling to 81 % were treated in and released from emergency automate the ROI analysis process. The proposed workflow departments, about 16 % were hospitalized and discharged, solution is applied to dynamic PET scanning with and approximately 3 % died [1]. However, these numbers underestimate the real prevalence of TBI as they do not K.-P. Wong (&) Á V. Kepe Á J. R. Barrio Á account for those people who did not seek for medical care, D. A. Hovda Á S.-C. Huang had non-fatal (mild or moderate) TBI and presented in Department of Molecular and Medical Pharmacology, David outpatient settings such as physician’s offices, or those who Geffen School of Medicine at UCLA, Los Angeles, CA, USA received medical care from federal, military, or Veterans e-mail: [email protected] Affairs hospitals [1]. M. Bergsneider Á T. C. Glenn Á D. A. Hovda Á P. M. Vespa Non-invasive neuroimaging technologies have been Department of Neurosurgery, David Geffen School of Medicine playing a pivotal role in the study of TBI, providing at UCLA, Los Angeles, CA, USA important information for anatomic localization, surgical P. M. Vespa planning, staging and monitoring the therapeutic responses, Department of Neurology, David Geffen School of Medicine at and predicting the recovery outcomes that could improve UCLA, Los Angeles, CA, USA the survival and change management in patients under acute and chronic conditions. Survivors of TBI typically S.-C. Huang Department of Biomathematics, David Geffen School of live with varying degrees of physical disability and suffer Medicine at UCLA, Los Angeles, CA, USA from significant cognitive deficits (e.g., impaired attention 123 2 K.-P. Wong et al. and poor executive function) and psychological health curves and surfaces as well as the boundaries and neigh- issues (e.g., depression and elevated impulsivity), all of borhoods between structures while allowing a large degree which require long-term or lifelong medical care and of transformation. Examples of this class of algorithms support. Advances in understanding the basic mechanisms include Demons [14], LDDMM [18], DARTEL [19] and pathophysiology of the brain following TBI are available in the SPM software (http://www.fil.ion.ucl.ac. somewhat limited due to lack of reliable image analysis uk/spm/), FNIRT [20] implemented in the FSL software methods that allow accurate quantitative assessment of (http://www.fmrib.ox.ac.uk/fsl/), and symmetric image TBI-induced structural and pathophysiological changes normalization (SyN) [21, 22] implemented in an open seen on anatomical and functional images obtained from source software package ANTs (Advanced Normalization multiple imaging modalities. The use of multimodal neu- Tools, [23]), which was built on an Insight Segmentation roimaging technologies has the advantages to overcome the and Registration ToolKit (ITK) framework (http://www. limitations of any individual imaging modality and to itk.org/). We used SyN in this study for non-linear brain aggregate clinical characteristics and features obtained warping as it can work with different similarity metrics and from different imaging techniques for knowledge mining regularization kernels [23] and has been extensively eval- and for guiding medical diagnosis and decision making [2]. uated with 8 different performance measures using 80 Some state-of-the-art methods of multimodal imaging and manually labeled MR brain images in a recent large-scale their uses in brain research can be found in the following comparative image registration study and was ranked the review articles [3, 4]. While conventional region-of-inter- overall best among 14 non-linear brain warping algorithms est (ROI) analysis enables quantitation of regional changes being assessed [24]. and serves as the basis for comparing data between indi- To address the technical difficulties in analyzing TBI viduals both within and across imaging modalities [5], imaging data, we propose and develop a workflow solution delineation of brain regions through manual ROI drawing framework that combined the use of non-linear brain is labor-intensive and time-consuming, and is also prone to warping of structural MR images and anatomical labeling reproducibility errors [2, 6, 7]. The complexity level of to automatically derive the regional cerebral blood flow ROI analysis increases tremendously by the complicated (CBF) parameters from multi-tracer PET studies. CBF is an nature of TBI that generally involves a combination of important physiological parameter for assessment of brain focal and diffuse injury mechanisms. Depending on the function in normal and pathological conditions. Brain tis- cause and severity of the brain injury, variability of indi- sues depend on CBF for delivery of nutrients and for vidual TBI brains is increased, particularly in the presence removal of metabolic products. Since Kety and Schmidt of focal lesions (e.g., contusion and hemorrhage) and large developed a theory of inert-gas exchange in 1940s [25], deformations within the brain (e.g., swelling and enlarge- many methods have become available for measuring CBF ment/shrinkage of ventricular space). This adds significant in human. Currently, PET imaging with 15O-water is con- difficulties to conduct group-level analysis (such as statis- sidered the gold standard for non-invasive quantification of tical parametric mapping, SPM [8]) and poses technical CBF [26, 27]. Using a hydrophobic tracer, 2-(1-{6-[(2- challenges to perform atlas-based anatomical labeling and [18F]fluoroethyl)(methyl)amino]-2-naphthyl}ethylidene)- ROI analysis [9, 10] with minimal or no human interven- malononitrile (18F-FDDNP), the initial uptake (0–6 min) of tion. The gist of the problem lies in the use of spatial which is perfusion-limited, it has been shown that regional normalization that integrates brain images obtained from perfusion can be inferred from the relative-delivery different modalities for the same individuals to establish a parameter derived by reference-tissue modeling and from one-to-one correspondence mapping between voxels of the early-summed image, and thus represent surrogate individual brains and a standard brain template in a com- indices of CBF [28]. The use of the proposed workflow mon stereotaxic space. solution is illustrated with neuroimaging data obtained A number of non-linear image registration methods have from T1-weighted magnetic resonance (MR) imaging and been proposed for spatial normalization. Some of these dynamic PET scanning of dual tracers (15O-water and methods use a linear combination of trigonometric func- 18F-FDDNP) on six TBI patients under acute condition. tions [11] or polynomials [12] as the transformation model. Because of an implicit assumption of small deformations, this class of methods would fail to normalize images with 2 Materials and methods large deformation, resulting in folding, shearing, and tearing of neighboring structures in the original image upon 2.1 Subjects and study protocol non-linear transformation. More recent research has been geared toward the development of a large deformation The study was approved by the UCLA Institutional Review framework [13–18] which preserves the continuity of Board and was conducted under the auspices of the UCLA 123 A semi-automated workflow solution for multimodal neuroimaging: application to patients 3 Brain Injury Research Center. Six patients with acute TBI 18F-FDDNP PET studies are
Recommended publications
  • A Riemannian Framework for Shape Analysis of Subcortical Brain Structures
    A Riemannian Framework for Shape Analysis of Subcortical Brain Structures A dissertation presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillment of the requirements for the degree Doctor of Philosophy Shuisheng Xie August 2013 © 2013 Shuisheng Xie. All Rights Reserved. 2 This dissertation titled A Riemannian Framework for Shape Analysis of Subcortical Brain Structures by SHUISHENG XIE has been approved for the School of Electrical Engineering and Computer Science and the Russ College of Engineering and Technology by Jundong Liu Associate Professor of Electrical Engineering and Computer Science Dennis Irwin Dean, Russ College of Engineering and Technology 3 Abstract XIE, SHUISHENG, Ph.D., August 2013, Electrical Engineering and Computer Science A Riemannian Framework for Shape Analysis of Subcortical Brain Structures (106 pp.) Director of Dissertation: Jundong Liu Measuring the volumetric and morphological changes of brain structures in MR images can provide a reliable basis for diagnosing and studying Alzheimer’s Disease (AD) and Mild Cognitive Impairment (MCI). Currently, the prevailing morphometric analysis frameworks are voxel-based morphometry (VBM), deformation-based morphometry (DBM) and tensor-based morphometry (TBM), which all involve a low-order spatial normalization procedure to project the computed voxel tissue map or motion field into a standard space. Shape and deformation information in local areas tend to be smoothed or even wiped out during the averaging process. In order to address the drawbacks of VBM/DBM/TBM, a novel local shape-based morphometry (SBM) framework on individual structures is developed in this dissertation, to measure the difference among a collection of anatomical images for individual subjects.
    [Show full text]
  • Spatial Normalization of Diffusion Tensor MRI Using Multiple Channels
    NeuroImage 20 (2003) 1995–2009 www.elsevier.com/locate/ynimg Spatial normalization of diffusion tensor MRI using multiple channels Hae-Jeong Park,a,b Marek Kubicki,a,b Martha E. Shenton,a,b Alexandre Guimond,c Robert W. McCarley,a Stephan E. Maier,d Ron Kikinis,b Ferenc A. Jolesz,d and Carl-Fredrik Westinb,* a Clinical Neuroscience Division, Laboratory of Neuroscience, Boston VA Health Care System-Brockton Division, Department of Psychiatry, Harvard Medical School, Boston, MA 02115, USA. b Surgical Planning Laboratory, MRI Division, Department of Radiology, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA. c Center for Neurological imaging, MRI Division, Department of Radiology, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA. d Department of Radiology, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA. Received 10 March 2003; revised 30 July 2003; accepted 5 August 2003 Abstract Diffusion Tensor MRI (DT-MRI) can provide important in vivo information for the detection of brain abnormalities in diseases characterized by compromised neural connectivity. To quantify diffusion tensor abnormalities based on voxel-based statistical analysis, spatial normalization is required to minimize the anatomical variability between studied brain structures. In this article, we used a multiple input channel registration algorithm based on a demons algorithm and evaluated the spatial normalization of diffusion tensor image in terms of the input information used for registration. Registration was performed on 16 DT-MRI data sets using different combinations of the channels, including a channel of T2-weighted intensity, a channel of the fractional anisotropy, a channel of the difference of the first and second eigenvalues, two channels of the fractional anisotropy and the trace of tensor, three channels of the eigenvalues of the tensor, and the six channel tensor components.
    [Show full text]
  • Handbook of Functional MRI Data Analysis
    Handbook of Functional MRI Data Analysis Functionalmagnetic resonance imaging (fMRI) has become the most popularmethod for imaging brain function. Handbook of Functional MRI Data Analysis provides a comprehensive and practical introduction to the methods used for fMRI data analysis. Using minimal jargon, this book explains the concepts behind processing fMRI data, focusing on the techniques that are most commonly used in the field. This book provides background about the methods employed by common data analysis packages including FSL, SPM, and AFNI. Some of the newest cutting-edge techniques, including pattern classification analysis, connectivity modeling, and resting state network analysis, are also discussed. Readers of this book, whether newcomers to the field or experienced researchers, will obtain a deep and effective knowledge of how to employ fMRI analysis to ask scientific questions and become more sophisticated users of fMRI analysis software. Dr. Russell A. Poldrack is the director of the Imaging Research Center and professor of Psychology and Neurobiology at the University of Texas at Austin. He has published more than 100 articles in the field of cognitive neuroscience, in journals including Science, Nature, Neuron, Nature Neuroscience, and PNAS. He is well known for his writings on how neuroimaging can be used to make inferences about psychologicalfunction, as wellasfor his research using fMRI and other imaging techniques to understand the brain systems that support learning and memory, decision making, and executive function. Dr. Jeanette A. Mumford is a research assistant professor in the Department of Psychology at the University of Texas at Austin. Trained in biostatistics, her research has focused on the development and characterization of new methods for statistical modeling and analysis of fMRI data.
    [Show full text]
  • Functional MRI User's Guide
    Functional MRI User's Guide Michael A. Yassa ● The Division of Psychiatric Neuroimaging ● ● Department of Psychiatry and Behavioral Sciences ● ● The Johns Hopkins School of Medicine ● ● Baltimore, MD ● 1 Document written in OpenOffice.org Writer 2.0 by Sun Microsystems Publication date: June 2005 (1st edition) Online versions available at http://pni.med.jhu.edu/intranet /fmriguide/ Acknowledgments: This document relies heavily on expertise and advice from the following individuals and/or groups: John Ashburner, Karl Friston, and Will Penny (FIL-UCL: London), Kalina Christoff (UBC: Canada), Matthew Brett (MRC-CBU: Cambridge), and Tom Nichols (SPH-UMichigan, Ann Arbor). Some portions of this document are adapted or copied verbatim from other sources, and are referenced as such. Supplemental Reading: Frackowiak RS, Friston K, Frith C, Dolan RJ, Price CJ, Zeki S, Ashburner J, & Perchey G (2004). Human Brain Function, 2nd edition, Elsevier Academic Press, San Diego, CA. Huettel SA, Song AW, McCarthy, G. (2004) Functional Magnetic Resonance Imaging. Sinaur Associates, Sunderland, MA. 2 Table of Contents Magnetic Resonance Physics..............................................................................6 How the MR Signal is Generated.............................................................................6 The BOLD Contrast Mechanism..............................................................................8 Hemodynamic Modeling.........................................................................................10 Signal and Noise
    [Show full text]
  • Manual Computational Anatomy Toolbox - CAT12
    Manual Computational Anatomy Toolbox - CAT12 Quick Start Guide 2 Version Information 4 Introduction and Overview 10 Getting Started 10 Download and Installation 10 Starting the Toolbox 11 Basic VBM Analysis (Overview) 11 Overview of CAT12 Processing 14 CAT12 Major Processing Steps 14 CAT12 Processing Steps in Detail 15 Basic VBM Analysis (detailed description) 17 Preprocessing Data 17 First Module: Segment Data 17 Second Module: Display slices (optionally) 18 Third Module: Estimate Total Intracranial Volume (TIV) 19 Fourth Module: Check sample 19 Fifth Module: Smooth 21 Building the Statistical Model 22 Two-sample T-Test 23 Full Factorial Model (for a 2x2 Anova) 24 Multiple Regression (Linear) 25 Multiple Regression (Polynomial) 26 Full Factorial Model (Interaction) 27 Full Factorial Model (Polynomial Interaction) 28 Estimating the Statistical Model 29 1 Checking for Design Orthogonality 29 Defining Contrasts 31 Special Cases 36 CAT12 for longitudinal data 36 Optional Change of Parameters for Preprocessing 39 Preprocessing of Longitudinal Data 39 Longitudinal Data in One Group 40 Longitudinal Data in Two Groups 42 Longitudinal Data in Two Groups with interaction of covariate by group 44 Adapting the CAT12 workflow for populations such as children 48 Customized Tissue Probability Maps 48 Customized Dartel- or Shooting-template 49 Other variants of computational morphometry 53 Deformation-based morphometry (DBM) 53 Surface-based morphometry (SBM) 54 Region of interest (ROI) analysis 57 Additional information on native, normalized and modulated volumes 58 Naming convention of output files 60 Calling CAT from the UNIX command line 62 Technical information 63 CAT12 Citation 66 References 68 2 Quick Start Guide VBM data ● Segment data using defaults (use Segment Longitudinal Data for longitudinal data).
    [Show full text]