Pitfalls and recommended strategies and metrics for suppressing motion artifacts in functional MRI

In resting-state functional magnetic resonance imaging (rs-fMRI), artefactual signals arising from subject motion can dwarf and obfuscate the neuronal activity signal. Typical motion correction approaches involve the generation of nuisance regressors, which are timeseries of non-brain signals regressed out of the fMRI timeseries to yield putatively artifact-free data. Recent work suggests that concatenating all regressors into a single regression model is more effective than the sequential application of individual regressors, which may reintroduce previously removed artifacts. This work compares 18 motion correction pipelines consisting of head motion, independent components analysis, and non-neuronal physiological signal regressors in sequential or concatenated combinations. The pipelines are evaluated on a dataset of cognitively normal individuals with repeat imaging and on datasets of studies of Autism Spectrum Disorder, Major Depressive Disorder, and Parkinsons Disease. Extensive metrics of motion artifact removal are measured, including resting state network recovery, Quality Control-Functional Connectivity (QC-FC) correlation, distance-dependent artifact, network modularity, and test-retest reliability of multiple rs-fMRI analyses. The results reveal limitations in previously proposed metrics, including the QC-FC correlation and modularity quality, and identify more robust motion correction metrics. The results also reveal limitations in the concatenated regression approach, which is outperformed by the sequential regression approach in the test-retest reliability metrics. Finally, pipelines are recommended that perform well based on quantitative and qualitative comparisons across multiple datasets and robust metrics. These new insights and recommendations help address the need for effective motion artifact correction to reduce noise and confounds in rs-fMRI.