Functional magnetic resonance imaging (fMRI) research has generally focused on drawing conclusions from average brain activation patterns. Importantly, the brain is inherently variable; growing literature has found that within-individual blood oxygen level-dependent (BOLD) signal variability may be meaningful, and not just “noise.” For example, recent research has identified increased BOLD signal variability in healthy younger and older adults during more effortful/complex task loads of n-back paradigms, commonly used tasks that involve important elements of executive function (e.g., attention, working memory, planning, inhibition, etc.). Verbal fluency is a complex cognitive domain that also involves similar processes to generate words given certain rules. As a result, the current study builds on existing literature to investigate within-individual BOLD signal variability patterns in peak coordinates of a verbal fluency network during different loads of a letter n-back task. Due to greater executive demands, greater variability was expected during more effortful/complex n-back task loads in regions of a verbal fluency network.

Forty-eight healthy young adults (Mage(SD) = 22.41(4.47), 25 females) from the Atlanta area completed a letter n-back task in an MRI scanner. After standard processing in AFNI, images were corrected for motion and physiological artifacts, which may be confounding sources of variability. Volumes associated with each load of the letter n-back task (0-back, 1-back, 2-back, 3-back, crosshair) were identified. Task runs were normalized and respective run means were subtracted prior to concatenating all runs for each load type. Standard deviations were calculated across this mean-run corrected time series. Ten peak regions of interest (ROIs) were identified from a verbal fluency network generated from 84 peer-reviewed publications for this domain gathered on NeuroSynth. Paired samples t-tests with Benjamini-Hochberg correction for multiple comparisons were conducted to explore differences in variability during n-back task loads.

In several of the verbal fluency network ROIs, within-individual BOLD signal variability was significantly greater for 2-back versus 0-back loads with medium to large effect sizes (p’s < .001 - < .01, Cohen’s d range: .53-.93). Variability was also significantly greater for 3-back versus 0-back loads with small to medium effect sizes (p’s < .001 - < .01, Cohen’s d range: .48-.74). Specific regions that evidenced this pattern included ROIs in the left inferior frontal gyrus, left cingulate, right inferior frontal gyrus, left middle frontal gyrus, and left superior parietal lobule. Only two regions demonstrated increased variability in the 1-back load versus crosshair (left middle frontal gyrus, p < .001, d = .63; left lentiform nucleus, p < .05, d = .42). No regions demonstrated a significant difference in variability in the 0-back load versus crosshair.

This study contributes to growing literature examining within-individual BOLD signal variability in healthy individuals by exploring variability patterns in a verbal fluency network. The observed pattern of results supports the hypothesis and is in line with previous research, demonstrating that greater variability occurs with greater executive task demands. Future research can use an inscanner task of verbal fluency and can extend variability findings during this in-scanner task to out-of-scanner measures of verbal fluency.