Statistical parametric network analysis of functional connectivity dynamics during a working memory task

Network analysis has become a tool of choice for the study of functional and structural Magnetic Resonance Imaging (MRI) data. Little research, however, has investigated connectivity dynamics in relation to varying cognitive load. In fMRI, correlations among slow (< 0.1 Hz) fluctuations of blood oxygen level dependent (BOLD) signal can be used to construct functional connectivity networks. Using an anatomical parcellation scheme, we produced undirected weighted graphs linking 90 regions of the brain representing major cortical gyri and subcortical nuclei, in a population of healthy adults (n = 43). Topological changes in these networks were investigated under different conditions of a classical working memory task — the N-back paradigm. A mass-univariate approach was adopted to construct statistical parametric networks (SPNs) that reflect significant modifications in functional connectivity between N-back conditions. Our proposed method allowed the extraction of 'lost' and 'gained' functional networks, providing concise graphical summaries of whole-brain network topological changes. Robust estimates of functional networks are obtained by pooling information about edges and vertices over subjects. Graph thresholding is therefore here supplanted by inference. The analysis proceeds by firstly considering changes in weighted cost (i.e. mean between-region correlation) over the different N-back conditions and secondly comparing small-world topological measures integrated over network cost, thereby controlling for differences in mean correlation between conditions. The results are threefold: (i) functional networks in the four conditions were all found to satisfy the small-world property and cost-integrated global and local efficiency levels were approximately preserved across the different experimental conditions; (ii) weighted cost considerably decreased as working memory load increased; and (iii) subject-specific weighted costs significantly predicted behavioral performances on the N-back task (Wald F = 13.39, df1 = 1, df2 = 83, p < 0.001), and therefore conferred predictive validity to functional connectivity strength, as measured by weighted cost. The results were found to be highly sensitive to the frequency band used for the computation of the between-region correlations, with the relationship between weighted cost and behavioral performance being most salient at very low frequencies (0.01–0.03 Hz). These findings are discussed in relation to the integration/specialization functional dichotomy. The pruning of functional networks under increasing cognitive load may permit greater modular specialization, thereby enhancing performance.