Abstract
The transition from childhood to adulthood is marked by pronounced functional and structural brain transformations that impact cognition and behavior. Here, we use a functional imaging approach to reveal dynamic changes in coupling strength between networks and the expression of discrete brain configurations over human development during rest and a cognitive control task. Although the brain's repertoire of functional states was generally preserved across ages, state-specific temporal features, such as the frequency of expression and the amount of time spent in select states, varied by age in ways that were dependent on condition. Increasing age was associated with greater variability of connection strengths across time at rest, while there was a selective inversion of this effect in higher-order networks during implementation of cognitive control. The results suggest that development is characterized by the modification of dynamic coupling to both maximize and constrain functional variability in response to ongoing cognitive and behavioral requirements.
Introduction
The abilities to maintain task goals and to filter distracting information are aspects of cognitive control (CC) that undergo profound change in early development (Zelazo, 2004; Luna et al., 2010; Somerville and Casey, 2010; Diamond, 2013). In the adult brain, these processes are rooted in the emergent dynamics of highly interconnected—and interacting—networks (Bullmore and Sporns, 2009). The spatiotemporal properties of the brain that constrain inter-regional dynamics, including anatomical coupling length, transmission time delays, and neurophysiological noise (Deco et al., 2009), change dramatically in early development as a consequence of a myriad of structural and functional alterations (for review, see Andersen, 2003; Brenhouse and Andersen, 2011) that impact both cognition and behavior. Indeed, numerical simulation and empirical electrophysiological investigations have highlighted the critical importance of internal dynamics for optimized brain function, and have shown, among other things, that age-related increases in behavioral speed and stability are associated with higher levels of neurophysiological noise and more complex brain responses to external stimulation (McIntosh et al., 2008). Thus, age-related differences in whole-brain emergent dynamics and their possible association with cognitive control deserve closer consideration as they could serve as important markers of developmental trajectories (indexing system integrity and function; Di Martino et al., 2014).
Age-related differences in inter-regional interaction—or functional connectivity (FC)—have been examined at a macroscopic, whole-brain level using neuroimaging techniques, and are linked to cognitive and behavioral maturation. Several intrinsic connectivity networks (ICNs) widely recognized from the study of adults—including prefrontal, default, sensory–motor, and salience networks—are evident early in development. These ICNs undergo topological reorganization over time (Supekar et al., 2009; Power et al., 2010), including the strengthening of long-range connections and the weakening of short-range connections within cognitive control networks (Fair et al., 2007), strengthening of the temporal segregation (anti-correlation) of task-positive and default networks (Marsh et al., 2006; Thomason et al., 2008), and a shift in location of cortical hubs from primary sensory-motor to heteromodal cortex (Fransson et al., 2011). While important, this picture of development overlooks the dynamic nature of the brain's functional connectivity profile, owing in part to the assumption that FC profiles are spatiotemporally invariant during the measurement period. Analysis strategies more sensitive to inter-regional coupling dynamics reveal short-lived meta-stable FC states that are reproducible across subjects and that deviate substantially from canonical connectivity patterns (Liu et al., 2013; Allen et al., 2014). The goal of the present study then was to examine age-related differences in FC dynamics and their possible association with behavior.
We had two predictions. First, given that structural constraints on emergent brain dynamics change over development, we predicted age-related differences in the temporal expression of whole-brain connectivity states that comprise the functional repertoire. Second, given age-related advances in the ability to respond to unpredicted challenges and hypothesized links between flexible cognition and the complexity of evoked brain dynamics, we predicted that inter-regional connectivity dynamics that occur in association with efficient stimulus/response selection and stable task-set maintenance would be more complex among older compared with younger individuals.
Materials and Methods
All procedures were approved by the University Research Ethics Board for Health Sciences Research at the University of Western Ontario (London, ON, Canada) and are in accordance with the 1964 Declaration of Helsinki. Portions of the data used in this study have been published previously (Wilk and Morton, 2012; Wilk et al., 2012). An overview of the methodological procedure is illustrated in Figure 1.
Illustration of task and analysis steps. A, Subjects were scanned in two different conditions, with eyes open at rest while fixating and during a size-congruency task administered in the form of a conflict adaptation paradigm. In the latter, on individual trials, participants were presented with two digits that differed in physical and numerical size, and pressed a button corresponding to the location of the numerically larger digit. On compatible trials, the numerically larger digit was also physically larger. On incompatible trials, the numerically larger digit was physically smaller. Individual trials were administered in conditions that differed in the proportion of overall trials that were compatible (25% or 75%). B, Data across subjects and conditions were decomposed using Group ICA into C = 90 components, 56 of which were identified as ICNs. GICA1 back-reconstruction was used to estimate the time courses (Ri) and spatial maps (Si) for each subject. C, Dynamic FC was estimated as the series of correlation matrices from windowed portions (W) of Ri, resulting in a concatenated data matrix of all ICN-to-ICN paired correlation values over time. D, Estimation of k was performed to select an optimal cluster number for clusters and then, following steps outlined in Allen et al., 2014, k-means was performed on the concatenated data matrix, resulting in a labeling of each window to a cluster and cluster centroids (state) that can be back-projected into the original matrix format. E, Variance in the correlation strength of each IC-to-IC coupling pair (IC2ICVar) was computed separately for the two conditions.
Subjects.
Fifty-one individuals (29 males) ranging in age from 9.4 to 32.4 years (mean, 16.9 years; SD, 5.8 years) participated, including 32 children (≤18 years of age) and 19 adults (>18 years of age). Participants <18 years of age were recruited from the University of Western Ontario Child Development Participant Pool. Parents of these participants gave written consent to their child's participation, and children provided verbal assent. Participants ≥18 years of age were recruited from the University of Western Ontario undergraduate and graduate student populations, and provided written consent to their participation. All participants were right handed and had normal or corrected-to-normal vision.
Resting-state protocol.
Participants were instructed to remain fixated on a centrally presented cross for 6 min and allow their thoughts to flow freely.
Task paradigm.
In each of three 7.8 min runs, participants were administered 76 trials of a standard stimulus–response compatibility task (Henik and Tzelgov, 1982; Fig. 1A). On each trial, two Arabic digits differing in physical size and numerical magnitude were simultaneously presented for 1950 ms in white on a black background. Digits ranged in value from 1 to 9. Physically larger and smaller digits were presented in 60- and 30-point font, respectively. On compatible trials, the numerically larger digit was physically larger; on incompatible trials, the numerically larger digit was physically smaller. Participants selected, by means of a right-handed button press, the numerically larger stimulus. Individual trials were randomly jittered by means of an intertrial interval (ITI; range, 1500–4500 ms; mean, 3000 ms) and were ordered to ensure that 50% of trials were compatible. To improve the signal-to-noise ratio, 11 8000 ms null events were added to each run. A fixed random run, trial, ITI, and null event ordering was used for all participants. The paradigm was presented using E-Prime software (Psychology Software Tools).
Data acquisition.
Data were collected using a 3 T Siemens Tim Trio MRI system fitted with a Siemens 32-channel head coil. Before scanning, participants <18 years of age received training in a mock-scanner facility to help limit motion and anxiety during the actual imaging procedure. Four functional runs were collected from each participant in a single scanning session, a single 180-volume resting-state run (eyes open, fixation) and three 234-volume task runs. Functional volumes consisted of 32 slices (thickness, 3 mm) oriented to the anterior commissure–posterior commissure axis and were collected in an ascending interleaved order using a T2*-weighted echoplanar imaging sequence (TR, 2000 ms; TE, 30 ms; flip angle, 78°), with an in-plane matrix of 64 × 64 pixels and a 21.1 × 21.1 cm field of view, yielding a 3 × 3 × 3 mm voxel resolution. There were no gaps between slices. To assist in the visualization of functional analyses, a T1-weighted anatomical scan consisting of 192 slices (thickness, 1 mm) with an in-plane 256 × 256 matrix and a 21.1 × 21.1 cm field of view (yielding a 1 × 1 × 1 mm voxel resolution) was also collected as part of the same scanning session.
Data preprocessing.
fMRI volumes were corrected for motion, smoothed with an 8 mm Gaussian kernel, and normalized to the MNI template using SPM8.
Independent component analysis.
Preprocessed volumes of concatenated resting-state and task data were submitted to high-model order group ICA using GIFT software (http://icatb.sourceforge.net; Calhoun et al., 2001; Fig. 1B). Two principal components analysis steps reduced the data to 90 components. Group data were then unmixed into 90 maximally independent components (ICs) using the Infomax algorithm. A high model order was used to improve the detection of noise sources and to limit the grouping of related subnetworks. Individual subject components were then created through back-reconstruction. Fifty-six neurophysiologically plausible sources, or ICNs, were identified by two trained viewers through independent visual (VIS) inspection of the spatial maps of the group ICs (Fig. 2). Plausible sources included auditory (AUD), somatomotor (SM), subcortical (SC), CC, cerebellar (CB), VIS, and default mode (DN) ICs. IC time courses were then filtered using a Butterworth filter and were despiked by replacing values greater than ±3 SDs with values equal to ±3 SDs. The group-averaged correlation matrix was calculated among all ICs following a Fisher z-transformation and the mean calculation of each subject's set of scans (Fig. 3).
Spatial profiles of the 56 neurophysiologically plausible ICNs revealed by a high-model order group ICA displayed on the medial and lateral views of the cortical surface. ICNs are grouped, and boxes are color coded based on their assignment, which is shown in Figure 3. In cases where a view or hemisphere is missing, no voxels above the lower-bound threshold (z > 0.75) were present.
Behavioral assessment.
Task performance was measured in terms of response time (RT) variability and error rate separately for each participant. RT variability was based on all trials with correct responses and was computed as the average squared deviation of each individual response time from the mean of all response times. The error rate was computed as the total number of errors divided by the total number of trials.
Mean FC estimation.
The mean FC was estimated by cross-correlating the full time courses of all 56 neurophysiologically plausible ICNs (Fig. 3). Age-related differences in the mean FC were estimated as the partial correlation of age and mean FC, with subject motion as a controlling variable (Fig. 4).
Dynamic matrix computation.
Internetwork dynamic FC analysis was performed on the 56 filtered and despiked IC time courses by means of a sliding window technique (Fig. 1C). Using a window of 30 s width (15 TRs), the time points of all 56 IC time courses within the window were cross-correlated using Pearson's correlation coefficient, yielding a total of 1540 unique IC-to-IC correlations. For each participant, these correlations were computed a total of 750 times (150 times from rest data, 600 times from task data) by sliding the window along subject-specific time courses using a time step of 2 s (1 TR).
Dynamic state detection.
Concatenated matrices across all 51 participants were submitted to a k-means cluster analysis (Fig. 1D). The k-means solution was evaluated across values of k ranging from 5 to 20 using the silhouette metric. The k value that returned the highest silhouette score without overfitting the data was 17 (Fig. 1D). Of the 17 total states, 5 were found to be unique to individual participants and were not considered, and the windows assigned to these states omitted from future analysis. Thus, our analysis focused on 12 connectivity states that dynamically reoccurred over time and were observable across multiple participants.
The following four metrics assessed connectivity state expression separately for each participant during rest and task: (1) frequency (F), measured as the proportion of all windows classified as instances of particular states, and computed separately for each of the 12 states; (2) mean dwell time (MDT), measured as the average number of consecutive windows classified as instances of the same state, and computed separately for each of the 12 states; (3) intertransition interval, measured as the number of consecutive instances before a transition to a different state and computed as an average across all states; and (4) the number of transitions (NT), measured as the number of state transitions across rest and task conditions.
To test the association between age and connectivity state expression while controlling for interindividual differences in motion, we modeled frequency, mean dwell time, intertransition interval, and number of transitions each separately as weighted linear sums of age and motion averaged across rest and task. Motion was computed separately for each participant as the average root mean square of the three translation and three rotation estimates generated during motion correction. To examine the effect of condition (rest vs task) and the interaction of age and condition, connectivity state metrics were computed from rest and task data separately. The effect of condition was then tested by means of paired-sample t tests. To test the interaction of age and condition, state expression metrics were modeled as weighted linear sums of age and motion separately for rest and task. We then tested whether the magnitude of the β-coefficient for age differed across rest and task conditions. To test whether connectivity state expression was further explainable by behavioral performance and sex, we added RT variability, error rate, and sex to the models, and tested whether these variables explained any additional variance in state expression variables. Scatterplots showing the association of age with the frequency and mean dwell time of each state separately for rest (blue circles) and task (red circles) are shown in Figure 5.
Assessment of internetwork correlation variability.
For each pairwise ICN-to-ICN coupling, variability in coupling strength was estimated as the variance of the ICN-to-ICN correlation (Fig. 1E), and was computed separately for rest and task data. Associations among internetwork correlation variability, age, and behavioral performance were examined in the context of multiple regression models in which internetwork correlation variability was modeled as the weighted sum of age, condition (task vs rest), age × condition, RT variability, and error rate. To control for the possible contaminating effects of subject motion known to differ across age (Power et al., 2012; Satterthwaite et al., 2012; Van Dijk et al., 2012; Yan et al., 2013), subject-wise root mean square motion was included as a nuisance regressor in all models.
Results
Behavior
RT variability decreased with increasing age (r(49) = −0.35, p < 0.05). The error rate across all trials was unrelated to age (r(49) = −0.17, not significant).
Mean FC
Analysis of the mean FC revealed a canonical pattern of inter-regional connectivity, including positive on-diagonal correlations among functionally related ICs; positive off-diagonal correlations between SM, CC, and VIS ICs; and negative off-diagonal correlations with DN ICs (Fig. 3). Static/mean IC-to-IC FC was associated with age (Fig. 4). The strongest associations were positive, and were evident largely among CC ICs, and between SM and VIS ICs. Weaker negative associations were particularly evident between CC and DN ICs.
Static (mean) functional network connectivity correlation matrix of ICNs across task and rest. Labels are color coded to match Figure 2.
Age-related correlations of static (mean) functional network connectivity. The bottom triangle of the matrix displays the correlation of ICNs across task and rest with age. The top matrix shows the statistically significant correlations at FDR thresholds of q = 0.01 (yellow) and q = 0.05 (orange).
Connectivity states
Twelve states dynamically reoccurred over time and were observable across multiple participants [Fig. 5 (for labels, see Fig. 3)]. Participants expressed an average of 4.15 (SD, 1.67) and 4.25 (SD, 1.96) states, respectively, during rest and task. The number of states expressed at rest was uncorrelated with age, but the number of states expressed during a task was positively associated with age (r(49) = 0.3, p < 0.03).
Age- and condition-dependent state expression. State solutions are shown for 12 states that were expressed by more than one subject. Below each ICN-to-ICN state matrix are distribution plots for frequency (left) and MDT (right) across ages for rest (blue circles) and task (red circles). Bold text marked by asterisks indicates significant main effects or interactions. a, Age; c, condition; mot, motion; r, rest; t, task.
Of the 12 states that were observable across participants, one state (state A) closely resembled the mean FC, with strong positive correlations along the diagonal (i.e., within functional groupings), and strong off-diagonal anti-correlations between CC and DN networks (Fig. 3). The remaining 11 states differed markedly from this pattern with differences between states both within and between functional groupings. Among SM networks, coupling ranged from being highly integrated (e.g., state C) to relatively disintegrated (e.g., state J), with some subclustering involving right and left motor cortex, superior parietal lobule, paracentral lobule, and postcentral gyrus evident in state E. Among CC networks, there was a departure away from one broadly integrated pattern evident in state A to two subclusters, one that included dorsolateral prefrontal cortex (dlPFC), anterior insula, and anterior cingulate cortex, and a second subcluster that was largely composed of lateral prefrontal and superior parietal ICNs. This was particularly evident in states F, H, E, and J. There was also striking variation in the coupling of the anterior insula to other CC networks, as such coupling ranged from being strongly and uniquely negative (states K and B) to being more moderate in strength (state J). Finally, among DN ICNs, coupling ranged from being highly integrated (state E) to being highly disintegrated (states C, L, I), with some evidence of subclustering among predominantly posterior DN ICNs, including posterior cingulate, right and left angular gyrus, and the MTL in states F and K. All 11 states also showed marked departures from the off-diagonal coupling evident in state A or the mean connectivity pattern. Anticorrelations between CC and DN ICNs, for example, were attenuated in states F, H, E, K, and I, although more so for DN coupling with lateral prefrontal/superior parietal ICNs than coupling with dlPFC and cinguloinsular ICNs. Furthermore, coupling among SM, SC, and CB networks, which were moderately strong in state A, were markedly weaker in state C, but were stronger in states D and E. Together, cluster centroids revealed marked and nonrandom departures from mean internetwork connectivity patterns.
Effects of condition and age on connectivity state expression
Scatterplots showing the association of age with the frequency and MDT of each state separately for rest (blue circles) and task (red circles) are shown in Figure 5.
Frequency
Age was negatively associated with the frequency of states I, F, and G, and positively associated with the frequency of states B and C. Condition was associated with the frequency of states K, B, and C, with states K and C more frequent during rest than task, and state B more frequent during task than rest. Age and condition had an interactive effect on the frequency of states H, B, and C. For state C, there was an age-related increase in frequency at rest, but no effect of age during task. For state B, there was an age-related increase in frequency during task, but no effect of age during rest. And for state H, there was an age-related increase in frequency during rest, but no effect of age during task. Adding behavioral and sex predictors to the models revealed a positive association between RT variability and the frequency of state K at rest, and a negative association between RT variability and the frequency of state E during task. Error rate and sex did not explain any additional variance in the frequency of any state. With the exception of state I, where motion was positively associated with frequency, and possibly state H, where motion was negatively associated with the frequency of expression during task, none of these associations were explainable by differences in motion between subjects.
Mean dwell time
Age was positively associated with the MDT of states H, D, and C, and negatively associated with the MDT of state I. Condition was also associated with MDT, with longer dwell times at task than rest for states L and B. There was an interaction of age and condition on the MDT of state B, with age positively associated with MDT during task but unrelated at rest. Adding behavioral and sex predictors to the models revealed a positive association between RT variability and the MDT of state I at rest, a positive association between RT variability and the MDT of state H during task, and a positive association of state E and error rate during task. Sex did not explain any additional variance in the MDT of any state. With the exception of state I, in which motion was positively related to the MDT, none of these associations or interactions was explainable by differences in motion.
Number of transitions
Age was positively associated with NT during task. Thus, older participants showed more state transitions during task than did younger participants (Fig. 6A).
Temporal properties of state expression. The distribution of the number of transitions (top) and the time between transitions (bottom) of the 12 states shown in Figure 5 are plotted across ages for rest (blue circles) and task (red circles).
Intertransition interval
Age was negatively associated with ITI during task. Thus, older participants showed more rapid state transitions during task than did younger participants (Fig. 6B).
Variability in internetwork connectivity
For the majority of ICN-to-ICN coupling pairs, correlation variability at rest was greater among older than younger participants (Fig. 7A). A few selected pairs, mostly among AUD, SM, and selected DN networks, showed the opposite pattern, with internetwork correlation variability greater among younger compared with older participants. There was, however, a marked reversal of this overall pattern during task. Under these conditions, a majority of ICN-to-ICN coupling pairs showed lower correlation variability among older compared with younger participants, save for several individual connections (Fig. 7B).
Age-related differences in internetwork connectivity variability. A, B, Variance in the ICN-to-ICN coupling strength was estimated for each pairwise ICN-to-ICN coupling, separately for rest (A) and task (B). Each coupling strength variability estimate was then correlated with age. Positive correlations (warm colors) indicate that connectivity variability was greater among older than younger participants; negative correlations (cool colors) indicate that variability was greater among younger than older participants. C, A statistical test of the interaction of age and condition, FDR corrected for multiple comparisons (q = 0.01), reveals coupling pairs whose correlation was more variable at rest than task for older relative to younger participants. No coupling pairs were more variable at rest than task for younger participants. D, Correlation of ICN-to-ICN variability and RT variability, FDR corrected (q = 0.05). E, Correlation of ICN-to-ICN variability and error rate, FDR corrected (q = 0.05). F, Conjunction of D and E.
To bring the difference in age-related ICN-to-ICN correlation variability across rest and task into greater relief and to test the statistical significance of this effect across all 1540 unique coupling pairs, we modeled the correlation variability of each coupling pair as a weighted sum of age, condition, and the interaction of age × condition, with subject motion included as a nuisance regressor. The β-coefficient of the age × condition interaction term was expressed as a t statistic, and a correction for multiple comparisons was implemented via the false discovery rate (FDR) (q = 0.01). The results are displayed in Figure 7C.
Coupling pairs that exhibited an age × condition interaction of their correlation variability were widely distributed, although they were clustered most densely among CC ICs, including right and left intraparietal sulcus, bilateral middle frontal gyrus, and cinguloinsular cortex, and between CC networks and SM and VIS networks. Dense clustering was also evident between DN ICs and SM and VIS networks. In summary, a widely distributed set of coupling pairs showed greater correlation variability at rest than at task for older relative to younger participants, with the most pronounced age-related increases in the dynamic range of ICN-to-ICN coupling strength evidenced within and between cognitive control networks.
To examine the implications of dynamic variation in internetwork connectivity for behavior, we examined, in two separate regression models, the linear association of ICN-to-ICN variability during task on the one hand and RT variability and error rate on the other, with participant motion included as a nuisance regressor in both models. Results were corrected for multiple comparisons by means of FDR (q = 0.05). They appear in Figure 7, D and E, respectively. Associations between ICN-to-ICN variability and RT variability were uniformly positive and evident in a widely distributed set of coupling pairs. Although there is no simple means of summarizing the distribution, noteworthy clusters include DN (medial temporal lobe)–CC (medial frontal gyrus, superior parietal cortex) coupling pairs, as well as DN–DN, and DN–VIS coupling pairs. Associations between ICN-to-ICN variability and error rate were also uniformly positive, but were more sparsely distributed. A conjunction of RT variability and error rate matrices (Fig. 7F) suggested some convergence in the association of ICN-to-ICN correlation variability and maladapted behavior (i.e., high response time variability and high error rate), most notably in DN–CC coupling pairs.
Discussion
The present study examined the effects of age and task condition on dynamic whole-brain FC. The analysis revealed 12 recurring FC states, only one of which was similar to the mean FC. Remaining FC states departed in marked, but nonrandom ways from the mean FC pattern, with substantial internetwork correlation variability. Age-dependent changes were evident in the transitional dynamics of the states across rest and task, and in internetwork correlation variability, suggesting that changes to large-scale brain dynamics, the “dynome” (Kopell et al., 2014) or the “chronnectome” (Calhoun et al., 2014), may be a fundamental feature of functional brain development.
Dynamic changes in network topology are evident both locally, within the perimeter of traditionally defined network boundaries (Chang and Glover, 2010; Hutchison et al., 2013b; Kucyi and Davis, 2014), and globally, at a whole-brain level (Liu et al., 2013; Allen et al., 2014). While, in principle, it is possible that connectivity states reflect noise-driven variation around a putative mean state, aspects of the present findings argue against this view. First, several connectivity states identified here bear a marked resemblance to connectivity states identified previously (Allen et al., 2014), including the following: (1) a “canonical” state akin to mean FC (state A); (2) a “somatomotor disconnection” state marked by strong coupling among SM ICNs, disconnection between SM and SC ICNs, and disintegration among DN ICNs (state C); and (3) two states that, despite varied topology, showed widespread disconnection of all ICNs with DN ICNs (states D and E). Evidence that select states are reproducible across studies suggests that they constitute a core feature of the adult state repertoire, and are robust against variation in methods, subject-specific brain morphology, and data quality. Second, principled changes in state expression occurred in association with changes in cognition and behavior. During the performance of a cognitive control task, states linked to drowsiness and inattention (state C) decreased in frequency relative to rest, whereas states marked by strong coupling among prefrontal cognitive control and visual ICNs (state B) increased in frequency relative to rest. The functional relevance of state B is particularly noteworthy, given the marked coupling among several cognitive control ICNs, including dorsolateral prefrontal, anterior cingulate, anterior insular, and superior parietal cortices. These regions have not only been linked to efficient stimulus/response selection and stable task-set maintenance in traditional studies of brain activation (Wilk and Morton, 2012; Wilk et al., 2012), but also undergo profound functional reorganization over development (Crone and Dahl, 2012). Note that strong cognitive control ICN-to-ICN coupling was not exclusive to this state (see states A, D, and K). This suggests that cognitive control is associated with a whole-brain state configuration that includes a unique connectivity fingerprint across all ICs, and not solely task-evoked activity in select “task-relevant” regions. Together then, the findings support the idea that dynamic state expression is not a reflection of random variance, but is a principled neurophysiological response to changing cognitive demands.
The findings motivate a reconceptualization of the link between cognitive control and FC. Current models assume that cognitive control networks are time-invariant structures, measureable by the mean FC assessed while participants are at rest, and whose enduring topological features somehow constrain executive processes (Cole and Schneider, 2007; Fair et al., 2007). Our results instead suggest that these networks are, like cognitive control process themselves, transient and dynamic. Furthermore, if, as emerging evidence suggests, dynamic connectivity states are nonrandom patterns (Hutchison et al., 2013a), then mean FC, on which static network characterizations are based, is a mixture of discrete connectivity states, each potentially evoked or exploited to some degree by cognition. From this perspective, the link between cognitive control and FC may be more direct than has been considered to date. Rather than indirectly correlated with static network topology, cognitive control may be directly coincident with dynamic changes in FC.
The implications for development are clear. Foremost perhaps is that alterations in the brain's spatial–temporal coupling structure that occur over the course of development (Deco et al., 2009) impact FC dynamics, with several lines of evidence consistent with this prediction. There were, for example, age-related differences in the frequency of state expression, with older participants expressing states C and D—states known to comprise the adult repertoire (Allen et al., 2014)—more frequently than younger participants, and younger participants expressing states F and G more frequently than older participants. Age-related differences were also evident in the effect of condition (rest vs task) on state expression. Older participants, for example, expressed the somatomotor disconnection state—state C—more frequently at rest than during task, and expressed state B less frequently at rest than during task. Younger participants, by contrast, expressed states B and C with comparable frequency across rest and task. Condition also affected the mean dwell time of states B and H for older, but not younger participants. Importantly, age-related differences in dynamic state expression were evident even though the states themselves were topologically similar for participants of different ages (Fig. 8). Together, our findings suggest that it is dynamic state transitions rather than a static topology traced by inter-regional connections that is the locus of change in the years spanning late childhood and early adulthood.
Average state fit of children (<216 months) and adults (≥216 months). A–L, Bars represent the mean spatial correlation of all windowed patterns with the canonical state configuration to which it was assigned by k-means clustering, derived separately for children (blue) and adults (red). Error bars represent 1 SD.
Whether changes in state expression are a consequence of, or can even be predicted from, changes in the brain's spatial–temporal coupling structure, remains unclear. Anatomical coupling length, transmission time delays, and neurophysiological noise (Deco et al., 2009), fundamentally constrain emergent neural dynamics and undergo dramatic change during the developmental period spanned by the current sample. Examination of the relationship between changes in the brain's coupling structure and changes in state dynamics deserves close consideration in future studies.
The noisy brain
A second prediction concerned age-related differences in inter-regional coupling variability and the complexity of evoked brain dynamics. Although major structural pathways and functional features within the human connectome develop before 2 years of age (LaMantia and Rakic, 1990; Luo and O'Leary, 2005; Low and Cheng, 2006; Gao et al., 2011; Yap et al., 2011), changes in the size (i.e., diameter and length), relative alignment, packing density, and myelination of axons continue until early adulthood, and promote the efficacy of neural communication (Innocenti, 1995; Chen et al., 2002; Bütefisch, 2004; Paus, 2010). Numerical simulations suggest that these changes, together with age-related increases in neurophysiological noise, lead to greater complexity in evoked brain dynamics. Much as an athlete is best prepared for action when in motion, noise lends computational capacity to the dynamic brain by facilitating transitions between different multistable connectivity states given external stimulation. Greater complexity in evoked dynamics, in turn, affords greater cognitive and behavioral flexibility.
The present findings were broadly convergent with these ideas. At rest, ICN-to-ICN coupling variability was typically higher among older than younger participants, and brain dynamics were of comparable complexity. However, with the administration of a cognitive control task, there was a marked reversal of these effects. ICN-to-ICN variability was typically lower among older participants, and brain dynamics were decidedly more complex—the number of states expressed, the number of transitions occurring between multistable connectivity states, and the rate of transition between states were all higher among older than younger participants. Thus, as predicted, there was an age-related increase in the complexity of brain dynamics evoked by a challenging cognitive control task.
Variability in ICN-to-ICN coupling also had implications for behavior. During the administration of a cognitive control task, ICN-to-ICN variability among selected DN-to-DN, DN-to-CC, and VIS-to-CC couplings were positively associated with RT variability and error rate. These associations remained largely unchanged after controlling for differences in age. Associations between DN function and behavior are well documented. DN activity, for example, increases in the moments preceding error commission (Eichele et al., 2008), and is associated with momentary lapses of attention (Weissman et al., 2006). Individual differences in static DN connectivity have been associated with the tendency to engage in mental time travel (Andrews-Hanna et al., 2010). And, in close alignment with the present findings, dynamic variability in the strength of DN connectivity has been linked to the frequency of mind wandering during task performance, with greater variability associated with more frequent instances of mind wandering (Kucyi and Davis, 2014).
Limitations
A major concern with using functional imaging across development are age-associated motion artifacts (Power et al., 2012). We instituted several procedures to mitigate the effects of motion, including mock-scanner pretraining, stringent motion-related data exclusion criteria, ICA-based denoising, and the use of motion regressors in all models. Beyond motion, cardiac and respiratory signals may contribute to age-related differences as normative rates differ by age. ICA is relatively effective at removing physiological signals; however, due to aliasing and possible task-specific effects, their contribution to dynamic changes cannot be ruled out without faster sampling and/or concurrent physiological recordings.
Conclusions
Much as the mind rapidly adapts to ever-changing experiences, functional connectivity dynamically varies on relatively short time scales. The present results suggest that how the brain maximizes and constrains connectivity variability in response to cognitive and behavioral challenge changes profoundly over the course of development.
Footnotes
This research was supported by grants from the National Science and Engineering Research Council, the Canadian Foundation for Innovation, and the Ontario Innovation Trust to J.B.M.; and a Canadian Institute of Health Research postdoctoral fellowship and a Brain & Behavior Research Foundation NARSAD Young Investigator Grant to R.M.H. We thank Elena Allen for sharing code used in the analysis, and Katie Insel for comments on an earlier version of the manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. J. Bruce Morton, Cognitive Development and Neuroimaging Laboratory, University of Western Ontario, London, ON N6A 3K7, Canada. bmorton3{at}uwo.ca or jbrucemorton{at}gmail.com