Resting state functional MRI data were collected on a 3T Siemens TIM Trio (Erlangen, Germany) scanner. Participants kept their eyes open during the 5-minute resting scan. Echo-planar EPI sequence images (TR = 2,000 ms, TE = 29 ms, flip angle = 75°) were acquired with an 8-channel head coil. Each volume consisted of 33 axial slices (64 × 64 matrix, 3.75 × 3.75 mm2, 3.5 mm thickness, 1 mm gap). Resting state fMRI data were preprocessed using the statistical parametric mapping software (SPM; http://www.fil.ion.ucl.ac.uk/spm) (Friston, 2003) including slice-timing correction, realignment, co-registration and spatial normalization. Images were transformed to the Montreal Neurological Institute (MNI) standard space. The first five scans were discarded to allow for T1 equilibration. The DVARS method (Power et al., 2012) was used to find spike regressors where the root mean square (RMS) head movement exceeded 3 standard deviations. Time courses, with a size of 145 time steps, were orthogonalized with respect to i) linear, quadratic and cubic trends; ii) the 6 realignment parameters and iii) realignment parameters derivatives. The decision to pre-process these nuisances at this point
