Initial removal of signal outliers, heart rate physiologic noise correction, slice timing correction, spatial smoothing, and registration to a T1 anatomical scan were performed. Susceptibility-induced off-resonance field correction was conducted by FSL “topup” software (www.fmrib.ox.ac.uk/fsl/). Quality control for head motion was performed by eliminating fMRI runs which did not meet the Parkes et al. (2018) stringent criteria [28]. Head motion re-alignment and signal correction was performed using the FSL MCFLIRT motion-correction program and ICA-AROMA, respectively [29–31]. Further denoising was performed using the aCompCor procedure implemented in CONN software (www.ni-trc.org/projects/conn, RRID:SCR_009550), and ICA components with possible motion-related structured noise are regressed out during the FSL dual regression procedure [32,33]. The denoised fMRI timeseries was transformed into MNI space. High pass filtering (cutoff=125 s), but not low pass filtering, was performed. Greater detail in the preprocessing steps is presented in Supplementary File.
