Key Brain Network Nodes Show Differential Cognitive Relevance and Developmental Trajectories during Childhood and Adolescence

Abstract Human adolescence is a period of rapid changes in cognition and goal-directed behavior, and it constitutes a major transitional phase towards adulthood. One of the mechanisms suggested to underlie the protracted maturation of functional brain networks, is the increased network integration and segregation enhancing neural efficiency. Importantly, the increasing coordinated network interplay throughout development is mediated through functional hubs, which are highly connected brain areas suggested to be pivotal nodes for the regulation of neural activity. To elucidate brain hub development during childhood and adolescence, we estimated voxel-wise eigenvector centrality (EC) using functional magnetic resonance imaging (fMRI) data from two different psychological contexts (resting state and a working memory task), in a large cross-sectional sample (n = 754) spanning the age from 8 to 22 years, and decomposed the maps using independent component analysis (ICA). Our results reveal significant age-related centrality differences in cingulo-opercular, visual, and sensorimotor network nodes during both rest and task performance, suggesting that common neurodevelopmental processes manifest across different mental states. Supporting the functional significance of these developmental patterns, the centrality of the cingulo-opercular node was positively associated with task performance. These findings provide evidence for protracted maturation of hub properties in specific nodes of the brain connectome during the course of childhood and adolescence and suggest that cingulo-opercular centrality is a key factor supporting neurocognitive development.


Introduction
The gradual transition from childhood to adulthood is characterized by profound physical, cognitive, and emotional changes. Individual adaptation to the dynamic and socially expanding environment throughout this period is enabled by the large potential for neuroplasticity, e.g., reflected in a transformation of cortical gray matter by a massive synaptic downscaling (Spear, 2000) and increases in brain white matter coherence and volume (Paus, 2005;Westlye et al., 2010;Alnaes et al., 2018). These protracted neurodevelopmental changes form the structural backbone allowing efficient neural signaling, enabling functional integration and adaptation of the brain networks underlying the substantial cognitive and functional maturation seen over this age span (Casey et al., 2008). Of note, while adolescence marks the successful transition from childhood to adulthood for the majority, this is also the period where most mental disorders emerge, supporting the critical role of adolescence as a sensitive period during ontogeny (Paus et al., 2008). Functional brain connectivity (FC) provides an intriguing window into the developing brain, and it has been extensively studied using resting-state functional magnetic resonance imaging (rs-fMRI). Resting-state networks (RSNs) show high spatial concordance with networks associated with task activations (Smith et al., 2009), and brain connectivity and activation are highly coordinated aspects of brain functioning (Gratton et al., 2016). In a large-scale study on children and adolescents from the Philadelphia Neurodevelopmental Cohort (PNC), Satterthwaite et al. (2014) analyzed fMRI data recorded during the performance of a cognitive task and found that a reciprocal pattern of executive-network activation and default-mode network (DMN) deactivation was associated with working memory performance. This pattern was more predictive of performance differences than activations and deactivations in isolation. Network dynamics and connectivity thus contains complementary information about the neural substrate supporting executive ability and its maturation Gu et al., 2015).
Whereas the hubs in the default mode and executivenetworks are fairly established from five years of age (Hwang et al., 2013;Wu et al., 2013;Cao et al., 2014), key dynamic and functional characteristics of brain nodes display maturational alterations throughout adolescence (Grayson et al., 2014). How these changes support the marked increase in executive abilities during this period is still unclear.
One commonly used metric for investigating hub properties is centrality, which reflects the number and strength of connections of a given node to the rest of the network.
Previous studies have reported a maturational increase in centrality of task positive networks like frontoparietal, cingulo-opercular and task negative DMN nodes (Fransson et al., 2011;Hwang et al., 2013;Cao et al., 2016). While most developmental studies on brain network connectivity have employed rs-fMRI, cognitive engagement modulates the connectivity and centrality of brain nodes (Chadick and Gazzaley, 2011;Cole et al., 2013;Alnaes et al., 2015a,b;Mišić and Sporns, 2016;Kaufmann et al., 2017b). Further, higher cognitive load has been shown to increase the sensitivity for age-related differences in FC (Dørum et al., 2017), motivating investigation of brain dynamics beyond the unconstrained resting state.
Here, to characterize age-related differences in functional connectivity in the adolescent brain, we estimated whole-brain voxel-wise eigenvector centrality (EC) of fMRI data collected during resting state and a working memory task (fractal n-back) for 754 children and adolescents aged 8 -22 years from the PNC . To allow for multiple levels of inference and to investigate the hierarchical organization of the brain network nodes, we decomposed the voxel-wise EC maps using independent component analysis (ICA) and tested for main effects of task, age, and their interactions both on the voxel-wise EC maps and the ICA subject weights. Based on previous studies and models of adolescent brain development, we hypothesized (1) increased centrality in task positive networks and decreased centrality in DMN, sensorimotor, and visual networks in response to task engagement; (2) and that these centrality changes would be associated with task performance. Furthermore, (3) we hypothesized that these effects on network centrality would parallel task-related activation through a characteristic pattern of increased activation in networks involved in task engagement and decreased activation in DMN. Lastly, we hypothesized (4) that the centrality of key nodes would show evidence of protracted maturation during the sampled age-span, with differential age-related differences across nodes and conditions, possibly indicating increasing agedependent differentiation between the two psychological contexts for both activation and centrality.
Philadelphia. After exclusion of participants with severe medical and psychiatric conditions, missing task responses, and poor normalization of fMRI data, the final sample comprised 754 individuals aged 8 -22 years (405 males, mean 15 years old, SD 3.3 years). All subjects were recruited through the Center for Applied Genomics at The Children's Hospital in Philadelphia.

fMRI task paradigm and behavior
The task has been previously described in detail (Satterthwaite et al., 2013). Briefly, we included data from two fMRI runs: One was a resting-state run (eyes open while fixating on central fixation cross), while the other was collected while participants performed a fractal version of the n-back task. The n-back task consisted of three load levels: 0-back, 1-back, and 2-back. During 0-back, participants were asked to press a button every time a target fractal was presented. During 1-and 2-back, participants were asked to respond whenever the presented fractal was identical to the one presented one or two trials prior, respectively. Each run consisted of three blocks per condition, and each block consisted of 20 trials, of which five were targets and 15 were non-target fractals. We computed the total hit-rate, false-positive rate, and D-prime. The latter was calculated by subtracting z score for false positives from the z score for hit-rate from the in-scanner responses during the n-back task (Stanislaw and Todorov, 1999). fMRI data processing fMRI data were processed using FMRI Expert Analysis Tool (FEAT) version 6.00, from FMRIB's Software Library (FSL; Smith et al., 2004;Jenkinson et al., 2012), and included the following steps: correction for motion using MCFLIRT (Jenkinson et al., 2002), linear trend removal and high-pass filtering (0.01 Hz), removal of non-brain tissue using BET , spatial smoothing with a Gaussian kernel of full width at half maximum (FWHM) of 6 mm (SUSAN; Smith and Brady, 1997). We employed automated procedures for data denoising, including (FIX; Salimi-Khorshidi et al., 2014) and ICA-AROMA (Pruim et al., 2015), which has been reported to reduce risk of inflating age-related effects due to higher in-scanner head-movement (Kaufmann et al., 2017a). Further, we applied nonlinear registration using FNIRT to Montreal Neurologic Institute (MNI) 152 standard space using the T1-weighted scan as an intermediate.

EC mapping (ECM)
Following previous studies (Alnaes et al., 2015b;Skåtun et al., 2016), we first created a group mask containing common voxels across all participants and sessions (nback and rest). To accommodate images to EC estimation, the pre-processed and normalized n-back, resting state runs and analysis mask were resampled to an isotropic voxel resolution of 3 mm, converted to the Vista image format (Pope and Lowe, 1994) using LIPSIA (Lohmann et al., 2001) and submitted to estimation of voxelwise EC (Lohmann et al., 2010) using the absolute of the full correlation coefficients (equal weight to positive and negative correlations), yielding one EC map per run for each participant. Individual EC maps and brain masks were converted to nifti format, resampled to 2-mm isotropic resolution, and submitted to further analysis (ICA and voxel-wise statistics).

ICA on EC maps
The 1508 EC maps were decomposed into a fixed set of 40 components, for optimal trade-off between spatial resolution of components and not over-fitting the model, using spatial ICA in MELODIC (Beckmann and Smith, 2004). To assess effects of condition on EC for each IC, the associated subject-weights from the ICA was used to calculate the difference in subject weights between rest and n-back runs for each IC. For visualization of the network structure, we used subject weights for both n-back and rest to form a component by component correlation matrix, which was submitted to hierarchical clustering using FSLNets (Smith et al., 2011). ICs covering white matter and CSF were excluded before clustering.

Dual regression
To assess the relationships between the centrality and task-related activation of the ECM-based ICs, we extracted individual time series from each of the component spatial maps from the n-back data using dual regression (Nickerson et al., 2017). GLMs were estimated for each participant and each component's time series, modeling 0-back, 1-back, and 2-back as well as responses and instructions. As a measure of BOLD task effect, ␤-estimates for 0-back was subtracted from the average of 1-back and 2-back ␤-estimates.

Voxel-wise analysis
For transparency, we compared our main results obtained using ICA-based decomposition of EC maps with voxel-wise effects of task and age and their interactions on EC using FSL randomize (Winkler et al., 2014). First, we computed difference maps by subtracting the individual resting-state EC maps from the n-back EC maps. We tested for main effect of condition by performing a onesample t test on the difference maps (n-back minus rest). We further tested for associations with age by performing linear regression on rest, n-back, and the difference maps, including sex as covariate. Statistical inference was done by permutation testing with 5000 iterations and threshold-free cluster enhancement (Smith and Nichols, 2009).

Statistical analysis
Analyses beyond the voxel-and imaging domain were performed in MATLAB 2014a (MathWorks). To test for effects of age on EC, we performed linear regressions with the subject weights from n-back or resting-state runs as dependent variable, and age and sex as independent variables for all ICs. To test for main effect of task on each component, we performed one-sample t tests on the difference scores between rest and n-back weights. To test for interactions between task and age on EC, we performed linear regressions with the difference score as dependent variable and age and sex as independent variables. To test for associations between EC and task performance (hit-rate, false positives, and D-prime), we used the same derived difference score as dependent variable, and hit-rate, false positives, or D-prime as independent variables in separate analyses, in addition to age and sex.
In the same manner, using the ␤-estimates derived from dual regression on the n-back data, we tested for effect of task by performing one-sample t tests. To test for effect of age, we performed linear regression with task effect as the dependent variable and age and sex as independent variables.
To test for associations between task-related effects on centrality and activation for each IC, we correlated the difference EC-score (n-back minus rest) with ␤-estimates for task engagement. To account for multiple comparisons, we adjusted the false discovery rate (FDR, q ϭ 0.05) across all tests and components (Genovese et al., 2002). For visualization of the association between centrality and activation, component-wise EC-difference-weights and ␤-estimates for the n-back task was imported to R (http:// www.r-project.org) and plotted with the ggplot2-package (Wickham, 2009).

ICA and hierarchical clustering
ICA of EC maps yielded 31 spatial maps corresponding to cortical gray-matter components, including frontoparietal, default mode, visual networks, and nine maps corresponding to white matter and CSF (for an anatomic description of the 31 maps included in the hierarchical clustering, see Table 1). Figure 1,1-1 shows the hierarchical clustering based on the correlation between ICweights across n-back and rest, and Figure 1,1-2 shows the associated spatial maps. Clustering revealed a structure coherent with known large-scale functional networks, largely reflecting visual, sensorimotor, frontoparietal, subcortical, cingulo-opercular, and DMN.

Component-wise EC task effects
Task engagement (n-back Ͼ rest) was associated with distributed differences in EC across most nodes. Figure  1,1-4 displays significant (FDR corrected, q ϭ 0.05) associations between EC, task engagement and task performance, and Table 2 summarizes corresponding statistics. All frontoparietal components displayed a significant increase in EC with task engagement, and for one of the components (IC 2) higher centrality was significantly associated with better performance. Three of the four cingulo-opercular nodes displayed significant increase of centrality with task engagement, and higher EC was associated with better performance (two of four ICs). The subcortical components displayed a mixed pattern of centrality changes, thalamus decreased, but cerebellum increased with task engagement. Higher EC for both cerebellum and thalamus was associated with better performance. DMN components displayed both increase (two of five ICs) and decrease (two of five ICs) in centrality with task engagement, but no significant associations with task performance. Sensorimotor components displayed reduced centrality (five of six ICs) with task engagement, and for two of the components, the degree of reduction was associated with task performance. Visual components displayed a mixed pattern with both increased (three of seven ICs) and decreased (three of seven ICs) centrality with task engagement, and centrality in two components showed negative correlations with task performance.

Component-wise activation
Significant effects (FDR corrected, q ϭ 0.05) of task engagement are displayed in Figure 1,1-5, lower panel. Table 2 summarizes the corresponding statistics. The componentwise activation pattern revealed several similarities with the pattern for EC. Task engagement was associated with significantly increased activation in the frontoparietal (six of seven ICs) and cingulo-opercular (four of four ICs) components. Among the subcortical components, we found increased activation in cerebellum and decreased activation in putamen. The visual, sensorimotor, and DMN components displayed reduced activation during task.

Associations between centrality and activation
Overall, we found a positive correlation between activation and EC (r ϭ 0.64, p ϭ 1.11e Ϫ4 ): components showing increased activity with task generally also showed taskrelated increases in EC. Figure 2 displays the componentwise association between task-related activation and differences in EC between n-back and rest. However, the results revealed some notable exceptions. For example, whereas components in the frontoparietal cluster showed both strong task-related activation and increased centrality, and visual and sensorimotor cluster components showed deactivation and decreased centrality, DMN components showed overall task related deactivations, but a mixed pattern of increased and decreased centrality. Figure 1, 1-3 and 1-5, upper panel, shows the significant (FDR corrected, q ϭ 0.05) associations between age centrality and task-related activation, whereas Figure 3 shows the age-associated trajectory for selected components. Table 3 summarizes the corresponding statistics. For frontoparietal components, we found significant positive associations between centrality and age in one component during n-back and two components during rest. One of the frontoparietal components displayed a significant interaction effect suggesting a stronger negative age effect during rest compared to n-back. In contrast, we found significantly increased activation with increasing age across all frontoparietal components, with the exception of one. Cingulo-opercular components displayed sig-nificant age-related increases in EC in two of four components, anterior cingulate and insula respectively. The increase in EC was paralleled by significant increase in activation in all cingulo-opercular components. Among the subcortical components, thalamus and cerebellum showed significant age-related decreased centrality, and the latter also showed age-related increased activation. Within the DMN cluster, one component showed significant age-related decrease in EC during n-back. In contrast, all DMN-components showed significant decrease in activation with increasing age, indicating stronger deactivation with higher age. Visual and sensorimotor components displayed a highly similar pattern: age-related increases in EC during both task and rest, but decreased activation with higher age. The exception was one component, showing age-related increases both in centrality and activation. Figure 4 shows the results from the full-brain voxel-wise analyses. Figure 4,4-1 displays the difference in centrality during n-back compared to rest, suggesting higher centrality in frontoparietal areas and posterior DMN and lower centrality in sensorimotor areas, mirroring the ICA-based analysis. Figure 4,4-2 and 4-3, shows associations between centrality and age during n-back and rest, respectively. We found significant age-related increases in centrality in both conditions in visual and sensorimotor regions as well as in the anterior cingulate. During n-back, we found significant age-related differences in the frontal pole, extending to the superior frontal gyrus. Negative associations with age were found in large white matter and CSF regions. However, there was no significant age effect in the derived difference maps, indicating no evidence of interactions between age and condition.

Discussion
Using graph-based metrics and data-driven decomposition of fMRI data obtained during an unconstrained resting state condition and performance of an n-back task, the current study yielded three main findings. First, we demonstrated the sensitivity of EC to differences in task demands. Next, we revealed age-dependent alterations in the centrality of key brain networks, with increasing centrality in sensorimotor, visual, and cinguloopercular network components and decreasing centrality in thalamus, cerebellum, and left temporal gyrus. Finally, the lack of condition by age interactions suggests that these age-related differences in network centrality are not strongly dependent on mental state, supporting that the maturation of the functional architecture of the brain is relatively pervasive across psychological contexts. These main findings will be discussed below.

Brain network centrality and task engagement
The most prominent finding in the current study was increased centrality of frontoparietal, cerebellar and cingulo-opercular components, and decreased centrality in sensorimotor and thalamic components with task engagement, which align with previous investigations of EC in adults (Alnaes et al., 2015b;Hearne et al., 2017). FC during task engagement is associated with increased cross-network connectivity, and decreased network modularity (Di et al., 2013;Cole et al., 2014;Shine et al., 2016). This dynamic shift in connectivity pattern in response to task demands supports increased global effectiveness, through shortening path lengths in the network, and increased efficiency of information transfer across the global network (Stanley et al., 2015). These alterations in brain network dynamics enable flexible reconfiguration in response to task demands, which supports executive control, including working memory performance (Braun et al., 2015;Stanley et al., 2015). Among the proposed neurocognitive mechanisms is an increased influence of executive nodes exerting top down control (Chadick and Gazzaley, 2011), as well as suppression of task-irrelevant nodes to minimize interference which may be detrimental to performance (Tomasi et al., 2014). Indeed, withinnetwork visual and sensorimotor connectivity have been reported to decrease with task demand, with simultaneous increases in connectivity with frontoparietal and cerebellar nodes (  (1-1) The hierarchical clustering of the components, based on the partial correlations between the difference weights (n-back minus rest). (1-2) Visualization of the IC spatial maps. (1-3) Visualization of statistics reflecting the effect of age on EC during both n-back and rest, and the interaction between age and task engagement respectively. (1-4) The effect of task engagement on EC for the ICs, and effect of EC on hit-rate, false positives and d-prime. (1-5) The significant effects of age on task activation, and the effect of task engagement. Summary of statistics for task effect is displayed in table 2, age effects in table 3. (1-6) The IC-subject weight correlation-matrix. Correlations across both conditions are displayed below the diagonal and was used for the hierarchical clustering displayed in 1-1. Above the diagonal are the correlations between IC difference-weights (n-back minus rest).
facilitates both suppression and biasing of sensory input. Thalamus plays a key role in the gating of sensory input toward higher cortical network structures (Halassa and Kastner, 2017;Hwang et al., 2017). Interestingly, reduced connectivity between specific subclusters of the thalamus and visual and motor cortices in response to task demands has been reported, also in absence of altered FC with executive or DMN networks (Fan et al., 2015). In light of these observations, our findings indicate that task engagement reduces thalamic FC with cortical areas, facilitating the suppression and filtering of sensory signal during task engagement.
Our results align with previous findings highlighting the importance of increased influence of executive and cerebellar nodes during task engagement, as higher cinguloopercular, cerebellar and frontoparietal centrality was associated with better task performance, and higher centrality in visual and sensorimotor nodes was associated with lower performance. Thalamus in particular was further associated with decrease in centrality during task engagement reflecting a drop in global FC, and the degree of differentiation between task and rest EC was positively associated with performance. Indeed, our results support both the interference hypothesis (D'Esposito et al., 1999), suggesting that performance during the n-back task is supported by increased influence from frontoparietal and cingulo-opercular nodes, providing top-down attentional control, as well as the sensory gating hypothesis (Cromwell et al., 2008), reflected through the decreased thalamic EC, which was associated with better performance.
Notably, to a large extent, the centrality pattern observed during task engagement paralleled the observed activation patterns. Studying activation and connectivity across several cognitive tasks, Gratton et al. (2016) reported coordinated changes in connectivity and activation, in which an increase in both measures were observed for brain regions involved in active task control. This supports the notion of connectivity and activation as coordinated aspects of brain functioning. Overall, our results comply with these findings of coordinated activity and network changes for the frontoparietal and the cingulo-opercular-network. However, while the DMN displayed a task-related increase in centrality, it also displayed the commonly observed deactivation in response to task demands. This could reflect a coordinated disengagement, as previous studies have reported task-driven increases in internal DMN connectivity and coherence during cognitive engagement (Anticevic et al., 2012;Goparaju et al., 2014). On the other hand, recent studies using a graph theoretical approach alternatively suggest that the DMN dynamically changes its cross talk with other networks during task execution (Vatansever et al., 2015a,b) and displays larger integration with the full connectome. In particular, variability in functional connectivity between executive networks and DMN during an executive task was shown to be positively associated with performance (Douw et al., 2016). Taken together, there is still a large degree of uncertainty regarding DMN involvement during task engagement, and further studies are needed to gain a better understanding of the DMN dynamics and its implications for brain development and cognitive functioning.

Age-related differences in activation
In line with our expectations and previous activation analysis in an overlapping sample (Satterthwaite et al., 2013), our analysis revealed age-related increases in the activation of frontoparietal components and increased deactivation of visual and DMN components, suggesting increased differentiation of these networks with age in response to task demands. We also observed additional age-related effects in cingulo-opercular and subcortical components, which were not reported earlier. This discrepancy may be related to higher sensitivity for the current ICA-based approach compared to the voxel-wise analysis, allowing for higher component-wise specificity and effectively reducing number of comparisons.

Age-related differences in centrality
For the connectivity and graph-based analysis, we found increasing centrality with increasing age in key cingulo-opercular components, as well as posterior cingulate, and superior frontal gyrus, reflecting increased centrality in task positive visual and sensorimotor components, along with reductions in thalamus, cerebellum and temporal gyrus during both n-back and rest. We found no significant interaction effects between task and age. Al-though one should be cautious when interpreting null findings, the lack of significant interactions suggests that age-related differences in network characteristics are not strongly dependent on cognitive load. A recent largescale developmental study demonstrated that maturation of the structure and organization of the brain functional connectome is characterized by increasing segregation of large-scale functional networks, along with differentiation of the various networks' influence on the full connectome (Gu et al., 2015). Indeed, converging evidence points toward an increasing global influence from DMN in general, and cingulo-opercular network nodes during cognitive processing, as hallmarks of healthy development (Gu et al., 2015;Marek et al., 2015). Specifically, Marek et al. (2015) reported that maturational increase in inhibitory control during early adolescence is accompanied by increased between-network connectivity, in particular within the sensorimotor, visual, and cingulo-opercular networks. Furthermore, temporal variability in cinguloopercular connectivity has been suggested to serve as a proxy for brain maturation (Sato et al., 2015b).
The importance of the cingulo-opercular network for cognitive functioning is further supported by studies investigating both healthy adults and clinical groups. Cingulo-opercular activity in healthy adults is associated with the ability to maintain sustained attention and integrate feedback (Cocchi et al., 2013), and within-network efficiency measured as degree of edge density has been linked to cognitive flexibility (Sheffield et al., 2016). Complementary, decreased functional connectivity between cingulo-opercular, sensorimotor, and DMN nodes with aging has been linked to deteriorated cognitive performance (Meier et al., 2012), and lesions targeting hub cingulo-opercular nodes are associated with greater cognitive impairment compared to non-hub nodes (Warren et al., 2014). Indeed, structural abnormalities in hubs has been found to be overrepresented across a range of brain disorders , and abnormalities in anterior cingulate of the cingulo-opercular network was proposed as one of the key brain predictors for schizophrenia. Identifying deviant hub development may show promise as a biomarker for risk of future psychopathology in the adolescent brain, as well as being a predictive marker for prognosis and the degree of symptom severity following brain insults.
Our results only displayed maturational effect in one of six frontoparietal components, and thus revealed sparse age associations with the centrality of the frontoparietal network. Recent studies suggest that maturation of task positive networks may be better reflected by dynamic rather than static measures of functional connectivity, and that development as well as aging are accompanied by changes in dynamic transitions between brain states, operationalized through temporal variability in betweennetwork connectivity (Ryali et al., 2016;Córdova-Palomera et al., 2017;Faghiri et al., 2017). Thus, graph and centrality measures based on a static apprehension of the brain connectome may fail to capture key aspects of maturational differences in the brain temporal dynamics.
The most prominent age-related reduction in centrality was seen in thalamus and cerebellum, converging with previous reports (Baker et al., 2015;Sato et al., 2015a;Petrican et al., 2017). Subcortical nodes have been reported to display enhanced connectivity with prefrontal nodes during adolescence (Ezekiel et al., 2013;Sato et al., 2015a). In conjunction, these findings suggest that subcortical networks reduce its influence throughout adolescence, through a drop in internal connectivity, shifting toward higher influence from frontal cortices and increased top-down control (Baker et al., 2015). Indeed, aberrant thalamo-cortical connectivity in adults have been linked to mental disorders (Wang et al., 2015;Skåtun et al., 2018), where the suggested mechanisms include altered inhibition of sensory signaling. In line with the notion that cognitive and emotional development during childhood and adolescence is accompanied and supported by substantial tuning and specialization of brain network activation and connectivity, it was recently demonstrated that the brain connectome develops into a more stable and individualized pattern, and that a delay in this development is associated with increased burden of mental health issues (Kaufmann et al., 2017a).The balance between maintaining a stable and idiosyncratic connectivity pattern across contexts on one hand and brain network flexibility in response to dynamic contextual changes on the other may reflect a key characteristic of a healthy mind. Delineating the maturational trajectories of dynamic shifts in hub properties between mental states may provide important clues for healthy and pathologic brain development, as well as shed light on adult neuronal pathology. Indeed, degree of damage to hub locations . Results from full-brain voxel-wise analysis. (4-1) Significant changes in EC associated with differentiation between task and rest across the age-span. (4-2) Age-effect in the n-back data and (4-3) Age-effect in the rest-data. All maps were corrected for multiple comparisons using permutation testing (n ϭ 5000) and TFCE with an alpha of 0.05 (FWE, two-tailed). Only significant results are displayed.
have been shown to be a key predictor of cognitive impairment (Warren et al., 2014).

Limitations
The current study employed a cross-sectional design. This limits the conclusions that can be made, as there is no measure of within-person changes, and cohort effects cannot be ruled out. To further differentiate healthy and pathologic development, longitudinal studies will be of importance to better identify markers for delayed and aberrant maturation. Furthermore, the mounting evidence for misinterpreted motion-driven age effects in previous developmental studies on FC (Power et al., 2012), calls for caution when interpreting maturational findings. Although we implemented stringent procedures for data denoising (Kaufmann et al., 2017a), it is not possible to completely rule out residual motion effects as a confounder regarding age-related findings.

Conclusion
Our study adds to the mounting evidence supporting the importance of hub development for healthy maturation during adolescence. The activation and centrality of cingulo-opercular components were associated with task engagement, task performance, as well as increasing age. This supports the importance of the hub properties of cerebellum and cingulo-opercular network as key to healthy development of executive capacities, as well as its vulnerability to injuries and deviant development.