Removing inter-subject technical variability in magnetic resonance imaging studies

Magnetic resonance imaging (MRI) intensities are acquired in arbitrary units, making scans non-comparable across sites and between subjects. Intensity normalization is a first step for the improvement of comparability of the images across subjects. However, we show that unwanted inter-scan variability associated with imaging site, scanner effect, and other technical artifacts is still present after standard intensity normalization in large multi-site neuroimaging studies. We propose RAVEL (Removal of Artificial Voxel Effect by Linear regression), a tool to remove residual technical variability after intensity normalization. As proposed by SVA and RUV [Leek and Storey, 2007, 2008, Gagnon-Bartsch and Speed, 2012], two batch effect correction tools largely used in genomics, we decompose the voxel intensities of images registered to a template into a biological component and an unwanted variation component. The unwanted variation component is estimated from a control region obtained from the cerebrospinal fluid (CSF), where intensities are known to be unassociated with disease status and other clinical covariates. We perform a singular value decomposition (SVD) of the control voxels to estimate factors of unwanted variation. We then estimate the unwanted factors using linear regression for every voxel of the brain and take the residuals as the RAVEL-corrected intensities. We assess the performance of RAVEL using T1-weighted (T1-w) images from more than 900 subjects with Alzheimer's disease (AD) and mild cognitive impairment (MCI), as well as healthy controls from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database. We compare RAVEL to two intensity-normalization-only methods: histogram matching and White Stripe. We show that RAVEL performs best at improving the replicability of the brain regions that are empirically found to be most associated with AD, and that these regions are significantly more present in structures impacted by AD (hippocampus, amygdala, parahippocampal gyrus, enthorinal area, and fornix stria terminals). In addition, we show that the RAVEL-corrected intensities have the best performance in distinguishing between MCI subjects and healthy subjects using the mean hippocampal intensity (AUC=67%), a marked improvement compared to results from intensity normalization alone (AUC=63% and 59% for histogram matching and White Stripe, respectively). RAVEL is promising for many other imaging modalities.

[1]  James J. Pekar,et al.  Regionally-specific diffusion tensor imaging in mild cognitive impairment and Alzheimer's disease , 2009, NeuroImage.

[2]  R. Luukkonen,et al.  Is the signal intensity of cerebrospinal fluid constant? Intensity measurements with high and low field magnetic resonance imagers. , 1993, Magnetic resonance imaging.

[3]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

[4]  R. Mayeux,et al.  Molecular drivers and cortical spread of lateral entorhinal cortex dysfunction in preclinical Alzheimer's disease , 2013, Nature Neuroscience.

[5]  C. Jack,et al.  Prediction of AD with MRI-based hippocampal volume in mild cognitive impairment , 1999, Neurology.

[6]  Nick C Fox,et al.  The Alzheimer's disease neuroimaging initiative (ADNI): MRI methods , 2008, Journal of magnetic resonance imaging : JMRI.

[7]  C. Guttmann,et al.  Deep Gray Matter Involvement on Brain MRI Scans Is Associated with Clinical Progression in Multiple Sclerosis , 2009, Journal of neuroimaging : official journal of the American Society of Neuroimaging.

[8]  Arthur W. Toga,et al.  Atlas-based whole brain white matter analysis using large deformation diffeomorphic metric mapping: Application to normal elderly and Alzheimer's disease participants , 2009, NeuroImage.

[9]  Michael Weiner,et al.  Robust atrophy rate measurement in Alzheimer's disease using multi-site serial MRI: Tissue-specific intensity normalization and parameter selection , 2010, NeuroImage.

[10]  C. Jack,et al.  Rates of hippocampal atrophy correlate with change in clinical status in aging and AD , 2000, Neurology.

[11]  Brigitte Landeau,et al.  Using voxel-based morphometry to map the structural changes associated with rapid conversion in MCI: A longitudinal MRI study , 2005, NeuroImage.

[12]  R. Eils,et al.  Comparison of normalization methods for Illumina BeadChip HumanHT-12 v3 , 2010, BMC Genomics.

[13]  Johann A. Gagnon-Bartsch,et al.  Using control genes to correct for unwanted variation in microarray data. , 2012, Biostatistics.

[14]  D. Arnold,et al.  Quantitative Measurement of tissue damage and recovery within new T2w lesions in pediatric- and adult-onset multiple sclerosis , 2015, Multiple sclerosis.

[15]  John Muschelli,et al.  fslr: Connecting the FSL Software with R , 2015, R J..

[16]  Alan C. Evans,et al.  A nonparametric method for automatic correction of intensity nonuniformity in MRI data , 1998, IEEE Transactions on Medical Imaging.

[17]  Clifford R. Jack,et al.  Alzheimer's disease diagnosis in individual subjects using structural MR images: Validation studies , 2008, NeuroImage.

[18]  Peter Dalgaard,et al.  R Development Core Team (2010): R: A language and environment for statistical computing , 2010 .

[19]  D. Arnold,et al.  Normalization of White Matter Intensity on T1‐Weighted Images of Patients with Acquired Central Nervous System Demyelination , 2015, Journal of neuroimaging : official journal of the American Society of Neuroimaging.

[20]  C. Crainiceanu,et al.  Statistical normalization techniques for magnetic resonance imaging , 2014, NeuroImage: Clinical.

[21]  M. Chong,et al.  Neuroimaging Biomarkers in Alzheimer's Disease , 2009 .

[22]  J. Pujol,et al.  Biological significance of iron-related magnetic resonance imaging changes in the brain. , 1992, Archives of neurology.

[23]  Jeffrey T. Leek,et al.  Opinion: Reproducible research can still be wrong: Adopting a prevention approach , 2015, Proceedings of the National Academy of Sciences.

[24]  Bilwaj Gaonkar,et al.  Analytic estimation of statistical significance maps for support vector machine based multi-variate image analysis and classification , 2013, NeuroImage.

[25]  C. Jack,et al.  Usefulness of MRI measures of entorhinal cortex versus hippocampus in AD , 2000, Neurology.

[26]  C. W. Tjoa,et al.  MRI T2 hypointensity of the dentate nucleus is related to ambulatory impairment in multiple sclerosis , 2005, Journal of the Neurological Sciences.

[27]  J. Leek.  svaseq: removing batch effects and other unwanted noise from sequencing data , 2014, bioRxiv.

[28]  R. Bakshi,et al.  Cognitive impairment is associated with subcortical magnetic resonance imaging grey matter T2 hypointensity in multiple sclerosis , 2006, Multiple sclerosis.

[29]  F. Wacker,et al.  A New Method for MRI Intensity Standardization with Application to Lesion Detection in the Brain , 2006 .

[30]  David M. Simcha,et al.  Tackling the widespread and critical impact of batch effects in high-throughput data , 2010, Nature Reviews Genetics.

[31]  C. Jack,et al.  3D maps from multiple MRI illustrate changing atrophy patterns as subjects progress from mild cognitive impairment to Alzheimer's disease. , 2007, Brain : a journal of neurology.

[32]  N. Schuff,et al.  Magnetic resonance imaging of the entorhinal cortex and hippocampus in mild cognitive impairment and Alzheimer's disease , 2001, Journal of neurology, neurosurgery, and psychiatry.

[33]  A. Noordhof,et al.  In the absence of a gold standard , 2010 .

[34]  S. Black,et al.  Beyond the hippocampus , 2001, Neurology.

[35]  D. Sparks,et al.  Amygdala cell loss and atrophy in Alzheimer's disease , 1992, Annals of neurology.

[36]  C. Bottino,et al.  Volumetric MRI Measurements Can Differentiate Alzheimer's Disease, Mild Cognitive Impairment, and Normal Aging , 2002, International Psychogeriatrics.

[37]  Jeffrey T Leek,et al.  A general framework for multiple testing dependence , 2008, Proceedings of the National Academy of Sciences.

[38]  H. Soininen,et al.  Hippocampus and entorhinal cortex in mild cognitive impairment and early AD , 2004, Neurobiology of Aging.

[39]  Peter A. Calabresi,et al.  OASIS is Automated Statistical Inference for Segmentation, with applications to multiple sclerosis lesion segmentation in MRI☆ , 2013, NeuroImage: Clinical.

[40]  Guido Gerig,et al.  Characterizing growth patterns in longitudinal MRI using image contrast , 2014, Medical Imaging.

[41]  Nick C Fox,et al.  Presymptomatic hippocampal atrophy in Alzheimer's disease. A longitudinal MRI study. , 1996, Brain : a journal of neurology.

[42]  J. Trojanowski,et al.  Prediction of MCI to AD conversion, via MRI, CSF biomarkers, and pattern classification , 2011, Neurobiology of Aging.

[43]  Russell T. Shinohara,et al.  Population-wide principal component-based quantification of blood–brain-barrier dynamics in multiple sclerosis , 2011, NeuroImage.

[44]  Otto W. Witte,et al.  Voxel-Based MRI Intensitometry Reveals Extent of Cerebral White Matter Pathology in Amyotrophic Lateral Sclerosis , 2014, PloS one.

[45]  Russell T. Shinohara,et al.  Relating multi-sequence longitudinal intensity profiles and clinical covariates in incident multiple sclerosis lesions , 2015, NeuroImage: Clinical.

[46]  Russell T. Shinohara,et al.  Statistical estimation of T 1 relaxation times using conventional magnetic resonance imaging , 2016, NeuroImage.

[47]  Brian B. Avants,et al.  N4ITK: Improved N3 Bias Correction , 2010, IEEE Transactions on Medical Imaging.

[48]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[49]  R Cameron Craddock,et al.  Disease state prediction from resting state functional connectivity , 2009, Magnetic resonance in medicine.

[50]  John D. Storey,et al.  Capturing Heterogeneity in Gene Expression Studies by Surrogate Variable Analysis , 2007, PLoS genetics.

[51]  Cees Jonker,et al.  Medial temporal lobe atrophy and memory dysfunction as predictors for dementia in subjects with mild cognitive impairment , 1999, Journal of Neurology.

[52]  Jayaram K. Udupa,et al.  New variants of a method of MRI scale standardization , 2000, IEEE Transactions on Medical Imaging.

[53]  Nick C Fox,et al.  Tracking atrophy progression in familial Alzheimer's disease: a serial MRI study , 2006, The Lancet Neurology.

[54]  E. Mori,et al.  Medial temporal structures relate to memory impairment in Alzheimer’s disease: an MRI volumetric study , 1997, Journal of neurology, neurosurgery, and psychiatry.

[55]  Dinggang Shen,et al.  Classifying spatial patterns of brain activity with machine learning methods: Application to lie detection , 2005, NeuroImage.

[56]  Stephen M. Smith,et al.  Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm , 2001, IEEE Transactions on Medical Imaging.

[57]  J. Morris,et al.  Amygdala atrophy is prominent in early Alzheimer's disease and relates to symptom severity , 2011, Psychiatry Research: Neuroimaging.

[58]  D. Rueckert,et al.  Multi-Method Analysis of MRI Images in Early Diagnostics of Alzheimer's Disease , 2011, PloS one.

[59]  Michael I. Miller,et al.  Amygdalar atrophy in symptomatic Alzheimer's disease based on diffeomorphometry: the BIOCARD cohort , 2015, Neurobiology of Aging.

[60]  Volker Schmid,et al.  Working with the DICOM and NIfTI Data Standards in R , 2011 .

[61]  R. Nieuwenhuys,et al.  Neuron loss and shrinkage in the amygdala in Alzheimer's disease , 1994, Neurobiology of Aging.

[62]  D. Louis Collins,et al.  Evaluating intensity normalization on MRIs of human brain with multiple sclerosis , 2011, Medical Image Anal..

[63]  L. Wahlund,et al.  Structural correlates of mild cognitive impairment , 2004, Neurobiology of Aging.

[64]  Hilkka Soininen,et al.  Diffusion tensor imaging and Tract-Based Spatial Statistics in Alzheimer's disease and mild cognitive impairment , 2011, Neurobiology of Aging.

[65]  Mathieu Lemire,et al.  Functional normalization of 450k methylation array data improves replication in large cancer studies , 2014 .

[66]  Stephen W. Scheff,et al.  Volumetric atrophy of the amygdala in Alzheimer's disease , 1991, Neurology.

[67]  L G Nyúl,et al.  On standardizing the MR image intensity scale , 1999, Magnetic resonance in medicine.

[68]  Simon K. Warfield,et al.  Normalization of joint image-intensity statistics in MRI using the Kullback-Leibler divergence , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[69]  S. Black,et al.  Beyond the hippocampus: MRI volumetry confirms widespread limbic atrophy in AD , 2001 .

[70]  D. Delis,et al.  Posterior cingulum white matter disruption and its associations with verbal memory and stroke risk in mild cognitive impairment. , 2012, Journal of Alzheimer's disease : JAD.

[71]  J. Morris,et al.  Profound Loss of Layer II Entorhinal Cortex Neurons Occurs in Very Mild Alzheimer’s Disease , 1996, The Journal of Neuroscience.

[72]  Giovanni B. Frisoni,et al.  Brain morphometry reproducibility in multi-center 3T MRI studies: A comparison of cross-sectional and longitudinal segmentations , 2013, NeuroImage.

[73]  Rainer Goebel,et al.  Combining multivariate voxel selection and support vector machines for mapping and classification of fMRI spatial patterns , 2008, NeuroImage.

[74]  J. Hort,et al.  Amygdalar volume and psychiatric symptoms in Alzheimer's disease: an MRI analysis , 2006, Acta neurologica Scandinavica.

[75]  D. Reich,et al.  Sample-size calculations for short-term proof-of-concept studies of tissue protection and repair in multiple sclerosis lesions via conventional clinical imaging , 2015, Multiple sclerosis.

[76]  Karl J. Friston,et al.  Voxel-Based Morphometry—The Methods , 2000, NeuroImage.

[77]  J. Baron,et al.  Mapping gray matter loss with voxel-based morphometry in mild cognitive impairment , 2002, Neuroreport.

[78]  R. Bakshi,et al.  T2 hypointensity in the deep gray matter of patients with multiple sclerosis: a quantitative magnetic resonance imaging study. , 2002, Archives of neurology.

[79]  A. Madabhushi,et al.  Comparing MR image intensity standardization against tissue characterizability of magnetization transfer ratio imaging , 2006, Journal of magnetic resonance imaging : JMRI.

[80]  Dominik S. Meier,et al.  Time-series modeling of multiple sclerosis disease activity: A promising window on disease progression and repair potential? , 2007, Neurotherapeutics.

[81]  John Quackenbush,et al.  Multiple-laboratory comparison of microarray platforms , 2005, Nature Methods.

[82]  H. Braak,et al.  Alzheimer’s disease: Pathogenesis and prevention , 2012, Alzheimer's & Dementia.

[83]  Brian B. Avants,et al.  Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain , 2008, Medical Image Anal..

[84]  Randolph W. Parks,et al.  Fronto-temporal-lobe atrophy in early-stage Alzheimer's disease identified using an improved detection methodology , 2007, Psychiatry Research: Neuroimaging.

[85]  Gianluca Bontempi,et al.  A comprehensive overview of Infinium HumanMethylation450 data processing , 2013, Briefings Bioinform..