True associations between resting fMRI time series based on innovations

We calculated voxel-by-voxel pairwise crosscorrelations between prewhitened resting-state BOLD fMRI time series recorded from 60 cortical areas (30 per hemisphere) in 18 human subjects (nine women and nine men). Altogether, more than a billion-and-a-quarter pairs of BOLD time series were analyzed. For each pair, a crosscorrelogram was computed by calculating 21 crosscorrelations, namely at zero lag +/- 10 lags of 2 s duration each. For each crosscorrelogram, in turn, the crosscorrelation with the highest absolute value was found and its sign, value, and lag were retained for further analysis. In addition, the crosscorrelations at zero lag (irrespective of the location of the peak) were also analyzed as a special case. Based on known varying density of anatomical connectivity, we distinguished four general brain groups for which we derived summary statistics of crosscorrelations between voxels within an area (group I), between voxels of paired homotopic areas across the two hemispheres (group II), between voxels of an area and all other voxels in the same (ipsilateral) hemisphere (group III), and voxels of an area and all voxels in the opposite (contralateral) hemisphere (except those in the homotopic area) (group IV). We found the following. (a) Most of the crosscorrelogram peaks occurred at zero lag, followed by +/- 1 lag; (b) over all groups, positive crosscorrelations were much more frequent than negative ones; (c) average crosscorrelation was highest for group I, and decreased progressively for groups II-IV; (d) the ratio of positive over negative crosscorrelations was highest for group I and progressively smaller for groups II-IV; (e) the highest proportion of positive crosscorrelations (with respect to all positive ones) was observed at zero lag; and (f) the highest proportion of negative crosscorrelations (with respect to all negative ones) was observed at lag = 2. These findings reveal a systematic pattern of crosscorrelations with respect to their sign, magnitude, lag and brain group, as defined above. Given that these groups were defined along a qualitative gradient of known overall anatomical connectivity, our results suggest that functional interactions between two voxels may simply reflect the density of such anatomical connectivity between the areas to which the voxels belong.