Activation in Right Dorsolateral Prefrontal Cortex Underlies Stuttering Anticipation

Abstract People who stutter learn to anticipate many of their overt stuttering events. Despite the critical role of anticipation, particularly how responses to anticipation shape stuttering behaviors, the neural bases associated with anticipation are unknown. We used a novel approach to identify anticipated and unanticipated words, which were produced by 22 adult stutterers in a delayed-response task while hemodynamic activity was measured using functional near infrared spectroscopy (fNIRS). Twenty-two control participants were included such that each individualized set of anticipated and unanticipated words was produced by one stutterer and one control participant. We conducted an analysis on the right dorsolateral prefrontal cortex (R-DLPFC) based on converging lines of evidence from the stuttering and cognitive control literatures. We also assessed connectivity between the R-DLPFC and right supramarginal gyrus (R-SMG), two key nodes of the frontoparietal network (FPN), to assess the role of cognitive control, and particularly error-likelihood monitoring, in stuttering anticipation. All analyses focused on the five-second anticipation phase preceding the go signal to produce speech. The results indicate that anticipated words are associated with elevated activation in the R-DLPFC, and that compared to non-stutterers, stutterers exhibit greater activity in the R-DLPFC, irrespective of anticipation. Further, anticipated words are associated with reduced connectivity between the R-DLPFC and R-SMG. These findings highlight the potential roles of the R-DLPFC and the greater FPN as a neural substrate of stuttering anticipation. The results also support previous accounts of error-likelihood monitoring and action-stopping in stuttering anticipation. Overall, this work offers numerous directions for future research with clinical implications for targeted neuromodulation.


INTRODUCTION
Stuttering is a complex neurodevelopmental communication disorder that often negatively impacts social, emotional, and professional opportunities for more than 50 million adults worldwide. The disorder manifests itself to listeners as intermittent interruptions in speech production including part-syllable repetitions and audible and inaudible prolongations of sounds. However, these behaviors do not always accompany stuttering events because most, if not all, stutterers develop the remarkable ability to anticipate stuttering. As a result, they can alter their speech plan prior to execution by, for example, avoiding, stalling, or using a speaking strategy (e.g., pull-out, easy onset; Jackson et al., 2015Jackson et al., , 2019. This potential discrepancy between the internal experience of the speaker (i.e., sensing upcoming speech breakdown) and how stuttering manifests itself to the listener (i.e., as something out of the ordinary) has important clinical implications because how individuals respond to anticipation shapes their communicative experiences. Behavioral and qualitative investigations of anticipation, especially in recent years, have improved our understanding of the phenomenon, but it is critical to augment this evidence with a neural account of stuttering anticipation, especially given the covert nature of anticipation. The purpose of this study was to initiate a brain-based understanding of stuttering anticipation by linking neural activation to self-reported anticipation and subsequent stuttering behaviors.

The Anticipation of Stuttering
Stuttering anticipation refers to the sense or prescience that upcoming speech will be stuttered, should the speaker execute their speech plan as originally intended without alterations (Jackson et al., 2015;Wingate, 1975). Anticipation occurs on a temporal continuum from a longer-term or looming sense of impending stuttering to a shorter-term immediate sense of upcoming stuttering. For example, a speaker may anticipate a word months in advance of saying it (a student knowing at the beginning of the semester that they have to say a certain word in a presentation at the end of the semester); minutes or seconds before (when they are about to introduce themselves); or immediately before executing speech. Anticipation is driven by error-likelihood monitoring whereby the speaker learns associations between "errors" (i.e., stuttered utterances) and listener reactions or other environmental consequences, thereby learning to predict the occurrences of these errors (Arenas, 2012(Arenas, , 2017Garcia-Barrera & Davidow, 2015). While adult stutterers, as a group, predict stuttering with high accuracy in experimental settings (greater than 90% accuracy; Knott et al., 1937;Milisen, 1938;Van Riper, 1936), there is a range in which speakers report anticipating stuttering, from, "sometimes" to "always" (Jackson et al., 2015). Anticipation is a relatively stable feature such that anticipated or feared words are stuttered in experiments even three months after they are identified by participants (Mersov et al., 2018). Arguably most important in the speaker's experience is how they learn or choose to respond to anticipation, whether by avoiding, approaching, or implementing physical speaking strategies that prevent stuttering from coming to the surface (Jackson et al., 2015(Jackson et al., , 2019. In this way, responding to anticipation is mediated by cognitive control.

Cognitive Control
Cognitive control refers to the ability to orchestrate brain functions to complete a given task or reach a certain goal (Miller, 2000). Cognitive control encompasses planning, initiating, and inhibiting actions or tasks, and being flexible and vigilant to tasks in response to environmental demands (Niendam et al., 2012). All of these processes are involved in responding to stuttering anticipation. For example, when a stutterer knows that they are going to stutter, they must initiate (or choose not to initiate) a response which may include avoidance or using a speaking strategy (Jackson et al., 2015(Jackson et al., , 2019; they may inhibit responses due to fear of negative reactions from the listener; they must be flexible with the challenge at hand (i.e., not being able to say what they want to say when they want to say it); and they must remain vigilant to their goal (i.e., producing speech). While numerous studies have examined cognitive control in children and adults who stutter (for review, see Anderson & Ofoe, 2019), no studies have assessed the relationship between cognitive control and anticipation directly.

A Potential Neural Substrate of Stuttering Anticipation
While anticipation is pervasive in the stuttering experience and contributes significantly to the negative impact on quality of life for stutterers (Jackson et al., 2015;Tichenor & Yaruss, 2019), the neural underpinnings of anticipation and related cognitive control processes are unknown. Neurofunctional investigations of stuttering have instead focused on the speech motor network, revealing atypical activation in left perisylvian and motor areas along the arcuate and superior longitudinal fasciculus, and their homologous regions in the right hemisphere, as well as atypical activity in basal ganglia and cerebellum (Braun et al., 1997;S.-E. Chang et al., 2011;De Nil et al., 2000;Fox et al., 2000;Kell et al., 2009;Neef et al., 2016;Neumann et al., 2004;Preibisch et al., 2003;Toyomura et al., 2018). Some of these studies reported significant findings outside of the speech motor network in areas related to cognitive control, even though these areas were not the focus of those investigations. For example, it is widely known that the right dorsolateral prefrontal cortex (R-DLPFC) plays a critical role in cognitive control processes (Koechlin et al., 2003;MacDonald et al., 2000;Miller, 2000;Ridderinkhof et al., 2004). Stutterers exhibit elevated activation in the R-DLPFC (Kell et al., 2009;Lu et al., 2009;Neef et al., 2018), and treatment temporarily reduces activity in R-DLPFC while reducing stuttering symptoms (De Nil et al., 2004;Neumann et al., 2005). Kell et al. (2009) also found that stutterers who reported recovering from stuttering without treatment did not show elevated activation in R-DLPFC, suggesting that these patterns reflect compensatory efforts not learned in therapy (e.g., avoiding, stalling, or using other self-learned speaking strategies). However, because anticipation and the R-DLPFC were not the focus of these studies, the potential relevance of the R-DLPFC to stuttering anticipation, and to stuttering more broadly, is unknown.
We focus primarily on the R-DLPFC in this study, but other regions and networks are likely recruited during anticipation and responding to anticipation. Modern imaging and computational approaches, including task-induced and task-free connectivity analysis, have identified multiple non-overlapping and distributed brain networks that underlie cognitive control, including the frontoparietal network (FPN), salience network, cingulo-opercular network, and dorsal and ventral attention networks (D'Esposito & Postle, 2002;Menon & D'Esposito, 2021;Niendam et al., 2012). The FPN, which includes the R-DLPFC and right supramarginal gyrus (R-SMG) (Menon & D'Esposito, 2021), is particularly relevant because it co-activates with multiple other networks (e.g., salience network, cingulo-opercular network) to carry out the diverse processes associated with cognitive control (Marek & Dosenbach, 2018). The FPN initiates and flexibly modulates interactions between the salience and cinguloopercular networks (Marek & Dosenbach, 2018), which both include the anterior cingulate cortex (ACC). ACC underlies error-likelihood monitoring (Brown & Braver, 2005): DLPFC co-activates with ACC, whereby ACC underlies the detection of errors in response to unintended outcomes and generates error signals, and DLPFC holds task-relevant information in working memory and initiates subsequent actions (Alexander & Brown, 2015;Holroyd & Yeung, 2012). It is reasonable to propose that the ACC underlies stuttering anticipation-the recognition of the breakdown or "glitch" in speech-language planning-and reasonable to predict that the R-DLPFC underlies initiating a response to this breakdown. Given the strong bidirectional connections within the FPN (Goldman-Rakic, 1988;Menon & D'Esposito, 2021;Mesulam, 1998), particularly between R-DLPFC and R-SMG, it is also reasonable to hypothesize that anticipation destabilizes these connections, resulting in altered connectivity.
While it was not possible in the current design to disentangle anticipation and responding to anticipation, the reader should be aware of these somewhat distinct but overlapping processes. We conceive of anticipation as an event-the point in time at which the speaker becomes aware that should they proceed as planned, they will overtly stutter-and responding to anticipation as the constellation of processes, involving the cognitive control system, that underlies how a speaker chooses to proceed in light of the knowledge that they are likely about to stutter. In the current design and given the hemodynamic lag associated with fNIRS, it was not possible to disentangle anticipation and responding to anticipation, and thus, the current study examined the processes involved in anticipation generally. Techniques with better temporal resolution than functional near-infrared spectroscopy (fNIRS; i.e., magnetoencephalography/electroencephalography (MEG/EEG)) are more suited to disentangle these two processes.
The current study examined the relationship between R-DLPFC and stuttering anticipation, which may clarify the significance of previous and seemingly incidental findings of elevated activation in R-DLPFC in stutterers (e.g., Kell et al., 2009;Neef et al., 2018). In a first visit, we used a clinical interview to determine individual-specific anticipated and unanticipated words (Jackson et al., 2020). In a second visit, which occurred between three and 10 days after the first visit, we used fNIRS to measure cortical activation immediately prior to participants producing the anticipated and unanticipated words in a delayed-response task. We focused on superficial cortical structures in part due to imaging depth restrictions associated with fNIRS. While we were not able to measure activation from deeper structures such as ACC, fNIRS offered several advantages compared to other techniques (e.g., functional magnetic resonance imaging (fMRI), EEG, MEG), including (1) robustness to speech movement artifact, and (2) allowing participants to produce speech while they sat upright and across from a communicative partner, which increased the likelihood of anticipation. fNIRS has also been validated as a tool to measure DLPFC activation associated with anticipation (Vassena et al., 2019). A matched control group was included to test whether stutterers recruit R-DLPFC differently than non-stutterers. Each anticipated and unanticipated word list was produced by a stutterer and a control speaker. We conducted a region of interest (ROI) analysis of R-DLPFC and hypothesized that (1) anticipation would be associated with greater activation in R-DLPFC, reflecting cognitive control processes associated with responding to anticipation, and (2) stutterers would exhibit greater activation than control speakers in R-DLPFC during this same time period. We also assessed functional connectivity to test whether anticipation was associated with reduced intrinsic connectivity within the FPN, specifically between R-DLPFC and R-SMG.

MATERIALS AND METHODS
This study was approved by the Institutional Review Boards at New York University and Yale University. Consent was obtained for all participants in accordance with the Declaration of Helsinki.

Participants
Twenty-seven adult stutterers were recruited through the first author's clinical network, mass emails distributed by Friends: The National Association of Young People Who Stutter and the National Stuttering Association, and by word of mouth. After the fNIRS screening (see below), 22 stutterers (9 female; mean age = 31.9, SD = 9.1; three left-handed) and 22 control speakers (10 female; mean = 27.4, SD = 8.0; three left-handed) participated in the study. Control participants were recruited after the stuttering participants so that they could be matched for age, gender, and stimuli (see below). Male-to-female ratio was lower than what is typically observed in the stuttering literature (59% vs. ∼75-80%). All participants were between the ages of 18 and 50, reported that American English was their primary language (multilingual was acceptable as long as English was learned during early childhood [younger than 6 years of age]), and reported negative histories of neurological, speech-language, psychological, learning, and hearing impairment. Participant characteristics, including age, gender, treatment history, and extent score from the Stuttering Anticipation Scale (SAS; Jackson et al., 2018) are included in Table 1. Note. Treatment history descriptions written as reported by participants. TS = Trials Stuttered; SAS = Stuttering Anticipation Scale (extent score out of 100, 0 = never, 100 = always). ID = participant; M = mean; SD = standard deviation.

First visit (stutterers only)
The first visit comprised (1) the stuttering assessment, (2) the clinical interview to determine participant-specific stimuli, and (3) the screening for fNIRS. Only the stuttering group participated in two visits; the fNIRS screening for the control group took place on the same day as fNIRS testing.

Diagnostic assessment
The stuttering group was identical to that in Jackson et al. (2020). That study validated the clinical interview used here, and the only overlapping data between that study and the current study are the speech classification data (stuttered/ambiguous/fluent, interrater reliability). Stuttering was diagnosed by the first author, an American Speech-Language-Hearing Association certified speech-language pathologist (SLP) with more than 10 years of expertise in stuttering intervention.
All stuttering participants (1) self-reported as a person who stutters; and (2) exhibited three or more stuttering-like disfluencies (Yairi & Ambrose, 1992) with temporally aligned physical concomitants (e.g., eye blinking, head movements) during a 5-10 min conversation. Participants also completed the SAS (Jackson et al., 2018), which provided self-report ratings of the extent of anticipation based on a 0-100 (never-always) visual analog scale ("How much do you anticipate stuttering?").

Clinical interview
The clinical interview was described previously in Jackson et al. (2020), and will be described briefly here. Interviews were conducted at New York University and Yale University, and for two participants remotely. The purpose of the interview was to identify 10 words that the participants anticipated they would stutter (hence: anticipated words) and 10 words that they did not anticipate they would stutter (hence: unanticipated words), resulting in individual lists for each participant. This approach extends previous methods for identifying anticipated/unanticipated words (Bowers et al., 2012;Mersov et al., 2018;Wymbs et al., 2013), by including clinical inference (e.g., asking participants whether stuttered words/sounds not immediately identified by participants as anticipated or "feared" words/sounds should be included as anticipated words), as well as using counseling techniques to create an environment in which participants were comfortable to identify/reveal feared words. The words were used to create the stimuli for fNIRS testing, which consisted of short questions or sentence completions that would require the participant to produce the words. Anticipated and unanticipated words were matched for length (number of syllables). Each stutterer was matched with a control participant who produced the same set of words.

fNIRS screening
The goal of the fNIRS screening was to determine whether a reliable hemodynamic signal could be acquired from each participant. This is because factors such as bone density and skull thickness weaken fNIRS signals (Krall & Dawson-Hughes, 1993;Okada & Delpy, 2003). Screening fNIRS participants limits acquisition of invalid data . The screening consisted of a finger tapping task during which participants were required to tap their fingers in an alternating pattern for 15 s then rest for 15 s, for a total of 3 min (right hand). Typically, this task elicits a robust response in left motor/premotor areas. Here, a response was determined to be reliable if there was sufficient separation between oxyhemoglobin (HbO) and deoxyhemoglobin (HbR) signals in any channel in left motor/premotor areas based on visual inspection of the event-triggered average. Importantly, potential participants were excluded before the study began. That is, no participants were excluded from the study after their signals were determined to be reliable. Twenty-seven stutterers were screened, and 22 participated in the study; twenty-eight controls were screened, and 22 participated in the study.

Second visit
The second visit (for the stuttering group) occurred between 3 and 10 days after the first visit, and included fNIRS testing in the Brain Function Laboratory at Yale. Participant-specific stimuli were created between the first and second visits based on the anticipated and unanticipated word lists established during the clinical interview (first visit). Each word list was used to create the stimuli for one stuttering and one control speaker. Stimuli included simple questions or sentence completions that required one-word responses (e.g., "You can fly in an ______" ➔ airplane. "What month comes after June?" ➔ July). All stimulus questions were approximately between 1 and 3 s. Participants were exposed to the stimuli before the experiment to minimize the potential impact of language retrieval and formulation processes (i.e., they became familiar with the questions and answers prior to the experiment). Figure 1 depicts the task timeline. The question was presented auditorily while participants viewed a cross on the monitor. Verbal responses were delayed by 5 s. Participants responded when the screen turned green. Participants were asked to look straight ahead and to try to remain still. The paradigm included interactive and alone conditions, but the condition contrast was not the focus of the current study. The two conditions were pooled in all analyses. During the interactive condition, participants responded to questions asked by the examiner, who was seated directly across from the participant and in full view. Participants were instructed to look at the examiner while responding (e.g., "make sure you're looking at me when you respond"). The experimenter was given a cue on a separate monitor (unseen by the participant) just before the go signal was presented, allowing him to look at the participants when they responded. During the alone condition, participants responded to prerecorded stimuli, the same questions as asked by the examiner as described above, while alone in the testing room. Questions during the interactive condition were matched, to the best of Figure 1. Task timeline. The question/sentence completion was presented auditorily. The anticipation period was 5 s; participants looked at a fixation cross during this time. The green screen signaled participants to produce the word. Participants again looked at the cross during the rest period. the examiner's ability, to the prerecorded questions in terms of duration and prosody. The paradigm included a total of 80 trials: 4 interactive runs, 4 alone runs; each run included 10 words. Each word was produced 4 times.

Data Acquisition
Behavioral Speech data were acquired via acoustic and video recordings. Acoustic signals were recorded using a head-mounted microphone with a pop-screen filter set at the same fixed distance for each participant. Video was captured using a Logitech c920 HD 1080p video camera mounted on the participant's monitor.

Neural
Data collection methods have been described previously (e.g., Hirsch et al., 2018Hirsch et al., , 2021, and are also described here. Hemodynamic signals were acquired using an 80-fiber continuouswave fNIRS system (Shimadzu LABNIRS, Kyoto, Japan) with a temporal resolution of 27 ms. Forty emitters and forty detectors were arranged in a 134-channel layout covering bilateral frontal, temporal, parietal, and occipital lobes. Depending on the size of the participant's head, caps with optode distances of either 2.75 cm or 3 cm were used. Three wavelengths of light (780, 805, and 830 nm) were delivered by each LABNIRS emitter. Absorption was converted to concentration changes for deoxyhemoglobin (HbR) and oxyhemoglobin (HbO) using the Beer-Lambert Law (Matcher et al., 1995). After the experiment, anatomical locations of optodes were determined based on standard head landmarks (inion, nasion, top center [Cz], and left and right tragi) using a Patriot 3D Digitizer (Polhemus, Colchester, VT) and linear transform techniques (Eggebrecht et al., 2012;Ferradal et al., 2014;Okamoto & Dan, 2005). Montreal Neurological Institute (MNI) coordinates for the channels were obtained using NIRS-SPM (Ye et al., 2009) in MATLAB (Mathworks, Natick, MA), and corresponding anatomical locations were determined for each channel. See Table S1 for group median coordinates, atlas-based probabilities, and anatomical regions for each channel; see Figure S1 for a visual representation. (Supporting Information can be found at https://doi.org/10.1162/nol_a _00073.) Channels were clustered into anatomical regions based on shared anatomy. The average number of channels per region was 2.69 ± 1.40.

Behavioral
Errors, which comprised non-productions or incorrect productions due to participants forgetting answers or producing erroneous speech, were not included in any analyses (Table 2). Reaction time was calculated as the time between the go signal (i.e., the green screen) and speech onset as defined by the first articulatory movement or accessory behavior (Table 3). Articulatory onset was marked as the first articulatory movement based on visual inspection using Davinci Resolve (Black Magic Design, Australia), which allowed for frame-by-frame scanning (29.97 frames per second) of the recordings of participants' faces. Interrater reliability between the first author and a SLP blind to the study yielded a Cohen's weighted kappa of 0.89 (p < 0.05), indicating strong agreement (McHugh, 2012). We used articulatory onset because (1) inaudible sound prolongations (blocks) typically included observable movement such as posturing; and (2) it appeared that for some participants, video and audio were not synchronized due to technical error. Although determining neural correlates of stuttered speech was not the primary goal of this study, we include a comparison of stuttered and fluent speech for completeness. A three-point rating system was used to classify stuttering response type: 0 indicated unambiguous fluency; 1 indicated ambiguity (unclear whether stuttered or fluent); and 2 indicated unambiguous stuttering (Jackson et al., 2020).

Neural
Data processing was similar to that previously reported (Hirsch et al., , 2018Zhang et al., 2017). Baseline drift was removed using a NIRS-SPM detrending procedure. Global components were removed using a principal component analysis spatial filter (Zhang et al., 2016), which is comparable to using short-source channels . Channels were rejected if the root mean square (RMS ) of the raw signal was 10 times greater than the group mean RMS, which resulted in a mean of 1.48 of 134 (1.1%) channels per participant rejected. Little to no motion artifact was observed in the data, based on visual inspection, likely because participants were instructed to remain still during the task. This is typical for compliant adult participants . Further, there was no speech movement observed during the anticipation phase. HbR and HbO signals were acquired. The fNIRS data sets for each subject were reshaped into 3D volume images for the general linear model (GLM) analysis using SPM8 ( Wellcome Trust, London, UK). Beta value coordinates were converted to standard MNI space using NIRS-SPM (Ye et al., 2009). Contrast images were rendered on a standardized MNI brain template using a p-value threshold of 0.05 and cluster size threshold of 50 voxels, for visual representation. Anatomical locations of peak voxel activity were identified using the Brodmann area Talairach atlas .

Behavioral
To examine reaction time and stuttering response type (stuttered/ambiguous/fluent), linear mixed effects models were fit using the lme4 (Bates et al., 2014) and lmerTest packages in R (Kuznetsova et al., 2017;R Core Team, 2014). The MuMIn package (Barton, 2020) was used to calculate estimated R 2 for model fit. For reaction time, the model included word type (anticipated/unanticipated), word length (number of syllables), trial, and stuttering response type as fixed factors, and participant as a random factor to account for expected variation due to individual differences. Wilcoxon rank sum tests were also used to assess reaction time differences for stuttered, ambiguous, and fluent speech. For stuttering response type, the model included word type, word length, and trial as fixed factors, and participant as a random factor.

Neural
All fNIRS analyses followed standard voxel-wise GLM techniques (Friston et al., 1994(Friston et al., , 1995 adapted for fNIRS (e.g., Descorbeth et al., 2020;Hirsch et al., 2021). Analyses targeted the anticipation phase, i.e., the 5 s time window between the end of the question and the go cue (see Figure 1). All included voxels were within 2 cm from the cortical surface. The primary ROI analyses focused on the R-DLPFC; secondary ROI analyses examined activation in the right inferior frontal gyrus (R-IFG) and right pre-supplementary motor area (R-preSMA). Connectivity between the R-DLPFC and R-SMG was also assessed.
ROI analyses. The mask for R-DLPFC was created by generating a 10 mm sphere using the MarsBar toolbox (Brett et al., 2002) and xyz coordinates [50 26 38] from Kell et al. (2009). Five ROI analyses of R-DLPFC were conducted.
(1) Anticipated vs. Unanticipated: Activation associated with anticipated vs. unanticipated words was compared by using SPM8 to convolve a 5 s block regressor during the anticipation phase, with a standard hemodynamic response function (HRF) that was fitted to the data. The first-level analysis yielded two beta values for each participant in each run (i.e., 10 trials per run), for anticipated and unanticipated words (five anticipated, five unanticipated). The results of the second-level contrast of anticipated vs. unanticipated words were projected onto a standardized MNI template image using SPM (p < 0.05 and cluster size of >50 voxels, uncorrected). Note that the "whole-brain" image was used only to test the ROIs; however, wholebrain results for this analysis, and all analyses below, are included as supplementary material. Significance for the ROI analysis was tested using a one-tailed t test in SPM (p < 0.05) by determining overlap between the second-level image (anticipated vs. unanticipated) and the mask. Two additional analyses were conducted to test whether the R-DLPFC ROI was significantly activated during the anticipation phase for anticipated and unanticipated words (separately). The second-level results, which compared activation for both anticipated and unanticipated words to rest, were projected onto the same MNI template, and the ROI analyses were carried out as described above.
(2) Stuttered vs. Fluent: To compare activation related to response type within stutterers (stuttered, ambiguous, fluent), the regressor was modulated in height (0 = fluent, 1 = ambiguous, 2 = stuttered ) in the first-level analysis. The assumption was that stuttered speech yielded more (or less) activity than ambiguous responses, which yielded more (or less) activity than fluent responses, effectively providing a "contrast" of stuttered and fluent speech. Beta values were compared to zero using a one-tailed t test in SPM. The ROI analysis was conducted as described above.
(3) Stutterers vs. Controls: The ROI group level comparison was similar to the anticipated vs. unanticipated contrast described above, except for the second-level contrast. The 5 s block regressor was convolved with the HRF, irrespective of whether the word was anticipated. Results of the second-level contrast for stutterers vs. controls were projected onto a standardized MNI template image, and the ROI analysis was carried out as described above.
(4) Task vs. Rest (controls only): We also conducted an analysis of the controls only, during the same 5 s window. The purpose of this analysis was to determine whether the R-DLPFC ROI was activated prior to speech execution, which could indicate whether R-DLPFC is involved in speech motor planning in unimpaired speakers. The second-level results, which compared activation of the controls to rest, were projected onto the same MNI template, and the ROI analysis was carried out as described above.
(5) Interactive vs. Alone: Finally, we compared activation between the interactive and alone conditions to justify our decision to exclude "condition" as a factor from all of the models. Confirmation of the null hypothesis would suggest that condition is not contributing to the anticipated and unanticipated responses.
We completed two additional ROI analyses. First, we used two control ROIs to confirm that our findings in R-DLPFC were specific to stuttering anticipation, and not, for example, general to stuttering or due to systemic artifact. These included L-DLPFC, the homologue of the R-DLPFC ROI, and right precentral gyrus (R-preCG) from Belyk et al. (2017), the most recent activation-likelihood meta-analysis of state stuttering (i.e., stuttered vs. fluent). ROIs were 10 mm spheres with centroids [−50 26 38] for L-DLPFC and [54 −14 34] for R-preCG. Second, we tested two key superficial cortical nodes of the action-stopping network (i.e., R-IFG and R-preSMA), as it has been proposed that action-stopping is associated with stuttering anticipation (Arenas, 2017;Hannah & Aron, 2021;Neef et al., 2018). Centroid coordinates were obtained from recent meta-analyses of stuttering: [46 23 −5] for R-IFG, included as a "state" finding in Belyk et al. (2015); and [15 13 59] for R-preSMA, included as a "trait" finding in Budde et al. (2014). Trait findings were used for R-preSMA because state findings were not reported in either meta-analysis. The coordinates were used to create 10 mm ROI spheres.
Functional connectivity. Psychophysiological Interaction (PPI) analysis (Friston et al., 1994(Friston et al., , 2003 was used to examine functional connectivity between R-DLPFC and R-SMG. PPI computations were performed on residual components of the modeled task (anticipated vs. unanticipated), after the GLM was effectively removed. Significant correlations are thought to reflect dynamic neural coupling, though not necessarily related to the task. PPI analysis was conducted using the gPPI toolbox (McLaren et al., 2012) with SPM8. The PPI analysis can be described by the following equations: in which the HRF is represented by H, and H(x) is the convolution of signal X with kernel H. The demeaned time course is represented by g p where 1 is task time and −1 is rest. β i is the PPI beta value, whereas β p is the beta value for the task, β k is the beta value of the time course of the seed, and Y k represents the fNIRS data collected at the seed region. Here, k is the functionally defined cluster for the seed region based on the GLM. x a represents the estimated neural activity for the seed region; the residual error is represented by e i . R-DLPFC had ROI coordinates [50 26 38], and R-SMG coordinates were determined using the atlas from NIRS-SPM (Rorden & Brett, 2000) and creating a 10 mm sphere that was projected onto the brain surface. Each area was used both as a seed and a target, resulting in two comparisons in total. Two-tailed t tests were used to compare residual activity for anticipated versus unanticipated words, and the Holm-Bonferroni method was used to correct for familywise error rates.

Behavioral
Two linear mixed effects models (not included in Jackson et al., 2020) were fit to examine reaction time and response type. Mean reaction time across all stuttering participants was 348.27 ms (SD = 132.34 ms). See Table 2. Reaction time was significantly impacted by word type (anticipated/unanticipated), such that anticipated words had longer reaction times (β = 11.45, t = 2.13, p < 0.05). Reaction time was also impacted by trial (β = −0.34, t = −2.91, p < 0.01), such that reaction time decreased as the experiment progressed. Reaction time was not impacted by word length (syllables) (β = −0.07, t = −0.73, p > 0.05) or stuttering response type (β = −5.45, t = −1.54, p > 0.05), indicating that reaction time for stuttered trials was not significantly longer than that for fluent trials. R 2 for the reaction time model was 0.24. In addition, post hoc Wilcoxon rank sum tests did not reveal differences in reaction time between stuttered and ambiguous trials ( W = 203, p > 1.10), ambiguous and fluent trials ( W = 190, p > 0.10), or stuttered and fluent trials ( W = 152, p > 0.10). Interrater reliability for reaction time between the first author and a SLP blind to the study yielded a Cohen's weighted kappa of 0.79 (p < 0.05), indicating moderate to strong agreement (McHugh, 2012). Table 3 shows the amount of stuttered, ambiguous, and fluent trials for each participant. For all trials, 43.6% were stuttered (2), 43.3% were fluent (0), and 13.1% were ambiguous (1), or not unambiguously stuttered or fluent (reported in Jackson et al., 2019; see Table 3). Interrater reliability between the first author and an SLP with 8 years of experience (blind to the study) yielded a Cohen's weighted kappa of 0.85 (p < 0.05), indicating high agreement. The remaining response type data were not reported in Jackson et al. (2020). 53.9% of anticipated words were unambiguously stuttered and 33.4% of unanticipated words were unambiguously stuttered. 35.4% of anticipated words were unambiguously fluent whereas 51.2% of unanticipated words were unambiguously fluent. In addition, 10.7% of anticipated words were ambiguous and 15.41% of unanticipated words were ambiguous. It is important to note that words characterized as ambiguous using the Jackson et al. (2020) approach would most likely have been categorized as fluent with a standard binary stuttered/fluent distinction that is most commonly applied clinically. There was more stuttering for anticipated than unanticipated words, as expected (β = −0.36, t = −9.76, p < 0.001). There was also more stuttering for

ROI analyses
(1) Anticipated vs. Unanticipated: The primary contrast compared activation in the R-DLPFC mask associated with anticipated vs. unanticipated words. Anticipated words were associated with greater activation than unanticipated words (p = 0.0217) based on HbR, but not HbO (p > 0.05). Figure 2 includes the ROI illustration (A) and event-related averages for anticipated vs. unanticipated words, for HbO (B) and HbR (C), respectively. The HbR difference between event-related averages at 5 s (i.e., the end of the anticipation phase) for anticipated vs. unanticipated words approached significance (t = 1.67, p = 0.0567), though for HbO, signals were similar (t = 0.57, p > 0.05). Note that the remaining analyses in this section focus on HbR due to the null HbO results, as well as evidence suggesting that HbR (1) is more strongly correlated to the blood oxygen level dependent signal (Boas et al., 2004(Boas et al., , 2014Ferrari & Quaresima, 2012); (2) has greater spatial specificity than HbO Zhang et al., 2016); (3) has been validated for speech-language tasks ; and (4) is less sensitive to systemic effects (e.g., heart rate, breathing, bloodflow) than HbO (Franceschini et al., 2003;Kirilina et al., 2012;Santosa et al., 2019;Scholkmann et al., 2013; Tachtsidis & Scholkmann,  (Brett et al., 2002) and the xjView toolbox (https://www.alivelearn.net/xjview). The image was then projected onto the template using BrainNet Viewer (Xia et al., 2013). (B-C) Event-related averages of activity for the stuttering group within the R-DLPFC ROI for anticipated vs. unanticipated words, for Oxyhemoglobin (HbO) and DeOxyhemoglobin (HbR), respectively. 0 seconds is the beginning of the anticipation phase, and 30 s is shown because the hemodynamic response could feasibly last this long. HbR signals prior to sign reversal such that decreasing values reflect increases in activation strength. Shading represents standard error of the mean.
2016; Zhang et al., 2016). Still, Figure S2 and Table S2 include the uncorrected whole-brain results for both chromophores.
The two within-group analyses that compared activation to rest provided some support that anticipated words are associated with increased activation (p = 0.02), whereas there was no change for unanticipated words (p = 0.82). Thus, it appears that prior to execution, anticipated words recruited the R-DLPFC ROI, whereas unanticipated words did not. Lastly, the correlation between reaction time and anticipation approached significance (r = 0.38, p = 0.08), indicating that for anticipated words, reaction time may increase as activation increases ( Figure 3).
(2) Stuttered vs. Fluent: Activation in the R-DLPFC mask, using the same 5 s time window during the anticipation phase, was also greater for stuttered vs. fluent trials (p = 0.0039). See Figure S3 and Table S3 for uncorrected whole-brain results for HbR as well as HbO. Given that anticipation is strongly associated with stuttering (above), it is possible that the previous results for the anticipated vs. unanticipated contrast may be in part due to atypical planning associated with stuttering that is not necessarily related to anticipation. It was not possible to compare activation associated with anticipated/stuttered, anticipated/fluent, unanticipated/ stuttered, and unanticipated/fluent, due to the unbalanced distribution of stuttering within participants (see Table 3) and the relatively limited number of trials. Importantly, however, 33.4% of unanticipated words were unambiguously stuttered, and 35.4% of anticipated words were unambiguously fluent (35.4%).
(3) Stutterers vs. Controls: Activation in the R-DLPFC mask, during the anticipation phase, was greater for stutterers compared to control speakers (p = 0.0442), irrespective of anticipation and stuttering. See Figure S4 and Table S4 for uncorrected whole-brain results for HbR and HbO. (4) Task vs. Rest (controls only): We did not find evidence that R-DLPFC was significantly activated during the "anticipation" phase in the control group (p = 0.1764). This may indicate that the area in the R-DLPFC mask was not recruited for speech planning, which is in line with previous accounts of speech motor control that do not attribute speech motor function to R-DLPFC.
(5) Interactive vs. Alone: Condition was not a variable of interest in the current study, but it was important to confirm that condition (interactive vs. alone) did not contribute to R-DLPFC overactivation. A significant difference between the interactive and alone conditions was not observed (p = 0.4826).
Significant differences in activation were not observed for anticipated vs. unanticipated words for the two control ROIs: L-DLPFC (p = 0.99); R-preCG (p = 0.95). The lack of significant differences in these areas provides evidence that the differences observed in the R-DLPFC ROI were due to the contrast (anticipated vs. unanticipated) and not to general differences associated with the stuttering brain or systemic artifact. In addition, significant activation differences were not observed between anticipated vs. unanticipated words in R-IFG (p = 0.75) or R-preSMA (p = 0.39).

Functional connectivity
Functional connectivity between the R-DLPFC and R-SMG was assessed using each node as a seed and target, and Holm-Bonferroni correction was applied for the two comparisons. Compared to unanticipated words, anticipated words were associated with lower intrinsic connectivity when R-DLPFC was the seed and R-SMG was the target (t = −2.89, p = 0.01), and also when R-SMG was the seed and R-DLPFC was the target (t = −2.06, p = 0.05). This reduction in functional connectivity between the R-DLPFC and R-SMG for anticipated relative to unanticipated words is taken as evidence of involvement of the FPN in stuttering anticipation.

DISCUSSION
That stutterers anticipate overt stuttering events is well known, but the neural substrates of anticipation had not been studied previously. In this study, we used a novel, clinically inspired approach to identify anticipated and unanticipated words in a relatively large sample of adults who stutter. The words were produced in a delayed-response paradigm while neural signals were recorded with fNIRS. We identified R-DLPFC as a neural substrate of stuttering anticipation. A connectivity analysis was also conducted to explore whether the FPN, specifically the R-DLPFC and R-SMG, is associated with stuttering anticipation. Results are discussed in the context of theoretical accounts of stuttering anticipation, error-likelihood monitoring, and action-stopping. Our findings and potential limitations, as well as possible clinical implications, are also discussed in the following sections.

Right Dorsolateral Prefrontal Cortex Underlies Stuttering Anticipation
The primary hypothesis was confirmed-anticipated words are associated with greater pre-execution activation in the R-DLPFC, compared to unanticipated words. This means that the production of words previously identified by participants as being difficult or likely to be stuttered, up-regulates activation in this area. It was also shown that stutterers exhibit greater activation than non-stutterers during this same time period, irrespective of anticipation and stuttering, and while anticipated words elicited activation in the R-DLPFC ROI in stutterers, unanticipated words of stutterers, and all words produced by controls, did not elicit activation in this area. Further, anticipated words were associated with longer reaction times, and there was some indication that as activation in R-DLPFC increased, so did reaction time (for anticipated words only). This extra time may be due to speakers delaying speech onset until the word can appear fluent to listeners, or "letting the stuttering pass," which could be a function of the R-DLPFC. Our results provide some clarification of Kell et al. (2009), which found that stutterers exhibit greater activation than controls in the R-DLPFC, but that after therapy this is not the case, suggesting that therapy down-regulates an overactive R-DLPFC. Kell et al. (2009) also found that stutterers who reportedly "recovered" from stuttering after early childhood exhibited similar activation compared to controls in the R-DLPFC, suggesting that elevated activation in the R-DLPFC is a maladaptive response to stuttering. Our results suggest that the Kell et al. (2009) result was due to stuttering anticipation. It is important to highlight that we focused on longer-(vs. shorter-) term anticipation in this study, as stimuli were identified between 3 and 10 days prior to the fNIRS experiment. While it is possible that the neural processes for longer-and shorter-term anticipation overlap, many words identified as anticipated before the experiment were not stuttered during the experiment, indicating that they may not have been anticipated either, suggesting that there may be differences in the underlying processes related to longer-vs. shorter-term anticipation. Future studies can attempt to disentangle these types of anticipation to provide clarity on the time scales associated with anticipation.
Stuttered speech was also associated with elevated activation in R-DLPFC, calling into question whether atypical activation in this area was due to anticipation, or aberrant planning associated with stuttering. Given our novel approach, it was not possible to completely differentiate processes associated with anticipation and aberrant planning because we did not have enough power to make comparisons between anticipated-stuttered, anticipated-fluent, unanticipated-stuttered, unanticipated-fluent speech. fNIRS, like fMRI, relies on the slow hemodynamic response, which requires longer trials. With 80 trials already, each lasting at least 20 s, the experiment lasted approximately 30 min, which pushed the comfort threshold for participants. Together with the unpredictability of stuttering, this resulted in too few trials for the aforementioned analysis. However, there is reason to believe that activation differences were in fact driven by anticipation. First, a significant portion of anticipated/unanticipated words were fluent/stuttered, respectively. 33% of anticipated words were unambiguously fluent, whereas 33% of unanticipated words were unambiguously stuttered. Second, it is likely that anticipation was present throughout the 5 s anticipation phase, whereas it is unlikely that speech planning processes would comprise the entire 5 s anticipation phase. Third, we did not find evidence that non-stutterers activate R-DLPFC during the anticipation phase, suggesting that R-DLPFC does not play a primary role in speech planning. This is consistent with work on speech motor control or production models that do not include DLPFC as a part of the speech motor network (e.g., Forseth et al., 2021;Guenther, 2016;Sörös et al., 2006). Future work can tease these processes apart by obtaining enough trials to make the necessary statistical comparisons. It should also be noted, however, that due to the dynamic nature of both anticipation and stuttering, it may not be possible to completely differentiate these processes, especially because speech execution, or the possibility of it, is required to elicit anticipation in speakers.

Error-Likelihood Monitoring
Anticipation is likely driven by error-likelihood monitoring, which refers to the ability to predict errors based on prior experience making those errors. In their original account, Brown and Braver (2005) showed that ACC learns to predict the likelihood of errors, generating a "warning signal" to heighten readiness or initiate cognitive control in response to predicted errors. Arenas (2012Arenas ( , 2017 and Garcia-Barrera and Davidow (2015) extended the error-likelihood account to stuttering, whereby associative learning is the basis for anticipation. The speaker learns that some words/sounds are difficult, and when the speaker next says these words/sounds, they are primed to respond to upcoming stuttering. DLPFC works in concert with ACC to detect and respond to anticipated errors, such that ACC underlies the detection of errors in response to unintended outcomes, and subsequently generates error signals, whereas DLPFC generates representations of these errors including holding task-relevant information in working memory, and initiating subsequent actions (Alexander & Brown, 2015;Botvinick et al., 2001;Donoso et al., 2014;Holroyd & Yeung, 2012;Kolling et al., 2012). Interestingly, while Brown and Braver (2005) primarily focused on ACC, they also found errorlikelihood effects in the R-DLPFC and R-SMG, which were less connected functionally for anticipated vs. unanticipated words in the current study. Thus, anticipation may destabilize the FPN, potentially reducing speakers' control in responding to anticipation. We propose a model that extends the Arenas and Garcia-Barrera and Davidow accounts to the right FPN (see Figure 4). The ACC detects the error (i.e., anticipates) and subsequently generates an error signal that is sent to the R-DLPFC, which coupled with R-SMG, holds this information in working memory and initiates a response. Ultimately, we could not measure ACC due to imaging depth restrictions associated with fNIRS. Future work can test this proposal using fMRI.

Action-stopping
A related interpretation involves the action-stopping network. Action-stopping has been studied primarily in the lab using the stop-signal task, in which participants are required to stop an initiated response when a signal occurs. Aron and colleagues (Aron, 2011;Wessel & Aron, 2017) proposed a prefrontal-basal ganglia-thalamocortical model of action-stopping in which prefrontal areas (R-IFG, R-preSMA) receive sensory information and subsequently produce a Figure 4. A model of stuttering anticipation. Anticipation develops via associative learning. Memory traces of associations form between overtly stuttered utterances and environmental/listener reactions. ACC detects upcoming stuttering via error-likelihood monitoring and signals R-DLPFC, which initiates cognitive control to respond to anticipation, thereby elevating activation in R-DLPFC. R-DLPFC initiates activation within the broader FPN to support the response. Atypical connection between R-DLPFC and areas within the FPN, such as SMG, impedes the speaker's ability to adaptively respond to anticipation. ACC = anterior cingulate cortex; R-DLPFC = right dorsolateral prefrontal cortex; FPN = frontoparietal control network; R-SMG = right supramarginal gyrus. stop command, which is then transmitted to the subthalamic nucleus. The subthalamic nucleus then relays this information to the globus pallidus interna of the basal ganglia, which then inhibits the excitatory drive to motor cortex, i.e., to stop the action (Aron, 2011;Wessel & Aron, 2017). It has been proposed that this mechanism underlies stuttering by inducing a global inhibition response (over-suppression) that impedes the execution of successive motor programs (Arenas, 2017;Aron et al., 2014;Neef et al., 2018). Our results did not support this hypothesis-we did not find evidence of heightened activity in R-IFG or R-preSMA. This may be because the global inhibition account attempts to explain stuttering events at a speech motor control level, whereas anticipation occurs at a cognitive control level.
Disentangling reactive and proactive control, two forms of action-stopping, may provide clarity regarding the distinction between speech motor and cognitive control processes associated with stuttering. Reactive control is stimulus-driven and habitual or automatic, whereas proactive control is prospective and goal-directed (Hannah & Aron, 2021). Reactive control is supported by the hyperdirect cortico-subthalamic-pallidal-thalamo-cortical pathway, and includes superficial cortical structures, the R-IFG, and R-preSMA (Jahanshahi et al., 2015). This pathway is likely associated with the global inhibition that has been proposed to prevent the execution of successive motor programs (e.g., Arenas, 2017;Aron et al., 2014;Neef et al., 2018), i.e., to cause or trigger stuttering events. Proactive control, on the other hand, is regulated by the indirect fronto-striato-pallido-thalamo-cortical pathway, which includes the R-DLPFC, R-SMG, caudate, and thalamus (A. Chang et al., 2017;Chikazoe et al., 2009;Jahanshahi et al., 2015;Jahfari et al., 2010;Marek & Dosenbach, 2018;Smittenaar et al., 2013). Neef et al. (2018) proposed that proactive control may underlie responses to stuttering anticipation, which is supported by the current findings of overactivation in the R-DLPFC and reduced connectivity between R-DLPFC and R-SMG. Thus, while the indirect and hyperdirect pathways share neural circuitry, the hyperdirect pathway may underlie global suppression to prevent the succession of speech gestures (i.e., the stuttering event), whereas the indirect pathway seems to be related to how the speaker responds to knowing that the stuttering event is going to happen. The temporal dynamics associated with fNIRS as well as the slow hemodynamic response makes fNIRS a sub-optimal tool to test this hypothesis. Future studies could use fMRI and MEG to tease apart the specific spatial and temporal contributions (respectively) of both pathways, as well as interactions with the error monitoring system (including the ACC), to the manifestation of overt stuttering events and the speaker's anticipation of them.

Clinical Implications
The current work may inform neuromodulation techniques that have recently been applied to stuttering in clinical trials. Chesters et al. (2018) tested the impact of anodal transcranial direct current stimulation (tDCS) of left inferior frontal cortex on overt severity in 30 adult stutterers. They found that the treatment significantly reduced overt stuttering severity at a 1-week follow-up assessment, and at 6 weeks for reading but not conversational speech (Chesters et al., 2018). Garnett et al. (2019) tested the impact of anodal tDCS on overt severity in 14 adult stutterers. They did not find significant effects on overt stuttering severity, though they found that the atypically strong association between overt severity and right thalamocortical activity was attenuated, especially in severe stutterers. It may be that the modest effects reported to date are due to the focus on single areas, specifically in the speech network. Stuttering is not simply a speech motor control problem. Our data show that the cognitive and sensorimotor processes that underlie anticipation and subsequent overt stuttering elicit elevated activation in R-DLPFC, as well as potentially other areas within the FPN. Down-regulating R-DLPFC using tDCS and concurrently up-regulating, for example, left speech-language areas (premotor cortex/inferior frontal gyrus), which are typically under-active in speakers who stutter, could be the basis for clinical trials of the effects of tDCS in stuttering therapy.
Behaviorally, the anticipation of stuttering allows stutterers to disguise stuttering, such that there will be a discrepancy between what listeners hear or see and what speakers experience. While anticipating stuttering, the speaker may experience anxiety, fear, shame, or other cognitive responses, but the listener may not be privy to this information, creating misunderstanding between the speaker and the listener that could lead to negative listener perceptions of stuttering. For example, a listener might judge a speaker for looking "nervous" or for not being intelligent because they do not respond in a timely manner (as perceived by the listener). Understanding how the brain functions during anticipation could mitigate this discrepancy through increased understanding of the phenomenon and public awareness. This could also improve how SLPs work with people who stutter, particularly with respect to helping clients develop adaptive responses to anticipation. Developing such adaptive responses during therapy can be challenging, primarily because of the unobservable or "hidden" nature of anticipation. A brain-based understanding of anticipation may provide an entry point to begin discussing anticipation with their clients.

Limitations
One limitation of this study is that the hypothesized effect in R-DLPFC-activation greater for anticipated vs. unanticipated words-was only evident based on the HbR signal, not HbO. For other motor tasks (e.g., finger tapping), effects are expected to be present for both chromophores. However, it is important to highlight that there is evidence that HbR may be a more reliable signal for speech tasks. For example, HbR is more strongly correlated to the blood oxygen level dependent signal compared to HbO (Boas et al., 2004(Boas et al., , 2014Ferrari & Quaresima, 2012), has greater spatial specificity than HbO Zhang et al., 2016), and has been validated for speech-language tasks . Further, HbR is less sensitive to systemic effects (e.g., heart rate, breathing, blood flow) than HbO (Franceschini et al., 2003;Kirilina et al., 2012;Santosa et al., 2019;Scholkmann et al., 2013;Tachtsidis & Scholkmann, 2016;Zhang et al., 2016), which may be relevant to stuttering because stuttering speakers may experience elevated anxiety/heart rate prior to speaking tasks. Neither anxiety nor autonomic data were collected in the current study, but either may have contributed to activation in R-DLPFC. Future studies could examine the impact of systemic effects on HbO and HbR, particularly for populations that may exhibit greater systemic effects such as anxiety. This would allow researchers both to determine whether systemic effects differentially impact HbO/HbR, and also to differentiate between anticipation and anxiety, to the extent that the latter is even possible.
Another limitation is that the design did not allow us to differentiate between the awareness of upcoming stuttering and the speaker's response to the knowledge that upcoming speech will be stuttered (should they initiate their speech plan as intended). Thus, anticipation and responding to anticipation are conflated. Future work can address this limitation, for example, by testing the model of anticipation proposed herein, which postulates that the initial awareness of upcoming speech breakdown may occur in ACC, whereas the response to anticipation may be supported by R-DLPFC and the broader FPN. A third limitation is that anticipation was not measured during the actual experiment (e.g., by using a button press), and therefore there was no indication of the extent to which participants anticipated, if they anticipated at all. There were several reasons for this. First, we did not want to increase the length of the experiment, which was already nearing the participant comfort threshold. Identifying whether anticipation occurred could be difficult for some, if not all, participants, which would require additional time. Second, the act of identifying anticipation may alter neural responses. If identification comes before the trial, the act of identifying itself could alter the neural response. If identification occurs after the trial, responses may be biased based on whether the trial was overtly stuttered. One compromise could be to ask participants to identify the extent of anticipation for all words that were previously identified as anticipated just prior to the experiment. Participants could be presented with stimuli as they would during the actual experiment, which may help with identifying whether a word would be anticipated.

Conclusion
This study determined that the R-DLPFC is a likely neural substrate of stuttering anticipation, and also that anticipation may be associated with reduced connectivity within the right hemisphere FPN. The results support accounts of error-likelihood monitoring and action-stopping and their association with stuttering events. Future investigations will benefit from adapting the current paradigm for use with fMRI and MEG to determine the relationships between error-likelihood monitoring, action-stopping, and stuttering, and whether functional connectivity within the FPN and related networks (salience, cingulo-opercular) is weaker for anticipated words. It will also be critical to study children who stutter so that the developmental bases of anticipation can be determined. Finally, results from this study may inform clinical trials that test the efficacy of neuromodulation in stuttering therapy, particularly by focusing not only on speech motor control networks, but also cognitive control and related networks.