INTRODUCTION
After mild traumatic brain injury (mTBI), neuropsychological functioning typically recovers to pre-injury levels (Belanger & Vanderploeg, Reference Belanger and Vanderploeg2005; Frencham, Fox, & Maybery, Reference Frencham, Fox and Maybery2005; McCrea, Reference McCrea2008). One notable exception, however, is the observation of residual behavioral (Bonnelle et al., Reference Bonnelle, Ham, Leech, Kinnunen, Mehta, Greenwood and Sharp2012; Pontifex, O’Connor, Broglio, & Hillman, Reference Pontifex, O’Connor, Broglio and Hillman2009; Seignourel et al., Reference Seignourel, Robins, Larson, Demery, Cole and Perlstein2005) and neural (Mayer et al., Reference Mayer, Yang, Yeo, Pena, Ling, Mannell, Stippler and Mojtahed2012; Pontifex et al., Reference Pontifex, O’Connor, Broglio and Hillman2009) alterations in cognitive control, the processes that allow for the flexible modulation of information processing in the service of goal-directed behavior (Botvinick, Braver, Barch, Carter, & Cohen, Reference Botvinick, Braver, Barch, Carter and Cohen2001). Cognitive control involves detecting salient events and errors as well as adjustment of attention in response to such events. Notably, neural alterations in cognitive control have been reported even in the absence of observable behavioral impairment in mTBI (Broglio, Pontifex, O’Connor, & Hillman, Reference Broglio, Pontifex, O’Connor and Hillman2009; Mayer et al., Reference Mayer, Yang, Yeo, Pena, Ling, Mannell, Stippler and Mojtahed2012, Reference Mayer, Hanlon, Dodd, Ling, Klimaj and Meier2015), suggesting that even when cognitive control is intact, its neural implementation may be altered by mTBI.
Although early work focused on residual cognitive control deficits in civilian mTBI, recent studies of blast-related mTBI have also shown neural alterations during the performance of cognitive control tasks, albeit in the absence of behavioral impairment. Scheibel and colleagues (2012) compared functional activation in a group of individuals with chronic blast-related mTBI and individuals without blast exposure or TBI in the context of a stimulus-response compatibility task. They found that mTBI was associated with increased activation in anterior regions such as the anterior cingulate cortex, insula, and medial frontal cortex as well as in posterior regions involved in visual processing. Fischer and colleagues (2014) administered a response inhibition task, and found that mTBI was associated with reduced activation during successful inhibition in regions including the medial frontal gyrus, middle frontal gyrus, middle temporal gyrus, and precuneus.
The apparent inconsistency in these findings (i.e., increased vs. decreased activation in mTBI) can be understood with reference to the involvement of distinct functional brain networks in cognitive control. Cognitive control strongly relies on several prefrontal areas, including the dorsal anterior cingulate cortex and dorsolateraleral prefrontal cortex, as well as the insula and posterior parietal cortex (Bunge, Hazeltine, Scanlon, Rosen, & Gabrieli, Reference Bunge, Hazeltine, Scanlon, Rosen and Gabrieli2002; Kerns et al., Reference Kerns, Cohen, MacDonald, Cho, Stenger and Carter2004; MacDonald, Cohen, Stenger, & Carter, Reference MacDonald, Cohen, Stenger and Carter2000). The insula and dorsal anterior cingulate cortex are two core hubs of the salience network, whereas the dorsolateral prefrontal cortex and posterior parietal cortex are key hubs of the central executive network (Seeley et al., Reference Seeley, Menon, Schatzberg, Keller, Glover, Kenna and Greicius2007). These networks typically co-activate during cognitive control tasks (Dosenbach, Fair, Cohen, Schlaggar, & Petersen, Reference Dosenbach, Fair, Cohen, Schlaggar and Petersen2008; Nee, Wager, & Jonides, Reference Nee, Wager and Jonides2007). Successful cognitive control also relies on deactivation of the default mode network (Bonnelle et al., Reference Bonnelle, Ham, Leech, Kinnunen, Mehta, Greenwood and Sharp2012; Kelly, Uddin, Biswal, Castellanos, & Milham, Reference Kelly, Uddin, Biswal, Castellanos and Milham2008; Singh & Fawcett, Reference Singh and Fawcett2008), a network that includes the posterior cingulate cortex/precuneus, medial frontal cortex, and temporal cortex and is typically deactivated during tasks. Thus, the results of Scheibel et al. (Reference Scheibel, Newsome, Troyanskaya, Lin, Steinberg, Radaideh and Levin2012), who reported increased activation in salience network regions in mTBI, and Fischer et al. (Reference Fischer, Parsons, Durgerian, Reece, Mourany, Lowe and Rao2014), who reported decreased activation in default mode network regions in mTBI, are consistent with the notion that cognitive control is associated with both heightened salience network activity and reduced default mode network activityFootnote 1 . In both studies, the pattern observed in mTBI represents an amplification of the pattern seen in control participants.
Moreover, cognitive control depends not only on these functional networks independently, but also on their interaction. Sridharan, Levitin, and Menon (Reference Sridharan, Levitin and Menon2008) demonstrated that the salience network switches between internally- and externally-directed cognitive processing by initiating control signals that upregulate the central executive network and deactivate the default mode network. Furthermore, disruption of the structural integrity of the salience network predicts reduced default mode network deactivation during a stop-signal task in moderate–severe TBI (Bonnelle et al., Reference Bonnelle, Ham, Leech, Kinnunen, Mehta, Greenwood and Sharp2012). However, it is unknown how mTBI impacts the interaction of these networks during cognitive control. In the present study, we extend previous work on cognitive control in blast-related mTBI by assessing the functioning of these networks as well as their interactions. To do this, we administered a functional magnetic resonance imaging (fMRI) flanker paradigm to Operation Enduring Freedom/Operation Iraqi Freedom (OEF/OIF) blast-exposed veterans with and without mTBI and used psychophysiological interaction (PPI) analysis (Friston et al., Reference Friston, Buechel, Fink, Morris, Rolls and Dolan1997; O’Reilly, Woolrich, Behrens, Smith, & Johansen-Berg, Reference O’Reilly, Woolrich, Behrens, Smith and Johansen-Berg2012) to assess task-based functional connectivity. Based on previous studies in blast-related mTBI, we hypothesized that mTBI would be associated with greater activation in salience network regions and/or greater deactivation in default mode network regions during increased demands on cognitive control. Furthermore, we hypothesized that the interaction between these networks would be altered in mTBI.
The flanker task requires participants to respond to a center arrow that is flanked by arrows in the same direction (congruent trials) or in the opposite direction (incongruent trials). Incongruent trials pose high response conflict because of the discrepancy in the direction of the flanking arrows compared to the center arrow, and require inhibition of task irrelevant information (i.e., flanking arrows) to correctly respond to task relevant stimuli (i.e., center arrow). The difference in response latency to incongruent versus congruent trials thus constitutes a measure of cognitive control. Additionally, processing during error trials can be examined as a measure of cognitive control, as error trials yield a response conflict resulting from competition between a correct response and a strong response tendency for an incorrect response. Moreover, several studies suggest that error processing concerns not only the identification of the error but also the correction of differences between the intended task goal and executed response (Menon, Adleman, White, Glover, & Reiss, Reference Menon, Adleman, White, Glover and Reiss2001; Taylor, Stern, & Gehring, Reference Taylor, Stern and Gehring2007; Ullsperger & von Cramon, Reference Ullsperger and von Cramon2001). Thus, the difference in response latency on trials immediately following an error (post-error) and trials immediately following a correct response (post-correct) constitutes an additional behavioral measure of cognitive control (i.e., post-error slowing). Using these measures, we examined behavioral and neural alterations in cognitive control associated with blast-related mTBI.
METHODS
Participants
Of 69 OEF/OIF veterans initially contacted for this study, 38 agreed to participate. Three were not enrolled because of exclusionary criteria (see below), leaving 35 participants who completed the protocol. The study sample consisted of 18 veterans with blast-related mTBI, as defined by the American Congress Rehabilitation Medicine (1993) criteria, and 17 blast-exposed veterans who reported no TBI from any mechanism of injury during deployment (control). TBI assessment was based on an extensive clinical interview described in detail in Verfaellie, Lafleche, Spiro, Tun, and Bousquet (Reference Verfaellie, Lafleche, Spiro, Tun and Bousquet2013). In brief, participants were queried about their blast exposure(s) to determine the index event, which they were then asked to describe in detail. Two investigators evaluated the interviews and sought consensus as to whether mTBI criteria had been met and whether any reported disorientation was the result of mTBI. The mTBI group consisted of seven individuals with loss of consciousness (LOC) and 11 without LOC. Study procedures were approved by the VA Boston Institutional Review Board and all participants provided written informed consent consistent with the Declaration of Helsinki.
Participants were excluded from the study if they reported a history of pre-deployment TBI with LOC or with symptoms persisting longer than three months post-injury, demonstrated questionable effort with raw scores below 45 on the retention trial of the Test of Memory Malingering (TOMM; Tombaugh & Tombaugh, Reference Tombaugh and Tombaugh1996), had structural brain abnormalities (e.g., hemorrhages, hematomas) on T2-FLAIR, susceptibility weighted imaging (SWI), or T1-weighted sequences as determined by a board-certified neuroradiologist, showed evidence of excessive alcohol use as reflected by scores above 20 on the Alcohol Use Disorder Identification Test (AUDIT), or reported a diagnosis of attention deficit hyperactivity disorder (ADHD) or medication use consistent with its treatment.
Two participants (one mTBI, one control) were unable to stay awake during the task, thus yielding a final sample of 33 participants. A summary of demographic characteristics can be found in Table 1.
Summary of demographic and clinical characteristics of participants
Note: For the control group, the scan interval indicates the time from blast exposure to scan, whereas for the mTBI group, the scan interval indicates the time from blast-related mTBI to scan. LOC=loss of consciousness; mTBI=mild traumatic brain injury; PCL-M=PTSD Checklist-Military version.
Clinical Assessment
The posttraumatic stress disorder (PTSD) Checklist – Military Version (PCL-M; Weathers, Huska, & Keane, Reference Weathers, Huska and Keane1991) was used to measure PTSD symptom severity within the last month preceding testing. The PCL has good convergent validity with the Clinician-Administered PTSD Scale (Wilkins, Lang, & Norman, Reference Wilkins, Lang and Norman2011), which is the gold standard for PTSD assessment (Blake et al., Reference Blake, Weathers, Nagy, Kaloupek, Gusman, Charney and Keane1995).
Flanker Task
The flanker task was administered in the scanner as an event-related fMRI paradigm. Visual stimuli were presented with E-Prime 2.10 Software (Psychology Software Tools, Pittsburg, PA) and were projected to a screen at the back of the scanner, which participants viewed with a mirror attached to the head coil. On each trial, participants viewed a string of arrows that was presented for 200 ms and was immediately followed by a crosshair that was randomly jittered between 1800 ms and 7800 ms (mean=4300 ms; see Figure 1). Participants’ task was to respond to the direction of the center arrow, which was surrounded by flanking arrows on either side. On half the trials, the flanking arrows pointed in the same direction as the center arrow (congruent condition), whereas on the other half of trials, the flanking arrows pointed in the opposite direction from the center arrow (incongruent condition). Direction of the center arrow was counterbalanced across trials. The order of stimulus presentation was pseudo-randomized to ensure that no more than three incongruent trials or no more than three trials with the center arrow pointing in the same direction appeared in a row. Responses were made with the index and middle fingers of the right hand. Responses were collected up to 2000 ms after stimulus onsetFootnote 2 . There were four runs with 80 trials per run, equaling a total of 320 trials. The order of runs was counterbalanced across participants. Instructions and practice were given outside the scanner a half hour before scanning took place. During fMRI data acquisition, response accuracy, onset time, and reaction time were recorded for each stimulus and only correct trials were included in analyses of response latency.
Flanker task. Stimuli were presented for 200 ms followed immediately by a cross hair with presentation duration randomly jittered between 1800 and 7800 ms.
Neuroimaging Acquisition
Data were acquired with a 32-channel head coil on a 3 Tesla Siemens Trio whole-body MRI scanner located at the VA Boston Healthcare System, Jamaica Plain campus. An auto align scout scan was acquired first. One T1-weighted three-dimensional magnetization-prepared rapid gradient-echo imaging (MP RAGE) scan was collected for each participant [field of view (FOV)=256; matrix=256×256×176 slices; 1×1×1 mm voxels; repetition time (TR)=2530 ms; echo time (TE)=3.32 ms; flip angle=7°]. A T2-FLAIR image was also collected (FOV=256; Matrix=512×512×160 slices; 0.49×0.49×1 mm voxels; TR=6000 ms; TE=388 ms; flip angle=120º) for each participant. Four blood oxygen level dependent (BOLD) T2*-weighted echo-planar imaging runs were acquired parallel to the anterior commissure-posterior commissure plane (FOV=192; TR=2000 ms; TE=30 ms; voxel size=2.67×2.67×3.75 mm; slice order=interleaved; flip angle=90º; matrix=722; volumes=185). The first five volumes of each run collected before stimulus presentation began were discarded to allow for signal magnetization equilibrium.
Behavioral Analysis
Statistical analyses were performed with SPSS, version 19 (IBM Corp., Armonk, NY). Demographic data were analyzed with independent samples t tests for linear variables and χ2 tests for categorical variables. Congruency effects were analyzed using 2 (group: mTBI, control)×2 (condition: incongruent, congruent) repeated measures analyses of covariance (ANCOVA) with reaction time and accuracy as the dependent measures, respectively. Effects of error processing were analyzed using an ANCOVA with post-error slowing scores as the dependent measure and group as the independent measure. PCL-M scores, number of pre-deployment TBIs, and age were entered as covariates in all analyses. However, because these variables did not contribute significant variance, they are not reported. Assumptions for ANCOVA were checked including normality, linear relationships between covariates for each group, outliers, homogeneity for regression slopes, homogeneity of covariance (Box’s M Test), and homogeneity of variance (Levene’s Test). No assumptions were violated, justifying the use of ANCOVA.
Neuroimaging Analysis
All preprocessing procedures and analyses were carried out using The Oxford Centre for FMRIB FSL software package (version 4.15; http://www.fmrib.ox.ac.uk/fsl). FMRI data were processed using fMRI Expert Analysis Tool (FEAT; Version 5.98). Data were preprocessed with the following pre-statistics: motion correction using MCFLIRT (Jenkinson, Bannister, Brady, & Smith, Reference Jenkinson, Bannister, Brady and Smith2002), slice-timing correction using Fourier-space time-series phase-shifting, non-brain removal using BET (Smith, Reference Smith2002), spatial smoothing using a Gaussian kernel of full-width/half-max 5 mm, and grand-mean intensity normalization of the entire 4D dataset by a single multiplicative factor. To remove head motion artifact, we used a data-driven independent component analysis (ICA) method to identify and remove motion-related components from the data (ICA-based strategy for Automatic Removal of Motion Artifact [ICA-AROMA]; Pruim, Mennes, van Rooij, et al., Reference Pruim, Mennes, van Rooij, Llera, Buitelaar and Beckmann2015). After removing these motion components, we removed signal from white matter and cerebral spinal fluid using nuisance regression to further minimize noise-related artifact in the functional data (Pruim, Mennes, Buitelaar, & Beckmann, Reference Pruim, Mennes, Buitelaar and Beckmann2015). Next, we applied linear detrending and a highpass temporal filter (σ=45.0 s). Registration to high-resolution structural and standard space images was carried out using FMRIB’s Linear Image Registration Tool (FLIRT) and further refined using FMRIB’s Nonlinear Image Registration Tool (FNIRT).
Functional analysis of incongruent trial processing was performed using an incongruent>congruent contrast for correct trials. The processing of error trials was analyzed using an incorrect>correct contrast for congruent and incongruent trials combined. Higher-level analyses were carried out using FMRIB’s Local Analysis of Mixed Effects (FLAME) stage 1 (Beckmann, Jenkinson, & Smith, Reference Beckmann, Jenkinson and Smith2003; Woolrich, Reference Woolrich2008; Woolrich, Behrens, Beckmann, Jenkinson, & Smith, Reference Woolrich, Behrens, Beckmann, Jenkinson and Smith2004). To examine each contrast, runs were combined for each participant. To determine activation differences across groups, group level activation maps were generated for each contrast using FLAME stage 1. Age, number of pre-deployment TBIs, and PTSD symptom severity were entered into the model as regressors. To examine group differences at the significance level of p=.005, Z statistic images were thresholded using clusters determined by Z>2.6 and a corrected cluster significance threshold of p=.05.
To examine the implications of observed group activation differences within the context of larger networks, we performed a PPI analysis (Friston et al., Reference Friston, Buechel, Fink, Morris, Rolls and Dolan1997; O’Reilly et al., Reference O’Reilly, Woolrich, Behrens, Smith and Johansen-Berg2012). PPI analysis involves a psychological regressor (i.e., incongruent vs. congruent or incorrect vs. correct responses), a physiological regressor, which is the time course of the seed region of interest (ROI), and the PPI term, which is the interaction between the psychological and physiological regressors. Group level PPI maps were generated to determine group differences in brain regions that were modulated by the interaction of the task and activation of the seed ROI. We extracted the time course of the seed region, defined as the significant region from the functional group analysis, from all participants. Higher-level analyses for the PPI regressor were carried out using FLAME stage 1. Runs were combined for all participants. To determine group differences, age, number of pre-deployment TBIs, and PTSD symptom severity were entered into the model as regressors and group PPI maps were generated. Z statistic images were thresholded using clusters determined by Z>2.6 and a corrected cluster significance threshold of p=.05. The maps of the PPI analysis represent the effects of the interaction that are over and above the main effects of the BOLD response to the task contrast and correlations with the seed region.
Finally, to identify brain regions that were associated with behavior, we performed a mixed effects (FLAME stage 1) group-level analysis using behavior as a regressor. Because significant group differences in neural activation were found only when examining the processing of errors, we performed this analysis for error-related performance only. We calculated a post-error slowing score, reflecting the difference in reaction time on trials immediately following an error (post-error) and trials immediately following a correct response (post-correct). Next, we identified brain regions that were positively and/or negatively associated with post-error slowing within the incorrect>correct contrast, by using the post-error slowing score as a regressor in the analysis. Additionally, to determine whether these brain-behavior associations differed across groups, we conducted a continuous covariate interaction analysis. The post-error slowing score was entered as a separate regressor for each group. Age, number of pre-deployment TBIs, and PTSD symptom severity were entered as regressors in all analyses. Z statistic images were thresholded using clusters determined by Z>2.6 and a corrected cluster significance threshold of p=.05.
RESULTS
Behavioral Performance
We examined congruency effects in both accuracy and reaction time data (see Table 2). A 2×2 repeated measures ANCOVA of accuracy data with group as the between subjects factor and congruency as the within subjects factor revealed that accuracy did not significantly differ as a function of group (F(1,28)<1; p>.4) or congruency (F(1,28)<1; p>.9). Moreover, the group by congruency interaction was also not significant (F(1,28)<1; p>.8). As a follow-up analysis of accuracy, we performed a signal detection analysis of discriminability and bias in which correct responses on congruent trials were considered hits and incorrect responses on incongruent trials were considered false alarms. Results revealed that neither discriminability nor response bias significantly differed as a function of group (d’: F(1,28)=1.4; p>.2; beta: F(1,28)<1; p>.4). Analysis of latency data revealed that performance again did not significantly differ as a function of group (F(1,28)<1; p>.5). Performance was numerically slower in the incongruent condition than congruent condition, but the difference was not significant (F(1,28)=2.1; p>.1). The group by congruency interaction was also not significant (F(1,28)<1; p>.6).
Performance as a function of congruency and post-error slowing
Note: Accuracy is reported as a percent. Reaction times are in milliseconds. mTBI=mild traumatic brain injury; RT=reaction time.
We next examined the behavioral measure of cognitive control as indexed by post-error slowing. An ANCOVA revealed that groups did not significantly differ in the magnitude of post-error slowing (F(1,28)<1; p>.9)
Neuroimaging Results
Processing of incongruent trials
The incongruent>congruent contrast was examined to determine if there were significant group differences in brain activation associated with the processing of incongruent information. There were no significant group differences in any brain region in the incongruent>congruent contrast. When we examined group maps separately, the control group showed increased activation in the right superior parietal lobe; the mTBI group showed increased activation in the left superior parietal lobe, left dorsal anterior cingulate cortex, right supramarginal gyrus, and right lateral occipital cortex (see Supplementary Table 1). Covariate effects are reported in the supplemental materials.
Processing of error trials
The incorrect>correct contrast was examined to determine if there were significant group differences in brain activation associated with the processing of errors. Both groups showed increased activation in the insula and dorsolateral prefrontal cortex (see Supplementary Table 2), regions that are part of the salience network and central executive network, respectively. There were no group differences in activation in these regions. By contrast, compared to controls, individuals with mTBI showed greater deactivation in areas of the default mode network including the left dorsomedial prefrontal cortex (dMPFC) and left posterior cingulate cortex (PCC)/precuneus (see Table 3; Figure 2). Covariate effects are reported in the supplemental materials.
Compared to controls, the mTBI group had significant deactivation in default mode network regions for the incorrect>correct contrast. In particular, the mTBI group had greater deactivation in the left dMPFC and left PCC. The hemodynamic response function for the dMPFC is plotted to the right of the figure. dMPFC = dorsomedial prefrontal cortex; mTBI = mild traumatic brain injury; PCC = posterior cingulate cortex.
Significant group differences in brain regions for the incorrect>correct contrast
Note: *=for these regions, control>mTBI for incorrect>correct represents greater deactivation in mTBI. Cluster size is number of voxels. Only clusters are reported; sub-clusters are not reported.
dMPFC=dorsomedial prefrontal cortex; mTBI=mild traumatic brain injury; PCC=posterior cingulate cortex.
We next examined whether there were regions where functional connectivity with the default mode network for incorrect (vs. correct) trials differed across groups. To do this, we extracted PCC and dMPFC ROIs based on significant group differences in the incorrect>correct analysis and used the PCC [peak Montreal Neurological Institute (MNI) coordinates=–4 –30 36] and the dMPFC (peak MNI coordinates=–4 60 28) ROIs as seeds in two separate PPI analyses. Using the PCC as a seed, there were no significant group differences in the functional connectivity between the PCC and any brain region for incorrect (vs. correct) trials. Using the dMPFC as a seed, we found that functional coupling with the left dorsal anterior cingulate cortex for incorrect versus correct trials was greater in the mTBI group than the control group. Upon further inspection, results showed that there was increased negative coupling between these regions in the mTBI group (see Table 4; Figure 3). There were no regions where the control group showed greater differential functional connectivity for incorrect versus correct trials than the mTBI group.
The mTBI group had enhanced functional connectivity between the dMPFC of the default mode network and a region of the salience network for the incorrect>correct contrast. Specifically, negative functional connectivity between the left dMPFC and left dorsal anterior cingulate cortex for incorrect versus correct trials was greater for the mTBI group than the control group. The hemodynamic response function is plotted to the right of the figure to show coupling. dACC = dorsal anterior cingulate cortex; dMPFC = dorsomedial prefrontal cortex; mTBI = mild traumatic brain injury.
Significant group differences in functional connectivity with dMPFC seed region in the incorrect>correct contrast as determined by PPI analysis.
Note: Cluster size is number of voxels. Only clusters are reported; sub-clusters are not reported.
dACC=dorsal anterior cingulate cortex; dMPFC=dorsomedial prefrontal cortex; mTBI=mild traumatic brain injury.
In a follow-up analysis, we examined the connectivity with the dMPFC for incorrect (vs. correct) trials within each group separately. Results revealed that within the mTBI group, the dMPFC was functionally coupled with the right insula and left postcentral gyrus, extending into the dorsal anterior cingulate cortex, as well as the right dorsolateral prefrontal cortex and right superior parietal lobe, regions of the salience and central executive networks, for incorrect (vs. correct) trials (see Table 5). This coupling was negative such that these regions are active when the dMPFC is deactivated. In the control group, there were no regions that showed greater coupling with the dMPFC for incorrect (vs. correct) trials.
Significant brain regions for PPI analysis of incorrect>correct contrast in each group
Note: Cluster size is number of voxels. Only clusters are reported; sub-clusters are not reported.
dACC=dorsal anterior cingulate cortex; DLPFC=dorsolateral prefrontal cortex; mTBI=mild traumatic brain injury.
Finally, to examine the functional significance of these neural alterations in cognitive control, we performed a whole-brain imaging analysis to determine whether brain activation associated with error processing was associated with a behavioral measure of cognitive control, as indexed by post-error slowing. The analysis revealed no significant associations between post-error slowing and activation in any brain region. However, group moderated the association between post-error slowing and activation in the right dorsolateral prefrontal cortex (ZMax=3.69; peak MNI coordinates: 46 38 26). Greater post-error slowing was associated with greater recruitment of the right dorsolateral prefrontal cortex in individuals with mTBI, but not in controls (Figure 4).
Increased post-error slowing was significantly associated with greater activity in the right dorsolateral prefrontal cortex, a central executive network region, in mTBI. There was no significant association in the control group. Activation is represented as contrast of parameter estimate values. DLPFC = dorsolateral prefrontal cortex; mTBI = mild traumatic brain injury; R = right.
DISCUSSION
We examined behavioral and neural indices of cognitive control in OEF/OIF veterans with blast-related mTBI in the context of a flanker task. Behavioral performance did not differ in individuals with and without mTBI, but the neural signature of cognitive control was amplified in the mTBI group. That is, with increased demands on cognitive control processes, the mTBI group showed greater deactivation of default mode network regions than the control group. Furthermore, there was enhanced negative connectivity between the dMPFC, a region within the default mode network, and regions of the salience network and central executive network. Taken together, these findings suggest that mTBI did not affect the ability to engage in cognitive control, but altered how such control was neurally implemented.
Although we evaluated activation associated with the processing of both incongruent information and errors as neural measures of cognitive control, group differences emerged for error processing only. The question arises whether the failure to observe group differences in the incongruent–congruent comparison may be due to the fact that this comparison does not provide a pure index of cognitive control. That is, it could be argued that this comparison reflects a combination of priming effects associated with the presence of congruent flankers and demands on cognitive control associated with the presence of incongruent flankers. However, studies that have included a baseline condition to disentangle these two effects have shown that reaction times for congruent and baseline trials do not differ, suggesting that priming effects are negligible and that differences between incongruent and congruent trials are largely due to demands on inhibitory cognitive control processes associated with incongruent trials (Bunge et al., Reference Bunge, Hazeltine, Scanlon, Rosen and Gabrieli2002; Hazeltine, Bunge, Scanlon, & Gabrieli, Reference Hazeltine, Bunge, Scanlon and Gabrieli2003). Such findings argue that the difference between incongruent and congruent trials is an appropriate measure of cognitive control. Thus, the fact that we did not observe group differences in the incongruent-congruent comparison is unlikely to reflect a measurement problem. Rather, the greater sensitivity of error processing in the current study may be due to the fact that the task was relatively easy, and cognitive control was more strongly taxed during the processing of errors than during the processing of incongruent flankers.
In the context of incorrect relative to correct responses and in comparison to controls, mTBI was associated with greater deactivation in the PCC and dMPFC, two regions within the default mode network. The default mode network is a well-established resting state network that is most active at rest and is involved in autobiographical memory retrieval, mind-wandering, and other self-generated thought (Gusnard, Akbudak, Shulman, & Raichle, Reference Gusnard, Akbudak, Shulman and Raichle2001; McKiernan, Kaufman, Kucera-Thompson, & Binder, Reference McKiernan, Kaufman, Kucera-Thompson and Binder2003; Raichle et al., Reference Raichle, MacLeod, Snyder, Powers, Gusnard and Shulman2001). The default mode network is deactivated during tasks that require externally oriented attention (Fox et al., Reference Fox, Snyder, Vincent, Corbetta, Van Essen and Raichle2005; Fransson, Reference Fransson2006; Uddin, Kelly, Biswal, Castellanos, & Milham, Reference Uddin, Kelly, Biswal, Castellanos and Milham2009). Previous work by Fischer et al. (Reference Fischer, Parsons, Durgerian, Reece, Mourany, Lowe and Rao2014) has shown that default mode network deactivation is exaggerated in blast-related mTBI, a finding that is replicated in the current study.
Using PPI analysis to examine functional connectivity, our findings go beyond Fischer et al. (Reference Fischer, Parsons, Durgerian, Reece, Mourany, Lowe and Rao2014) by demonstrating that during the processing of errors, individuals with mTBI, in comparison to controls, show enhanced negative coupling between the dMPFC and the dorsal anterior cingulate cortex, a region of the salience network. When we examined the groups separately, we found that the mTBI group in particular displayed significantly enhanced negative coupling between the dMPFC and regions of both the salience and central executive networks associated with error processing. Previous studies suggest that this coupling between the default mode, salience, and central executive networks facilitates successful task performance (Fransson, Reference Fransson2006; McKiernan et al., Reference McKiernan, Kaufman, Kucera-Thompson and Binder2003). The salience network has been hypothesized to be a “switching” network in cognitive control and may be especially important in switching between the default mode and central executive networks to accomplish task goals (Menon & Uddin, Reference Menon and Uddin2010). Accordingly, Sridharan and colleagues (2008) showed in healthy individuals that with increased demands on cognitive control, the salience network is engaged to suppress default mode network regions and to amplify the response of central executive network regions. In the present study, these network dynamics were up-regulated in the mTBI group. Given that this occurred in the context of intact behavioral performance in the mTBI group, we postulate that the enhanced negative coupling between the default mode network and salience and central executive networks in mTBI serves as a compensatory mechanism for successful task performance. Consistent with this notion, we found an association between post-error slowing and dorsolateral prefrontal cortex recruitment in mTBI, suggesting that recruitment of the central executive network may facilitate the implementation of cognitive control and contribute to adjustments in behavioral performance after an error.
By focusing not only on activation in distinct brain regions, but also on the interaction of the networks they are a part of, our study sheds light on the findings of Scheibel et al. (Reference Scheibel, Newsome, Troyanskaya, Lin, Steinberg, Radaideh and Levin2012) and Fischer et al. (Reference Fischer, Parsons, Durgerian, Reece, Mourany, Lowe and Rao2014) and suggests that the neural alterations reported in those studies can be understood with reference to larger network dynamics. These findings point to the importance of examining networks and their interactions in cognitive control and emphasize that cognitive control is not a singular process but instead relies on the modulation of multiple processes through network interactions.
Altered communication between the functional networks involved in cognitive control has also been observed in civilian TBI, albeit in the context of more severe injury that led to behavioral impairment (Bonnelle et al., Reference Bonnelle, Ham, Leech, Kinnunen, Mehta, Greenwood and Sharp2012). In a recent study focused on cognitive control in civilian mTBI, enhanced activation in the mTBI group was also interpreted as reflecting a compensatory mechanism, but this activation was observed in inferior parietal cortex, a region not part of the functional networks discussed in the current study (Mayer et al., Reference Mayer, Hanlon, Dodd, Ling, Klimaj and Meier2015). Moreover, this study did not examine functional network interactions. Thus, it is unknown whether our findings of compensatory functional network interactions would generalize to civilian mTBI, especially given that the pathology in blast and non-blast mTBI is somewhat different. For example, white matter alterations associated with blast-related mTBI tend to be spatially variable (Hayes, Miller, Lafleche, Salat, & Verfaellie, Reference Hayes, Miller, Lafleche, Salat and Verfaellie2015; Miller, Hayes, Lafleche, Salat, & Verfaellie, Reference Miller, Hayes, Lafleche, Salat and Verfaellie2016), which is in contrast to the more consistent findings of white matter alterations in specific long fiber pathways in non-blast mTBI (Aoki, Inokuchi, Gunshin, Yahagi, & Suwa, Reference Aoki, Inokuchi, Gunshin, Yahagi and Suwa2012; Hayes, Bigler, & Verfaellie, Reference Hayes, Bigler and Verfaellie2016). Thus, it remains an open question whether the mechanism of injury impacts how cognitive control is implemented in mTBI.
Given that the focus of this study was on changes in cognitive control associated with blast-related mTBI, we included a control group of individuals who had been exposed to blast, but did not suffer TBI. The inclusion of a blast-exposed control group in this study thus helped to isolate the contribution of mTBI. However, recent reports suggest that blast exposure itself is associated with neural changes (Robinson et al., Reference Robinson, Lindemer, Fonda, Milberg, McGlinchey and Salat2015; Taber et al., Reference Taber, Hurley, Haswell, Rowland, Hurt, Lamar and Morey2015), leaving open the possibility that blast-exposure might be associated with altered dynamics in the functional networks mediating cognitive control. Future studies examining the effects of blast exposure on cognitive control are needed to evaluate this possibility.
The results reported in this study should be considered within the context of the limitation that mTBI group assignment was based on self-report. However, mTBI assessment was conducted with an in-depth structured clinical interview, which is currently the gold standard of diagnosis (Corrigan & Bogner, Reference Corrigan and Bogner2007). Another limitation is the small sample size of the control and patient groups of this study. The exclusion of individuals with ADHD as well as those with current alcohol abuse limited potential enrollment in this study. It will be important to replicate these findings in studies with larger samples. A third limitation is the inability to examine LOC-associated effects on cognitive control. As shown in Matthews, Simmons, and Strigo (Reference Matthews, Simmons and Strigo2011), LOC may moderate the neural changes associated with cognitive control in blast-related mTBI, but our sample size was too small to allow for a direct comparison between individuals who suffered mTBI with and without LOC. Future studies will need to examine whether there are differences in cognitive control in individuals with mTBI as a function of presence of LOC.
In summary, we report robust neural differences associated with error processing in individuals with chronic blast-related mTBI compared to blast-exposed controls. In particular, individuals with mTBI exhibited increased deactivation in regions of the default mode network. Furthermore, these regions showed greater negative functional connectivity with regions of the salience and central executive networks during the processing of errors. Importantly, these brain changes in mTBI occurred in the context of intact behavioral performance. Taken together, these results suggest that with increased demands on cognitive control, greater deactivation of regions of the default mode network, and enhanced negative coupling between the default mode network and regions of the salience and central executive networks may act as a compensatory mechanism for successful task performance in mTBI.
ACKNOWLEDGMENTS
The contents do not represent the views of the U.S. Department of Veterans Affairs or the United States Government. The authors thank Michael Esterman, Ph.D., for valuable input on this work. The authors declare no conflicts of interest. This work was supported by the VA Rehabilitation Research & Development Service (M.V., grant number I01RX000216 and J.P.H., grant number I21RX001594); the VA Clinical Science Research and Development Service; and the National Center for PTSD. This work was further supported with resources and the use of facilities at the Neuroimaging Research for Veterans Center, VA Boston Healthcare System.
SUPPLEMENTARY MATERIAL
To view supplementary material for this article, please visit https://doi.org/10.1017/S1355617718000279
