Dissociable multi-scale patterns of development in personalized brain networks

Functional connectome through the human life span

The lifespan growth of the functional connectome remains unknown. Here, we assemble task-free functional and structural magnetic resonance imaging data from 33,250 individuals aged 32 postmenstrual weeks to 80 years from 132 global sites. We report critical inflection points in the nonlinear growth curves of the global mean and variance of the connectome, peaking in the late fourth and late third decades of life, respectively. After constructing a fine-grained, lifespan-wide suite of system-level brain atlases, we show distinct maturation timelines for functional segregation within different systems. Lifespan growth of regional connectivity is organized along a primary-to-association cortical axis. These connectome-based normative models reveal substantial individual heterogeneities in functional brain networks in patients with autism spectrum disorder, major depressive disorder, and Alzheimer's disease. These findings elucidate the lifespan evolution of the functional connectome and can serve as a normative reference for quantifying individual variation in development, aging, and neuropsychiatric disorders.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

pNet: A toolbox for personalized functional networks modeling

Personalized functional networks (FNs) derived from functional magnetic resonance imaging (fMRI) data are useful for characterizing individual variations in the brain functional topography associated with the brain development, aging, and disorders. To facilitate applications of the personalized FNs with enhanced reliability and reproducibility, we develop an open-source toolbox that is user-friendly, extendable, and includes rigorous quality control (QC), featuring multiple user interfaces (graphics, command line, and a step-by-step guideline) and job-scheduling for high performance computing (HPC) clusters. Particularly, the toolbox, named personalized functional network modeling (pNet), takes fMRI inputs in either volumetric or surface type, ensuring compatibility with multiple fMRI data formats, and computes personalized FNs using two distinct modeling methods: one method optimizes the functional coherence of FNs, while the other enhances their independence. Additionally, the toolbox provides HTML-based reports for QC and visualization of personalized FNs. The toolbox is developed in both MATLAB and Python platforms with a modular design to facilitate extension and modification by users familiar with either programming language. We have evaluated the toolbox on two fMRI datasets and demonstrated its effectiveness and user-friendliness with interactive and scripting examples. pNet is publicly available at https://github.com/MLDataAnalytics/pNet.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Release the Krakencoder: A unified brain connectome translation and fusion tool

Brain connectivity can be estimated in many ways, depending on modality and processing strategy. Here we present the Krakencoder, a joint connectome mapping tool that simultaneously, bidirectionally translates between structural (SC) and functional connectivity (FC), and across different atlases and processing choices via a common latent representation. These mappings demonstrate unprecedented accuracy and individual-level identifiability; the mapping between SC and FC has identifiability 42-54% higher than existing models. The Krakencoder combines all connectome flavors via a shared low-dimensional latent space. This "fusion" representation i) better reflects familial relatedness, ii) preserves age- and sex-relevant information and iii) enhances cognition-relevant information. The Krakencoder can be applied without retraining to new, out-of-age-distribution data while still preserving inter-individual differences in the connectome predictions and familial relationships in the latent representations. The Krakencoder is a significant leap forward in capturing the relationship between multi-modal brain connectomes in an individualized, behaviorally- and demographically-relevant way.

Functional brain connectivity in children with focal epilepsy: A systematic review of functional MRI studies

Epilepsy is increasingly recognised as a brain network disorder and many studies have investigated functional connectivity (FC) in children with epilepsy using functional MRI (fMRI). This systematic review of fMRI studies, published up to November 2023, investigated profiles of FC changes and their clinical relevance in children with focal epilepsy compared to healthy controls. A literature search in PubMed and Web of Science yielded 62 articles. We categorised the results into three groups: 1) differences in correlation-based FC between patients and controls; 2) differences in other FC measures between patients and controls; and 3) associations between FC and disease variables (for example, age of onset), cognitive and seizure outcomes. Studies revealed either increased or decreased FC across multiple brain regions in children with focal epilepsy. However, findings lacked consistency: conflicting FC alterations (decreased and increased FC) co-existed within or between brain regions across all focal epilepsy groups. The studies demonstrated overall that 1) interhemispheric connections often displayed abnormal connectivity and 2) connectivity within and between canonical functional networks was decreased, particularly for the default mode network. Focal epilepsy disrupted FC in children both locally (e.g., seizure-onset zones, or within-brain subnetworks) and globally (e.g., whole-brain network architecture). The wide variety of FC study methodologies limits clinical application of the results. Future research should employ longitudinal designs to understand the evolution of brain networks during the disease course and explore the potential of FC biomarkers for predicting cognitive and postsurgical seizure outcomes.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Development of white matter fiber covariance networks supports executive function in youth

During adolescence, the brain undergoes extensive changes in white matter structure that support cognition. Data-driven approaches applied to cortical surface properties have led the field to understand brain development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. Although white matter development also appears asynchronous, previous studies have relied largely on anatomical tract-based atlases, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Harnessing advances in diffusion modeling and machine learning, we identified 14 data-driven patterns of covarying white matter structure in a large sample of youth. Fiber covariance networks aligned with known major tracts, while also capturing distinct patterns of spatial covariance across distributed white matter locations. Most networks showed age-related increases in fiber network properties, which were also related to developmental changes in executive function. This study delineates data-driven patterns of white matter development that support cognition.

Personalized functional brain network topography is associated with individual differences in youth cognition

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

Socioeconomic resources in youth are linked to divergent patterns of network integration and segregation across the brain's transmodal axis

Socioeconomic resources (SER) calibrate the developing brain to the current context, which can confer or attenuate risk for psychopathology across the lifespan. Recent multivariate work indicates that SER levels powerfully influence intrinsic functional connectivity patterns across the entire brain. Nevertheless, the neurobiological meaning of these widespread alterations remains poorly understood, despite its translational promise for early risk identification, targeted intervention, and policy reform. In the present study, we leverage the resources of graph theory to precisely characterize multivariate and univariate associations between household SER and the functional integration and segregation (i.e., participation coefficient, within-module degree) of brain regions across major cognitive, affective, and sensorimotor systems during the resting state in 5,821 youth (ages 9-10 years) from the Adolescent Brain Cognitive Development (ABCD) Study. First, we establish that decomposing the brain into profiles of integration and segregation captures more than half of the multivariate association between SER and functional connectivity with greater parsimony (100-fold reduction in number of features) and interpretability. Second, we show that the topological effects of SER are not uniform across the brain; rather, higher SER levels are related to greater integration of somatomotor and subcortical systems, but greater segregation of default mode, orbitofrontal, and cerebellar systems. Finally, we demonstrate that the effects of SER are spatially patterned along the unimodal-transmodal gradient of brain organization. These findings provide critical interpretive context for the established and widespread effects of SER on brain organization, indicating that SER levels differentially configure the intrinsic functional architecture of developing unimodal and transmodal systems. This study highlights both sensorimotor and higher-order networks that may serve as neural markers of environmental stress and opportunity, and which may guide efforts to scaffold healthy neurobehavioral development among disadvantaged communities of youth.

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

A critical period plasticity framework for the sensorimotor-association axis of cortical neurodevelopment

To understand human brain development it is necessary to describe not only the spatiotemporal patterns of neurodevelopment but also the neurobiological mechanisms that underlie them. Human neuroimaging studies have provided evidence for a hierarchical sensorimotor-to-association (S-A) axis of cortical neurodevelopment. Understanding the biological mechanisms that underlie this program of development using traditional neuroimaging approaches has been challenging. Animal models have been used to identify periods of enhanced experience-dependent plasticity - 'critical periods' - that progress along cortical hierarchies and are governed by a conserved set of neurobiological mechanisms that promote and then restrict plasticity. In this review we hypothesize that the S-A axis of cortical development in humans is partly driven by the cascading maturation of critical period plasticity mechanisms. We then describe how recent advances in in vivo neuroimaging approaches provide a promising path toward testing this hypothesis by linking signals derived from non-invasive imaging to critical period mechanisms.

Different patterns of intrinsic functional connectivity at the default mode and attentional networks predict crystalized and fluid abilities in childhood

Crystallized abilities are skills used to solve problems based on experience, while fluid abilities are linked to reasoning without evoke prior knowledge. To what extent crystallized and fluid abilities involve dissociated or overlapping neural systems is debatable. Due to often deployed small sample sizes or different study settings in prior work, the neural basis of crystallized and fluid abilities in childhood remains largely unknown. Here we analyzed within and between network connectivity patterns from resting-state functional MRI of 2707 children between 9 and 10 years from the ABCD study. We hypothesized that differences in functional connectivity at the default mode network (DMN), ventral, and dorsal attentional networks (VAN, DAN) explain differences in fluid and crystallized abilities. We found that stronger between-network connectivity of the DMN and VAN, DMN and DAN, and VAN and DAN predicted crystallized abilities. Within-network connectivity of the DAN predicted both crystallized and fluid abilities. Our findings reveal that crystallized abilities rely on the functional coupling between attentional networks and the DMN, whereas fluid abilities are associated with a focal connectivity configuration at the DAN. Our study provides new evidence into the neural basis of child intelligence and calls for future comparative research in adulthood during neuropsychiatric diseases.

Atypical functional connectivity hierarchy in Rolandic epilepsy

Functional connectivity hierarchy is an important principle in the process of brain functional organization and an important feature reflecting brain development. However, atypical brain network hierarchy organization in Rolandic epilepsy have not been systematically investigated. We examined connectivity alteration with age and its relation to epileptic incidence, cognition, or underlying genetic factors in 162 cases of Rolandic epilepsy and 117 typically developing children, by measuring fMRI multi-axis functional connectivity gradients. Rolandic epilepsy is characterized by contracting and slowing expansion of the functional connectivity gradients, highlighting the atypical age-related change of the connectivity hierarchy in segregation properties. The gradient alterations are relevant to seizure incidence, cognition, and connectivity deficit, and development-associated genetic basis. Collectively, our approach provides converging evidence for atypical connectivity hierarchy as a system-level substrate of Rolandic epilepsy, suggesting this is a disorder of information processing across multiple functional domains, and established a framework for large-scale brain hierarchical research.

Age-related differences in resting-state functional connectivity from childhood to adolescence

The human brain is active at rest, and spontaneous fluctuations in functional MRI BOLD signals reveal an intrinsic functional architecture. During childhood and adolescence, functional networks undergo varying patterns of maturation, and measures of functional connectivity within and between networks differ as a function of age. However, many aspects of these developmental patterns (e.g. trajectory shape and directionality) remain unresolved. In the present study, we characterised age-related differences in within- and between-network resting-state functional connectivity (rsFC) and integration (i.e. participation coefficient, PC) in a large cross-sectional sample of children and adolescents (n = 628) aged 8-21 years from the Lifespan Human Connectome Project in Development. We found evidence for both linear and non-linear differences in cortical, subcortical, and cerebellar rsFC, as well as integration, that varied by age. Additionally, we found that sex moderated the relationship between age and putamen integration where males displayed significant age-related increases in putamen PC compared with females. Taken together, these results provide evidence for complex, non-linear differences in some brain systems during development.

ModelArray: An R package for statistical analysis of fixel-wise data

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Computing personalized brain functional networks from fMRI using self-supervised deep learning

A novel self-supervised deep learning (DL) method is developed to compute personalized brain functional networks (FNs) for characterizing brain functional neuroanatomy based on functional MRI (fMRI). Specifically, a DL model of convolutional neural networks with an encoder-decoder architecture is developed to compute personalized FNs directly from fMRI data. The DL model is trained to optimize functional homogeneity of personalized FNs without utilizing any external supervision in an end-to-end fashion. We demonstrate that a DL model trained on fMRI scans from the Human Connectome Project can identify personalized FNs and generalizes well across four different datasets. We further demonstrate that the identified personalized FNs are informative for predicting individual differences in behavior, brain development, and schizophrenia status. Taken together, the self-supervised DL allows for rapid, generalizable computation of personalized FNs.

Multiscale functional connectivity patterns of the aging brain learned from harmonized rsfMRI data of the multi-cohort iSTAGING study

To learn multiscale functional connectivity patterns of the aging brain, we built a brain age prediction model of functional connectivity measures at seven scales on a large fMRI dataset, consisting of resting-state fMRI scans of 4186 individuals with a wide age range (22 to 97 years, with an average of 63) from five cohorts. We computed multiscale functional connectivity measures of individual subjects using a personalized functional network computational method, harmonized the functional connectivity measures of subjects from multiple datasets in order to build a functional brain age model, and finally evaluated how functional brain age gap correlated with cognitive measures of individual subjects. Our study has revealed that functional connectivity measures at multiple scales were more informative than those at any single scale for the brain age prediction, the data harmonization significantly improved the brain age prediction performance, and the data harmonization in the functional connectivity measures' tangent space worked better than in their original space. Moreover, brain age gap scores of individual subjects derived from the brain age prediction model were significantly correlated with clinical and cognitive measures. Overall, these results demonstrated that multiscale functional connectivity patterns learned from a large-scale multi-site rsfMRI dataset were informative for characterizing the aging brain and the derived brain age gap was associated with cognitive and clinical measures.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Omnipresence of the sensorimotor-association axis topography in the human connectome

Low-dimensional representations are increasingly used to study meaningful organizational principles within the human brain. Most notably, the sensorimotor-association axis consistently explains the most variance in the human connectome as its so-called principal gradient, suggesting that it represents a fundamental organizational principle. While recent work indicates these low dimensional representations are relatively robust, they are limited by modeling only certain aspects of the functional connectivity structure. To date, the majority of studies have restricted these approaches to the strongest connections in the brain, treating weaker or negative connections as noise despite evidence of meaningful structure among them. The present work examines connectivity gradients of the human connectome across a full range of connectivity strengths and explores the implications for outcomes of individual differences, identifying potential dependencies on thresholds and opportunities to improve prediction tasks. Interestingly, the sensorimotor-association axis emerged as the principal gradient of the human connectome across the entire range of connectivity levels. Moreover, the principal gradient of connections at intermediate strengths encoded individual differences, better followed individual-specific anatomical features, and was also more predictive of intelligence. Taken together, our results add to evidence of the sensorimotor-association axis as a fundamental principle of the brain's functional organization, since it is evident even in the connectivity structure of more lenient connectivity thresholds. These more loosely coupled connections further appear to contain valuable and potentially important information that could be used to improve our understanding of individual differences, diagnosis, and the prediction of treatment outcomes.

Development of White Matter Fiber Covariance Networks Supports Executive Function in Youth

The white matter architecture of the human brain undergoes substantial development throughout childhood and adolescence, allowing for more efficient signaling between brain regions that support executive function. Increasingly, the field understands grey matter development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. While white matter development also appears asynchronous, previous studies have largely relied on anatomical atlases to characterize white matter tracts, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Here, we leveraged advances in diffusion modeling and unsupervised machine learning to delineate white matter fiber covariance networks comprised of structurally similar areas of white matter in a cross-sectional sample of 939 youth aged 8-22 years. We then evaluated associations between fiber covariance network structural properties with both age and executive function using generalized additive models. The identified fiber covariance networks aligned with the known architecture of white matter while simultaneously capturing novel spatial patterns of coordinated maturation. Fiber covariance networks showed heterochronous increases in fiber density and cross section that generally followed hierarchically organized temporal patterns of cortical development, with the greatest increases in unimodal sensorimotor networks and the most prolonged increases in superior and anterior transmodal networks. Notably, we found that executive function was associated with structural features of limbic and association networks. Taken together, this study delineates data-driven patterns of white matter network development that support cognition and align with major axes of brain maturation.

Hierarchical functional system development supports executive function

In this perspective, we describe how developmental improvements in youth executive function (EF) are supported by hierarchically organized maturational changes in functional brain systems. We first highlight evidence that functional brain systems are embedded within a hierarchical sensorimotor-association axis of cortical organization. We then review data showing that functional system developmental profiles vary along this axis: systems near the associative end become more functionally segregated, while those in the middle become more integrative. Developmental changes that strengthen the hierarchical organization of the cortex may support EF by facilitating top-down information flow and balancing within- and between-system communication. We propose a central role for attention and frontoparietal control systems in the maturation of healthy EF and suggest that reduced functional system differentiation across the sensorimotor-association axis contributes to transdiagnostic EF deficits.

Fiber-specific structural properties relate to reading skills in children and adolescents

Recent studies suggest that the cross-sectional relationship between reading skills and white matter microstructure, as indexed by fractional anisotropy, is not as robust as previously thought. Fixel-based analyses yield fiber-specific micro- and macrostructural measures, overcoming several shortcomings of the traditional diffusion tensor model. We ran a whole-brain analysis investigating whether the product of fiber density and cross-section (FDC) related to single-word reading skills in a large, open, quality-controlled dataset of 983 children and adolescents ages 6-18. We also compared FDC between participants with (n = 102) and without (n = 570) reading disabilities. We found that FDC positively related to reading skills throughout the brain, especially in left temporoparietal and cerebellar white matter, but did not differ between reading proficiency groups. Exploratory analyses revealed that among metrics from other diffusion models - diffusion tensor imaging, diffusion kurtosis imaging, and neurite orientation dispersion and density imaging - only the orientation dispersion and neurite density indexes from NODDI were associated (inversely) with reading skills. The present findings further support the importance of left-hemisphere dorsal temporoparietal white matter tracts in reading. Additionally, these results suggest that future DWI studies of reading and dyslexia should be designed to benefit from advanced diffusion models, include cerebellar coverage, and consider continuous analyses that account for individual differences in reading skill.