A network correspondence toolbox for quantitative evaluation of novel neuroimaging results.

Decades of neuroscience research has shown that macroscale brain dynamics can be reliably decomposed into a subset of large-scale functional networks, but the specific spatial topographies of these networks and the names used to describe them can vary across studies. Such discordance has hampered interpretation and convergence of research findings across the field. To address this problem, we have developed the Network Correspondence Toolbox (NCT) to permit researchers to examine and report spatial correspondence between their novel neuroimaging results and sixteen widely used functional brain atlases, consistent with recommended reporting standards developed by the Organization for Human Brain Mapping. The atlases included in the toolbox show some topographical convergence for specific networks, such as those labeled as default or visual. Network naming varies across atlases, particularly for networks spanning frontoparietal association cortices. For this reason, quantitative comparison with multiple atlases is recommended to benchmark novel neuroimaging findings. We provide several exemplar demonstrations using the Human Connectome Project task fMRI results and UK Biobank independent component analysis maps to illustrate how researchers can use the NCT to report their own findings through quantitative evaluation against multiple published atlases. The NCT provides a convenient means for computing Dice coefficients with spin test permutations to determine the magnitude and statistical significance of correspondence among user-defined maps and existing atlas labels. The NCT also includes functionality to incorporate additional atlases in the future. The adoption of the NCT will make it easier for network neuroscience researchers to report their findings in a standardized manner, thus aiding reproducibility and facilitating comparisons between studies to produce interdisciplinary insights.

Engaging dystonia networks with subthalamic stimulation.

Deep brain stimulation is a viable and efficacious treatment option for dystonia. While the internal pallidum serves as the primary target, more recently, stimulation of the subthalamic nucleus (STN) has been investigated. However, optimal targeting within this structure and its complex surroundings have not been studied in depth. Indeed, multiple historical targets that have been used for surgical treatment of dystonia are directly adjacent to the STN. Further, multiple types of dystonia exist, and outcomes are variable, suggesting that not all types would profit maximally from the exact same target. Therefore, a thorough investigation of the neural substrates underlying effects on dystonia symptoms is warranted. Here, we analyze a multi-center cohort of isolated dystonia patients with subthalamic implantations (N = 58) and relate their stimulation sites to improvement of appendicular and cervical symptoms as well as blepharospasm. Stimulation of the ventral oral posterior nucleus of thalamus and surrounding regions was associated with improvement in cervical dystonia, while stimulation of the dorsolateral STN was associated with improvement in limb dystonia and blepharospasm. This dissociation was also evident for structural connectivity, where the cerebellothalamic, corticospinal and pallidosubthalamic tracts were associated with improvement of cervical dystonia, while hyperdirect and subthalamopallidal pathways were associated with alleviation of limb dystonia and blepharospasm. Importantly, a single well-placed electrode may reach the three optimal target sites. On the level of functional networks, improvement of limb dystonia was correlated with connectivity to the corresponding somatotopic regions in primary motor cortex, while alleviation of cervical dystonia was correlated with connectivity to the recently described 'action-mode' network that involves supplementary motor and premotor cortex. Our findings suggest that different types of dystonia symptoms are modulated via distinct networks. Namely, appendicular dystonia and blepharospasm are improved with modulation of the basal ganglia, and, in particular, the subthalamic circuitry, including projections from the primary motor cortex. In contrast, cervical dystonia was more responsive when engaging the cerebello-thalamo-cortical circuit, including direct stimulation of ventral thalamic nuclei. These findings may inform DBS targeting and image-based programming strategies for patient-specific treatment of dystonia.

Principles of cortical areas and their implications for neuroimaging.

Cortical organization should constrain the study of how the brain performs behavior and cognition. A fundamental concept in cortical organization is that of arealization: that the cortex is parceled into discrete areas. In part one of this report, we review how non-human animal studies have illuminated principles of cortical arealization by revealing: (1) what defines a cortical area, (2) how cortical areas are formed, (3) how cortical areas interact with one another, and (4) what "computations" or "functions" areas perform. In part two, we discuss how these principles apply to neuroimaging research. In doing so, we highlight several examples where the commonly accepted interpretation of neuroimaging observations requires assumptions that violate the principles of arealization, including nonstationary areas that move on short time scales, large-scale gradients as organizing features, and cortical areas with singular functionality that perfectly map psychological constructs. Our belief is that principles of neurobiology should strongly guide the nature of computational explanations.

Proceedings of the 11th Annual Deep Brain Stimulation Think Tank: pushing the forefront of neuromodulation with functional network mapping, biomarkers for adaptive DBS, bioethical dilemmas, AI-guided neuromodulation, and translational advancements.

The Deep Brain Stimulation (DBS) Think Tank XI was held on August 9-11, 2023 in Gainesville, Florida with the theme of "Pushing the Forefront of Neuromodulation". The keynote speaker was Dr. Nico Dosenbach from Washington University in St. Louis, Missouri. He presented his research recently published in Nature inn a collaboration with Dr. Evan Gordon to identify and characterize the somato-cognitive action network (SCAN), which has redefined the motor homunculus and has led to new hypotheses about the integrative networks underpinning therapeutic DBS. The DBS Think Tank was founded in 2012 and provides an open platform where clinicians, engineers, and researchers (from industry and academia) can freely discuss current and emerging DBS technologies, as well as logistical and ethical issues facing the field. The group estimated that globally more than 263,000 DBS devices have been implanted for neurological and neuropsychiatric disorders. This year's meeting was focused on advances in the following areas: cutting-edge translational neuromodulation, cutting-edge physiology, advances in neuromodulation from Europe and Asia, neuroethical dilemmas, artificial intelligence and computational modeling, time scales in DBS for mood disorders, and advances in future neuromodulation devices.

Intrusive Traumatic Re-Experiencing Domain: Functional Connectivity Feature Classification by the ENIGMA PTSD Consortium.

Background: Intrusive traumatic re-experiencing domain (ITRED) was recently introduced as a novel perspective on posttraumatic psychopathology, proposing to focus research of posttraumatic stress disorder (PTSD) on the unique symptoms of intrusive and involuntary re-experiencing of the trauma, namely, intrusive memories, nightmares, and flashbacks. The aim of the present study was to explore ITRED from a neural network connectivity perspective. Methods: Data were collected from 9 sites taking part in the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) PTSD Consortium (n= 584) and included itemized PTSD symptom scores and resting-state functional connectivity (rsFC) data. We assessed the utility of rsFC in classifying PTSD, ITRED-only (no PTSD diagnosis), and trauma-exposed (TE)-only (no PTSD or ITRED) groups using a machine learning approach, examining well-known networks implicated in PTSD. A random forest classification model was built on a training set using cross-validation, and the averaged cross-validation model performance for classification was evaluated using the area under the curve. The model was tested using a fully independent portion of the data (test dataset), and the test area under the curve was evaluated. Results: rsFC signatures differentiated TE-only participants from PTSD and ITRED-only participants at about 60% accuracy. Conversely, rsFC signatures did not differentiate PTSD from ITRED-only individuals (45% accuracy). Common features differentiating TE-only participants from PTSD and ITRED-only participants mainly involved default mode network-related pathways. Some unique features, such as connectivity within the frontoparietal network, differentiated TE-only participants from one group (PTSD or ITRED-only) but to a lesser extent from the other group. Conclusions: Neural network connectivity supports ITRED as a novel neurobiologically based approach to classifying posttrauma psychopathology.

MRI economics: Balancing sample size and scan duration in brain wide association studies.

A pervasive dilemma in neuroimaging is whether to prioritize sample size or scan duration given fixed resources. Here, we systematically investigate this trade-off in the context of brain-wide association studies (BWAS) using resting-state functional magnetic resonance imaging (fMRI). We find that total scan duration (sample size × scan duration per participant) robustly explains individual-level phenotypic prediction accuracy via a logarithmic model, suggesting that sample size and scan duration are broadly interchangeable. The returns of scan duration eventually diminish relative to sample size, which we explain with principled theoretical derivations. When accounting for fixed costs associated with each participant (e.g., recruitment, non-imaging measures), we find that prediction accuracy in small-scale BWAS might benefit from much longer scan durations (>50 min) than typically assumed. Most existing large-scale studies might also have benefited from smaller sample sizes with longer scan durations. Both logarithmic and theoretical models of the relationships among sample size, scan duration and prediction accuracy explain well-predicted phenotypes better than poorly-predicted phenotypes. The logarithmic and theoretical models are also undermined by individual differences in brain states. These results replicate across phenotypic domains (e.g., cognition and mental health) from two large-scale datasets with different algorithms and metrics. Overall, our study emphasizes the importance of scan time, which is ignored in standard power calculations. Standard power calculations inevitably maximize sample size at the expense of scan duration. The resulting prediction accuracies are likely lower than would be produced with alternate designs, thus impeding scientific discovery. Our empirically informed reference is available for future study design: WEB_APPLICATION_LINK.

Functional parcellation of the neonatal cortical surface.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that existing parcellations, including surface-based parcels derived from older samples as well as volume-based neonatal parcels, are a poor fit for neonatal surface data. We next derive a set of 283 cortical surface parcels from a sample of n = 261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

The brainstem's red nucleus was evolutionarily upgraded to support goal-directed action.

The red nucleus is a large brainstem structure that coordinates limb movement for locomotion in quadrupedal animals (Basile et al., 2021). The humans red nucleus has a different pattern of anatomical connectivity compared to quadrupeds, suggesting a unique purpose (Hatschek, 1907). Previously the function of the human red nucleus remained unclear at least partly due to methodological limitations with brainstem functional neuroimaging (Sclocco et al., 2018). Here, we used our most advanced resting-state functional connectivity (RSFC) based precision functional mapping (PFM) in highly sampled individuals (n = 5) and large group-averaged datasets (combined N ~ 45,000), to precisely examine red nucleus functional connectivity. Notably, red nucleus functional connectivity to motor-effector networks (somatomotor hand, foot, and mouth) was minimal. Instead, red nucleus functional connectivity along the central sulcus was specific to regions of the recently discovered somato-cognitive action network (SCAN; (Gordon et al., 2023)). Outside of primary motor cortex, red nucleus connectivity was strongest to the cingulo-opercular (CON) and salience networks, involved in action/cognitive control (Dosenbach et al., 2007; Newbold et al., 2021) and reward/motivated behavior (Seeley, 2019), respectively. Functional connectivity to these two networks was organized into discrete dorsal-medial and ventral-lateral zones. Red nucleus functional connectivity to the thalamus recapitulated known structural connectivity of the dento-rubral thalamic tract (DRTT) and could prove clinically useful in functionally targeting the ventral intermediate (VIM) nucleus. In total, our results indicate that far from being a 'motor' structure, the red nucleus is better understood as a brainstem nucleus for implementing goal-directed behavior, integrating behavioral valence and action plans.

Disuse-driven plasticity in the human thalamus and putamen.

Motor adaptation in cortico-striato-thalamo-cortical loops has been studied mainly in animals using invasive electrophysiology. Here, we leverage functional neuroimaging in humans to study motor circuit plasticity in the human subcortex. We employed an experimental paradigm that combined two weeks of upper-extremity immobilization with daily resting-state and motor task fMRI before, during, and after the casting period. We previously showed that limb disuse leads to decreased functional connectivity (FC) of the contralateral somatomotor cortex (SM1) with the ipsilateral somatomotor cortex, increased FC with the cingulo-opercular network (CON) as well as the emergence of high amplitude, fMRI signal pulses localized in the contralateral SM1, supplementary motor area and the cerebellum. From our prior observations, it remains unclear whether the disuse plasticity affects the thalamus and striatum. We extended our analysis to include these subcortical regions and found that both exhibit strengthened cortical FC and spontaneous fMRI signal pulses induced by limb disuse. The dorsal posterior putamen and the central thalamus, mainly CM, VLP and VIM nuclei, showed disuse pulses and FC changes that lined up with fmri task activations from the Human connectome project motor system localizer, acquired before casting for each participant. Our findings provide a novel understanding of the role of the cortico-striato-thalamo-cortical loops in human motor plasticity and a potential link with the physiology of sleep regulation. Additionally, similarities with FC observation from Parkinson Disease (PD) questions a pathophysiological link with limb disuse.

Framewise multi-echo distortion correction for superior functional MRI.

Functional MRI (fMRI) data are severely distorted by magnetic field (B0) inhomogeneities which currently must be corrected using separately acquired field map data. However, changes in the head position of a scanning participant across fMRI frames can cause changes in the B0 field, preventing accurate correction of geometric distortions. Additionally, field maps can be corrupted by movement during their acquisition, preventing distortion correction altogether. In this study, we use phase information from multi-echo (ME) fMRI data to dynamically sample distortion due to fluctuating B0 field inhomogeneity across frames by acquiring multiple echoes during a single EPI readout. Our distortion correction approach, MEDIC (Multi-Echo DIstortion Correction), accurately estimates B0 related distortions for each frame of multi-echo fMRI data. Here, we demonstrate that MEDIC's framewise distortion correction produces improved alignment to anatomy and decreases the impact of head motion on resting-state functional connectivity (RSFC) maps, in higher motion data, when compared to the prior gold standard approach (i.e., TOPUP). Enhanced framewise distortion correction with MEDIC, without the requirement for field map collection, furthers the advantage of multi-echo over single-echo fMRI.

Functional parcellation of the neonatal brain.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that adult- and older infant-derived parcels are a poor fit with neonatal data, emphasizing the need for neonatal-specific parcels. We next derive a set of 283 cortical surface parcels from a sample of n=261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

Substructure of the brain's Cingulo-Opercular network.

The Cingulo-Opercular network (CON) is an executive network of the human brain that regulates actions. CON is composed of many widely distributed cortical regions that are involved in top-down control over both lower-level (i.e., motor) and higher-level (i.e., cognitive) functions, as well as in processing of painful stimuli. Given the topographical and functional heterogeneity of the CON, we investigated whether subnetworks within the CON support separable aspects of action control. Using precision functional mapping (PFM) in 15 participants with > 5 hours of resting state functional connectivity (RSFC) and task data, we identified three anatomically and functionally distinct CON subnetworks within each individual. These three distinct subnetworks were linked to Decisions, Actions, and Feedback (including pain processing), respectively, in convergence with a meta-analytic task database. These Decision, Action and Feedback subnetworks represent pathways by which the brain establishes top-down goals, transforms those goals into actions, implemented as movements, and processes critical action feedback such as pain.

Structure-Function Coupling in Highly Sampled Individual Brains.

Structural connections (SC) between distant regions of the brain support synchronized function known as functional connectivity (FC) and give rise to the large-scale brain networks that enable cognition and behavior. Understanding how SC enables FC is important to understand how injuries to structural connections may alter brain function and cognition. Previous work evaluating whole-brain SC-FC relationships showed that SC explained FC well in unimodal visual and motor areas, but only weakly in association areas, suggesting a unimodal-heteromodal gradient organization of SC-FC coupling. However, this work was conducted in group-averaged SC/FC data. Thus, it could not account for inter-individual variability in the locations of cortical areas and white matter tracts. We evaluated the correspondence of SC and FC within three highly sampled healthy participants. For each participant, we collected 78 minutes of diffusion-weighted MRI for SC and 360 minutes of resting state fMRI for FC. We found that FC was best explained by SC in visual and motor systems, as well as in anterior and posterior cingulate regions. A unimodal-to-heteromodal gradient could not fully explain SC-FC coupling. We conclude that the SC-FC coupling of the anterior-posterior cingulate circuit is more similar to unimodal areas than to heteromodal areas. SIGNIFICANCE STATEMENT: Structural connections between distant regions of the human brain support networked function that enables cognition and behavior. Improving our understanding of how structure enables function could allow better insight into how brain disconnection injuries impair brain function.Previous work using neuroimaging suggested that structure-function relationships vary systematically across the brain, with structure better explaining function in basic visual/motor areas than in higher-order areas. However, this work was conducted in group-averaged data, which may obscure details of individual-specific structure-function relationships.Using individual-specific densely sampled neuroimaging data, we found that in addition to visual/motor regions, structure strongly predicts function in specific circuits of the higher-order cingulate gyrus. The cingulate's structure-function relationship suggests that its organization may be unique among higher-order cortical regions.

Neuroimaging-based classification of PTSD using data-driven computational approaches: A multisite big data study from the ENIGMA-PGC PTSD consortium.

BACKGROUND: Recent advances in data-driven computational approaches have been helpful in devising tools to objectively diagnose psychiatric disorders. However, current machine learning studies limited to small homogeneous samples, different methodologies, and different imaging collection protocols, limit the ability to directly compare and generalize their results. Here we aimed to classify individuals with PTSD versus controls and assess the generalizability using a large heterogeneous brain datasets from the ENIGMA-PGC PTSD Working group. METHODS: We analyzed brain MRI data from 3,477 structural-MRI; 2,495 resting state-fMRI; and 1,952 diffusion-MRI. First, we identified the brain features that best distinguish individuals with PTSD from controls using traditional machine learning methods. Second, we assessed the utility of the denoising variational autoencoder (DVAE) and evaluated its classification performance. Third, we assessed the generalizability and reproducibility of both models using leave-one-site-out cross-validation procedure for each modality. RESULTS: We found lower performance in classifying PTSD vs. controls with data from over 20 sites (60 % test AUC for s-MRI, 59 % for rs-fMRI and 56 % for d-MRI), as compared to other studies run on single-site data. The performance increased when classifying PTSD from HC without trauma history in each modality (75 % AUC). The classification performance remained intact when applying the DVAE framework, which reduced the number of features. Finally, we found that the DVAE framework achieved better generalization to unseen datasets compared with the traditional machine learning frameworks, albeit performance was slightly above chance. CONCLUSION: These results have the potential to provide a baseline classification performance for PTSD when using large scale neuroimaging datasets. Our findings show that the control group used can heavily affect classification performance. The DVAE framework provided better generalizability for the multi-site data. This may be more significant in clinical practice since the neuroimaging-based diagnostic DVAE classification models are much less site-specific, rendering them more generalizable.

Controversies and progress on standardization of large-scale brain network nomenclature.

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macroscale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field toward standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including (1) network scale, resolution, and hierarchies; (2) interindividual variability of networks; (3) dynamics and nonstationarity of networks; (4) consideration of network affiliations of subcortical structures; and (5) consideration of multimodal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

Psilocybin desynchronizes brain networks.

1The relationship between the acute effects of psychedelics and their persisting neurobiological and psychological effects is poorly understood. Here, we tracked brain changes with longitudinal precision functional mapping in healthy adults before, during, and for up to 3 weeks after oral psilocybin and methylphenidate (17 MRI visits per participant) and again 6+ months later. Psilocybin disrupted connectivity across cortical networks and subcortical structures, producing more than 3-fold greater acute changes in functional networks than methylphenidate. These changes were driven by desynchronization of brain activity across spatial scales (area, network, whole brain). Psilocybin-driven desynchronization was observed across association cortex but strongest in the default mode network (DMN), which is connected to the anterior hippocampus and thought to create our sense of self. Performing a perceptual task reduced psilocybin-induced network changes, suggesting a neurobiological basis for grounding, connecting with physical reality during psychedelic therapy. The acute brain effects of psilocybin are consistent with distortions of space-time and the self. Psilocybin induced persistent decrease in functional connectivity between the anterior hippocampus and cortex (and DMN in particular), lasting for weeks but normalizing after 6 months. Persistent suppression of hippocampal-DMN connectivity represents a candidate neuroanatomical and mechanistic correlate for psilocybin's pro-plasticity and anti-depressant effects.

Expansion of a frontostriatal salience network in individuals with depression.

Hundreds of neuroimaging studies spanning two decades have revealed differences in brain structure and functional connectivity in depression, but with modest effect sizes, complicating efforts to derive mechanistic pathophysiologic insights or develop biomarkers. 1 Furthermore, although depression is a fundamentally episodic condition, few neuroimaging studies have taken a longitudinal approach, which is critical for understanding cause and effect and delineating mechanisms that drive mood state transitions over time. The emerging field of precision functional mapping using densely-sampled longitudinal neuroimaging data has revealed unexpected, functionally meaningful individual differences in brain network topology in healthy individuals, 2-5 but these approaches have never been applied to individuals with depression. Here, using precision functional mapping techniques and 11 datasets comprising n=187 repeatedly sampled individuals and >21,000 minutes of fMRI data, we show that the frontostriatal salience network is expanded two-fold in most individuals with depression. This effect was replicable in multiple samples, including large-scale, group-average data (N=1,231 subjects), and caused primarily by network border shifts affecting specific functional systems, with three distinct modes of encroachment occurring in different individuals. Salience network expansion was unexpectedly stable over time, unaffected by changes in mood state, and detectable in children before the subsequent onset of depressive symptoms in adolescence. Longitudinal analyses of individuals scanned up to 62 times over 1.5 years identified connectivity changes in specific frontostriatal circuits that tracked fluctuations in specific symptom domains and predicted future anhedonia symptoms before they emerged. Together, these findings identify a stable trait-like brain network topology that may confer risk for depression and mood-state dependent connectivity changes in frontostriatal circuits that predict the emergence and remission of depressive symptoms over time.

Brain hubs defined in the group do not overlap with regions of high inter-individual variability.

Connector 'hubs' are brain regions with links to multiple networks. These regions are hypothesized to play a critical role in brain function. While hubs are often identified based on group-average functional magnetic resonance imaging (fMRI) data, there is considerable inter-subject variation in the functional connectivity profiles of the brain, especially in association regions where hubs tend to be located. Here we investigated how group hubs are related to locations of inter-individual variability. To answer this question, we examined inter-individual variation at group-level hubs in both the Midnight Scan Club and Human Connectome Project datasets. The top group hubs defined based on the participation coefficient did not overlap strongly with the most prominent regions of inter-individual variation (termed 'variants' in prior work). These hubs have relatively strong similarity across participants and consistent cross-network profiles, similar to what was seen for many other areas of cortex. Consistency across participants was further improved when these hubs were allowed to shift slightly in local position. Thus, our results demonstrate that the top group hubs defined with the participation coefficient are generally consistent across people, suggesting they may represent conserved cross-network bridges. More caution is warranted with alternative hub measures, such as community density (which are based on spatial proximity to network borders) and intermediate hub regions which show higher correspondence to locations of individual variability.

Thalamo-cortical and cerebello-cortical functional connectivity in development.

The thalamus is a critical relay center for neural pathways involving sensory, motor, and cognitive functions, including cortico-striato-thalamo-cortical and cortico-ponto-cerebello-thalamo-cortical loops. Despite the importance of these circuits, their development has been understudied. One way to investigate these pathways in human development in vivo is with functional connectivity MRI, yet few studies have examined thalamo-cortical and cerebello-cortical functional connectivity in development. Here, we used resting-state functional connectivity to measure functional connectivity in the thalamus and cerebellum with previously defined cortical functional networks in 2 separate data sets of children (7-12 years old) and adults (19-40 years old). In both data sets, we found stronger functional connectivity between the ventral thalamus and the somatomotor face cortical functional network in children compared with adults, extending previous cortico-striatal functional connectivity findings. In addition, there was more cortical network integration (i.e. strongest functional connectivity with multiple networks) in the thalamus in children than in adults. We found no developmental differences in cerebello-cortical functional connectivity. Together, these results suggest different maturation patterns in cortico-striato-thalamo-cortical and cortico-ponto-cerebellar-thalamo-cortical pathways.

Resting-state cortical hubs in youth organize into four categories.

During childhood, neural systems supporting high-level cognitive processes undergo periods of rapid growth and refinement, which rely on the successful coordination of activation across the brain. Some coordination occurs via cortical hubs-brain regions that coactivate with functional networks other than their own. Adult cortical hubs map into three distinct profiles, but less is known about hub categories during development, when critical improvement in cognition occurs. We identify four distinct hub categories in a large youth sample (n = 567, ages 8.5-17.2), each exhibiting more diverse connectivity profiles than adults. Youth hubs integrating control-sensory processing split into two distinct categories (visual control and auditory/motor control), whereas adult hubs unite under one. This split suggests a need for segregating sensory stimuli while functional networks are experiencing rapid development. Functional coactivation strength for youth control-processing hubs are associated with task performance, suggesting a specialized role in routing sensory information to and from the brain's control system.