The nuisance of nuisance regression: Spectral misspecification in a common approach to resting-state fMRI preprocessing reintroduces noise and obscures functional connectivity

Recent resting-state functional connectivity fMRI (RS-fcMRI) research has demonstrated that head motion during fMRI acquisition systematically influences connectivity estimates despite bandpass filtering and nuisance regression, which are intended to reduce such nuisance variability. We provide evidence that the effects of head motion and other nuisance signals are poorly controlled when the fMRI time series are bandpass-filtered but the regressors are unfiltered, resulting in the inadvertent reintroduction of nuisance-related variation into frequencies previously suppressed by the bandpass filter, as well as suboptimal correction for noise signals in the frequencies of interest. This is important because many RS-fcMRI studies, including some focusing on motion-related artifacts, have applied this approach. In two cohorts of individuals (n=117 and 22) who completed resting-state fMRI scans, we found that the bandpass-regress approach consistently overestimated functional connectivity across the brain, typically on the order of r=.10-.35, relative to a simultaneous bandpass filtering and nuisance regression approach. Inflated correlations under the bandpass-regress approach were associated with head motion and cardiac artifacts. Furthermore, distance-related differences in the association of head motion and connectivity estimates were much weaker for the simultaneous filtering approach. We recommend that future RS-fcMRI studies ensure that the frequencies of nuisance regressors and fMRI data match prior to nuisance regression, and we advocate a simultaneous bandpass filtering and nuisance regression strategy that better controls nuisance-related variability.

[1]  Joshua Carp,et al.  Optimizing the order of operations for movement scrubbing: Comment on Power et al. , 2013, NeuroImage.

[2]  Marisa O. Hollinshead,et al.  The organization of the human cerebral cortex estimated by intrinsic functional connectivity. , 2011, Journal of neurophysiology.

[3]  Thomas E. Nichols,et al.  Non-white noise in fMRI: Does modelling have an impact? , 2006, NeuroImage.

[4]  V. Calhoun,et al.  Changes in the interaction of resting‐state neural networks from adolescence to adulthood , 2009, Human brain mapping.

[5]  Robert S. Stawski,et al.  Multilevel Analysis: An Introduction to Basic and Advanced Multilevel Modeling (2nd Edition) , 2013 .

[6]  Karsten Mueller,et al.  The General Linear Model , 2015 .

[7]  Thomas T. Liu,et al.  A component based noise correction method (CompCor) for BOLD and perfusion based fMRI , 2007, NeuroImage.

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

[9]  R W Cox,et al.  AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. , 1996, Computers and biomedical research, an international journal.

[10]  Simon B. Eickhoff,et al.  An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data , 2013, NeuroImage.

[11]  Stephen M. Smith,et al.  Advances and Pitfalls in the Analysis and Interpretation of Resting-State FMRI Data , 2010, Front. Syst. Neurosci..

[12]  Abraham Z. Snyder,et al.  Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion , 2012, NeuroImage.

[13]  Christopher L. Asplund,et al.  The organization of the human cerebellum estimated by intrinsic functional connectivity. , 2011, Journal of neurophysiology.

[14]  S. Petersen,et al.  Development of distinct control networks through segregation and integration , 2007, Proceedings of the National Academy of Sciences.

[15]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

[16]  Hang Joon Jo,et al.  Mapping sources of correlation in resting state FMRI, with artifact detection and removal , 2010, NeuroImage.

[17]  D. Bates,et al.  Linear Mixed-Effects Models using 'Eigen' and S4 , 2015 .

[18]  Donald B. Percival,et al.  Spectral Analysis for Physical Applications , 1993 .

[19]  Biyu J. He,et al.  Electrophysiological correlates of the brain's intrinsic large-scale functional architecture , 2008, Proceedings of the National Academy of Sciences.

[20]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[21]  R. Baayen,et al.  Mixed-effects modeling with crossed random effects for subjects and items , 2008 .

[22]  Kaustubh Supekar,et al.  Development of Large-Scale Functional Brain Networks in Children , 2009, NeuroImage.

[23]  Archana Venkataraman,et al.  Intrinsic functional connectivity as a tool for human connectomics: theory, properties, and optimization. , 2010, Journal of neurophysiology.

[24]  S. Petersen,et al.  Control networks in paediatric Tourette syndrome show immature and anomalous patterns of functional connectivity , 2008, Brain : a journal of neurology.

[25]  Jonathan D. Power,et al.  Prediction of Individual Brain Maturity Using fMRI , 2010, Science.

[26]  Peter H. Westfall,et al.  Multiple Testing of General Contrasts Using Logical Constraints and Correlations , 1997 .

[27]  Ewald Moser,et al.  On the origin of respiratory artifacts in BOLD-EPI of the human brain. , 2002, Magnetic resonance imaging.

[28]  Maurizio Corbetta,et al.  The human brain is intrinsically organized into dynamic, anticorrelated functional networks. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[29]  Mert R. Sabuncu,et al.  The influence of head motion on intrinsic functional connectivity MRI , 2012, NeuroImage.

[30]  Vince D. Calhoun,et al.  A review of group ICA for fMRI data and ICA for joint inference of imaging, genetic, and ERP data , 2009, NeuroImage.

[31]  Peter A. Bandettini,et al.  Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI , 2006, NeuroImage.

[32]  Trenton A. Jerde,et al.  True associations between resting fMRI time series based on innovations , 2011, Journal of neural engineering.

[33]  M. Schölvinck,et al.  Neural basis of global resting-state fMRI activity , 2010, Proceedings of the National Academy of Sciences.

[34]  Joshua Carp,et al.  The secret lives of experiments: Methods reporting in the fMRI literature , 2012, NeuroImage.

[35]  B. Miller,et al.  Neurodegenerative Diseases Target Large-Scale Human Brain Networks , 2009, Neuron.

[36]  Peter A. Bandettini,et al.  The respiration response function: The temporal dynamics of fMRI signal fluctuations related to changes in respiration , 2008, NeuroImage.

[37]  R. Engle Band Spectrum Regression , 1974 .

[38]  Mark A. Elliott,et al.  Impact of in-scanner head motion on multiple measures of functional connectivity: Relevance for studies of neurodevelopment in youth , 2012, NeuroImage.

[39]  Kevin Murphy,et al.  The impact of global signal regression on resting state correlations: Are anti-correlated networks introduced? , 2009, NeuroImage.

[40]  T. Hothorn,et al.  Simultaneous Inference in General Parametric Models , 2008, Biometrical journal. Biometrische Zeitschrift.

[41]  Karl J. Friston,et al.  Movement‐Related effects in fMRI time‐series , 1996, Magnetic resonance in medicine.

[42]  D. G. Watts,et al.  Spectral analysis and its applications , 1968 .

[43]  David C. Van Essen,et al.  The future of the human connectome , 2012, NeuroImage.

[44]  G H Glover,et al.  Image‐based method for retrospective correction of physiological motion effects in fMRI: RETROICOR , 2000, Magnetic resonance in medicine.

[45]  A. Dale,et al.  Whole Brain Segmentation Automated Labeling of Neuroanatomical Structures in the Human Brain , 2002, Neuron.

[46]  Andrea Lancichinetti,et al.  Community detection algorithms: a comparative analysis: invited presentation, extended abstract , 2009, VALUETOOLS.

[47]  Michael W. Cole,et al.  Variable Global Dysconnectivity and Individual Differences in Schizophrenia , 2011, Biological Psychiatry.

[48]  Rupert Lanzenberger,et al.  Correlations and anticorrelations in resting-state functional connectivity MRI: A quantitative comparison of preprocessing strategies , 2009, NeuroImage.

[49]  N. Logothetis,et al.  The Amplitude and Timing of the BOLD Signal Reflects the Relationship between Local Field Potential Power at Different Frequencies , 2012, The Journal of Neuroscience.

[50]  Jonathan D. Power,et al.  The Development of Human Functional Brain Networks , 2010, Neuron.

[51]  Catie Chang,et al.  Influence of heart rate on the BOLD signal: The cardiac response function , 2009, NeuroImage.

[52]  Abraham Z. Snyder,et al.  Steps toward optimizing motion artifact removal in functional connectivity MRI; a reply to Carp , 2013, NeuroImage.

[53]  Jacob Cohen Statistical Power Analysis for the Behavioral Sciences , 1969, The SAGE Encyclopedia of Research Design.