The preprocessed task-free functional MRI datasets were parcellated into 90 regions using anatomical templates defined by Tzourio-Mazoyer et al. [45]. A task-free fMRI timeseries was computed for each of the 90 regions by averaging all voxels within each region at each time point in the time series, resulting in 172 time points for each of the 90 anatomical regions of interest. These regional fMRI time series were then used to construct a 90 node whole-brain task-free functional connectivity network for each subject.