The same data preprocessing steps were applied to both the Discovery and the Replication cohort. Briefly, images for each subject were realigned to the first volume in the time series to correct for head motion, spatially normalized into a standard stereotaxic space (Montreal Neurological Institute template) using a 12-parameter affine model (final resolution of functional images=2 mm isotropic voxels), and smoothed to minimize noise and residual difference in gyral anatomy with a Gaussian filter, set at 6-mm full-width at half-maximum. Voxel-wise signal intensities were ratio normalized to the whole-brain global mean.
