ORIGINAL RESEARCH article
Dorsal anterior cingulate cortex modulates supplementary motor area in coordinated unimanual motor behavior
- 1Center for Complex Systems and Brain Sciences, Florida Atlantic University, Boca Raton, FL, USA
- 2Brain Imaging Research Division, Department of Psychiatry and Behavioral Neuroscience, Wayne State University School of Medicine, Detroit, MI, USA
- 3Department of Psychology, Florida Atlantic University, Boca Raton, FL, USA
Motor control is integral to all types of human behavior, and the dorsal Anterior Cingulate Cortex (dACC) is thought to play an important role in the brain network underlying motor control. Yet the role of the dACC in motor control is under-characterized. Here we aimed to characterize the dACC’s role in adolescent brain network interactions during a simple motor control task involving visually coordinated unimanual finger movements. Network interactions were assessed using both undirected and directed functional connectivity analysis of functional Magnetic Resonance Imaging (fMRI) Blood-Oxygen-Level-Dependent (BOLD) signals, comparing the task with a rest condition. The relation between the dACC and Supplementary Motor Area (SMA) was compared to that between the dACC and Primary Motor Cortex (M1). The directed signal from dACC to SMA was significantly elevated during motor control in the task. By contrast, the directed signal from SMA to dACC, both directed signals between dACC and M1, and the undirected functional connections of dACC with SMA and M1, all did not differ between task and rest. Undirected coupling of dACC with both SMA and dACC, and only the dACC-to-SMA directed signal, were significantly greater for a proactive than a reactive task condition, suggesting that dACC plays a role in motor control by maintaining stimulus timing expectancy. Overall, these results suggest that the dACC selectively modulates the SMA during visually coordinated unimanual behavior in adolescence. The role of the dACC as an important brain area for the mediation of task-related motor control may be in place in adolescence, continuing into adulthood. The task and analytic approach described here should be extended to the study of healthy adults to examine network profiles of the dACC during basic motor behavior.
Executive control is a central organizing principal underlying human brain network interactions and a ubiquitous element of higher human function (Royall et al., 2002; Banich, 2009). Motor control is an important and specific form of control that subserves multiple complex behaviors, especially the integration of motor responses with cognitive decisions (Paus, 2001). The dorsal anterior cingulate cortex (dACC) has been strongly implicated in adult motor control. Tasks requiring motor control (including both the inhibition and excitation of motor responses) consistently demonstrate robust activation of the adult human dACC.
Cognitive-control-related brain networks, typically studied using complex response inhibition tasks, appear to be functionally articulated by the time of adolescence (Mennigen et al., 2014). Although the dACC continues to develop into the second decade of life (Gogtay et al., 2004), its role in mediating basic motor control in adolescence has been relatively under-studied. Our goal here was to investigate network interactions of the adolescent dACC: its undirected and directed interactions with the supplementary motor area (SMA), a high-level constituent of the motor system that has been associated with the control of simple finger movements, and with the primary motor cortex (M1). We attempted this using a simple unimanual visuomotor integration task with basic motor control demands. The task is highly suitable for characterizing brain network effects in adolescence as it is readily performed with no developmental changes in performance. Moreover, evidence of systematic network effects of the dACC in adolescence provides suggestive evidence for the early establishment of brain network interactions subserving this basic domain in the adult. Our results indicate that the adolescent dACC undergoes task-related undirected functional coupling with both the SMA and M1, but sends task-specific directed signals specifically to the SMA. These results suggest hierarchical control-related interactions between the dACC and SMA as part of a large-scale executive brain network underlying simple motor control in adolescence that most likely continues into adulthood.
The anterior cingulate cortex is a structurally and functionally complex frontal region (Bozkurt et al., 2005; Margulies et al., 2007). It has been proposed that the dorsal part of the anterior cingulate cortex (dACC, also called the anterior midcingulate cortex) participates in the control of behavior by acting as a major interface between sensorimotor and cognitive processing (Paus, 2001). The anterior cingulate cortex is perhaps best known as a mediator of cognitive conflict (Botvinick et al., 1999; Cohen et al., 2000, 2005; Milham and Banich, 2005; Aarts et al., 2008; Brown, 2009). However, the dACC is connected with the frontal motor system in both humans (Barbas, 1986; Luppino et al., 1990, 1991; Paus et al., 1993; Koski and Paus, 2000; Schulz et al., 2011; Amiez and Petrides, 2014) and non-human primates (Matelli et al., 1991; Morecraft and Van Hoesen, 1992; Vogt et al., 1992, 2005; Devinsky et al., 1995; Picard and Strick, 1996, 2001; Dum and Strick, 2002; Morecraft et al., 2012). That the dACC is also involved in motor function is seen by evidence from experimental and clinical lesion, electrical stimulation, and PET studies (Paus, 2001). In monkeys, unilateral cingulate lesions produce contralateral motor neglect (Watson et al., 1973). In humans, excessive activity in the dACC in the cingulate epilepsy syndrome is associated with motor impairment (Mesulam, 1981; Devinsky et al., 1995). Other human neurological syndromes that impact the dACC also involve disruption of motor control or cognitive processes associated with motor events (Devinsky et al., 1995; Paus, 2001). Electrical stimulation of the dACC elicits simple and complex contralateral and bilateral movements in both monkeys and humans (Luppino et al., 1991). PET studies confirm that motor regions of the human dACC are activated during manual movement (Paus et al., 1993, 1998). The dACC remains neurodevelopmentally dynamic in adolescence (Fjell et al., 2012) though this does not preclude highly articulated interactions with other brain regions, in relatively simple sensori-motor domains such as basic motor function (Ordaz et al., 2013).
The dACC is ideally suited to serve as a motor control area of both the adolescent and adult human brain: its dense projections to various regions of motor cortex suggest that the structure is in a privileged topological position to modulate their activities. Nonetheless, relatively few functional Magnetic Resonance Imaging (fMRI) studies have questioned whether the dACC modulates other brain regions during tasks with simple motor control demands. Whereas this question has been studied in the context of response conflict (Fan et al., 2008), the closest precedent comes from the recent work of Schulz and colleagues (Schulz et al., 2011). With a seed in BA 32, they used the psychophysiological interaction (PPI) measure in a go-no-go task (Friston et al., 1997) to assess dACC modulation of other brain regions during response preparation. Their results provide evidence for the task-related involvement of regions of the dACC, the dPFC and the basal ganglia (but not the SMA). These investigations were conducted in adults, and whether a similar functional network structure is present in adolescence is not yet known.
Here we predicted a motor control function for the dACC in a highly simple behavioral paradigm: unimanual finger-movement in response to exogenous visual stimuli. The task itself does not contain explicit excitation or inhibition demands and (as our results show) is robustly performed by adolescents. To test involvement of the dACC in motor control, we investigated its functional interactions with the SMA, itself a prominent cortical motor region, and with M1. As is well known, M1 is the primary motor outflow region of the cerebral cortex. The SMA is hypothesized to be involved in controlling simple unimanual finger-movements that are coordinated by sensory stimuli (Romo and Schultz, 1987; Thaler et al., 1988; Picard and Strick, 2003; Grefkes et al., 2008; Witt et al., 2008). However, it is unknown whether SMA acts alone to control simple coordinated manual behavior, or whether it requires modulatory signals from other brain regions.
We investigated functional relations of the dACC with the SMA and M1 using both undirected and directed functional connectivity analysis of fMRI Blood-Oxygen-Level-Dependent (BOLD) time series data. Undirected functional connectivity analysis measures correlated fMRI BOLD activity between different brain areas (Biswal et al., 1995; Maldjian, 2001; Martuzzi et al., 2010), and was assessed for the dACC with the SMA and M1 in the task and at rest. To assess modulatory signaling from the dACC to the SMA and M1, we used directed functional connectivity analysis, which measures the predictability of fMRI BOLD activity in one brain area from that in another area by multivariate autoregression analysis (Bressler et al., 2008; Deshpande and Hu, 2012; Tang et al., 2012). By comparing directed functional connectivity in the task and at rest, we tested whether the dACC exerts task-specific modulation of the SMA and M1.
Our results indicate that the adolescent dACC exerts task-specific modulation of the SMA, but not M1, in visually coordinated unimanual finger movement. We suggest that this modulation reflects motor control mechanisms associated with the adolescent dACC, consistent with other studies supporting a role for the dACC in motor control (Bush et al., 1999, 2000; Ridderinkhof et al., 2004), and that its role continues into adulthood. More generally, the results directly support the concept of the dACC being crucially involved in the intentional motor control of behavior (Paus, 2001). Our focus on adolescence allows us to address the working hypothesis that the adolescent motor system evinces complex network interactions during relatively simple visuo-motor coordination.
Materials and Methods
Eleven (7 male, 4 female) healthy adolescent participants (age: 8–18 yrs., mean = 14 yrs.) participated in the study. The participants, all from the metro Detroit area, provided informed consent or assent under a protocol approved by the Human Investigative Committee at the Wayne State University School of Medicine, and were monetarily compensated for their participation. Participants were screened with Schedule for Affective Disorders and Schizophrenia for School-Age Children-Present and Lifetime Version (K-SADS-PL; Kaufman et al., 1997) to rule out the presence of psychiatric diagnoses. All participants were predominantly right handed as evaluated using the structured Neurological Evaluation Scale (Buchanan and Heinrichs, 1989).
Participants were instructed to tap the forefinger of their right hand as quickly as possible in response to a flashing white stimulus in the center of the display panel (RGB:255, 255, 255; extent: 34 × 32 mm; subtended visual angle: ~17º; duration: 100 ms). Four behavioral paradigms were employed, in which the stimulus was presented with frequency of either 1 Hz or 0.5 Hz, and the presentation had either Periodic or Pseudo-random intervals between stimuli. Inter-stimulus intervals (in s) for the Pseudo-random epochs (either 1 Hz or 0.5 Hz) were created by pseudo-randomly sampling values from Gaussian distributions (μ = 1.0 sand σ = 0.5 s or μ = 2.0 sand σ = 1.0 s). The lower bound on inter-stimulus intervals was 300 ms (exceeding typical lower limits in response latency). Stimulus onsets during Pseudo-random epochs were adjusted so that the average frequency of the elicited response (and therefore the number of elicited responses) was equal to the periodic counterpart. Finger responses were collected from the receptive surface (extent: 33 × 33 mm) of a fiber-optic response touchpad (Current Design Systems, Inc.) interfaced with the Presentation software package (Neurobehavioral Systems, Inc.). During the resting control condition, participants were instructed to fixate on a cross hair in the center of the field of vision. In the course of a single scan, participants alternated between epochs involving one of the four behavioral paradigms (30 s duration each, 8 epochs total) and short rest (10 s duration each, 11 epochs total). In addition, epochs 10 and 18 had long rest (30 s duration). Eight behavioral paradigm epochs (involving finger movement) formed the Task condition and 13 rest epochs (11 short and 2 long) formed the Rest condition.
Functional Magnetic Resonance Imaging (fMRI) BOLD data collection was conducted on a Bruker MedSpec 4T system running the Siemens Syngo console at the Vaetkiveckius Imaging Institute in the Wayne State University School of Medicine. The scanner conformed to quality control standards exceeding clinical systems by virtue of daily quality control and SNR assessment. Whole-head gradient echo-planar images were acquired using an 8-channel head coil (TR: 2 s, TE: 30 ms, matrix: 64 × 64, 24 slices, voxels: 3.8 × 3.8 × 4.0 mm). fMRI data were collected continuously across all the conditions. An entire scan lasted 6.83 min.
Whole-head functional volumetric images were preprocessed with SPM5 using a standard protocol (Friston et al., 1995), including realignment to the first image in the series, correction for susceptibility-by-movement interactions, normalization to a standard EPI template, and smoothing with an 8 mm full width at half maximum (FWHM) isotropic Gaussian kernel.
In first-level analyses on individual participant data, windows of interest treated as boxcar waveforms (30 s) were convolved with the canonical hemodynamic response function (HRF) to produce reference waveforms. These waveforms were used for later contrast assessment within the General Linear Model framework. Motion effects were modeled using six movement parameters (for translation and rotation) as covariates of no interest.
Three Regions of Interest (ROIs) were considered for connectivity analysis: (1) Primary Motor Cortex, M1 (BA 4); (2) Supplementary Motor Cortex, SMA (medial BA 6) and (3) dACC (BA 24 and 32). Masks for the ROIs were generated based on the MNI atlas in the Talairach Daemon database using the WFU_PickAtlas software (ANSIR Laboratory) and applied to the whole-head functional volumetric image of each participant. As the Pickatlas mask does not specifically identify the region we label as the dACC, we identified those portions of the structure that have been most closely associated with motor control (Paus, 2001), and in separate studies have been referred to as the “midcingulate” cortex (Hoffstaedter et al., 2014, 2015). This region is also presumed to be functionally parcellated into at least two sub-divisions associated with motor function: a caudal cingulate zone and a rostral cingulate zone. The latter has been further divided into the anterior and posterior aspects, of which the latter has been more strongly associated with motor function (Milham and Banich, 2005). The dACC ROI, created in stereotactic space using BA 24 and the supra-genual aspects of BA 32 subsumed these divisions, and was anatomically distinct from, and non-overlapping with, the SMA.
Within each ROI, an effects of interest contrast was used to extract eigenvariate time series from voxels involved in the task. This contrast was based on the significance of the F statistic from the comparison of all task conditions, representing the weighted means of modeled effects, and was robust under the heterogeneity of the task conditions. Values were averaged over voxels within a 5 mm radial sphere (covering a volume of ~524 mm3) centered on the peak of the F contrast significance in each ROI. The average provides a relatively stable estimate of the representative signal around the significance peak. The sphere was located in the left hemisphere (contralateral to the finger movement behavior) for the SMA and M1, and in the hemisphere showing the higher peak of F contrast for the dACC. The resulting average eigenvariate intensity over a sequence of images represented the ROI time series that were then subjected to time series analysis.
Time Series Preprocessing
Each ROI time series was first z-normalized at each time point by subtracting the mean value (over all time points) and dividing by the standard deviation (over all time points). Images with outlier intensity values were rejected. Consistent with standard practice, the first four images in the fMRI acquisition were discarded to allow signal stabilization, and did not enter in subsequent analyses. The use of motion parameters as covariates in the first level analyses removed motion effects from subsequent analyses, ensuring that our data were free of outliers from scanner inhomogeneities or motion.
Undirected Functional Connectivity Analysis
The Pearson product-moment correlation coefficient was computed between the fMRI BOLD time series of the dACC with that of the SMA and of M1, in each participant separately for the Task (combined Periodic and Pseudo-random) and Rest conditions, and then also for the two task conditions (Periodic and Pseudo-random). All correlation values were computed using the cor function in the R software environment.
Directed Functional Connectivity Analysis
Directed functional connectivity was estimated from MultiVariate AutoRegressive (MVAR) models. For both the dACC with the SMA and with M1, separate MVAR models were estimated from the fMRI BOLD time series of the ROIs for the Task (combined Periodic and Pseudo-random) and Rest conditions in each participant. The MVAR model order (Bressler and Seth, 2011), indicating the number of previous time points in the model used to estimate a current time point, was one (Tang et al., 2012).
The strength of directed functional connectivity from one ROI to another was estimated by the magnitude of a t-statistic from significance testing of the corresponding MVAR model coefficient implemented in the glm function in the R software environment. This strength thus served as a metric for the dynamic causal relationship between the time series, and is similar to Granger Causality in being derived from the MVAR model (Seth and Edelman, 2007; Bressler and Seth, 2011). For each participant, directed functional connectivity was estimated in both directions between the two ROIs of each model, in the Task (combined Periodic and Pseudo-random) and Rest conditions, and then also for the two task conditions (Periodic and Pseudo-random).
Statistical Analysis, Effect Sizes, and Post hoc Tests
In the undirected analysis, we compared correlations between Task and Rest conditions using both the parametric paired t-test and the nonparametric Mann-Whitney-Wilcoxon signed-rank test. We also compared correlations between the Periodic and Pseudo-random task conditions using both the paired t-test and the signed-rank test. In the directed analysis, Task vs. Rest, and Periodic vs. Pseudo-random comparisons were performed using two-way, repeated-measures, within-participant Analysis of Variance (ANOVA), with Condition (Task and Rest, or Periodic and Pseudo-random) and Direction (dACC → SMA or M1, and SMA or M1 → dACC) as independent variables, and MVAR coefficient t-values as dependent variables. ANOVA was performed using the aov function in the R software environment. The α level for determining the significance of F-values was the generally accepted value of 5% (p < 0.05).
The effect size of a significant F-value for a factor or interaction was determined by computing η2 (effect sum of squares/ total sum of squares), which represents the fraction of the total variance in the dependent variable of the ANOVA that is attributable to the factor or interaction independent variable. We interpreted η2 as representing the proportion of variance in the dependent variable accounted for by a main factor or interaction (Tabachnick and Fidell, 2001).
Post hoc paired t-tests were used to analyze the results when significant main effects or interactions were found by ANOVA. The significances of these post hoc tests were corrected for multiple comparisons using Dunn’s multiple comparison procedure (based on the Bonferroni inequality). To control for the family-wise error rate, the alpha level for corrected t-values was the generally accepted value of 5% (p < 0.05). The effect size of a significant t-value was determined by computing the effect size (Pearson) correlation (r) (Cohen, 1988).
Response data were analyzed to assess patterns of missed responses and latencies for made responses. The percentage of missed responses and latencies for made responses were analyzed in separate repeated measures analyses of covariance with frequency and periodicity as within-participant factors, and age as covariate. No effects reached significance for the percentage of missed responses, indicating that these responses were not driven by frequency or periodicity. Separate regression analyses were used to assess age-related dependencies in the percentage of missed response. No significant relationship was observed for responses during Periodic epochs. During Pseudo-random epochs, there was a significant negative relationship between the percentage of missed responses and age, [F(1,8) = 3.99, p < 0.05, one-tailed, r2 = 0.33], suggesting that the frequency of missed responses during Pseudo-random epochs decreased with age. For the analyses of latencies, there was a significant main effect of periodicity [F(1,7) = 18.69, p < 0.003, MSe = 785.45]. Age-related dependencies assessed with regression analyses did not reveal significant effects. Response latencies during Pseudo-random epochs were on average 88 ms longer than during Periodic epochs.
Figure 1 depicts a conjunction analysis (Nichols et al., 2005) performed over the Periodic and Pseudo-random epochs and thresholded in core ROIs in the motor cortex including the dACC, the SMA and M1 (p < 0.05, cluster level). These clusters represent overlapping patterns of activation across each of the task periods. The peak activations in the dACC converge with regions previously labeled the posterior rostral anterior cingulate (Schulz et al., 2011), and associated with motor response related activation and control. More generally these ROIs are consistent with other studies of exogenously guided finger responses (Witt et al., 2008), and are notably present in contra- and ipsi-lateral M1 consistent with other evidence of mutually supplementary activations in ipsi-lateral M1 during unimanual responses (Newton et al., 2005). The substantial overlap across the ROIs highlighted in this activation-based analysis, emphasizes the value of deriving more specific (and differentiable) network interactions driving these activations (Friston, 2005).
Figure 1. Task-related BOLD activation maps. Activation patterns are shown for the conjunction of the two task variants (Periodic and Pseudo-random) averaged across participants. Significant activation clusters (p < 0.01, cluster level) are overlaid on a mosaic of axial views that reveal task-related activations in the dACC, SMA, M1 and the cerebellum.
Undirected Functional Connectivity
Comparisons of Task vs. Rest conditions on the undirected functional connectivity (correlation) between the dACC and the SMA were not significant by either the parametric paired t-test or the nonparametric Mann-Whitney-Wilcoxon signed-rank test. The same was true between the dACC and M1. By contrast, comparison of Periodic and Pseudo-random task conditions revealed that the correlations both between dACC and SMA, and between dACC and M1, were significantly greater in the Periodic than the Pseudo-random condition by both paired t-test [dACC and SMA: t(10) = 3.08, p < 0.05 (corrected); dACC and M1: t(10) = 3.41, p < 0.05 (corrected)] and by signed-rank test [dACC and SMA: V = 58, p < 0.05 (corrected); dACC and M1: V = 64, p < 0.05 (corrected)]. The effect size for dACC with SMA was r = 0.32. For dACC with M1, it was r = 0.47.
The finding of a significantly larger undirected functional connectivity between the dACC with both the SMA and M1 in the Periodic than the Pseudo-random condition suggests that the dACC activity was more tightly coupled with both these areas when the stimulus was more predictable as compared to when it was less predictable, but does not indicate the direction of influence underlying this coupling.
Directed Functional Connectivity
To test for a significant difference between Task and Rest conditions in directed functional connectivity between the dACC and the SMA or M1, two-way Analyses of Variance (ANOVAs) were performed on the MVAR coefficient values, with Condition (Task or Rest) and Direction (dACC → SMA, SMA → dACC or dACC → M1, M1 → dACC) as factors. For the dACC with the SMA, neither the Condition nor Direction main factors were significant, but their interaction was significant [F(1,10) = 11.28, p < 0.01, η2 = 0.07]. This value of η2 indicates that 7% of the total variance in the dependent variable was attributable to the interaction. For the dACC with M1, no significant main factor or interaction was observed.
The finding of a significant interaction for the dACC with the SMA between Condition and Direction factors suggests that the difference between Task and Rest conditions significantly varied with direction. To examine this further, post hoc paired t-tests were run to compare the two conditions separately for each direction. These tests revealed that the influence from the dACC to the SMA significantly differed between Task and Rest conditions [t(10) = 3.69, p < 0.01 corrected, effect size r = 0.41], but the influence from the SMA to the dACC did not. The effect size r indicates that there was a medium to strong effect between Task and Rest conditions for the influence from the dACC to the SMA (Cohen, 1988). The positive t-value indicates that this influence was greater for the Task than the Rest condition. This result is illustrated by the group means in Figure 2 and by the box plots in Figure 3.
Figure 2. The group mean and standard error of influence in both directions between the 2 ROIs (SMA to dACC, dACC to SMA) for Task (blue) and Rest (red) conditions. Post hoc paired t-tests revealed that the influence from dACC to SMA was significantly greater for Task than Rest, but that the influence from SMA to dACC did not significantly differ (**p < 0.01).
Figure 3. The distribution across subjects for both directions between the 2 ROIs (SMA to dACC, dACC to SMA) and for Task (blue) and Rest (red) conditions, as box plots corresponding to the bar plots in Figure 2. Note that the distribution for the influence from dACC to SMA during the Task is both more compact than and elevated above the distribution during Rest.
We then tested for a significant difference between Periodic and Pseudo-random conditions in directed functional connectivity between the dACC and the SMA, and between the dACC and M1. Again, two-way Analyses of Variance (ANOVAs) were performed, with Condition (Periodic or Pseudo-random) and Direction (dACC → SMA, SMA → dACC or dACC → M1, M1 → dACC) as factors. For the dACC with the SMA, both the Condition factor [F = 14.13, p < 0.005, η2 = 0.20] and the Condition x Direction interaction [F = 5.64, p < 0.05, η2 = 0.026] were significant. Post hoc tests revealed that the Periodic condition was significantly greater than the Pseudo-random condition for the influence from the dACC to the SMA [t(10) = 4.85, p < 0.001, effect size r = 0.54], but that the difference was not significant in the reverse direction. This result is illustrated by the group means in Figure 4 and by the box plots in Figure 5.
Figure 4. The group mean and standard error of influence in both directions between the 2 ROIs (SMA to dACC, dACC to SMA) for Periodic (dark blue) and Pseudo-random (light blue) task conditions. Post hoc paired t-tests revealed that the influence from dACC to SMA was significantly greater for Periodic than Pseudo-random conditions, but that the influence from SMA to dACC did not significantly differ (**p < 0.01).
Figure 5. The distribution across subjects for both directions between the 2 ROIs (SMA to dACC, dACC to SMA) and for Periodic (dark blue) and Pseudo-random (light blue) task conditions, as box plots corresponding to the bar plots in Figure 4. Note that the distribution for the influence from dACC to SMA during the Periodic condition is both more compact than and elevated above the distribution during the Pseudo-random condition.
The finding of a significant Condition factor, a significant interaction between Condition and Direction factors, and a significantly greater influence for Periodic than Pseudo-random conditions from dACC to SMA but not from SMA to dACC, suggest that the elevated undirected coupling observed between these two regions in the Periodic condition was due to the influence from dACC to SMA. The results also suggest that the observed Task vs. Rest difference in the dACC to SMA influence was due more to the Periodic than the Pseudo-random task trials.
We next used regression analysis to test whether the influence from dACC to SMA was significantly dependent on age. The F-value resulting from the regression analysis was not significant for either the Task or Rest condition.
For the 2-way ANOVA between the dACC and M1, the Condition factor was significant [F = 8.72, p < 0.05, η2 = 0.16], but neither the Direction factor nor the Condition x Direction interaction was significant. This finding indicates that the Periodic condition was significantly different from the Pseudo-random condition, regardless of direction. In fact, the mean influence, regardless of direction, was greater for the Periodic than Pseudo-random condition (Figures 4, 5). This finding that the directional functional connectivity between dACC and M1 was greater in the Periodic than the Pseudo-random condition, regardless of direction, is consistent with that from the undirected functional connectivity between dACC and M1. It suggests that the Direction factor was not important for the functional relation between these two regions.
In this study, we used time series analysis to investigate undirected and directed functional interactions of the adolescent dACC with the Supplementary Motor Area (SMA) and the Primary Motor Cortex (M1) during a simple unimanual visuomotor task. Undirected functional connectivity analysis revealed significantly greater coupling of dACC with both SMA and M1 for a proactive than a reactive task condition. Directed functional connectivity analysis demonstrated a significant task-specific modulatory influence from the dACC to the SMA, but not in the reverse direction, and in neither direction between the dACC and M1. Although our results were obtained in adolescents, they motivate the search for comparable analyses in adults. The employed behavioral paradigm is a very basic sensorimotor coordination task well performed by the adolescent subjects and un-confounded by development effects on task-proficiency.
The participants in this study were healthy adolescents. Even though cognitive control remains dynamic across adolescence, basic mechanisms of sensorimotor integration are relatively mature by then (Witt and Stevens, 2013). We demonstrate that the task-specific modulatory influence from the dACC to the SMA was independent of age, suggesting that motor-control-related network signaling by the dACC during simple tasks is sophisticated in adolescence. The maturing of dACC-related motor control thus may occur before maturation of higher order cognitive control. Undoubtedly further testing will be required in the future with adults performing the same task (and this is a focus of our ongoing studies).
Our analyses demonstrate a task-specific directed functional relation from the dACC to the SMA in unimanual visuomotor control, thus suggesting a specialized role for the dACC in human motor control. Although it is generally accepted that executive control depends on the coordination of sensory and motor operations by frontal regions of the brain (Funahashi, 2001; Miller and Cohen, 2001), different frontal regions are responsible for different aspects of executive control.
The finding that the dACC→SMA influence was significantly larger for the Periodic than Pseudo-random task conditions suggests that the strength of this influence is greater when the subject is in a proactive (Periodic) as compared to a reactive (Pseudo-random) state. This difference may reflect the maintained expectancy of stimulus timing that is inherent to proactive states. Maintaining the system in a proactive, as opposed to a reactive response state, presumably requires ongoing adjustments in control. Control and behavioral adjustments, particularly during conflict robustly activate the anterior cingulate cortex (Botvinick et al., 1999; Kerns et al., 2004; Sohn et al., 2007; Fan et al., 2008). Whereas the specific demands induced by response conflict differ from the demands of maintaining a response set during finger tapping, we surmise that our results reveal that the dACC is involved in the continual fine tuning of finger responses during a proactive response set. This result is an extension of previous work (Kerns et al., 2004), while providing evidence of network profiles of the dACC to the SMA.
Our results provide a measure of continuity with the work of Grefkes and colleagues (Grefkes et al., 2008). Using a task with manual (left or right) or bimanual responses to a flashing circle (1.5 Hz), they employed dynamic causal modeling (DCM) to assess connectivity and dynamics in a sub-network involving bilateral SMA, M1 and pre-motor cortex. They showed that the intra hemispheric SMA → M1 pathway was positively modulated when responding with the contra-lateral response hand, whereas the bilateral SMA → M1 pathways were modulated during bimanual responses. Although the dACC was not included in their network analyses, their results suggest a hierarchical organization within the motor system. By comparison, our results (though achieved using a complementary analytic method) suggest that the dACC, a region that is associated with more general mechanisms of cognitive control, exercises substantial modulation of the SMA as a function of motor task demand. These observations emphasize the value of effective connectivity techniques in revealing the complex and “causal” dynamics of the motor system.
By employing both undirected and directed functional connectivity, our study demonstrates the relative value of these two types of measure. Using undirected functional connectivity analysis, elevated coupling was observed between dACC and both SMA and M1 in the proactive, as compared to the reactive, condition. Directed functional connectivity analysis was required to disambiguate the pattern of directed interactions underlying the undirected functional coupling: that between dACC and SMA was supported by a unidirectional influence from dACC to SMA, whereas that between dACC and M1 was supported by bidirectional influences between the two areas.
The statistical bases of time series analyses do not allow us to determine the anatomical pathway over which modulatory signals might be conveyed between areas. This, and other limitations associated with the BOLD signal, means that we are unable to determine whether the presumed unidirectional modulation of the SMA by the dACC is monosynaptic or polysynaptic. Nonetheless, our findings are consistent with other evidence that the dACC provides supplementary functionality to the SMA (Hatanaka et al., 2003).
The application of autoregressive time series analysis to fMRI BOLD data has been controversial because of the possibility of inter-regional variation in the HRF (David et al., 2008; Friston, 2009, 2011; Roebroeck et al., 2011a,b; Smith et al., 2011; Valdes-Sosa et al., 2011; Stephan and Roebroeck, 2012; Friston et al., 2013), or measurement noise (Seth et al., 2013). However, neither of these potential problems can account for the dependence of the dACC-to-SMA influence on experimental condition in our analyses, since neither HRF nor measurement noise is expected to differ with experimental condition.
The observation that only the directional influence from the dACC to the SMA significantly distinguished task from rest suggests that the dACC may reside at a hierarchically higher organizational level than the SMA in motor control. Since the SMA is a premotor region, this result supports a triune cortical hierarchy of frontal executive networks (Fuster, 2003). In this view, the dACC is hierarchically above the SMA, and selectively modulates it as part of the SMA’s function in controlling simple coordinated manual behavior.
Our results inextricably link the dACC with the motor system, suggesting that the SMA is not alone in driving motor behavior. They open the possibility for differential contributions by the dACC and the SMA to motor behavior, such that the dACC might provide high-level control of the SMA. However, our analyses do not preclude the possibility of other “higher-order” areas being involved. For example, to maintain visually coordinated motor behavior in the present task, the dACC would need to be informed by the visual system. Our analyses did not allow us to infer whether the dACC receives direct projections from visual cortex (Sepulcre, 2014), or whether visual information is first passed to another “higher-order” area, which then informs the dACC. More detailed experimentation is needed to determine all the areas involved in motor control and the various conditions of their involvement.
We recorded and analyzed fMRI BOLD data from healthy adolescent participants performing a visually coordinated unimanual finger-movement task. The observation of task-specific modulatory signals to the SMA from the dACC indicates that the dACC is systematically involved in motor control in adolescence, consistent with its predicted role in motor control. Adolescent neuro- and behavioral-development is heterochronous, and brain interactions in the adolescent state are not necessarily predictive of the adult state. Nevertheless, the lack of age-related effects in our analyses suggests that both the behavioral domain, and the functional brain responses it evoked were invariant. It is plausible, though as yet unknown, that the combination of task and time series analyses we applied will reveal similar principles of motor organization in adults. We suggest that control mechanisms in core and extended motor circuitry may be hierarchically organized, with the dACC operating at a higher hierarchical level of motor control than the SMA.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The work was supported by the National Institute of Mental Health (MH 68680), the Children’s Research Center of Michigan (CRCM), the Children’s Hospital of Michigan Foundation, the Cohen Neuroscience Endowment, a Career Development Chair from the Office of the President (Wayne State University), and the Lyckaki-Young Fund from the State of Michigan. The authors declare no competing financial interests.
Aarts, E., Roelofs, A., and van Turennout, M. (2008). Anticipatory activity in anterior cingulate cortex can be independent of conflict and error likelihood. J. Neurosci. 28, 4671–4678. doi: 10.1523/jneurosci.4400-07.2008
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Botvinick, M., Nystrom, L. E., Fissell, K., Carter, C. S., and Cohen, J. D. (1999). Conflict monitoring versus selection-for-action in anterior cingulate cortex. Nature 402, 179–181. doi: 10.1038/46035
Bozkurt, A., Zilles, K., Schleicher, A., Kamper, L., Arigita, E. S., Uylings, H. B. M., et al. (2005). Distributions of transmitter receptors in the macaque cingulate cortex. Neuroimage 25, 219–229. doi: 10.1016/j.neuroimage.2004.10.040
Bressler, S. L., Tang, W., Sylvester, C. M., Shulman, G. L., and Corbetta, M. (2008). Top-down control of human visual cortex by frontal and parietal cortex in anticipatory visual spatial attention. J. Neurosci. 28, 10056–10061. doi: 10.1523/jneurosci.1776-08.2008
Brown, J. W. (2009). Conflict effects without conflict in anterior cingulate cortex: multiple response effects and context specific representations. Neuroimage 47, 334–341. doi: 10.1016/j.neuroimage.2009.04.034
Buchanan, R. W., and Heinrichs, D. W. (1989). The neurological evaluation scale (NES): a structured instrument for the assessment of neurological signs in schizophrenia. Psychiatry Res. 27, 335–350. doi: 10.1016/0165-1781(89)90148-0
Bush, G., Frazier, J. A., Rauch, S. L., Seidman, L. J., Whalen, P. J., Jenike, M. A., et al. (1999). Anterior cingulate cortex dysfunction in attention-deficit/hyperactivity disorder revealed by fMRI and the counting stroop. Biol. Psychiatry 45, 1542–1552. doi: 10.1016/s0006-3223(99)00083-9
Cohen, M. X., Heller, A. S., and Ranganath, C. (2005). Functional connectivity with anterior cingulate and orbitofrontal cortices during decision-making. Brain Res. Cogn. Brain Res. 23, 61–70. doi: 10.1016/j.cogbrainres.2005.01.010
David, O., Guillemain, I., Saillet, S., Reyt, S., Deransart, C., Segebarth, C., et al. (2008). Identifying neural drivers with functional MRI: an electrophysiological validation. PLoS Biol. 6, 2683–2697. doi: 10.1371/journal.pbio.0060315
Deshpande, G., and Hu, X. (2012). Investigating effective brain connectivity from FMRI data: past findings and current issues with reference to Granger causality analysis. Brain Connect. 2, 235–245. doi: 10.1089/brain.2012.0091
Fan, J., Hof, P. R., Guise, K. G., Fossella, J. A., and Posner, M. I. (2008). The functional integration of the anterior cingulate cortex during conflict processing. Cereb. Cortex 18, 796–805. doi: 10.1093/cercor/bhm125
Fjell, A. M., Walhovd, K. B., Brown, T. T., Kuperman, J. M., Chung, Y., Hagler, D. J., et al. (2012). Multimodal imaging of the self-regulating developing brain. Proc. Natl. Acad. Sci. U S A 109, 19620–19625. doi: 10.1073/pnas.1208243109
Friston, K. (2011). Dynamic causal modeling and Granger causality Comments on: the identification of interacting networks in the brain using fMRI: model selection, causality and deconvolution. Neuroimage 58, 303–305; author reply 310–311. doi: 10.1016/j.neuroimage.2009.09.031
Friston, K. J., Buechel, C., Fink, G. R., Morris, J., Rolls, E., and Dolan, R. J. (1997). Psychophysiological and modulatory interactions in neuroimaging. Neuroimage 6, 218–229. doi: 10.1006/nimg.1997.0291
Friston, K. J., Holmes, A. P., Poline, J. B., Grasby, P. J., Williams, S. C., Frackowiak, R. S., et al. (1995). Analysis of fMRI time-series revisited. Neuroimage 2, 45–53. doi: 10.1006/nimg.1995.1007
Gogtay, N., Giedd, J. N., Lusk, L., Hayashi, K. M., Greenstein, D., Vaituzis, A. C., et al. (2004). Dynamic mapping of human cortical development during childhood through early adulthood. Proc. Natl. Acad. Sci. U S A 101, 8174–8179. doi: 10.1073/pnas.0402680101
Grefkes, C., Eickhoff, S. B., Nowak, D. A., Dafotakis, M., and Fink, G. R. (2008). Dynamic intra- and interhemispheric interactions during unilateral and bilateral hand movements assessed with fMRI and DCM. Neuroimage 41, 1382–1394. doi: 10.1016/j.neuroimage.2008.03.048
Hatanaka, N., Tokuno, H., Hamada, I., Inase, M., Ito, Y., Imanishi, M., et al. (2003). Thalamocortical and intracortical connections of monkey cingulate motor areas. J. Comp. Neurol. 462, 121–138. doi: 10.1002/cne.10720
Hoffstaedter, F., Grefkes, C., Caspers, S., Roski, C., Palomero-Gallagher, N., Laird, A. R., et al. (2014). The role of anterior midcingulate cortex in cognitive motor control: evidence from functional connectivity analyses. Hum. Brain Mapp. 35, 2741–2753. doi: 10.1002/hbm.22363
Hoffstaedter, F., Grefkes, C., Roski, C., Caspers, S., Zilles, K., and Eikhoff, S. B. (2015). Age-related decrease of functional connectivity additional to gray matter atrophy in a network for movement initiation. Brain Struct. Funct. 220, 999–1012. doi: 10.1007/s00429-013-0696-2
Kaufman, J., Birmaher, B., Brent, D., Rao, U., Flynn, C., Moreci, P., et al. (1997). Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version (K-SADS-PL): initial reliability and validity data. J. Am. Acad. Child Adolesc. Psychiatry 36, 980–988. doi: 10.1097/00004583-199707000-00021
Kerns, J. G., Cohen, J. D., MacDonald, A. W. 3rd, Cho, R. Y., Stenger, V. A., and Carter, C. S. (2004). Anterior cingulate conflict monitoring and adjustments in control. Science 303, 1023–1026. doi: 10.1126/science.1089910
Luppino, G., Matelli, M., Camarda, R. M., Gallese, V., and Rizzolatti, G. (1991). Multiple representations of body movements in mesial area 6 and the adjacent cingulate cortex: an intracortical microstimulation study in the macaque monkey. J. Comp. Neurol. 311, 463–482. doi: 10.1002/cne.903110403
Luppino, G., Matelli, M., and Rizzolatti, G. (1990). Cortico-cortical connections of two electrophysiologically identified arm representations in the mesial agranular frontal cortex. Exp. Brain Res. 82, 214–218. doi: 10.1007/bf00230855
Margulies, D. S., Kelly, A. M. C., Uddin, L. Q., Biswal, B. B., Castellanos, F. X., and Milham, M. P. (2007). Mapping the functional connectivity of anterior cingulate cortex. Neuroimage 37, 579–588. doi: 10.1016/j.neuroimage.2007.05.019
Martuzzi, R., Ramani, R., Qiu, M., Rajeevan, N., and Constable, R. T. (2010). Functional connectivity and alterations in baseline brain state in humans. Neuroimage 49, 823–834. doi: 10.1016/j.neuroimage.2009.07.028
Matelli, M., Luppino, G., and Rizzolatti, G. (1991). Architecture of superior and mesial area 6 and the adjacent cingulate cortex in the macaque monkey. J. Comp. Neurol. 311, 445–462. doi: 10.1002/cne.903110402
Mennigen, E., Rodehacke, S., Müller, K. U., Ripke, S., Goschke, T., and Smolka, M. N. (2014). Exploring adolescent cognitive control in a combined interference switching task. Neuropsychologia 61, 175–189. doi: 10.1016/j.neuropsychologia.2014.06.022
Morecraft, R. J., Stilwell-Morecraft, K. S., Cipolloni, P. B., Ge, J., McNeal, D. W., and Pandya, D. N. (2012). Cytoarchitecture and cortical connections of the anterior cingulate and adjacent somatomotor fields in the rhesus monkey. Brain Res. Bull. 87, 457–497. doi: 10.1016/j.brainresbull.2011.12.005
Morecraft, R. J., and Van Hoesen, G. W. (1992). Cingulate input to the primary and supplementary motor cortices in the rhesus monkey: evidence for somatotopy in areas 24c and 23c. J. Comp. Neurol. 322, 471–489. doi: 10.1002/cne.903220403
Newton, J. M., Sunderland, A., and Gowland, P. A. (2005). fMRI signal decreases in ipsilateral primary motor cortex during unilateral hand movements are related to duration and side of movement. Neuroimage 24, 1080–1087. doi: 10.1016/j.neuroimage.2004.10.003
Ordaz, S. J., Foran, W., Velanova, K., and Luna, B. (2013). Longitudinal growth curves of brain function underlying inhibitory control through adolescence. J. Neurosci. 33, 18109–18124. doi: 10.1523/jneurosci.1741-13.2013
Paus, T., Koski, L., Caramanos, Z., and Westbury, C. (1998). Regional differences in the effects of task difficulty and motor output on blood flow response in the human anterior cingulate cortex: a review of 107 PET activation studies. Neuroreport 9, R37–R47. doi: 10.1097/00001756-199806220-00001
Paus, T., Petrides, M., Evans, A. C., and Meyer, E. (1993). Role of the human anterior cingulate cortex in the control of oculomotor, manual and speech responses: a positron emission tomography study. J. Neurophysiol. 70, 453–469.
Ridderinkhof, K. R., van den Wildenberg, W. P., Segalowitz, S. J., and Carter, C. S. (2004). Neurocognitive mechanisms of cognitive control: the role of prefrontal cortex in action selection, response inhibition, performance monitoring and reward-based learning. Brain Cogn. 56, 129–140. doi: 10.1016/j.bandc.2004.09.016
Roebroeck, A., Formisano, E., and Goebel, R. (2011a). The identification of interacting networks in the brain using fMRI: model selection, causality and deconvolution. Neuroimage 58, 296–302. doi: 10.1016/j.neuroimage.2009.09.036
Roebroeck, A., Formisano, E., and Goebel, R. (2011b). Reply to Friston and David: after comments on: the identification of interacting networks in the brain using fMRI: model selection, causality and deconvolution. Neuroimage 58, 310–311. doi: 10.1016/j.neuroimage.2009.10.077
Royall, D. R., Lauterbach, E. C., Cummings, J. L., Reeve, A., Rummans, T. A., Kaufer, D. I., et al. (2002). Executive control function: a review of its promise and challenges for clinical research. A report from the committee on research of the american neuropsychiatric association. J. Neuropsychiatry Clin. Neurosci. 14, 377–405. doi: 10.1176/appi.neuropsych.14.4.377
Schulz, K. P., Bedard, A. C., Czarnecki, R., and Fan, J. (2011). Preparatory activity and connectivity in dorsal anterior cingulate cortex for cognitive control. Neuroimage 57, 242–250. doi: 10.1016/j.neuroimage.2011.04.023
Seth, A. K., Chorley, P., and Barnett, L. C. (2013). Granger causality analysis of fMRI BOLD signals is invariant to hemodynamic convolution but not downsampling. Neuroimage 65, 540–555. doi: 10.1016/j.neuroimage.2012.09.049
Smith, S. M., Miller, K. L., Salimi-Khorshidi, G., Webster, M., Beckmann, C. F., Nichols, T. E., et al. (2011). Network modelling methods for FMRI. Neuroimage 54, 875–891. doi: 10.1016/j.neuroimage.2010.08.063
Sohn, M. H., Albert, M. V., Jung, K., Carter, C. S., and Anderson, J. R. (2007). Anticipation of conflict monitoring in the anterior cingulate cortex and the prefrontal cortex. Proc. Natl. Acad. Sci. U S A 104, 10330–10334. doi: 10.1073/pnas.0703225104
Tang, W., Bressler, S. L., Sylvester, C. M., Shulman, G. L., and Corbetta, M. (2012). Measuring Granger causality between cortical regions from voxelwise fMRI BOLD signals with LASSO. PLoS Comput. Biol. 8:e1002513. doi: 10.1371/journal.pcbi.1002513
Thaler, D. E., Rolls, E. T., and Passingham, R. E. (1988). Neuronal activity of the supplementary motor area (SMA) during internally and externally triggered wrist movements. Neurosci. Lett. 93, 264–269. doi: 10.1016/0304-3940(88)90093-6
Valdes-Sosa, P. A., Roebroeck, A., Daunizeau, J., and Friston, K. (2011). Effective connectivity: influence, causality and biophysical modeling. Neuroimage 58, 339–361. doi: 10.1016/j.neuroimage.2011.03.058
Vogt, B. A., Finch, D. M., and Olson, C. R. (1992). Functional heterogeneity in cingulate cortex: the anterior executive and posterior evaluative regions. Cereb. Cortex 2, 435–443. doi: 10.1093/cercor/2.6.435-a
Witt, S. T., Laird, A. R., and Meyerand, M. E. (2008). Functional neuroimaging correlates of finger-tapping task variations: an ALE meta-analysis. Neuroimage 42, 343–356. doi: 10.1016/j.neuroimage.2008.04.025
Keywords: dorsal anterior cingulate cortex (dACC), supplementary motor area (SMA), primary motor cortex (M1), motor control, connectivity analysis, adolescence
Citation: Asemi A, Ramaseshan K, Burgess A, Diwadkar VA and Bressler SL (2015) Dorsal anterior cingulate cortex modulates supplementary motor area in coordinated unimanual motor behavior. Front. Hum. Neurosci. 9:309. doi: 10.3389/fnhum.2015.00309
Received: 04 December 2014; Accepted: 14 May 2015;
Published online: 04 June 2015.
Edited by:Joydeep Bhattacharya, Goldsmiths, University of London, UK
Reviewed by:Theodore Zanto, University of California, San Francisco, USA
Lars Michels, University Hospital Zurich, Switzerland
Copyright © 2015 Asemi, Ramaseshan, Burgess, Diwadkar and Bressler. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vaibhav A. Diwadkar, Brain Imaging Research Division, Department of Psychiatry and Behavioral Neuroscience, Wayne State University School of Medicine, Suite 5A, Tolan Park Medical Building, 3901 Chrysler Service Drive, Detroit, MI 48201, USA, firstname.lastname@example.org